首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 42 毫秒
1.
An idealized numerical study of the influence of a tidal flow around an island has been undertaken with ROMS. The study focusses on coastal island wakes which are mainly controlled by elliptical tidal current flows on shallow shelves. This model is typical of some isolated continental shelf islands. The model is forced by a semi-diurnal barotropic inertia gravity wave imposed on the four open boundaries of a rectangular domain and its propagation results in an elliptical tidal flow within the domain in which the circular island lies. The influence of the surrounding island bathymetry and of the ellipse shape has been studied both in two and three dimensions. In the island vicinity, the residual circulation patterns over a tidal period show alongshore flow divergence along the major axis and convergence along the minor axis. A thin tidal ellipse (i.e. with a large ratio between major and minor axes) leads to strong eddy activity periods in the lee of the island during the flood and ebb phases, with eddy dissipation phases in between. By contrast, an almost round ellipse (axis ratio nearly 1) leads to vorticity filaments which continuously progress around the island without eddy shedding. The presence of a topographic slope in the vicinity of the island strengthens the eddy activity. This study suggests that the tidal current rotation favors the development of the eddy rotating in the same direction and weakens the development of the second eddy. In three dimensions with a surrounding bathymetry, an intense upwelling occurs in a large area in the lee of the island and the vertical velocities are stronger with thinner ellipses. With a flat bottom the vertical motions are almost fully generated by convergence and divergence of the secondary flow. With a varying bottom topography, the vertical motions come from a combination of this mechanism with convergence and divergence of the depth averaged flow.  相似文献   

2.
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.

  相似文献   

3.
Lagrangian drifters, moored acoustic Doppler current meters and hydrographic observations are combined with wind observations to describe the mean and variable nature of flow around Kapiti Island, New Zealand. Thirteen day-long deployments of up to six Lagrangian drifters show the mean flow is to the southwest, with evidence of stronger flows in the channel separating the island from the mainland, and an island wake in the lee of the island. Vortices in this island wake may be tidally driven. Scaling considerations suggest the flow is strong enough that tidal-generated vortices are shed on each tidal cycle. Both the drifters and mooring data suggest that the d’Urville Current around Kapiti Island has a significant wind-driven component. During north-westerlies, the drifters tend to hug the coast, and south-eastwards flows in the Rauoterangi Channel are accelerated. We suggest the observed correlation is the local expression of a South Taranaki basin scale response to the winds.  相似文献   

4.
A three-dimensional finite-element model is used to investigate the tidal flow around Rattray Island, Great Barrier Reef, Australia. Field measurements and visual observations show both stable eddies developing at rising and falling tide in the wake of the island. The water turbidity suggests intense upwelling able to carry bed sediments upwards. Based on previous numerical studies, it remains unclear at this point whether the most intense upwelling occurs near the centre of the eddies or off the island's tips, closer to the island. All these studies resorted to a very simple turbulence closure, with a zero-equation model whereby the coefficient of vertical viscosity is computed via an algebraic expression. In this work, we aim at studying the influence of the turbulence closure on model results, with emphasis on the prediction of vertical motions. The Mellor and Yamada level 2.5 closure scheme is used and an increase in the intensity of vertical transport is observed. This increase is partly explained by the fact that the Mellor and Yamada model takes into account the hysteresis effect in the time variation of turbulence variables. The influence of the advection of turbulence variables is estimated to be negligible. By a better representation of transient coastal phenomena, the Mellor and Yamada level 2.5 turbulence closure improves the model to a significant degree.  相似文献   

