首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
利用1958—2014年47个CMIP6模式输出资料和NCEP/NCAR再分析资料,研究了模式大气中南北涛动(InterHemispheric Oscillation,IHO)的季节变化特征,且评估了CMIP6对IHO季节特征的模拟能力。结果表明:47个CMIP6模式都能模拟出IHO的季节演变特征,但模式间存在一定差异。通过比较,筛选出模拟IHO季节循环较好的16个模式,它们能成功模拟出半球大气质量的时间演变和空间结构。进一步分析表明,水汽对IHO季节变化有抵消作用且半球内部水汽质量变化可驱动越赤道质量流的产生;地表净短波辐射夏高冬低,其加热造成的水汽蒸发在水汽质量变化中起到重要作用;地表净长波辐射在春秋变化幅度较大,与大气质量逐月变化吻合。对比再分析资料表明,CMIP6模式模拟的半球大气质量的峰谷值变化有明显的月份偏差,且CMIP6模式模拟的地表气压异常值的偏差主要出现在北太平洋、欧亚大陆、南半球中纬度地区和两极极区,模拟的南北半球的蒸发和降水量、赤道风场、地表净长波和短波辐射通量等均存在明显的偏差。  相似文献   

2.
赵丹  张丽霞  周天军 《大气科学》2022,46(3):557-572
本文基于观测和再分析资料,采用Brubaker二元模型评估了第六次国际耦合模式比较计划(CMIP6)中19个模式对中国东部季风区气候态水循环过程的模拟能力,并分析了模拟误差来源。结果表明,CMIP6模式集合平均(MME)能够合理再现观测降水和蒸发的年平均气候态空间分布及年循环特征,与观测值的空间相关系数分别为0.92和0.87。较之观测,MME高估了华北地区降水(0.55 mm d?1),低估了华南沿海地区降水(?0.3 mm d?1)。所有CMIP6模式均高估蒸发强度(偏差0.03~0.98 mm d?1),使得模拟的降水与蒸发之差偏少。模式整体能够模拟出我国东部季风区降水再循环率及不同边界水汽来源的贡献率,但低估了由南边界进入季风区的水汽贡献,导致东亚季风区偏干。通过分析模式对影响水汽通量的两个气象要素(风速和大气比湿)的模拟能力,发现研究区南边界的风速大小决定了模式间水汽输送差异。南边界风速越大的模式,由南边界进入的水汽通量越大,模式模拟的降水越多。西北太平洋辐合带的东西位置是影响南边界南风速的重要系统之一,辐合带位置偏东的模式模拟的南风强度较弱,使得水汽输送偏弱、降水偏少;反之,南边界水汽输送偏强、降水偏多。本文通过评估最新一代CMIP6模式在东亚水循环方面的模拟性能,指出了当前气候模式在模拟西太平洋辐合带位置方面存在的偏差及其对东亚水循环的影响。  相似文献   

3.
In this study, we thoroughly analyzed abrupt behaviors, trends, and periodicity properties of water vapor flux and moisture budget entering and exiting the four edges of the Pearl River basin based on the NCAR/NCEP reanalysis dataset by using the continuous wavelet transform and the simple two-phase linear regression technique. Possible implications for hydrological cycle and water resource management of these changes are also discussed. The results indicate that: (1) the water vapor propagating through the four edges of the Pearl River basin is decreasing, and it is particularly true for the changes of the water vapor flux exiting from the north edge of the study river basin. The transition point from increase to decrease occurs in the early 1960s; (2) The wavelet transform spectra indicate that the monthly water vapor flux through the north edge decreases and this decrease is mainly reflected by intermittent distribution of the wavelet power spectra after early 1980s. The periodicity properties of the water vapor flux through the north edge imply that the northward propagation of water vapor flux decreases after the 1980s; (3) close relations between water vapor flux, precipitation and streamflow implies that the altered hydrological cycle in the Pearl River basin is mainly manifested by seasonal shifts of water vapor flux after early 1960s. One of the direct consequences of these changes of water vapor flux is the seasonal transition of wet and dry conditions across the Pearl River basin. Regional responses of hydrological cycle to climate variation/change could be different from one river basin to another. Hydrological responses of the Pearl River basin to the global warming are mainly demonstrated by seasonal shifts of precipitation changes: winter comes to be wetter and summer tends to be dryer. The finding of the seasonal transition of precipitation in the Pearl River basin is of great scientific and practical merits in basin scale water resource management in the Pearl River basin under the changing climate and global warming in particular.  相似文献   

