首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.

Global Climate Models (GCMs) generally exhibit significant biases in the representation of large-scale atmospheric circulation. Even after a sensible bias adjustment these errors remain and are inherited to some extent by the derived downscaling products, impairing the credibility of future regional projections. In this study we perform a process-based evaluation of state-of-the-art GCMs from CMIP5 and CMIP6, with a focus on the simulation of the synoptic climatological patterns having a most prominent effect on the European climate. To this aim, we use the Lamb Weather Type Classification (LWT, Lamb British isles weather types and a register of the daily sequence 736 of circulation patterns 1861-1971. METEOROL OFF, GEOPHYS MEM; 737 GB; DA 1972; NO 116; PP 1-85; BIBL 2P1/2, 1972), a subjective classification of circulation weather types constructed upon historical simulations of daily mean sea level pressure. Observational uncertainty has been taken into account by considering four different reanalysis products of varying characteristics. Our evaluation unveils an overall improvement of salient atmospheric circulation features consistent across observational references, although this is uneven across models and large frequency biases still remain for the main LWTs. Some CMIP6 models attain similar or even worse results than their CMIP5 counterparts, although in most cases consistent improvements have been found, demonstrating the ability of the new models to better capture key synoptic conditions. In light of the large differences found across models, we advocate for a careful selection of driving GCMs in downscaling experiments with a special focus on large-scale atmospheric circulation aspects.

  相似文献   

2.
One of the main sources of uncertainty in estimating climate projections affected by global warming is the choice of the global climate model (GCM). The aim of this study is to evaluate the skill of GCMs from CMIP3 and CMIP5 databases in the north-east Atlantic Ocean region. It is well known that the seasonal and interannual variability of surface inland variables (e.g. precipitation and snow) and ocean variables (e.g. wave height and storm surge) are linked to the atmospheric circulation patterns. Thus, an automatic synoptic classification, based on weather types, has been used to assess whether GCMs are able to reproduce spatial patterns and climate variability. Three important factors have been analyzed: the skill of GCMs to reproduce the synoptic situations, the skill of GCMs to reproduce the historical inter-annual variability and the consistency of GCMs experiments during twenty-first century projections. The results of this analysis indicate that the most skilled GCMs in the study region are UKMO-HadGEM2, ECHAM5/MPI-OM and MIROC3.2(hires) for CMIP3 scenarios and ACCESS1.0, EC-EARTH, HadGEM2-CC, HadGEM2-ES and CMCC-CM for CMIP5 scenarios. These models are therefore recommended for the estimation of future regional multi-model projections of surface variables driven by the atmospheric circulation in the north-east Atlantic Ocean region.  相似文献   

3.
动力降尺度CMIP5的2006—2035新疆夏季降水变化   总被引:1,自引:0,他引:1  
腊孟珂  周洋  朱汉辰  董轩 《气象科学》2019,39(3):413-420
评估了CMIP5 9个模式对新疆夏季降水的模拟效果,从中选取5个较好的模式结果,利用RegCM4动力降尺度再集合平均。通过对CMIP5 9个模式的评估,发现CanCM4、CMCC-CM、CNRM-CM5、HadCM3、MIROC4h对新疆夏季降水的模拟较好。进一步分析动力降尺度和集合平均结果,发现动力降尺度能更好地描述新疆复杂地形造成的降水,同时发现2006—2035年比1976—2005年新疆夏季的总体降水有所减少,存在明显的减少趋势。在空间分布上主要表现为天山山脉和其南侧降水显著减少,昆仑山山脉北缘降水显著增加,这与2006—2035年新疆东北侧500 hPa出现的位势高度正异常和与之相对应的反气旋式环流异常有关。这一环流异常造成了天山山脉上空水汽的显著辐散和昆仑山山脉北缘上空水汽的显著辐合。  相似文献   

4.
潘延  张洋  李舒婷 《气象科学》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中有更多模式可以模拟出极地—高原之间的遥相关波列,且对波列结构的模拟更完整。  相似文献   

