首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 250 毫秒
1.
Seasonal and interannual variability of ocean bottom pressure(OBP) in the Southern Ocean was investigated using Gravity Recovery and Climate Experiment(GRACE) data and a Pressure Coordinate Ocean Model(PCOM)based on mass conservation. By comparing OBP, steric sea level, and sea level, it is found that at high latitudes the OBP variability dominates the sea level variability at seasonal-to-decadal time scales. The diagnostic OBP based on barotropic vorticity equation has a good correlation with t...  相似文献   

2.
The Southern Ocean plays an important role in the global overturning circulation as a significant proportion of deep water is converted into intermediate and deeper water masses in this region. Recently, a secular trend has been reported in wind stress around the Southern Ocean and it is thought theoretically that the strength of the ACC is closely related to wind stress, so one consequence should be a corresponding increase in ACC transport and hence changes in the rate of the global overturning. There are no long-term data sets of ACC transport and so we must examine other data that may also respond to changing wind stress. Here we calculate surface currents in Drake Passage every seven days over 11.25 years from 1992 to 2004. We combine surface velocity anomalies calculated from satellite altimeter sea surface heights with measured surface currents. Since 1992, the UK has regularly occupied WOCE hydrographic section SR1b across the ACC in Drake Passage. From seven hydrographic sections surface currents are estimated by referencing relative geostrophic velocities from CTD sections with current measurements made by shipboard and lowered acoustic Doppler current profilers. Combining the seven estimates of surface currents with the altimeter data reduces bias in the estimates of average currents over time through Drake Passage and we show that surface current anomalies estimated by satellite and in situ observations are in good agreement. The strongest surface currents are found in the Subantarctic and Polar Fronts with average speeds of 50 cm/s and 35 cm/s, respectively and are inversely correlated, so that maximum westward flow in one corresponds to minimum westward flow in the other. The average cross-sectional weighted surface velocity from 1992 to 2004 is 16.7 ± 0.2 cm/s. A spectral analysis of the average surface current has only weakly increasing energy at higher frequencies and there is no dominant mode of variability. The standard deviation of the seven day currents is 0.68 cm/s and a running 12 month average has only a slightly smaller standard deviation of 0.52 ± 0.16 cm/s. The southern annular mode (SAM) measures the circumpolar average of wind stress and like the surface currents its spectrum has slightly increased energy at frequencies greater than 1 cpy. A cospectral analysis of these, averaging cospectra of five slightly overlapping 36 month segments improve statistical reliability, suggests that there is coherence between them at 1 cpy with the currents leading changes in the Southern annular mode. We conclude that the SAM and average Drake Passage surface currents are weakly correlated with no dominant co-varying modes, and hence predicting Southern Ocean transport variability from the SAM is not likely to give significant results and that secular trends in surface currents are likely to be masked by weekly and interannual variability.  相似文献   

3.
Observed subsurface signature of Southern Ocean sea level rise   总被引:1,自引:0,他引:1  
Satellite altimetry data show a strong increase in sea level in various parts of the Southern Ocean over the 1990s. In this paper, we examine the causes of the observed sea level rise in the region south of Australia, using 13 years of repeat hydrographic data from the WOCE SR3 sections, and the SURVOSTRAL XBT and surface salinity data. The hydrographic data show a poleward shift in the position of the Subtropical and the Subantarctic Fronts over the period. In the Antarctic Zone, the Antarctic Surface Water has become warmer and fresher, and the Winter Water tongue has become warmer, fresher, thinner and shallower. Increased freshening south of the Polar Front is linked to increased precipitation over the 1990s. Temperature changes over the upper 500 m account for only part of the altimetric sea level rise. The CTD sections show that the deeper layers are also warmer and slightly saltier and the observed sea level can be explained by steric expansion over the upper 2000 m. ENSO variability impacts on the northern part of the section, and a simple Sverdrup transport model shows how large-scale changes in the wind forcing, related to the Southern Annular Mode, may contribute to the deeper warming to the south.  相似文献   

4.
Using observations and numerical simulations, this study examines the intraseasonal variability of the surface zonal current(u ISV) over the equatorial Indian Ocean, highlighting the seasonal and spatial differences, and the causes of the differences. Large-amplitude u ISV occurs in the eastern basin at around 80°–90°E and near the western boundary at 45°–55°E. In the eastern basin, the u ISV is mainly caused by the atmospheric intraseasonal oscillations(ISOs), which explains 91% of the standard...  相似文献   

