首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
This study examines seasonal circulation, hydrography, and associated spatial variability over the inner shelf of the northern South China Sea (NSCS) using a nested-grid coastal ocean circulation model. The model external forcing consists of tides, atmospheric forcing, and open boundary conditions based on the global ocean circulation and hydrography reanalysis produced by the Hybrid Coordinate Ocean model. Five numerical experiments are conducted with different combinations of external forcing functions to examine main physical processes affecting the seasonal circulation in the study region. Model results demonstrate that the monthly mean circulation in the study region features the Guangdong Coastal Current (GCC) over coastal waters and the South China Sea Warm Current (SCSWC) in the offshore deep waters. The GCC produced by the model flows nearly southwestward in winter months and northwestward in summer months, which agrees with previous studies. The SCSWC flows roughly northeastward and is well defined in summer months. In winter months, by comparison, the SCSWC is superseded by the southwestward strong wind-driven currents. Analysis of model results in five different experiments demonstrates that the monthly mean circulation over coastal and inner shelf waters of the NSCS can be approximated by barotropic currents forced by the southwestward monsoon winds in winter months. In summer months, by comparison, the monthly mean circulation in the study region is affected significantly by baroclinic dynamics associated with freshwater runoff from the Pearl River and advection of warm and saline waters carried by the SCSWC over the NSCS.  相似文献   

2.
A new approach to understand the physical processes that govern internal variability of the large scale North Atlantic ocean circulation is outlined and current methods and results are reviewed. In this approach, based on the theory of dynamical systems, internal variability is viewed as arising through successive transitions when parameters are changed. The potential of the approach is demonstrated through analysesof solutions of intermediate complexity models of the wind-driven ocean circulation in the North Atlantic. In a quasi-geostrophic modelfor the flow in a rectangular basin with idealized wind forcing, the basic transitions are already found and physical mechanisms at work can be described in detail. Qualitatively, this transition behavior remains robust in more realistic models, having shallow water dynamics, realistic wind forcingand continental geometry, although patterns and time scales changethrough the model hierarchy. The relevance of the results is discussed inrelation to those of observations and of ocean general circulation models.  相似文献   

3.
The Pearl River Estuary (PRE) in South China's Guangdong Province is a subtropical estuary with highly irregular topography and dynamically complicated circulations. A nested-grid coastal circulation modelling system is used in this study to examine dynamic responses of the PRE to tides, meteorological forcing and buoyancy forcing. The nested-grid modelling system is based on the Princeton Ocean Model and consists of three downscaling subcomponents: including an outer-most model with a coarse horizontal resolution of ~7 km for simulating tidally forced and wind-driven surface elevations and depth-mean currents over the China Seas from Bohai Sea to the northern South China Sea and an innermost model with a fine resolution of ~1.2 km for simulating the 3D coastal circulation and hydrography over the PRE and adjacent coastal waters. Model results during the winter northeast monsoon surge in January and super typhoon Koryn in June of 1993 are used to demonstrate that the 3D coastal circulation and hydrographic distributions in the PRE are affected by tides, winds and buoyancy forcing associated with river discharge from the Pearl River with significant seasonal and synoptic variabilities.  相似文献   

