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

2.
基于1961~2013年观测资料和第六次国际耦合模式比较计划(CMIP6)中32个气候模式的最低气温数据,评估了各模式及多模式集合平均方法对中国冬季寒潮频次的模拟能力,并筛选了中国不同区域的最优模式,为中国未来寒潮频次预估和气候模式改进提供理论支撑。结果表明,CMIP6全球气候模式可以很好地再现出中国冬季寒潮频次由北向南逐级递减的空间特征,EC-Earth3-Veg对中国寒潮频次模拟能力最好;绝大多数模式可以模拟出中国冬季寒潮频次下降的趋势,但对趋势变化幅度的模拟能力则相对有限。多模式等权重集合平均和多模式中位数集合平均较单一模式分别对中国北方地区和南方地区冬季寒潮频次空间格局的模拟效果改进显著,但是对趋势值的模拟普遍偏低。  相似文献   

3.
青藏高原中东部夏季降水主要表现为东北和东南反位相变化的双极型特征。采用经验正交函数(empirical orthogonal function,EOF)分解方法,系统性地评估参与第五次耦合模式比较计划 (Coupled Model Intercomparison Project Phase 5,CMIP5)历史模拟试验的 47 个模式对青藏高原中东部夏季降水双极型变化特征的模拟能力。结果表明,大多数模式基本可以反映青藏高原中东部夏季降水东北部和东南部反位相的变化特征。模式间 EOF 分析结果表明在35°N 以南的东西向模拟偏差是 CMIP5 模式模拟降水空间型态的主要偏差,且大多数模式对时间系数的模拟效果差于空间型态。文中定义了一个综合评估指标 Snew 来定量描述模式对空间型态、时间系数以及方差贡献的综合模拟效果。由定量评估结果来看,MIROC-ESM、HadGEM2-CC 和 ACCESS1-0 (FIO-ESM、 HadGEM2-AO 和 MIROC-ESM-CHEM)模式对观测降水的 EOF1(EOF2)模态的综合模拟能力相对较好,而 GISS 系列模式、CESM1-CAM5 和 MPI-ESM-LR (CMCC-CESM、MPI-ESM-MR 和 GFDL- CM3)模式对观测降水的 EOF1(EOF2)模态的综合模拟效果较差。由 EOF1 和 EOF2 的综合评估结果来看,MIROC-ESM-CHEM模式对观测降水的 EOF1 和 EOF2 模态的综合模拟效果最好。  相似文献   

4.
中国冬季多种积雪参数的时空特征及差异性   总被引:6,自引:2,他引:4  
利用1979~2006年冬季中国站点最大雪深和站点雪日、卫星遥感雪深、积雪覆盖率和雪水当量5种积雪资料,从多角度深入细致地分析了我国冬季积雪的时空变化特征。结果表明:5种积雪资料的经验正交分解第一模态都表现为中国南、北方反位相的特征,即当新疆和东北三省-内蒙古地区积雪偏多(少)时,青藏高原和南方地区积雪偏少(多)。新疆和东北三省-内蒙古地区的雪深、积雪覆盖率和雪日随时间有逐渐增多的趋势,而其中边缘山区的雪水当量表现出减少的趋势,青藏高原地区的积雪表现出与其完全相反的特征。南方地区站点最大雪深和雪日表现出随时间减少的趋势,卫星遥感难以监测到该区积雪。相比较而言,卫星遥感资料比较适合高原和山区缺少气象站的地区及北半球更大区域积雪的研究,而站点资料更适用于中国中东部和平原地区积雪的区域研究。雪深、雪日、积雪覆盖率和雪水当量这些多样性积雪参数存在一定的差异性,因此5种积雪资料结合使用才能得到更准确的结论。  相似文献   

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

6.
基于卫星观测数据,评估了23个CMIP5耦合模式对北半球3—4月积雪面积的模拟能力,在此基础上应用多模式集合平均结果,预估了未来不同温室气体排放情景下北半球3—4月积雪面积的变化情况。结果表明:整体上看,CMIP5耦合模式对北半球3—4月积雪面积具有一定的模拟能力,模式基本能再现北半球3—4月积雪面积的分布特征,但对高原等复杂地形地区积雪的模拟偏差较大并且低估了北半球积雪的减少趋势,这些可能是由卫星资料本身的缺陷以及模式参数化方案的不同造成的。多模式集合预估结果表明,未来几十年北半球3—4月积雪将继续减少并且集中发生在欧亚大陆中西部地区。温室气体排放将会对未来北半球积雪的变化产生显著影响。在RCP8.5情景下,未来北半球积雪减少最显著;在RCP4.5和RCP6.0情景下,在21世纪前半叶北半球积雪减少趋势与RCP8.5情景相当,但是在21世纪后半叶积雪的减少趋势明显小于RCP8.5情景;在RCP2.6情景下,北半球积雪减少趋势最小。所以,控制温室气体排放对于未来北半球积雪的生存至关重要。  相似文献   

