首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
中高纬度大气低频模态研究进展   总被引:1,自引:0,他引:1  
大气低频模态是导致中高纬地区季节及气候变化的主要因素,加强低频模态动力机制的认识对于提高短期气候预测水平具有十分重要的科研意义和应用价值。由于中高纬大气低频环流本身的复杂性,关于其动力机制的研究仍然是大气动力学中重点难点问题。本文简单回顾了激发维持低频模态的前两种机制,即:外源强迫、气流的纬向非均匀性对大气低频模态的影响。从观测事实和模式结果出发,着重介绍了瞬变波与基本气流相互作用激发维持低频模态的第三种机制,且详细阐述了线性假设条件下,瞬变波与基本气流相互作用的动力机制。本文最后针对瞬变波与基本气流相互作用研究中一些尚未解决的问题进行了讨论。  相似文献   

2.
Interaction between synoptic eddy and low-frequency flow (SELF) has been the subject of many studies. In this study, we further examine the interaction by introducing a transformed eddy-potential-vorticity (TEPV) flux that is obtained from eddy-potential-vorticity flux through a quasi-geostrophic potential-vorticity inversion. The main advantage of using the TEPV flux is that it combines the effects of the eddy-vorticity and heat fluxes into the net acceleration of the low-frequency flow in such a way that the TEPV flux tends to be analogous to the eddy-vorticity fluxes in the barotropic framework. We show that the anomalous TEPV fluxes are preferentially directed to the left-hand side of the low-frequency flow in all vertical levels throughout the troposphere for monthly flow anomalies and for climate modes such as the Arctic Oscillation (AO). Furthermore, this left-hand preference of the TEPV flux direction is a convenient three-dimensional indicator of the positive reinforcement of the low-frequency flow by net eddy-induced acceleration. By projecting the eddy-induced net accelerations onto the low-frequency flow anomalies, we estimate the eddy-induced growth rates for the low frequency flow anomalies. This positive eddy-induced growth rate is larger (smaller) in the lower (upper) troposphere. The stronger positive eddy feedback in the lower troposphere may play an important role in maintaining an equivalent barotropic structure of the low-frequency atmospheric flow by balancing some of the strong damping effect of surface friction.  相似文献   

3.
不同水平分辨率BCC_CSM模式对中亚地面气温模拟能力评估   总被引:1,自引:1,他引:1  
本研究基于IPCC AR5(Intergovernmental Panel on Climate Change:Fifth Assessment Report)中BCC_CSM1.1(Beijing Climate Center Climate System Model version 1.1)和BCC_CSM1.1(m)(Beijing Climate Center Climate System Model version 1.1 with a Moderate Resolution)气候模式的历史试验结果和CRU(Climatic Research Unit)资料, 采用趋势分析和滑动平均等方法检验了两个版本BCC_CSM模式对中亚地区1948~2011年平均地表气温、各热通量及其趋势的模拟性能, 并讨论了不同模式水平分辨率的影响, 结果显示:BCC_CSM1.1和BCC_CSM1.1(m)两个模式均能够模拟出中亚地区显著增温以及感热通量、长/短波净辐射等要素由南向北递减的总体趋势。其中, BCC_CSM1.1(m) 在对中亚地面年平均气温、感热通量和长/短波净辐射空间分布的模拟结果好于BCC_CSM1.1, 但对于气温标准差的模拟, BCC-CSM1.1模式略好于BCC-CSM1.1(m)。模式分辨率的提高, 能够更好地表现出地形的影响, 对气温和各热通量模拟性能改善较大, 在中亚地区年平均气温的模拟中表现出了一定的优势。  相似文献   