4.
In this paper, we use the unstructured grid model SCHISM to simulate the thermohydrodynamics in a chain of baroclinic, interconnected basins. The model shows a good skill in simulating the horizontal circulation and vertical profiles of temperature, salinity, and currents. The magnitude and phases of the seasonal changes of circulation are consistent with earlier observations. Among the mesoscale and subbasin-scale circulation features that are realistically simulated are the anticyclonic coastal eddies, the Sebastopol and Batumi eddies, the Marmara Sea outflow around the southern coast of the Limnos Island, and the pathway of the cold water originating from the shelf. The superiority of the simulations compared to earlier numerical studies is demonstrated with the example of model capabilities to resolve the strait dynamics, gravity currents originating from the straits, high-salinity bottom layer on the shallow shelf, as well as the multiple intrusions from the Bosporus Strait down to 700 m depth. The warm temperature intrusions from the strait produce the warm water mass in the intermediate layers of the Black Sea. One novel result is that the seasonal intensification of circulation affects the interbasin exchange, thus allowing us to formulate the concept of circulation-controlled interbasin exchange. To the best of our knowledge, the present numerical simulations, for the first time, suggest that the sea level in the interior part of the Black Sea can be lower than the sea level in the Marmara Sea and even in some parts of the Aegean Sea. The comparison with observations shows that the timings and magnitude of exchange flows are also realistically simulated, along with the blocking events. The short-term variability of the strait transports is largely controlled by the anomalies of wind. The simulations demonstrate the crucial role of the narrow and shallow strait of Bosporus in separating the two pairs of basins: Aegean-Marmara Seas from one side and Azov-Black Seas from the other side. The straits of Kerch and Dardanelles provide sufficient interbasin connectivity that prevents large phase lags of the sea levels in the neighboring basins. The two-layer flows in the three straits considered here show different dependencies upon the net transport, and the spatial variability of this dependence is also quite pronounced. We show that the blocking of the surface flow can occur at different net transports, thus casting doubt on a previous approach of using simple relationships to prescribe (steady) outflow and inflow. Specific attention is paid to the role of synoptic atmospheric forcing for the basin-wide circulation and redistribution of mass in the Black Sea. An important controlling process is the propagation of coastal waves. One major conclusion from this research is that modeling the individual basins separately could result in large inaccuracies because of the critical importance of the cascading character of these interconnected basins.  相似文献   

5.
The western Baltic Sea infront of the German coast is a highly variable dynamical system, dominated by a complex and small-scale morphometry, the water exchange between the Baltic and North Seas, and driven by local wind. Neither data collection, nor satellite images or model simulations alone were able to explain the observed spatial patterns and transport processes. Therefore, all these methods were combined to explain the dynamical features and to systematise them according to the typical local wind pattern and time series. The aim was to develop an instrument for regional authorities which supports the interpretation of coastal water monitoring data and forms a basis for an improved monitoring strategy. Satellite data of sea surface temperature and ocean colour from the sensors NOAA-AVHRR and SeaWiFS were applied for synoptic investigations in the entire region and Landsat-7-ETM+ for regional studies. Model simulations were performed for the western Baltic using a 3D model MOM-3 and for the Szczecin Lagoon using 2D model FEMFLOW. For the first time, regional particularities in the coastal dynamical features and processes are derived for the main wind directions and for transitions between dominant wind situations west and east as derived from wind statistics. The simulated transport of particles released from different coastal and open sea sources indicate the affected areas during changing forcing conditions. The results support the interpretation of acquired coastal monitoring data as well as the assessment and optimisation of the monitoring programme.  相似文献   

6.
The impact of the choice of high-resolution atmospheric forcing on ocean summertime circulation in the Gulf of Lions (GoL; Mediterranean Sea) is evaluated using three different datasets: AROME (2.5 km, 1 h), ALADIN (9.5 km, 3 h), and MM5 (9 km, 3 h). A short-term ocean simulation covering a 3-month summer period was performed on a 400-m configuration of the GoL. The main regional features of both wind and oceanic dynamics were well-reproduced by all three atmospheric models. Yet, at smaller scales and for specific hydrodynamic processes, some differences became apparent. Inertial oscillations and mesoscale variability were accentuated when high-resolution forcing was used. Sensitivity tests suggest a predominant role for spatial rather than temporal resolution of wind. The determinant influence of wind stress curl was evidenced, both in the representation of a mesoscale eddy structure and in the generation of a specific upwelling cell in the north-western part of the gulf.  相似文献   

7.
We assess the performance of an eddy-recognizing numerical ocean model in simulating the pattern and variability of the hydrography in the Skagerrak/northern North Sea area. The model we use is a version of the widely used Princeton ocean model employing a terrain-following vertical coordinate. Results from a series of five multi-year simulations of the mesoscale response are described. The simulations differ in their representation of the lateral freshwater supply to the model ocean of which the first is a reference simulation. The next four are variations in which the river discharges and/or the Baltic outflow are given more realistic representations. For validation, we have used in situ hydrographic data. A novelty is that we use the concepts of freshwater height and potential energy anomaly as objective validation tools. We find that, in general, the model faithfully reproduces many of the observed hydrographic features including their mean patterns and their variance. Not surprisingly, we find that the Baltic outflow is by far the most significant freshwater source in terms of its influence on the hydrography in the area, a result corroborating earlier findings. The best validation is obtained when all freshwater supply is made as realistic as possible, in particular the Baltic outflow. We also find that the large scale cyclonic circulation and the location of fronts are robust characteristics of the Skagerrak/northern North Sea circulation given the impact changes in the freshwater input has on the hydrography. Finally, we find that a further exploration of the impact of the lateral open boundary forcing, e.g., the input of Atlantic water, is needed.  相似文献   

