首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 234 毫秒
1.
The South Atlantic Convergence Zone (SACZ) is an intrinsic characteristic of the South American Summer Monsoon. In a recent study, we verified that the main mode of coupled variability over the South Atlantic (South Atlantic Dipole (SAD)) plays a role in modulating the position of extratropical cyclones that affect the SACZ precipitation. In this study, we perform numerical experiments to further investigate the mechanisms between SAD and the SACZ. Numerical experiments forced with prescribed SST anomalies showed that, even though the Atlantic SST affects the position of the cyclone associated with the SACZ, the atmospheric response and precipitation patterns over land are opposed to the observations. On the other hand, experiments forced with prescribed anomalous driving fields showed that the atmospheric component of SAD plays a significant role for the right position and intensity of precipitation associated with the SACZ. SAD negative anomalies provide the low-level and upper-level atmospheric support for the intensification of the cyclone at surface and for the increase in precipitation over the land portion of the SACZ. Therefore, the numerical experiments suggest that, during El Niño Southern Oscillation neutral conditions, the SACZ precipitation variability associated with SAD is largely dependent on the atmospheric variability rather than the underlying SST.  相似文献   

2.
This paper examines moisture transport on intraseasonal timescales over the continent and over the South Atlantic convergence zone (SACZ) during the South America (SA) summer monsoon. Combined Empirical Orthogonal Function analysis (EOFc) of Global Precipitation Climatology Project pentad precipitation, specific humidity, air temperature, zonal and meridional winds at 850?hPa (NCEP/NCAR reanalysis) are performed to identify the large-scale variability of the South America monsoon system and the SACZ. The first EOFc was used as a large-scale index for the South American monsoon (LISAM), whereas the second EOFc characterized the SACZ. LISAM (SACZ) index showed spectral variance on 30?C90 (15?C20) days and were both band filtered (10?C100?days). Intraseasonal wet anomalies were defined when LISAM and SACZ anomalies were above the 75th percentile of their respective distribution. LISAM and SACZ wet events were examined independently of each other and when they occur simultaneously. LISAM wet events were observed with the amplification of wave activity in the Northern Hemisphere and the enhancement of northwesterly cross-equatorial moisture transport over tropical continental SA. Enhanced SACZ was observed with moisture transport from the extratropics of the Southern Hemisphere. Simultaneous LISAM and SACZ wet events are associated with cross-equatorial moisture transport along with moisture transport from Subtropical Southwestern Atlantic.  相似文献   

3.
The regional influence of the Madden–Julian oscillation (MJO) on South America is described. Maps of probability of weekly-averaged rainfall exceeding the upper tercile were computed for all seasons and related statistically with the phase of the MJO as characterized by the Wheeler–Hendon real-time multivariate MJO (RMM) index and with the OLR MJO Index. The accompanying surface air temperature and circulation anomalies were also calculated. The influence of the MJO on regional scales along with their marked seasonal variations was documented. During December–February when the South American monsoon system is active, chances of enhanced rainfall are observed in southeastern South America (SESA) region mainly during RMM phases 3 and 4, accompanied by cold anomalies in the extratropics, while enhanced rainfall in the South Atlantic Convergence Zone (SACZ) region is observed in phases 8 and 1. The SESA (SACZ) signal is characterized by upper-level convergence (divergence) over tropical South America and a cyclonic (anticyclonic) anomaly near the southern tip of the continent. Impacts during March–May are similar, but attenuated in the extratropics. Conversely, in June–November, reduced rainfall and cold anomalies are observed near the coast of the SACZ region during phases 4 and 5, favored by upper-level convergence over tropical South America and an anticyclonic anomaly over southern South America. In September–November, enhanced rainfall and upper-level divergence are observed in the SACZ region during phases 7 and 8. These signals are generated primarily through the propagation of Rossby wave energy generated in the region of anomalous heating associated with the MJO.  相似文献   