5.
Theoretical and Applied Climatology - This study compared precipitation projections of CMIP5 and CMIP6 GCMs over Yulin City, China. The performance of CMIP5 and CMIP6 GCMs in replicating Global...  相似文献   

6.
We assess the ability of Global Climate Models participating in phases 3 and 5 of the Coupled Model Intercomparison Project (CMIP3 and CMIP5) to simulate observed annual precipitation cycles over the Caribbean. Compared to weather station records and gridded observations, we find that both CMIP3 and CMIP5 models can be grouped into three categories: (1) models that correctly simulate a bimodal distribution with two rainfall maxima in May–June and September–October, punctuated by a mid-summer drought (MSD) in July–August; (2) models that reproduce the MSD and the second precipitation maxima only; and (3) models that simulate only one precipitation maxima, beginning in early summer. These categories appear related to model simulation of the North Atlantic Subtropical High (NASH) and sea surface temperature (SST) in the Caribbean Sea and Gulf of Mexico. Specifically, models in category 2 tend to anticipate the westward expansion of the NASH into the Caribbean in early summer. Early onset of NASH results in strong moisture divergence and MSD-like conditions at the time of the May–June observed precipitation maxima. Models in category 3 tend to have cooler SST across the region, particularly over the central Caribbean and the Gulf of Mexico, as well as a weaker Caribbean low-level jet accompanying a weaker NASH. In these models, observed June-like patterns of moisture convergence in the central Caribbean and the Central America and divergence in the east Caribbean and the Gulf of Mexico persist through September. This analysis suggests systematic biases in model structure may be responsible for biases in observed precipitation variability over the Caribbean and more confidence may be placed in the precipitation simulated by the GCMs that are able to correctly simulate seasonal cycles of SST and NASH.  相似文献   

7.
European temperatures and their projected changes under the 8.5 W/m2 Representative Concentration Pathway scenario are evaluated in an ensemble of 33 global climate models participating in the fifth phase of the Coupled Model Intercomparison Project (CMIP5). Respective contributions of large-scale dynamics and local processes to both biases and changes in temperatures, and to the inter-model spread, are then investigated from a recently proposed methodology based on weather regimes. On average, CMIP5 models exhibit a cold bias in winter, especially in Northern Europe. They overestimate summer temperatures in Central Europe, in association with a greater diurnal range than observed. The projected temperature increase is stronger in summer than in winter, with the highest summer warming occurring over Mediterranean regions. Links between biases and sensitivities are evidenced in winter, suggesting a potential influence of snow cover biases on the projected surface warming. A brief analysis of daily temperature extremes suggests that the intra-seasonal variability is projected to decrease (slightly increase) in winter (summer). Then, in order to understand model discrepancies in both present-day and future climates, we disentangle effects of large-scale atmospheric dynamics and regional physical processes. In particular, in winter, CMIP5 models simulate a stronger North-Atlantic jet stream than observed and, in contrast with CMIP3 results, the majority of them suggests an increased frequency of the negative phase of the North-Atlantic Oscillation under future warming. While large-scale circulation only has a minor contribution to ensemble-mean biases or changes, which are primarily dominated by non-dynamical processes, it substantially affects the inter-model spread. Finally, other sources of uncertainties, including the North-Atlantic warming and local radiative feedbacks related to snow cover and clouds, are briefly discussed.  相似文献   

8.
CMIP3气候模式对东亚冬季大气环流模拟能力的评估   总被引:1,自引:0,他引:1  
利用1960—1999年ECMWF月平均再分析资料(ERA40)和耦合模式比较计划(Phase 3 of the CoupledModel Intercomparison Project,简称CMIP3)21个气候耦合模式对20世纪气候模拟试验的模式结果,从气候态和年际变化两个方面,评估了CMIP3气候模式对东亚冬季大气环流的模拟能力。结果表明:(1)模式对东亚地区冬季海平面气压、850 hPa纬向风、经向风和500 hPa高度场气候态的模拟存在不同程度的偏差,但均能较好模拟出上述要素气候态的空间分布特征。总体而言,模式对500 hPa高度场气候态的模拟效果最好,而对850 hPa经向风的模拟效果较差。(2)模式基本上能抓住近40年来东亚地区冬季500 hPa高度场的主要变化特征,但基本上不能模拟出冬季海平面气压、850 hPa纬向风和经向风的变化特征。此外,模式对阿留申低压、蒙古高压和东亚冬季风强度的变化特征几乎没有模拟能力。  相似文献   

