首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
In this paper SST imagery and a three-dimensional numerical model of a river plume were employed to detect upwelling induced by tidal straining in the Rhine ROFI (region of fresh water influence). Previous studies have shown that the Rhine ROFI in the North Sea exhibits strong cross-shore density gradients that compete with tidal and wind mixing to establish stratification. During neap periods with low mixing energy an area measuring 30 km offshore by 100 km alongshore becomes stratified. When the ROFI is stratified strong cross-shore currents are observed, with surface currents rotating anti-cyclonically and bottom currents rotating cyclonically. The cross-shore currents interact with the cross-shore density gradients to produce a semi-diurnal cycle of stratification. Due to continuity requirements imposed by the proximity of the coast, the offshore-directed surface currents and onshore-directed bottom currents should lead to coastal upwelling.  相似文献   

2.
The Basque coastal area, in the southeastern Bay of Biscay, can be characterised as being more influenced by land climate and inputs, than other typically ‘open sea’ areas. The influence of coastal processes, together with the presence of irregular and steep topography, complicate greatly the water circulation patterns. Water movement along the Basque coastal area is not well understood; observations are scarce and long-term current records are lacking. The knowledge available is confined to the surface currents: the surface water circulation is controlled mainly by wind forcing, with tidal and density currents being weak. However, there is a lack of knowledge available on currents within the lower levels of the water column; likewise, on the main time-scales involved in the water circulation. This study quantifies the contribution of the tidal and wind-induced currents, to the overall water circulation; it identifies the main time-scales involved within the tidal and wind-induced flows, investigating difference in such currents, throughout the water column, within Pasaia Bay (Basque coast). Within this context, extensive oceanographic and meteorological data have been obtained, in order to describe the circulation. The present investigation reveals that the circulation, within the surface and the sub-surface waters, is controlled mainly by wind forcing fluctuations, over a wide range of meteorological frequencies: third-diurnal, semidiurnal and diurnal land–sea breezes; synoptic variability; frequencies, near fortnightly periods; and seasonal. At the lower levels of the water column, the main contribution to the water circulation arises from residual currents, followed by wind-induced currents on synoptic time-scales. In contrast, tidal currents contribute minimally to the overall circulation throughout the water column.  相似文献   

3.
This article characterizes the spatial and temporal current variations, in the subtidal and tidal ranges, during the rainy and dry seasons, at the continental shelf off the Jaguaribe River, through measurements of continuous current field data from an acoustic Doppler current profiler (ADCP) mooring during 124 days, from June 12 to October 14, 2009. To support this dataset, we collected corresponding data from a meteorological station located at the estuary. The spatial variation showed that highest current speeds occur near the coast, with an offset of a NNW coastal jet, decreasing intensity, monotonically, towards offshore up to 0.1 ms?1. In the rainy season, small inversions of the wind field were observed, lasting 2 to 3 days on average and were accompanied by the direction of surface currents only. In the dry season, the period of reversal of wind fields and currents lasted 14 and 35 h, respectively. The analysis of empirical orthogonal functions in rainy and dry seasons showed that the continental shelf is predominantly barotropic, where the second and third modes explained only 7% of the total variance, during the dry season. The tidal currents are more intense in the direction normal to the coast, showing a semidiurnal tidal regime. Energy distribution between tidal currents and currents of longer periods showed that for the component parallel to bathymetry, subtidal frequency currents are dominant, contributing to more than 70% of the variance. For the normal component to the coastline bathymetry, there is a significant increase of power concerning tidal currents, at all depths, so they contribute with about 55% of the total variance.  相似文献   

