首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The variability and predictability of the surface wind field at the regional scale is explored over a complex terrain region in the northeastern Iberian Peninsula by means of a downscaling technique based on Canonical Correlation Analysis. More than a decade of observations (1992–2005) allows for calibrating and validating a statistical method that elicits the main associations between the large scale atmospheric circulation over the North Atlantic and Mediterranean areas and the regional wind field. In an initial step the downscaling model is designed by selecting parameter values from practise. To a large extent, the variability of the wind at monthly timescales is found to be governed by the large scale circulation modulated by the particular orographic features of the area. The sensitivity of the downscaling methodology to the selection of the model parameter values is explored, in a second step, by performing a systematic sampling of the parameters space, avoiding a heuristic selection. This provides a metric for the uncertainty associated with the various possible model configurations. The uncertainties associated with the model configuration are considerably dependent on the spatial variability of the wind. While the sampling of the parameters space in the model set up moderately impact estimations during the calibration period, the regional wind variability is very sensitive to the parameters selection at longer timescales. This fact illustrates that downscaling exercises based on a single configuration of parameters should be interpreted with extreme caution. The downscaling model is used to extend the estimations several centuries to the past using long datasets of sea level pressure, thereby illustrating the large temporal variability of the regional wind field from interannual to multicentennial timescales. The analysis does not evidence long term trends throughout the twentieth century, however anomalous episodes of high/low wind speeds are identified.  相似文献   

2.
 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  相似文献   

3.
The key aspect of the ocean circulation off Peru?CChile is the wind-driven upwelling of deep, cold, nutrient-rich waters that promote a rich marine ecosystem. It has been suggested that global warming may be associated with an intensification of upwelling-favorable winds. However, the lack of high-resolution long-term observations has been a limitation for a quantitative analysis of this process. In this study, we use a statistical downscaling method to assess the regional impact of climate change on the sea-surface wind over the Peru?CChile upwelling region as simulated by the global coupled general circulation model IPSL-CM4. Taking advantage of the high-resolution QuikSCAT wind product and of the NCEP reanalysis data, a statistical model based on multiple linear regressions is built for the daily mean meridional and zonal wind at 10?m for the period 2000?C2008. The large-scale 10?m wind components and sea level pressure are used as regional circulation predictors. The skill of the downscaling method is assessed by comparing with the surface wind derived from the ERS satellite measurements, with in situ wind observations collected by ICOADS and through cross-validation. It is then applied to the outputs of the IPSL-CM4 model over stabilized periods of the pre-industrial, 2?×?CO2 and 4?×?CO2 IPCC climate scenarios. The results indicate that surface along-shore winds off central Chile (off central Peru) experience a significant intensification (weakening) during Austral winter (summer) in warmer climates. This is associated with a general decrease in intra-seasonal variability.  相似文献   

4.
Statistical downscaling of 14 coupled atmosphere-ocean general circulation models (AOGCM) is presented to assess potential changes of the 10 m wind speeds in France. First, a statistical downscaling method is introduced to estimate daily mean 10 m wind speed at specific sites using general circulation model output. Daily 850 hPa wind field has been selected as the large scale circulation predictor. The method is based on a classification of the daily wind fields into a few synoptic weather types and multiple linear regressions. Years are divided into an extended winter season from October to March and an extended summer season from April to September, and the procedure is conducted separately for each season. ERA40 reanalysis and observed station data have been used to build and validate the downscaling algorithm over France for the period 1974–2002. The method is then applied to 14 AOGCMs of the coupled model intercomparison project phase 3 (CMIP3) multi-model dataset. Three time periods are focused on: a historical period (1971–2000) from the climate of the twentieth century experiment and two climate projection periods (2046–2065 and 2081–2100) from the IPCC SRES A1B experiment. Evolution of the 10 m wind speed in France and associated uncertainties are discussed. Significant changes are depicted, in particular a decrease of the wind speed in the Mediterranean area. Sources of those changes are investigated by quantifying the effects of changes in the weather type occurrences, and modifications of the distribution of the days within the weather types.  相似文献   

