首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 343 毫秒
1.
Samples of surface snow were collected for stable isotope analysis along the traverse route from Zhongshan to Dome A (East Antarctica) from Dec 28th, 2007 to Feb. 8th, 2008. The local relationship between δD and surface temperature is established to be 6.4 ± 0.2 ‰ per °C, very similar to the average for East Antarctic. The deuterium excess shows a pattern of high values over Antarctica, particularly at Dome A. We compare our data with an atmospheric general circulation model which includes stable water isotopes (ECHAM5-wiso). The model simulation captures the right levels of δD, but overestimates δ18O. This study provides support for the ongoing deep ice core project at Dome A.  相似文献   

2.
The Hamburg atmospheric general circulation model (AGCM) ECHAM‐4 is used to identify the main source regions of precipitation falling on Greenland and Antarctica. Both water isotopes H218O and HDO are explicitly built into the water cycle of the AGCM, and in addition the capability to trace water from different source regions was added to the model. Present and LGM climate simulations show that water from the most important source regions has an isotopic signature similar to the mean isotope values of the total precipitation amount. But water from other source regions (with very different isotopic signatures) contributes an additional, non‐negligible part of the total precipitation amount on both Greenland and Antarctica. Analyses of the temperature‐isotope‐relations for both polar regions reveal a solely bias of the glacial isotope signal on Greenland, which is caused by a strong change in the seasonal deposition of precipitation originating from nearby polar seas and the northern Atlantic. Although the performed simulations under LGM boundary conditions show a decrease of the δ 18O values in precipitation in agreement with ice core measurements, the AGCM fails to reproduce the observed simultaneous decrease of the deuterium excess signal.  相似文献   

3.
ISOTOPE RECORDS FROM MONGOLIAN AND ALPINE ICE CORES AS CLIMATE INDICATORS   总被引：3，自引：0，他引：3  
The link between long term changes in the isotopic composition of precipitation and surface air temperature at a given location is of exceptional importance for paleoclimatic studies, as ahs been demonstrated by many recent publications based on the isotope records from polar ice cores. By means of direct comparison with instrumental data, this paper evaluates the potential of the deuterium and oxygen-18 records from two continental glaciers for monitoring climatic trends. The isotopic data presented characterize climatically contrasted enviroments. The records from the Swiss glacier show distinct seasonal variations. Oxygen-18 is fairly well correlated with the instrumental record of atmospheric temperature; the seasonal differences in deuterium excess reflect nearness to the oceanic moisture source. By contrast, the isotope data from the Mongolian site show poor correlation with atmospheric temperature. The seasonal variations in deuterium excess, with higher values during summer time, indicate that precipitation largely originates from re-evaporated continental moisture sources. In both cases however, the correlation with temperature is significantly improved by the elimination of values derived from years where major changes in seasonal distribution and/or snow loss obviously have occurred, thereby distoring the isotopic ratios for that particular year. Depending on the site selected for study, the stable isotope composition of ice cores should therefore be viewed not only as a proxy for atmospheric temperature, but also as an additional hydrometeorological parameter and source indicator for atmospheric moisture.  相似文献   

4.
Four Holocene-long East Antarctic deuterium excess records are used to study past changes of the hydrological cycle in the Southern Hemisphere. We combine simple and complex isotopic models to quantify the relationships between Antarctic deuterium excess fluctuations and the sea surface temperature (SST) integrated over the moisture source areas for Antarctic snow. The common deuterium excess increasing trend during the first half of the Holocene is therefore interpreted in terms of a warming of the average ocean moisture source regions over this time. Available Southern Hemisphere SST records exhibit opposite trends at low latitudes (warming) and at high latitudes (cooling) during the Holocene. The agreement between the Antarctic deuterium excess and low-latitude SST trends supports the idea that the tropics dominate in providing moisture for Antarctic precipitation. The opposite trends in SSTs at low and high latitudes can potentially be explained by the decreasing obliquity during the Holocene inducing opposite trends in the local mean annual insolation between low and high latitudes. It also implies an increased latitudinal insolation gradient that in turn can maintain a stronger atmospheric circulation transporting more tropical moisture to Antarctica. This mechanism is supported by results from a mid-Holocene climate simulation performed using a coupled ocean-atmosphere model. Received: 7 July 1999 / Accepted: 21 July 2000  相似文献   