4.
The Tibetan Plateau (TP) surfaces have been experiencing an overall rapid warming and wetting while wind speed and solar radiation have been declining in the last three decades. This study investigated how climate changes influenced the hydrological cycle on the TP during 1984??2006. To facilitate the analysis, a land surface model was used to simulate surface water budget at all CMA (China Meteorological Administration) stations on the TP. The simulated results were first validated against observed ground temperature and observation-derived heat flux on the western TP and observed discharge trends on the eastern TP. The response of evaporation and runoff to the climate changes was then analyzed. Major finding are as follows. (1) Surface water balance has been changed in recent decades. Observed precipitation shows insignificant increasing trends in central TP and decreasing trends along the TP periphery while evaporation shows overall increasing trends, leading to decreased discharge at major TP water resource areas (semi-humid and humid zones in the eastern and southern TP). (2) At the annual scale, evaporation is water-limited in dry areas and energy-limited (radiation and air temperature) in wet areas; these constraints can be interpreted by the Budyko-curve. Evaporation in autumns and winters was strongly controlled by soil water storage in summers, weakening the dependence of evaporation on precipitation at seasonal scales. (3) There is a complementary effect between the simulated actual evaporation and potential evaporation, but this complementary relationship may deviate from Bouchet??s hypothesis when vapor pressure deficit (or air temperature) is too low, which suppresses the power of vapor transfer.  相似文献   

5.
华文剑  陈海山 《大气科学》2011,35(1):121-133
利用“国际耦合模式比较计划” (Phase 3 of the Coupled Model Intercomparison Project, CMIP3) 12个模式对20世纪 (The Twentieth-Century Climate in Coupled Models, 20C3M) 和21世纪SRES (Special Report on Emissions Scenarios) A1B 情景下的模拟结果, 通过21世纪 (2001~2099年) 与20世纪 (1901~1999年) 陆面能量和水文变量的对比分析, 揭示了陆面过程对全球变暖响应的基本特征, 并探讨了其可能的响应机制。结果表明, 与20世纪相比, 21世纪全球陆面平均的表面温度、 地表净辐射、 潜热通量明显增加; 而感热通量有所减小。降水、 径流、 蒸发等地表水循环分量也表现出不同程度的增加, 而土壤含水量有减小趋势。通过分析近地层主要大气强迫变量与陆面变量之间的联系, 发现陆面能量平衡过程对全球变暖的响应主要受向下长波辐射和气温变化的影响, 而温度的变化对陆面水文过程的影响起决定性的作用。进一步分析表明, 陆面过程对全球变暖的响应存在明显的区域性差异, 陆面温度和感热对全球变暖响应最显著的区域位于北半球中高纬, 而净辐射和潜热对全球变暖的响应在亚洲中部和非洲大陆最显著。相对于20世纪, 21世纪主要是长波辐射和温度对陆面能量平衡过程的贡献重要。对于陆面水文过程, 径流和土壤含水量对全球变暖的响应在亚洲中部以及北美最显著。在全球变暖背景下, 21世纪相对于20世纪, 温度对陆面水循环的影响更加显著, 主要体现在北半球中纬度地区。  相似文献   

6.
The atmospheric water holding capacity will increase with temperature according to Clausius-Clapeyron scaling and affects precipitation.The rates of change in future precipitation extremes are quantified with changes in surface air temperature.Precipitation extremes in China are determined for the 21st century in six simulations using a regional climate model,RegCM4,and 17 global climate models that participated in CMIP5.First,we assess the performance of the CMIP5 models and RCM runs in their simulation of extreme precipitation for the current period(RF:1982-2001).The CMIP5 models and RCM results can capture the spatial variations of precipitation extremes,as well as those based on observations:OBS and XPP.Precipitation extremes over four subregions in China are predicted to increase in the mid-future(MF:2039-58)and far-future(FF:2079-98)relative to those for the RF period based on both the CMIP5 ensemble mean and RCM ensemble mean.The secular trends in the extremes of the CMIP5 models are predicted to increase from 2008 to 2058,and the RCM results show higher interannual variability relative to that of the CMIP5 models.Then,we quantify the increasing rates of change in precipitation extremes in the MF and FF periods in the subregions of China with the changes in surface air temperature.Finally,based on the water vapor equation,changes in precipitation extremes in China for the MF and FF periods are found to correlate positively with changes in the atmospheric vertical wind multiplied by changes in surface specific humidity(significant at the p<0.1 level).  相似文献   

