首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Climate changes during the next 100 years caused by anthropogenic emissions of greenhouse gases have been simulated for the Intergovernmental Panel on Climate Change Scenarios A (business as usual) and D (accelerated policies) using a coupled ocean-atmosphere general circulation model. In the global average, the near-surface temperature rises by 2.6 K in Scenario A and by 0.6 K in Scenario D. The global patterns of climate change for both IPCC scenarios and for a third step-function 2 x CO2 experiment were found to be very similar. The warming delay over the oceans is larger than found in simulations with atmospheric general circulation models coupled to mixed-layer models, leading to a more pronounced land-sea contrast and a weaker warming (and in some regions even an initial cooling) in the Southern Ocean. During the first forty years, the global warming and sea level rise due to the thermal expansion of the ocean are significantly slower than estimated previously from box-diffusion-upwelling models, but the major part of this delay can be attributed to the previous warming history prior to the start of present coupled ocean-atmosphere model integration (cold start).  相似文献   

2.
A new, biogeochemical model of ice age cycles is developed and applied which explains major features of climate variations in the late Quaternary —rapid ice age terminations, large glacial-interglacial amplitudes and 100-kyr cycles — in a way consistent with the paleorecord. Existing models which invoke non-linear, ice-sheet-earth-crust dynamics to explain ice age cycles are not consistent with simultaneous terminations in both hemispheres and other phase relationships implied by the paleorecord. The present model relates climate change to oscillations of oceanic primary (new) production controlled by the availability of inorganic nitrogen. Large oscillations follow shelf erosion events triggered by small sea-level drops. These drops are due to glacial buildup associated with a minimum in Northern Hemisphere insolation. Rapid global warming at terminations is initiated by open ocean denitrification events leading to new production crashes and rapid modification of atmospheric trace gas concentrations (CO2, DMS, N2O). Other feedbacks of the land-ice-atmosphere-ocean system control the rest of the climate cycle. 100-kyr cycles derive from orbital pacemaking of the strong, low-frequency model response. Results suggest that the climate regime transition near 800 kyr B.P. may be related to changes in the continental shelf slope, that existing chronologies based on orbital tuning may need to be revised and that temporary increases in atmospheric N2O concentrations at terminations, due to the denitrification events, may have caused significant greenhouse warming. A spike of elevated N2O concentration at terminations may be recorded in polar ice.  相似文献   

3.
Progress in understanding how terrestrial ice volume is linked to Earths orbital configuration has been impeded by the cost of simulating climate system processes relevant to glaciation over orbital time scales (103–105 years). A compromise is usually made to represent the climate system by models that are averaged over one or more spatial dimensions or by three-dimensional models that are limited to simulating particular snapshots in time. We take advantage of the short equilibration time (10 years) of a climate model consisting of a three-dimensional atmosphere coupled to a simple slab ocean to derive the equilibrium climate response to accelerated variations in Earths orbital configuration over the past 165,000 years. Prominent decreases in ice melt and increases in snowfall are simulated during three time intervals near 26, 73, and 117 thousand years ago (ka) when aphelion was in late spring and obliquity was low. There were also significant decreases in ice melt and increases in snowfall near 97 and 142 ka when eccentricity was relatively large, aphelion was in late spring, and obliquity was high or near its long term mean. These glaciation-friendly time intervals correspond to prominent and secondary phases of terrestrial ice growth seen within the marine 18O record. Both dynamical and thermal effects contribute to the increases in snowfall during these periods, through increases in storm activity and the fraction of precipitation falling as snow. The majority of the mid- to high latitude response to orbital forcing is organized by the properties of sea ice, through its influence on radiative feedbacks that nearly double the size of the orbital forcing as well as its influence on the seasonal evolution of the latitudinal temperature gradient.  相似文献   