5.
A layered-aquifer model of groundwater occurrence in an atoll island was tested with a solute-transport numerical model. The computer model used, SUTRA, incorporates density-dependent flow. This can be significant in freshwater-saltwater interactions associated with the freshwater lens of an atoll island. Boundary conditions for the model included ocean and lagoon tidal variations. The model was calibrated to field data from Enjebi Island, Enewetak Atoll, and tested for sensitivity to a variety of parameters. This resulted in a hydraulic conductivity of 10 m day−1 for the surficial aquifer and 1000 m day−1 for the deeper aquifer; this combination of values gave an excellent reproduction of the tidal response data from test wells. The average salinity distribution was closely reproduced using a dispersivity of 0.02m. The computer simulation quantitatively supports the layered-aquifer model, including under conditions of density-dependent flow, and shows that tidal variations are the predominant driving force for flow beneath the island. The oscillating, vertical flow produced by the tidal variations creates an extensive mixing zone of brackish water. The layered-aquifer model with tidally driven flow is a significant improvement over the Ghyben-Herzberg-Dupuit model as it is conventionally applied to groundwater studies for many Pacific reef islands.  相似文献   

6.
Complex geophysical wake flows   总被引：1，自引：1，他引：0  
Idealized studies of island wakes often use a cylinder-like island to generate the wake, whereas most realistic studies use a close representation of the oceanic bathymetry immersed in a complex representation of the “ambient” geophysical flows. Here, a system of multiple islands was placed into numerical and experimental channels, in order to focus on the complexity of the archipelago wake, including (a) the influence of small neighboring islands and (b) the role of the island-shelf. The numerical geostrophic and stratified channel was built using a three-dimensional primitive equation model, considering a realistic representation of the Madeira archipelago bathymetry, with prescribed initial and boundary conditions. Results from the simulations show that the neighboring islands alter the near-field wake. Small eddies generated by the neighboring islands lead to destabilization of the shear layers of the larger island. Laboratory experiments carried out in the Coriolis rotating tank corroborated this near-field disruptive mechanism. The neighboring island perturbation effect was present whatever the direction of the incoming flow, but under different regimes. North–south wakes produced geostrophic eddies (≥ R d), whereas west–east wakes produced (exclusively) ageostrophic submesoscale eddies (< < R d) which traveled offshore with wave-like motion. The archipelago shelf contributed to the asymmetric vertical migration of oceanic vorticity. Cyclonic vorticity dominated the surface dynamics, whereas anticyclonic circulation prevailed at the bottom part of the linearly stratified upper layer. This study identifies several likely wake scenarios induced by the Madeira archipelago, and may serve as guide for future multiscale numerical studies and in situ campaigns.  相似文献   

7.
The problem of flow separation around islands is investigated using a dynamically adaptive finite element model to allow for resolution of the shear layers that form in the advent of separation. The changes in secondary circulation and vertical motion that occur in both attached and separated flows are documented, as is the degree of closure of the wake eddies. In the numerical experiments presented, the strongest motion always takes place at the sides of the idealised island, where flow curvature and shear act together to induce ascent. In contrast, it is the slower motion within the wake eddies that allow streamlines to extend from the bottom to the surface. We find no evidence for closure of the wake eddies. Rather, all of our separated experiments show that streamlines that pass through the eddies originate outside of the shear layers and frictional boundary layers on the upstream side of the idealised island. The numerical experiments demonstrate the potential for dynamically adaptive, unstructured meshes to resolve the separated shear layers that occur downstream of the idealised island, as well as the narrow boundary layers that form on the island itself.  相似文献   

8.
China's island and reef coasts are mainly distributed along the coasts of Zhejiang and Fujian provinces. These areas have favorable deep-water conditions and great development potential. The Zhoushan Archipelago is a typical island and reef coast. In the midst of the islands, there are many tidal channels with high water depth and high current velocity. Daishan Island is located on the northern part of the Zhoushan Archipelago. The sea north of Daishan Island can be characterized as having great water depth (generally between 10 and 20 m), high current velocity (the maximum measured value is about 2 m/s), and hyperconcentrated suspended sediment (annual average concentration under normal weather conditions is about 0.6 kg/m3). Sediment transport is influenced by the co-action of waves and currents. Land reclamation can ease the contradiction between the desire for development and the land resource shortage on island and reef coasts. Land reclamation generally starts from deposition-promoting projects on island and reef coasts with hyperconcentration of suspended sediment with the aim of reducing the reclamation cost. Based on the morphodynamic characteristics of the sea area near north of Daishan Island, a two-dimensional mathematical model was developed to analyze the co-action of current and waves and the model was verified for spring as well as neap tides in 2007 and 2012. The simulated values of tidal stages, tidal current velocities, tidal current directions, and suspended sediment concentration at 15 stations in 2007 and 9 stations in 2012 were in good agreement with the in-situ measured values. Furthermore, the simulated magnitude and pattern of erosion and sedimentation agreed well with the observations. This model was used to evaluate the effects of a proposed siltation scheme on the hydrodynamic environment and morphological change of the study area. The results of this study can inform the protection, development, and utilization of shoal-channels on island and reef coasts, as well provide a scientific basis for comprehensive maritime protection and development.  相似文献   

