首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The study examines future scenarios of precipitation extremes over Central Europe in an ensemble of 12 regional climate model (RCM) simulations with the 25-km resolution, carried out within the European project ENSEMBLES. We apply the region-of-influence method as a pooling scheme when estimating distributions of extremes, which consists in incorporating data from a ‘region’ (set of gridboxes) when fitting an extreme value distribution in any single gridbox. The method reduces random variations in the estimates of parameters of the extreme value distribution that result from large spatial variability of heavy precipitation. Although spatial patterns differ among the models, most RCMs simulate increases in high quantiles of precipitation amounts when averaged over the area for the late-twenty-first century (2070–2099) climate in both winter and summer. The sign as well as the magnitude of the projected change vary only little for individual parts of the distribution of daily precipitation in winter. In summer, on the other hand, the projected changes increase with the quantile of the distribution in all RCMs, and they are negative (positive) for parts of the distribution below (above) the 98% quantile if averaged over the RCMs. The increases in precipitation extremes in summer are projected in spite of a pronounced drying in most RCMs. Although a rather general qualitative agreement of the models concerning the projected changes of precipitation extremes is found in both winter and summer, the uncertainties in climate change scenarios remain large and would likely further increase considerably if a more complete ensemble of RCM simulations driven by a larger suite of global models and with a range of possible scenarios of the radiative forcing is available.  相似文献   

2.
This paper combines the climatological and societal perspectives for assessing future climatic extremes over Kangasabati River basin in India using an ensemble of four high resolution (25 km) regional climate model (RCM) simulations from 1970 to 2050. The relevant extreme indices and their thresholds are defined in consultation with stakeholders and are then compared using RCM simulations. To evaluate the performance of RCM in realistically representing atmospheric processes in the basin, model simulations driven with ERAInterim global re-analysis data from 1989 to 2008 are compared with observations. The models perform well in simulating seasonality, interannual variability and climatic extremes. Future climatic extremes are evaluated based on RCM simulations driven by GCMs, for present (1970–1999) and for the SRES A1B scenario for future (2021–2050) period. The analysis shows an intensification of majority of extremes as projected by future ensemble mean. The study suggests that there is a marked consistency in stakeholder observed changes in climate extremes and future predicted trends.  相似文献   

3.
The ability of a large ensemble of regional climate models to accurately simulate heat waves at the regional scale of Europe was evaluated. Within the EURO-CORDEX project, several state-of-the art models, including non-hydrostatic meso-scale models, were run for an extended time period (20 years) at high resolution (12 km), over a large domain allowing for the first time the simultaneous representation of atmospheric phenomena over a large range of spatial scales. Eight models were run in this configuration, and thirteen models were run at a classical resolution of 50 km. The models were driven with the same boundary conditions, the ERA-Interim re-analysis, and except for one simulation, no observations were assimilated in the inner domain. Results, which are compared with daily temperature and precipitation observations (ECA&D and E-OBS data sets) show that, even forced by the same re-analysis, the ensemble exhibits a large spread. A preliminary analysis of the sources of spread, using in particular simulations of the same model with different parameterizations, shows that the simulation of hot temperature is primarily sensitive to the convection and the microphysics schemes, which affect incoming energy and the Bowen ratio. Further, most models exhibit an overestimation of summertime temperature extremes in Mediterranean regions and an underestimation over Scandinavia. Even after bias removal, the simulated heat wave events were found to be too persistent, but a higher resolution reduced this deficiency. The amplitude of events as well as the variability beyond the 90th percentile threshold were found to be too strong in almost all simulations and increasing resolution did not generally improve this deficiency. Resolution increase was also shown to induce large-scale 90th percentile warming or cooling for some models, with beneficial or detrimental effects on the overall biases. Even though full causality cannot be established on the basis of this evaluation work, the drivers of such regional differences were shown to be linked to changes in precipitation due to resolution changes, affecting the energy partitioning. Finally, the inter-annual sequence of hot summers over central/southern Europe was found to be fairly well simulated in most experiments despite an overestimation of the number of hot days and of the variability. The accurate simulation of inter-annual variability for a few models is independent of the model bias. This indicates that internal variability of high summer temperatures should not play a major role in controlling inter-annual variability. Despite some improvements, especially along coastlines, the analyses conducted here did not allow us to generally conclude that a higher resolution is clearly beneficial for a correct representation of heat waves by regional climate models. Even though local-scale feedbacks should be better represented at high resolution, combinations of parameterizations have to be improved or adapted accordingly.  相似文献   