4.
Summary  In subtropical Argentina, Paraguay and southern Brazil, precipitation is most abundant during summer but its interannual variability is large. At this time a zone of low-level convergence, upper-level divergence, and intense convection is developed to the north of this area. This feature is known as the South Atlantic convergence zone (SACZ) and seems to be related to the interannual variability of summer rainfall to its south. The aim of this work is to document this relationship. Reduced (increased) precipitation in southern Brazil, most of Uruguay and northeastern Argentina is associated with a strong (weak) SACZ and a northward (southward) displacement of it, while increased (reduced) rainfall occurs further south in subtropical Argentina. Also, warm (cold) SST in the region 20° S–40° S and west of 30° W is likely accompanied by a southward (northward) shift of the SACZ. Aside of this relation with the SACZ that affect on the precipitation field of Southeastern South America, the proximate Atlantic Ocean SST seems to force the precipitation over this region by other mechanisms as well. The result of this additional SST forcing is to enhance the signal of the SACZ in northeastern Argentina, Uruguay and southern Brazil and to oppose the SACZ effect in southern subtropical Argentina. Received July 24, 1999 Revised July 5, 2000  相似文献   

5.
S. B. Cerne  C. S. Vera 《Climate Dynamics》2011,36(11-12):2265-2277
The influence of the intraseasonal variability on heat wave development over subtropical South America during austral summer is analyzed. The role of the South Atlantic Convergence Zone (SACZ) on this development is documented. Results show that intraseasonal variability can explain on average at least 32% of summer temperature variance. Moreover, 73% of the heat waves in subtropical South America develop in association with an active SACZ. The analysis of pentad maps shows that warm conditions in the region under study develop in association with the strengthening of an anticyclonic anomaly, which is discernible over the subtropical regions at least 15?days before temperature peak occurrence. That circulation anomaly is embedded in a large-scale Rossby wave train extending along the South Pacific Ocean that is linked to convection anomalies at the equatorial western and central Pacific Ocean. In addition, the development of the anticyclonic circulation over subtropical South America appears to be strengthened by the subsidence conditions promoted by the active SACZ, which result in temperature rise in the subtropical region under relatively dry conditions. On the other hand, during the last 2?days of evolution, SACZ activity weakens and the progressive temperature rise in the region is dominated by warmer and moister air being anomalously advected from the north. Results confirm the important role that SACZ activity on intraseasonal time scales has in inducing persistent circulation anomalies at the subtropical regions that can result in the development of persistent heat waves, and very extreme daily temperature.  相似文献   

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

7.
This study explores the ocean–atmosphere interaction in the formation and dynamics of the South Atlantic Convergence Zone (SACZ), through the analysis of the heat sources estimated through the outgoing longwave radiation. The results obtained with this study show that the coupled variability between SACZ and the South Atlantic Ocean indicates that in northern positioned SACZ cases (over Southeastern Brazil), westerly anomalies are verified in the low level continental tropical circulation, consistent with the active phase of the South America Monsoon System (SAMS). In these cases, cold anomalies in the subtropical Atlantic Ocean cause an increase in the continent–ocean temperature gradient, favoring an easterly flow in this region, and blocking the SACZ at a northerly position. Easterly anomalies in the tropical continent were verified in the low level circulation in southern positioned cases (over Southern Brazil), consistent with the SAMS break phase. The SST anomaly patterns indicate cold anomalies in the tropics and warm anomalies in the subtropics, which do not favor the development of an easterly flow at low levels over the western tropical Atlantic. In these cases, two situations may occur: the strengthening of the Low Level Jet (LLJ), which prevails in the eastern subtropical South America and convergence with the South Atlantic Subtropical High at its southern position; or the atmospheric unstable conditions caused by ocean warm SST anomalies (in this case the LLJ may be weaker than its climatological intensity).  相似文献   