8.
Large-scale redistribution of sand by hydrodynamical processes in shelf seas is important for basin and coastal evolution on time scales of a thousand to tens of thousands of years. The influence of tides on the large-scale net sand-transport patterns in the North Sea has received much attention, but the influence of wind-driven flow and wind waves has hardly been investigated. Here, to establish the present-day situation and to develop a method that can also be used for palaeo-situations and forecasts for different sea levels, this influence is assessed for the present southern North Sea using a numerical flow model, a parametric wave model and a wave-averaged sand-transport formulation. Various forcing combinations are used to identify the dominant transport mechanisms: tides only, tides and wind, tides and waves, and combined tides, wind and waves. Wind forcing is applied in two ways to find an efficient, but still representative, method of incorporating this stochastic process: a statistical wind climatology and an observed time series. The results show that (i) the wind climatology yields a good approximation of the sand transport computed using the time series; (ii) wind-driven flow and waves only contribute significantly to the net sand transport by tides when acting together where tidal currents are small; and (iii) various combinations of forcings dominate the net sand transport in different regions of the southern North Sea: (a) tides dominate in the southern, middle and northwestern parts of the Southern Bight and in the region of The Wash; (b) tides, wind-driven flow and waves all are important in the northeastern part of the Southern Bight; and (c) wind-driven flow and waves dominate north of the Friesian Islands, in the German Bight and on the Dogger Bank. Qualitative comparison with observations shows good agreement.  相似文献   

9.
The baroclinic response of a stratified coastal embayment (Lunenburg Bay of Nova Scotia) to the observed wind forcing is examined using two numerical models. A linear baroclinic model based on the normal mode approach shows skill at reproducing the observed isotherm movements and sub-surface currents during a time of strong stratification in the bay. The linear model also shows that the isotherm movement in Lunenburg Bay is influenced by the wind forcing and propagation of baroclinic Kelvin waves from neighbouring Mahone Bay. The effects of nonlinearity and topography are investigated using a three-dimensional nonlinear coastal circulation model. The nonlinear model results demonstrate that the nonlinear advection terms generate a gyre circulation at the entrance of Lunenburg Bay, and the slope bottom topography at the mouth of the bay strengthens the sub-surface time-mean inflow on the southern side of the bay. A comparison of model-calculated currents in different numerical experiments clearly shows that baroclinicity plays a dominant role in the dynamics of wind-driven circulation in Lunenburg Bay.  相似文献   

10.
Coastal polynyas around Antarctica are the place of intense air–sea exchanges which eventually lead to the formation of high-salinity shelf waters (HSSW) over continental shelves. Here, the influence of atmospheric forcing on coastal polynyas in the Ross Sea is studied by contrasting the response of a regional ocean/sea-ice circulation model to two different atmospheric forcing sets. A first forcing (DFS3) is based on ERA40 atmospheric surface variables and satellite products. A second forcing (MAR) is produced on the basis of ERA40 with a dynamical downscaling procedure. As compared to DFS3, MAR forcing is shown to improve substantially the representation of small-scale patterns of coastal winds with stronger katabatic winds along the coast. The response of the ocean/sea-ice model to the two forcing sets shows that the MAR forcing improves substantially the geographical distribution of polynyas in the Ross Sea. With the MAR forcing, the polynya season is also shown to last longer with a greater ice-production rate. As a consequence, a greater flow of dense water out of the polynyas is found with the MAR forcing and the properties of HSSW are notably improved as compared to the DFS3 forcing. The factors contributing to the activity of Terra Nova Bay and Ross Ice Shelf polynyas in the model are studied in detail. The general picture that emerges from our simulations is that the properties of HSSW are mostly set by brine rejection when the polynya season resume. We found that coastal polynyas in the Ross Sea export about 0.4 Sv of HSSW which then flows along three separate channels over the Ross Shelf. A 6-month time lag is observed between the peak of activity of polynyas and the maximum transport across the sills in the channels with a maximum transport of about 1 Sv in February. This lag corresponds to the time it takes to the newly formed HSSW to spread from the polynya to the sills (at a speed of nearly 2 cm s−1).  相似文献   