5.
Coal-fired power stations in Collie, Western Australia emit 10 million tonnes of CO2 per year. This study assesses the potential opportunities of geological storage of CO2 both within the Collie Basin and the onshore part of the adjacent Southern Perth Basin of Western Australia within 50 km of Collie town site through a desktop evaluation of existing data. The aquifers and coal formations within both basins have been evaluated for their suitability for storage based on geological, geographical and environmental criteria related to storage capacity, injectivity, proximity to sources of CO2, location of other natural resources and containment security. The study has concluded that there is limited scope for large-scale storage of CO2 within the Collie Basin. In addition the potential for storage within coals of either basin is not a viable solution. This assessment is based on published criteria for CO2 storage in sedimentary basins and coal-bearing formations.  相似文献   

6.
《Oceanologica Acta》1999,22(1):1-17
Sea surface elevation in the South China Sea is examined in the Topex/Poseidon altimeter data from 1992 to 1995. Sea level anomalies are smoothed along satellite tracks and in time with tidal errors reduced by harmonic analysis. The smoothed data are sampled every ten days with an along-track separation of about 40 km. The data reveal significant annual variations in sea level. In winter, low sea level is over the entire deep basin with two local lows centred off Luzon and the Sunda Shelf. In summer, sea level is high off Luzon and off the Sunda Shelf, and a low off Vietnam separates the two highs. The boundary between the Vietnam low and Sunda high coincides with the location of a jet leaving the coast of Vietnam described in earlier studies. Principal component analysis shows that the sea level variation consists mainly of two modes, corresponding well to the first two modes of the wind stress curl. Mode 1 represents the oscillation in the southern basin and shows little inter-annual variation. The mode 2 oscillation is weak in the southern basin and is strongest off central Vietnam. During the winters of 1992–1993 and 1994–1995 and the following summers, the wind stress curl is weak, and the mode 2 sea level variation in the northern basin is reduced, resulting in weaker winter and summer gyres. Weakening of the Vietnam low in summer implies diminishing of the eastward jet leaving the coast of Vietnam. The results are consistent with model simulations.  相似文献   

7.
A long term simulation of the barotropic circulation in the Río de la Plata estuary was performed with the aim of identifying the characteristic patterns of sea surface height (SSH) variability from synoptic to inter-annual time scales and their forcing mechanisms. Hamburg Shelf Ocean Model (HamSOM), forced by tides, monthly mean runoffs and 4-daily 10 m winds and surface atmospheric pressure from the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis was run. The solution was analyzed for the period 1965–2004. Inter-annual variability accounts for almost 10% of the variance. The first EOF mode of SSH variability on this time scale is associated with a mean anomaly of approximately 0.25 m at the upper estuary forced by both runoff and winds, which seems to be strongly associated with the ENSO cycles. Other two modes, with periodicities around 2.5 and 10 years were also found. Even though they are linked to weaker SSH anomalies, they are consistent with inter-annual modes of wind variability reported by other authors. Those modes are important, particularly if they act in phase, because they can provide a background for stronger surges. In contrast with the salinity field, SSH variability on seasonal time scales accounts for a very small percentage of variance and it is the combination of an annual and a semi-annual signal forced by winds and runoff, respectively. Approximately 90% of the variance is due to wind driven variability on sub-annual time scales. The most significant SSH anomalies in this band are associated with cyclogenetic events in the atmosphere, occurring either over Uruguay or over the Patagonian Shelf, whereas the strengthening or weakening of the semi-permanent South Atlantic anticyclone displays a relatively smaller influence. In agreement with previous publications, the estuary's spatial patterns in response to short-scale wind variability seem to be determined by wind direction more than by wind speed.  相似文献   