8.
We perform a systematic study of the predictability of surface air temperature and precipitation in Southeastern South America (SESA) using ensembles of AGCM simulations, focusing on the role of the South Atlantic and its interaction with the El Niño-Southern Oscillation (ENSO). It is found that the interannual predictability of climate over SESA is strongly tied to ENSO showing high predictability during the seasons and periods when there is ENSO influence. The most robust ENSO signal during the whole period of study (1949–2006) is during spring when warm events tend to increase the precipitation over Southeastern South America. Moreover, the predictability shows large inter-decadal changes: for the period 1949–1977, the surface temperature shows high predictability during late fall and early winter. On the other hand, for the period 1978–2006, the temperature shows (low) predictability only during winter, while the precipitation shows not only high predictability in spring but also in fall. Furthermore, it is found that the Atlantic does not directly affect the climate over SESA. However, the experiments where air–sea coupling is allowed in the south Atlantic suggest that this ocean can act as a moderator of the ENSO influence. During warm ENSO events the ocean off Brazil and Uruguay tends to warm up through changes in the atmospheric heat fluxes, altering the atmospheric anomalies and the predictability of climate over SESA. The main effect of the air–sea coupling is to strengthen the surface temperature anomalies over SESA; changes in precipitation are more subtle. We further found that the thermodynamic coupling can increase or decrease the predictability. For example, the air–sea coupling significantly increases the skill of the model in simulating the surface air temperature anomalies for most seasons during period 1949–1977, but tends to decrease the skill in late fall during period 1978–2006. This decrease in skill during late fall in 1978–2006 is found to be due to a wrong simulation of the remote ENSO signal that is further intensified by the local air–sea coupling in the south Atlantic. Thus, our results suggest that climate models used for seasonal prediction should simulate correctly not only the remote ENSO signal, but also the local air–sea thermodynamic coupling.  相似文献   

9.
Summary ?This paper presents an objective analysis of the structure of daily rainfall variability over the South American/South Atlantic region (15°–60° W and 0°–40° S) during individual austral summer months of November to March. From EOF analysis of satellite derived daily rainfall we find that the leading mode of variability is represented by a highly coherent meridional dipole structure, organised into 2 extensive bands, oriented northwest to southeast across the continent and Atlantic Ocean. We argue that this dipole structure represents variability in the meridional position of the South Atlantic Convergence Zone (SACZ). During early and later summer, in the positive (negative) phase of the dipole, enhanced (suppressed) rainfall over eastern tropical Brazil links with that over the subtropical and extra-tropical Atlantic and is associated with suppressed (enhanced) rainfall over the sub-tropical plains and adjacent Atlantic Ocean. This structure is indicative of interaction between the tropical, subtropical and temperate zones. Composite fields from NCEP reanalysis products (associated with the major positive and negative events) show that in early and late summer the position of the SACZ is associated with variability in: (a) the midlatitude wave structure, (b) the position of the continental low, and (c) the zonal position of the South Atlantic Subtropical High. Harmonic analysis of the 200 hPa geopotential anomaly structure in the midlatitudes indicates that reversals in the rainfall dipole structure are associated primarily with variability in zonal wave 4. There is evidence of a wave train extending throughout the midlatitudes from the western Pacific into the SACZ region. During positive (negative) events the largest anomalous moisture advection occurs within westerlies (easterlies) primarily from Amazonia (the South Atlantic). In both phases a convergent poleward flow results along the leading edge of the low-level trough extending from the tropics into temperate latitudes. High summer events differ from those in early and late summer in that the rainfall dipole is primarily associated with variability in the phase of zonal wave 3, and that tropical-temperate link is not clearly evident in positive events. Received May 31, 2001; revised October 17, 2001; accepted June 13, 2002  相似文献   

10.
The interdecadal change in seasonal predictability and numerical models’ seasonal forecast skill in the Northern Hemisphere are examined using both observations and the seasonal hindcast from six coupled atmosphere-ocean climate models from the 21 period of 1960–1980 (P1) to that of 1981–2001 (P2). It is shown that the one-month lead seasonal forecast skill of the six models’ multi-model ensemble is significantly increased from P1 to P2 for all four seasons. We identify four possible reasons accounting for the interdecadal change of the seasonal forecast skill. Firstly, the numerical model’s ability to simulate the mean state, the time variability and the spatial structures of the sea surface temperature and precipitation over the tropical Pacific is improved in P2 compared to P1. Secondly, an examination of the potential predictability of the atmosphere, estimated by the ratio of the total variance to the variance due to the internal dynamics of the model atmosphere, reveals that the atmospheric potential predictability is significantly increased after 1980s which is mainly due to an increased influence of El Niño-Southern Oscillation signal over the North Pacific and North American regions. Thirdly, the long-term climate trends in the atmosphere are found to contribute, to some extent, to the increased seasonal forecast skill especially over the Eurasian regions. Finally, the improved ocean observations in P2 may provide better initial conditions for the coupled models’ seasonal forecast.  相似文献   