4.
Summary ¶In order to better understand land-atmosphere interactions and increase the predictability of climate models, it is important to investigate the role of forest representation in climate modeling. Corresponding to the big-leaf model commonly employed in land surface schemes to represent the effects of a forest, a so called big-tree model, which uses multi-layer vegetation to represent the vertical canopy heterogeneity, was introduced and incorporated into the National Center for Atmospheric Research (NCAR) regional climate model RegCM2, to make the vegetation model more physically based. Using this augmented RegCM2 and station data for China during 1991 Meiyu season, we performed 10 experiments to investigate the effects of the application of the big-tree model on the summer monsoon climate.With the big-tree model incorporated into the regional climate model, some climate characteristics, e.g. the 3-month-mean surface temperature, circulation, and precipitation, are significantly and systematically changed over the model domain, and the change of the characteristics differs depending on the area. Due to the better representation of the shading effect in the big-tree model, the temperature of the lower layer atmosphere above the plant canopy is increased, which further influences the 850hPa temperature. In addition, there are significant decreases in the mean latent heat fluxes (within 20–30W/m2) in the three areas of the model domain.The application of the big-tree model influences not only the simulated climate of the forested area, but also that of the whole model domain, and its impact is greater on the lower atmosphere than on the upper atmosphere. The simulated rainfall and surface temperature deviate from the originally simulated result and are (or seem to be) closer to the observations, which implies that an appropriate representation of the big-tree model may improve the simulation of the summer monsoon climate.We also find that the simulated climate is sensitive to some big-tree parameter values and schemes, such as the shape, height, zero-plane displacement height and mixing-length scheme. The simulated local/grid differences may be very large although the simulated areal-average differences may be much lower. The area-average differences in the monthly-mean surface temperature and heat fluxes can amount to 0.5°C and 4W/m2, respectively, which correspond to maximum local/grid differences of 3.0°C and 40W/m2 respectively. It seems that the simulated climate is most sensitive to the parameter of the zero-plane displacement among the parameters studied.  相似文献   

5.
The primary goal of this investigation is to focus on a realistic scenario for simulating impacts on regional African climate of future deforestation in a greenhouse-warmed world. Combined effects of plausible land-cover change and greenhouse warming are assessed by time-slice simulations with an atmospheric general circulation model (AGCM) for the middle of the twenty first century. Three time-slice integrations have been performed with the ARPEGE-Climat AGCM incorporating a zooming technique to achieve a resolution of about 100 km over Africa. A control run for the current climate is forced by observed climatological sea surface temperatures (SSTs) and the observed vegetation distribution is specified from a new vegetation database, in order to improve the geographical distribution and properties of the vegetation cover. Future SST changes are derived from a transient coupled atmosphere–ocean simulation for scenario B2 of the International Panel on Climate Change (IPCC). Future vegetation changes are specified from a simulation of scenario B2 with the Integrated Model to Assess the Global Environment (IMAGE) developed at the National Institute of Public Health and the Environment in the Netherlands (RIVM). The results show that land surface processes can locally modulate greenhouse warming effects for African climate, with reductions of surface transpiration and small increases of surface temperature. Deforestation of tropical Africa has overall only a marginal effect on precipitation because of a compensatory increase in moisture convergence. Energy budget analyses show that increases in surface temperature are produced both by increases of greenhouse gases (GHG) concentration from the increase in downward atmospheric longwave radiation, and by African tropical deforestation from the resulting reduction in transpiration. This study indicates that realistic land-use changes, though of smaller amplitude than greenhouse gas forcing, may have a small regional effect in projections of future climate.  相似文献   

6.
Due to restrictions in the available computing resources and a lack of suitable observational data, transient climate change experiments with global coupled ocean-atmosphere models have been started from an initial state at equilibrium with the present day forcing. The historical development of greenhouse gas forcing from the onset of industrialization until the present has therefore been neglected. Studies with simplified models have shown that this cold start error leads to a serious underestimation of the anthropogenic global warming. In the present study, a 150-year integration has been carried out with a global coupled ocean-atmosphere model starting from the greenhouse gas concentration observed in 1935, i.e., at an early time of industrialization. The model was forced with observed greenhouse gas concentrations up to 1985, and with the equivalent C02 concentrations stipulated in Scenario A (Business as Usual) of the Intergovernmental Panel on Climate Change from 1985 to 2085. The early starting date alleviates some of the cold start problems. The global mean near surface temperature change in 2085 is about 0.3 K (ca. 10%) higher in the early industrialization experiment than in an integration with the same model and identical Scenario A greenhouse gas forcing, but with a start date in 1985. Comparisons between the experiments with early and late start dates show considerable differences in the amplitude of the regional climate change patterns, particularly for sea level. The early industrialization experiment can be used to obtain a first estimate of the detection time for a greenhouse-gas-induced near-surface temperature signal. Detection time estimates are obtained using globally and zonally averaged data from the experiment and a long control run, as well as principal component time series describing the evolution of the dominant signal and noise modes. The latter approach yields the earliest detection time (in the decade 1990–2000) for the time-evolving near-surface temperature signal. For global-mean temperatures or for temperatures averaged between 45°N and 45°S, the signal detection times are in the decades 2015–2025 and 2005–2015, respectively. The reduction of the cold start error in the early industrialization experiment makes it possible to separate the near-surface temperature signal from the noise about one decade earlier than in the experiment starting in 1985. We stress that these detection times are only valid in the context of the coupled model's internally-generated natural variability, which possibly underestimates low frequency fluctuations and does not incorporate the variance associated with changes in external forcing factors, such as anthropogenic sulfate aerosols, solar variability or volcanic dust.  相似文献   