5.
Snow accumulation and its moisture origin over Dome Argus, Antarctica   总被引：1，自引：0，他引：1  
The spatial and temporal variability of snow accumulation near Dome Argus, Antarctica, is assessed using new snow pit and stake measurement data together with existing snow pit, ice core and automatic weather station records. Snow accumulation rate shows large inter-annual variations, but stable multi-decadal levels over the last seven centuries. Spatial variations in snow accumulation within the space of 50 km of Dome Argus are relatively small, probably thanks to the smooth topography. A comparison of theses accumulation observations with ECMWF reanalyses (ERA-40 and ERA-Interim) suggests ECMWF reanalysis captures the seasonal variations, but underestimates the overall snow accumulation at Dome Argus by ~50 %. The moisture sources for precipitation over Dome Argus are examined by means of a Lagrangian moisture source diagnostic, based on the tracing of specific humidity changes along air parcel trajectories, for the period 2000–2004 using operational ECMWF analysis data. Dome Argus mainly receives moisture from the mid-latitude (46 ± 4°S) South Indian Ocean, with a seasonal latitudinal shift of about 6°. Compared to other central East Antarctic deep ice core sites such as Dome F, Dome C, Vostok, and EPICA Dronning Maud Land, Dome Argus has a more southerly moisture origin, probably due to topographic influences on the moisture transport paths. These results have important implications for the interpretation of future ice cores at Dome Argus.  相似文献   

6.
A detailed climatic study of the Antarctic Peninsula from 1850 to 1980 has been carried out through the analysis of deuterium content in the snow layers of Dalinger Dome (James Ross Island, Antarctic Peninsula). It is based on the high correlation found between mean deuterium contents at this site and temperature data from stations within this region going back as far as April 1903 for the Argentine Orcadas station.The strong correlation between isotopes and temperatures first reveals a 1956 isotope reference for the region considered. Secondly, the isotope-temperature gradient is estimated at 4.5%. °C–1 for deuterium.After checking that the major temperature anomalies on the Antarctic Peninsula recorded since 1904 (according to available data) correspond to annual mean stable isotope peaks at Dalinger Dome, the amplitude of four prior anomalies are estimated in °C. Finally, a cooling of about 2 °C since 1850 is suggested for the region.  相似文献   

7.
Three east Antarctic ice cores (Dome B, EPICA-Dome C and Komsomolskaia) give evidence for a uniform dust input to the polar plateau during the last glacial maximum (LGM)/Holocene transition (20 to 10 kyr BP) and the 87Sr/86Sr versus 143Nd/144Nd isotopic signature of the mineral particles highlights a common provenance from southern South America at that time. However, the size distribution of dust from the three ice cores highlights important differences within the east Antarctic during the LGM and shows clearly opposite regional trends during the climatic transition. Between Dome B and Dome C the timing of these changes is also different. A geographical diversity also arises from the different phasing of the short-term (multi-secular scale) dust size oscillations that are superposed at all sites on the main trends of glacial to interglacial changes. We hypothesize the dust grading is controlled by size fractionation inresponse to its atmospheric pathway, either in terms of horizontal trajectory or in altitude of transport. Such mechanism is supported also by the dust size changes observed during a volcanic event recorded in Vostok ice. Ice core dust size data suggest preferential upper air subsidence over the EDC-KMS region and easier penetration of relatively lower air masses to the DB area during the LGM. At the end of the last glacial period and during the climatic transition the region of relatively higher subsidence progressively moved southward. The scenario proposed, supported also by the LGM/Holocene regional changes of snow accumulation, likely operates even at sub-millennial time scale.  相似文献   