11.
This study explores the decadal potential predictability of the Atlantic Meridional Overturning Circulation (AMOC) as represented in the IPSL-CM5A-LR model, along with the predictability of associated oceanic and atmospheric fields. Using a 1000-year control run, we analyze the prognostic potential predictability (PPP) of the AMOC through ensembles of simulations with perturbed initial conditions. Based on a measure of the ensemble spread, the modelled AMOC has an average predictive skill of 8 years, with some degree of dependence on the AMOC initial state. Diagnostic potential predictability of surface temperature and precipitation is also identified in the control run and compared to the PPP. Both approaches clearly bring out the same regions exhibiting the highest predictive skill. Generally, surface temperature has the highest skill up to 2 decades in the far North Atlantic ocean. There are also weak signals over a few oceanic areas in the tropics and subtropics. Predictability over land is restricted to the coastal areas bordering oceanic predictable regions. Potential predictability at interannual and longer timescales is largely absent for precipitation in spite of weak signals identified mainly in the Nordic Seas. Regions of weak signals show some dependence on AMOC initial state. All the identified regions are closely linked to decadal AMOC fluctuations suggesting that the potential predictability of climate arises from the mechanisms controlling these fluctuations. Evidence for dependence on AMOC initial state also suggests that studying skills from case studies may prove more useful to understand predictability mechanisms than computing average skill from numerous start dates.  相似文献   

12.
Summary This study investigates the capabilities of two regional models (the ICTP RegCM3 and the climate version of the CPTEC Eta model – EtaClim) in simulating the summer quasi-stationary circulations over South America during two extreme cases: the 1997–1998 El Ni?o and 1998–1999 La Ni?a. The results showed that both the models are successful in simulating the interannual variability of summer quasi-stationary circulation over South America. Both the models simulated the intensification of subtropical jet stream during the El Ni?o event, which favoured the blocking of transient systems and increased the precipitation over south Brazil. The models simulated the increase (decrease) of precipitation over north (west) Amazonia during the La Ni?a (El Ni?o) event. The upper level circulation is in agreement with the simulated distribution of precipitation. In general, the results showed that both the models are capable of capturing the main changes of the summer climate over South America during these two extreme cases and consequently they have potential to predict climate anomalies.  相似文献   

13.
C. Junquas  C. Vera  L. Li  H. Le Treut 《Climate Dynamics》2012,38(9-10):1867-1883
December–January–February (DJF) rainfall variability in southeastern South America (SESA) is studied in 18 coupled general circulation models from the WCRP/CMIP3 dataset, for present climate and the SRES-A1B climate change scenario. The analysis is made in terms of properties of the first leading pattern of rainfall variability in the region, characterized by a dipole-like structure with centers of action in the SESA and South Atlantic Convergence Zone (SACZ) regions. The study was performed to address two issues: how rainfall variability in SESA would change in a future climate and how much of that change explains the projected increasing trends in the summer mean rainfall in SESA identified in previous works. Positive (negative) dipole events were identified as those DJF seasons with above (below) normal rainfall in SESA and below (above) normal rainfall in the SACZ region. Results obtained from the multi-model ensemble confirm that future rainfall variability in SESA has a strong projection on the changes of seasonal dipole pattern activity, associated with an increase of the frequency of the positive phase. In addition, the frequency increase of positive dipole phase in the twenty first century seems to be associated with an increase of both frequency and intensity of positive SST anomalies in the equatorial Pacific, and with a Rossby wave train-like anomaly pattern linking that ocean basin to South America, which regionally induces favorable conditions for moisture transport convergence and rainfall increase in SESA.  相似文献   