8.
Parameterizations of normal atmospheric modes (NAMs) and orographic gravity waves (OGWs) are implemented into the mechanistic general circulation model of the middle and upper atmosphere (MUA). Numerical experiments of sudden stratospheric warming (SSW) events are performed for climatological conditions typical for January and February using meteorological reanalysis data from the UK MET Office in the MUA model averaged over the years 1992–2011 with the easterly phase of quasi-biennial oscillation (QBO). The simulation shows that an increase in the OGW amplitudes occurs at altitudes higher than 30 km in the Northern Hemisphere after SSW. The OGW amplitudes have maximums at altitudes of about 50 km over the North American and European mountain systems before and during SSW, as well as over the Himalayas after SSW. At high latitudes of the Northern Hemisphere, significant (up to 50–70%) variations in the amplitudes of stationary planetary waves (SPWs) are observed during and after the SSW. Westward travelling NAMs have local amplitude maximums not only in the Northern Hemisphere, but also in the Southern Hemisphere, where there are waveguides for the propagation of these modes. Calculated variations of SPW and NAM amplitudes correspond to changes in the mean temperature and wind fields, as well as the Eliassen-Palm flux and atmospheric refractive index for the planetary waves, during SSW. Including OGW thermal and dynamical effects leads to an increase in amplitude (by 30–70%) of almost all SPWs before and during SSW and to a decrease (up to 20–100%) after the SSW at middle and high latitudes of the Northern Hemisphere.  相似文献   

9.
Limitations in sea surface salinity (SSS) observations and timescale separation methods have led to an incomplete picture of the mechanisms of SSS decadal variability in the tropical Pacific Ocean, where the El Niño Southern Oscillation (ENSO) dominates. Little is known regarding the roles of the North Pacific Gyre Oscillation (NPGO) and the Pacific Decadal Oscillation (PDO) in the large-scale SSS variability over the tropical basin. A self-organizing map (SOM) clustering analysis is performed on the intrinsic mode function (IMF) maps, which are decomposed from SSS and other hydrological fields by ensemble empirical mode decomposition (EEMD), to extract their asymmetric features on decadal timescales over the tropical Pacific. For SSS, an anomalous pattern appeared during 1997 to 2004, a period referred to as the anomalous late 1990s, when strong freshening prevailed in large areas over the southwestern basin and moderate salinization occurred in the western equatorial Pacific. During this period, the precipitation and surface currents were simultaneously subjected to anomalous fluctuations: the precipitation dipole and zonal current divergence along the equator coincided with the SSS increase in the far western equatorial Pacific, while the weak zonal current convergence in the southwestern basin and large-scale southward meridional currents tended to induce SSS decreases there. The dominant decadal modes of SSS and sea surface temperature (SST) in the tropical Pacific both resemble the NPGO but occur predominantly during the negative and positive NPGO phases, respectively. The similarities between the NPGO and Central Pacific ENSO (CP-ENSO) in their power spectra and associated spatial patterns in the tropics imply their dynamical links; the correspondence between the NPGO-like patterns during negative (positive) phases and the CP La Niña (CP El Niño) patterns for SSS is also discussed.  相似文献   

10.
This study uses the Climate Forecast System Reanalysis(CFSR) to investigate the responses of the Southern Hemisphere(SH) extratropical climate to two types of El Ni?o–Southern Oscillation(ENSO)—the eastern Pacific(EP) type and the central Pacific(CP) type in different seasons. The responses are denoted by the anomalies of climate variables associated with one-standard-deviation increase in the Ni?o3 or Ni?o4 index. The results show that in austral spring the differences in the ENSO-related anomaly(ERA) patterns of atmospheric circulation between the EP ENSO period(1979–1998) and CP ENSO period(1999–2010) are mainly associated with the change in the ENSO-PSA2 relationship. Such differences affect the ERA fields of surface air temperature and mixed layer temperature, and finally result in significant differences in sea-ice concentration anomalies in the Atlantic sector. In austral summer, significant correlation exists between the variations of SAM and both of the variations of Ni?o3 and Ni?o4 in 1979–1998, while the correlation between SAM and Ni?o4 disappears in 1999–2010. For all seasons, the strength of the climate ERAs depend on if there are close relationship between ENSO and the major climate variation modes of the SH extratropics. For the climate variables, the ERA patterns of surface air temperature are generally controlled by surface wind anomalies and mirrored by the mixed layer temperature anomalies. The mixed layer depth anomalies are primarily modulated by surface heat flux anomalies and occasionally by anomalous wind. There are strikingly strong anomalies of surface heat flux in the autumn of 1979–1998 related to the Ni?o3 variation, the period when there is only significant correlation between ENSO and PSA2. There are no evidence that the SH extratropical climate variability induced by Ni?o3 variations are stronger in the EP-ENSO period, and that variability induced by Ni?o4 variations are stronger in the CP-ENSO period.  相似文献   