11.
This paper examines the subtidal circulation and associated variability in the Gulf of St. Lawrence, the Scotian Shelf, and the Gulf of Maine (GSL-SS-GOM) in 1988–2004 based on results produced by a nested-grid shelf circulation model. The model has a fine-resolution child model (～ (1/12)°) embedded inside a coarse-resolution parent model (～ (1/4)°) of the northwest Atlantic. The combination of the semi-prognostic method and the spectral nudging method is used to reduce the model seasonal bias and drift. The child model reproduces the general features of the observed circulation and hydrography over the study region during the study period. The child model results demonstrate that the time-mean circulation in the GSL is affected by the time-mean atmospheric forcing and inflow through the Strait of Belle Isle. The temporal variability in atmospheric forcing affects the outflow through western Cabot Strait, which in turn affects the transport of the Nova Scotian Current and the gulf-wide cyclonic circulation in the GOM. The simulated seasonal variability of salinity in the top 30 m of the GSL-SS-GOM is mainly affected by the equatorward advection of low-salinity waters from the lower St. Lawrence Estuary to the GOM through the Scotian Shelf. The simulated intraseasonal variability of circulation in the GSL is affected by the variability in the estuarine circulation in response to the temporal variability in atmospheric forcing. On the Scotian Shelf, the intraseasonal variability is mainly driven by the variability of wind forcing and mesoscale and nonlinear dynamics over the shelf break and slope region. The interannual variability in the simulated temperature and salinity are spatially coherent in the intermediate waters in the GSL, which is caused partially by the local response to atmospheric variability and partially by variabilities over the southern Newfoundland Shelf that enter the GSL through the eastern Cabot Strait. By comparison, on the Scotian Shelf, the interannual variability of simulated circulation is affected by anomalies produced by the nonlinear dynamics which are advected equatorwards by the shelf break jet.  相似文献   

12.
We present an idealized network model for storm surges in the Wadden Sea, specifically including a time-dependent wind forcing (wind speed and direction). This extends the classical work by H.A. Lorentz who only considered the equilibrium response to a steady wind forcing. The solutions obtained in the frequency domain for the linearized shallow-water equations in a channel are combined in an algebraic system for the network. The velocity scale that is used for the linearized friction coefficient is determined iteratively. The hindcast of the storm surge of 5 December 2013 produces credible time-varying results. The effects of storm and basin parameters on the peak surge elevation are the subject of a sensitivity analysis. The formulation in the frequency domain reveals which modes in the external forcing lead to the largest surge response at coastal stations. There appears to be a minimum storm duration, of about 3–4 h, that is required for a surge to attain its maximum elevation. The influence of the water levels at the North Sea inlets on the Wadden Sea surges decreases towards the shore. In contrast, the wind shearing generates its largest response near the shore, where the fetch length is at its maximum.  相似文献   

13.
Using in situ, continuous, high frequency (8–16 Hz) measurements of velocity, suspended sediment concentration (SSC), and salinity, we investigate the factors affecting near-bed sediment flux during and after a meteorological event (cold front) on an intertidal flat in central San Francisco Bay. Hydrodynamic forcing occurs over many frequency bands including wind wave, ocean swell, seiching (500–1000 s), tidal, and infra-tidal frequencies, and varies greatly over the time scale of hours and days. Sediment fluxes occur primarily due to variations in flow and SSC at three different scales: residual (tidally averaged), tidal, and seiching. During the meteorological event, sediment fluxes are dominated by increases in tidally averaged SSC and flow. Runoff and wind-induced circulation contribute to an order of magnitude increase in tidally averaged offshore flow, while waves and seiching motions from wind forcing cause an order of magnitude increase in tidally averaged SSC. Sediment fluxes during calm periods are dominated by asymmetries in SSC over a tidal cycle. Freshwater forcing produces sharp salinity fronts which trap sediment and sweep by the sensors over short (∼30 min) time scales, and occur primarily during the flood. The resulting flood dominance in SSC is magnified or reversed by variations in wind forcing between the flood and ebb. Long-term records show that more than half of wind events (sustained speeds of greater than 5 m/s) occur for 3 h or less, suggesting that asymmetric wind forcing over a tidal cycle commonly occurs. Seiching associated with wind and its variation produces onshore sediment transport. Overall, the changing hydrodynamic and meteorological forcing influence sediment flux at both short (minutes) and long (days) time scales.  相似文献   