4.
BCC_CSM模式对热带降水年循环模态的模拟   总被引:5,自引:2,他引:5  
本文评估了国家气候中心发展的两个不同分辨率海—陆—气—冰多圈层耦合气候系统模式BCC_CSM1.1和BCC_CSM1.1 (m) 对热带降水两个年循环模态——揭示降水冬夏季节差异的季风模态和揭示过渡季节春季和秋季非对称特征的春秋非对称模态的模拟能力,讨论了模拟偏差产生的可能原因。分析结果表明,BCC_CSM1.1和BCC_CSM1.1 (m) 均能合理再现全球年平均降水的基本分布特征,也能较合理再现热带降水年循环模态的基本分布特征,尤其季风模态中降水与环流关于赤道反对称的特征;能够较合理再现春秋非对称模态与热带海洋表面温度(SST)年循环之间的关系。大气温度场、环流场以及热带SST的模拟偏差对降水季风模态有影响;热带SST年循环的偏差对降水春秋非对称模态的模拟偏差有贡献;模式分辨率对降水年循环模态的模拟也有一定影响。对比分析显示,大气模式和陆面模式水平分辨率提高之后模式在某些模拟性能上有所提高,这表现在:BCC_CSM1.1 (m) 模拟的1~12月降水气候态的空间变率更接近观测;热带海表温度年循环总体上更接近观测;模拟的热带降水年循环模态的部分特征更合理。但BCC_CSM1.1 (m) 的模拟结果相对观测仍存在较大偏差,有待进一步改进。  相似文献   

5.
较全面地介绍了北京气候中心气候系统模式(BCC_CSM)研发所取得的一些进展及其在气候变化研究中的应用,重点介绍了全球近280 km较低分辨率的全球海-陆-气-冰-生物多圈层耦合的气候系统模式BCC_CSM1.1和110 km中等大气分辨率的BCC_CSM1.1(m),以及大气、陆面、海洋、海冰各分量模式的发展。BCC_CSM1.1和BCC_CSM1.1(m)气候系统模式均包含了全球碳循环和动态植被过程。当给定全球人类活动导致的碳源排放后,就可以模拟和预估人类活动对气候变化的影响。BCC_CSM1.1和BCC_CSM1.1(m)已应用于IPCC AR5模式比较,为中外开展气候变化机理分析和未来气候变化预估提供了大量的试验数据。还介绍了BCC_CSM1.1和BCC_CSM1.1(m)参与国际耦合模式比较计划(CMIP5)的大量试验分析评估结果,BCC_CSM能够较好地模拟20世纪气温和降水等气候平均态和季节变化特征,以及近1000年的历史气候变化,所预估的未来100年气候变化与国际上其他模式的CMIP5试验预估结果相当。初步的分析表明,分辨率相对高的BCC_CSM1.1(m)在区域气候平均态的模拟上优于分辨率较低的BCC_CSM1.1。  相似文献   

6.
Long-lead precipitation forecasts for 1–4 seasons ahead are usually difficult in dynamical climate models due to the model deficiencies and the limited persistence of initial signals. But, these forecasts could be empirically improved by statistical approaches. In this study, to improve the seasonal precipitation forecast over the southern China (SC), the statistical downscaling (SD) models are built by using the predictors of atmospheric circulation and sea surface temperature (SST) simulated by the Beijing Climate Center Climate System Model version 1.1 m (BCC_CSM1.1 m). The different predictors involved in each SD model is selected based on both its close relationship with the target seasonal precipitation and its reasonable prediction skill in the BCC_CSM1.1 m. Cross and independent validations show the superior performance of the SD models, relative to the BCC_CSM1.1 m. The temporal correlation coefficient of SD models could reach > 0.4, exceeding the 95 % confidence level. The SC precipitation index can be much better forecasted by the SD models than by the BCC_CSM1.1 m in terms of the interannual variability. In addition, the errors of the precipitation forecast in all four seasons are significantly reduced over most of SC in the SD models. For the 2015/2016 strong El Niño event, the SD models outperform the dynamical BCC_CSM1.1 m model on the spatial and regional-average precipitation anomalies, mostly due to the effective SST predictor in the SD models and the weak response of the SC precipitation to El Niño-related SST anomalies in the BCC_CSM1.1 m.  相似文献   