7.
环青海湖地区气候变化及其对湖泊水位的影响   总被引:32,自引:3,他引:32       下载免费PDF全文
分析青海湖地区1961~2000年气象观测资料得出:年和四季的气温、地表蒸发以及年和夏季、冬季降水变化的气候倾向率均为正值;而春季、秋季降水变化的气候倾向率为负值。气温升高、地表蒸发加大的趋势比较显著,而降水增多的趋势不显著且年代际变化比较大,气温、地表蒸发等气象要素有向暖干化过渡的趋势,这种暖干化趋势是造成青海湖水位下降的主要原因之一。  相似文献   

8.
Future change of global monsoon in the CMIP5   总被引:5,自引:1,他引:4  
This study investigates future changes of Global Monsoon (GM) under anthropogenic global warming using 20 coupled models that participated in the phase five of Coupled Model Intercomparison Project (CMIP5) by comparing two runs: the historical run for 1850–2005 and the Representative Concentration Pathway (RCP) 4.5 run for 2006–2100. A metrics for evaluation of models’ performance on GM is designed to document performance for 1980–2005 and best four models are selected. The four best models’ multi-model ensemble (B4MME) projects the following changes in the twenty-first century under the RCP4.5 scenario. (1) Monsoon domain will not change appreciably but land monsoon domain over Asia tends to expand westward by 10.6 %. (2) The annual mean and range of GM precipitation and the percentage of local summer rainfall will all amplify at a significant level over most of the global region, both over land and over ocean. (3) There will be a more prominent northern-southern hemispheric asymmetry and eastern-western hemispheric asymmetry. (4) Northern Hemisphere (NH) monsoon onset will be advanced and withdrawal will be delayed. (5) Changes in monsoon precipitation exhibits huge differences between the NH and the Southern hemisphere (SH). The NH monsoon precipitation will increase significantly due to increase in temperature difference between the NH and SH, significant enhancement of the Hadley circulation, and atmospheric moistening, against stabilization of troposphere. There is a slight decrease of the Walker circulation but not significant against the inter-model spread. There are important differences between the CMIP 3 and CMIP5 results which are discussed in detail.  相似文献   

9.
The impact of mid- and upper-level dry air, represented by low relative humidity (RH) values, on the genesis of tropical cyclone (TC) Durian (2001) in the South China Sea was investigated by a series of numerical experiments using the Weather Research and Forecasting model. The mid-level RH was lowered in different regions relative to TC Durian (2001)'s genesis location. Results suggest that the location of dry air was important to Durian (2001)'s genesis and intensification. The rapid development of the TC was accompanied by sustained near-saturated mid- and upper-level air, whereas low humidity decelerated its development. Water vapor budget analysis showed that moisture at mid and upper levels was mainly supplied by the vertical convergence of moisture flux and the divergence terms, and consumed by the condensation process. The horizontal convergence of moisture flux term supplied moisture in the air moistening process but consumed moisture in the air drying process. With a dryer mid- and upper-level environment, convective and stratiform precipitation were both inhibited. The upward mass fluxes and the diabatic heating rates associated with these two precipitation types were also suppressed. Generally, convection played the dominant role, since the impact of the stratiform process on vertical mass transportation and diabatic heating was much weaker. The vorticity budget showed that the negative vorticity convergence term, which was closely related to the inhibited convection, caused the vorticity to decrease above the lower troposphere in a dryer environment. The negative vorticity tendency is suggested to slow down the vertical coherence and the development rate of TCs.  相似文献   

