首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 169 毫秒
1.
陆面蒸散对气候变化的影响   总被引:24,自引:0,他引:24  
利用含有较真实的陆面过程的GOALS/LASG陆气耦合模式, 分别进行亚洲/北美洲陆面蒸散的敏感性试验来研究陆地与大气环流的相互作用. 模拟结果表明:模式气候对地表蒸散的变化是极其敏感的. 尤其是亚洲地表蒸散的变化将引起极为显著的气候效应, 若地表无蒸散将使气候在一定程度上变暖变干; 此外, 陆面蒸散的异常还通过季风降水的变化和 β 效应进一步影响副热带高压的形成和变异; 进而造成北半球甚至全球大气环流发生显著变化. 因此, 除了传统观点使人们很重视副热带高压活动对我国东部大陆夏季降水的影响以外, 陆面蒸散的异常通过季风降水的变化也会对副热带高压的活动产生明显影响. 由此, 夏季陆面蒸散及其水汽相变所致的大气内热源的变化是影响天气和气候的一个重要外强迫.  相似文献   

2.
植被恢复对黄土高原局地降水的反馈效应研究   总被引:2,自引:0,他引:2  
自开展大规模植被恢复以来,黄土高原大量坡耕地被转化为草地和林地,其植被覆盖状况得到了显著改善.然而,大规模植被恢复通过改变下垫面性质导致大气与地表间的能量和水分循环过程发生变化,进而对区域气候和局地降水产生反馈.如何定量评估该反馈效应是水文气象学领域亟需解决的重要科学问题.本研究通过构建考虑植被动态变化的陆气双向耦合模式,分析黄土高原不同土地利用类型和地表参数的时空变化趋势,在考虑植被恢复对流域水热循环和陆气交互作用影响的基础上,定量评估植被建设对黄土高原局地降水的反馈效应,并阐明其作用机理.通过设立考虑植被恢复的真实情景和不考虑植被恢复的假设情景两组模拟试验,对比分析植被建设对局地降水过程的净影响.结果表明,植被恢复对黄土高原局地降水量的增加具有积极效应.观测数据显示,2000~2015年黄土高原降水量以7.84mm a~(-2)的速率显著增加,模拟结果显示大规模植被恢复对研究区降水量增加速率的贡献约为37.4%,而外部水汽环流变化对研究区降水量增加速率的贡献约占62.6%.考虑植被恢复情景比不考虑植被恢复情景下的黄土高原年均降水量高12.4%.上述研究结果对黄土高原生态保护和高质量发展,以及植被建设的可持续性管理具有重要的理论与现实意义.  相似文献   

3.
陆面过程模式中有关土壤水热传输、植被冠层和空气动力学等过程的关键参数的不确定性严重制约着地表-大气相互作用模拟能力的提高.本文利用架设在我国吉林通榆和甘肃榆中典型半干旱区陆-气相互作用野外试验站的观测资料,结合大气边界层理论,利用多种方法系统估算了上述试验站地表空气动力学粗糙度(z0m)以及热传输附加阻尼(?B?1)的量值.结果表明,z0m在半干旱区具有明显的季节变化和年际变化特征,且在植被低矮的下垫面与现行通用的陆面模式中的默认值相差较大;而?B?1的日变化和季节变化特征明显.将修正后的z0m和?B?1参数化方案引入陆面过程模式,发现能够明显改善模式对于半干旱区地表感热通量的模拟能力.这些结果说明有必要进一步结合半干旱区的野外观测试验对陆面模式在该地区的缺省参数设置进行更广泛的评估,而基于外场观测试验和大气边界层理论估算的关键地表参数对于改进陆气相互作用的模拟体现出较大的应用潜力.  相似文献   