7.
The sea surface temperature anomalies(SSTAs) in the tropical Indian Ocean(TIO) show two dominant modes at interannual time scales,referred to as the Indian Ocean basin mode(IOBM) and dipole mode(IOD).Recent studies have shown that the IOBM and IOD not only affect the local climate,but also induce remarkable influences in East Asia via teleconnections.In this study,we assess simulations of the IOBM and IOD,as well as their teleconnections,using the operational seasonal prediction models from the Met Office(Had GEM3) and Beijing Climate Center [BCC CSM1.1(m)].It is demonstrated that the spatial patterns and seasonal cycles are generally reproduced by the control simulations of BCC CSM1.1(m) and Had GEM3,although spectra biases exist.The relationship between the TIO SSTA and El Nio is successfully simulated by both models,including the persistent IOBM warming following El Nio and the IOD–El Nio interactions.BCC CSM1.1(m) and Had GEM3 are capable of simulating the observed local impact of the IOBM,such as the strengthening of the South Asian high.The influences of the IOBM on Yangtze River rainfall are also captured well by both models,although this teleconnection is slightly weaker in BCC CSM1.1(m) due to the underestimation of the northwestern Pacific subtropical high.The local effect of the IOD on East African rainfall is reproduced by both models.However,the remote control of the IOD on rainfall over southwestern China is not clear in either model.It is shown that the realistic simulations of TIO SST modes and their teleconnections give rise to the source of skillful seasonal predictions over China.  相似文献   

8.
利用1991-2017年BCC_CSM1.1(m)模式模拟数据和NCEP/NCAR逐月再分析资料,评估了BCC_CSM1.1(m)对初夏东北冷涡的模拟能力。结果表明:BCC_CSM1.1(m)模式可以对500 hPa位势高度场气候态进行模拟,均方根误差显示该模式对中国东北南部地区500 hPa位势高度场的模拟要优于东北北部地区。EOF第一模态结果显示,该模式可以较好地模拟出500 hPa位势高度场的主要时空变化特征。BCC_CSM1.1(m)模式能够模拟出近27 a东北冷涡指数的上升趋势和年际变化,但模拟的上升趋势较实况偏强,年际变率较实况偏弱。BCC_CSM1.1(m)模式能够模拟出东北冷涡指数的年代际突变,但是对突变开始时间的模拟较实况偏晚。BCC_CSM1.1(m)模式能够模拟出东北冷涡指数和500 hPa位势高度场在东北及其附近地区的显著正相关,不能模拟出东北冷涡指数与东北初夏降水之间的显著负相关。此外,模式东北冷涡指数对东北初夏降水的预测能力十分有限。  相似文献   

9.
This paper reviews recent progress in the development of the Beijing Climate Center Climate System Model(BCC-CSM) and its four component models(atmosphere,land surface,ocean,and sea ice).Two recent versions are described:BCC-CSM1.1 with coarse resolution(approximately 2.8125°×2.8125°) and BCC-CSM1.1(m) with moderate resolution(approximately 1.125°×1.125°).Both versions are fully coupled climate-carbon cycle models that simulate the global terrestrial and oceanic carbon cycles and include dynamic vegetation.Both models well simulate the concentration and temporal evolution of atmospheric CO_2 during the 20th century with anthropogenic CO2 emissions prescribed.Simulations using these two versions of the BCC-CSM model have been contributed to the Coupled Model Intercomparison Project phase five(CMIP5) in support of the Intergovernmental Panel on Climate Change(IPCC) Fifth Assessment Report(AR5).These simulations are available for use by both national and international communities for investigating global climate change and for future climate projections.Simulations of the 20th century climate using BCC-CSMl.l and BCC-CSMl.l(m) are presented and validated,with particular focus on the spatial pattern and seasonal evolution of precipitation and surface air temperature on global and continental scales.Simulations of climate during the last millennium and projections of climate change during the next century are also presented and discussed.Both BCC-CSMl.l and BCC-CSMl.l(m) perform well when compared with other CMIP5 models.Preliminary analyses indicate that the higher resolution in BCC-CSM1.1(m) improves the simulation of mean climate relative to BCC-CSMl.l,particularly on regional scales.  相似文献   