7.
潘延  张洋  李舒婷 《气象科学》2022,42(4):440-456
本文评估了36个CMIP5模式和39个CMIP6模式对近期观测中揭示的北半球冬季大气环流与高原冬春气温之间的相关关系的模拟能力。利用最大协方差(MCA)分析方法,计算并比较了观测和模式中冬季北半球200 hPa位势高度场与同后期青藏高原近地面气温的耦合关系。整体而言,大部分CMIP模式能够模拟出显著的冬季北半球大气环流与青藏高原气温之间的相关关系,且CMIP6模式模拟相关特征和作用机制的能力较CMIP5均有所提升。与观测相比,历史情景下36个CMIP5模式中有26个能够模拟出显著的大气环流与同后期高原气温之间的相关关系,其中对于相关的位势高度场空间模态的模拟明显好于对高原气温异常场空间模态的模拟。同情景下39个CMIP6模式中有37个能模拟出显著相关关系,且CMIP6模式更能模拟出观测中MCA模态的位势高度场上北极涛动(AO)和西太平洋遥相关型(WP)反相位叠加的大气环流特征。在对MCA模态时间变率的模拟上,大部分模式都能重现青藏高原整体变暖的趋势,部分模式能够模拟出观测中位势高度场时间主成分的年际变率,并且CMIP6表现要优于CMIP5。对耦合环流型的动力诊断显示,相比CMIP5模式,CMIP6中有更多模式可以模拟出极地—高原之间的遥相关波列,且对波列结构的模拟更完整。  相似文献   

8.
CMIP5模式对冬季北极涛动的模拟和预估   总被引:1,自引:0,他引:1  
基于NCEP/NCAR再分析资料和CMIP5的19个模式结果,从异常模态、年代际趋势和周期特征等方面评估了CMIP5耦合模式对冬季北极涛动(Arctic Oscillation,AO)的模拟能力,并对未来RCP4.5、RCP8.5两种浓度路径下AO的可能变化趋势给出了定性的预估。CMIP5模式历史试验结果显示,大多数模式都能够模拟出AO模态的基本结构,但是对中心位置、强度的模拟存在较大的偏差,其中MPI-ESM-LR和Had GEM2-AO能较好地模拟出AO整体模态来。在历史演变和周期特征的刻画方面,模式的冬季海平面气压经验正交函数分解第一模态时间序列(Principal Component,PC1)基本能够反映出1950~1970年以来的减弱趋势,但对1970年以后的增长趋势模拟并不明显,而北半球环状模指数(Zonal Index,ZI)序列对两个阶段的趋势均可模拟出来,模式的PC1和ZI序列总体表现为正的变化趋势。有一半以上的模式对2~3 a高频周期模拟较好,但对20 a左右的周期模拟较差,其中仅有Can ESM2、CNRM-CM5、GFDL-ESM2G这3个模式对ZI指数的两个周期变化模拟较好。在RCP4.5和RCP8.5两种浓度路径下,ZI序列有显著的上升趋势,从长期趋势系数看RCP4.5路径下有14个模式呈现正的变化趋势,其中有10个模式通过了检验。RCP8.5浓度路径下,16个模式为正变化趋势,有11个模式通过了检验,集合平均结果正变化趋势较为显著。两种浓度路径下不同时段的海平面气压变化趋势表明,ZI序列的年代际变化明显,存在3个不同的变化阶段——2006~2039年、2070~2100年为两个上升阶段,2040~2069年为缓慢下降阶段。  相似文献   

9.
青藏高原东坡陡峭地形区是气候模式陆地降水模拟偏差的大值区,且这一偏差长期未得到有效改善.基于17个参加国际耦合模式比较计划第六阶段(CMIP6)的全球气候模式的日降水结果,评估了当前最新一代的气候模式对青藏高原东坡地区2000—2014年暖季(5—9月)降水气候态及其季节内演变的模拟能力.结果表明:高原东坡降水正偏差存...  相似文献   