8.
Mid-latitude winter atmospheric variability in the South Indian Ocean and southwest Pacific Ocean regions of the circum-Antarctic are reconstructed using sea-salt aerosol concentrations measured in the high resolution Law Dome (DSS) ice core from East Antarctica. The sea-salt aerosol concentration data, as sodium (Na), were measured at approximately monthly resolution spanning the past 700 years. Analyses of covariations between Na concentrations in Law Dome ice, and mean sea-level pressure (MSLP) and wind field data were conducted to define the mid-latitude and sub-Antarctic atmospheric circulation patterns associated with variations in Na delivery. High Na concentrations in Law Dome snow are associated with increased meridional aerosol transport from mid-latitude sources. The seasonal average Na concentration for early winter (May, June, July (MJJ)) is strongly correlated to the mid-latitude MSLP field in the South Indian and southwest Pacific Oceans, and southern Australian regions. In addition, the average MJJ Na concentrations display a strong association with the stationary Rossby wave number 3 circulation, and are anti-correlated to the Southern Annular Mode (SAM) index of climate variability: high (low) Na concentrations occurring during negative (positive) SAM phases. This observed relationship is used to derive a proxy record for early-winter MSLP anomalies and the SAM in the South Indian and southwest Pacific Ocean regions over the period 1300–1995 AD. The proxy SAM index from 1300 to 1995 AD shows pronounced decadal-scale variability throughout. The period after 1500 AD is marked by a tendency toward slower variations and a weakly-positive mean SAM (enhanced westerlies in the 50° to 65°S zone) compared to the early part of the record.  相似文献   

9.
Snow precipitation is the primary mass input to the Antarctic ice sheet and is one of the most direct climatic indicators, with important implications for paleoclimatic reconstruction from ice cores. Provenance of precipitation and the dynamic conditions that force these precipitation events at four deep ice core sites (Dome C, Law Dome, Talos Dome, and Taylor Dome) in East Antarctica were analysed with air mass back trajectories calculated using the Lagrangian model and the mean composite data for precipitation, geopotential height and wind speed field data from the European Centre for Medium Range Weather Forecast from 1980 to 2001. On an annual basis, back trajectories showed that the Atlantic-Indian and Ross-Pacific Oceans were the main provenances of precipitation in Wilkes Land (80%) and Victoria Land (40%), respectively, whereas the greatest influence of the ice sheet was on the interior near the Vostok site (80%) and in the Southwest Ross Sea (50%), an effect that decreased towards the coast and along the Antarctic slope. Victoria Land received snowfall atypically with respect to other Antarctica areas in terms of pathway (eastern instead of western), seasonality (summer instead of winter) and velocity (old air age). Geopotential height patterns at 500 hPa at low (>10 days) and high (2–6 days) frequencies during snowfall cycles at two core sites showed large positive anomalies at low frequencies developing in the Tasman Sea-Eastern Indian Ocean at higher latitudes (60–70°S) than normal. This could be considered part of an atmospheric blocking event, with transient eddies acting to decelerate westerlies in a split region area and accelerate the flow on the flanks of the low-frequency positive anomalies.  相似文献   

10.
The interpretation of stable isotopes in speleothems in terms of past temperature variability or precipitation rates requires a comprehensive understanding of the climatic factors and processes that influence the δ18O signal in the way through the atmosphere to the cave, where carbonate precipitates acquiring its final isotopic composition. This study presents for the first time in the Iberia Peninsula an integrated analysis of the isotopic composition of rainfall (δ18Op) during 2010–2012 years and, through a detailed monitoring survey, the transference of the primary isotopic signal throughout the soil and epikarst into the Molinos cave (Teruel, NE Spain). Both air temperature and amount of precipitation have an important effect on δ18Op values, clearly imprinting a seasonal variability modulated by an amount effect when rainfall events are more frequent or intense. Air mass history and atmospheric circulation influences are considered through the study of weather types, synoptic-scale climate patterns and large-scale atmospheric circulation indexes (North Atlantic Oscillation and Western Mediterranean Oscillation) revealing a dominant source effect on δ18Op values in this region where tropical North Atlantic and Western Mediterranean are the two moisture source regions. A delay of 2–3 months occurs between the dripwater oxygen isotopic composition (δ18Od) respect to δ18Op values as a consequence of large residence time in the epikarst. Limited calcite precipitates are found from winter to spring when δ18Od values are less negative and dripwater rates are constant. This study suggests that NE Iberian δ18Ocalcite proxy records are best interpreted as reflecting a combination of parameters, not just paleotemperature or paleorainfall and, if extending present-day situation towards the recent past, a biased signal towards winter values should be expected in Molinos speleothem records.  相似文献   