10.
In this paper, major features of the upper-tropospheric jet streams simulated by a coupled Climate System Model BCC_CSM1.1 are evaluated through comparison with the NCEP/NCAR reanalysis. The jet streams consist of the East Asian subtropical jet (EASJ) and the East Asian polar-front jet (EAPJ). Associated stationary wave and synoptic-scale transient eddy activities (STEA) are also examined. The results show that the climatological positions of the westerly jet streams are well captured by BCC_CSM1.1, but with slight intensity biases. Statistics from the 6-h model outputs reveal that the jet core number (JCN) of ESPJ is significantly underestimated. Examination of the simulated seasonal evolution of the westerly jet stream indicates that the model has produced a westward movement of the EASJ core in May, one month earlier than that in the reanalysis. Analysis of stationary wave activities shows that the overestimated meridional wind component may have caused considerable enhancement of meridional momentum and heat transport. The stationary Rossby wave represented by the wave activity flux at the southern flank of the Tibetan Plateau is favorable to the growth of asymmetric zonal wind and the multiple-center pattern of JCN. Unlike the stationary wave heat flux transport, the model tends to systematically generate weaker transient heat flux over East Asia. Further analysis of STEA exhibits a general consistent pattern between the simulation and the reanalysis, while the intensity of the northern STEA branch associated with the EAPJ is greatly reduced. The deficiencies of eddy momentum and heat flux transport and accompanied eddy forcing may contribute to the biases of the simulated upper-tropospheric jet streams, suggesting the potential importance of midlatitude internal atmospheric dynamics in shaping the tropospheric general circulation, which is not yet fully and accurately resolved in the current BCC_CSM1.1.  相似文献   

11.
利用1979—2005年OAFlux (Objectively Analyzed air-sea Fluxes) 观测资料以及CMIP5的15个耦合模式的模拟结果,评估了BCC_CSM1.1(m) 模式对热带太平洋年平均潜热通量气候态和变化趋势的模拟能力,并分析造成趋势偏差的可能原因。结果表明:BCC_CSM1.1(m) 模式模拟热带太平洋年平均潜热通量气候态在各纬度上差异较大, 其中在赤道的模拟能力较佳,而在10°N和8°S附近模拟偏差较大;BCC_CSM1.1(m) 模式对热带太平洋年平均潜热通量趋势的模拟能力一般,造成趋势偏差的主要原因是该模式低估了风速对潜热通量的局地贡献以及它对风速的非局地贡献的模拟存在较大偏差。此外,该模式未能较好地模拟出风速对全球变暖响应。因此,BCC_CSM1.1(m) 模式对热带太平洋年平均潜热通量趋势模拟的改进需加强其对风速模拟的改进。  相似文献   

12.
The summer Asian–Pacific Oscillation(APO) is a major teleconnection pattern that reflects the zonal thermal contrast between East Asia and the North Pacific in the upper troposphere. The performance of Beijing Climate Center Climate System Models(BCC CSMs) with different horizontal resolutions, i.e., BCC CSM1.1 and BCC CSM1.1(m), in reproducing APO interannual variability, APO-related precipitation anomalies, and associated atmospheric circulation anomalies, is evaluated.The results show that BCC CSM1.1(m) can successfully capture the interannual variability of the summer APO index. It is also more capable in reproducing the APO's spatial pattern, compared to BCC CSM1.1, due to its higher horizontal resolution. Associated with a positive APO index, the northward-shifted and intensified South Asian high, strengthened extratropical westerly jet, and tropical easterly jet in the upper troposphere, as well as the southwesterly monsoonal flow over North Africa and the Indian Ocean in the lower troposphere, are realistically represented by BCC CSM1.1(m), leading to an improvement in reproducing the increased precipitation over tropical North Africa, South Asia, and East Asia, as well as the decreased precipitation over subtropical North Africa, Japan, and North America. In contrast, these features are less consistent with observations when simulated by BCC CSM1.1. Regression analysis further indicates that surface temperature anomalies over the North Pacific and the southern and western flanks of the Tibetan Plateau are reasonably reproduced by BCC CSM1.1(m), which contributes to the substantial improvement in the simulation of the characteristics of summer APO compared to that of BCC CSM1.1.  相似文献   

13.
针对江苏夏季旱涝和高温热浪等异常气候的预测难题,以江苏夏季站点降水和气温为预测目标,建立了一种基于全球动力模式BCC_CSM1.1(m)和最优可预测气候模态和异常相对倾向(SMART)原理结合的统计降尺度季节气候预测方法。利用历史观测资料和SVD方法提取决定中国夏季降水异常相对倾向的同期热带地区向外长波辐射(Outgoing Longwave Radiation,OLR)和北半球中高纬500 hPa位势高度场异常相对倾向的最优可预测气候模态,并利用逐步回归法构建其与同期江苏站点降水和气温异常相对倾向同期关系的统计降尺度模型;将动力模式对最优可预测气候模态的预测带入统计降尺度模型,实现对区域降水和气温异常相对倾向的预测;最后通过引入观测的近期背景异常实现对降尺度的降水和气温总距平的预测。通过对1991—2019年江苏夏季降水和气温的回报检验表明,本文建立的统计降尺度模型效果较BCC_CSM1.1(m)动力模式的直接预测效果有显著提高,为区域精细化季节气候预测提供了一种有效的手段。  相似文献   