4.
Climate scenarios for the Netherlands are constructed by combining information from global and regional climate models employing a simplified, conceptual framework of three sources (levels) of uncertainty impacting on predictions of the local climate. In this framework, the first level of uncertainty is determined by the global radiation balance, resulting in a range of the projected changes in the global mean temperature. On the regional (1,000–5,000 km) scale, the response of the atmospheric circulation determines the second important level of uncertainty. The third level of uncertainty, acting mainly on a local scale of 10 (and less) to 1,000 km, is related to the small-scale processes, like for example those acting in atmospheric convection, clouds and atmospheric meso-scale circulations—processes that play an important role in extreme events which are highly relevant for society. Global climate models (GCMs) are the main tools to quantify the first two levels of uncertainty, while high resolution regional climate models (RCMs) are more suitable to quantify the third level. Along these lines, results of an ensemble of RCMs, driven by only two GCM boundaries and therefore spanning only a rather narrow range in future climate predictions, are rescaled to obtain a broader uncertainty range. The rescaling is done by first disentangling the climate change response in the RCM simulations into a part related to the circulation, and a residual part which is related to the global temperature rise. Second, these responses are rescaled using the range of the predictions of global temperature change and circulation change from five GCMs. These GCMs have been selected on their ability to simulate the present-day circulation, in particular over Europe. For the seasonal means, the rescaled RCM results obey the range in the GCM ensemble using a high and low emission scenario. Thus, the rescaled RCM results are consistent with the GCM results for the means, while adding information on the small scales and the extremes. The method can be interpreted as a combined statistical–dynamical downscaling approach, with the statistical relations based on regional model output.  相似文献   

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

6.
This study aims at exploring potential impacts of land-use vegetation change (LUC) on regional climate variability and extremes. Results from a pair of Australian Bureau of Meteorology Research Centre (BMRC) climate model 54-yr (1949-2002) integrations have been analysed. In the model experiments, two vegetation datasets are used, with one representing current vegetation coverage in China and the other approximating its potential coverage without human intervention. The model results show potential impacts ...  相似文献   

7.
The behaviour of precipitation and maximum temperature extremes in the Mediterranean area under climate change conditions is analysed in the present study. In this context, the ability of synoptic downscaling techniques in combination with extreme value statistics for dealing with extremes is investigated. Analyses are based upon a set of long-term station time series in the whole Mediterranean area. At first, a station-specific ensemble approach for model validation was developed which includes (1) the downscaling of daily precipitation and maximum temperature values from the large-scale atmospheric circulation via analogue method and (2) the fitting of extremes by generalized Pareto distribution (GPD). Model uncertainties are quantified as confidence intervals derived from the ensemble distributions of GPD-related return values and described by a new metric called “ratio of overlapping”. Model performance for extreme precipitation is highest in winter, whereas the best models for maximum temperature extremes are set up in autumn. Valid models are applied to a 30-year period at the end of the twenty-first century (2070–2099) by means of ECHAM5/MPI-OM general circulation model data for IPCC SRES B1 scenario. The most distinctive future changes are observed in autumn in terms of a strong reduction of precipitation extremes in Northwest Iberia and the Northern Central Mediterranean area as well as a simultaneous distinct increase of maximum temperature extremes in Southwestern Iberia and the Central and Southeastern Mediterranean regions. These signals are checked for changes in the underlying dynamical processes using extreme-related circulation classifications. The most important finding connected to future changes of precipitation extremes in the Northwestern Mediterranean area is a reduction of southerly displaced deep North Atlantic cyclones in 2070–2099 as associated with a strengthened North Atlantic Oscillation. Thus, the here estimated future changes of extreme precipitation are in line with the discourse about the influence of North Atlantic circulation variability on the changing climate in Europe.  相似文献   

