首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 890 毫秒
1.
Interactions between the tropical and subtropical northern Pacific at decadal time scales are examined using uncoupled oceanic and atmospheric simulations. An atmospheric model is forced with observed Pacific sea surface temperatures (SST) decadal anomalies, computed as the difference between the 2000–2009 and the 1990–1999 period. The resulting pattern has negative SST anomalies at the equator, with a global pattern reminiscent of the Pacific decadal oscillation. The tropical SST anomalies are responsible for driving a weakening of the Hadley cell and atmospheric meridional heat transport. The atmosphere is then shown to produce a significant response in the subtropics, with wind-stress-curl anomalies having the opposite sign from the climatological mean, consistent with a weakening of the oceanic subtropical gyre (STG). A global ocean model is then forced with the decadal anomalies from the atmospheric model. In the North Pacific, the shallow subtropical cell (STC) spins down and the meridional heat transport is reduced, resulting in positive tropical SST anomalies. The final tropical response is reached after the first 10 years of the experiment, consistent with the Rossby-wave adjustment time for both the STG and the STC. The STC provides the connection between subtropical wind stress anomalies and tropical SSTs. In fact, targeted simulations show the importance of off-equatorial wind stress anomalies in driving the oceanic response, whereas anomalous tropical winds have no role in the SST signal reversal. We further explore the connection between STG, STC and tropical SST with the help of an idealized model. We argue that, in our models, tropical SST decadal variability stems from the forcing of the Pacific subtropical gyre through the atmospheric response to ENSO. The resulting Ekman pumping anomaly alters the STC and oceanic heat transport, providing a negative feedback on the SST. We thus suggest that extratropical atmospheric responses to tropical forcing have feedbacks onto the ocean dynamics that lead to a time-delayed response of the tropical oceans, giving rise to a possible mechanism for multidecadal ocean-atmosphere coupled variability.  相似文献   

2.
Summary Hindcasts for the Indian summer monsoons (ISMs) of 2002 and 2003 have been produced from an ensemble of numerical simulations performed with a global model by changing SST. Two sets of ensemble simulations have been produced without vegetation: (i) by prescribing the weekly observed SST from ECMWF (European Centre for Medium Range Weather Forecasting) analyses, and (ii) by adding weekly SST anomalies (SSTA) of April to the climatological SST during the simulation period from May to August. For each ensemble, 10 simulations have been realized with different initial conditions that are prepared from ECMWF data with five each from April and May analyses of both the years. The predicted June–July monsoon rainfall over the Indian region shows good agreement with the GPCP (observed) pentad rainfall distribution when 5 member ensemble is taken from May initial conditions. The All-India June–July simulated rainfall time series matches favourably with the observed time series in both the years for the five member ensemble from May initial condition but drifts away from observation with April initial conditions. This underscores the role of initial conditions in the seasonal forecasting. But the model has failed to capture the strong intra-seasonal oscillation in July 2002. Heating over equatorial Indian Ocean for June 2002 in a particular experiment using 29th May 12 GMT as initial conditions shows some intra-seasonal oscillation in July 2002 rainfall, as in observation. Further evaluation of the seasonal simulations from this model is done by calculating the empirical orthogonal functions (EOFs) of the GPCP rainfall over India. The first four EOFs explain more than 80% of the total variance of the observed rainfall. The time series of expansion coefficients (principal components), obtained by projecting on the observed EOFs, provide a better framework for inter-comparing model simulations and their evaluation with observed data. The main finding of this study is that the All-India rainfall from various experiments with prescribed SST is better predicted on seasonal scale as compares to prescribed SST anomalies. This is indicative of a possible useful seasonal forecasts from a GCM at least for the case when monsoon is going to be good. The model responses do not differ much for 2002 and 2003 since the evolution of SST during these years was very similar, hence July rainfall seems to be largely modulated by the other feedbacks on the overall circulation.  相似文献   