14.
This study evaluates the ability of four versions BCC (Beijing Climate Center or National Climate Center) models (BCC_AGCM2.1, BCC_AGCM2.2, BCC_CSM1.1 and BCC_CSM1.1m) in simulating the MJO phenomenon using the outputs of the AMIP (Atmospheric Model Intercomparison Project) and historical runs. In general, the models can simulate some major characteristics of the MJO, such as the intensity, the periodicity, the propagation, and the temporal/spatial evolution of the MJO signals in the tropics. There are still some biases between the models and the observation/reanalysis data, such as the overestimated total intraseasonal variability, but underestimated MJO intensity, shorter significant periodicity, and excessive westward propagation. The differences in the ability of simulating the MJO between AMIP and historical experiments are also significant. Compared to the AMIP runs, the total intraseasonal variability is reduced and more realistic, however the ratio between the MJO and its westward counterpart decreases in the historical runs. This unrealistic simulation of the zonal propagation might have been associated with the greater mean precipitation over the Pacific and corresponded to the exaggeration of the South Pacific Convergence Zone structure in precipitation mean state. In contrast to the T42 versions, the improvement of model resolution demonstrate more elaborate topography, but the enhanced westward propagation signals over the Arabia Sea followed. The underestimated (overestimated) MJO variability over eastern Indian Ocean (Pacific) was assumed to be associated with the mean state. Three sets of sensitive experiments using BCC_CSM1.1m turn out to support this argument.  相似文献   

15.
黄淮地区夏季降水的统计降尺度预测   总被引:1,自引:2,他引:1       下载免费PDF全文
利用1991-2011年黄淮地区夏季降水、NCEP/NCAR再分析资料和国家气候中心第2代动力气候模式(BCC_CSM1.1m)夏季回报结果,研究黄淮地区夏季降水降尺度预测模型和可预报性来源。诊断发现,黄淮地区夏季降水与同期南亚高压、乌拉尔山附近阻塞高压、西风急流、西太平洋赤道上空200 hPa纬向风场呈明显正相关。分析BCC_CSM1.1m对夏季环流的回报结果发现,模式对200 hPa和500 hPa位势高度场、200 hPa纬向风场和850 hPa经向风场上影响黄淮地区夏季降水的部分关键区域有较好的模拟能力。利用模式预报技巧较高且对黄淮地区夏季降水的影响有物理含义的环流特征作为预测因子,对比预测因子进行独立性筛选前后分别建立的降尺度预测模型发现,黄淮地区夏季降水预测与实况的距平符号一致率由61%提高到72%。预测技巧来源分析发现,降尺度预测能力与BCC_CSM1.1m对影响黄淮地区夏季降水的3个关键因子乌拉尔山附近环流、南亚高压、西太平洋赤道上空西风强弱的预测技巧密切相关,尤其是模式对西太平洋赤道上空西风的模拟能力起到决定性作用。  相似文献   

16.
利用大气能量循环框图,对比分析中国科学院大气物理研究所大气科学与地球流体力学数值模拟国家重点实验室(LASG/IAP)全球海-陆-气耦合系统模式(GOALS)两个版本(GOALS-2和GOALS-4),以及观测的全球平均大气能量循环的主要特征,并从能量循环贮蓄和转换项的纬向平均贡献去解释全球积分值改善和转坏的原因,以及诊断分析参数化方案变化后产生的影响.结果表明:模式的两个版本基本上能正确地模拟出全球能量循环的主要特征.旧版本GOALS-2能较好模拟全球积分值,常常是不同符号局地误差的相互抵消结果.新版本GOALS-4中某种局地过程的改善在一些情况下导致了全球积分值的转坏.引入辐射日变化参数化方案可能对能量循环各参数的局地贡献有着明显的影响.如纬向平均有效位能向瞬变涡动有效位能的斜压转换率、瞬变涡动有效位能向瞬变涡动动能的斜压转换率以及定常涡动动能的局地贡献有明显改善.南极地区不合实际的上升运动,是模拟的纬向平均有效位能与纬向平均动能之间的转换项全球积分值为负数的主要原因.  相似文献   