4.
大气和地表之间热辐射交换引起的地气温度耦合(即大气温度反馈)是影响地表能量收支平衡的重要因子.文章旨在阐述大气温度反馈机理,讨论影响其强度和空间分布的主要因子;并以全球变暖为例,论述大气温度反馈如何与外强迫和气候反馈过程耦合最终对全球增暖产生贡献.基于ERA-Interim再分析资料,利用地表反馈响应分析方法,计算大气温度反馈核,以此来阐述大气温度反馈的物理机制及其强度的空间分布与气候态温度、水汽和云水含量空间分布的关系,以及全球增暖加速期间大气温度反馈对全球平均表面温度增加的贡献.分析表明大气温度反馈过程主要通过与气候系统外强迫和内部过程的耦合作用,将各独立过程引起的地表能量收支异常信号放大.研究结果表明大气温度反馈显著放大了CO_2浓度升高、冰雪融化、水汽含量增加和海洋热量吸收减缓引起的地表增暖,削弱了云量增加引起的地表降温效应.同时,也放大了地表潜热通量增加造成的地表冷却效应.从全球平均结果来看,全球快速变暖前后,尽管外强迫和气候系统内部过程引起的全球平均总地面直接能量通量扰动为负,但大气温度反馈造成的全球平均总地面能量通量扰动却为正,且后者幅度远大于前者,这导致全球平均总地面净能量通量扰动正异常.由此可见,大气温度反馈对全球变暖起到了至关重要的作用.  相似文献   

5.
南京城市化进程对大气边界层的影响研究   总被引:11,自引:2,他引:11       下载免费PDF全文
为了研究城市化进程对城市边界层结构的影响,通过高分辨率的卫星遥感资料获得土地利用类型,以及地表反照率、植被指数等地表参数,以南京为例,运用数值模拟手段进行了研究,结果显示这是一种可行的手段.数值模拟结果表明,随着城市的发展,城市反照率减小、植被减少、地表湿度减小,蒸发耗热减少,感热通量增加140 W/m2,城市波恩比增加.地表和大气之间热量交换加强,湍流热量通量增大了一倍,湍流交换发展加剧,14:00混合层高度抬高了500 m.地表湍流水汽通量和空气中水汽含量都减小,这也使更多的热量用于加热地表和大气,使地温、气温的日变化幅度增加.  相似文献   

6.
基于卫星遥感资料的中国区域土壤湿度EnKF数据同化   总被引:6,自引:0,他引:6       下载免费PDF全文
土壤湿度在陆气相互作用过程中扮演着重要的角色,是气候、水文、农业、林业等研究中重要的地球物理参数之一.土壤湿度影响地面蒸散,径流、地表反射率、地表发射率以及地表感热和潜热通量,从而对气候有重要影响,它对大气的影响在全球尺度上仅次于海面温度,在陆地尺度其影响甚至超过海面温度.本文介绍了基于EnKF及陆面过程模型的中国区域陆面土壤湿度同化系统(CLSMDAS,China Land Soil Moisture Data Assimilation System),以及该系统应用于中国区域陆面土壤湿度同化试验的结果.CLSMDAS包括以下几个部分:1)陆面模式采用美国国家大气研究中心NCAR的陆面过程模型Community Land Model Version3.0(简写为CLM3.0);2)大气驱动场数据中的降水和地面入射太阳辐射数据来自FY2静止气象卫星每小时产品;3)陆面数据同化方法采用EnKF(Ensemble Kalman Filter)同化方法;4)观测数据包括AMSR-E卫星反演土壤湿度产品以及地面土壤湿度观测资料.利用CLSMDAS对2006年6~9月的土壤湿度同化试验结果的分析表明:陆面模式模拟和同化结果都能比较合理地反映出土壤湿度时空分布,同化的土壤湿度分布与2006年8月重庆、四川发生建国以来最严重的夏伏旱有非常好的对应关系,与发生在9月的湖北东部、广西南部等地的干旱区也有非常好的对应关系.  相似文献   

7.
应用定西地区的陆面过程野外观测资料和四次EOS_MODIS卫星资料,对典型的西北半干旱雨养农业区的基本地表特征参数进行反演,在此基础上对各能量通量和蒸散量进行估算,进而分别分析其季节和空间变化.结果表明:应用卫星数据估算的各参量的相对误差基本都在20%以内,其空间和时间变化基本反映了当地的实际情况;区域各通量和日蒸散量分布极不均匀,表现为四季相差比较大的单峰型变化特征,冬季最小、其次是秋季、夏春两季最大;与面积较大的中、低植被覆盖度区域相比,小范围的高植被覆盖度区域的地表净辐射、潜热和日蒸散量相对较高,土壤热通量和显热通量相对较低,并且都在春、夏季节表现得更加显著.低、中植被覆盖度区的各能量通量的季节变化不显著的特征反映了半干旱西北雨养农业区土壤的干土层相对比较厚的特征.  相似文献   