11.
As part of the United States’ contribution to the International Trans-Antarctic Scientific Expedition (ITASE), a network of precisely dated and highly resolved ice cores was retrieved from West Antarctica. The ITASE dataset provides a unique record of spatial and temporal variations of stable water isotopes (δ18O and δD) across West Antarctica. We demonstrate that, after accounting for water vapor diffusion, seasonal information can be successfully extracted from the ITASE cores. We use meteorological reanalysis, weather station, and sea ice data to assess the role of temperature, sea ice, and the state of the large-scale atmospheric circulation in controlling seasonal average water isotope variations in West Antarctica. The strongest relationships for all variables are found in the cores on and west of the West Antarctic Ice Sheet Divide and during austral fall. During this season positive isotope anomalies in the westernmost ITASE cores are strongly related to a positive pressure anomaly over West Antarctica, low sea ice concentrations in the Ross and Amundsen Seas, and above normal temperatures. Analyses suggest that this seasonally distinct climate signal is due to the pronounced meridional oriented circulation and its linkage to enhanced sea ice variations in the adjacent Southern Ocean during fall, both of which also influence local to regional temperatures.  相似文献   

12.
A regional atmospheric climate model with multi-layer snow module (RACMO2) is forced at the lateral boundaries by global climate model (GCM) data to assess the future climate and surface mass balance (SMB) of the Antarctic ice sheet (AIS). Two different GCMs (ECHAM5 until 2100 and HadCM3 until 2200) and two different emission scenarios (A1B and E1) are used as forcing to capture a realistic range in future climate states. Simulated ice sheet averaged 2 m air temperature (T2m) increases (1.8–3.0 K in 2100 and 2.4–5.3 K in 2200), simultaneously and with the same magnitude as GCM simulated T2m. The SMB and its components increase in magnitude, as they are directly influenced by the temperature increase. Changes in atmospheric circulation around Antarctica play a minor role in future SMB changes. During the next two centuries, the projected increase in liquid water flux from rainfall and snowmelt, together 60–200 Gt year?1, will mostly refreeze in the snow pack, so runoff remains small (10–40 Gt year?1). Sublimation increases by 25–50 %, but remains an order of magnitude smaller than snowfall. The increase in snowfall mainly determines future changes in SMB on the AIS: 6–16 % in 2100 and 8–25 % in 2200. Without any ice dynamical response, this would result in an eustatic sea level drop of 20–43 mm in 2100 and 73–163 mm in 2200, compared to the twentieth century. Averaged over the AIS, a strong relation between $\Updelta$ SMB and $\Updelta\hbox{T}_{2{\rm m}}$ of 98 ± 5 Gt w.e. year?1 K?1 is found.  相似文献   

13.
Three ice cores and a set of snow pit samples collected on James Ross Island, Antarctic Peninsula, in 1979, 1981 and 1991 have been analyzed for water stable isotope content D or 18O (isotopic temperature) and major chemical species. A reliable and detailed chronological scale has been established first for the upper 24.5 m of water equivalent (1990–1943) where various data sets can be compared, then extended down to 59.5 m of water equivalent (1847) with the aid of seasonal variations and the sulphate peak reflecting the 1883 Krakatoa volcanic eruption. At James Ross Island, sea-salt aerosol is generally produced by ice-free marine surfaces during the summer months, although some winter sea-salt events have been observed. For the upper part of the core (1990–1943), correlations (positive or negative) were calculated between isotopic temperature, chloride content (a sea-salt indicator), sea-ice extent, regional atmospheric temperature changes and atmospheric circulation. The D and chloride content correlation was then extended back to 1847, making it possible to estimate decadal sea-ice cover fluctuations over the study period. Our findings suggest that ice-core records from James Ross Island reflect the recent warming and sea-ice decrease trends observed in the Antarctic Peninsula area from the mid-1940s.  相似文献   

14.
The NASA/GISS GCM is used to estimate the evaporative contributions of several oceanic regions (defined by temperature) to Antarctica's July precipitation. Tracer diagnostics in the GCM suggest that the weighted average evaporative source temperature for Antarctic precipitation as a whole is about 12°C. The average source temperature for local precipitation there varies from 9° C to 14° C. To examine the effect of evaporative source on water isotope concentration, the GCM also follows a global deuterium (HDO) tracer and deuterium tracers evaporating from each oceanic region. The results suggest that although evaporative source temperature does affect the concentrations of the individual HDO tracers, differences in evaporative source do not explain the scatter in the roughly linear relationship between condensation temperature and isotope concentration. Offprint requests to: R Koster  相似文献   