5.
High-resolution sea wind hindcasts over the Mediterranean area   总被引：1，自引：1，他引：0  
The goal of this study is to develop a high-resolution atmospheric hindcast over the Mediterranean area using the WRF-ARW model, focusing on offshore surface wind fields. In order to choose the most adequate model configuration, the study provides details on the calibration of the experimental saet-up through a sensitivity test considering the October–December 2001 period (the 2001 super-storm event in the West Mediterranean). A daily forecast outperforms the spectral technique of previous products and the boundary data from ERA-Interim reanalysis produces the most accurate estimates in terms of wind variability and hour-to-hour correspondence. According to the sensitivity test, two data sets of wind hindcast are produced: the SeaWind I (30-km horizontal resolution for a period of 60 years) and the SeaWind II (15-km horizontal resolution for 20 years). The validation of the resulting surface winds is undertaken considering two offshore observational datasets. On the one hand, hourly surface buoy stations are used to validate wind time series at specific locations; on the other hand, wind altimeter satellite observations are considered for spatial validation in the whole Mediterranean Sea. The results obtained from this validation process show a very good agreement with observations for the southern Europe region. Finally, SeaWind I and II are used to characterize offshore wind fields in the Mediterranean Sea. The statistical structure of sea surface wind is analyzed and the agreement with Weibull probability distribution is discussed. In addition, wind persistence and extreme wind speed (50 year return period) are characterized and relevant areas of wind power generation are described by estimating wind energy quantities.  相似文献   

6.
The study examines two methods for classification of daily pressure-pattern sequences (extended PCA and extended cluster analysis) and evaluates whether sequence classification is more suitable to describe surface air-temperature conditions in Europe. For this purpose, sequences of daily sea-level pressure fields are classified for each month of the year using a wide range of class numbers (2–40) and sequence lengths from 1 to 6 days. In each case, the classification is used to reconstruct temperature variability in 5°×5° grid-boxes throughout Europe for the period 1850–2003 in order to compare the circulation type based reconstruction to monthly observations for skill estimation. Results show that extended PCA leads to more dynamic pattern sequences for subordinated classes, but also that it is limited in skill, because classes of very low frequencies are generated. In contrast, extended cluster analysis offers higher skill levels for temperature downscaling. However, both methods can benefit from using sequences instead of single days for classification, especially for continental regions.  相似文献   

7.
This study evaluates how statistical and dynamical downscaling models as well as combined approach perform in retrieving the space–time variability of near-surface temperature and rainfall, as well as their extremes, over the whole Mediterranean region. The dynamical downscaling model used in this study is the Weather Research and Forecasting (WRF) model with varying land-surface models and resolutions (20 and 50 km) and the statistical tool is the Cumulative Distribution Function-transform (CDF-t). To achieve a spatially resolved downscaling over the Mediterranean basin, the European Climate Assessment and Dataset (ECA&D) gridded dataset is used for calibration and evaluation of the downscaling models. In the frame of HyMeX and MED-CORDEX international programs, the downscaling is performed on ERA-I reanalysis over the 1989–2008 period. The results show that despite local calibration, CDF-t produces more accurate spatial variability of near-surface temperature and rainfall with respect to ECA&D than WRF which solves the three-dimensional equation of conservation. This first suggests that at 20–50 km resolutions, these three-dimensional processes only weakly contribute to the local value of temperature and precipitation with respect to local one-dimensional processes. Calibration of CDF-t at each individual grid point is thus sufficient to reproduce accurately the spatial pattern. A second explanation is the use of gridded data such as ECA&D which smoothes in part the horizontal variability after data interpolation and damps the added value of dynamical downscaling. This explains partly the absence of added-value of the 2-stage downscaling approach which combines statistical and dynamical downscaling models. The temporal variability of statistically downscaled temperature and rainfall is finally strongly driven by the temporal variability of its forcing (here ERA-Interim or WRF simulations). CDF-t is thus efficient as a bias correction tool but does not show any added-value regarding the time variability of the downscaled field. Finally, the quality of the reference observation dataset is a key issue. Comparison of CDF-t calibrated with ECA&D dataset and WRF simulations to local measurements from weather stations not assimilated in ECA&D, shows that the temporal variability of the downscaled data with respect to the local observations is closer to the local measurements than to ECA&D data. This highlights the strong added-value of dynamical downscaling which improves the temporal variability of the atmospheric dynamics with regard to the driving model. This article highlights the benefits and inconveniences emerging from the use of both downscaling techniques for climate research. Our goal is to contribute to the discussion on the use of downscaling tools to assess the impact of climate change on regional scales.  相似文献   