7.
Grid transformation for incorporating the Arctic in a global ocean model   总被引：2，自引：0，他引：2  
A grid transformation is described which isolates the Arctic and North Atlantic, rotates the spherical grid to pass an equator up the Atlantic through the north pole and remaps the Coriolis parameter. Boundary condition information is exchanged along the equatorial Atlantic so that the Arctic-Atlantic model is dynamically coupled to a model of the rest of the global ocean (which remains on the geographic spherical grid). The transform produces a more regular grid over the Arctic and eliminates the need for filtering or special treatment at the pole. The transform has been implemented in the GFDL Modular Ocean Model. After testing with idealized geometry, a 300 y global integration is compared to an integration using the geographic spherical grid and Fourier filtering. Results are similar, with differences in the Arctic and western North Atlantic regions leading to smaller air-sea heat flux near the Gulf Stream separation latitude for the transform case. Use of the transform also leads to a reduction in computation time.  相似文献   

8.
We have studied the response of the Atlantic meridional overturning circulation to surface freshwater forcing using an ocean GCM coupled to an energy-moisture-balance atmosphere model. The overturning collapses rapidly when a slowly increasing forcing applied to the North Atlantic passes a positive threshold, and spins up equally quickly when the forcing falls below a negative threshold. This well-known behaviour is referred to as hysteresis because the thresholds in forcing are different for the transitions in opposite directions. However, we argue that the behaviour of the Atlantic salinity is more fundamental than the forcing. Hysteresis as a function of freshwater forcing occurs because the states with North Atlantic overturning on and off each tend to reinforce their associated salinity distributions and inhibit the transition to the other state. During the collapse, the Atlantic becomes less saline because of the import of 80 Sv year of freshwater by ocean transports across 30°S; during the spin-up this freshwater is exported again. We show that qualitatively similar hysteresis behaviour can be produced by perturbing the system without any net freshwater forcing. The salinity flip-flop is associated with the appearance and disappearance of a shallow reverse overturning circulation south of the Equator, which is present while the northern overturning is absent, and may provide the mechanism for the ocean freshwater influx during collapse.  相似文献   

9.
Summary Regional climate model (RegCM2) and sulfur transport model (NJUADMS) were combined to simulate the distribution of anthropogenic sulfate aerosol burden over China, where a look up table method was applied to illustrate sulfate formation from SO2-oxidation. Direct radiative forcing of sulfate aerosol was further estimated using the scheme suggested by Charlson et al (1991). Investigations show that the annual average total sulfate column over mainland China is 2.01mg/m2 with high value in East and Central areas (more than 7mg/m2). The annual average direct radiative forcing of China is about –0.85W/m2. The forcing can reach –7W/m2 in Central and East China during the winter season. Total sulfate column shows significant seasonal variations with winter maximum-summer minimum in the Southern part of China and spring maximum-autumn minimum in the northern part of China. Strong seasonal cycles of direct radiative forcing are also found due to the influence of total sulfate column, cloud, relative humidity and the reflectivity of underlying surfaceReceived May 16, 2001; accepted August 5, 2002 Published online: May 8, 2003  相似文献   