3.
Sea surface temperature (SST) variations include negative feedbacks from the atmosphere, whereas SST anomalies are specified in stand-alone atmospheric general circulation simulations. Is the SST forced response the same as the coupled response? In this study, the importance of air–sea coupling in the Indian and Pacific Oceans for tropical atmospheric variability is investigated through numerical experiments with a coupled atmosphere-ocean general circulation model. The local and remote impacts of the Indian and Pacific Ocean coupling are obtained by comparing a coupled simulation with an experiment in which the SST forcing from the coupled simulation is specified in either the Indian or the Pacific Ocean. It is found that the Indian Ocean coupling is critical for atmospheric variability over the Pacific Ocean. Without the Indian Ocean coupling, the rainfall and SST variations are completely different throughout most of the Pacific Ocean basin. Without the Pacific Ocean coupling, part of the rainfall and SST variations in the Indian Ocean are reproduced in the forced run. In regions of large mean rainfall where the atmospheric negative feedback is strong, such as the North Indian Ocean and the western North Pacific in boreal summer, the atmospheric variability is significantly enhanced when air–sea coupling is replaced by specified SST forcing. This enhancement is due to the lack of the negative feedback in the forced SST simulation. In these regions, erroneous atmospheric anomalies could be induced by specified SST anomalies derived from the coupled model. The ENSO variability is reduced by about 20% when the Indian Ocean air–sea coupling is replaced by specified SST forcing. This change is attributed to the interfering roles of the Indian Ocean SST and Indian monsoon in western and central equatorial Pacific surface wind variations.  相似文献   

4.
Boreal winter North Atlantic climate change since 1950 is well described by a trend in the leading spatial structure of variability, known as the North Atlantic Oscillation (NAO). Through diagnoses of ensembles of atmospheric general circulation model (AGCM) experiments, we demonstrate that this climate change is a response to the temporal history of sea surface temperatures (SSTs). Specifically, 58 of 67 multi-model ensemble members (87%), forced with observed global SSTs since 1950, simulate a positive trend in a winter index of the NAO, and the spatial pattern of the multi-model ensemble mean trend agrees with that observed. An ensemble of AGCM simulations with only tropical SST forcing further suggests that variations in these SSTs are of primary importance. The probability distribution function (PDF) of 50-year NAO index trends from the forced simulations are, moreover, appreciably different from the PDF of a control simulation with no interannual SST variability, although chaotic atmospheric variations are shown to yield substantial 50-year trends. Our results thus advance the view that the observed linear trend in the winter NAO index is a combination of a strong tropically forced signal and an appreciable noise component of the same phase. The changes in tropical rainfall of greatest relevance include increased rainfall over the equatorial Indian Ocean, a change that has likely occurred in nature and is physically consistent with the observed, significant warming trend of the underlying sea surface.  相似文献   

5.
Three ensembles of AMIP-type simulations using the Arpege-climat coupled land–atmosphere model have been designed to assess the relative influence of SST and soil moisture (SM) on climate variability and predictability. The study takes advantage of the GSWP2 land surface reanalysis covering the 1986–1995 period. The GSWP2 forcings have been used to derive a global SM climatology that is fully consistent with the model used in this study. One ensemble of ten simulations has been forced by climatological SST and the simulated SM is relaxed toward the GSWP2 reanalysis. Another ensemble has been forced by observed SST and SM is evolving freely. The last ensemble combines the observed SST forcing and the relaxation toward GSWP2. Two complementary aspects of the predictability have been explored, the potential predictability (analysis of variance) and the effective predictability (skill score). An analysis of variance has revealed the effects of the SST and SM boundary forcings on the variability and potential predictability of near-surface temperature, precipitation and surface evaporation. While in the tropics SST anomalies clearly maintain a potentially predictable variability throughout the annual cycle, in the mid-latitudes the SST forced variability is only dominant in winter and SM plays a leading role in summer. In a similar fashion, the annual cycle of the hindcast skill (evaluated as the anomalous correlation coefficient of the three ensemble means with respect to the “observations”) indicates that the SST forcing is the dominant contributor over the tropical continents and in the winter mid-latitudes but that SM is supporting a significant part of the skill in the summer mid-latitudes. Focusing on boreal summer, we have then investigated different aspects of the SM and SST contribution to climate variations in terms of spatial distribution and time-evolution. Our experiments suggest that SM is potentially an additional source of climate predictability. A realistic initialization of SM and a proper representation of the land–atmosphere feedbacks seem necessary to improve state-of-the-art dynamical seasonal predictions, but will be actually efficient only in the areas where SM anomalies are themselves predictable at the monthly to seasonal timescale (since remote effects of SM are probably much more limited than SST teleconnections).  相似文献   

