首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
In this paper, lead-time and spatial dependence in skill for prediction of monthly mean climate variability is analyzed. The analysis is based on a set of extensive hindcasts from the Climate Forecast System at the National Centers for Environmental Prediction. The skill characteristics of initialized predictions is also compared with the AMIP simulations forced with the observed sea surface temperature (SST) to quantify the role of initial versus boundary conditions in the prediction of monthly means. The analysis is for prediction of monthly mean SST, precipitation, and 200-hPa height. The results show a rapid decay in skill with lead time for the atmospheric variables in the extratropical latitudes. Further, after a lead-time of approximately 30?C40?days, the skill of monthly mean prediction is essentially a boundary forced problem, with SST anomalies in the tropical central/eastern Pacific playing a dominant role. Because of the larger contribution from the atmospheric internal variability to monthly time-averages (compared to seasonal averages), skill for monthly mean prediction associated with boundary forcing is also lower. The analysis indicates that the prospects of skillful prediction of monthly means may remain a challenging problem, and may be limited by inherent limits in predictability.  相似文献   

2.
Seasonal predictions of Arctic sea ice have typically been based on statistical regression models or on results from ensemble ice model forecasts driven by historical atmospheric forcing. However, in the rapidly changing Arctic environment, the predictability characteristics of summer ice cover could undergo important transformations. Here global coupled climate model simulations are used to assess the inherent predictability of Arctic sea ice conditions on seasonal to interannual timescales within the Community Climate System Model, version 3. The role of preconditioning of the ice cover versus intrinsic variations in determining sea ice conditions is examined using ensemble experiments initialized in January with identical ice?Cocean?Cterrestrial conditions. Assessing the divergence among the ensemble members reveals that sea ice area exhibits potential predictability during the first summer and for winter conditions after a year. The ice area exhibits little potential predictability during the spring transition season. Comparing experiments initialized with different mean ice conditions indicates that ice area in a thicker sea ice regime generally exhibits higher potential predictability for a longer period of time. In a thinner sea ice regime, winter ice conditions provide little ice area predictive capability after approximately 1?year. In all regimes, ice thickness has high potential predictability for at least 2?years.  相似文献   

3.
A verification framework for interannual-to-decadal predictions experiments   总被引:2,自引:1,他引:1  
Decadal predictions have a high profile in the climate science community and beyond, yet very little is known about their skill. Nor is there any agreed protocol for estimating their skill. This paper proposes a sound and coordinated framework for verification of decadal hindcast experiments. The framework is illustrated for decadal hindcasts tailored to meet the requirements and specifications of CMIP5 (Coupled Model Intercomparison Project phase 5). The chosen metrics address key questions about the information content in initialized decadal hindcasts. These questions are: (1) Do the initial conditions in the hindcasts lead to more accurate predictions of the climate, compared to un-initialized climate change projections? and (2) Is the prediction model’s ensemble spread an appropriate representation of forecast uncertainty on average? The first question is addressed through deterministic metrics that compare the initialized and uninitialized hindcasts. The second question is addressed through a probabilistic metric applied to the initialized hindcasts and comparing different ways to ascribe forecast uncertainty. Verification is advocated at smoothed regional scales that can illuminate broad areas of predictability, as well as at the grid scale, since many users of the decadal prediction experiments who feed the climate data into applications or decision models will use the data at grid scale, or downscale it to even higher resolution. An overall statement on skill of CMIP5 decadal hindcasts is not the aim of this paper. The results presented are only illustrative of the framework, which would enable such studies. However, broad conclusions that are beginning to emerge from the CMIP5 results include (1) Most predictability at the interannual-to-decadal scale, relative to climatological averages, comes from external forcing, particularly for temperature; (2) though moderate, additional skill is added by the initial conditions over what is imparted by external forcing alone; however, the impact of initialization may result in overall worse predictions in some regions than provided by uninitialized climate change projections; (3) limited hindcast records and the dearth of climate-quality observational data impede our ability to quantify expected skill as well as model biases; and (4) as is common to seasonal-to-interannual model predictions, the spread of the ensemble members is not necessarily a good representation of forecast uncertainty. The authors recommend that this framework be adopted to serve as a starting point to compare prediction quality across prediction systems. The framework can provide a baseline against which future improvements can be quantified. The framework also provides guidance on the use of these model predictions, which differ in fundamental ways from the climate change projections that much of the community has become familiar with, including adjustment of mean and conditional biases, and consideration of how to best approach forecast uncertainty.  相似文献   