4.
Tidal residual eddies and their effect on water exchange in Puget Sound   总被引:1,自引:0,他引:1  
Tidal residual eddies are one of the important hydrodynamic features in tidally dominant estuaries and coastal bays, and they could have significant effects on water exchange in a tidal system. This paper presents a modeling study of tides and tidal residual eddies in Puget Sound, a tidally dominant fjord-like estuary in the Pacific Northwest coast, using a three-dimensional finite-volume coastal ocean model. Mechanisms of vorticity generation and asymmetric distribution patterns around an island/headland were analyzed using the dynamic vorticity transfer approach and numerical experiments. Model results of Puget Sound show that a number of large twin tidal residual eddies exist in the Admiralty Inlet because of the presence of major headlands in the inlet. Simulated residual vorticities near the major headlands indicate that the clockwise tidal residual eddy (negative vorticity) is generally stronger than the anticlockwise eddy (positive vorticity) because of the effect of Coriolis force. The effect of tidal residual eddies on water exchange in Puget Sound and its subbasins was evaluated by simulations of dye transport. It was found that the strong transverse variability of residual currents in the Admiralty Inlet results in a dominant seaward transport along the eastern shore and a dominant landward transport along the western shore of the inlet. A similar transport pattern in Hood Canal is caused by the presence of tidal residual eddies near the entrance of the canal. Model results show that tidal residual currents in Whidbey Basin are small in comparison to other subbasins. A large clockwise residual circulation is formed around Vashon Island near entrance of South Sound, which can potentially constrain the water exchange between the Central Basin and South Sound.  相似文献   

5.
The south-flowing waters of the Kamchatka and Oyashio currents are key components of the western subarctic Pacific gyre. The dissipation of tidal energy in shallow and coastal regions of these currents and the attendant mixing are the important processes that affect the upper layer temperature and salinity. Examples of the impact of tidal currents on water temperature and salinity are the persistent tide-driven mixing around the Kashevarov and Kruzenshtern banks. The Kruzenshtern Bank is a shallow submarine bank stretching along the eastern continental slope of the Kuril Islands with the minimum depth of 86 m. Surface drifters observations are used to determine the characteristics of tidal currents and the circulation over these banks. New software that allows more versatility in the harmonic analysis is used for drifter’s data. The two banks have similar features. The variations in current velocities are dominated by the diurnal signals. The K1 and O1 tidal ellipses over the banks are the largest and clockwise. The enhanced tidal currents suggest that the formation of cold and saline water in summer is due to mixing of water column over the banks with intermediate waters. Variations of tidal ellipses over the bank may explain the formation of polynya at the western end of the Kashevarov Bank. We found that the 18.6-year lunar nodal cycle is a significant characteristic of salinity variation at the temperature minimum in the eddies eastward of the Boussole Strait over the period 1990–2015.  相似文献   

6.
In this paper, high-resolution wave, current and water depth fields derived by marine X-Band radar are presented for a coastal region of extreme tidal currents in the presence of inhomogeneous bathymetry at the south coast of New Zealand’s North Island. The current and water depth information for the presented location covers an area of approximately 13 km2 with a spatial resolution of 225 m and an update rate of 3 min. The sea state data provides a spatial representation of coastal effects like wave shoaling and refraction forced by bathymetry and current interaction. The near-surface current measurements about 3 km off the coast show expected tidal current pattern with maximum northwest/southeast current of 1.5–2 m/s alongshore. This is in agreement with currents from the RiCOM hydrodynamic model. The spatial resolution of the observed current field exhibits in addition small-scale current features caused by the influence of the local bathymetry. These data demonstrate the insight to be gained in complex, high-energy coastal situations through the use of high-resolution remote sensing techniques.  相似文献   

7.
This paper presents the development of a multiple‐station neural network for predicting tidal currents across a coastal inlet. Unlike traditional hydrodynamic models, the neural network model does not need inputs of coastal topography and bathymetry, grids, surface and bottom frictions, and turbulent eddy viscosity. Without solving hydrodynamic equations, the neural network model applies an interconnected neural network to correlate the inputs of boundary forcing of water levels at a remote station to the outputs of tidal currents at multiple stations across a local coastal inlet. Coefficients in the neural network model are trained using a continuous dataset consisting of inputs of water levels at a remote station and outputs of tidal currents at the inlet, and verified using another independent input and output dataset. Once the neural network model has been satisfactorily trained and verified, it can be used to predict tidal currents at a coastal inlet from the inputs of water levels at a remote station. For the case study at Shinnecock Inlet in the southern shore of New York, tidal currents at nine stations across the inlet were predicted by the neural network model using water level data located from a station about 70 km away from the inlet. A continuous dataset in May 2000 was used for the training, and another dataset in July 2000 was used for the verification of the neural network model. Comparing model predictions and observations indicates correlation coefficients range from 0·95 to 0·98, and the root‐mean‐square error ranges from 0·04 to 0·08 m s?1 at the nine current locations across the inlet. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