15.
The isotope enabled atmospheric water balance model is applied to examine the spatial and temporal variations of δ18O in precipitation, amount effect and meteoric water lines (MWL) under four scenarios with different fractionation nature and surface evaporation inputs. The experiments are conducted under the same weather forcing in the framework of the water balance and stable water isotope balance. Globally, the spatial patterns of mean δ18O and global MWLs simulated by four simulation tests are in reasonably good agreement with the Global Network of Isotopes in Precipitation observations. The results indicate that the assumptions of equilibrium fractionation for simulating spatial distribution in mean annual δ18O and the global MWL, and kinetic fractionation in simulating δ18O seasonality are acceptable. In Changsha, four simulation tests all reproduce the observed seasonal variations of δ18O in precipitation. Compared with equilibrium fractionation, the depleted degree of stable isotopes in precipitation is enhanced under kinetic fractionation, in company with a decrease of isotopic seasonality and inter-event variability. The alteration of stable isotopes in precipitation caused by the seasonal variation of stable isotopes in vapour evaporated from the surface is opposite between cold and warm seasons. Four simulations all produce the amount effect commonly observed in monsoon areas. Under kinetic fractionation, the slope of simulated amount effect is closer to the observed one than other scenarios. The MWL for warm and humid climate in monsoon areas are well simulated too. The slopes and intercepts of the simulated MWLs decrease under kinetic fractionation.  相似文献   

16.
Scaling analysis shows that the mean thickness of an ice sheet depends on the product of two poorly known quantities, the ice viscosity and the net snow accumulation rate. We adjust the viscosity of an ice sheet in order to get a consistent value of this product for the present-day ice sheet volume and area given the net snow accumulation rate calculated by an atmospheric general circulation model (GCM). We then hold this artificial rheology constant in further numerical experiments. We hope that in doing so we can partially compensate for systematic GCM errors in simulating the snow accumulation rate, and, therefore, thickening/thinning of ice sheets will depend mostly on the tendency in the net accumulation change rather than on its absolute value. Using this approach, the response of the Greenland and Antarctic ice sheets to doubling CO2 concentration is simulated and the horizontal distribution of possible thickening/thinning of polar ice obtained. We find that, initially, the region of thickening ice is close to the area of increased snowfall rate, but later it significantly changes under the influence of internal ice flow dynamics. The sea-level changes predicted by our experiments agree with some empirical estimates. The sensitivity experiment with assigned basal sliding does not show significant changes in the large-scale ice topography, meaning, for example, that there is no indication of a possible disintegration of the West Antarctic ice sheet. At the same time, the regional thickening/thinning of ice (and consequently the sea-level change) depends strongly on processes at the ice sheet bottom.  相似文献   

17.
The South American Summer Monsoon (SASM) is a prominent feature of summertime climate over South America and has been identified in a number of paleoclimatic records from across the continent, including records based on stable isotopes. The relationship between the stable isotopic composition of precipitation and interannual variations in monsoon strength, however, has received little attention so far. Here we investigate how variations in the intensity of the SASM influence δ18O in precipitation based on both observational data and Atmospheric General Circulation Model (AGCM) simulations. An index of vertical wind shear over the SASM entrance (low level) and exit (upper level) region over the western equatorial Atlantic is used to define interannual variations in summer monsoon strength. This index is closely correlated with variations in deep convection over tropical and subtropical South America during the mature stage of the SASM. Observational data from the International Atomic Energy Agency-Global Network of Isotopes in Precipitation (IAEA-GNIP) and from tropical ice cores show a significant negative association between δ18O and SASM strength over the Amazon basin, SE South America and the central Andes. The more depleted stable isotopic values during intense monsoon seasons are consistent with the so-called ’‘amount effect‘’, often observed in tropical regions. In many locations, however, our results indicate that the moisture transport history and the degree of rainout upstream may be more important factors explaining interannual variations in δ18O. In many locations the stable isotopic composition is closely related to El Niño-Southern Oscillation (ENSO), even though the moisture source is located over the tropical Atlantic and precipitation is the result of the southward expansion and intensification of the SASM during austral summer. ENSO induces significant atmospheric circulation anomalies over tropical South America, which affect both SASM precipitation and δ18O variability. Therefore many regions show a weakened relationship between SASM and δ18O, once the SASM signal is decomposed into its ENSO-, and non-ENSO-related variance.  相似文献   