9.
We analyzed gravity data obtained in Juneau and global positioning system (GPS) data obtained from three PBO sites in southeastern Alaska (SE-AK), which are part of a US research facility called ‘EarthScope’, and we compared the obtained tidal amplitudes and phases with those estimated from the predicted tides including both effects of the body tide and ocean tide. Global tide models predict the ocean tides in this region of complex coastline and bathymetry. To improve the accuracy of prediction, we developed a regional ocean tide model in SE-AK.Our comparison results suggest: (1) by taking into account the ocean tide effect, the amplitude differences between the observation and the predicted body tide is remarkably reduced for both the gravity and displacement (e.g. for the M2 constituent, 8.5–0.3 μGal, and 2.4–0.1 cm at the AB50 GPS site in Juneau in terms of the vector sum of three components of the north–south, east–west and up–down), even though the ocean tide loading is large in SE-AK. (2) We have confirmed the precise point positioning (PPP) method, which was used to extract the tidal signals from the original GPS time series, works well to recover the tidal signals. Although the GPS analysis results still contain noise due to the atmosphere and multipath, we may conclude that the GPS observation surely detects the tidal signals with the sub-centimeter accuracy or better for some of the tidal constituents. (3) In order to increase the accuracy of the tidal prediction in SE-AK, it is indispensable to improve the regional ocean tide model developed in this study, especially for the phase.  相似文献   

10.
Two interdisciplinary cruises aimed at relating the ecology of marine fish populations to oceanographic conditions were fielded during the late summer and late winter seasons near Palmyra Atoll (5.9°N, 162.1°W) in the Line Islands. Ocean current and hydrographic measurements revealed interaction of the flow with the steep topography. During the first cruise (August/September 1990) satellite-tracked surface drifters and acoustic Doppler current profiler (ADCP) measurements showed a strong eastward setting North Equatorial Counter Current (NECC) with maximum speeds exceeding 1 m s–1 at 80 m depth approximately. This current turned southeastward on closer approach to Palmyra. The drifter paths exhibited excursions with zonal wavelength of approximately 250 km, meridional amplitude of 25 km and period of approximately 5 days. During the second cruise (February/March 1992), the ADCP-derived speeds of the NECC were weaker (maxima approximately 33 cm s–1) while the relative geostrophic flow component was of magnitude similar to 1990 and the signal of zonal geostrophic currents reached much deeper to approximately 650 m depth (150 m in 1990). Doming isopycnals beneath the surface mixed layer as well as thick (10–25 m) internal mixed layers were found near Palmyra during both cruises, with slightly different positions relative to the island. The discontinuous vertical temperature profiles may have been a result of strong boundary mixing due to breaking internal waves on Palmyras steep slopes. In the immediate vicinity of the island variations in flow speed, stratification and mixing in both the alongshore and cross-isobath directions were observed. Overall, the current speeds were reduced during February/March 1992, the peak time of the 1991–1993 warm event in the tropical Pacific. While parameters of turbulent two-dimensional wake theory are suggestive of formation and shedding of eddies in the lee of the island, no direct observations of circular motions were made in either expedition.Responsible Editor: Hans Burchard  相似文献   