14.
The SST-precipitation relationship in the intraseasonal variability (ISV) over the Asian monsoon region is examined using recent high quality satellite data and simulations from a state of the art coupled model, the climate forecast system version 2 (CFSv2). CFSv2 demonstrates high skill in reproducing the spatial distribution of the observed climatological mean summer monsoon precipitation along with its interannual variability, a task which has been a conundrum for many recent climate coupled models. The model also exhibits reasonable skill in simulating coherent northward propagating monsoon intraseasonal anomalies including SST and precipitation, which are generally consistent with observed ISV characteristics. Results from the observations and the model establish the existence of spatial variability in the atmospheric convective response to SST anomalies, over the Asian monsoon domain on intraseasonal timescales. The response is fast over the Arabian Sea, where precipitation lags SST by ~5 days; whereas it is slow over the Bay of Bengal and South China Sea, with a lag of ~12 days. The intraseasonal SST anomalies result in a similar atmospheric response across the basins, which consists of a destabilization of the bottom of the atmospheric column, as observed from the equivalent potential temperature anomalies near the surface. However, the presence of a relatively strong surface convergence over the Arabian Sea, due to the presence of a strong zonal gradient in SST, which accelerates the upward motion of the moist air, results in a relatively faster response in terms of the local precipitation anomalies over the Arabian Sea than over the Bay of Bengal and South China Sea. With respect to the observations, the ocean–atmosphere coupling is well simulated in the model, though with an overestimation of the intraseasonal SST anomalies, leading to an exaggerated SST-precipitation relationship. A detailed examination points to a systematic bias in the thickness of the mixed layer of the ocean model, which needs to be rectified. A too shallow (deep) mixed layer enhances (suppress) the amplitude of the intraseasonal SST anomalies, thereby amplifying (lessening) the ISV and the active-break phases of the monsoon in the model.  相似文献   

15.
Using the Flexible Global Ocean--Atmosphere--Land System model (FGOALS) version g1.11, a group of seasonal hindcasting experiments were carried out. In order to investigate the potential predictability of sea surface temperature (SST), singular value decomposition (SVD) analyses were applied to extract dominant coupled modes between observed and predicated SST from the hindcasting experiments in this study. The fields discussed are sea surface temperature anomalies over the tropical Pacific basin (20oS--20oN, 120oE--80oW), respectively starting in four seasons from 1982 to 2005. On the basis of SVD analysis, the simulated pattern was replaced with the corresponding observed pattern to reconstruct SST anomaly fields to improve the ability of the simulation. The predictive skill, anomaly correlation coefficients (ACC), after systematic error correction using the first five modes was regarded as potential predictability. Results showed that: 1) the statistical postprocessing approach was effective for systematic error correction; 2) model error sources mainly arose from mode 2 extracted from the SVD analysis---that is, during the transition phase of ENSO, the model encountered the spring predictability barrier; and 3) potential predictability (upper limits of predictability) could be high over most of the tropical Pacific basin, including the tropical western Pacific and an extra 10-degrees region of the mid and eastern Pacific.  相似文献   

16.
Heavy rains frequently occur over the Serra do Mar, in the southeast coastal mountain region in Brazil, particularly during the summer season. These rains can cause landslides and loss of life. The objective of this work is to produce a synoptic climatology of heavy rainfall episodes of the period from November to April using 10-year reanalyses data. The identification of the synoptic pattern of these events should provide guidance to forecasters. The landslide events were mostly related to cold front passages and the South Atlantic Convergence Zone (SACZ). These systems differed from the climatology by exhibiting more intense characteristics for the frontal situations and a blocking circulation configuration for the SACZ situations. In both cases, the composite fields showed that the 250-hPa mass divergence was strikingly more intense than the climatology and had a preferred location in the 24 h prior to landslide events. Anomalies of this ten-year event climatology showed above-normal moisture anomalies, which are more evident in the SACZ than in the frontal cases.  相似文献   