8.
植被是重要的城市地表覆盖类型之一,它通过蒸散降温和阻挡强冷空气,对城市风温湿大气微环境和污染物扩散特征会产生显著影响.准确预测植被环境流动和标量输运特性,理解植被与大气之间湍流交换过程,对城市环境改善具有重要意义.文章建立和发展了适用于森林植被环境湍流流动和标量场演化仿真计算的大涡模拟方法,将中尺度气象预报模式与微尺度精细大涡模拟方法耦合,利用地表能量平衡关系考虑大气辐射等多物理过程,研究了典型天气条件下复杂森林下垫面大气流动问题.通过考察不同大气稳定度条件下森林植被环境流动发现,在不稳定和稳定大气情况,浮力分别起到增强和抑制大气湍流混合作用,风切变也相应减小和增加.在北京奥林匹克森林公园的多尺度模拟试验中,通过与外场观测结果比较,文章验证了耦合模型可以较好地模拟城市复杂地表上空风速、温度和相对湿度的日变化.尤其是对于风场的模拟,耦合方案明显优于传统中尺度数值模拟,这主要是由于耦合方案既考虑了外部中尺度流动的影响,又通过降尺度方法可以精细分辨城市非均匀地表粗糙元素分布,从而较为精确的复现了城市粗糙子层内复杂大气流动和湍流通量输运.  相似文献   

9.
陆面过程的研究对于更好地认识气候和天气系统的演变规律、陆地-大气水热交换过程、人类活动对气候和环境的影响等具有重要意义. 建立了综合考虑土壤冻融、土壤水汽通量、植被覆盖和陆面-大气近地层水热交换的一维冻土-植被-大气连续体模型, 模拟了固液相变、汽态水迁移、土壤水、汽、热耦合迁移等过程, 反映了液态水从未冻区向冻结区迁移、冻结及其引起的潜热迁移的冻土物理本质, 也反映了汽态水分从高温区向低温区迁移所引起的温度及水分场的变化, 并对模型进行了检验. 水分运动方程采用混合Richards方程, 可适应各种边界条件. 土壤水热传输模型求解引入了修正的Picard迭代法, 不仅使计算迭代收敛更快, 而且能更好地保证数值计算过程中的水量平衡. 结合GAME/Tibet实验1998年5月份、7月份的观测数据, 应用该模型对青藏高原安多观测点的水热交换过程进行了模拟分析. 模拟结果表明: 土壤的冻融过程对地温变化会产生负反馈作用; 若净辐射相同, 土壤表层含水量较高的情况下考虑冻结时其地热通量在冰融化时明显增加, 显热通量减少, 而潜热通量变化不大, 但是冻结时各通量的变化不明显; 而土壤发生融化时, 尽管地热通量增加, 但是地表温度仍然减小; 土壤发生冻结时, 尽管土壤负温要比不考虑冻结时高, 但整体上热通量变化不大.  相似文献   

10.
基于遥感的鄱阳湖湖区蒸散特征及环境要素影响   总被引:1,自引:0,他引:1  
蒸散是湖泊湿地生态系统水循环的重要组成部分,研究湖区地表蒸散量的时空变化对了解鄱阳湖湖区水量平衡关系具有重要意义.本研究基于MODIS数据,应用地面温度-植被指数三角关系法反演2000-2009年鄱阳湖湖区的实际蒸散量,分析湖区蒸散的时空分布特征及主要气象因子对流域蒸散的影响.结果表明:2000-2009年鄱阳湖湖区年蒸散量在685~921 mm之间,平均年蒸散量为797 mm,最大蒸散量出现在2004年.2000-2009年多年平均水体蒸发量为1107 mm,高于湖区植被蒸散量(774 mm).湖区汇水区域中蒸散量占降水的平均比例为55%,是水量平衡的主要支出项,径流系数约为0.45.湖区蒸散主要受辐射和气温的影响,月蒸散量与气温呈显著的指数相关,2007年蒸散量对温度的关系最为敏感.降水量距平与蒸散量距平的关系除2007年呈显著负相关外,其他年份相关性不显著.鄱阳湖湿地蒸散与湖泊水域面积总体呈正相关,但在水文干旱严重的2006年,当水域面积<30%时,蒸散速率随水域面积增加而减小.  相似文献   