11.
With the existence of eight substantial islands in the Southern California Bight, the oceanic circulation is significantly affected by island wakes. In this paper a high-resolution numerical model (on a 1 km grid), forced by a high-resolution wind (2 km), is used to study the wakes. Island wakes arise due both to currents moving past islands and to wind wakes that force lee currents in response. A comparison between simulations with and without islands shows the surface enstrophy (i.e., area-integrated square of the vertical component of vorticity at the surface) decreases substantially when the islands in the oceanic model are removed, and the enstrophy decrease mainly takes place in the areas around the islands. Three cases of wake formation and evolution are analyzed for the Channel Islands, San Nicolas Island, and Santa Catalina Island. When flows squeeze through gaps between the Channel Islands, current shears arise, and the bottom drag makes a significant contribution to the vorticity generation. Downstream the vorticity rolls up into submesoscale eddies. When the California Current passes San Nicolas Island from the northwest, a relatively strong flow forms over the shelf break on the northeastern coast and gives rise to a locally large bottom stress that generates anticyclonic vorticity, while on the southwestern side, with an adverse flow pushing the main wake current away from the island, positive vorticity has been generated and a cyclonic eddy detaches into the wake. When the northward Southern California Countercurrent passes the irregular shape of Santa Catalina Island, cyclonic eddies form on the southeastern coast of the island, due primarily to lateral stress rather than bottom stress; they remain coherent as they detach and propagate downstream, and thus they are plausible candidates for the submesoscale “spirals on the sea” seen in many satellite images. Finally, the oceanic response to wind wakes is analyzed in a spin-up experiment with a time-invariant wind that exhibits strips of both positive and negative curl in the island lee. Corresponding vorticity strips in the ocean develop through the mechanism of Ekman pumping.  相似文献   

12.
This work presents results from two complementary and interconnected approaches to study water temperature and salinity patterns in an estuarine tidal channel. This channel is one of the four main branches of the Ria de Aveiro, a shallow lagoon located in the Northwest coast of the Iberian Peninsula. Longitudinal and cross-sectional fields of water temperature and salinity were determined by spatial interpolation of field measurements. A numerical model (Mohid) was used in a 2D depth-integrated mode in order to compute water temperature and salinity patterns. The main purpose of this work was to determine the horizontal patterns of water temperature and salinity in the study area, evaluating the effects of the main forcing factors. The field results were depth-integrated and compared to numerical model results. These results obtained using extreme tidal and river runoff forcing, are also presented. The field results reveal that, when the river flow is weak, the tidal intrusion is the main forcing mechanism, generating saline and thermal fronts which migrate with the neap/spring tidal cycle. When the river flow increases, the influence of the freshwater extends almost as far as the mouth of the lagoon and vertical stratification is established. Results of numerical modelling reveal that the implemented model reproduces quite well the observed horizontal patterns. The model was also used to study the hydrology of the study area under extreme forcing conditions. When the model is forced with a low river flow (1 m3 s−1) the results confirm that the hydrology is tidally dominated. When the model is forced with a high river flow (1,000 m3 s−1) the hydrology is dominated by freshwater, as would be expected in such an area.  相似文献   

13.
This paper describes the results of numerical experiments with a three-dimensional hydrodynamic and mud transport model in which sediment–fluid interaction is taken into account through the effects of hindered settling, buoyancy destruction in the turbulence k– model and sediment-induced barocline pressure gradients in the momentum equations. The model was applied to a schematic case representing a coastal area with a tidal river, navigation channel and harbour basin, and a real-world case, viz. Rotterdam harbour area in The Netherlands. The results show that the sediment transport into the harbour area, and subsequent siltation rates, increase by a factor 3 to 5 due to the sediment–fluid interaction. It is shown that the larger contribution stems from an increase in vertical gradients in suspended sediment.Responsible Editor: Jens Kappenberg  相似文献   