4.
Given observed initial conditions, how well do coupled atmosphere–ocean models predict precipitation climatology with 1-month lead forecast? And how do the models’ biases in climatology in turn affect prediction of seasonal anomalies? We address these questions based on analysis of 1-month lead retrospective predictions for 21 years of 1981–2001 made by 13 state-of-the-art coupled climate models and their multi-model ensemble (MME). The evaluation of the precipitation climatology is based on a newly designed metrics that consists of the annual mean, the solstitial mode and equinoctial asymmetric mode of the annual cycle, and the rainy season characteristics. We find that the 1-month lead seasonal prediction made by the 13-model ensemble has skills that are much higher than those in individual model ensemble predictions and approached to those in the ERA-40 and NCEP-2 reanalysis in terms of both the precipitation climatology and seasonal anomalies. We also demonstrate that the skill for individual coupled models in predicting seasonal precipitation anomalies is positively correlated with its performances on prediction of the annual mean and annual cycle of precipitation. In addition, the seasonal prediction skill for the tropical SST anomalies, which are the major predictability source of monsoon precipitation in the current coupled models, is closely link to the models’ ability in simulating the SST mean state. Correction of the inherent bias in the mean state is critical for improving the long-lead seasonal prediction. Most individual coupled models reproduce realistically the long-term annual mean precipitation and the first annual cycle (solstitial mode), but they have difficulty in capturing the second annual (equinoctial asymmetric) mode faithfully, especially over the Indian Ocean (IO) and Western North Pacific (WNP) where the seasonal cycle in SST has significant biases. The coupled models replicate the monsoon rain domains very well except in the East Asian subtropical monsoon and the tropical WNP summer monsoon regions. The models also capture the gross features of the seasonal march of the rainy season including onset and withdraw of the Asian–Australian monsoon system over four major sub-domains, but striking deficiencies in the coupled model predictions are observed over the South China Sea and WNP region, where considerable biases exist in both the amplitude and phase of the annual cycle and the summer precipitation amount and its interannual variability are underestimated.  相似文献   

5.
Conflicting results have been presented regarding the link between Arctic sea-ice loss and midlatitude cooling, particularly over Eurasia. This study analyzes uncoupled(atmosphere-only) and coupled(ocean–atmosphere) simulations by the Climate Forecast System, version 2(CFSv2), to examine this linkage during the Northern Hemisphere winter, focusing on the simulation of the observed surface cooling trend over Eurasia during the last three decades. The uncoupled simulations are Atmospheric Model Intercomparison Project(AMIP) runs forced with mean seasonal cycles of sea surface temperature(SST)and sea ice, using combinations of SST and sea ice from different time periods to assess the role that each plays individually,and to assess the role of atmospheric internal variability. Coupled runs are used to further investigate the role of internal variability via the analysis of initialized predictions and the evolution of the forecast with lead time.The AMIP simulations show a mean warming response over Eurasia due to SST changes, but little response to changes in sea ice. Individual runs simulate cooler periods over Eurasia, and this is shown to be concurrent with a stronger Siberian high and warming over Greenland. No substantial differences in the variability of Eurasian surface temperatures are found between the different model configurations. In the coupled runs, the region of significant warming over Eurasia is small at short leads, but increases at longer leads. It is concluded that, although the models have some capability in highlighting the temperature variability over Eurasia, the observed cooling may still be a consequence of internal variability.  相似文献   

6.
Advanced warning of extreme sea level events is an invaluable tool for coastal communities, allowing the implementation of management policies and strategies to minimise loss of life and infrastructure damage. This study is an initial attempt to apply a dynamical coupled ocean–atmosphere model to the prediction of seasonal sea level anomalies (SLA) globally for up to 7 months in advance. We assess the ability of the Australian Bureau of Meteorology’s operational seasonal dynamical forecast system, the Predictive Ocean Atmosphere Model for Australia (POAMA), to predict seasonal SLA, using gridded satellite altimeter observation-based analyses over the period 1993–2010 and model reanalysis over 1981–2010. Hindcasts from POAMA are based on a 33-member ensemble of seasonal forecasts that are initialised once per month for the period 1981–2010. Our results show POAMA demonstrates high skill in the equatorial Pacific basin and consistently exhibits more skill globally than a forecast based on persistence. Model predictability estimates indicate there is scope for improvement in the higher latitudes and in the Atlantic and Southern Oceans. Most characteristics of the asymmetric SLA fields generated by El-Nino/La Nina events are well represented by POAMA, although the forecast amplitude weakens with increasing lead-time.  相似文献   

