首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.  相似文献   

2.

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.

  相似文献   

3.
Three prominent features of ocean surface turbulent heat fluxes (THF) trends during 1958–2013 are identified based on the Objectively Analyzed air-sea Fluxes (OAFlux) data set. The associated ocean-atmosphere dynamics changes are further investigated based on the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis. First, the THF are enhanced over the mid-latitude expansions of the subtropical western boundary currents (WBCs). An intensified oceanic heat transport, forced by stronger near-surface zonal wind, is likely to be the cause of such THF tendency. Second, the THF are reduced over the tropical eastern Pacific Ocean, which is primarily caused by the decreasing near-surface wind speed and sea surface temperature (SST), associated with a local coupled ocean-atmosphere cooling mode. Finally, the THF are reduced over the northern tropical Atlantic Ocean, which is attributed to the decreasing air-sea humidity and temperature differences as a result of the convergence of near-surface air and the divergence of ocean currents (upwelling).  相似文献   

4.
A high-resolution, regional coupled atmosphere–ocean model is used to investigate strong air–sea interactions during a rapidly developing extratropical cyclone (ETC) off the east coast of the USA. In this two-way coupled system, surface momentum and heat fluxes derived from the Weather Research and Forecasting model and sea surface temperature (SST) from the Regional Ocean Modeling System are exchanged via the Model Coupling Toolkit. Comparisons are made between the modeled and observed wind velocity, sea level pressure, 10 m air temperature, and sea surface temperature time series, as well as a comparison between the model and one glider transect. Vertical profiles of modeled air temperature and winds in the marine atmospheric boundary layer and temperature variations in the upper ocean during a 3-day storm period are examined at various cross-shelf transects along the eastern seaboard. It is found that the air–sea interactions near the Gulf Stream are important for generating and sustaining the ETC. In particular, locally enhanced winds over a warm sea (relative to the land temperature) induce large surface heat fluxes which cool the upper ocean by up to 2 °C, mainly during the cold air outbreak period after the storm passage. Detailed heat budget analyses show the ocean-to-atmosphere heat flux dominates the upper ocean heat content variations. Results clearly show that dynamic air–sea interactions affecting momentum and buoyancy flux exchanges in ETCs need to be resolved accurately in a coupled atmosphere–ocean modeling framework.  相似文献   

5.
利用2002年4月24日至6月20日在西沙海区进行的第三次南海海-气通量观测试验资料,采用涡相关法和TOGA COARE25b版本通量计算方案,计算了西南季风爆发前后海洋-大气间的通量交换,讨论了辐射、动量、感热通量、潜热通量、海洋热量净收支的时间变化特征及其与气象要素变化的关系.结果表明:西南季风爆发前后,太阳短波辐射、海面净辐射、潜热通量和海洋热量净收支变化特别强烈;通量变化受不同环境要素的影响:感热通量与海-气温差呈正相关关系,与气温呈明显的负相关关系.潜热通量与风速、海-气温差及海面水温均有正相关关系,其中与风速的关系最密切.动量通量(τ)主要随风速变化,它与风速(V)的关系可以表示为τ=000185V2-000559V+001248.  相似文献   

6.
A coupled ocean and boundary layer flux numerical modeling system is used to study the upper ocean response to surface heat and momentum fluxes associated with a major hurricane, namely, Hurricane Dennis (July 2005) in the Gulf of Mexico. A suite of experiments is run using this modeling system, constructed by coupling a Navy Coastal Ocean Model simulation of the Gulf of Mexico to an atmospheric flux model. The modeling system is forced by wind fields produced from satellite scatterometer and atmospheric model wind data, and by numerical weather prediction air temperature data. The experiments are initialized from a data assimilative hindcast model run and then forced by surface fluxes with no assimilation for the time during which Hurricane Dennis impacted the region. Four experiments are run to aid in the analysis: one is forced by heat and momentum fluxes, one by only momentum fluxes, one by only heat fluxes, and one with no surface forcing. An equation describing the change in the upper ocean hurricane heat potential due to the storm is developed. Analysis of the model results show that surface heat fluxes are primarily responsible for widespread reduction (0.5°–1.5°C) of sea surface temperature over the inner West Florida Shelf 100–300 km away from the storm center. Momentum fluxes are responsible for stronger surface cooling (2°C) near the center of the storm. The upper ocean heat loss near the storm center of more than 200 MJ/m2 is primarily due to the vertical flux of thermal energy between the surface layer and deep ocean. Heat loss to the atmosphere during the storm’s passage is approximately 100–150 MJ/m2. The upper ocean cooling is enhanced where the preexisting mixed layer is shallow, e.g., within a cyclonic circulation feature, although the heat flux to the atmosphere in these locations is markedly reduced.  相似文献   