10.
The exact Eulerian velocity probability density function (pdf) of a turbulent field is generally unknown, and one normally has available only partial information in the form of low order moments. We compare two alternative Lagrangian Stochastic (LS) approaches formed from this partial information, (i) the moments approximation approach (Kaplan and Dinar, 1993); and (ii) the well-mixed model (Thomson, 1987) that corresponds to the maximum missing information pdf formed from the available information. We show that the moments approximation model does not in general satisfy the well-mixed constraint, and can give an inferior prediction of dispersion.  相似文献   

11.
Climate drift is a common and serious problem in most state-of-the-art coupled atmosphere-ocean-sea ice models. We consider the nature of climate drift in such a model, and in particular address the question of whether or not climate drift is inherent to the model, or whether the drift can be averted by a suitable choice of initial conditions or coupling procedure. The synchronous approach to coupling was adopted in which the ocean, atmosphere and sea ice models were spun-up independently to equilibrium using climatological forcing fields. The models were then coupled and integrated forward in time. Several experiments were performed which were designed to assess the impact of different coupling methodologies and changes in the initial conditions of the component models on the climate drift of the system. The results of our experiments indicate that climate drift is a problem inherent to the coupled model in that systematic errors in the components lead to incompatibilities in the surface fluxes required by the component models to maintain realistic climatologies. We conclude that climate drift can be averted only if the parameterizations of certain important physical processes are improved which should have the effect of reducing or eliminating these incompatibilities.  相似文献   

12.
Radiosondes releases during the NOPEX-WINTEX experiment carried out in late winter in Northern Finland were analysed for the determination of the height h of the atmospheric boundary layer. We investigate various possible scaling approaches, based on length scales using micrometeorological turbulence surface measurements and the background atmospheric stratification above h. Under stable conditions, the three previously observed turbulence regimes delineated by values of z/L (L is the Obukhov length) appears as a blueprint for understanding the departures found for the suitability of the Ekman scaling based on LE = u/f (u is the friction velocity and f the Coriolis parameter). The length scale LN = u/N (where N is the Brunt–Väisälä frequency) appears to be a useful scale under most stable conditions, especially in association with L. Under unstable conditions, shear production of turbulence is still significant, so that the three scales L, LN and LE are again relevant and the dimensionless ratios N = LN/L and LN/LE = N/f describe well the WINTEX data. Furthermore, in the classical scaling framework, the unstable domain may also be divided into three regimes as reflected by the dependence ofu/f on instability (z/L).  相似文献   

13.
Summary The first GCM climate change projections to include dynamic vegetation and an interactive carbon cycle produced a very significant amplification of global warming over the 21st century. Under the IS92a business as usual emissions scenario CO2 concentrations reached about 980ppmv by 2100, which is about 280ppmv higher than when these feedbacks were ignored. The major contribution to the increased CO2 arose from reductions in soil carbon because global warming is assumed to accelerate respiration. However, there was also a lesser contribution from an alarming loss of the Amazonian rainforest. This paper describes the phenomenon of Amazonian forest dieback under elevated CO2 in the Hadley Centre climate-carbon cycle model.  相似文献   

14.
The effect of changes in zonal and meridional atmospheric moisture transports on Atlantic overturning is investigated. Zonal transports are considered in terms of net moisture export from the Atlantic sector. Meridional transports are related to the vigour of the global hydrological cycle. The equilibrium thermohaline circulation (THC) simulated with an efficient climate model is strongly dependent on two key parameters that control these transports: an anomaly in the specified Atlantic–Pacific moisture flux (Fa) and atmospheric moisture diffusivity (Kq). In a large ensemble of spinup experiments, the values of Fa and Kq are varied by small increments across wide ranges, to identify sharp transitions of equilibrium THC strength in a 2-parameter space (between Conveyor On and Off states). Final states from this ensemble of simulations are then used as the initial states for further such ensembles. Large differences in THC strength between ensembles, for identical combinations of Fa and Kq, reveal the co-existence of two stable THC states (Conveyor On and Off)—i.e. a bistable regime. In further sensitivity experiments, the model is forced with small, temporary freshwater perturbations to the mid-latitude North Atlantic, to establish the minimum perturbation necessary for irreversible THC collapse in this bistable regime. A threshold is identified in terms of the forcing duration required. The model THC, in a Conveyor On state, irreversibly collapses to a Conveyor Off state under additional freshwater forcing of just 0.1 Sv applied for around 100 years. The irreversible collapse is primarily due to a positive feedback associated with suppressed convection and reduced surface heat loss in the sinking region. Increased atmosphere-to-ocean freshwater flux, under a collapsed Conveyor, plays a secondary role.  相似文献   