8.
The cooling water discharged from a power station may contain dissolved material, in addition to the heat rejected from the station turbines, which must be dispersed away from the power station in a manner which is environmentally acceptable. Mathematical methods are presented for use in the calculation of the rises in background concentration arising from the discharge of such material into lakes, rivers, estuaries or the sea. The calculations take account of factors including the replacement of the receiving water by freshwater, the available volume of water in the system, the rate at which the material may be degraded in the receiving water and some aspects of the design of the cooling water system. Particular calculations emphasize the different levels of concentration which may be expected from similar discharges into the different types of receiving water body. These examples suggest that for conditions representative of some present power station sites the concentration in a lake may be ten times that in an estuary, and that in an estuary may be ten times that in the sea. Of course, this is not a general result and the concentration level must be calculated for each particular case. The analyses also suggest what hydrographic factors must be measured during field surveys at a particular site if the mathematical methods presented here are to be used for the calculation of concentrations of any discharged material.  相似文献   

9.
To model currents in a fjord accurate tidal forcing is of extreme importance. Due to complex topography with narrow and shallow straits, the tides in the innermost parts of a fjord are both shifted in phase and altered in amplitude compared to the tides in the open water outside the fjord. Commonly, coastal tide information extracted from global or regional models is used on the boundary of the fjord model. Since tides vary over short distances in shallower waters close to the coast, the global and regional tidal forcings are usually too coarse to achieve sufficiently accurate tides in fjords. We present a straightforward method to remedy this problem by simply adjusting the tides to fit the observed tides at the entrance of the fjord. To evaluate the method, we present results from the Oslofjord, Norway. A model for the fjord is first run using raw tidal forcing on its open boundary. By comparing modelled and observed time series of water level at a tidal gauge station close to the open boundary of the model, a factor for the amplitude and a shift in phase are computed. The amplitude factor and the phase shift are then applied to produce adjusted tidal forcing at the open boundary. Next, we rerun the fjord model using the adjusted tidal forcing. The results from the two runs are then compared to independent observations inside the fjord in terms of amplitude and phases of the various tidal components, the total tidal water level, and the depth integrated tidal currents. The results show improvements in the modelled tides in both the outer, and more importantly, the inner parts of the fjord.  相似文献   

10.
An experimental study of temperature cycles and the heat budget in the Duplin River, a tidal creek bordered by extensive intertidal salt marshes, was carried out in late summer of 2003 and spring of 2004 near Sapelo Island on the central Georgia coast in the southeastern US. Three water masses are identified with differing temperature and salinity regimes, the characteristics of which are dictated by channel morphology, tidal communication with the neighboring sound, ground water hydrology, the extent of local intertidal salt marshes and side channels and the spring–neap tidal cycle (which controls both energetic mixing and, presumably, ground water input). For the first experiment, heat budgets are constructed for the upper (warmer) and lower (cooler) areas of the Duplin River showing the diminishing importance of tidal advection away from the mouth of the creek along with the concomitant increase in the importance of both direct atmospheric fluxes and of interactions with the marsh and side creeks. The second experiment, in the spring of 2004, reexamines the heat budget on seasonal and daily averaged scales revealing the decreased importance of advective fluxes relative to direct atmospheric fluxes on this scale but the constant importance of marsh/creek interactions regardless of time scale or season. Short period temperature fluctuations which affect larval development are examined and analogies are drawn to use heat to understand the marsh as a source of sediment, carbon and other nutrients.  相似文献   

11.
Kodaira  Tsubasa  Waseda  Takuji 《Ocean Dynamics》2019,69(11):1373-1385

The Kuroshio current is well known for generating cold wakes behind islands over Izu Ridge in Northwestern Pacific. Observational data from the geostationary Himawari-8 satellite for 2015–2017 revealed the occurrence of cold waters during the period when the Kuroshio current flows away from the islands. With a focus on tidal currents, this study presents an investigation of dynamical processes responsible for the formation of areas with low sea surface temperature (SST) through the adoption of a high-resolution numerical ocean model for an event that happened in July 2017. Areas with cold water emerged only when tidal currents are included in the numerical model. The model results indicate the cold surface waters are formed in the vicinity of the islands because of upwelling and vertical mixing. Qualitative features of the cold water formation for each island are found to depend on its size, topography, and ambient currents. Near Kozu Island, the tidal excursion is large enough to cause eddy shedding. These shed eddies are stirred by tidal currents to extend the surface cooling effect to wider areas. Near Hachijo Island, a persistent wake is formed by the ambient northward current. Inclusion of tidal currents destabilizes the wake, and consequently leads to the formation of a low SST area, although no clear eddy shedding is detected. The flow patterns around the islands are classified using an additional non-dimensional parameter, defined as the ratio between tidal excursion and island diameter.

  相似文献   