8.
A new method is proposed to estimate future net basin supplies and lake levels for the Laurentian Great Lakes based on GCM projections of global climate change. The method first dynamically downscales the GCM simulation with a regional climate model, and then bias—corrects the simulated net basin supply in order to be used directly in a river—routing/lake level scheme. This technique addresses two weaknesses in the traditional approach, whereby observed sequences of climate variables are perturbed with fixed ratios or differences derived directly from GCMs in order to run evaporation and runoff models. Specifically, (1) land surface—atmosphere feedback processes are represented, and (2) changes in variability can be analyzed with the new approach. The method is demonstrated with a single, high resolution simulation, where small changes in future mean lake levels for all the upper Great Lakes are found, and an increase in seasonal range—especially for Lake Superior—is indicated. Analysis of a small ensemble of eight lower resolution regional climate model simulations supports these findings. In addition, a direct comparison with the traditional approach based on the same GCM projections used as the driving simulations in this ensemble shows that the new method indicates smaller declines in level for all the upper Great Lakes than has been reported previously based on the traditional method, though median differences are only a few centimetres in each case.  相似文献   

9.
It is increasingly accepted that any possible climate change will not only have an influence on mean climate but may also significantly alter climatic variability. A change in the distribution and magnitude of extreme rainfall events (associated with changing variability), such as droughts or flooding, may have a far greater impact on human and natural systems than a changing mean. This issue is of particular importance for environmentally vulnerable regions such as southern Africa. The sub-continent is considered especially vulnerable to and ill-equipped (in terms of adaptation) for extreme events, due to a number of factors including extensive poverty, famine, disease and political instability. Rainfall variability and the identification of rainfall extremes is a function of scale, so high spatial and temporal resolution data are preferred to identify extreme events and accurately predict future variability. The majority of previous climate model verification studies have compared model output with observational data at monthly timescales. In this research, the assessment of ability of a state of the art climate model to simulate climate at daily timescales is carried out using satellite-derived rainfall data from the Microwave Infrared Rainfall Algorithm (MIRA). This dataset covers the period from 1993 to 2002 and the whole of southern Africa at a spatial resolution of 0.1° longitude/latitude. This paper concentrates primarily on the ability of the model to simulate the spatial and temporal patterns of present-day rainfall variability over southern Africa and is not intended to discuss possible future changes in climate as these have been documented elsewhere. Simulations of current climate from the UK Meteorological Office Hadley Centre’s climate model, in both regional and global mode, are firstly compared to the MIRA dataset at daily timescales. Secondly, the ability of the model to reproduce daily rainfall extremes is assessed, again by a comparison with extremes from the MIRA dataset. The results suggest that the model reproduces the number and spatial distribution of rainfall extremes with some accuracy, but that mean rainfall and rainfall variability is under-estimated (over-estimated) over wet (dry) regions of southern Africa.  相似文献   

10.
This study presents a performance-based comprehensive weighting factor that accounts for the skill of different regional climate models (RCMs), including the effect of the driving lateral boundary condition coming from either atmosphere–ocean global climate models (AOGCMs) or reanalyses. A differential evolution algorithm is employed to identify the optimal relative importance of five performance metrics, and corresponding weighting factors, that include the relative absolute mean error (RAME), annual cycle, spatial pattern, extremes and multi-decadal trend. Based on cumulative density functions built by weighting factors of various RCMs/AOGCMs ensemble simulations, current and future climate projections were then generated to identify the level of uncertainty in the climate scenarios. This study selected the areas of southern Ontario and Québec in Canada as a case study. The main conclusions are as follows: (1) Three performance metrics were found essential, having the greater relative importance: the RAME, annual variability and multi-decadal trend. (2) The choice of driving conditions from the AOGCM had impacts on the comprehensive weighting factor, particularly for the winter season. (3) Combining climate projections based on the weighting factors significantly increased the consistency and reduced the spread among models in the future climate changes. These results imply that the weighting factors play a more important role in reducing the effects of outliers on plausible future climate conditions in regions where there is a higher level of variability in RCM/AOGCM simulations. As a result of weighting, substantial increases in the projected warming were found in the southern part of the study area during summer, and the whole region during winter, compared to the simple equal weighting scheme from RCM runs. This study is an initial step toward developing a likelihood procedure for climate scenarios on a regional scale using equal or different probabilities for all models.  相似文献   