10.
利用CMIP5中20个模式的历史模拟结果和英国东英格利亚大学CRU观测数据,采用泰勒图、趋势分析、滑动平均和EOF等方法从气候态和气候变率两方面检验各个模式对中亚年平均气温的模拟能力。结果表明:各模式能较好地模拟1951—2005年中亚地区显著增温趋势和年平均气温的空间分布特点,尤其是高、低值中心和等值线数值分布。泰勒图分析显示,大部分模式的均方根误差在0.5左右,空间相关系数在0.85~0.90之间,标准差在0.5~1.0之间。EOF分析结果表明,模式集合平均能够较好地表现出中亚气温一致升高以及南北反位相波动这两个主要模态的时空变化特征。  相似文献   

11.
Based on the simulations of 31 global models in CMIP5, the performance of the models in simulating the Hadley and Walker circulations is evaluated. In addition, their change in intensity by the end of the 21 st century(2080–2099) under the RCP4.5 and RCP8.5 scenarios, relative to 1986–2005, is analyzed from the perspective of 200 h Pa velocity potential.Validation shows good performance of the individual CMIP5 models and the multi-model ensemble mean(MME) in reproducing the meridional(zonal) structure and magnitude of Hadley(Walker) circulation. The MME can also capture the observed strengthening tendency of the winter Hadley circulation and weakening tendency of the Walker circulation. Such secular trends can be simulated by 39% and 74% of the models, respectively. The MME projection indicates that the winter Hadley circulation and the Walker circulation will weaken under both scenarios by the end of the 21 st century. The weakening amplitude is larger under RCP8.5 than RCP4.5, due to stronger external forcing. The majority of the CMIP5 models show the same projection as the MME. However, for the summer Hadley circulation, the MME shows little change under RCP4.5 and large intermodel spread is apparent. Around half of the models project an increase, and the other half project a decrease. Under the RCP8.5 scenario, the MME and 65% of the models project a weakening of the summer southern Hadley circulation.  相似文献   

12.
本文基于第五次耦合模式比较计划的23个全球气候模式所提供的最高气温与最低气温在RCP4.5情景下的逐日格点资料,根据模式对5个极端气温指数的模拟能力,使用秩加权方法研究了中国未来极端气温变化的概率预估及其不确定性。结果表明,21世纪中期(2046—2065年)中国区域平均最高气温和平均最低气温的增加幅度相对于历史时期(1986—2005年)可能超过2.0℃(概率>66%),增加的大值区主要位于青藏高原南部。暖夜指数在中国大部分地区增加超过15%,西南和东南部沿海是增加的大值区,增幅超过20%。霜冻日数在全国范围内减少,减少的大值区位于青藏高原周围,减少日数超过了20 d。热浪指数在整个中国区域可能增加10 d以上,大值区位于西藏西南部,可达30 d。不确定性的结果表明,除热浪指数的可信度较低外,其余指数都有较高的可信度。到21世纪末期(2081—2100年),中国区域极端气温增加幅度超过前期,平均最高气温和平均最低气温很可能增加超过2.0℃(概率>90%),大值区除中国西部地区外,还扩展到了东北和青藏高原西南地区。中国大部分地区的暖夜指数增加超过15%,西南和南部沿海可能超过25%。大部分地区的霜冻日数减少20 d,青藏高原周围减少超过40 d。热浪指数在中国范围内增加20 d,青藏高原西南部增加40 d以上。除霜冻指数的信噪比略比21世纪中期大外,其余指数的信噪比与中期基本一致。  相似文献   

13.
The trends and fluctuations of observed and CMIP5-simulated yearly mean surface air temperature over China were analyzed.In general,the historical simulations replicate the observed increase of temperature,but the multi-model ensemble (MME) mean does not accurately reproduce the drastic interannual fluctuations.The correlation coefficient of the MME mean with the observations over all runs and all models was 0.77,which was larger than the largest value (0.65) from any single model ensemble.The results showed that winter temperatures are increasing at a higher rate than summer temperatures,and that winter temperatures exhibit stronger interannual variations.It was also found that the models underestimate the differences between winter and summer rates.The ensemble empirical mode decomposition technique was used to obtain six intrinsic mode functions (IMFs) for the modeled temperature and observations.The periods of the first two IMFs of the MME mean were 3.2 and 7.2,which represented the cycle of 2-7-yr oscillations.The periods of the third and fourth IMFs were 14.7 and 35.2,which reflected a multi-decadal oscillation of climate change.The corresponding periods of the first four IMFs were 2.69,7.24,16.15 and 52.5 in the observed data.The models overestimate the period of low frequency oscillation of temperature,but underestimate the period of high frequency variation.The warming rates from different representative concentration pathways (RCPs) were calculated,and the results showed that the temperature will increase by approximately 0.9℃,2.4℃,3.2℃ and 6.1℃ in the next century under the RCP2.6,RCP4.5,RCP6.0 and RCP8.5 scenarios,respectively.  相似文献   