12.
The use of low levels of oxidants as antifouling agents in coastal and estuarine power station cooling water circuits can give rise to concerns over the likely production of chlorination byproducts (CBPs) including organohalogens, of largely cryptic identity, loading and impact. In order to review the current understanding of the degree of impact involved this paper describes a recent collaborative programme of work on CBPs within the wider context of a number of allied studies on power station cooling water entrainment and discharge effects. Bromoform was the single most important CBP found in coastal power station effluents. Although bromoform was found to accumulate in the liver of the sea bass Dicentrarchus labrax (L.) continuously exposed to a chlorinated effluent stream, no otherwise untoward impact on continuously exposed individuals was discerned. Populations of the common blue mussel Mytilus edulis L., likewise exposed, exhibited two responses: a high degree of genoplasticity and the evolution of stress proteins--although both may have been due to thermal rather than toxic stress. Experiments on plankton that simulated the conditions of entrainment within cooling water systems isolated the influence of oxidant usage, with mortality varying between taxa and life stage. The results suggest a very limited impact of oxidant use and the associated CBPs on receiving waters both in terms of plume toxicity or any more widespread ecotoxicological influence. Instead, it is in combination with applied temperature (DeltaT), pressure and the other stressors involved in entrainment within the power station cooling water circuits themselves that an impact is most marked.  相似文献   

13.
MASNUM wave-tide-circulation coupled numerical model (MASNUM coupled model, hereinafter) is de-veloped based on the Princeton Ocean Model (POM). Both POM and MASNUM coupled model are ap-plied in the numerical simulation of the upwelling off Yangtze River estuary and in Zhejiang coastal waters in summer. The upwelling mechanisms are analyzed from the viewpoint of tide, and a new mechanism is proposed. The study suggests that the tidally inducing mechanism of the upwelling in-cludes two dynamic aspects: the barotropic and the baroclinic process. On the one hand, the residual currents induced by barotropic tides converge near the seabed, and upwelling is generated to maintain mass conservation. The climbing of the residual currents along the sea bottom slope also contributes to the upwelling. On the other hand, tidal mixing plays a very important role in inducing the upwelling in the baroclinic sea circumstances. Strong tidal mixing leads to conspicuous front in the coastal waters. The considerable horizontal density gradient across the front elicits a secondary circulation clinging to the tidal front, and the upwelling branch appears near the frontal zone. Numerical experiments are de-signed to determine the importance of tide in inducing the upwelling. The results indicate that tide is a key and dominant inducement of the upwelling. Experiments also show that coupling calculation of the four main tidal constituents (M2, S2, K1, and O1), rather than dealing with the single M2 constituent, im-proves the modeling precision of the barotropic tide-induced upwelling.  相似文献   

14.
The turbulent kinetic energy dissipation rate, ε, in tidal seas is maximum at the bottom during full flood and during full ebb, i.e. when tidal currents are strongest. In coastal regions with tides similar to a Kelvin wave, this coincides with high water and low water. If there is a freshwater source at the coast, stratification in such a region will be most stable at high water and least at low water. Measurements of ε in the Rhine region of freshwater influence performed by previous studies have revealed bottom maxima at both high and low water. In addition, a maximum in the upper half of the water column was found around high water, which cannot be explained by tidal shear at the bottom, convective instabilities or wind mixing. This study investigates the dissipation rate and relevant physical properties in the Rhine region of freshwater influence by means of three-dimensional numerical simulations using the General Estuarine Transport Model and idealised conditions. The measurements are well reproduced; two distinct peaks of ε are evident in the upper layer shortly before and after high water. These maxima turn out to be due to strong peaks in the alongshore shear occurring when the fore- and the back-front of the plume transit the water column.  相似文献   