17.
通过分析BCC_CSM1.1m、ECMWF_SYSTEM5、NCEP_CFSv2、FGOALS_f、FGOALS_s2、PCCSM4气候模式的模式初始误差、空间相关系数和时间距平相关系数,评估国内外6种气候预测模式在西南地区的预测能力.结果表明:FGOALS_s2和ECMWF_SYSTEM5模式的初始误差最大,FGOA...  相似文献   

18.
This paper explores the role of synoptic eddy feedback in the air-sea interaction in the North Atlantic region, particularly the interaction between the North Atlantic Oscillation (NAO) and the North Atlantic sea surface temperature anomalies (SSTA) tripole. A linearized five-layer primitive equation atmospheric model with synoptic eddy and low-frequency flow (SELF) interaction is coupled with a linearized oceanic mixed-layer model to investigate this issue. In this model, the “climatological” storm track/activity (or synoptic eddy activity) is characterized in terms of spatial structures, variances, decay time scales and propagation speeds through the complex empirical orthogonal function (CEOF) analysis on the observed data, which provides a unique tool to investigate the role of synoptic eddy feedback in the North Atlantic air–sea coupling. Model experiments show that the NAO-like atmospheric circulation anomalies can produce tripole-like SSTA in the North Atlantic Ocean, and the tripole-like SSTA can excite a NAO-like dipole with an equivalent barotropic structure in the atmospheric circulation, which suggests a positive feedback between the NAO and the SSTA tripole. This positive feedback makes the NAO/SSTA tripole-like mode be the leading mode of the coupled dynamical system. The synoptic eddy feedback plays an essential role in the origin of the NAO/SSTA tripole-like leading mode and the equivalent barotropic structure in the atmosphere. Without synoptic eddy feedback, the atmosphere has a baroclinic structure in the response field to the tripole-like SSTA forcing, and the leading mode of the dynamic system does not resemble NAO/SSTA tripole pattern.  相似文献   

19.
Response of the Kuroshio Current to Eddies in the Luzon Strait   总被引:1,自引:0,他引:1       下载免费PDF全文
The impact of eddies on the Kuroshio Current in the Luzon Strait (LS) area is investigated by using the sea surface height anomaly (SSHA) satellite observation data and the sea surface height (SSH) assimilation data. The influence of the eddies on the mean current depends upon the type of eddies and their relative position. The mean current is enhanced (weakened) as the cyclonic (anticyclonic) eddy becomes slightly far from it, whereas it is weakened (enhanced) as the cyclonic (anticyclonic) eddy moves near or within the position of the mean current; this is explained as the eddy-induced meridional velocity and geostrophic flow relationship. The anticyclonic (cyclonic) eddy can increase (decrease) the mean meridional flow due to superimposition of the eddy-induced meridional flow when the eddy is within the region of the mean current. However, when the eddy is slightly far from the mean current region, the anticyclonic (cyclonic) eddy tends to decrease (increase) the zonal gradient of the SSH, which thus results in weakening (strengthening) of the mean current in the LS region.  相似文献   

20.
东亚夏季风模式跨季预测的EOF-相似误差订正   总被引:1,自引:1,他引:0       下载免费PDF全文
利用国家气候中心第2代季节气候预测模式BCC_CSM1.1(m) 的1991—2010年每年2月起报的历史回算资料集,考察模式对于5个夏季风指数的预测能力,并通过发展基于经验正交函数分解与相似分析的EOF-相似误差订正方法,对5个夏季风指数的模式预测进行再修正。交叉检验和独立样本检验结果表明:该模式对1991—2010年东亚夏季风指数与西北太平洋夏季风指数预测技巧较高;EOF-相似误差订正方法适用于模式预测技巧较低的指数,这些指数经订正后预测效果均有不同程度改进,而预测技巧较高的夏季风指数经订正后改进效果不明显;在交叉检验中,线性部分订正多优于非线性部分订正效果,而对于独立试报的年份,非线性部分订正多优于线性部分订正效果,显示出良好的应用前景。  相似文献   

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

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