8.
Statistical Downscaling of Wind Variability from Meteorological Fields   总被引：1，自引：0，他引：1  
Measurements show that on numerous occasions the low-level wind is highly variable across a large portion of south-eastern Australia. Under such conditions the risk of a large rapid change in total wind power is increased. While variability tends to increase with mean wind speed, a large component of wind variability is not explained by wind speed alone. In this work, reanalysis fields from the US National Centers for Environmental Prediction (NCEP) are statistically downscaled to model wind variability at a coastal location in Victoria, Australia. In order to reduce the dimensionality of the problem, the NCEP fields are each decomposed using empirical orthogonal function (EOF) techniques. The downscaling technique is applied to two periods in the seasonal cycle, namely (i) winter to early spring, and (ii) summer. In each case, data representing 2 years are used to form a model that is then validated using independent data from another year. The EOFs that best predict wind variability are examined. To allow for non-linearity and complex interaction between variables, all empirical models are built using random forests. Quantitatively, the model compares favourably with a simple regression of wind variability against wind speed, as well as multiple linear regression models.  相似文献   

9.
Summary Regional climate model and statistical downscaling procedures are used to generate winter precipitation changes over Romania for the period 2071–2100 (compared to 1961–1990), under the IPCC A2 and B2 emission scenarios. For this purpose, the ICTP regional climate model RegCM is nested within the Hadley Centre global atmospheric model HadAM3H. The statistical downscaling method is based on the use of canonical correlation analysis (CCA) to construct climate change scenarios for winter precipitation over Romania from two predictors, sea level pressure and specific humidity (either used individually or together). A technique to select the most skillful model separately for each station is proposed to optimise the statistical downscaling signal. Climate fields from the A2 and B2 scenario simulations with the HadAM3H and RegCM models are used as input to the statistical downscaling model. First, the capability of the climate models to reproduce the observed link between winter precipitation over Romania and atmospheric circulation at the European scale is analysed, showing that the RegCM is more accurate than HadAM3H in the simulation of Romanian precipitation variability and its connection with large-scale circulations. Both models overestimate winter precipitation in the eastern regions of Romania due to an overestimation of the intensity and frequency of cyclonic systems over Europe. Climate changes derived directly from the RegCM and HadAM3H show an increase of precipitation during the 2071–2100 period compared to 1961–1990, especially over northwest and northeast Romania. Similar climate change patterns are obtained through the statistical downscaling method when the technique of optimum model selected separately for each station is used. This adds confidence to the simulated climate change signal over this region. The uncertainty of results is higher for the eastern and southeastern regions of Romania due to the lower HadAM3H and RegCM performance in simulating winter precipitation variability there as well as the reduced skill of the statistical downscaling model.  相似文献   

10.
About 75 % of the Antarctic surface mass gain occurs over areas below 2,000 m asl, which cover 40 % of the grounded ice-sheet. As the topography is complex in many of these regions, surface mass balance modelling is highly dependent on horizontal resolution, and studying the impact of Antarctica on the future rise in sea level requires physical approaches. We have developed a computationally efficient, physical downscaling model for high-resolution (15 km) long-term surface mass balance (SMB) projections. Here, we present results of this model, called SMHiL (surface mass balance high-resolution downscaling), which was forced with the LMDZ4 atmospheric general circulation model to assess Antarctic SMB variability in the twenty first and the twenty second centuries under two different scenarios. The higher resolution of SMHiL better reproduces the geographical patterns of SMB and increase significantly the averaged SMB over the grounded ice-sheet for the end of the twentieth century. A comparison with more than 3200 quality-controlled field data shows that LMDZ4 and SMHiL reproduce the observed values equally well. Nevertheless, field data below 2,000 m asl are too scarce to efficiently show the added value of SMHiL and measuring the SMB in these undocumented areas should be a future scientific priority. Our results suggest that running LMDZ4 at a finer resolution (15 km) may give a future increase in SMB in Antarctica that is about 30 % higher than by using its standard resolution (60 km) due to the higher increase in precipitation in coastal areas at 15 km. However, a part (～15 %) of these discrepancies could be an artefact from SMHiL since it neglects the foehn effect and likely overestimates the precipitation increase. Future changes in the Antarctic SMB at low elevations will result from the competition between higher snow accumulation and runoff. For this reason, developing downscaling models is crucial to represent processes in sufficient detail and correctly model the SMB in coastal areas.  相似文献   