6.
The main goal of this study is to determine the oceanic regions corresponding to variability in African rainfall and seasonal differences in the atmospheric teleconnections. Canonical correlation analysis (CCA) has been applied in order to extract the dominant patterns of linear covariability. An ensemble of six simulations with the global atmospheric general circulation model ECHAM4, forced with observed sea surface temperatures (SSTs) and sea ice boundary variability, is used in order to focus on the SST-related part of African rainfall variability. Our main finding is that the boreal summer rainfall (June–September mean) over Africa is more affected by SST changes than in boreal winter (December–March mean). In winter, there is a highly significant link between tropical African rainfall and Indian Ocean and eastern tropical Pacific SST anomalies, which is closely related to El Niño-Southern Oscillation (ENSO). However, long-term changes are found to be associated with SST changes in the Indian and tropical Atlantic Oceans, thus, showing that the tropical Atlantic plays a critical role in determining the position of the intertropical convergence zone (ITCZ). Since ENSO is less in summer, the tropical Pacific and the Indian Oceans are less important for African rainfall. The African summer monsoon is strongly influenced by SST variations in the Gulf of Guinea, with a response of opposite sign over the Sahelian zone and the Guinean coast region. SST changes in the subtropical and extratropical oceans mostly take place on decadal time scales and are responsible for low-frequency rainfall fluctuations over West Africa. The modelled teleconnections are highly consistent with the observations. The agreement for most of the teleconnection patterns is remarkable and suggests that the modelled rainfall anomalies serve as suitable predictors for the observed changes.  相似文献   

7.
In this study we investigate the impact of large-scale oceanic forcing and local vegetation feedback on the variability of the Sahel rainfall using a global biosphere-atmosphere model, the coupled GENESIS-IBIS model, running at two different resolutions. The observed global sea surface temperature in the twentieth century is used as the primary model forcing. Using this coupled global model, we experiment on treating vegetation as a static boundary condition and as a dynamic component of the Earth climate system. When vegetation is dynamic, the R30-resolution model realistically reproduces the multi-decadal scale fluctuation of rainfall in the Sahel region; keeping vegetation static in the same model results in a rainfall regime characterized by fluctuations at much shorter time scales, indicating that vegetation dynamics act as a mechanism for persistence of the regional climate. Even when vegetation dynamics is included, the R15 model fails to capture the main characteristics of the long-term rainfall variability due to the exaggerated atmospheric internal variability in the coarse resolution model. Regardless how vegetation is treated and what model resolution is used, conditions in the last three decades of the twentieth century are always drier than normal in the Sahel, suggesting that global oceanic forcing during that period favors the occurrence of a drought. Vegetation dynamics is found to enhance the severity of this drought. However, with both the observed global SST forcing and feedback from dynamic vegetation in the model, the simulated drought is still not as persistent as that observed. This indicates that anthropogenic land cover changes, a mechanism missing in the model, may have contributed to the occurrence of the twentieth century drought in the Sahel.  相似文献   

8.
9.
Decadal climate predictability is examined in hindcast experiments by a multi-model ensemble using three versions of the coupled atmosphere-ocean model MIROC. In these hindcast experiments, initial conditions are obtained from an anomaly assimilation procedure using the observed oceanic temperature and salinity with prescribed natural and anthropogenic forcings on the basis of the historical data and future emission scenarios in the Intergovernmental Panel of Climate Change. Results of the multi-model ensemble in our hindcast experiments show that predictability of surface air temperature (SAT) anomalies on decadal timescales mostly originates from externally forced variability. Although the predictable component of internally generated variability has considerably smaller SAT variance than that of externally forced variability, ocean subsurface temperature variability has predictive skills over almost a decade, particularly in the North Pacific and the North Atlantic where dominant signals associated with Pacific decadal oscillation (PDO) and the Atlantic multidecadal oscillation (AMO) are observed. Initialization enhances the predictive skills of AMO and PDO indices and slightly improves those of global mean temperature anomalies. Improvement of these predictive skills in the multi-model ensemble is higher than that in a single-model ensemble.  相似文献   