14.
A geological and a geophysical survey have been carried out at Ponza Island, Tyrrhenian sea, Italy. Geological and structural data allowed to identify three main tectonic systems: NW-SE, NE-SW and E-W trending. The first one is related to a pre-volcanic tectonic event, probably linked to the Pliocene extensional activity of the Tyrrhenian evolution; the other two systems affected the volcanic units in two different stages of the Lower Pleistocene, the earlier one after the rhyolitic hyaloclastic formation (HF) emplacement and the later one after the emplacement of older trachytic pyroclastic deposits (Lower Pyroclastic Units—LPU). The latter event was followed by the emersion of the whole Ponza area, as testified by a marked erosional surface and marine terrace deposits cropping out at the top of LPU. The Upper Pyroclastic Units (UPU) represent the younger trachytic activity of the island (1.3 Myr) and do not show evidence of tectonic activity.The NW-SE-trending tectonic system probably assisted the rhyolitic magma rise, while the NE-SW- and E-W-trending systems mainly assisted the trachytic magma rise, responsible for the explosive and effusive activity in the southern area and for the hydrothermal fluids that caused alteration processes in the northern area.A 161-station gravimetric survey was carried out on the island and surrounding islets. The geological data and the gravimetric survey have been used to propose a 2.5 D model in which rhyolitic hyaloclastic deposits (ρ = 1.7 g cm−3) overlay an articulated Meso-Cenozoic sedimentary substratum (ρ = 2.6 g cm−3) laying at a depth to 300 m below sea level. Both formations are crossed by rhyolitic dykes (ρ = 2.4 g cm−3) which mark feeder fractures. In the M. Guardia area, where a maximum is present, this model accounts for the presence of a horst of the rigid basement, a shallow trachytic lava flow and its feeder (ρ = 2.8 g cm−3).  相似文献   

15.
The MECO hydrodynamic model (MECO Technical Report No. OMR-118/120, CSIRO Marine Research, 1998) was adapted for the Torres Strait–Gulf of Papua region at 0.05° resolution. Validation of the hydrodynamic model was carried out against observed current meter data and calculated tidal sea levels. Dispersal pathways of sediments derived from the Fly River, and from a resuspension event along the northern Great Barrier Reef were investigated using an Eulerian approach. Sediment input into Torres Strait is found to be greater during the Trade season by approximately 10%. Wave data were also obtained, and together with hydrodynamic model output, sediment mobility due to currents, waves and wave–current interactions was considered for both the Trade and Monsoon seasons. Sediment mobility in the Gulf of Papua is dominated by wave motion, whereas Torres Strait is a mixed environment of waves and tidal currents.  相似文献   

16.
Tidal water table fluctuations in a coastal aquifer are driven by tides on a moving boundary that varies with the beach slope. One-dimensional models based on the Boussinesq equation are often used to analyse tidal signals in coastal aquifers. The moving boundary condition hinders analytical solutions to even the linearised Boussinesq equation. This paper presents a new perturbation approach to the problem that maintains the simplicity of the linearised one-dimensional Boussinesq model. Our method involves transforming the Boussinesq equation to an ADE (advection–diffusion equation) with an oscillating velocity. The perturbation method is applied to the propagation of spring–neap tides (a bichromatic tidal system with the fundamental frequencies ω1andω2) in the aquifer. The results demonstrate analytically, for the first time, that the moving boundary induces interactions between the two primary tidal oscillations, generating a slowly damped water table fluctuation of frequency ω1−ω2, i.e., the spring–neap tidal water table fluctuation. The analytical predictions are found to be consistent with recently published field observations.  相似文献   

17.
The validation and subsequent application of the current three-dimensional numerical hydrodynamic model of Chesapeake Bay is presented. The numerical model solves conservation equations for water mass, momentum, salinity, and heat on a boundary-fitted grid in the horizontal plane and a Cartesian z-grid in the vertical. A generalized ADI finite difference scheme is employed in conjunction with mode splitting technique, solving external and the internal modes. The 10-year boundary conditions including tide, slinity, temperature, wind, heat exchange coefficient, river and non-point source flows were constructed. Model validation was accomplished by demonstrating the model's ability to reproduce observed data over time scales ranging from tidal to seasonal periods. The major parameters compared include tidal elevation, intra-tidal and residual velocities, salinity, temperature, stratification, and flux calculated through the Bay mouth.After validation, the model was applied to simulate bay hydrodynamics for the 10 years of 1985–94. These results were used to drive the three-dimensional water quality model of Chesapeake Bay, which is discussed in a companion paper.  相似文献   