15.
Summary The evolving modes of the sea-surface temperature (SST) in the Tropical Atlantic on the short interannual (IA) timescale were obtained by performing the extended empirical orthogonal function (EEOF) analyses on this variable separately for the 106-year (1871–1976) and 20-year (1881–1900; 1901–1920; 1921–1940; 1941–1960) periods. The equatorial and inter-hemispheric patterns manifest in the first EEOF mode of each analysis as part of the short IA evolution of the SST anomalies in the Tropical Atlantic. Another outstanding feature of the first EEOF mode of each analysis concerns the propagations of the SST anomalies in the meridional direction within the 20°N–20°S band and in the zonal direction in the sector 40°W–20°W. For all analyses, the SST anomalies propagate northward from the equator to 15°N and southward from 20°N to 15°N, with the same sign anomalies merging approximately at 15°N. On the other hand, the SST anomalies propagate westward in the sector 40°W–20°W with a propagation rate close to that of the phase speed of the fastest baroclinic Rossby wave in the ocean. So, the observed propagations of the SST anomalies in the 20°N–20°S band might result from the combined effect of the surface oceanic currents in this band and the baroclinic Rossby waves in the ocean.  相似文献   

16.
Parameterization of evaporation from a non-plant-covered surface is very important in the hierarchy strategy of modelling land surface processes. One of the representations frequently used in its computation is the resistance formulation. The performance of the evaporation schemes using the , , and their combination resistance approaches to parameterize evaporation from bare soil surfaces is discussed. For that purpose, the nine schemes, based on a different dependence of and on volumetric soil moisture content and its saturated value, are used.The tests of performances of the considered schemes are based on time integrations by the land surface module (BARESOIL) using observed data. The 23 data sets at a bare surface experimental site in Rimski anevi, Yugoslavia on chernozem soil, were used for the resistance algorithm evaluation. The quality of the schemes was compared with the observed values of the latent heat flux using several statistical parameters.  相似文献   

17.
The impact of climate change on the river rhine: A scenario study   总被引：3，自引：0，他引：3  
This paper concerns the impact of human-induced global climate change on the River Rhine discharge. For this purpose a model for climate assessment, named ESCAPE, is coupled to a water balance model, named RHINEFLOW. From climate scenarios, changes in regional annual water availability and seasonal discharge in the River Rhine Basin are estimated. The climate scenarios are based on greenhouse gases emissions scenarios. An assessment is made for best guess seasonal discharge changes and for changes in frequencies of low and high discharges in the downstream reaches of the river. In addition, a quantitative estimation of the uncertainties associated with this guess is arrived at.The results show that the extent and range of uncertainty is large with respect to the best guess changes. The uncertainty range is 2–3 times larger for the Business-as-Usual than for the Accelerated Policies scenarios. This large range stems from the doubtful precipitation simulations from the present General Circulation Models. This scenario study showed the precipitation scenarios to be the key-elements within the present range of reliable climate change scenarios.For the River Rhine best guess changes for annual water availability are small according to both scenarios. The river changes from a present combined snow-melt-rain fed river to an almost entirely rain fed river. The difference between present-day large average discharge in winter and the small average discharge in autumn should increase for all scenarios. This trend is largest in the Alpine part of the basin. Here, winter discharges should increase even for scenarios forecasting annual precipitation decreases. Summer discharge should decrease. Best guess scenarios should lead to increased frequencies of both low and high flow events in the downstream (Dutch) part of the river. The results indicate changes could be larger than presently assumed in worst case scenarios used by the Dutch water management authorities.  相似文献   