11.
Freshwater (FW) leaves the Arctic Ocean through sea-ice export and the outflow of low-salinity upper ocean water. Whereas the variability of the sea-ice export is known to be mainly caused by changes in the local wind and the thickness of the exported sea ice, the mechanisms that regulate the variability of the liquid FW export are still under investigation. To better understand these mechanisms, we present an analysis of the variability of the liquid FW export from the Arctic Ocean for the period 1950–2007, using a simulation from an energy and mass conserving global ocean–sea ice model, coupled to an Energy Moisture Balance Model of the atmosphere, and forced with daily winds from the NCEP reanalysis. Our results show that the simulated liquid FW exports through the Canadian Arctic Archipelago (CAA) and the Fram Strait lag changes in the large-scale atmospheric circulation over the Arctic by 1 and 6 years, respectively. The variability of the liquid FW exports is caused by changes in the cyclonicity of the atmospheric forcing, which cause a FW redistribution in the Arctic through changes in Ekman transport in the Beaufort Gyre. This in turn causes changes in the sea surface height (SSH) and salinity upstream of the CAA and Fram Strait, which affect the velocity and salinity of the outflow. The SSH changes induced by the large-scale atmospheric circulation are found to explain a large part of the variance of the liquid FW export, while the local wind plays a much smaller role. We also show that during periods of increased liquid FW export from the Arctic, the strength of the simulated Atlantic meridional overturning circulation is reduced and the ocean heat transport into the Arctic is increased. These results are particularly relevant in the context of global warming, as climate simulations predict an increase in the liquid FW export from the Arctic during the twenty-first century.  相似文献   

12.
This study documents the temporal and spatial variability of surface wind conditions over the Norwegian county of Finnmark and the coupling of local surface winds to the larger-scale atmospheric circulation, represented by the mean sea level pressure field. At locations along the northern coast, thermally driven offshore winds from the south dominate, especially during the cold season. During the warm season, downward mixing of westerly overlying winds becomes more important as the stability of the boundary layer stratification decreases. In the western part, locations are situated in valleys, resulting in two opposing along-channel dominant wind directions. Seasonal changes in atmospheric conditions are reflected in a reversal of the dominant wind direction, with a component towards the coast during the cold season. At all locations, a clear separation between different prevailing surface wind directions in each season can be achieved based exclusively on local mean sea level geostrophic wind direction. This allows statistical downscaling of the prevailing surface wind conditions from lower-resolution simulations of the surface pressure field and may improve local wind forecasts over complex terrain.  相似文献   

13.
黑河流域大气资料尺度转换的对比分析   总被引：8，自引：2，他引：6  
尺度转换是解决尺度差异问题的有效方法之一。尺度转换统计模式的可靠性和正确性如何,在很大程度上取决于建立统计模式时所用的实际气象资料的处理,尤其在下垫面分布极不均匀地区。本文研究了黑河流域这一复杂下垫面条件下的大气资料的尺度转换问题,对比分析了由NCEP再分析资料直接内插到较小尺度而得到的局地气候变化与由尺度转换统计模式,经过客观分析所得到的局地气候变化。结果表明,直接插值所得的细网格值只能用于对平坦地区的温、湿状况的描述,而对较高山区的区域特征描述不够;进行观测资料的客观分析可以更准确地反映黑河流域复杂下垫面背景下近地层大气要素场的变化。研究表明,客观分析过程是尺度转换过程的重要环节。  相似文献   