15.
An irregular mesh model of the west coast of Britain is used to examine the sensitivity of tidal residuals to mesh resolution in the region. Computed residuals are compared with earlier published results determined with a high resolution (1 km grid) finite difference model of the eastern Irish Sea. Initial calculations show that tidal residuals are largest in nearshore regions particularly in the vicinity of headlands. Local refinement of the mesh in these regions leads to a more detailed picture of the flow field, particularly adjacent to the coast. Although large scale offshore features of the flow can be resolved using the high resolution finite difference model, such an approach leads to a “stair case” representation of the coastal boundary with an adjacent near coastal region of spurious tidal residuals. By using an irregular mesh that follows the coast, this effect is removed. In the Mersey river region the tidal residual is resolved with a mesh resolution of 120 m, although calculations show that its distribution is particularly sensitive to small scale features of the topography. A variable mesh that can accurately represent the lateral variations in river width and details of topography in both the nearshore and estuarine environment appears essential in modelling the coastal spread of freshwater plumes from rivers and pollutants discharged into the near coastal environment.  相似文献   

16.
In the summer of 2001, a coastal ocean measurement program in the southeastern portion of San Pedro Bay, CA, was designed and carried out. One aim of the program was to determine the strength and effectiveness of local cross-shelf transport processes. A particular objective was to assess the ability of semidiurnal internal tidal currents to move suspended material a net distance across the shelf. Hence, a dense array of moorings was deployed across the shelf to monitor the transport patterns associated with fluctuations in currents, temperature and salinity. An associated hydrographic program periodically monitored synoptic changes in the spatial patterns of temperature, salinity, nutrients and bacteria. This set of measurements show that a series of energetic internal tides can, but do not always, transport subthermocline water, dissolved and suspended material from the middle of the shelf into the surfzone. Effective cross-shelf transport occurs only when (1) internal tides at the shelf break are strong and (2) subtidal currents flow strongly downcoast. The subtidal downcoast flow causes isotherms to tilt upward toward the coast, which allows energetic, nonlinear internal tidal currents to carry subthermocline waters into the surfzone. During these events, which may last for several days, the transported water remains in the surfzone until the internal tidal current pulses and/or the downcoast subtidal currents disappear. This nonlinear internal tide cross-shelf transport process was capable of carrying water and the associated suspended or dissolved material from the mid-shelf into the surfzone, but there were no observation of transport from the shelf break into the surfzone. Dissolved nutrients and suspended particulates (such as phytoplankton) transported from the mid-shelf into the nearshore region by nonlinear internal tides may contribute to nearshore algal blooms, including harmful algal blooms that occur off local beaches.  相似文献   

17.
A numerical simulation of circulation in the Columbia River estuary and plume during the summer of 2004 is used to explore the mixing involved as river water is transformed into shelf water. The model is forced with realistic river flow, tides, wind stress, surface heat flux, and ocean boundary conditions. Simulated currents and water properties on the shelf near the mouth are compared with records from three moorings (all in 72 m of water) and five CTD sections. The model is found to have reasonable skill; statistically significant correlations between observed and modeled surface currents, temperature, and salinity are all 0.42–0.72 for the mooring records. Equations for the tidally averaged, volume-integrated mechanical energy budget (kinetic and potential) are derived, with attention to the effects of: (i) Reynolds averaging, (ii) a time varying volume due to the free surface, and (iii) dissipation very close to the bottom. It is found that convergence of tidal pressure work is the most important forcing term in the estuary. In the far field plume (which has a volume 15 times greater than that of the estuary), the net forcing is weaker than that in the estuary, and may be due to either tidal currents or wind stress depending on the time period considered. These forcings lead to irreversible mixing of the stratification (buoyancy flux) that turns river water into shelf water. This occurs in both the plume and estuary, but appears to be more efficient (17% vs. 5%), and somewhat greater (4.2 MW vs. 3.3 MW), in plume vs. estuary. This demonstrates the importance of both wind and tidal forcing to watermass transformation, and the need to consider the estuary and plume as part of a single system.  相似文献   