11.
Data on bottom-water potential temperature, turbidity and current indications show that in the Southern Ocean west of the Kerguelen Plateau, Antarctic Bottom Water (AABW) of Weddell Sea origin spreads northwards from the Atlantic—Indian Basin in two directions: (1) AABW enters the Agulhas Basin through relatively deep areas in the Mid-Indian Ridge at 20–25°E and possibly at 35°E, and flows northwards into the Mozambique Basin as far as its northern limits; (2) a more easterly spreading path extends from the Atlantic—Indian Basin through the Crozet into the Madagascar, Mascarene, Somali and Arabian Basins. The passage in the western branch of the Indian Ridge for the AABW spreading from the Crozet into the Madagascar Basin appears to be at 29-26°S and 60–64°E.East of the Kerguelen Plateau in the South Indian Basin, the bottom water formed mainly along the Adélie Coast and Ross Sea travels west towards the Kerguelen Plateau and then parallel to it. This water finally flows eastwards hugging the Southeast Indian Ridge. Significant deviations from this general circulation pattern occur due to local topographic effects. Some AABW in the South Indian Basin exits through a passage at 120–125°E in the region of the Australian—Antarctic discordance in the Southeast Indian Ridge and enters the South Australian Basin and subsequently the Wharton Basin. This passage is clearly indicated by the northward extension of a cold, bottom-water tongue as shown by the temperature distribution in the region; the bottom-water effects in the passage are reflected in the high turbidity and current lineations on the sea floor.In the Southern Ocean basins, bottom-water turbidity is generally high, reflecting in part the strong bottom-water activity. The effects of AABW circulation on the sea floor—in the form of well-developed small- or large-scale current ripples and erosional/depositional features, manganese-nodule formations, and unconformities and reworking of sediments observed in cores — are also marked in these basins. Even though the AABW in the Wharton Basin is cold, its spreading effects on the sea floor are minimal in this basin in contrast to the basins west of the Mid-Indian Ridge at comparable latitudes.  相似文献   

12.
Interannual variability of the sea surface height (SSH) over the northeast Pacific Ocean is hindcast with a reduced-gravity, quasi-geostrophic model that includes linear damping. The model is forced with monthly Ekman pumping fields derived from the NCEP reanalysis wind stresses. The numerical solution is compared with SSH observations derived from satellite altimeter data and gridded at a lateral resolution of 1 degree. Provided that the reduced gravity parameter is chosen appropriately, the results demonstrate that the model has significant hindcast skill over interior regions of the basin, away from continental boundaries. A damping time scale of 2 to 3 years is close to optimal, although the hindcast skill is not strongly dependent on this parameter.A simplification of the quasi-geostrophic model is considered in which Rossby waves are eliminated, yielding a Markov model driven by local Ekman pumping. The results approximately reproduce the hindcast skill of the more complete quasi-geostrophic model and indicate that the interannual SSH variability is dominated by the local response to wind forcing. There is a close correspondence the two leading empirical orthogonal modes of the local model and those of the observed SSH anomalies. The latter account for over half of the variance of the interannual signal over the region.  相似文献   