11.
A number of uncertainties exist in climate simulation because the results of climate models are influenced by factors such as their dynamic framework, physical processes, initial and driving fields, and horizontal and vertical resolution. The uncertainties of the model results may be reduced, and the credibility can be improved by employing multi-model ensembles. In this paper, multi-model ensemble results using 10-year simulations of five regional climate models (RCMs) from December 1988 to November 1998 over Asia are presented and compared. The simulation results are derived from phase II of the Regional Climate Model Inter-comparison Project (RMIP) for Asia. Using the methods of the arithmetic mean, the weighted mean, multivariate linear regression, and singular value decomposition, the ensembles for temperature, precipitation, and sea level pressure are carried out. The results show that the multi-RCM ensembles outperform the single RCMs in many aspects. Among the four ensemble methods used, the multivariate linear regression, based on the minimization of the root mean square errors, significantly improved the ensemble results. With regard to the spatial distribution of the mean climate, the ensemble result for temperature was better than that for precipitation. With an increasing number of models used in the ensembles, the ensemble results were more accurate. Therefore, a multi-model ensemble is an efficient approach to improve the results of regional climate simulations.  相似文献   

12.
We investigate how well the variability of extreme daily precipitation events across the United Kingdom is represented in a set of regional climate models and the E-OBS gridded data set. Instead of simply evaluating the climatologies of extreme precipitation measures, we develop an approach to validate the representation of physical mechanisms controlling extreme precipitation variability. In part I of this study we applied a statistical model to investigate the influence of the synoptic scale atmospheric circulation on extreme precipitation using observational rain gauge data. More specifically, airflow strength, direction and vorticity are used as predictors for the parameters of the generalised extreme value (GEV) distribution of local precipitation extremes. Here we employ this statistical model for our validation study. In a first step, the statistical model is calibrated against a gridded precipitation data set provided by the UK Met Office. In a second step, the same statistical model is calibrated against 14 ERA40 driven 25?km resolution RCMs from the ENSEMBLES project and the E-OBS gridded data set. Validation indices describing relevant physical mechanisms are derived from the statistical models for observations and RCMs and are compared using pattern standard deviation, pattern correlation and centered pattern root mean squared error as validation measures. The results for the different RCMs and E-OBS are visualised using Taylor diagrams. We show that the RCMs adequately simulate moderately extreme precipitation and the influence of airflow strength and vorticity on precipitation extremes, but show deficits in representing the influence of airflow direction. Also very rare extremes are misrepresented, but this result is afflicted with a high uncertainty. E-OBS shows considerable biases, in particular in regions of sparse data. The proposed approach might be used to validate other physical relationships in regional as well as global climate models.  相似文献   

13.
Observations as well as most climate model simulations are generally in accord with the hypothesis that the hydrologic cycle should intensify and become highly volatile with the greenhouse-gas-induced climate change, although uncertainties of these projections as well as the spatial and seasonal variability of the changes are much larger than for temperature extremes. In this study, we examine scenarios of changes in extreme precipitation events in 24 future climate runs of ten regional climate models, focusing on a specific area of the Czech Republic (central Europe) where complex orography and an interaction of other factors governing the occurrence of heavy precipitation events result in patterns that cannot be captured by global models. The peaks-over-threshold analysis with increasing threshold censoring is applied to estimate multi-year return levels of daily rainfall amounts. Uncertainties in scenarios of changes for the late 21st century related to the inter-model and within-ensemble variability and the use of the SRES-A2 and SRES-B2 greenhouse gas emission scenarios are evaluated. The results show that heavy precipitation events are likely to increase in severity in winter and (with less agreement among models) also in summer. The inter-model and intra-model variability and related uncertainties in the pattern and magnitude of the change is large, but the scenarios tend to agree with precipitation trends recently observed in the area, which may strengthen their credibility. In most scenario runs, the projected change in extreme precipitation in summer is of the opposite sign than a change in mean seasonal totals, the latter pointing towards generally drier conditions in summer. A combination of enhanced heavy precipitation amounts and reduced water infiltration capabilities of a dry soil may severely increase peak river discharges and flood-related risks in this region.  相似文献   