14.
The wind-driven circulation in the northwestern Pacific and the South China Sea (SCS) is simplified as a two-layer, quasi-geostrophic model in two rectangular basins connected by an idealized strait. This model is used to investigate the impact of the western boundary current (WBC) on the adjacent marginal sea. The variability of the circulation in the two basins is investigated with a high resolution and at low viscosity, which allows the numerical solution to resolve mesoscale eddy forcing. The model ocean is driven by the time-independent asymmetric wind stress acting on the idealized Pacific (large basin) only. Under the reference parameters used in this study, the WBC can intrude into the idealized SCS (small basin) in the form of a loop current, shedding eddies regularly. The rate of eddy shedding is nearly constant throughout the entire integration time of the model; however, the intensity of the eddy-shedding exhibits multiple timescale variability ranging from quasi-biennial to decadal timescale. A set of sensitivity experiments demonstrate that our results are robust against changes to model parameters and geometry. Multivariate spectral analysis is used to extract the spatiotemporal feature of the variability. Joint analysis for the two basins shows that the circulation in the idealized SCS is significantly impacted by the variability at decadal (15-year), interannual (5–7-year and quasi-biennial), and quasi-annual timescales. The spatial structures of the modes of variablility suggest that the variance in position of the WBC, combined with mesoscale activity, act to influence the low-frequency modes of the idealized SCS. The structural differences between the modes strongly impacting the idealized SCS and those having weak influence are also presented.  相似文献   

15.
The spatial resolution of wind forcing fields is critical for modeling ocean surface waves. We analyze here the performance of the non-hydrostatic numerical weather prediction system WRF-ARW (Weather Research and Forecasting) run with a 14-km resolution for hindcasting wind waves in the North Atlantic. The regional atmospheric model was run in the domain from 20° N to 70° N in the North Atlantic and was forced with ERA-Interim reanalysis as initial and boundary conditions in a spectral nudging mode. Here, we present the analysis of the impact of spectral nudging formulation (cutoff wavelengths and depth through which full weighting from reanalysis data is applied) onto the performance of the modeled 10-m wind speed and wind wave fields for 1 year (2010). For modeling waves, we use the third-generation spectral wave model WAVEWATCH III. The sensitivity of the atmospheric and wave models to the spectral nudging formulation is investigated via the comparison with reanalysis and observational data. The results reveal strong and persistent agreement with reanalysis data during all seasons within the year with well-simulated annual cycle and regional patterns independently of the nudging parameters that were tested. Thus, the proposed formulation of the nudging provides a reliable framework for future long-term experiments aiming at hindcasting climate variability in the North Atlantic wave field. At the same time, dynamical downscaling allows for simulation of higher waves in coastal regions, specifically near the Greenland east coast likely due to a better representation of the mesoscale atmospheric dynamics in this area.  相似文献   

16.
Rationalized by the observational circulation pattern in the upper ocean of the North Pacific, meridional friction term is first incorporated in a barotropic theoretical model of the wind-driven circulation. The governing potential vorticity equation thence has β term and wind stress curl term (the two of the Sverdrup balance), zonal friction term and meridional friction term. The analytical solution satisfactorily captures many important features of the wind-driven circulation in the North Pacific: Kuroshio, Oyashio, Kuroshio extension, North Equatorial Current, and especially the eastern boundary currents in the North Pacific, i.e. California current and Alaska current.  相似文献   