10.
While most models project large increases in agricultural drought frequency and severity in the 21st century, significant uncertainties exist in these projections. Here, we compare the model-simulated changes with observation-based estimates since 1900 and examine model projections from both the Coupled Model Inter-comparison Project Phase 3 (CMIP3) and Phase 5 (CMIP5). We use the self-calibrated Palmer Drought Severity Index with the Penman-Monteith potential evapotranspiration (PET) (sc_PDSI_pm) as a measure of agricultural drought. Results show that estimated long-term changes in global and hemispheric drought areas from 1900 to 2014 are consistent with the CMIP3 and CMIP5 model-simulated response to historical greenhouse gases and other external forcing, with the short-term variations within the model spread of internal variability, despite that regional changes are still dominated by internal variability. Both the CMIP3 and CMIP5 models project continued increases (by 50–200 % in a relative sense) in the 21st century in global agricultural drought frequency and area even under low-moderate emissions scenarios, resulting from a decrease in the mean and flattening of the probability distribution functions (PDFs) of the sc_PDSI_pm. This flattening is especially pronounced over the Northern Hemisphere land, leading to increased drought frequency even over areas with increasing sc_PDSI_pm. Large differences exist in the CMIP3 and CMIP5 model-projected precipitation and drought changes over the Sahel and northern Australia due to uncertainties in simulating the African Inter-tropical convergence zone (ITCZ) and the subsidence zone over northern Australia, while the wetting trend over East Africa reflects a robust response of the Indian Ocean ITCZ seen in both the CMIP3 and CMIP5 models. While warming-induced PET increases over all latitudes and precipitation decreases over subtropical land are responsible for mean sc_PDSI_pm decreases, the exact cause of its PDF flattening needs further investigation.  相似文献   

11.
The climatological characteristics of precipitation and the water vapor budget in the Haihe River basin (HRB) are analyzed using daily observations at 740 stations in China in 1951-2007 and the 4-time daily ERA40 reanalysis data in 1958-2001. The results show that precipitation and surface air temperature present significant interannual and interdecadal variability, with cold and wet conditions before the 1970s but warm and dry conditions after the 1980s. Precipitation has reduced substantially since the 1990s, with a continued increase of surface air temperature. The total column water vapor has also reduced remarkably since the late 1970s. The multi-model ensemble from the Fourth Assessment Report (AR4) of the Intergovernmental Panel on Climate Change (IPCC) has capably simulated the 20th century climate features and successfully reproduced the spatial patterns of precipitation and temperature. Unfortunately, the models do not reproduce the interdecadal changes. Based on these results, future projections of the climate in the HRB are discussed under the IPCC Special Report on Emissions Scenarios (SRES) B1, A1B, and A2. The results show that precipitation is expected to increase in the 21st century, with substantial interannual fluctuations relative to the models’ baseline climatology. A weak increasing trend in precipitation is projected before the 2040s, followed by an abrupt increase after the 2040s, especially in winter. Precipitation is projected to increase by 10%-18% by the end of the 21st century. Due to the persistent warming of surface air temperature, water vapor content in the lower troposphere is projected to increase. Relative humidity will decrease in the mid-lower troposphere but increase in the upper troposphere. On the other hand, precipitation minus evaporation remains positive throughout the 21st century. Based on these projection results, the HRB region is expected to get wetter in the 21st century due to global warming.  相似文献   

12.
使用多种观测资料和43个参加耦合模式比较计划第五阶段(CMIP5)的全球气候模式模拟数据,评估分析了全球气候模式对中国地区1980-2005年降水特征的模拟能力。结果表明:多数CMIP5模式能够模拟出中国降水由西北向东南递增的分布特点,这与耦合模式比较计划第三阶段(CMIP3)的模式模拟结果类似,但华南地区降水模拟偏少,西部高原地区降水模拟偏多。模式能够较好地模拟出降水冬弱夏强的季节变化特征,但降水模拟系统性偏多。从EOF分析结果来看,多数CMIP5模式可以再现中国地区年平均降水的时空变化特征,集合平均的表现优于CMIP3。多模式集合在月、季、年时间尺度下模拟的平均值优于大部分单个模式的结果。CMIP5中6个中国模式的模拟能力与其他模式相当,其中FGOALS-g2、BCC-CSM1-1-m的模拟能力相对较好。  相似文献   