18.
On the vertical structure of the Rhine region of freshwater influence   总被引：1，自引：0，他引：1  
An idealised three-dimensional numerical model of the Rhine region of fresh water influence (ROFI) was set up to explore the effect of stratification on the vertical structure of the tidal currents. Prandle’s dynamic Ekman layer model, in the case of zero-depth-averaged, cross-shore velocities, was first used to validate the response of the numerical model in the case of barotropic tidal flow. Prandle’s model predicted rectilinear tidal currents with an ellipse veering of up to 2%. The behaviour of the Rhine ROFI in response to both a neap and a spring tide was then investigated. For the given numerical specifications, the Rhine plume region was well mixed over the vertical on spring tide and stratified on neap tide. During spring conditions, rectilinear tidal surface currents were found along the Dutch coast. In contrast, during neap conditions, significant cross-shore currents and tidal straining were observed. Prandle’s model predicted ellipse veering of 50%, and was found to be a good indicator of ellipticity magnitude as a function of bulk vertical eddy viscosity. The modelled tidal ellipses showed that surface currents rotated anti-cyclonically whereas bottom currents rotated cyclonically. This caused a semi-diurnal cross-shore velocity shearing which was 90° out of phase with the alongshore currents. This cross-shore shear subsequently acted on the horizontal density gradient in the plume, thereby causing a semi-diurnal stratification pattern, with maximum stratification around high water. The same behaviour was exhibited in simulations of a complete spring–neap tidal cycle. This showed a pattern of recurring stratification on neaps and de-stratification on springs, in accordance with observations collected from field campaigns in the 1990’s. To understand the increase in ellipticities to 30% during neaps and the precise shape of the vertical ellipse structure, stratification has to be taken into account. Here, a full three-dimensional numerical model was employed, and was found to represent the effect of de-coupling of the upper and lower layers due to a reduction of mixing at the pycnocline.  相似文献   

19.
Profiles were analysed in conjunction with wave climate to assess offshore island influences on an embayed beach at Tenby, Wales. Time series analyses showed medium and short‐term beach oscillation, with volume exchanges between zones lagging by up to six months. Dominant southerly and southwesterly waves caused sub and low tidal longshore drift from south towards north, while less frequent southeasterly waves generated counter drift. Modelled inshore breaking waves had less energy than offshore ones and the former behaved differently between the low and high tidal zones (spring tidal range of 7 · 5 m). Variations in wave direction from directly behind the islands resulted in reduced wave heights and statistical analyses agreed with wave model results. These were correlated to morphological change and it was concluded that offshore islands change wave dynamics and modify the morphology of embayed beaches in their lee. Consequently, this work provides significant new insights into offshore island influences, shoreline behaviour and especially tidal setting Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

20.
While recent studies have revealed that tidal fluctuations in an estuary significantly affect groundwater flows and salt transport in the riparian zone, only seawater salinity in the estuary has been considered. A numerical study is conducted to investigate the influence of estuarine salinity variations on the groundwater flow and salt dynamics in the adjacent aquifer to extend our understanding of these complex and dynamic systems. Tidal salinity fluctuations (synchronous with estuary stage) were found to alter the magnitude and distribution of groundwater discharge to the estuary, which subsequently impacted on groundwater salinity patterns and residence times, especially in the riparian zone. The effects of salinity fluctuations were not fully captured by adopting a constant, time-averaged estuarine salinity. The modelling analysis also included an assessment of the impact of a seasonal freshwater flush in the estuary, similar to that expected in tropical climates (e.g. mean estuary level during flood significantly greater than average), on adjacent groundwater flow and salinity conditions. The three-month freshwater flushing event temporarily disrupted the salt distribution and re-circulation patterns predicted to occur under conditions of constant salinity and tidal water level fluctuations in the estuary. The results indicate that the salinity variations in tidal estuaries impact significantly on estuary–aquifer interaction and need to be accounted for to properly assess salinity and flow dynamics and groundwater residence times of riparian zones.  相似文献   

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

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