13.
From late 1995 through early 2001, three major interannual climate events occurred in the tropical Pacific; the 1995–97 La Niña (LN), 1997–98 El Niño (EN), and 1998–2001 LN. We analyze atmospheric and upper oceanic anomalies in the northeast Pacific (NEP) during these events, and compare them to anomalies both elsewhere in the north and tropical Pacific, and to typical EN and LN anomaly patterns. The atmospheric and oceanic anomalies varied strongly on intraseasonal and interannual scales. During the 1995–97 LN and 1997–98 EN, the Northeast Pacific was dominated by negative SLP and cyclonic wind anomalies, and by upper ocean temperature and sea surface height (SSH) anomalies. The latter were positive along the North American west coast and in the NEP thermal anomaly pool (between Hawaii, Vancouver Island, and Baja California), and negative in the central north Pacific. This atmospheric/oceanic anomaly pattern is typical of EN. An eastward shift in the atmospheric teleconnection from east Asia created EN-like anomalies in the NEP during the 1995–97 LN, well before the 1997–98 EN had begun. The persistence of negative sea-level pressure (SLP) and cyclonic wind anomalies in the NEP during the 1997–98 EN intensified pre-existing upper oceanic anomalies. Atmospheric anomalies were shifted eastward during late 1996–early 1998, leading to a similar onshore shift of oceanic anomalies. This produced exceptionally strong positive upper ocean temperature and SSH anomalies along the west coast during the 1997–98 EN, and explains the unusual coastal occurrences of several species of large pelagic warm-water fishes. The growth and eastward shift of these pre-existing anomalies does not appear to have been linked to tropical Pacific EN anomalies until late 1997, when a clear atmospheric teleconnection between the two regions developed. Prior to this, remote atmospheric impacts on the NEP were primarily from east Asia. As the 1998–2001 LN developed, NEP anomalies began reversing toward the typical LN pattern. This led to predominantly negative SLP and cyclonic wind anomalies in the NEP, and upper ocean temperature and SSH anomalies that were mainly negative along the west coast and positive in the central north Pacific. The persistence of these anomalies into mid-2001, and a number of concurrent biological changes in the NEP, suggest that a decadal climate shift may have occurred in late 1998.During 1995–2001, NEP oceanic anomalies tracked the overlying atmospheric anomalies, as indicated by the maintenance of a characteristic spatial relationship between these anomalies. In particular, wind stress curl and SSH anomalies in the NEP maintained an inverse relationship that strengthened and shifted eastward toward the west coast during late 1996–early 1998. This consistent relationship indicates that anomalous Ekman transport driven by regional atmospheric forcing was an important contributor to temperature and SSH anomalies in the NEP and CCS during the 1997–98 EN. Other studies have shown that coastal propagations originating from the tropical Pacific also may have contributed to coastal NEP anomalies during this EN. Our results indicate that at least some of this coastal anomaly signal may have been generated by regional atmospheric forcing within the NEP.  相似文献   

14.
An ocean general circulation model (OGCM) is used to identify a Southern Ocean southeast Pacific intrinsic mode of low frequency variability. Using CORE data a comprehensive suite of experiments were carried out to elucidate excitation and amplification responses of this intrinsic mode to low frequency forcing (ENSO, SAM) and stochastic forcing due to high frequency winds. Subsurface anomalies were found to teleconnect the Pacific and Atlantic regions of the Antarctic Circumpolar Current (ACC) thermocline. The Pacific region of the ACC is characterised by intrinsic baroclinic disturbances that respond to both SAM and ENSO, while the Atlantic sector of the ACC is sensitive to higher frequency winds that act to amplify thermocline anomalies propagating downstream from the Pacific. Non-stationary cluster analysis was used to identify the system’s dynamical regimes and characterise meta-stability, persistence and transitions between the respective states. This analysis reveals significant trends, indicating fundamental changes to the meta-stability of the ocean dynamics in response to changes in atmospheric forcing. Intrinsic variability in sea-ice concentration was found to be coupled to thermocline processes. Sea-ice variability localised in the Atlantic was most closely associated with high frequency weather forcing. The SAM was associated with a circumpolar sea-ice response whereas ENSO was found to be a major driver of sea-ice variability only in the Pacific. This simulation study identifies plausible mechanisms that determine the predictability of the Southern Ocean climate on multi-decadal timescales.  相似文献   

15.
We here investigate the frequency and intensity of oscillations in oceanographic data within intraseasonal time scales using spectral analysis of surface wind and wave time-series data collected at off-island weather stations or moored buoys around Taiwan. Data from marine weather stations were used to trace atmospheric conditions, while we used buoy data to examine sea states. The spectra and wavelet scalogram of the wind fields revealed oscillations with a period of around 20–33 days, and the energy density of the wind field at the off-island stations was stronger than that at the data buoy stations. However, the wavelet scalogram of the wave height measured at the buoy stations was stronger than its associated wind field. This long-period oscillation is consistent with the wavelet scalogram of the wind field calculated from the off-island weather stations. About 20–33 day oscillations exist within intraseasonal variations, which are closely linked to the atmospheric environment and to wind and ocean wave fields. Oscillations with a period of 5–10 days are a pronounced feature over northeastern Taiwan waters during the winter season and can be interpreted as the wave pattern following synoptic weather systems.  相似文献   