7.
The unprecedented Zhengzhou heavy rainfall in July 2021 occurred under the background of a northward shift of the western Pacific subtropical high(WPSH). Although the occurrence of this extreme event could not be captured by seasonal predictions, a skillful prediction of the WPSH variation might have warned us of the increased probability of extreme weather events in Central and Northern China. However, the mechanism for the WPSH variation in July2021 and its seasonal predictability are still un...  相似文献   

8.
The latest operational version of the ECMWF seasonal forecasting system is described. It shows noticeably improved skill for sea surface temperature (SST) prediction compared with previous versions, particularly with respect to El Nino related variability. Substantial skill is shown for lead times up to 1?year, although at this range the spread in the ensemble forecast implies a loss of predictability large enough to account for most of the forecast error variance, suggesting only moderate scope for improving long range El Nino forecasts. At shorter ranges, particularly 3?C6?months, skill is still substantially below the model-estimated predictability limit. SST forecast skill is higher for more recent periods than earlier ones. Analysis shows that although various factors can affect scores in particular periods, the improvement from 1994 onwards seems to be robust, and is most plausibly due to improvements in the observing system made at that time. The improvement in forecast skill is most evident for 3-month forecasts starting in February, where predictions of NINO3.4 SST from 1994 to present have been almost without fault. It is argued that in situations where the impact of model error is small, the value of improved observational data can be seen most clearly. Significant skill is also shown in the equatorial Indian Ocean, although predictive skill in parts of the tropical Atlantic are relatively poor. SST forecast errors can be especially high in the Southern Ocean.  相似文献   

9.
There are two main approaches for dealing with model biases in forecasts made with initialized climate models. In full-field initialization, model biases are removed during the assimilation process by constraining the model to be close to observations. Forecasts drift back towards the model’s preferred state, thereby re-establishing biases which are then removed with an a posterior lead-time dependent correction diagnosed from a set of historical tests (hindcasts). In anomaly initialization, the model is constrained by observed anomalies and deviates from its preferred climatology only by the observed variability. In theory, the forecasts do not drift, and biases may be removed based on the difference between observations and independent model simulations of a given period. Both approaches are currently in use, but their relative merits are unclear. Here we compare the skill of each approach in comprehensive decadal hindcasts starting each year from 1960 to 2009, made using the Met Office decadal prediction system. Both approaches are more skilful than climatology in most regions for temperature and some regions for precipitation. On seasonal timescales, full-field initialized hindcasts of regional temperature and precipitation are significantly more skilful on average than anomaly initialized hindcasts. Teleconnections associated with the El Niño Southern Oscillation are stronger with the full-field approach, providing a physical basis for the improved precipitation skill. Differences in skill on multi-year timescales are generally not significant. However, anomaly initialization provides a better estimate of forecast skill from a limited hindcast set.  相似文献   

10.
The impact of realistic atmospheric initialisation on the seasonal prediction of tropical Pacific sea surface temperatures is explored with the Predictive Ocean–Atmosphere Model for Australia (POAMA) dynamical seasonal forecast system. Previous versions of POAMA used data from an Atmospheric Model Intercomparison Project (AMIP)-style simulation to initialise the atmosphere for the hindcast simulations. The initial conditions for the hindcasts did not, therefore, capture the true intra-seasonal atmospheric state. The most recent version of POAMA has a new Atmosphere and Land Initialisation scheme (ALI), which captures the observed intra-seasonal atmospheric state. We present the ALI scheme and then compare the forecast skill of two hindcast datasets, one with AMIP-type initialisation and one with realistic initial conditions from ALI, focussing on the prediction of El Niño. For eastern Pacific (Niño3) sea surface temperature anomalies (SSTAs), both experiments beat persistence and have useful SSTA prediction skill (anomaly correlations above 0.6) at all lead times (forecasts are 9 months duration). However, the experiment with realistic atmospheric initial conditions from ALI is an improvement over the AMIP-type initialisation experiment out to about 6 months lead time. The improvements in skill are related to improved initial atmospheric anomalies rather than an improved initial mean state (the forecast drift is worse in the ALI hindcast dataset). Since we are dealing with a coupled system, initial atmospheric errors (or differences between experiments) are amplified though coupled processes which can then lead to long lasting errors (or differences).  相似文献   