14.
Based on integrated simulations of 26 global climate models provided by the Coupled Model Intercomparison Project(CMIP), this study predicts changes in temperature and precipitation across China in the 21 st century under different representative concentration pathways(RCPs), and analyzes uncertainties of the predictions using Taylor diagrams. Results show that increases of average annual temperature in China using three RCPs(RCP2.6, RCP4.5,RCP8.5) are 1.87 ℃, 2.88 ℃ and 5.51 ℃, respectively. Increases in average annual precipitation are 0.124, 0.214, and 0.323 mm/day, respectively. The increased temperature and precipitation in the 21 st century are mainly contributed by the Tibetan Plateau and Northeast China. Uncertainty analysis shows that most CMIP5 models could predict temperature well, but had a relatively large deviation in predicting precipitation in China in the 21 st century. Deviation analysis shows that more than 80% of the area of China had stronger signals than noise for temperature prediction;however, the area proportion that had meaningful signals for precipitation prediction was less than 20%. Thus, the multi-model ensemble was more reliable in predicting temperature than precipitation because of large uncertainties of precipitation.  相似文献   

15.
中国降水季节性的预估   总被引:2,自引:1,他引:1  
姚世博  姜大膀  范广洲 《大气科学》2018,42(6):1378-1392
本文使用国际耦合模式比较计划第五阶段(CMIP5)中共46个全球气候模式的数值试验数据,通过评估择优选取了14个模式来预估21世纪中国各季节降水百分率及其变率。结果表明,模式集合平均能够较好地模拟各季节降水百分率及其变率,但模式与观测间、各模式间都存在一定不同,空间上西部差异较大,季节上夏季差异明显。21世纪中国降水百分率整体表现为夏季大冬季小,但存在区域性,如华南春季降水百分率大于夏季。与1986~2004年相比,中国降水百分率整体表现为在夏季显著减少,冬春季显著增加,但高原则与之相反。此外,模式对于长江中下游地区降水百分率的预估存在较大不确定性。RCP8.5情景下降水季节性变幅要大于RCP4.5情景。降水季节性的变率在四季均表现出一定的增加趋势,但21世纪早、中和末期与1986~2004年相比并无显著差异(置信水平为95%)。  相似文献   

16.
Climate changes in future 21 st century China and their uncertainties are evaluated based on 22 climate models from the Coupled Model Intercomparison Project Phase 5(CMIP5). By 2081–2100, the annual mean surface air temperature(SAT) is predicted to increase by 1.3℃± 0.7℃, 2.6℃± 0.8℃ and 5.2℃± 1.2℃ under the Representative Concentration Pathway(RCP) scenarios RCP2.6, RCP4.5 and RCP8.5, relative to 1986–2005, respectively. The future change in SAT averaged over China increases the most in autumn/winter and the least in spring, while the uncertainty shows little seasonal variation.Spatially, the annual and seasonal mean SAT both show a homogeneous warming pattern across China, with a warming rate increasing from southeastern China to the Tibetan Plateau and northern China, invariant with time and emissions scenario.The associated uncertainty in SAT decreases from northern to southern China. Meanwhile, by 2081–2100, the annual mean precipitation increases by 5% ± 5%, 8% ± 6% and 12% ± 8% under RCP2.6, RCP4.5 and RCP8.5, respectively. The national average precipitation anomaly percentage, largest in spring and smallest in winter, and its uncertainty, largest in winter and smallest in autumn, show visible seasonal variations. Although at a low confidence level, a homogeneous wetting pattern is projected across China on the annual mean scale, with a larger increasing percentage in northern China and a weak drying in southern China in the early 21 st century. The associated uncertainty is also generally larger in northern China and smaller in southwestern China. In addition, both SAT and precipitation usually show larger seasonal variability on the sub-regional scale compared with the national average.  相似文献   

17.
利用4个海气耦合模式对1960~2005年的多年代际回报结果,评估模式对中国区域年代际气候变化(温度和降水)的预测潜力,并初步给出2005~2015年的气候预测结果。与CMIP3多模式集合1960~2000年结果以及观测实况比较的结果表明:融入观测资料进行同化的年代际气候预测模式,对中国区域温度和降水的模拟能力总体好于CMIP3模式。年代际气候预测模式对温度气候场的模拟仍以"冷偏差"为主,但较之CMIP3模式已有显著改进,中国区域平均的冷偏差减少1.3°C;对降水气候场的模拟仍以"湿偏差"为主,但在华南沿海和西北内陆降水的模拟能力优于CMIP3模式,中国区域平均的湿偏差降低了20%。年代际模式和CMIP3模式都能较好地模拟出中国区域尤其是北方20世纪后期的增暖信号;但CMIP3模式对20世纪后期中国东部降水的旱涝结构演变的模拟与观测相反;而年代际气候预测模式未能再现华北偏旱的变化,但能成功地模拟江淮流域和华南沿海的旱涝演变。2005~2015年的10年预测表明中国区域将继续增暖0.3~0.7°C,且增温幅度北方大于南方,增幅中心位于西北内陆和青藏高原;而降水的变化趋势不显著,黄淮地区、西北内陆和青藏高原的降水略有增加,而西南地区降水将减少。但需要指出的是,这种预测的不确定性是相当大的。  相似文献   