14.
We study the influence of synoptic scale atmospheric circulation on extreme daily precipitation across the United Kingdom, using observed time series from 689 rain gauges. To this end we employ a statistical model, that uses airflow strength, direction and vorticity as predictors for the generalised extreme value distribution of monthly precipitation maxima. The inferred relationships are connected with the dominant westerly flow, the orography, and the moisture supply from surrounding seas. We aggregated the results for individual rain gauges to regional scales to investigate the temporal variability of extreme precipitation. Airflow explains a significant fraction of the variability on subannual to decadal time scales. A large fraction of the especially heavy winter precipitation during the 1980s and 1990s in north Scotland can be attributed to a prevailing positive phase of the North Atlantic Oscillation. Our statistical model can be used for statistical downscaling and to validate regional climate model output.  相似文献   

15.
The present and twenty-first century near-surface wind climate of Greenland is presented using output from the regional atmospheric climate model RACMO2. The modelled wind variability and wind distribution compare favourably to observations from three automatic weather stations in the ablation zone of southwest Greenland. The Weibull shape parameter is used to classify the wind climate. High values (κ > 4) are found in northern Greenland, indicative of uniform winds and a dominant katabatic forcing, while lower values (κ < 3) are found over the ocean and southern Greenland, where the synoptic forcing dominates. Very high values of the shape parameter are found over concave topography where confluence strengthens the katabatic circulation, while very low values are found in a narrow band along the coast due to barrier winds. To simulate the future (2081–2098) wind climate RACMO2 was forced with the HadGEM2-ES general circulation model using a scenario of mid-range radiative forcing of +4.5 W m?2 by 2100. For the future simulated climate, the near-surface potential temperature deficit reduces in all seasons in regions where the surface temperature is below the freezing point, indicating a reduction in strength of the near-surface temperature inversion layer. This leads to a wind speed reduction over the central ice sheet where katabatic forcing dominates, and a wind speed increase over steep coastal topography due to counteracting effects of thermal and katabatic forcing. Thermally forced winds over the seasonally sea ice covered region of the Greenland Sea are reduced by up to 2.5 m s?1.  相似文献   

16.
This work uses a WRF numerical simulation from 1960 to 2005 performed at a high horizontal resolution (2 km) to analyze the surface wind variability over a complex terrain region located in northern Iberia. A shorter slice of this simulation has been used in a previous study to demonstrate the ability of the WRF model in reproducing the observed wind variability during the period 1992–2005. Learning from that validation exercise, the extended simulation is herein used to inspect the wind behavior where and when observations are not available and to determine the main synoptic mechanisms responsible for the surface wind variability. A principal component analysis was applied to the daily mean wind. Two principal modes of variation accumulate a large percentage of the wind variability (83.7%). The first mode reflects the channeling of the flow between the large mountain systems in northern Iberia modulated by the smaller topographic features of the region. The second mode further contributes to stress the differentiated wind behavior over the mountains and valleys. Both modes show significant contributions at the higher frequencies during the whole analyzed period, with different contributions at lower frequencies during the different decades. A strong relationship was found between these two modes and the zonal and meridional large scale pressure gradients over the area. This relationship is described in the context of the influence of standard circulation modes relevant in the European region like the North Atlantic Oscillation, the East Atlantic pattern, East Atlantic/Western Russia pattern, and the Scandinavian pattern.  相似文献   

17.
马文通  朱蓉  李泽椿  龚玺 《气象学报》2016,74(1):89-102
复杂地形导致近地层风场时空变化大,是影响风电场短期风电功率预测准确率的重要因素。为此,基于中尺度数值预报模式和微尺度计算流体力学模式,建立了风电场短期风电功率动力降尺度预测系统。该系统由中尺度数值预报模式、微尺度风场基础数据库、风电功率预测集成系统组成,能够预测复杂地形风电场中每台风电机组未来72 h逐15 min的发电量。提高了复杂地形风场发电功率预测准确率,同时还可以在上报电网的风电功率预测结果中考虑运行维护计划和限电等因素对实际并网功率的影响。2014年7月-2015年1月的业务预测试验表明,风电场短期风电功率动力降尺度预测系统的月预测相对误差均小于0.2,满足中国国家电网对风电功率预测误差和时效性的业务要求。动力降尺度技术不受具体项目地形复杂程度和历史观测数据样本量的限制,可以在新建风电场中推广应用,具备实际的可操作性。   相似文献   