18.
An analytical and a numerical model are used to understand the response of velocity and sediment distributions over Gaussian-shaped estuarine cross-sections to changes in tidal forcing and water depth. The estuaries considered here are characterized by strong mixing and a relatively weak along-channel density gradient. It is also examined under what conditions the fast, two-dimensional analytical flow model yields results that agree with those obtained with the more complex three-dimensional numerical model. The analytical model reproduces and explains the main velocity and sediment characteristics in large parts of the parameter space considered (average tidal velocity amplitude, 0.1–1 m s − 1 and maximum water depth, 10–60 m). Its skills are lower for along-channel residual flows if nonlinearities are moderate to high (strong tides in deep estuaries) and for transverse flows and residual sediment concentrations if the Ekman number is small (weak tides in deep estuaries). An important new aspect of the analytical model is the incorporation of tidal variations in the across-channel density gradient, causing a double circulation pattern in the transverse flow during slack tides. The gradient also leads to a new tidally rectified residual flow component via net advection of along-channel tidal momentum by the density-induced transverse tidal flow. The component features landward currents in the channel and seaward currents over the slopes and is particularly effective in deeper water. It acts jointly with components induced by horizontal density differences, Coriolis-induced tidal rectification and Stokes discharge, resulting in different along-channel residual flow regimes. The residual across-channel density gradient is crucial for the residual transverse circulation and for the residual sediment concentration. The clockwise density-induced circulation traps sediment in the fresher water over the left slope (looking up-estuary in the northern hemisphere). Model results are largely consistent with available field data of well-mixed estuaries.  相似文献   

19.
Two very high-frequency radars (VHFR) operating on the Opal coast of eastern English Channel provided a nearly continuous 35-day long dataset of surface currents over a 500 km2 area at 0.6–1.8 km resolution. Argo drifter tracking and CTD soundings complemented the VHFR observations, which extended approximately 25 km offshore. The radar data resolve three basic modes of the surface velocity variation in the area, that are driven by tides, winds and freshwater fluxes associated with seasonal river discharge. The first mode, accounting for 90% of variability, is characterized by an along-shore flow pattern, whereas the second and third modes exhibit cross-shore, and eddy-like structures in the current velocity field. All the three modes show the dominant semi-diurnal variability and low-frequency modulation by the neap-spring tidal cycle. Although tidal forcing provides the major contribution to variability of local currents, baroclinicity plays an important role in shaping the 3D velocity field averaged over the tidal cycle and may strongly affect tracer dynamics on larger time scales. An empirical orthogonal function (EOF) decomposition and a spectral rotary analysis of the VHFR data reveal a discontinuity in the velocity field occurring approximately 10 km offshore which was caused by the reversal in the sign of rotation of the current vector. This feature of local circulation is responsible for surface current convergence on ebb, divergence on flood and strong oscillatory vertical motion. Spectral analysis of the observed currents and the results of the Agro drifter tracking indicate that the line of convergence approximately follows the 30-m isobath. The most pronounced feature of the radar-derived residual circulation is the along-coast intensification of surface currents with velocity magnitude of 0.25 m/s typical for the Regions of Freshwater Influence (ROFI). The analysis has provided a useful, exploratory examination of surface currents, suggesting that the circulation off the Opal coast is governed by ROFI dynamics on the hypertidal background.  相似文献   

20.
A three-dimensional hydrodynamic model is used to investigate intra-tidal and spring–neap variations of turbulent mixing, stratification and residual circulation in the Chesapeake Bay estuary. Vertical profiles of salinity, velocity and eddy diffusivity show a marked asymmetry between the flood and ebb tides. Tidal mixing in the bottom boundary layer is stronger and penetrates higher on flood than on ebb. This flood–ebb asymmetry results in a north–south asymmetry in turbulent mixing because tidal currents vary out of phase between the lower and upper regions of Chesapeake Bay. The asymmetric tidal mixing causes significant variation of salinity distribution over the flood–ebb tidal cycle but insignificant changes in the residual circulation. Due to the modulation of tidal currents over the spring–neap cycle, turbulent mixing and vertical stratification show large fortnightly and monthly fluctuations. The stratification is not a linear function of the tidal-current amplitude. Strong stratification is only established during those neap tides when low turbulence intensity persists for several days. Residual circulation also shows large variations over the spring–neap cycle. The tidally averaged residual currents are about 50% stronger during the neap tides than during the spring tides.  相似文献   

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

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