7.
A regressive correction method is presented with the primary goal of improving ENSO simulation in regional coupled GCM. It focuses on the correction of ocean-atmosphere exchanged fluxes. On the basis of numerical experiments and analysis, the method can be described as follows: first, driving the ocean model with heat and momentum flux computed from a long-term observation data set; the pro-duced SST is then applied to force the AGCM as its boundary condition; after that the AGCM’s simula-tion and the corresponding observation can be correlated by a linear regressive formula. Thus the re-gressive correction coefficients for the simulation with spatial and temporal variation could be obtained by linear fitting. Finally the coefficients are applied to redressing the variables used for the calculation of the exchanged air-sea flux in the coupled model when it starts integration. This method together with the anomaly coupling method is tested in a regional coupled model, which is composed of a global grid-point atmospheric general circulation model and a high-resolution tropical Pacific Ocean model. The comparison of the results shows that it is superior to the anomaly coupling both in reducing the coupled model ‘climate drift’ and in improving the ENSO simulation in the tropical Pacific Ocean.  相似文献   

8.
An important part of the influence of the oceans on the atmosphere is through direct radiation, sensible heat flux and release of latent heat of evaporation, whereby all of these processes are directly related to the surface temperature of the oceans. A main effect of the atmosphere on the oceans is through momentum exchange at the air-ocean interface, and this process is directly related to the surface wind stress. The sea surface temperature (SST) and the surface wind stress are the two important components in the air-ocean system. If SST is given, a thermally forced boundary layer atmospheric circulation can be simulated. On the other hand, if the surface wind stress is given, the wind-driven ocean waves and ocean currents can be computed.The relationship between SST and surface wind is a coupling of the atmosphere and the oceans. It changes a one-way effect (ocean mechanically driven by atmosphere, or atmosphere thermally forced by oceans) into two-way air-sea interactions. Through this coupling the SST distribution, being an output from an ocean model, leads to the thermally forced surface winds, which feeds back into the ocean model as an additional forcing.Based on Kuo's planetary boundary layer model a linear algebraic equation is established to link the SST gradient with the thermally forced surface wind. The surface wind blows across the isotherms from cold to warm region with some deflection angle to the right (left) in the Northern (Southern) Hemisphere. Results from this study show that the atmospheric stratification reduces both the speed and the deflection angle of the thermally forced wind, however, the Coriolis' effect increases the wind speed in stable atmosphere (Ri>10–4) and increases the deflection angle.  相似文献   

9.
The effect of wind waves on water level and currents during two storms in the North Sea is investigated using a high-resolution Nucleus for European Modelling of the Ocean (NEMO) model forced with fluxes and fields from a high-resolution wave model. The additional terms accounting for wave-current interaction that are considered in this study are the Stokes-Coriolis force, the sea-state-dependent energy and momentum fluxes. The individual and collective role of these processes is quantified and the results are compared with a control run without wave effects as well as against current and water-level measurements from coastal stations. We find a better agreement with observations when the circulation model is forced by sea-state-dependent fluxes, especially in extreme events. The two extreme events, the storm Christian (25–27 October 2013), and about a month later, the storm Xaver (5–7 December 2013), induce different wave and surge conditions over the North Sea. Including the wave effects in the circulation model for the storm Xaver raises the modelled surge by more than 40 cm compared with the control run in the German Bight area. For the storm Christian, a difference of 20–30 cm in the surge level between the wave-forced and the stand-alone ocean model is found over the whole southern part of the North Sea. Moreover, the modelled vertical velocity profile fits the observations very well when the wave forcing is accounted for. The contribution of wave-induced forcing has been quantified indicating that this represents an important mechanism for improving water-level and current predictions.  相似文献   

10.
A sea spray generation function(SSGF)for bubble-derived droplets that takes into account the impact of wave state on whitecap coverage was presented in this study.By combining the new SSGF with a previous wave-state-dependent SSGF for spume droplets,an SSGF applicable to both bubble-derived and spume droplets that includes the impacts of wave state was obtained.The produced SSGF varies with surface wind as well as with wave development.As sea surface wind increases,more sea spray droplets are produced,resulting in larger SSGFs and volume fluxes.Meanwhile,under the same wind conditions,the SSGF is mediated by wave state,with larger SSGFs corresponding to older waves and larger windsea Reynolds numbers.The impact of wave state on sea spray heat flux was then estimated by applying this SSGF while considering the thermodynamic feedback process.Under given atmospheric and oceanic conditions,the estimated sea spray heat flux increases with wind speed,wave age,and windsea Reynolds number.  相似文献   