10.
The influence of prescribed changes in vegetation on the climate of the North American monsoon region is examined using the National Center for Atmospheric Research Community Climate System Model Version 3.5 (NCAR CCSM3.5). Initial value ensemble experiments are performed in which the vegetation cover fraction over the North American monsoon region is reduced by 0.2 and the intra-annual climatic response is assessed probabilistically in each one-year ensemble experiment. Changes in the surface radiation budget include decreases in sensible and latent heat fluxes and increases in upward longwave and downward shortwave radiation fluxes, with small net changes in surface albedo. The climatic responses to reduced vegetation cover fraction include year-round increases in ground and surface air temperature, a dampened hydrologic cycle with decreased springtime evaporation, springtime and autumnal precipitation, and autumnal cloud cover, and enhanced atmospheric subsidence in late autumn. Decreased vegetation shifts the monsoon season over the Southwest United States earlier in the year. Within the North American monsoon region, the most robust vegetation feedbacks to climate are found over woody landscapes.  相似文献   

11.
Influence of SST biases on future climate change projections   总被引:1,自引:0,他引:1  
We use a quantile-based bias correction technique and a multi-member ensemble of the atmospheric component of NCAR CCSM3 (CAM3) simulations to investigate the influence of sea surface temperature (SST) biases on future climate change projections. The simulations, which cover 1977?C1999 in the historical period and 2077?C2099 in the future (A1B) period, use the CCSM3-generated SSTs as prescribed boundary conditions. Bias correction is applied to the monthly time-series of SSTs so that the simulated changes in SST mean and variability are preserved. Our comparison of CAM3 simulations with and without SST correction shows that the SST biases affect the precipitation distribution in CAM3 over many regions by introducing errors in atmospheric moisture content and upper-level (lower-level) divergence (convergence). Also, bias correction leads to significantly different precipitation and surface temperature changes over many oceanic and terrestrial regions (predominantly in the tropics) in response to the future anthropogenic increases in greenhouse forcing. The differences in the precipitation response from SST bias correction occur both in the mean and the percent change, and are independent of the ocean?Catmosphere coupling. Many of these differences are comparable to or larger than the spread of future precipitation changes across the CMIP3 ensemble. Such biases can affect the simulated terrestrial feedbacks and thermohaline circulations in coupled climate model integrations through changes in the hydrological cycle and ocean salinity. Moreover, biases in CCSM3-generated SSTs are generally similar to the biases in CMIP3 ensemble mean SSTs, suggesting that other GCMs may display a similar sensitivity of projected climate change to SST errors. These results help to quantify the influence of climate model biases on the simulated climate change, and therefore should inform the effort to further develop approaches for reliable climate change projection.  相似文献   

12.
We quantify the feedbacks from the physical climate system on the radiative forcing for idealized climate simulations using four different methods. The results differ between the methods and differences are largest for the cloud feedback. The spatial and temporal variability of each feedback is used to estimate the averaging scale necessary to satisfy the feedback concept of one constant global mean value. We find that the year-to-year variability, combined with the methodological differences, in estimates of the feedback strength from a single model is comparable to the model-to-model spread in feedback strength of the CMIP3 ensemble. The strongest spatial and temporal variability is in the short-wave component of the cloud feedback. In our simulations, where many sources of natural variability are neglected, long-term averages are necessary to get reliable feedback estimates. Considering the large natural variability and relatively small forcing present in the real world, as compared to the forcing imposed by doubling CO2 concentrations in the simulations, implies that using observations to constrain feedbacks is a challenging task and requires reliable long-term measurements.  相似文献   

13.
The impact of internal atmospheric variability on North Pacific sea surface temperature (SST) variability is examined based on three coupled general circulation model simulations. The three simulations differ only in the level of atmospheric noise occuring over the ocean at the air-sea interface. The amplitude of atmospheric noise is controlled by use of the interactive ensemble technique. This technique simultaneously couples multiple realizations of a single atmospheric model to a single realization of an ocean model. The atmospheric component models all experience the same SST, but the ocean component is forced by the ensemble averaged fluxes thereby reducing the impact of internal atmospheric dynamics at the air-sea interface. The ensemble averaging is only applied at the air-sea interface so that the internal atmospheric dynamics (i.e., transients) of each atmospheric ensemble member is unaffected. This interactive ensemble technique significantly reduces the SST variance throughout the North Pacific. The reduction in SST variance is proportional to the number of ensemble members indicating that most of the variability can simply be explained as the response to atmospheric stochastic forcing. In addition, the impact of the internal atmospheric dynamics at the air-sea interface masks out much of the tropical-midlatitude SST teleconnections on interannual time scales. Once this interference is reduced (i.e., applying the interactive ensemble technique), tropical-midlatitude SST teleconnections are easily detected.  相似文献   