9.
Experimental outputs of 11 Atmospheric Model Intercomparison Project (AMIP) models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) are analyzed to assess the atmospheric circulation anomaly over Northern Hemisphere induced by the anomalous rainfall over tropical Pacific and Indian Ocean during boreal winter.The analysis shows that the main features of the interannual variation of tropical rainfall anomalies,especially over the Central Pacific (CP) (5°S-5°N,175°E-135°W) and Indo-western Pacific (IWP) (20°S-20°N,110°-150°E) are well captured in all the CMIP5/AMIP models.For the IWP and western Indian Ocean (WIO) (10°S-10°N,45°-75°E),the anomalous rainfall is weaker in the 11 CMIP5/AMIP models than in the observation.During El Ni(n)o/La Ni(n)a mature phases in boreal winter,consistent with observations,there are geopotential height anomalies known as the Pacific North American (PNA) pattern and Indo-western Pacific and East Asia (IWPEA) pattern in the upper troposphere,and the northwestern Pacific anticyclone (cyclone) (NWPA) in the lower troposphere in the models.Comparison between the models and observations shows that the ability to simulate the PNA and NWPA pattern depends on the ability to simulate the anomalous rainfall over the CP,while the ability to simulate the IWPEA pattern is related to the ability to simulate the rainfall anomaly in the IWP and WIO,as the SST anomaly is same in AMIP experiments.It is found that the tropical rainfall anomaly is important in modeling the impact of the tropical Indo-Pacific Ocean on the extratropical atmospheric circulation anomaly.  相似文献   

10.
For the construction of regional climate change scenarios spanning a relevant fraction of the spread in climate model projections, an inventory of major drivers of regional climate change is needed. For the Netherlands, a previous set of regional climate change scenarios was based on the decomposition of local temperature/precipitation changes into components directly linked to the level of global warming, and components related to changes in the regional atmospheric circulation. In this study this decomposition is revisited utilizing the extensive modelling results from the CMIP5 model ensemble in support for the 5th IPCC assessment. Rather than selecting a number of GCMs based on performance metrics or relevant response features, a regression technique was developed to utilize all available model projections. The large number of projections allows a quantification of the separate contributions of emission scenarios, systematic model responses and natural variability to the total likelihood range. Natural variability plays a minor role in modelled differences in the global mean temperature response, but contributes for up to 50 % to the range of mean sea level pressure responses and local precipitation. Using key indicators (“steering variables”) for the temperature and circulation response, the range in local seasonal mean temperature and precipitation responses can be fairly well reproduced.  相似文献   

11.
A parallel comparison is made of the circulation climatology and the leading oscillation mode of the northern winter stratosphere among six reanalysis products and 24 CMIP5(Coupled Model Intercomparison Project Phase 5) models. The results reveal that the NCEP/NCAR, NECP/DOE, ERA40, ERA-Interim and JRA25 reanalyses are quite consistent in describing the climatology and annual cycle of the stratospheric circulation. The 20 CR reanalysis, however, exhibits a remarkable"cold pole" bias accompanied by a much stronger stratospheric polar jet, similar as in some CMIP5 models. Compared to the1–2 month seasonal drift in most coupled general circulation models(GCMs), the seasonal cycle of the stratospheric zonal wind in most earth system models(ESMs) agrees very well with reanalysis. Similar to the climatology, the amplitude of Polar Vortex Oscillation(PVO) events also varies among CMIP5 models. The PVO amplitude in most GCMs is relatively weaker than in reanalysis, while that in most of the ESMs is more realistic. In relation to the "cold pole" bias and the weaker oscillation in some CMIP5 GCMs, the frequency of PVO events is significantly underestimated by CMIP5 GCMs; while in most ESMs, it is comparable to that in reanalysis. The PVO events in reanalysis(except in 20CR) mainly occur from mid-winter to early spring(January–March); but in some of the CMIP5 models, a 1–2 month delay exists, especially in most of the CMIP5 GCMs. The long-term trend of the PVO time series does not correspond to long-term changes in the frequency of PVO events in most of the CMIP5 models.  相似文献   

