首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 375 毫秒
1.
In North American Land Data Assimilation System Phase 2 (NLDAS‐2) Noah simulation, the NLDAS team introduced an intermediate ‘fix’ to constrain the surface exchange coefficient when the atmospheric boundary layer is stable. In the current NLDAS‐2 Noah version, this fix is used for all stable cases including snow‐free grid cells. In this study, we simply apply this fix to the grid cells in which both stable atmospheric boundary layer and snow exist simultaneously, excluding the snow‐free grid cells as we recognize that the fix in NLDAS‐2 is too strong. We conduct a 31‐year (1979–2009) NLDAS‐2 Noah interim (Noah‐I) run and use observed streamflow, evapotranspiration, land surface temperature, soil temperature, and ground heat flux to evaluate the results, including comparisons with the original NLDAS‐2 Noah run. The results show that Noah‐I has the same performance as NLDAS‐2 Noah for snow water equivalent; however, Noah‐I significantly improved the simulation of other hydrometeorological products as noted earlier when compared with NLDAS‐2 Noah and the observations. This simple modification is being included in the next Noah version used in NLDAS. The hydrometeorological products from the improved NLDAS‐2 Noah‐I are being staged on the National Centers for Environmental Prediction public server. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

2.
The North American Land Data Assimilation System project phase 2 (NLDAS‐2) has run four land surface models for a 30‐year (1979–2008) retrospective period. Land surface evapotranspiration (ET) is one of the most important model outputs from NLDAS‐2 for investigating land–atmosphere interaction or to monitor agricultural drought. Here, we evaluate hourly ET using in situ observations over the Southern Great Plains (Atmospheric Radiation Measurement/Cloud and Radiation Testbed network) for 1 January 1997–30 September 1999 and daily ET u‐sing in situ observations at the AmeriFlux network over the conterminous USA for an 8‐year period (2000–2007). The NLDAS‐2 models compare well against observations, with the National Centers for Environmental Prediction's Noah land surface model performing best, followed, in order, by the Variable Infiltration Capacity, Sacramento Soil Moisture Accounting, and Mosaic models. Daily evaluation across the AmeriFlux network shows that for all models, performance depends on season and vegetation type; they do better in spring and fall than in winter or summer and better for deciduous broadleaf forest and grasslands than for croplands or evergreen needleleaf forest. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

3.
This study reports results from an analysis of the relationship between atmospheric forcing and model‐simulated water and energy fluxes for the North American Land Data Assimilation System Project Phase 2 (NLDAS‐2). The relationships between mean monthly precipitation and total runoff are stronger in the Sacramento (SAC) and variable infiltration capacity (VIC) models, which grew out of the hydrological community, than in the Noah and Mosaic models, which grew out of the soil‐vegetation‐atmosphere transfer (SVAT) community. The reverse is true for the relationship between mean monthly precipitation and evapotranspiration. In addition, surface energy fluxes in VIC are less sensitive to model forcing (except for air temperature) than those in the Noah and Mosaic model. Notwithstanding these general conclusions, the relationships between forcings and model‐simulated water and energy fluxes for all models vary for different seasons, variables, and regions. These findings will ultimately inspire a combination of SVAT‐type model energy components with hydrological model water components to develop a SVAT‐hydrology model to improve both evapotranspiration and total runoff simulations. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
Ocean Dynamics - Forecasts of 10-m wind (U10) and significant wave height (Hs) from the National Centers for Environmental Prediction (NCEP) Ensemble Forecast System are evaluated using altimeter...  相似文献   