18.
Specific studies about the stable isotope composition (18O/16O and D/H) of atmospheric icy conglomerations are still scarce. The present work offers, for the first time, a very detailed analysis of oxygen and hydrogen isotopic signatures of unusually large ice conglomerations, or “megacryometeors”, that fell to the ground in Spain during January 2000. The hydrochemical analysis is based on the bulk isotopic composition and systematic selective sampling (deuterium isotopic mapping) of eleven selected specimens. δ18O and δD (V-SMOW) of all samples fall into the Meteoric Water Line matching well with typical tropospheric values. The distribution of the samples on Craig's line suggests either a variation in condensation temperature and/or different residual fractions of water vapour (Rayleigh processes). Three of the largest megacryometeors exhibited unequivocally distinctive negative values (δ18O = −17.2%0 and δD = −127 %0 V-SMOW), (δ18O = −15.6%0 and δD = −112%0 V-SMOW) and (δ18O = −14.4%0 and δD = −100%0 V-SMOW), suggesting an atmospheric origin typical of the upper troposphere. Theoretical calculations indicate that the vertical trajectory of growth was lower than 3.2 km. During the period in which the fall of megacryometeors occurred, anomalous atmospheric conditions were observed to exist: a substantial lowering of the tropopause with a deep layer of saturated air below, ozone depression and strong wind shear. Moreover, these large ice conglomerations occurred during non-thunderstorm conditions, suggesting an alternative process of ice growth was responsible for their formation.  相似文献   

19.
The mechanisms involved in the glacial inception are still poorly constrained due to a lack of high resolution and cross-dated climate records at various locations. Using air isotopic measurements in the recently drilled NorthGRIP ice core, we show that no evidence exists for stratigraphic disturbance of the climate record of the last glacial inception (∼123–100 kyears BP) encompassing Dansgaard–Oeschger events (DO) 25, 24 and 23, even if we lack sufficient resolution to completely rule out disturbance over DO 25. We quantify the rapid surface temperature variability over DO 23 and 24 with associated warmings of 10±2.5 and 16±2.5°C, amplitudes which mimic those observed in full glacial conditions. We use records of δ18O of O2 to propose a common timescale for the NorthGRIP and the Antarctic Vostok ice cores, with a maximum uncertainty of 2,500 years, and to examine the interhemispheric sequence of events over this period. After a synchronous North–South temperature decrease, the onset of rapid events is triggered in the North through DO 25. As for later events, DO 24 and 23 have a clear Antarctic counterpart which does not seem to be the case for the very first abrupt warming (DO 25). This information, when added to intermediate levels of CO2 and to the absence of clear ice rafting associated with DO 25, highlights the uniqueness of this first event, while DO 24 and 23 appear similar to typical full glacial DO events.  相似文献   

20.
In order to investigate Last Glacial Maximum and future climate, we “precalibrate” the intermediate complexity model GENIE-1 by applying a rejection sampling approach to deterministic emulations of the model. We develop ~1,000 parameter sets which reproduce the main features of modern climate, but not precise observations. This allows a wide range of large-scale feedback response strengths which generally encompass the range of GCM behaviour. We build a deterministic emulator of climate sensitivity and quantify the contributions of atmospheric (±0.93°C, 1σ) vegetation (±0.32°C), ocean (±0.24°C) and sea–ice (±0.14°C) parameterisations to the total uncertainty. We then perform an LGM-constrained Bayesian calibration, incorporating data-driven priors and formally accounting for structural error. We estimate climate sensitivity as likely (66% confidence) to lie in the range 2.6–4.4°C, with a peak probability at 3.6°C. We estimate LGM cooling likely to lie in the range 5.3–7.5°C, with a peak probability at 6.2°C. In addition to estimates of global temperature change, we apply our ensembles to derive LGM and 2xCO2 probability distributions for land carbon storage, Atlantic overturning and sea–ice coverage. Notably, under 2xCO2 we calculate a probability of 37% that equilibrium terrestrial carbon storage is reduced from modern values, so the land sink has become a net source of atmospheric CO2.  相似文献   

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

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