11.
Decadal predictability and forecast skill   总被引:2,自引:1,他引:1  
The “potential predictability” of the climate system is the upper limit of available forecast skill and can be characterized by the ratio p of the predictable variance to the total variance. While the potential predictability of the actual climate system is unknown its analog q may be obtained for a model of the climate system. The usual correlation skill score r and the mean square skill score M are functions of p in the case of actual forecasts and potential correlation ρ and potential mean square skill score $\mathcal{M}$ are the same functions of q in the idealized model context. In the large ensemble limit the connection between model-based potential predictability and skill scores is particularly straightforward with $q=\rho^{2}=\mathcal{M}.$ Decadal predictions of annual mean temperature produced with the Canadian Centre for Climate Modelling and Analysis coupled climate model are analyzed for information on decadal climate predictability and actual forecast skill. Initialized forecast results are compared with the results of uninitialized climate simulations. Model-based values of potential predictability q and potential correlation skill ρ are obtained and ρ is compared with the actual forecast correlation skill r. The skill of externally forced and internally generated components of the variability are separately estimated. As expected, ρ > r and both decline with forecast range τ, at least for the first five years. The decline of skill is associated mainly with the decline of the skill of the internally generated component. The potential and actual skill of a forecast of time-averaged temperature depends on the averaging period. The skill of uninitialized simulations is low for short averaging times and increases as averaging time increases. By contrast, skill is high at short averaging times for forecasts initialized from observations and declines as averaging times increase to about three years, then increases somewhat at longer averaging times. The skills of the initialized forecasts and uninitialized simulations begin to converge for longer averaging times. The potential correlation skill ρ of the externally forced component of temperature is largest at tropical latitudes and the skill of the internally generated component is largest over the North Atlantic, parts of the Southern Ocean and to some extent the North Pacific. Potential skill over extratropical land is somewhat weaker than over oceans. The distribution of actual correlation skill r is broadly similar to that of potential skill for the externally forced component but less so for the internally generated component. Differences in potential and actual skill suggest where improvements in the forecast system might be found.  相似文献   

12.
Ma  Youwei  Li  Jianping  Zhang  Shaoqing  Zhao  Haoran 《Climate Dynamics》2021,56(11):3489-3509

Of great importance for guiding numerical weather and climate predictions, understanding predictability of the atmosphere in the ocean − atmosphere coupled system is the first and critical step to understand predictability of the Earth system. However, previous predictability studies based on prefect model assumption usually depend on a certain model. Here we apply the predictability study with the Nonlinear Local Lyapunov Exponent and Attractor Radius to the products of multiple re-analyses and forecast models in several operational centers to realize general predictability of the atmosphere in the Earth system. We first investigated the predictability characteristics of the atmosphere in NCEP, ECMWF and UKMO coupled systems and some of their uncoupled counterparts and other uncoupled systems. Although the ECMWF Integrated Forecast System shows higher skills in geopotential height over the tropics, there is no certain model providing the most precise forecast for all variables on all levels and the multi-model ensemble not always outperforms a single model. Improved low-frequency signals from the air − sea and stratosphere − troposphere interactions that extend predictability of the atmosphere in coupled system suggests the significance of air − sea coupling and stratosphere simulation in practical forecast development, although uncertainties exist in the model representation for physical processes in air − sea interactions and upper troposphere. These inspire further exploration on predictability of ocean and stratosphere as well as sea − ice and land processes to advance our understanding of interactions of Earth system components, thus enhancing weather − climate prediction skills.

  相似文献   