11.
Ocean–atmosphere coupling in the Humboldt Current System (HCS) of the Southeast Pacific is studied using the Scripps Coupled Ocean–atmosphere Regional (SCOAR) model, which is used to downscale the National Center for Environmental Prediction (NCEP) Reanalysis-2 (RA2) product for the period 2000–2007 at 20-km resolution. An interactive 2-D spatial smoother within the sea-surface temperature (SST)–flux coupler is invoked in a separate run to isolate the impact of the mesoscale (~50–200 km, in the oceanic sense) SST field felt by the atmosphere in the fully coupled run. For the HCS, SCOAR produces seasonal wind stress and wind stress curl patterns that agree better with QuikSCAT winds than those from RA2. The SCOAR downscaled wind stress distribution has substantially different impacts on the magnitude and structure of wind-driven upwelling processes along the coast compared to RA2. Along coastal locations such as Arica and Taltal, SCOAR and RA2 produce seasonally opposite signs in the total wind-driven upwelling transport. At San Juan, SCOAR shows that upwelling is mainly due to coastal Ekman upwelling transport, while in RA2 upwelling is mostly attributed to Ekman pumping. Fully coupled SCOAR shows significant SST–wind stress coupling during fall and winter, while smoothed SCOAR shows insignificant coupling throughout, indicating the important role of ocean mesoscale eddies on air–sea coupling in HCS. Coupling between SST, wind speed, and latent heat flux is incoherent in large-scale coupling and full coupling mode. In contrast, coupling between these three variables is clearly identified for oceanic mesoscales, which suggests that mesoscale SST affects latent heat directly through the bulk formulation, as well as indirectly through stability changes on the overlying atmosphere, which affects surface wind speeds. The SST–wind stress and SST–heat-flux couplings, however, fail to produce a strong change in the ocean eddy statistics. No rectified effects of ocean–atmosphere coupling were identified for either the atmospheric or oceanic mean conditions, suggesting that mesoscale coupling is too weak in this region to strongly alter the basic climate state.  相似文献   

12.
Coupled assimilation for an intermediated coupled ENSO prediction model   总被引:4,自引:0,他引:4  
Fei Zheng  Jiang Zhu 《Ocean Dynamics》2010,60(5):1061-1073
The value of coupled assimilation is discussed using an intermediate coupled model in which the wind stress is the only atmospheric state which is slavery to model sea surface temperature (SST). In the coupled assimilation analysis, based on the coupled wind–ocean state covariance calculated from the coupled state ensemble, the ocean state is adjusted by assimilating wind data using the ensemble Kalman filter. As revealed by a series of assimilation experiments using simulated observations, the coupled assimilation of wind observations yields better results than the assimilation of SST observations. Specifically, the coupled assimilation of wind observations can help to improve the accuracy of the surface and subsurface currents because the correlation between the wind and ocean currents is stronger than that between SST and ocean currents in the equatorial Pacific. Thus, the coupled assimilation of wind data can decrease the initial condition errors in the surface/subsurface currents that can significantly contribute to SST forecast errors. The value of the coupled assimilation of wind observations is further demonstrated by comparing the prediction skills of three 12-year (1997–2008) hindcast experiments initialized by the ocean-only assimilation scheme that assimilates SST observations, the coupled assimilation scheme that assimilates wind observations, and a nudging scheme that nudges the observed wind stress data, respectively. The prediction skills of two assimilation schemes are significantly better than those of the nudging scheme. The prediction skills of assimilating wind observations are better than assimilating SST observations. Assimilating wind observations for the 2007/2008 La Niña event triggers better predictions, while assimilating SST observations fails to provide an early warning for that event.  相似文献   

13.
利用中等复杂程度热带大气和海洋模式研究了热带太平洋和大西洋SST通过风应力桥梁的相互作用.利用1958~1998年NCEP分析的海表面温度场(SST)强迫大气模式得到的表面风应力与NCEP分析的同期热通量共同驱动海洋模式,作为控制试验;和控制试验平行,但强迫大气模式的SST在某一海盆取为多年气候平均值的试验作为敏感性试验;比较控制试验与敏感性试验模拟,则可反映风应力桥梁作用下热带某海盆SST异常对其他海盆的影响.结果表明,热带某一海盆SST暖(冷)异常总是引起局地海盆表面西部西(东)风异常和东部东(西)风异常;热带太平洋SST暖(冷)异常导致的该海盆东部表面东(西)风异常可以扩展到热带大西洋,从而导致热带大西洋SST冷(暖)异常;热带大西洋SST暖(冷)异常导致的该海盆西部表面西(东)风异常可以扩展到热带太平洋,从而导致热带太平洋SST暖(冷)异常.  相似文献   