5.
The long‐term and large‐scale soil moisture (SM) record is important for understanding land atmosphere interactions and their impacts on the weather, climate, and regional ecosystem. SM products are one of the parameters used in some Earth system models, but these records require evaluation before use. The water resources on the Qinghai–Tibet Plateau (QTP) are important to the water security of billions of people in Asia. Therefore, it is necessary to know the SM conditions on the QTP. In this study, the evaluation metrics of multilayer (0–10, 10–40, and 40–100 cm) SM in different reanalysis datasets of the European Centre for Medium‐Range Weather Forecasts interim reanalysis (ERA‐Interim [ERA]), National Centers for Environmental Prediction Climate Forecast System and the Climate Forecast System version 2 (CFSv2), and China Meteorological Administration Land Data Assimilation System (CLDAS) are compared with in situ observations at 5 observation sites, which represent alpine meadow, alpine swamp meadow, alpine grassy meadow, alpine desert steppe, and alpine steppe environments during the thawing season from January 1, 2011, to December 31, 2013, on the QTP. The ERA SM remains constant at approximately 0.2 m3?m?3 at all observation sites during the entire thawing season. The CLDAS and CFSv2 SM products show similar patterns with those of the in situ SM observations during the thawing season. The CLDAS SM product performs better than the CFSv2 and ERA for all vegetation types except the alpine swamp meadow. The results indicate that the soil texture and land cover types play a more important role than the precipitation to increase the biases of the CLDAS SM product on the QTP.  相似文献   

6.
To improve simulations of regional‐scale snow processes and related cold‐season hydroclimate, the Community Land Model version 3 (CLM3), developed by the National Center for Atmospheric Research (NCAR), was coupled with the Pennsylvania State University/NCAR fifth‐generation Mesoscale Model (MM5). CLM3 physically describes the mass and heat transfer within the snowpack using five snow layers that include liquid water and solid ice. The coupled MM5–CLM3 model performance was evaluated for the snowmelt season in the Columbia River Basin in the Pacific Northwestern United States using gridded temperature and precipitation observations, along with station observations. The results from MM5–CLM3 show a significant improvement in the SWE simulation, which has been underestimated in the original version of MM5 coupled with the Noah land‐surface model. One important cause for the underestimated SWE in Noah is its unrealistic land‐surface structure configuration where vegetation, snow and the topsoil layer are blended when snow is present. This study demonstrates the importance of the sheltering effects of the forest canopy on snow surface energy budgets, which is included in CLM3. Such effects are further seen in the simulations of surface air temperature and precipitation in regional weather and climate models such as MM5. In addition, the snow‐season surface albedo overestimated by MM5–Noah is now more accurately predicted by MM5–CLM3 using a more realistic albedo algorithm that intensifies the solar radiation absorption on the land surface, reducing the strong near‐surface cold bias in MM5–Noah. The cold bias is further alleviated due to a slower snowmelt rate in MM5–CLM3 during the early snowmelt stage, which is closer to observations than the comparable components of MM5–Noah. In addition, the over‐predicted precipitation in the Pacific Northwest as shown in MM5–Noah is significantly decreased in MM5–CLM3 due to the lower evaporation resulting from the longer snow duration. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

7.
陆面过程模式是气候模式和天气模式的核心组成部分之一.在土壤—植被—大气耦合模式(Soil-PlantAtmosphere Model,SPAM)的基础上,发展了新一代北京大学陆面过程模式PKULM(Peking University Land Model).本文首先介绍了PKULM的辐射传输、湍流输送、光合作用、土壤水热输送等过程的参数化方案;采用隐式迭代计算框架,发展并应用了一个快速的线性方程组求解算法,提高了模式计算稳定性;提出并使用了二分搜索算法计算气孔阻抗,避免了CLM(Community Land Model)等使用的迭代方法在干旱区不稳定的情况,提高了模式的适用性;采用水势为基础的土壤水分扩散方程,使模式能够模拟土壤饱和区的水分输送过程,为进一步与水文过程模式耦合奠定了基础;还发展了一个地表积水与径流过程的机理模型,提高了模式对地表水分平衡过程的模拟能力;最后,使用"中国西北干旱区陆—气相互作用观测试验"平凉站的资料对模式进行了检验并与NOAH(National Center for Environmental Prediction,Oregon State University,Air Force,and Hydrology Lab model)陆面过程模式的模拟结果进行了比较,结果表明PKULM能够较好地模拟西北半干旱区农田下垫面地气交换过程.  相似文献   