13.
年代际气候预测计划(DCPP)是第六次国际耦合模式比较计划(CMIP6)的子计划之一,其目标是利用多模式开展气候系统年代际预测、可预测性和变率机制研究。DCPP设计了3组试验,即年代际回报试验、预报试验以及理解年代际变率机制和可预测性的敏感性试验。目前有21个模式拟参与DCPP计划,其中包括5个来自中国的模式。DCPP将推动解决气候系统从年际到年代际尺度预测相关的多项科学问题,评估当前气候预测系统预报技巧,挖掘潜在可预报性,研究长时间尺度气候变率形成机制,提供对科学和社会有用的预测产品。  相似文献   

14.
基于乌鲁木齐区域数值预报业务系统,运用Ts和Bias评分方法,对2012年9月1日—2015年8月31日逐日2个起报时次的逐6 h累积降水量的年与季节预报性能进行检验,并从空间上分析了2015年全疆站点逐6 h累积降水量在4个预报时段的评分特征。结果表明:(1)2个起报时次的降水评分相差较小,00 UTC起报略优于12 UTC起报,2015年系统改进了白天大量级降水的空报现象。(2)系统对晴雨预报较为准确,Bias接近1,空报、漏报率很小;随着降水阈值的升高,Ts评分减小,Bias变幅增大,空、漏报率也随之增加。系统对强降水过程以漏报为主。(3)系统的降水预报能力存在季节差异,夏季Ts评分最高,秋季次之,冬季最小;随时间模式对四季降水预报能力均有提高,降低了冬季大量级降水的漏报率和夏季大量级降水的空报率。(4)在新疆地区,08—14 BT(Beijing Time)、14—20 BT、20—次日02 BT空报站点数多于漏报,14—20 BT空报率最高;在02—08 BT整体呈漏报。(5)各站点整体来看,白天Ts评分高于夜间,山区及邻近地区评分高于平原地区;西天山评分略优于东天山,夜间晴雨预报有天山北坡漏报、南坡空报的趋势。  相似文献   

15.
The impact of transient eddies on extratropical seasonal-mean prediction and predictability was examined using DEMETER seasonal prediction data. Two distinct groups were found among the seven DEMETER models based on the simulated properties of their climatological state: (1) models of a strong jet stream and strong transient activity (strong transient models), which is close to the observed intensity, and (2) models of a weak jet stream and weak transient activity (weak transient models). In addition to climatology, the strong transient models tend to predict strong Pacific North American (PNA) patterns, whereas the weak transient models predict weak PNA patterns. Here we demonstrate that these differences mainly result from differences in the eddy feedback intensity. Due to synoptic eddy feedback, the strong transient models exhibit not only strong signal variance but also strong noise variance compared with those of the weak transient models. Interestingly two groups of models show the potential predictability of deterministic forecast, measured by the signal to noise ratio, which is similar to each other. However, the strong transient models produce the error to spread ratio smaller than that of the weak transient models, implying that the former models produce a more reliable spread for the probabilistic forecast. This study implies that a better representation of transient statistics is needed to improve the extratropical predictability of the dynamical seasonal prediction.  相似文献   

16.
Through the analysis of ensembles of coupled model simulations and projections collected from CMIP3 and CMIP5, we demonstrate that a fundamental spatial scale limit might exist below which useful additional refinement of climate model predictions and projections may not be possible. That limit varies among climate variables and from region to region. We show that the uncertainty (noise) in surface temperature predictions (represented by the spread among an ensemble of global climate model simulations) generally exceeds the ensemble mean (signal) at horizontal scales below 1000 km throughout North America, implying poor predictability at those scales. More limited skill is shown for the predictability of regional precipitation. The ensemble spread in this case tends to exceed or equal the ensemble mean for scales below 2000 km. These findings highlight the challenges in predicting regionally specific future climate anomalies, especially for hydroclimatic impacts such as drought and wetness.  相似文献   

17.
对月平均大气环流预报试验、季度预测和中国汛期降水预测进行了总结。结果表明气候预测的对象必须是要素的时间平均场。利用数值模拟进行气候预测是今后的主要发展方向,而季度预测技巧的提高依赖于对物理参数化和物理机制的研究。最后,讨论了季平均气温和季总降水的可预报性问题,即时效性和准确率。  相似文献   