12.
Based on 15 Coupled Model Intercomparison Project (CMIP) phase 3 (CMIP3) and 32 CMIP phase 5 (CMIP5) models, a detailed diagnosis was carried out to understand what compose the biases in simulation of the Indian Ocean basin mode (IOBM) and its capacitor effect. Cloud-radiation-SST (CRS) feedback and wind-evaporation-SST (WES) feedback are the two major atmospheric processes for SST changes. Most CMIP models simulate a stronger CRS feedback and a weaker WES feedback. During boreal fall of the El Niño/Southern Oscillation developing year and the following spring, there are weak biases of suppressed rainfall anomalies over the Maritime Continent and anomalous anticyclone over South Indian Ocean. Most CMIP models simulate reasonable short wave radiation (SWR) and weaker latent heat flux (LHF) anomalies. This leads to a weak bias of atmospheric processes. During winter, however, the rainfall anomalies are stronger due to west bias, and the anomalous anticyclone is comparable to observations. As such, most models simulate stronger SWR and reasonable LHF anomalies, leading to a strong bias of atmospheric processes. The thermocline feedback is stronger in most models. Though there is a deep bias of climatology thermocline, most models capture reasonable sea surface height-induced SST anomalies. Therefore, the effect of oceanic processes offset the weak bias of atmospheric processes in spring, and the tropical Indian Ocean warming persists into summer. However, anomalous northwest Pacific (NWP) anticyclone is weaker due to weak and west bias of the capacitor effect. The unrealistic western Pacific SST anomalies in models favor the westward extension of Rossby wave from the Pacific, weakening the effect of Kelvin wave from the Indian Ocean. Moreover, the western Pacific warming forces the NWP anticyclone move farther north than observations, suggesting a major forcing from the Pacific. Compared to CMIP3, CMIP5 models simulate the feedbacks more realistically and display less diversity. Thus, the overall performance of CMIP5 models is better than that of CMIP3 models.  相似文献   

13.
 The study seeks to describe one method of deriving information about local daily temperature extremes from larger scale atmospheric flow patterns using statistical tools. This is considered to be one step towards downscaling coarsely gridded climate data from global climate models (GCMs) to finer spatial scales. Downscaling is necessary in order to bridge the spatial mismatch between GCMs and climate impact models which need information on spatial scales that the GCMs cannot provide. The method of statistical downscaling is based on physical interaction between atmospheric processes with different spatial scales, in this case between synoptic scale mean sea level pressure (MSLP) fields and local temperature extremes at several stations in southeast Australia. In this study it was found that most of the day-to-day spatial variability of the synoptic circulation over the Australian region can be captured by six principal components. Using the scores of these PCs as multivariate indicators of the circulation a substantial part of the local daily temperature variability could be explained. The inclusion of temperature persistence noticeably improved the performance of the statistical model. The model established and tested with observations is thought to be finally applied to GCM-simulated pressure fields in order to estimate pressure-related changes in local temperature extremes under altered CO2 conditions. Received: 26 March 1996 / Accepted: 20 September 1996  相似文献   