14.
Using an ensemble of four high resolution (~25 km) regional climate models, this study analyses the future (2021–2050) spatial distribution of seasonal temperature and precipitation extremes in the Ganges river basin based on the SRES A1B emissions scenario. The model validation results (1989–2008) show that the models simulate seasonality and spatial distribution of extreme temperature events better than precipitation. The models are able to capture fine topographical detail in the spatial distribution of indices based on their ability to resolve processes at a higher regional resolution. Future simulations of extreme temperature indices generally agree with expected warming in the Ganges basin, with considerable seasonal and spatial variation. Significantly warmer summers in the central part of the basin along with basin-wide increase in night temperature are expected during the summer and monsoon months. An increase in heavy precipitation indices during monsoon, coupled with extended periods without precipitation during the winter months; indicates an increase in the incidence of extreme events.  相似文献   

15.
Going to the Extremes   总被引:8,自引:1,他引:8  
Projections of changes in climate extremes are critical to assessing the potential impacts of climate change on human and natural systems. Modeling advances now provide the opportunity of utilizing global general circulation models (GCMs) for projections of extreme temperature and precipitation indicators. We analyze historical and future simulations of ten such indicators as derived from an ensemble of 9 GCMs contributing to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC-AR4), under a range of emissions scenarios. Our focus is on the consensus from the GCM ensemble, in terms of direction and significance of the changes, at the global average and geographical scale. The climate extremes described by the ten indices range from heat-wave frequency to frost-day occurrence, from dry-spell length to heavy rainfall amounts. Historical trends generally agree with previous observational studies, providing a basic sense of reliability for the GCM simulations. Individual model projections for the 21st century across the three scenarios examined are in agreement in showing greater temperature extremes consistent with a warmer climate. For any specific temperature index, minor differences appear in the spatial distribution of the changes across models and across scenarios, while substantial differences appear in the relative magnitude of the trends under different emissions rates. Depictions of a wetter world and greater precipitation intensity emerge unequivocally in the global averages of most of the precipitation indices. However, consensus and significance are less strong when regional patterns are considered. This analysis provides a first overview of projected changes in climate extremes from the IPCC-AR4 model ensemble, and has significant implications with regard to climate projections for impact assessments. An erratum to this article is available at . An erratum to this article can be found at  相似文献   

16.
Assessments of the impacts of uncertainties in parameters on mean climate and climate change in complex climate models have, to date, largely focussed on perturbations to parameters in the atmosphere component of the model. Here we expand on a previously published study which found the global impacts of perturbed ocean parameters on the rate of transient climate change to be small compared to perturbed atmosphere parameters. By separating the climate-change-induced ocean vertical heat transport in each perturbed member into components associated with the resolved flow and each parameterisation scheme, we show that variations in global mean heat uptake in different perturbed versions are an order of magnitude smaller than the average heat uptake. The lack of impact of the perturbations is attributed to (1) the relatively small impact of the perturbation on the direct vertical heat transport associated with the perturbed process and (2) a compensation between those direct changes and indirect changes in heat transport from other processes. Interactions between processes and changes appear to combine in complex ways to limit ensemble spread and uncertainty in the rate of warming. We also investigate regional impacts of the perturbations that may be important for climate change predictions. We find variations across the ensemble that are significant when measured against natural variability. In terms of the experimental set-up used here (models without flux adjustments) we conclude that perturbed physics ensembles with ocean parameter perturbations are an important component of any probabilistic estimate of future climate change, despite the low spread in global mean quantities. Hence, careful consideration should be given to assessing uncertainty in ocean processes in future probabilistic assessments of regional climate change.  相似文献   