17.
基于NCEP/NCAR再分析资料,利用信噪比(SNR)方法研究了南海夏季风强度潜在可预报性的年代际变化,结果表明:南海夏季风强度潜在可预报性在20世纪80年代出现了由潜在可预报性偏低位相向偏高位相的年代际转折。进一步的研究发现,南海夏季风强度潜在可预报性和东印度洋—西太平洋(EIOWP)的海表温度(SST)存在明显的正相关。当EIOWP区域SST年际变率较大时,对南海夏季风影响较强,使得南海夏季风的外部信号增强,从而提高了潜在可预报性;当EIOWP区域SST年际变率较小时,对南海夏季风影响较弱,南海夏季风的外部信号进而减弱,潜在可预报性降低。  相似文献   

18.
The climate and variability of seasonal ensemble integrations, made with a recent version of ECMWF model (used for ERA-40 production) at relatively high horizontal resolution (TL159), have been studied for the 10-year period, 1980–1989. The model systematic error over the Atlantic-European region has been substantially improved when compared with the earlier model versions (e.g. from the PROVOST and AMIP-2 projects). However, it has worsened over the Pacific-North American region. This systematic error reduces the amplitude of planetary waves and has a negative impact on intraseasonal variability and predictability of the PNA mode. The signal-to-noise analysis yields results similar to earlier model versions: only during relatively strong ENSO events do some parts of the extratropics exhibit potential predictability. For precipitation, there is more disagreement between observed and model climatologies over sea than over land, but interannual variations over many parts of the tropical ocean are reasonably well represented. The south Asian monsoon in the model is severely weakened when compared to observations; this is seen in both poor climatology and interannual variability. Overall, comparing the ERA-40 model with earlier versions, there seems to be a balance between model improvements and deteriorations due to systematic errors. For the seasonal time-scale predictability, it is not clear that this model cycle constitutes an advantage over the earlier versions. Therefore, since it is not always possible to achieve distinct improvements in model climate and variability, a careful and detailed strategy ought to be considered when introducing a new model version for operational seasonal forecasting.  相似文献   

19.
青藏高原冬春积雪和地表热源影响亚洲夏季风的研究进展   总被引:2,自引:0,他引:2  
青藏高原冬春积雪和地表热源的气候效应是青藏高原气候动力学的两个重要内容。大量资料分析和数值试验研究均表明这两个因子对亚洲季风有一定的预测意义,本文对此做了比较系统的回顾和总结,并进一步比较了青藏高原积雪和地表热源影响东亚和南亚夏季降水的异同。结果表明,东亚夏季降水在年际和年代际尺度上均存在"三极型"和"南北反相"型的空间分布特征,高原春季地表热源在年代际和年际尺度上主要影响东亚夏季降水"三极型"模态;在年代际尺度上它是中国东部出现"南涝北旱"格局的重要原因,而高原冬季积雪的作用相反。另一方面,高原冬季积雪在年际和年代际尺度上对印度夏季风降水的预测效果均要优于高原地表热源。无论是空间分布还是时间演变特征,高原冬季积雪与春季地表热源整体上均无统计意义上的显著联系。不断完善高原地面观测网和改进模式在高原地区的模拟性能,将是进一步深入理解高原积雪和地表热源影响亚洲季风物理过程和机制的关键所在。  相似文献   

20.
Winter rainfall over South China shows strong interannual variability,which accounts for about half of the total winter rainfall over South China.This study investigated the predictability of winter (December-January-February; DJF) rainfall over South China using the retrospective forecasts of five state-of-the-art coupled models included in the ENSEMBLES project for the period 1961-2006.It was found that the ENSEMBLES models predicted the interannual variation of rainfall over South China well,with the correlation coefficient between the observed/station-averaged rainfall and predicted/areaaveraged rainfall being 0.46.In particular,above-normal South China rainfall was better predicted,and the correlation coefficient between the predicted and observed anomalies was 0.64 for these wetter winters.In addition,the models captured well the main features of SST and atmospheric circulation anomalies related to South China rainfall variation in the observation.It was further found that South China rainfall,when predicted according to predicted DJF Nifio3.4 index and the ENSO-South China rainfall relationship,shows a prediction skill almost as high as that directly predicted,indicating that ENSO is the source for the predictability of South China rainfall.  相似文献   

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

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