11.
River basins in mountainous regions are characterized by strong variations in topography, vegetation, soils, climatic conditions and snow cover conditions, and all are strongly related to altitude. The high spatial variation needs to be considered when modelling hydrological processes in such catchments. A complex hydrological model, with a great potential to account for spatial variability, was developed and applied for the hourly simulation of evapotranspiration, soil moisture, water balance and the runoff components for the period 1993 and 1994 in 12 subcatchments of the alpine/pre‐alpine basin of the River Thur (area 1703 km2). The basin is located in the north‐east of the Swiss part of the Rhine Basin and has an elevation range from 350 to 2500 m a.s.l. A considerable part of the Thur Basin is high mountain area, some of it above the tree‐line and a great part of the basin is snow covered during the winter season. In the distributed hydrological model, the 12 sub‐basins of the Thur catchment were spatially subdivided into sub‐areas (hydrologically similar response units—HRUs or hydrotopes) using a GIS. Within the HRUs a hydrologically similar behaviour was assumed. Spatial interpolations of the meteorological input variables wereemployed for each altitudinal zone. The structure of the model components for snow accumulation and melt, interception, soil water storage and uptake by evapotranspiration, runoff generation and flow routing are briefly outlined. The results of the simulated potential evapotranspiration reflect the dominant role of altitudinal change in radiation and albedo of exposure, followed by the influence of slope. The actual evapotranspiration shows, in comparison with the potential evapotranspiration, a greater variability in the lower and medium altitudinal zones and a smaller variability in the upper elevation zones, which was associated with limitations of available moisture in soil and surface depression storages as well as with the evaporative demand of the local vegetation. The higher altitudinal dependency and variability of runoff results from the strong increase in precipitation and the decrease in evaporation with increased altitude. An increasing influence of snow cover on runoff as well as evapotranspiration with altitude is obvious. The computed actual evapotranspiration and runoff were evaluated against the observed values of a weighting lysimeter and against runoff hydrographs. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

12.
As an alternative to geostatistical modeling, we characterized the hydrology of a semi-arid landscape in southeastern Washington state, USA, by coupling spatial patterns identified in the distributions of relative relief and vegetation with the influence each has on soil moisture storage and evapotranspiration at the appropriate scale. Gauging precipitation, soil moisture, and evapotranspiration over a two-year period while concurrently mapping relative relief and vegetation distributions at three scales ranging from centimeters to 90 m, we determined that soil moisture and soil moisture storage are significantly greater in topographic concavities than in convexities at the microrelief (20–50 cm) scale but are not significantly different in relief features at larger scales. A generalized microrelief surface produced using a two-dimensional Fourier transformation provided a good representation of the distribution of soil moisture within microrelief when scaled to soil moisture values. Applying a spatial point process analysis we determined that big sage are randomly distributed across the landscape at all scales, suggesting that lysimeter-derived sage evapotranspiration rates also be distributed randomly across the landscape. Where sage were not present, we applied an autoregressive moving-average model conditioned on grass lysimeter measurements to derive evapotranspiration rates. Combining these hydrologic spatial patterns derived from distributions in relief and vegetation with measured precipitation inputs and evapotranspiration outputs, we created a spatially distributed model of soil moisture which we tested against measured values over an eight-week period. The model provides accurate characterization of soil moisture, allows estimates of soil moisture between measurement points, permits extrapolation of soil moisture distributions outside the gauged area, and maintains small-scale variability when aggregating soil moisture to successively larger scales.  相似文献   

13.
Using the defined sensitivity index, the sensitivity of streamflow, evapotranspiration and soil moisture to climate change was investigated in four catchments in the Haihe River basin. Climate change contained three parts: annual precipitation and temperature change and the change of the percentage of precipitation in the flood season (Pf). With satisfying monthly streamflow simulation using the variable infiltration capacity model, the sensitivity was estimated by the change of simulated hydrological variables with hypothetical climatic scenarios and observed climatic data. The results indicated that (i) the sensitivity of streamflow would increase as precipitation or Pf increased but would decrease as temperature increased; (ii) the sensitivity of evapotranspiration and soil moisture would decrease as precipitation or temperature increased, but it to Pf varied in different catchments; and (iii) hydrological variables were more sensitive to precipitation, followed by Pf, and then temperature. The nonlinear response of streamflow, evapotranspiration and soil moisture to climate change could provide a reference for water resources planning and management under future climate change scenarios in the Haihe River basin. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