18.
基于国际耦合模式比较计划第五阶段(CMIP5)历史模拟试验(historical run)的模式输出结果以及遥感数据,采用相关分析、均方根误差、标准差等统计方法,评估了13个气候(或地球)系统模式对欧亚大陆积雪覆盖率的模拟能力,在此基础上,采用多模式集合平均的方法对未来不同温室气体排放情景下(rcp2.6、rcp4.5和rcp8.5)欧亚大陆积雪覆盖率的变化进行预估。结果显示:尽管各模式模拟的积雪覆盖率在高原地区与观测差异较大,但总体看来模式能够对欧亚大陆积雪覆盖率的空间形态、季节变化及年际变化特征做出较好地模拟。未来预估结果表明,多模式集合平均预估的欧亚大陆积雪覆盖率从2006年到2040年左右减少趋势非常明显,且不同排放情景下模式模拟的积雪减少速率非常接近;然而,大约从2040年之后,不同排放情景下的积雪覆盖率减小趋势的差异越来越大,rcp2.6和rcp4.5下积雪覆盖率的变化趋于平缓,而rcp8.5情景下,积雪覆盖率一直减少,冬季、春季和秋季都明显减少,减少最显著的区域位于西欧和青藏高原地区。由此可见,控制温室气体的排放对于未来欧亚大陆积雪的变化是至关重要的。  相似文献   

19.
Snow cover changes in the middle (2040–2059) and end (2080–2099) of the twenty-first century over China were investigated with a regional climate model, nested within the global model BCC_CSM1.1. The simulations had been conducted for the period of 1950–2099 under the RCP4.5 and RCP8.5 scenarios. Results show that the model perform well in representing contemporary (1986–2005) spatial distributions of snow cover days (SCDs) and snow water equivalent (SWE). However, some differences between observation and simulation were detected. Under the RCP4.5 scenarios, SCDs are shortened by 10–20 and 20–40 days during the middle and end of the twenty-first century, respectively. Whereas simulated SWE is lowered by 0.1–10 mm in most areas over the Tibetan Plateau (TP). On the other hand, the spatial distributions of SWE are reversed between the middle and end terms in the northeast China. Furthermore, compared with the changes of RCP4.5 scenario, SCDs are reduced by 5–20 days in the middle period under RCP8.5 scenario with even larger decreasing amplitude in the end term. SWE was lowered by 0.1–2.5 mm in most areas except the northeast of China in middle term under RCP8.5 scenario. The great center of SCDs and SWE changes are always located over TP. The regional mean of SCDs and SWE for the TP and for China display a declining trend from 2006 to 2099 with more pronounced changes in the TP than in China as a whole. Under the RCP8.5 scenario, the changes are enhanced compared to those under RCP4.5.  相似文献   

20.
利用国际耦合模式比较计划第5阶段(CMIP5)中的21个气候模式的RCP4.5和RCP8.5情景预估结果,分析了全球变暖1.5℃和2℃阈值时青藏高原气温年和季节的变化特征。结果表明,对应1.5℃和2℃全球变暖,青藏高原变暖幅度明显更大,就整体而言,在RCP4.5/RCP8.5情景下,高原区域平均的平均、最高、最低气温变暖分别为2.11℃/2.10℃和2.96℃/2.85℃、2.02℃/2.02℃和2.89℃/2.77℃、2.34℃/2.34℃和3.20℃/3.14℃,冬季平均气温的变暖幅度(2.19℃/2.31℃和3.13℃/3.05℃)较其他季节更大;从空间分布形势上看,年变暖呈西南高东北低的分布,而春、冬变暖呈南高北低的分布,夏、秋变暖则呈西高东低的分布。到达同一温升阈值时,RCP4.5与RCP8.5情景下高原气温的响应也存在区域差异。高原年与各季平均气温对全球变暖1.5℃与2℃的响应差异均>0.5℃,其中冬季最明显,区域平均差异可达0.94℃,局地差异超过1.1℃。  相似文献   

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

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