18.
The behavior of the water cycle in the Coupled Forecast System version 2 reforecasts and reanalysis is examined. Attention is focused on the evolution of forecast biases as the lead-time changes, and how the lead-time dependent model climatology differs from the reanalysis. Precipitation biases are evident in both reanalysis and reforecasts, while biases in soil moisture grow throughout the duration of the forecasts. Locally, the soil moisture biases may shrink or reverse sign. These biases are reflected in evaporation and runoff. The Noah land surface scheme shows the necessary relationships between evaporation and soil moisture for land-driven climate predictability. There is evidence that the atmospheric model cannot maintain the link between precipitation and antecedent soil moisture as strongly as in the real atmosphere, potentially hampering prediction skill, although there is better precipitation forecast skill over most locations when initial soil moisture anomalies are large. Bias change with lead-time, measured as the variance across ten monthly forecast leads, is often comparable to or larger than the interannual variance. Skill scores when forecast anomalies are calculated relative to reanalysis are seriously reduced over most locations when compared to validation against anomalies based on the forecast model climate at the corresponding lead-time. When all anomalies are calculated relative to the 0-month forecast, some skill is recovered over some regions, but the complex manner in which biases evolve indicates that a complete suite of reforecasts would be necessary whenever a new version of a climate model is implemented. The utility of reforecast programs is evident for operational forecast systems.  相似文献   

19.
Decadal prediction is one focus of the upcoming 5th IPCC Assessment report. To be able to interpret the results and to further improve the decadal predictions it is important to investigate the potential predictability in the participating climate models. This study analyzes the upper limit of climate predictability on decadal time scales and its dependency on sea ice albedo parameterization by performing two perfect ensemble experiments with the global coupled climate model EC-Earth. In the first experiment, the standard albedo formulation of EC-Earth is used, in the second experiment sea ice albedo is reduced. The potential prognostic predictability is analyzed for a set of oceanic and atmospheric parameters. The decadal predictability of the atmospheric circulation is small. The highest potential predictability was found in air temperature at 2?m height over the northern North Atlantic and the southern South Atlantic. Over land, only a few areas are significantly predictable. The predictability for continental size averages of air temperature is relatively good in all northern hemisphere regions. Sea ice thickness is highly predictable along the ice edges in the North Atlantic Arctic Sector. The meridional overturning circulation is highly predictable in both experiments and governs most of the decadal climate predictability in the northern hemisphere. The experiments using reduced sea ice albedo show some important differences like a generally higher predictability of atmospheric variables in the Arctic or higher predictability of air temperature in Europe. Furthermore, decadal variations are substantially smaller in the simulations with reduced ice albedo, which can be explained by reduced sea ice thickness in these simulations.  相似文献   

20.
Identifying regions sensitive to external radiative changes, including anthropogenic (sulphate aerosols and greenhouse gases) and natural (volcanoes and solar variations) forcings, is important to formulate actionable information at multi-year time-scales. Internally-generated climate variability can overcome this radiative forcing, especially at regional level, so that detecting the areas for this potential dominance is likewise critical for decadal prediction. This work aims to clarify where each contribution has the largest effect on North Atlantic sea surface temperature (SST) predictions in relation to the Atlantic multi-decadal variability (AMV). Initialized decadal hindcasts and radiatively-forced historical simulations from the fifth phase of the Climate Model Intercomparison Project are analysed to assess multi-year skill of the AMV. The initialized hindcasts reproduce better the phase of the AMV index fluctuations. The radiatively-forced component consists of a residual positive trend, although its identification is ambiguous. Initialization reduces the inter-model spread when estimating the level of AMV skill, thus reducing its uncertainty. Our results show a skilful performance of the initialized hindcasts in capturing the AMV-related SST anomalies over the subpolar gyre and Labrador Sea regions, as well as in the eastern subtropical basin, and the inability of the radiatively-forced historical runs to simulate the horseshoe-like AMV signature over the North Atlantic. Initialization outperforms empirical predictions based on persistence beyond 1–4 years ahead, suggesting that ocean dynamics play a role in the AMV predictability beyond the thermal inertia. The initialized hindcasts are also more skilful at reproducing the observed AMV teleconnection to the West African monsoon. The impact of the start date frequency is also described, showing that the standard of 5-year interval between start dates yields the main features of the AMV skill that are robustly detected in hindcasts with yearly start date sampling. This work updates previous studies, complementing them, and concludes that skilful initialized multi-model forecasts of the AMV-related climate variability can be formulated, improving uninitialized projections, until 3–6 years ahead.  相似文献   

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

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