8.
M. Rahman  M. Sulis  S. J. Kollet 《水文研究》2016,30(10):1563-1573
Subsurface and land surface processes (e.g. groundwater flow, evapotranspiration) of the hydrological cycle are connected via complex feedback mechanisms, which are difficult to analyze and quantify. In this study, the dual‐boundary forcing concept that reveals space–time coherence between groundwater dynamics and land surface processes is evaluated. The underlying hypothesis is that a simplified representation of groundwater dynamics may alter the variability of land surface processes, which may eventually affect the prognostic capability of a numerical model. A coupled subsurface–land surface model ParFlow.CLM is applied over the Rur catchment, Germany, and the mass and energy fluxes of the coupled water and energy cycles are simulated over three consecutive years considering three different lower boundary conditions (dynamic, constant, and free‐drainage) based on groundwater dynamics to substantiate the aforementioned hypothesis. Continuous wavelet transform technique is applied to analyze scale‐dependent variability of the simulated mass and energy fluxes. The results show clear differences in temporal variability of latent heat flux simulated by the model configurations with different lower boundary conditions at monthly to multi‐month time scales (~32–91 days) especially under soil moisture limited conditions. The results also suggest that temporal variability of latent heat flux is affected at even smaller time scales (~1–3 days) if a simple gravity drainage lower boundary condition is considered in the coupled model. This study demonstrates the importance of a physically consistent representation of groundwater dynamics in a numerical model, which may be important to consider in local weather prediction models and water resources assessments, e.g. drought prediction. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

9.
A statistical downscaling model, based on the outputs of general circulation models (GCMs) as predictors, was proposed to simulate the daily precipitations in the Shih‐Men reservoir catchment in Taiwan. The structure of the proposed downscaling model is composed of two parts: classification and regression. Predictors of classification and regression models were selected from the large‐scale weather variables in the National Centers for Environmental Prediction (NECP) reanalysis data based on statistical tests. Discriminant analysis and support vector regression (SVR) were applied to build the classification and regression models. The outputs of five atmosphere‐ocean GCMs, which are reported to have properly considered tropical cyclone information and East Asian Monsoon modelling, were used for projecting future precipitations. Data from four grids covering Taiwan were used for developing the downscaling model. The potential of the downscaling models in simulating local precipitations was evaluated, and downscaling results reveal that the proposed downscaling model can reproduce local daily precipitations from large‐scale weather variables. Projected local precipitations under two emission scenarios show that the precipitations in the study area tend to decrease. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

10.
The operational prediction of climatic variables in monthly-to-seasonal scales has been issued by National Centers for Environmental Prediction (NCEP) through Climate Forecast System model (CFSv1) since 2004. After incorporating significant changes, a new version of this model (CFSv2) was released in 2011. The present study is based on the comparative evaluation of performances of CFSv2 and CFSv1 for the southwest monsoon season (June-July-August-September, JJAS) over India with May initial condition during 1982–2009. It was observed that CFSv2 has improved over CFSv1 in simulating the observed monsoon rainfall climatology and inter annual variability. The movement of the cell of Walker circulation in years of excessive and deficient rainfall is better captured in CFSv2, as well. The observed teleconnection pattern between ISMR-sea surface temperature (SST) is also better captured in CFSv2. The overall results suggest that the changes incorporated in CFSv1 through the development of CFSv2 have resulted in an improved prediction of ISMR.21  相似文献   

11.
WRF模式不同陆面方案对一次暴雨事件模拟的影响   总被引：5，自引：1，他引：4       下载免费PDF全文
本文利用中尺度模式Weather Research and Forecasting Model (WRF) 3.1版本及National Centers for Environmental Prediction (NCEP)分析资料,就2003年6月下旬我国江淮及南方地区的强降水事件, 以24 h短期天气模拟的方式,研究了模式中四个不同陆面方案对降水模拟的影响.结果表明,此次暴雨事件模拟对不同陆面方案是比较敏感的,模拟区域内雨量级别越高,不同方案的TS评分差异就越大,较大范围雨量可存在30%的差异,四种方案的暴雨中心值可存在100%～150%的较大差别;不同陆面方案还导致了模拟平均感热通量及潜热通量的系统性差异,这些差异的分布具有地域特点;陆面方案通过两种机理对模拟降水产生重要影响,即主要影响地表蒸发量,以及主要影响低层环流及水汽辐合,从而分别影响模拟的较大范围降水(如,平均约7%、最大约30%的较大范围雨量差异)及包含模拟降水中心的较小范围暴雨(如,方案间暴雨中心雨量可存在100%～150%的较大差别).可见,不同陆面过程可从不同空间尺度、不同程度上影响暴雨天气,改进陆面方案可以提高WRF模式对暴雨的模拟能力.  相似文献   