18.
Parameterizing turbulent diffusion through the joint probability density   总被引：3，自引：3，他引：0  
The convective mass flux parameterization often used in meteorological modeling expresses the vertical flux of a transported scalar as proportional to the product of the difference in mean values of the scalar in updrafts and downdrafts and their characteristic velocity. The proportionality factor is a constant to be specified. We show that this proportionality factor also appears in the relaxed eddy accumulation technique of Businger and Oncley. That associates the surface-layer flux of a scalar with the product of the standard deviation of vertical velocity and the mean concentration difference between updrafts and downdrafts.We show that this constant (b) is determined uniquely by the joint probability density (jpd) of vertical velocity and the scalar. Using large-eddy simulation, we generate this jpd for a conservative scalar diffusing through a convective boundary layer. It has quite different forms in top-down and bottom-up diffusion geometries. The bottom-up jpd is fairly well represented by a jointly Gaussian form and implies b ~ 0.6, in good agreement with the surface-layer value reported by Businger and Oncley. The top-down jpd is strikingly non-Gaussian and gives b ~ 0.47. Updrafts carry the bulk of the scalar flux - 70% in the bottom-up case, 60% in the top-down case.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

19.
The development of a theory of the evolution of the climate of the earth over millions of years can be subdivided into three fundamental, nested, problems: Firstly, to establish by equilibrium climate models (e.g., general circulation models) the diagnostic relations, valid at any time, between the fast-response climate variables (i.e., the weather statistics) and both the prescribed external radiative forcing and the prescribed distribution of the slow-response variables (e.g., the ice sheets and shelves, the deep ocean state, and the atmospheric CO2 concentration). Secondly, to construct, by an essentially inductive process, a model of the time-dependent evolution of the slow-response climatic variables over time scales longer than the damping times of these variables but shorter than the time scale of ultra-slow tectonic and astronomical changes in the boundary conditions (e.g., altered geography and elevation of the continents, slow outgassing and weathering and and solar radiative output). Thirdly, to determine the nature of these ultra-slow processes and their effects on the evolution of the equilibrium state of the climatic system about which the previously mentioned time-dependent variations occur. In this review we discuss the basis for this resolution, and give a broad overview of the contributions that have been made thus far in each area, emphasizing the work of the Yale climate group.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil  相似文献   

20.
Results from a global coupled ocean-atmosphere general circulation model (GCM) are used to perform the first in a series of studies of the various time and space scales of climate anomalies in an environment of gradually increasing carbon dioxide (CO2) (a linear transient increase of 1% per year in the coupled model). Since observed climate anomaly patterns often are computed as time-averaged differences between two periods, climate-change signals in the coupled model are defined using differences of various averaging intervals between the transient and control integrations. Annual mean surface air temperature differences for several regions show that the Northern Hemisphere warms faster than the Southern Hemisphere and that land areas warm faster than ocean. The high northern latitudes outside the North Atlantic contribute most to global warming but also exhibit great variability, while the high southern latitudes contribute the least. The equatorial tropics warm more slowly than the subtropics due to strong upwelling and mixing in the ocean. The globally averaged surface air temperature trend computed from annual mean differences for years 23–60 is 0.03 C per year. Projecting this trend to the time of CO2 doubling in year 100 produces a warming of 2.3° C. By chance, one particular northern winter five-year average geographical difference pattern in the Northern Hemisphere from the coupled model resembles the recent observed pattern of surface temperature and sea-level pressure anomalies. This pattern is not consistent from one five-year period to the next in any season in the model. However, multidecadal averages in the coupled model show that the North Atlantic warms less than the rest of the high northern latitudes, and recent observations may be a manifestation of this phenomenon. Consistent geographic patterns of climate anomalies forced by increased CO2 in the model are more evident with a longer averaging interval. There is also the possibility that the CO2 climate-change signal may itself be a function of time and space. The general pattern of zonal mean temperature anomalies for all periods in the model shows warming in the troposphere and cooling in the stratosphere. This pattern (or one similar to it taking into account the rest of the trace gases) could be looked for in observations to verify the enhanced greenhouse effect. A zonal mean pattern, however, could prove scientifically satisfactory but of little value to policymakers seeking regional climate-change forecasts. These results from the coupled model underscore the difficulty in identifying a time- and space-dependent fingerprint of greenhouse warming that has some practical use from short climatic records and point to the need to understand the mechanisms of decadal-scale variability.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

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

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