16.
A three-dimensional numerical model is used to simulate sea level and velocity variations in the South China Sea for 1992–1995. The model is driven by daily wind and daily sea surface temperature fields derived from the NCEP/NCAR 40-year reanalysis project. The four-year model outputs are analyzed using time-domain Empirical Orthogonal Functions (EOF). Spatial and temporal variations of the first two modes from the simulation compare favorably with those derived from satellite altimetry. Mode 1, which is associated with a southern gyre, shows symmetric seasonal reversal. Mode 2, which contributes to a northern gyre, is responsible for the asymmetric seasonal and interannual variations. In winter, the southern and northern cyclonic gyres combine into a strong basin-wide cyclonic gyre. In summer, a cyclonic northern gyre and an anticyclonic southern gyre form a dipole with a jet leaving the coast of Vietnam. Interannual variations are particularly noticeable during El Niño. The winter gyre is generally weakened and confined to the southern basin, and the summer dipole structure does not form. Vertical motions weaken accordingly with the basin-wide circulation. Variations of the wind stress curl in the first two EOF modes coincide with those of the model-derived sea level and horizontal velocities. The mode 1 wind stress curl, significant in the southern basin, coincides with the reversal of the southern gyre. The mode 2 curl, large in the central basin, is responsible for the asymmetry in the winter and summer gyres. Lack of the mode 2 contribution during El Niño events weakens the circulation. The agreement indicates that changes in the wind stress curl contribute to the seasonal and interannual variations in the South China Sea.  相似文献   

17.
南海海面高度季节变化的数值模拟   总被引:8,自引:1,他引:8  
比较POM模式模拟与观测(TOPEX/Poseidon高度计资料)的南海海面高度(SSH)的季节变化在空间分布上的一致性和差异.结果表明:本文使用的POM模式能较好地模拟南海SSH的季节变化;冬季与夏季,春季与秋季南海海面异常场形式完全相反,冬季Ekman输运造成在西海岸的堆积要比夏季在东海岸堆积更明显,而吕宋冷涡中心附近和吕宋海峡海面季节变化振幅最大;除春季以外,在南海绝大部分海域,海面高度的季节变化主要受风力的控制,南海海面热量通量对SSH的季节变化贡献约为20%,风应力对SSH的季节变化的贡献约为80%.  相似文献   

18.
地形与热源强迫下的南方涛动   总被引:1,自引:0,他引:1  
用El Nino和La Nina位相时的海温异常和地形作为大气下垫面的异常强迫,引入IAP的两层原始方程大气环流模式,模拟出了南方涛动的典型结构.当去掉地形后,仅仅由海温异常也能模拟出太平洋东西部的气压异常振荡,但太平洋东部振荡中心的位置并不与观测的一致.由此可见,观测到的南方涛动是在实际地形下对全球海温异常的响应.  相似文献   

19.
对海洋中起伏运动(heaving)信号的时空分布研究能够帮助我们更好地了解气候系统中的年际和年代际变率。文章通过再分析资料和模式对太平洋区域的heaving主要模态进行了研究。研究结果表明: 太平洋区域主要存在两种heaving模态: 第一模态主要表现为赤道东西两侧的温跃层异常信号反位相; 第二模态表现为赤道区域和副热带区域的温跃层异常信号呈现反位相变化的规律。本文对这两个主要heaving模态所涉及的物理过程进行详细讨论, 结果表明: 东西反位相模态主要是受赤道波动调节的结果; 而经向结构模态则主要是由赤道地区的波动和副热带区域的风应力旋度异常作用共同导致。此外, 我们还讨论了heaving模态可以通过海洋波动以及Ekman输送等过程对海盆尺度的热输送(振幅约为5×1014W)以及海洋热含量(振幅约为1.5×1020J)的再分配起到了关键的调制作用, 进一步表明heaving模态对全球气候变化有着重要的作用。  相似文献   

20.
SouthernOscillationforcedbyheatsourceandtopographyQianWeihong,YouXintian(ReceivedJune5,1996,acceptedDecember9,1996)Abstract──...  相似文献   

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

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