14.
In this research, tropopause temperature (TT) and tropopause geopotential height (TGH) over the inner-core and environmental regions of all tropical cyclones (TCs) over the northwest of the Indian Ocean (NWIO) from 1990 to 2019 were investigated. To this aim, observational/analysis/reanalysis data and also simulated data from both historical and Representative Concentration Pathway 8.5 (RCP8.5) experiments of some global climate models (GCMs) from the Coupled Model Intercomparision Project (CMIP5) were used. Dynamical and thermo-dynamical environmental factors controlling TC, together with their correlation with different phases of some climatic indices were considered. Results indicated that the eastern part of the NWIO was more favorable for TC genesis and intensification.Lower-level stratospheric (upper-level tropospheric) cooling (warming) was detected over the NWIO during 1990−2019. Over the both inner-core and environmental regions of the NWIO TCs, the coldest tropopause occurred at a CS-Category and the warmest tropopause happened at the first stage of a VSCS event. Over the inner-core (environmental) region, the highest tropopause was detected at the first stage of a CS event (at the end of a VSCS life cycle). A significant majority of the used CMIP5 GCMs produced stratospheric cooling and tropospheric warming trends over the NWIO, similar to those obtained using ERA5 reanalysis dataset. Finally, the decreasing trend of TT over the both inner-core and environmental regions of NWIO TCs together with temperature decreasing trend obtained from the CMIP5 GCMs simulations suggest that the NWIO is prone to experience more TCs, especially the intense ones, in the future.  相似文献   

15.
In order to assess the impact of the mid-tropospheric circulation over the Greenland ice sheet (GrIS) on surface melt, as simulated by the regional climate model MAR, an automatic Circulation type classification (CTC) based on 500?hPa geopotential height from reanalyses is developed. General circulation correlates significantly with the surface melt anomalies for the summers in the period 1958?C2009. The record surface melt events observed during the summers of 2007?C2009 are linked to the exceptional persistence of atmospheric circulations favouring warm air advection. The CTC emphasizes that summer 500?hPa circulation patterns have changed since the beginning of the 2000s; this process is partly responsible for the recent warming observed over the GrIS.  相似文献   

16.
We evaluate the representation of dynamic sea surface height (SSH) fields of 33 global coupled models (GCMs) contributed to the fifth phase of the Coupled Model Intercomparison Project (CMIP5). We use observations from satellite altimetry and basic performance metrics to quantify the ability of the GCMs to replicate observed SSH of the time-mean, seasonal cycle, and inter-annual variability patterns. The time-mean SSH representation has markedly improved from CMIP3 to CMIP5, both in terms of overall reduction in root-mean square differences, and in terms of reduced GCM ensemble spread. Biases of the time-mean SSH field in the Indian and Pacific Ocean equatorial regions are consistent with biases in the zonal surface wind stress fields identified with independent measurements. In the Southern Ocean, the latitude of the maximum meridional gradient of the zonal mean SSH CMIP5 models is shifted equatorward, consistent with the GCMs’ spatial biases in the maximum of the zonal mean westerly surface wind stress fields. However, while the Southern Ocean SSH gradients correlate well with the maximum Antarctic circumpolar current transports, there is no significant correlation with the maximum zonal mean wind stress amplitudes, consistent with recent findings that the eddy parameterisations in GCMs dominate over wind stress amplitudes in this region. There is considerable spread across the CMIP5 ensemble for the seasonal and interannual SSH variability patterns. Because of the short observational period, the interannual variability patterns depend on the time-period over which they are derived, while no such dependency is found for the time-mean patterns. The model performance metrics for SSH presented here provide insight into GCM shortcoming due to inadequate model physics or processes. While the diagnostics of CMIP5 GCM performance relative to observations reveal that some models are clearly better than others, model performance is sensitive to the spatio-temporal scales chosen.  相似文献   

17.
CMIP5模式对中国年平均气温模拟及其与CMIP3模式的比较   总被引:5,自引:0,他引:5  
利用CRUT3v和CN05两套观测资料,评估25个CMIP5模式对1906-2005年中国年平均气温变化的模拟能力,并与CMIP3模式对比。结果表明:1906-2005年中国平均温升速率为0.84℃/100a,CMIP5多模式集合平均模拟的增温率为0.77℃/100a。模式对20世纪后期温升模拟好于前期,仅有两个模式能模拟中国20世纪40年代异常增暖。模式对气温气候态空间分布模拟较好,但在中国西部地区存在最大模拟冷偏差和不确定性。1961-1999年,中国北方增暖大于南方。多模式集合平均可以较好地模拟气温变化线性趋势的空间分布,但对南北气温变化趋势的差异模拟过小。总体说来,在中国平均气温变化趋势、气温气候态空间分布和气温变化趋势空间分布三方面,CMIP5模式都较CMIP3模式有所提高。  相似文献   