14.
Abstract

In an ocean with a horizontal bottom where no wind is blowing it is shown that the spin (angular momentum) of the ocean is conserved. Thus, when energy is dissipated, at least one of three things will happen: i) Wave spectra may move towards lower frequencies. ii) The directional distribution may be changed towards long-crested waves. iii) Shear currents may be generated. By neglecting ii) and iii), the frequency shift of a spectrum is calculated due to molecular dissipation. When all energy transforming phenomena as e.g. wave breaking and turbulence generation are taken into account, the conservation of spin seems to be able to explain the frequency shift of wave spectra. In shallow water it is shown that there is energy transfer from the waves to shear currents.  相似文献   

15.
Here, we address the sediment dynamics in the Black Sea based on analysis of remote sensing data from the Medium Resolution Imaging Spectrometer and numerical simulations with Nucleus for European Modelling of the Ocean model. Boundary conditions consist of realistic meteorological forcing, including significant wave height generated by wave prediction model. A number of sensitivity runs was analysed with the aim to find the most suitable parameters governing sediment fluxes. The comparison between numerical simulations and remote sensing data gives credibility to the quality of simulations. The combined effect of wind waves and currents in the bed layer controls the sediment resuspension that appears to be the major basin-wide source of sediment. Sensitivity experiments included or excluded different forcing terms, e.g. sediment flux from rivers enable to determine the spatial extensions of different point sources. It is concluded that wind-wave forcing is manifested in the sediment dynamics through episodic high energy events contributing to the increase of horizontal sediment fluxes over the northwestern shelf. Both satellite images and numerical model simulations demonstrated that the penetration of suspended sediment into the basin interior was governed by the dynamics of coastal and open-ocean eddies. While fine sediment at sea surface could cross the continental slope propagating into the open ocean, coarser fractions follow the bottom and their penetration into the open ocean is limited. The conclusion is thus that the deposition patterns correlate with the specific shape of Black Sea topography, and the largest depositions are observed in the area of continental slope.  相似文献   

16.
As large-scale ocean circulation is a key regulator in the redistribution of oceanic energy, evaluating the multi-decadal trends in the western Pacific Ocean circulation under global warming is essential for not only understanding the basic physical processes but also predicting future climate change in the western Pacific. Employing the hydrological observations of World Ocean Atlas 2018(WOA18) from 1955 to 2017, this study calculated the geostrophic currents, volume transport and multidecadal trends for the North Equatorial Current(NEC), the North Equatorial Countercurrent(NECC), the Mindanao Current(MC), the Kuroshio Current(KC) in the origin and the New Guinea Coastal Undercurrent(NGCUC) within tropical western Pacific Ocean over multi-decades. Furthermore, this study examined the contributions of temperature and salinity variations. The results showed significant strengthening trends in NEC, MC and NGCUC over the past six decades, which is mainly contributed by temperature variations and consistent with the tendency in the dynamic height pattern. Zonal wind stress averaged over the western Pacific Ocean in the same latitude of each current represents the decadal variation and multi-decadal trends in corresponding ocean currents, indicating that the trade wind forcing plays an important role in the decadal trend in the tropical western Pacific circulation. Uncertainties in the observed hydrological data and trends in the currents over the tropical western Pacific are also discussed. Given that the WOA18 dataset covers most of the historical hydrological sampling data for the tropical western Pacific, this paper provides important observational information on the multi-decadal trend of the large-scale ocean circulation in the western Pacific.  相似文献   