12.
Ocean–atmosphere coupling in the Humboldt Current System (HCS) of the Southeast Pacific is studied using the Scripps Coupled Ocean–atmosphere Regional (SCOAR) model, which is used to downscale the National Center for Environmental Prediction (NCEP) Reanalysis-2 (RA2) product for the period 2000–2007 at 20-km resolution. An interactive 2-D spatial smoother within the sea-surface temperature (SST)–flux coupler is invoked in a separate run to isolate the impact of the mesoscale (～50–200 km, in the oceanic sense) SST field felt by the atmosphere in the fully coupled run. For the HCS, SCOAR produces seasonal wind stress and wind stress curl patterns that agree better with QuikSCAT winds than those from RA2. The SCOAR downscaled wind stress distribution has substantially different impacts on the magnitude and structure of wind-driven upwelling processes along the coast compared to RA2. Along coastal locations such as Arica and Taltal, SCOAR and RA2 produce seasonally opposite signs in the total wind-driven upwelling transport. At San Juan, SCOAR shows that upwelling is mainly due to coastal Ekman upwelling transport, while in RA2 upwelling is mostly attributed to Ekman pumping. Fully coupled SCOAR shows significant SST–wind stress coupling during fall and winter, while smoothed SCOAR shows insignificant coupling throughout, indicating the important role of ocean mesoscale eddies on air–sea coupling in HCS. Coupling between SST, wind speed, and latent heat flux is incoherent in large-scale coupling and full coupling mode. In contrast, coupling between these three variables is clearly identified for oceanic mesoscales, which suggests that mesoscale SST affects latent heat directly through the bulk formulation, as well as indirectly through stability changes on the overlying atmosphere, which affects surface wind speeds. The SST–wind stress and SST–heat-flux couplings, however, fail to produce a strong change in the ocean eddy statistics. No rectified effects of ocean–atmosphere coupling were identified for either the atmospheric or oceanic mean conditions, suggesting that mesoscale coupling is too weak in this region to strongly alter the basic climate state.  相似文献   

13.
The numerical modeling of the impacts of urban buildings in mesoscale meteorological models has gradually improved in recent years. Correctly representing the latent heat flux from urban surfaces is a key issue in urban land-atmosphere coupling studies but is a common weakness in current urban canopy models. Using the surface energy balance data at a height of 140 m from a 325 m meteorological tower in Beijing, we conducted a 1-year continuous off-line simulation by using a coupled land surface model and a single-layer urban canopy model and found that this model has a relatively large systematic error for simulated latent heat flux. To improve the numerical method for modeling latent heat flux from urban surfaces, we combined observational analysis and urban land surface model to derive an oasis effect coefficient for urban green areas; to develop a temporal variation formula for water availability in urban impervious surfaces; and to specify a diurnal profile and the maximum values of anthropogenic latent heat release for four seasons. These results are directly incorporated into the urban land surface model to improve model performance. In addition, this method serves as a reference for studies in other urban areas.  相似文献   

14.
Forecast of Low Visibility and Fog from NCEP: Current Status and Efforts   总被引：2，自引：0，他引：2  
Based on the visibility analysis data during November 2009 through April 2010 over North America from the Aviation Digital Database Service (ADDS), the performance of low visibility/fog predictions from the current operational 12?km-NAM, 13?km-RUC and 32?km-WRF-NMM models at the National Centers for Environmental Prediction (NCEP) was evaluated. The evaluation shows that the performance of the low visibility/fog forecasts from these models is still poor in comparison to those of precipitation forecasts from the same models. In order to improve the skill of the low visibility/fog prediction, three efforts have been made at NCEP, including application of a rule-based fog detection scheme, extension of the NCEP Short Range Ensemble Forecast System (SREF) to fog ensemble probabilistic forecasts, and a combination of these two applications. How to apply these techniques in fog prediction is described and evaluated with the same visibility analysis data over the same period of time. The evaluation results demonstrate that using the multi-rule-based fog detection scheme significantly improves the fog forecast skill for all three models relative to visibility-diagnosed fog prediction, and with a combination of both rule-based fog detection and the ensemble technique, the performance skill of fog forecasting can be further raised.  相似文献   