17.
A three-dimensional sigma coordinate numerical model with wetting and drying (WAD) and a Mellor–Yamada turbulence closure scheme has been used in an idealized island configuration to evaluate how tidally driven dynamics and mixing are affected by inundation processes. Comprehensive sensitivity experiments evaluate the influence of various factors, including tidal amplitudes (from 1- to 9-m range), model grid size (from 2 to 16 km), stratification, wind, rotation, and the impact of WAD on the mixing. The dynamics of the system involves tidally driven basin-scale waves (propagating anticlockwise in the northern hemisphere) and coastally trapped waves propagating around the island in an opposite direction. The evolutions of the surface mixed layer (SML) and the bottom boundary layer (BBL) under different forcing have been studied. With small amplitude tides, wind-driven mixing dominates and the thickness of the SML increases with time, while with large-amplitude tides, tidal mixing dominates and the thickness of the BBL increases with time. The inclusion of WAD in the simulations increases bottom stress and impacts the velocities, the coastal waves, and the mixing. However, the impact of WAD is complex and non-linear. For example, WAD reduces near-coast currents during flood but increases currents during ebb as water drains from the island back to the sea. The impacts of WAD, forcing, and model parameters on the dynamics are summarized by an analysis of the vorticity balance for the different sensitivity experiments.  相似文献   

18.
The long-term variability of the non-tidal circulation in Southampton Water, a partially mixed estuary, was investigated using 71-day acoustic Doppler current profiler (ADCP) time series. The data show evidence that the spring–neap tidal variability of the turbulent mixing modulates the strength of the non-tidal residual circulation, with subtidal neap tide surface flows reaching 0.12 m s–1 compared to <0.05 m s–1 at spring tides. The amplitude of the neap-tide events in this non-tidal circulation is shown to be related to a critical value of the tidal currents, illustrating the strong dependence on tidal mixing. The results suggest that the dominant mechanism for generating these neap-tide circulation events is the baroclinic forcing of the horizontal density gradient, rather than barotropic forcing associated with ebb-induced periodic stratification. While tidal turbulence is thought to be the dominant control on this gravitational circulation, there is evidence of the additional effect of wind-driven mixing, including the effects of wind fetch and possibly wave development with along-estuary winds being more efficient at mixing the estuary than across-estuary winds. Rapid changes in atmospheric pressure also coincided with fluctuations in the gravitational circulation. The observed subtidal flows are shown to be capable of rapidly flushing buoyant material out of the estuary and into the coastal sea at neap tides.Responsible Editor: Iris Grabemann  相似文献   

19.
This study examines the circulation and associated monthly-to-seasonal variability in the Caribbean Sea using a regional ocean circulation model. The model domain covers the region between 99.0 and 54.0°W and between 8.0 and 30.3°N, with a horizontal resolution of 1/6°. The ocean circulation model is driven by 6-hourly atmospheric reanalysis data from the National Center for Environmental Prediction and boundary forcing extracted from 5-day global ocean reanalysis data produced by Smith et al. (Mercator Newsletter 36:39–49, 2010), and integrated for 7 years. A comparison of model results with observations demonstrates that the regional ocean circulation model has skill in simulating circulation and associated variability in the study region. Analysis of the model results, as well as a companion model run that uses steady annual mean forcing, illustrates the role of Caribbean eddies for driving monthly-to-seasonal circulation variability in the model. It is found that vertically integrated transport between Nicaragua and Jamaica is influenced by the interaction between the density perturbations associated with Caribbean eddies and the Nicaraguan Ridge. The impact of Caribbean eddies squeezing through the Yucatan Channel is also discussed.  相似文献   

20.
An event based, long-term, climatological analysis is presented that allows the creation of coastal ocean atmospheric forcing on the coastal ocean that preserves both frequency of occurrence and event time history. An algorithm is developed that identifies individual storm event (cold fronts, warm fronts, and tropical storms) from meteorological records. The algorithm has been applied to a location along the South Atlantic Bight, off South Carolina, an area prone to cyclogenesis occurrence and passages of atmospheric fronts. Comparison against daily weather maps confirms that the algorithm is efficient in identifying cold fronts and warm fronts, while the identification of tropical storms is less successful. The average state of the storm events and their variability are represented by the temporal evolution of atmospheric pressure, air temperature, wind velocity, and wave directional spectral energy. The use of uncorrected algorithm-detected events provides climatologies that show a little deviation from those derived using corrected events. The effectiveness of this analysis method is further verified by numerically simulating the wave conditions driven by the characteristic wind forcing and comparing the results with the wave climatology that corresponds to each storm type. A high level of consistency found in the comparison indicates that this analysis method can be used for accurately characterizing event-based oceanic processes and long-term storm-induced morphodynamic processes on wind-dominated coasts.  相似文献   

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

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