17.
Considered are the changes in the climate impact on the objects of construction and infrastructure on the territory of Russia. The focus is on the changes in the characteristics of daily air temperature and precipitation expected by the middle of the 21st century, which are of high importance in terms of the building design. The assessment of expected changes is based on the results of ensemble computations using the MGO gobal climate model and the embedded regional model (MGO RCM) with the horizontal resolution of 25 km. Along with the ensemble-averaged estimates of changes in applied climate parameters, an uncertainty of estimates associated with the natural climate variability is analyzed using the data of numerical experiments. The attention is drawn to the most significant effects of climate changes which should be taken into account when developing the measures for adapting the construction sector in Russia.  相似文献   

18.
In this study, an ensemble of four multi-year climate simulations is performed with the regional climate model ALADIN to evaluate its ability to simulate the climate over North America in the CORDEX framework. The simulations differ in their driving fields (ERA-40 or ERA-Interim) and the nudging technique (with or without large-scale nudging). The validation of the simulated 2-m temperature and precipitation with observationally-based gridded data sets shows that ALADIN performs similarly to other regional climate models that are commonly used over North America. Large-scale nudging improves the temporal correlation of the atmospheric circulation between ALADIN and its driving field, and also reduces the warm and dry summer biases in central North America. The differences between the simulations driven with different reanalyses are small and are likely related to the regional climate model’s induced internal variability. In general, the impact of different driving fields on ALADIN is smaller than that of large-scale nudging. The analysis of the multi-year simulations over the prairie and the east taiga indicates that the ALADIN 2-m temperature and precipitation interannual variability is similar or larger than that observed. Finally, a comparison of the simulations with observations for the summer 1993 shows that ALADIN underestimates the flood in central North America mainly due to its systematic dry bias in this region. Overall, the results indicate that ALADIN can produce a valuable contribution to CORDEX over North America.  相似文献   

19.
A high resolution regional climate model (RCM) is used to simulate climate of the recent past and to project future climate change across the northeastern US. Different types of uncertainties in climate simulations are examined by driving the RCM with different boundary data, applying different emissions scenarios, and running an ensemble of simulations with different initial conditions. Empirical orthogonal functions analysis and K-means clustering analysis are applied to divide the northeastern US region into four climatologically different zones based on the surface air temperature (SAT) and precipitation variability. The RCM simulations tend to overestimate SAT, especially over the northern part of the domain in winter and over the western part in summer. Statistically significant increases in seasonal SAT under both higher and lower emissions scenarios over the whole RCM domain suggest the robustness of future warming. Most parts of the northeastern US region will experience increasing winter precipitation and decreasing summer precipitation, though the changes are not statistically significant. The greater magnitude of the projected temperature increase by the end of the twenty-first century under the higher emissions scenario emphasizes the essential role of emissions choices in determining the potential future climate change.  相似文献   

20.
Validation results of the MGO regional climate model (RCM) with 50-km resolution for Siberia are discussed. For the specification of side boundary conditions, the reanalysis data are used. It is shown that the model satisfactorily simulates the sea-level pressure and temperature fields for all seasons and the year as a whole. The lowest computational errors in the simulation of regional surface temperature arise in the fall and winter; in spring and summer, the temperature errors are slightly higher. The model slightly underestimates the variability of daily mean temperature in winter relative to the reanalysis data. In summer, on the contrary, the RCM-simulated variability exceeds the variability in reanalysis. In winter, the space distribution of model precipitation is in qualitative agreement with the data of observational analysis; in summer, the space variability of model precipitation is significantly higher than that of precipitation in the reanalysis, especially in the mountains. Agreement between time changes in precipitation and temperature anomalies in RCM and in the reanalysis is better in the areas with a relatively large number of weather stations. The model can be used for estimation of future climate changes in the above-mentioned region.  相似文献   

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

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