15.
ABSTRACT

Ballona Creek watershed in Los Angeles, California provides a unique combination of heterogeneous urban land cover, a semi-arid environment, and a large outdoor water-use flux that presents a challenge for physically-based models. We ran simulations using the Noah Land Surface Model and Parflow-Community Land Model and compared to observations of evapotranspiration (ET), runoff, and land surface temperature (LST) for the entire 11-year study period. Both models were systematically adjusted to test the impact of land cover and urban irrigation on simulation results. Monthly total runoff and ET results are greatly improved when compared to an in-situ stream gauge and meteorological tower data: from 0.64 to 0.81 for the Nash–Sutcliffe efficiency (NSE) for runoff and from a negative NSE to 0.82 for ET. The inclusion of urban irrigation in semi-arid urban environments is found to be vital, but not sufficient, for the accurate simulation of variables in the studied models.  相似文献   

16.
Abstract

The hydrological regime of a mountainous catchment, in this instance the Mesochora catchment in Central Greece, was simulated for altered climates resulting when using the Goddard Institute for Space Studies (GISS) model for carbon dioxide doubling. The catchment snow water equivalent was predicted on the basis of the snow accumulation and ablation model of the US National Weather Service River Forecast System (NWSRFS), while the catchment runoff, as well as actual evapotranspiration and soil moisture storages, were simulated through application of the soil moisture accounting model of NWSRFS. Two scenarios of monthly climate change were drawn from the GISS model, one associated with temperature and precipitation changes, while the other referred to temperature changes alone. A third hypothetical scenario with temperature and precipitation changes similar to those corresponding to the mean monthly GISS scenarios was used to test the sensitivity of the monthly climate change of the hypothetical case on catchment hydrology. All three scenarios projected decreases in average snow accumulations and in spring and summer runoff and soil moisture, as well as increases in winter runoff and soil moisture storage and spring evapotranspiration.  相似文献   

17.
The default fractional vegetation cover and terrain height were replaced by the estimated fractional vegetation cover, which was calculated by the Normalized Difference Vegetation Index(NDVI) of Earth Observing System Moderate-Resolution Imaging Spectroradiometer(EOS-MODIS) and the Digital Elevation Model of the Shuttle Radar Topography Mission(SRTM) system. The near-surface meteorological elements over northeastern China were assimilated into the three-dimensional variational data assimilation system(3DVar) module in the Weather Research and Forecasting(WRF) model. The structure and daily variations of air temperature, humidity, wind and energy fields over northeastern China were simulated using the WRF model. Four groups of numerical experiments were performed, and the simulation results were analyzed of latent heat flux, sensible heat flux, and their relationships with changes in the surface energy flux due to soil moisture and precipitation over different surfaces. The simulations were compared with observations of the stations Tongyu, Naiman, Jinzhou, and Miyun from June to August, 2009. The results showed that the WRF model achieves high-quality simulations of the diurnal characteristics of the surface layer temperature, wind direction, net radiation, sensible heat flux, and latent heat flux over semiarid northeastern China in the summer. The simulated near-surface temperature, relative humidity, and wind speed were improved in the data assimilation case(Case 2) compared with control case(Case 1). The simulated sensible heat fluxes and surface heat fluxes were improved by the land surface parameterization case(Case 3) and the combined case(Case 4). The simulated temporal variations in soil moisture over the northeastern arid areas agree well with observations in Case 4, but the simulated precipitation should be improved in the WRF model. This study could improve the land surface parameters by utilizing remote sensing data and could further improve atmospheric elements with a data assimilation system. This work provides an effective attempt at combining multi-source data with different spatial and temporal scales into numerical simulations. The assimilation datasets generated by this work can be applied to research on climate change and environmental monitoring of arid lands, as well as research on the formation and stability of climate over semiarid areas.  相似文献   