18.
Coupled ocean–atmosphere general circulation models (GCMs) lack sufficient resolution to model the regional detail of changes to mean circulation and rainfall with projected climate warming. In this paper, changes in mean circulation and rainfall in GCMs are compared to those in a variable resolution regional climate model, the Conformal Cubic Atmospheric Model (CCAM), under a high greenhouse gas emissions scenario. The study site is Tasmania, Australia, which is positioned within the mid-latitude westerlies of the southern hemisphere. CCAM projects a different response in mean sea level pressure and mid-latitude westerly circulation to climate warming to the GCMs used as input, and shows greater regional detail of the boundaries between regions of increasing and decreasing rainfall. Changes in mean circulation dominate the mean rainfall response in western Tasmania, whereas changes to rainfall in the East Coast are less related to mean circulation changes. CCAM projects an amplification of the dominant westerly circulation over Tasmania and this amplifies the seasonal cycle of wet winters and dry summers in the west. There is a larger change in the strength than in the incidence of westerly circulation and rainfall events. We propose the regional climate model displays a more sensitive atmospheric response to the different rates of warming of land and sea than the GCMs as input. The regional variation in these results highlight the need for dynamical downscaling of coupled general circulation models to finely resolve the influence of mean circulation and boundaries between regions of projected increases and decreases in rainfall.  相似文献   

19.
Joint variable spatial downscaling   总被引:1,自引:0,他引:1  
Joint Variable Spatial Downscaling (JVSD), a new statistical technique for downscaling gridded climatic variables, is developed to generate high resolution gridded datasets for regional watershed modeling and assessments. The proposed approach differs from previous statistical downscaling methods in that multiple climatic variables are downscaled simultaneously and consistently to produce realistic climate projections. In the bias correction step, JVSD uses a differencing process to create stationary joint cumulative frequency statistics of the variables being downscaled. The functional relationship between these statistics and those of the historical observation period is subsequently used to remove GCM bias. The original variables are recovered through summation of bias corrected differenced sequences. In the spatial disaggregation step, JVSD uses a historical analogue approach, with historical analogues identified simultaneously for all atmospheric fields and over all areas of the basin under study. Analysis and comparisons are performed for 20th Century Climate in Coupled Models (20C3M), broadly available for most GCMs. The results show that the proposed downscaling method is able to reproduce the sub-grid climatic features as well as their temporal/spatial variability in the historical periods. Comparisons are also performed for precipitation and temperature with other statistical and dynamic downscaling methods over the southeastern US and show that JVSD performs favorably. The downscaled sequences are used to assess the implications of GCM scenarios for the Apalachicola-Chattahoochee-Flint river basin as part of a comprehensive climate change impact assessment.  相似文献   

20.
Dynamical downscaling has been recognized as a useful tool not only for the climate community, but also for associated application communities such as the environmental and hydrological societies. Although climate projection data are available in lower-resolution general circulation models (GCMs), higher-resolution climate projections using regional climate models (RCMs) have been obtained over various regions of the globe. Various model outputs from RCMs with a high resolution of even as high as a few km have become available with heavy weight on applications. However, from a scientific point of view in numerical atmospheric modeling, it is not clear how to objectively judge the degree of added value in the RCM output against the corresponding GCM results. A key factor responsible for skepticism is based on the fundamental limitations in the nesting approach between GCMs and RCMs. In this article, we review the current status of the dynamical downscaling for climate prediction, focusing on basic assumptions that are scrutinized from a numerical weather prediction (NWP) point of view. Uncertainties in downscaling due to the inconsistencies in the physics packages between GCMs and RCMs were revealed. Recommendations on how to tackle the ultimate goal of dynamical downscaling were also described.  相似文献   

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

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