14.
Wang Yanhui 《水文研究》1992,6(2):241-251
Black locust (Robina pseudoacacia) has become one of the most important shelter species in the loess area of northwest China. This paper summarizes recent research concerning its hydrological influence, including canopy interception, litter absorption capacity, its effect on rainfall kinetic energy, infiltration rates, surface runoff, soil moisture, and evapotranspiration, and its role in soil conservation. Several predictive models are listed. on the basis of existing results, optimum characteristics for an effective plantation are defined, and problems requiring further research are identified.  相似文献   

15.
This study investigates the impact of climate change on rainfall, evapotranspiration, and discharge in northern Taiwan. The upstream catchment of the Shihmen reservoir in northern Taiwan was chosen as the study area. Both observed discharge and soil moisture were simultaneously adopted to optimize the HBV‐based hydrological model, clearly improving the simulation of the soil moisture. The delta change of monthly temperature and precipitation from the grid cell of GCMs (General Circulation Models) that is closest to the study area were utilized to generate the daily rainfall and temperature series based on a weather generating model. The daily rainfall and temperature series were further inputted into the calibrated hydrological model to project the hydrological variables. The studies show that rainfall and discharge will be increased during the wet season (May to October) and decreased during the dry season (November to April of the following year). Evapotranspiration will be increased in the whole year except in November and December. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

16.
Accurate soil moisture information is useful in agricultural practice, weather forecasting, and various hydrological applications. Although land surface modeling provides a viable approach to simulating soil moisture, many factors such as errors in the precipitation can affect the accuracy of soil moisture simulations. This paper examined how precipitation rate and evapotranspiration rate affect the accuracy of soil moisture simulation using simple biosphere model with and without data assimilation through ensemble Kalman filter (EnKF). For each of the two variables, seven levels of relative errors (?20, ?10, ?5, 0, 5, 10 and 20 %) were introduced independently, thus a total of 49 combined cases were investigated. Observations from Wudaogou Hydrology Experimental site in the Huaihe River basin, China, were used to drive and verify the simulations. Results indicate that when the error of precipitation rate is within 10 % of the observations, the resulting error in soil moisture simulations is less significant and manageable, thus the simulated precipitation can be used to drive hydrological models in poorly gauged catchments when observations are not available. When the error of evapotranspiration rate is within 20 % of the observations, which is partly caused by model structural and parameterization errors, its impact on soil moisture simulation is less significant and can be acceptable. This study also demonstrated that the EnKF can perform consistently well to improve soil moisture simulation with less sensitivity to precipitation errors.  相似文献   

17.
A conceptual water‐balance model was modified from a point application to be distributed for evaluating the spatial distribution of watershed water balance based on daily precipitation, temperature and other hydrological parameters. The model was calibrated by comparing simulated daily variation in soil moisture with field observed data and results of another model that simulates the vertical soil moisture flow by numerically solving Richards' equation. The impacts of soil and land use on the hydrological components of the water balance, such as evapotranspiration, soil moisture deficit, runoff and subsurface drainage, were evaluated with the calibrated model in this study. Given the same meteorological conditions and land use, the soil moisture deficit, evapotranspiration and surface runoff increase, and subsurface drainage decreases, as the available water capacity of soil increases. Among various land uses, alfalfa produced high soil moisture deficit and evapotranspiration and lower surface runoff and subsurface drainage, whereas soybeans produced an opposite trend. The simulated distribution of various hydrological components shows the combined effect of soil and land use. Simulated hydrological components compare well with observed data. The study demonstrated that the distributed water balance approach is efficient and has advantages over the use of single average value of hydrological variables and the application at a single point in the traditional practice. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

18.
Abstract