14.
This is a study of the annual and interannual variability of regional rainfall produced by the Center for Weather Forecasts and Climate Studies/Center for Ocean, Land and Atmospheric Studies (CPTEC/COLA) atmospheric global climate model. An evaluation is made of a 9-member ensemble run of the model forced by observed global sea surface temperature (SST) anomalies for the 10-year period 1982–1991. The Brier skill score and, Relative Operating Characteristics (ROC) are used to assess the predictability of rainfall and to validate rainfall simulations, in several regions world wide. In general, the annual cycle of precipitation is well simulated by the model for several continental and oceanic regions in the tropics and mid latitudes. Interannual variability of rainfall during the peak rainy season is realistically simulated in Northeast Brazil, Amazonia, central Chile, and southern Argentina–Uruguay, Eastern Africa, and tropical Pacific regions, where the model shows good skill. Some regions, such as northwest Peru–Ecuador, and southern Brazil exhibit a realistic simulation of rainfall anomalies associated with extreme El Niño warming conditions, while in years with neutral or La Niña conditions, the agreement between observed and simulated rainfall anomalies is not always present. In the monsoon regions of the world and in southern Africa, even though the model reproduces the annual cycle of rainfall, the skill of the model is low for the simulation of the interannual variability. This is indicative of mechanisms other than the external SST forcing, such as the effect of land–surface moisture and snow feedbacks or the representation of sub-grid scale processes, indicating the important role of factors other than external boundary forcing. The model captures the well-known signatures of rainfall anomalies of El Niño in 1982–83 and 1986–87, indicating its sensitivity to strong external forcing. In normal years, internal climate variability can affect the predictability of climate in some regions, especially in monsoon areas of the world.  相似文献   

15.
A nine-member ensemble of simulations with a state-of-the-art atmospheric model forced only by the observed record of sea surface temperature (SST) over 1930–2000 is shown to capture the dominant patterns of variability of boreal summer African rainfall. One pattern represents variability along the Gulf of Guinea, between the equator and 10°N. It connects rainfall over Africa to the Atlantic marine Intertropical Convergence Zone, is controlled by local, i.e., eastern equatorial Atlantic, SSTs, and is interannual in time scale. The other represents variability in the semi-arid Sahel, between 10°N and 20°N. It is a continental pattern, capturing the essence of the African summer monsoon, while at the same time displaying high sensitivity to SSTs in the global tropics. A land–atmosphere feedback associated with this pattern translates precipitation anomalies into coherent surface temperature and evaporation anomalies, as highlighted by a simulation where soil moisture is held fixed to climatology. As a consequence of such feedback, it is shown that the recent positive trend in surface temperature is consistent with the ocean-forced negative trend in precipitation, without the need to invoke the direct effect of the observed increase in anthropogenic greenhouse gases. We advance plausible mechanisms by which the balance between land–ocean temperature contrast and moisture availability that defines the monsoon could have been altered in recent decades, resulting in persistent drought. This discussion also serves to illustrate ways in which the monsoon may be perturbed, or may already have been perturbed, by anthropogenic climate change.  相似文献   

16.
The present study examined the major features of the interdecadal variation of the summer rainfall over eastern China (IVRC) and the possible association with sea surface temperature (SST). We noted that the first leading mode of IVRC (accounting for nearly half of the total variance and with maximum loading for the summer rainfall anomalies over South China) may be not forced by SST. On the other hand, the second and third leading modes [accounting for 17.1 and 13.6 % of the total variance and mainly associated with the summer rainfall anomalies over the Yangtze River valley (YRV) and North China, respectively] in some extent are forced by SST anomalies. These observational results are confirmed by atmospheric general circulation model (AGCM) simulations forced by observed SST. By eliminating the internal dynamical process driven rainfall though ensemble mean, the simulations further suggest an overall enhancement of the intensity of IVRC in the corresponding ensemble mean, especially in the YRV and North China regions, but not in South China. That implies the different role of SST in driving IVRC over different regions.  相似文献   