18.
Observations show a multidecadal signal in the North Atlantic ocean, but the underlying mechanism and cause of its timescale remain unknown. Previous studies have suggested that it may be driven by the North Atlantic Oscillation (NAO), which is the dominant pattern of winter atmospheric variability. To further address this issue, the global ocean general circulation model, Nucleus for European Modelling of the Ocean (NEMO), is driven using a 2,000 years long white noise forcing associated with the NAO. Focusing on key ocean circulation patterns, we show that the Atlantic Meridional Overturning Circulation (AMOC) and Sub-polar gyre (SPG) strength both have enhanced power at low frequencies but no dominant timescale, and thus provide no evidence for a oscillatory ocean-only mode of variability. Instead, both indices respond linearly to the NAO forcing, but with different response times. The variability of the AMOC at 30°N is strongly enhanced on timescales longer than 90 years, while that of the SPG strength starts increasing at 15 years. The different response characteristics are confirmed by constructing simple statistical models that show AMOC and SPG variability can be related to the NAO variability of the previous 53 and 10 winters, respectively. Alternatively, the AMOC and the SPG strength can be reconstructed with Auto-regressive (AR) models of order seven and five, respectively. Both statistical models reconstruct interannual and multidecadal AMOC variability well, while on the other hand, the AR(5) reconstruction of the SPG strength only captures multidecadal variability. Using these methods to reconstruct ocean variables can be useful for prediction and model intercomparision.  相似文献   

19.
Local to regional climate anomalies are to a large extent determined by the state of the atmospheric circulation. The knowledge of large-scale sea level pressure (SLP) variations in former times is therefore crucial when addressing past climate changes across Europe and the Mediterranean. However, currently available SLP reconstructions lack data from the ocean, particularly in the pre-1850 period. Here we present a new statistically-derived 5° × 5° resolved gridded seasonal SLP dataset covering the eastern North Atlantic, Europe and the Mediterranean area (40°W–50°E; 20°N–70°N) back to 1750 using terrestrial instrumental pressure series and marine wind information from ship logbooks. For the period 1750–1850, the new SLP reconstruction provides a more accurate representation of the strength of the winter westerlies as well as the location and variability of the Azores High than currently available multiproxy pressure field reconstructions. These findings strongly support the potential of ship logbooks as an important source to determine past circulation variations especially for the pre-1850 period. This new dataset can be further used for dynamical studies relating large-scale atmospheric circulation to temperature and precipitation variability over the Mediterranean and Eurasia, for the comparison with outputs from GCMs as well as for detection and attribution studies.  相似文献   

20.

Drought over the Greek region is characterized by a strong seasonal cycle and large spatial variability. Dry spells longer than 10 consecutive days mainly characterize the duration and the intensity of Greek drought. Moreover, an increasing trend of the frequency of drought episodes has been observed, especially during the last 20 years of the 20th century. Moreover, the most recent regional circulation models (RCMs) present discrepancies compared to observed precipitation, while they are able to reproduce the main patterns of atmospheric circulation. In this study, both a statistical and a dynamical downscaling approach are used to quantify drought episodes over Greece by simulating the Standardized Precipitation Index (SPI) for different time steps (3, 6, and 12 months). A statistical downscaling technique based on artificial neural network is employed for the estimation of SPI over Greece, while this drought index is also estimated using the RCM precipitation for the time period of 1961–1990. Overall, it was found that the drought characteristics (intensity, duration, and spatial extent) were well reproduced by the regional climate models for long term drought indices (SPI12) while ANN simulations are better for the short-term drought indices (SPI3).

  相似文献   

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

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