18.
Changes in climate and land use can significantly influence the hydrological cycle and hence affect water resources. Understanding the impacts of climate and land‐use changes on streamflow can facilitate development of sustainable water resources strategies. This study investigates the flow variation of the Zamu River, an inland river in the arid area of northwest China, using the Soil and Water Assessment Tool distributed hydrological model. Three different land‐use and climate‐change scenarios were considered on the basis of measured climate data and land‐use cover, and then these data were input into the hydrological model. Based on the sensitivity analysis, model calibration and verification, the hydrological response to different land‐use and climate‐change scenarios was simulated. The results indicate that the runoff varied with different land‐use type, and the runoff of the mountain reaches of the catchment increased when grassland area increased and forestland decreased. The simulated runoff increased with increased precipitation, but the mean temperature increase decreased the runoff under the same precipitation condition. Application of grey correlation analysis showed that precipitation and temperature play a critical role in the runoff of the Zamu River basin. Sensitivity analysis of runoff to precipitation and temperature by considering the 1990s land use and climate conditions was also undertaken. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

19.
Land surface energy fluxes are required in many environmental studies, including hydrology, agronomy and meteorology. Surface energy balance models simulate microscale energy exchange processes between the ground surface and the atmospheric layer near ground level. Spatial variability of energy fluxes limits point measurements to be used for larger areas. Remote sensing provides the basis for spatial mapping of energy fluxes. Remote‐sensing‐based surface energy flux‐mapping was conducted using seven Landsat images from 1997 to 2002 at four contiguous crop fields located in Polk County, northwestern Minnesota. Spatially distributed surface energy fluxes were estimated and mapped at 30 m pixel level from Landsat Thematic Mapper and Enhanced Thematic Mapper images and weather information. Net radiation was determined using the surface energy balance algorithm for land (SEBAL) procedure. Applying the two‐source energy balance (TSEB) model, the surface temperature and the latent and sensible heat fluxes were partitioned into vegetation and soil components and estimated at the pixel level. Yield data for wheat and soybean from 1997 to 2002 were mapped and compared with latent heat (evapotranspiration) for four of the fields at pixel level. The spatial distribution and the relation of latent heat flux and Bowen ratio (ratio of sensible heat to latent heat) to crop yield were studied. The root‐mean‐square error and the mean absolute percentage of error between the observed and predicted energy fluxes were between 7 and 22 W m−2 and 12 and 24% respectively. Results show that latent heat flux and Bowen ratio were correlated (positive and negative) to the yield data. Wheat and soybean yields were predicted using latent heat flux with mean R2 = 0·67 and 0·70 respectively, average residual means of −4·2 bushels/acre and 0·11 bushels/acre respectively, and average residual standard deviations of 16·2 bushels/acre and 16·6 bushels/acre respectively (1 bushel/acre ≈ 0·087 m3 ha−1). The flux estimation procedure from the SEBAL‐TSEB model was useful and applicable to agricultural fields. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

20.
The European Centre for Medium-Range Weather Forecasts Reanalysis ERA40,National Centers for Environmental Prediction(NCEP) 20th-century reanalysis,and three station observations along an Antarctic traverse from Zhongshan to Dome-A stations are used to assess 2-m temperature simulation skill of a regional climate model.This model(HIRHAM) is from the Alfred Wegener Institute for Polar and Marine Research in Germany.Results show:(1) The simulated multiyear averaged 2-m temperature field pattern is close to that of ERA40 and NCEP;(2) the cold bias relative to ERA40 over all of Antarctic regions is 1.8℃,and that to NCEP reaches 5.1℃;(3) bias of HIRHAM relative to ERA40 has seasonal variation,with a cold bias mainly in the summer,as much as 3.4℃.There is a small inland warm bias in autumn of 0.3℃.Further analysis reveals that the reason for the cold bias of 2-m temperature is that physical conditions of the near-surface boundary layer simulated by HIRHAM are different from observations:(1) During the summer,observations show that near-surface atmospheric stability conditions have both inversions and non-inversions,which is due to the existence of both positive and negative sensible heat fluxes,but HIRHAM almost always simulates a situation of inversion and negative sensible heat flux;(2) during autumn and winter,observed near-surface stability is almost always that of inversions,consistent with HIRHAM simulations.This partially explains the small bias during autumn and winter.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号