An integrated model, combining a surface energy balance system, an LAI-based interception model and a distributed monthly water balance model, was developed to predict hydrological impacts of land-use/land-cover change (LUCC) in the East River basin, China, with the aid of GIS/RS. The integrated model is a distributed model that not only accounts for spatial variations in basin terrain, rainfall and soil moisture, but also considers spatial and temporal variation of vegetation cover and evapotranspiration (ET), in particular, thus providing a powerful tool for investigating the hydrological impact of LUCC. The model was constructed using spatial data on topography, soil types and vegetation characteristics together with time series of precipitation from 170 stations in the basin. The model was calibrated and validated based on river discharge data from three stations in the basin for 21 years. The calibration and validation results suggested that the model is suitable for application in the basin. The results show that ET has a positive relationship with LAI (leaf area index), while runoff has a negative relationship with LAI in the same climatic zone that can be described by the surface energy balance and water balance equation. It was found that deforestation would cause an increase in annual runoff and a decrease in annual ET in southern China. Monthly runoff for different land-cover types was found to be inversely related to ET. Also, for most of the scenarios, and particularly for grassland and cropland, the most significant changes occurred in the rainy season, indicating that deforestation would cause a significant increase in monthly runoff in that season in the East River basin. These results are important for water resources management and environmental change monitoring.
Editor Z.W. Kundzewicz  相似文献   

19.
ABSTRACT

Climate change alters hydrological processes and results in more extreme hydrological events, e.g. flooding and drought, which threaten human livelihoods. In this study, the large-scale distributed variable infiltration capacity (VIC) model was used to simulate future hydrological processes in the Yarlung Zangbo River basin (YZRB), China, with a combination of the CMIP5 (Coupled Model Intercomparison Project, fifth phase) and MIROC5 (Model for Interdisciplinary Research on Climate, fifth version) datasets. The results indicate that the performance of the VIC model is suitable for the case study, and the variation in runoff is remarkably consistent with that of precipitation, which exhibits a decreasing trend for the period 2046–2060 and an increasing trend for 2086–2100. The seasonality of runoff is evident, and substantial increases are projected for spring runoff, which might result from the increase in precipitation as well as the increase in the warming-induced melting of snow, glaciers and frozen soil. Moreover, evapotranspiration exhibits an increase between 2006–2020 and 2046–2060 over the entire basin, and soil moisture decreases in upstream areas and increases in midstream and downstream areas. For 2086–2100, both evapotranspiration and soil moisture increase slightly in the upstream and midstream areas and decrease slightly in the downstream area. The findings of this study could provide references for runoff forecasting and ecological protection for similar studies in the future.  相似文献   

20.
The implementation of large-scale vegetation restoration over the Chinese Loess Plateau has achieved clear improvements in vegetation fraction, as evidenced by large areas of slopes and plains being restored to grassland or forest.However, such large-scale vegetation restoration has altered land-atmosphere exchanges of water and energy, as the land surface characteristics have changed. These variations could affect regional climate, especially local precipitation. Quantitatively evaluating this feedback is an important scientific question in hydrometeorology. This study constructs a coupled land-atmosphere model incorporating vegetation dynamics, and analyzes the spatio-temporal changes of different land use types and land surface parameters over the Loess Plateau. By considering the impacts of vegetation restoration on the water-energy cycle and on land-atmosphere interactions, we quantified the feedback effect of vegetation restoration on local precipitation across the Loess Plateau, and discussed the important underlying processes. To achieve a quantitative evaluation, we designed two simulation experiments, comprising a real scenario with vegetation restoration and a hypothetical scenario without vegetation restoration.These enabled a comparison and analysis of the net impact of vegetation restoration on local precipitation. The results show that vegetation restoration had a positive effect on local precipitation over the Loess Plateau. Observations show that precipitation on the Loess Plateau increased significantly, at a rate of 7.84 mm yr~(-2), from 2000 to 2015. The simulations show that the contribution of large-scale vegetation restoration to the precipitation increase was about 37.4%, while external atmospheric circulation changes beyond the Loess Plateau contributed the other 62.6%. The average annual precipitation under the vegetation restoration scenario over the Loess Plateau was 12.4% higher than that under the scenario without vegetation restoration. The above research results have important theoretical and practical significance for the ecological protection and optimal development of the Loess Plateau, as well as the sustainable management of vegetation restoration.  相似文献   

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

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