17.
A regional climate model is used to investigate the mechanism of interdecadal rainfall variability, specifically the drought of the 1970s and 1980s, in the Sahel region of Africa. The model is the National Center for Environmental Prediction’s (NCEPs) Regional Spectral Model (RSM97), with a horizontal resolution of approximately equivalent to a grid spacing of 50 km, nested within the ECHAM4.5 atmospheric general circulation model (AGCM), which in turn was forced by observed sea surface temperature (SST). Simulations for the July–September season of the individual years 1955 and 1986 produced wet conditions in 1955 and dry conditions in 1986 in the Sahel, as observed. Additional July–September simulations were run forced by SSTs averaged for each month over the periods 1950–1959 and the 1978–1987. These simulations yielded wet conditions in the 1950–1959 case and dry conditions in the 1978–1987 case, confirming the role of SST forcing in decadal variability in particular. To test the hypothesis that the SST influences Sahel rainfall via stabilization of the tropospheric sounding, simulations were performed in which the temperature field from the AGCM was artificially modified before it was used to force the regional model. We modified the original 1955 ECHAM4.5 temperature profiles by adding a horizontally uniform, vertically varying temperature increase, taken from the 1986–1955 tropical mean warming in either the AGCM or the NCEP/National Center for Atmospheric Research Reanalysis. When compared to the 1955 simulations without the added tropospheric warming, these simulations show a drying in the Sahel similar to that in the 1986–1955 difference and to the decadal difference between the 1980s and 1950s. This suggests that the tropospheric warming may have been, at least in part, the agent by which the SST increases led to the Sahel drought of the 1970s and 1980s.  相似文献   

18.
Observations indicate that recent tropical Pacific decadal climate variability tends to be associated with the extratropical North Pacific through a relay teleconnection of a fast coupled ocean-atmosphere bridge and a slow oceanic tunnel. A coupled ocean-atmosphere model, forced by the observed decadal wind in the extratropical North Pacific, explicitly demonstrates that extratropical decadal sea surface temperature (SST) anomalies may propagate to the tropics through a coupled wind-evaporative-SST (WES) feedback. The WES feedback cannot only lead to a nearly synchronous change of tropical SST, but also force a delayed adjustment of the meridional overturning circulation in the upper ocean to further sustain the tropical SST change. The study further suggests that the extratropical–tropical teleconnection provides a positive feedback to sustain the decadal changes in both the tropical and extratropical North Pacific.  相似文献   

19.
This study presents the spatial and temporal structures of the decadal variability of the Pacific from an extended control run of a coupled global climate model (GCM).The GCM used was version-g2.0 of the Flexible Global Ocean Atmosphere Land System (FGOALS-g2.0) developed at LASG/IAP.The GCM FGOALS-g2.0 re-produces similar spatial-temporal structures of sea surface temperature (SST) as observed in the Pacific decadal os-cillation (PDO) with a significant period of approximately 14 years.Correspondingly,the PDO signals were closely related to the decadal change both in the upper-ocean temperature anomalies and in the atmospheric circulation.The present results suggest that warm SST anomalies along the equator relax the trade winds,causing the SSTs to warm even more in the eastern equatorial Pacific,which is a positive feedback.Meanwhile,warm SST anomalies along the equator force characteristic off-equa-torial wind stress curl anomalies,inducing much more poleward transport of heat,which is a negative feedback.The upper-ocean meridional heat transport,which is asso-ciated with the PDO phase transition,links the equatorial to the off-equatorial Pacific Ocean,acting as a major mechanism responsible for the tropical Pacific decadal variations.Therefore,the positive and negative feedbacks working together eventually result in the decadal oscilla-tion in the Pacific.  相似文献   

20.
The sea surface temperature (SST) of the tropical western Pacific Ocean (TWPO) showed a pronounced warming in the late 1990s. Using numerical experiments of a regional climate model (RegCM), we analyzed the impact of this warming on rainfall over the Yangtze-Huaihe River valley of China during the Meiyu period (June-July). The model results revealed that the observed decadal changes in Meiyu rainfall since the late 1990s can be reproduced by a control experiment forced by the observed SST. Additionally, the sensitivity experiments suggested that the warming trend in the TWPO played a substantial role in the northward shift of the Meiyu belt in the late 1990s.  相似文献   

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

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