17.
Shear- and convection-driven turbulence coexists with wind-generated surface gravity waves in the upper ocean. The turbulent Reynolds stresses in the oceanic mixed layer can therefore interact with the shear of the wave-generated Stokes drift velocity to extract energy from the surface waves and inject it into turbulence, thus augmenting the mean shear-driven turbulence. Stokes production of turbulence kinetic energy (TKE) is difficult to measure in the field, since it requires simultaneous measurement of the turbulent stress and the Stokes drift profiles in the water column. However, it is readily inferred using second moment closure models of the oceanic mixed layer provided: (1) wave properties are available, along with the usual water mass properties, and radiative and air–sea fluxes needed to drive the mixed layer model and (2) the model skill can be assessed by comparing the model results against the observed dissipation rates of TKE. Comprehensive measurements made during the Reynolds 2002 campaign in the Baltic Sea have made the estimation of Stokes production possible, and in this paper, we report on the effort and the conclusions reached. Measurements of air–sea exchange parameters and water mass properties during the campaign allowed a mixed layer model to be run and the turbulent stress in the water column to be inferred. Simultaneous wave spectrum measurements enabled Stokes drift profile to be deduced and wave breaking to be included in the model run, and the Stokes production of TKE in the water column estimated. Direct measurements of the TKE dissipation rate from an upward traversing microstructure profiler were used to assure that the model could reproduce the turbulent dissipation rate in the water column. The model results indicate that the Stokes production of TKE in the mixed layer is of the same order of magnitude as the shear production and must therefore be included in mixed layer models.  相似文献   

18.
Using reanalysis data, the role of initial signals in the tropical Pacific Ocean in predictions of negative Indian Ocean Dipole (IOD) events were analyzed. It was found that the summer predictability barrier (SPB) phenomenon exists in predictions, which is closely related to initial sea temperature errors in the tropical Pacific Ocean, with type-1 initial errors presenting a significant west-east dipole pattern in the tropical Pacific Ocean, and type-2 initial errors showing the opposite spatial pattern. In contrast, SPB-related initial sea temperature errors in the tropical Indian Ocean are relatively small. The initial errors in the tropical Pacific Ocean induce anomalous winds in the tropical Indian Ocean by modulating the Walker circulation in the tropical oceans. In the first half of the prediction year, the anomalous winds, combined with the climatological winds in the tropical Indian Ocean, induce a basin-wide mode of sea surface temperature (SST) errors in the tropical Indian Ocean. With the reversal of the climatological wind in the second half of the prediction year, a west-east dipole pattern of SST errors appears in the tropical Indian Ocean, which is further strengthened under the Bjerknes feedback, yielding a significant SPB. Moreover, two types of precursors were also identified: a significant west-east dipole pattern in the tropical Pacific Ocean and relatively small temperature anomalies in the tropical Indian Ocean. Under the combined effects of temperature anomalies in the tropical Indian and Pacific oceans, northwest wind anomalies appear in the tropical Indian Ocean, which induce a significant west-east dipole pattern of SST anomalies, and yield a negative IOD event.  相似文献   

19.
Both the tropical Indian and tropical Pacific Oceans are active atmosphere-ocean interactive regions with robust interannual variability, which also constitutes a linkage between the two basins in the mode of variability. Using a global atmosphereocean coupled model, we conducted two experiments(CTRL and PC) to explore the contributions of Indian Ocean interannual sea surface temperature(SST) modes to the occurrence of El Ni?o events. The results show that interannual variability of the SST in the Indian Ocean induces a rapid growth of El Ni?o events during the boreal autumn in an El Ni?o developing year. However, it weakens El Ni?o events or even promotes cold phase conversions in an El Ni?o decaying year. Therefore, the entire period of the El Ni?o is shortened by the interannual variations of the Indian Ocean SST. Specifically, during the El Ni?o developing years, the positive Indian Ocean Dipole(IOD) events force an anomalous Walker circulation, which then enhances the existing westerly wind anomalies over the west Pacific. This will cause a warmer El Ni?o event, with some modulations by ocean advection and oceanic Rossby and Kelvin waves. However, with the onset of the South Asian monsoon, the Indian Ocean Basin(IOB) warming SST anomalies excite low level easterly wind anomalies over the west tropical Pacific during the El Ni?o decaying years. As a result, the El Ni?o event is prompted to change from a warm phase to a cold phase. At the same time, an associated atmospheric anticyclone anomaly appears and leads to a decreasing precipitation anomaly over the northwest Pacific. In summary, with remote forcing in the atmospheric circulation, the IOD mode usually affects the El Ni?o during the developing years, whereas the IOB mode affects the El Ni?o during the decaying years.  相似文献   

20.
Cheng  Shukui  Cao  Anzhou  He  Hailun  Li  Peiliang  Song  Jinbao 《Ocean Dynamics》2023,73(7):449-461
Ocean Dynamics - Sea surface temperature (SST) cooling is a typical ocean response to tropical cyclones (TCs). Previous studies have shown that TC-induced SST cooling is influenced by the TC...  相似文献   

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

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