13.
Surface-level moisture transport over the Indian Ocean has been computed using NOAA/HIRS data for the years 1980, 1981 and 1984. The global relation between monthly mean surface-level humidity and precipitable water (Liu, 1986) has been applied for the computation of surface-level humidity using monthly mean satellite-derived water vapour. The monthly mean surface wind fields over the Indian ocean provided by Florida State University have been used for the surface-level moisture flux computations. Our analysis indicates net positive surface-level moisture flux divergence over the Arabian Sea and negative moisture flux divergence over the Bay of Bengal. It has also been found that evaporation over the Arabian Sea is a variable quantity and forms a significant part of the net moisture budget over the Arabian Sea. The relative contribution of cross-equatorial flux and evaporation from the Arabian Sea has been studied for all three years.  相似文献   

14.
The East Asian Monsoon Simulation with IAP AGCMs-A Composite StudyWangHuijunandBiXunqiang(InstituteofAtmosphericPhysics(IAP),...  相似文献   

15.
The water vapor budget and the cloud microphysical processes associated with a heavy rainfall system in the Dabie Mountain area in June 2008 were analyzed using mesoscale reanalysis data(grid resolution 0.03 × 0.03,22 vertical layers,1-h intervals),generated by amalgamating the local analysis and prediction system(LAPS).The contribution of each term in the water vapor budget formula to precipitation was evaluated.The characteristics of water vapor budget and water substances in various phase states were evaluated and their differences in heavy and weak rainfall areas were compared.The precipitation calculated from the total water vapor budget accounted for 77% of actual precipitation;surface evaporation is another important source of water vapor.Water vapor within the domain of interest mainly came from the lower level along the southern boundary and the lower-middle level along the western boundary.This altitude difference for water vapor flux was caused by different weather systems.The decrease of local water vapor in the middle-lower layer in the troposphere during the system development stage also contributed to precipitation.The strength and the layer thickness of water vapor convergence and the content of various water substances in the heavy rainfall areas were obviously larger than in the weak rainfall areas.The peak values of lower-level water vapor convergence,local water vapor income,and the concentration of cloud ice all preceded the heaviest surface rainfall by a few hours.  相似文献   

16.
CMIP5全球气候模式对青藏高原地区气候模拟能力评估   总被引:9,自引:4,他引:5  
胡芩  姜大膀  范广洲 《大气科学》2014,38(5):924-938
青藏高原是气候变化的敏感和脆弱区,全球气候模式对于这一地区气候态的模拟能力如何尚不清楚。为此,本文使用国际耦合模式比较计划第五阶段(CMIP5)的历史模拟试验数据,评估了44 个全球气候模式对1986~2005 年青藏高原地区地表气温和降水两个基本气象要素的模拟能力。结果表明,CMIP5 模式低估了青藏高原地区年和季节平均地表气温,年均平均偏低2.3℃,秋季和冬季冷偏差相对更大;模式可较好地模拟年和季节平均地表气温分布型,但模拟的空间变率总体偏大;地形效应校正能够有效订正地表气温结果。CMIP5 模式对青藏高原地区降水模拟能力较差。尽管它们能够模拟出年均降水自西北向东南渐增的分布型,但模拟的年和季节降水量普遍偏大,年均降水平均偏多1.3 mm d-1,这主要是源于春季和夏季降水被高估。同时,模式模拟的年和季节降水空间变率也普遍大于观测值,尤其表现在春季和冬季。相比较而言,44 个模式集合平均性能总体上要优于大多数单个模式;等权重集合平均方案要优于中位数平均;对择优挑选的模式进行集合平均能够提高总体的模拟能力,其中对降水模拟的改进更为显著。  相似文献   

17.
Various ocean reanalysis data reveal that the subarctic Atlantic sea surface temperature (SST) has been cooling during the twentieth century. A similar cooling pattern is found in the doubling CO2 experiment obtained from the CMIP3 (coupled model intercomparison project third phase) compared to the pre-industrial experiment. Here, in order to investigate the main driver of this cooling, we perform the heat budget analysis on the subarctic Atlantic upper ocean temperature. The net surface heat flux associated with the increased concentration of greenhouse gases heats the subarctic ocean surface. In the most of models, the longwave radiation, latent heat flux, and sensible heat flux exert a warming effect, and the shortwave radiation exerts a cooling effect. On the other hand, the thermal advection by the meridional current reduces the subarctic upper ocean temperature in all models. This cold advection is attributed to the weakening of the meridional overturning circulation, which is related to the reduction in the ocean surface density. In particular, greater warming of the surface air than of the sea surface results in the reduction of surface evaporation and thereby enhanced freshening of the ocean surface water, while precipitation change was smaller than evaporation change. The thermal advections by both the wind-driven Ekman current and the density-driven geostrophic current contribute to cooling in most of the models, where the heat transport by the geostrophic current tends to be larger than that by the Ekman current.  相似文献   

18.
The temperature biases of 28 CMIP5 AGCMs are evaluated over the Tibetan Plateau(TP) for the period 1979–2005. The results demonstrate that the majority of CMIP5 models underestimate annual and seasonal mean surface 2-m air temperatures(T_(as)) over the TP. In addition, the ensemble of the 28 AGCMs and half of the individual models underestimate annual mean skin temperatures(T_s) over the TP. The cold biases are larger in T_(as) than in T_s, and are larger over the western TP. By decomposing the T_s bias using the surface energy budget equation, we investigate the contributions to the cold surface temperature bias on the TP from various factors, including the surface albedo-induced bias, surface cloud radiative forcing, clear-sky shortwave radiation, clear-sky downward longwave radiation, surface sensible heat flux, latent heat flux,and heat storage. The results show a suite of physically interlinked processes contributing to the cold surface temperature bias.Strong negative surface albedo-induced bias associated with excessive snow cover and the surface heat fluxes are highly anticorrelated, and the cancelling out of these two terms leads to a relatively weak contribution to the cold bias. Smaller surface turbulent fluxes lead to colder lower-tropospheric temperature and lower water vapor content, which in turn cause negative clear-sky downward longwave radiation and cold bias. The results suggest that improvements in the parameterization of the area of snow cover, as well as the boundary layer, and hence surface turbulent fluxes, may help to reduce the cold bias over the TP in the models.  相似文献   

19.
林志强  乔钰  唐叔乙 《气象科技》2014,42(1):110-117
利用西藏高原39个气象站自1961年或建站开始至2010年汛期(5—9月)的逐月降水资料,根据Z指数和区域旱涝指数分析了近50年西藏高原汛期旱涝的变化特征。在此基础上,选取了典型旱、涝年,采用NCEP/NCAR再分析月平均资料,通过合成分析方法分析了旱、涝年的水汽输送特征的差异。结果表明:近50年以来,西藏高原以旱为主,涝灾较少,干湿状况具有明显年代际变化;西藏高原汛期旱涝与水汽输送异常联系紧密,典型旱(涝)年北印度洋、西北太平洋和亚洲中纬度地区的整层水汽通量经向、纬向及整层水汽通量和水汽通量散度均存在显著差异,呈反位相分布。越赤道急流区的差异最为显著,索马里越赤道气流区的水汽输送差异主要表现在经向上,而菲律宾越赤道气流区的水汽输送差异主要表现在纬向上。  相似文献   

20.
一次西南涡特大暴雨的中尺度诊断分析   总被引:1,自引:0,他引:1       下载免费PDF全文
采用LAPS中尺度分析模式大气资料,对2008年7月一次西南涡暴雨过程进行天气学降水运动的中尺度诊断计算与分析。诊断计算包括:可降水量、层结不稳定能量、对流可降水量、水汽权重平均风速、水汽通量散度、云水、云冰总量及其通量散度和垂直速度与凝结函数降水率等。结果表明:“西南涡-切变线”系统的暴雨发生在暖湿气团与变性冷气团之间的中尺度风场辐合上升运动区,中尺度雨团发生在层结不稳定的暖湿气团一侧。计算的中尺度垂直运动与凝结函数降水率场,降水率为暴雨到特大暴雨。计算的水汽通量辐合降水率与凝结函数降水率不会完全重合,且水汽通量辐合既可致中尺度“雨”,又可成大尺度“云”,并且云水、云冰通量辐合/辐散,可解释为它们的“正”/“负”碰并增长,而碰并增长产生水凝物增量(降水率)也促成大暴雨。因此,在凝结函数降水率场中产生的中、小尺度对流雨团,加上水汽与云水、云冰通量辐合及其碰并增长,并且借助层结不稳定能量释放和可能产生的强迫“次级环流”及水汽与云水、云冰输送,是这次“西南涡-切变线”系统造成襄樊特大暴雨的天气学成因。  相似文献   

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

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