首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

2.
Groundwater flow and chemical transport in subterranean estuaries are poorly understood despite their potentially important implications for chemical fluxes from aquifers to coastal waters. Here, a numerical study of the dynamics in a subterranean estuary subject to tidal forcing is presented. Simulations show that salt transport associated with tidally driven seawater recirculation leads to the formation of an upper saline plume in the intertidal region. Computed transit times and flow velocities indicate that this plume represents a more active zone for mixing and reaction than the dispersion zone of the lower, classical salt wedge. Proper conceptualisation of this surficial mixing zone extends our understanding of processes within the subterranean estuary. Numerical tracer simulations reveal that tidal forcing may reduce the threat of a land-derived contaminant discharging to the marine environment by modifying the subsurface transport pathway and local geochemical conditions. Mixing and stratification in the subterranean estuary are strongly affected by both inland and tidal forcing. Based on the estuarine analogy we present a systematic classification of subterranean estuaries.  相似文献   

3.
The variation of seawater level resulting from tidal fluctuations is usually neglected in regional groundwater flow studies. Although the tidal oscillation is damped near the shoreline, there is a quasi‐steady‐state rise in the mean water‐table position, which may have an influence on regional groundwater flow. In this paper the effects of tidal fluctuations on groundwater hydraulics are investigated using a variably saturated numerical model that includes the effects of a realistic mild beach slope, seepage face and the unsaturated zone. In particular the impact of these factors on the velocity field in the aquifer is assessed. Simulations show that the tidal fluctuation has substantial consequences for the local velocity field in the vicinity of the exit face, which affects the nearshore migration of contaminant in coastal aquifers. An overheight in the water table as a result of the tidal fluctuation is observed and this has a significant effect on groundwater discharge to the sea when the landward boundary condition is a constant water level. The effect of beach slope is very significant and simplifying the problem by considering a vertical beach face causes serious errors in predicting the water‐table position and the groundwater flux. For media with a high effective capillary fringe, the moisture retained above the water table is important in determining the effects of the tidal fluctuations. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

4.
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.  相似文献   

5.
Sediment flux in marsh tidal creeks is commonly used to gauge sediment supply to marshes. We conducted a field investigation of temporal variability in sediment flux in tidal creeks in the accreting tidal marsh at China Camp State Park adjacent to northern San Francisco Bay. Suspended‐sediment concentration (SSC), velocity and depth were measured near the mouths of two tidal creeks during three 6‐ to 10‐week deployments: two in winter and one in summer. Currents, wave properties and SSC were measured in the adjacent shallows. All deployments spanned the largest spring tides of the season. Results show that tidally averaged suspended‐sediment flux (SSF) in the tidal creeks varied from slightly landward to strongly bayward with increasing tidal energy. SSF was negative (bayward) for tidal cycles with maximum water surface elevation above the marsh plain. Export during the largest spring tides dominated the cumulative SSF for each deployment. During ebb tides following the highest tides, velocities exceeded 1 m s?1 in the narrow tidal creeks, resulting in negative tidally averaged water flux, and mobilizing sediment from the creek banks or bed. Storm surge also produced negative SSF. Tidally averaged SSF was positive in wavy conditions with moderate tides. Spring tide sediment export at the creek mouth was about twice that at a station 130 m further up the tidal creek. The negative tidally averaged water flux near the creek mouth during spring tides indicates that in the lower marsh some of the water flooding directly across the bay–marsh interface drains through the tidal creeks, and suggests that this interface may be a pathway for sediment supply to the lower marsh as well. Copyright © 2018 John Wiley & Sons, Ltd.  相似文献   

6.
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.  相似文献   

7.
The outflow of brackish waters from the Tees estuary on an ebb tide forms a plume which is initially spread by the action of gravity. Under quiescent conditions, the rate of vertical mixing of the brackish field is relatively slow and an area of low salinity may form in the surface waters as a result of the superimposition of the outflow plumes from successive tides. A field area of some 70 km2 has been observed to be enclosed by the 34.2 × 10?3 contour.Averaged over a tidal period, the estuary outflow from the Tees forms a broad plume of low-salinity water which extends to seaward of the 20-m contour. The movement of the field thereafter depends on the residual coastal drift which runs either to the northwest or to the southeast with a typical speed of 0.015 m s?1. The tidally averaged vertical mixing coefficient in this offshore field was estimated to be about 10?4 m?2 s?1, which is similar to that off the Wear estuary and an order of magnitude smaller than that calculated for the waters off the Humber.  相似文献   

8.
A numerical study was conducted to investigate the influence of tides on the fate of terrestrially derived BTEX discharging through an unconfined aquifer to coastal waters. Previous studies have revealed that tide-induced seawater circulations create an active salt–freshwater mixing zone in the near-shore aquifer and alter the specific subsurface pathway for contaminants discharging to the coastal environment. Here the coupled density-dependent flow and multi-species reactive transport code PHWAT was used to examine the impact of these tidal effects on the aerobic biodegradation of BTEX released in a coastal aquifer and its subsequent loading to coastal waters. Simulations indicated that tides significantly enhance BTEX attenuation in the near-shore aquifer. They also reduce the rate of chemical transfer from the aquifer to the ocean and exit concentrations at the beach face. For the base case consisting of toluene transport and biodegradation, 79% of toluene initially released in the aquifer was attenuated prior to discharge with tides present, compared to only 1.8% for the non-tidal case. The magnitude of tidal forcing relative to the fresh groundwater flow rate was shown to influence significantly the extent of biodegradation as it controls the intensity of salt–freshwater mixing, period of exposure of the contaminant to the mixing zone and rate of oxygen delivery to the aquifer. The oxygen available for biodegradation also depends on the rate at which oxygen is consumed by natural processes such as organic matter decomposition. While simulations conducted with heterogeneous conductivity fields highlighted the uncertainties associated with predicting contaminant loadings, the study revealed overall that BTEX may undergo significant attenuation in tidally influenced aquifers prior to discharge.  相似文献   

9.
Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary’s entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429–35) for two-dimensional non-interacting tidal waves.  相似文献   

10.
Low-permeability layer (LPL), formed by natural deposit or artificial reclamation and commonly found below the intertidal zone of coastal groundwater system, can retard the ingress of seawater and contaminants, and shorten the travel time of the land-sourced contaminant to the marine environment compared with a homogenous sandy coastal aquifer. However, there is limited understanding on how an intertidal LPL, a condition occurred in a coastal aquifer at Moreton Bay, Australia, influences the groundwater and contaminant transport across the shallow beach aquifer system. We characterized the aquifer hydrological parameters, monitored the in situ groundwater heads, and constructed a 2-D numerical model to analyses the cross-shore hydrological processes in this stratified system. The calibrated model suggests that in the lower aquifer, the inland-source fresh groundwater flowed horizontally towards the sea, upwelled along the freshwater–saltwater interface, and exited the aquifer at the shore below the LPL. Whereas in the upper aquifer, the tidally driven seawater circulation formed a barrier that prevented fresh groundwater from horizontal transport and discharge to the beach above the LPL, thereby directing its leakage to the lower aquifer. A contaminant represented by a conservative tracer was ‘released’ the upper aquifer in the model and results showed that the spreading extent of the contaminant plume, the maximum rate of contaminant discharge to the ocean, and its plume length decreased compared with a simulation case in a homogenous sandy aquifer. Sensitivity analysis was also conducted to investigate the characteristics of the LPL, including its continuity and hydraulic conductivity, which were found to vary along the beach at Moreton Bay. The result shows that with a lower hydraulic conductivity and continuous layer of LPL reduced the groundwater exchange and contaminant transport between upper and lower aquifer. The findings from the combined field and modelling investigations on the impact of an intertidal LPL on coastal aquifer systems highlight its significant implications to alter the groundwater and mass transport across the land–ocean interface.  相似文献   

11.
Sediment deposition is one of the key mechanisms to counteract the impact of sea level rise in tidal freshwater wetlands (TFWs). However, information about sediment deposition rates in TFWs is limited, especially for those located in the transition zone between the fluvially dominated and tidally dominated sections of a river delta where sedimentation rates are affected by the combined impact of river discharge, wind, and tides. Using a combined hydrodynamic–morphological model, we examined how hydrometeorological boundary conditions control sedimentation rates and patterns in a TFW located in the Rhine–Meuse estuary in the Netherlands. The modelling results show that net sedimentation rate increases with the magnitude of the river discharge, whereas stronger wind increasingly prevents sedimentation. Sediment trapping efficiency decreases for both increasing river discharge and wind magnitude. The impact of wind storms on the trapping efficiency becomes smaller for higher water discharge. The spatial sedimentation patterns are affected by all controls. Our study illustrates the importance of evaluating both the separate and the joint impact of discharge, wind, and tides when estimating sedimentation rates in a TFW affected by these controls. Such insights are relevant to design measures to reactivate the sedimentation process in these areas.  相似文献   

12.
Monte Carlo simulations are conducted to evaluate microbial-mediated contaminant reactions in an aquifer comprised of spatially variable microbial biomass concentrations, aquifer hydraulic conductivities, and initial electron donor/acceptor concentrations. A finite element simulation model is used that incorporates advection, dispersion, and Monod kinetic expressions to describe biological processes. Comparisons between Monte Carlo simulations of heterogeneous systems and simulations using homogeneous formulation of the same two-dimensional transport problem are presented. For the assumed set of parameters, physical aquifer heterogeneity is found to have a minor effect on the mass of contaminant biodegraded/transformed when compared to a homogeneous system; however, it noticeably changes the dispersion, skewness, and peakness of contaminant concentration distributions. Similarly, for low microbial growth rate, given favorable microbial growth characteristics, biological heterogeneity has minor effect on the mass of contaminant biodegraded/transformed when compared to a homogeneous system. On the other hand, when higher effective growth rates are assumed, biological heterogeneity and spatial heterogeneities in essential electron donor/acceptors reduce the efficiency of biotic contaminant reactions; consequently, model simulations derived from heterogeneous biomass distributions predict remediation time scales that are longer than those simulated for homogeneous systems. When correlations between physical aquifer and biological heterogeneities are considered, the assumed correlation affects predicted mean and variance of contaminant concentration and biomass distributions. For example, an assumed negative correlation between hydraulic conductivity and the initial biomass distribution produces a plume where less efficient biotic contaminant reactions occur at the leading edge of the plume; this is consistent with less degradation/transformation occurring over regions of higher groundwater velocities. However, the presence and absence of these correlations do not appear to affect the efficiency of microbial-mediated contaminant attenuation.  相似文献   

13.
Observations of the flow field over an elongated hollow (bathymetric depression) in the lower Chesapeake Bay showed tidally asymmetric distributions. Current speed increased over the landward side of the hole during flood tides and decreased in the deepest part of the hollow during ebb tides. A simple conceptual analysis indicated that the presence of a horizontal density gradient can generate the asymmetric spatial variations of flow structure depending on the sign of the horizontal density gradient. When water density decreases downstream, the velocity increases over the downstream edge of the hollow. Conversely when water density increases downstream, the flow decreases over the hollow more than a case without a horizontal density gradient. The conceptual analysis is confirmed by numerical experiments of simplified hollows in steady open channel flows and of an idealized tidal estuary. These hollows also alter the local current field of tidally averaged estuarine exchange flows. The residual depth-averaged currents over a hollow show a two-cell circulation when Coriolis forcing is neglected and an asymmetric two-cell circulation, with a stronger cyclonic eddy, when Coriolis forcing is included.  相似文献   

14.
Biased monitoring of fresh water-salt water mixing zone in coastal aquifers   总被引：2，自引：0，他引：2  
In coastal aquifers, significant vertical hydraulic gradients are formed where fresh water and underlying salt water discharge together upward to the seafloor. Monitoring boreholes may act as "short circuits" along these vertical gradients, connecting between the higher and the lower hydraulic head zones. When a sea tide is introduced, the fluctuations of both the water table and the depth of the mixing zone are also biased due to this effect. This problem is intensified in places of long-screen monitoring boreholes, which are common in many places in the world. For example, all approximately 500 boreholes of the fresh water-salt water mixing zone in the coastal aquifer of Israel are installed with 10 to 50 m long screens. We present field measurements of these fluctuations, along with a three-dimensional numerical model. We find that the in-well fluctuation magnitude of the mixing zone is an order of magnitude larger than that in the porous media of the actual aquifer. The primary parameters that affect the magnitude of this bias are the anisotropy of the aquifer conductivity and the borehole hydraulic parameters. With no sea tide, borehole interference is higher for the anisotropic case because the vertical hydraulic gradients are high. When tides are introduced, the amplitude of the mixing zone fluctuation is higher for the isotropic case because the overall effective hydraulic conductivity is greater than the conductivity in the anisotropic case. In the aquifer, the fresh water-salt water mixing zone fluctuations are dampened, and tens of meters inland from the shoreline, the fluctuations are on the order of few centimeters.  相似文献   

15.
This paper addresses the impact of atmospheric variability on ocean circulation in tidal and non-tidal basins. The data are generated by an unstructured-grid numerical model resolving the dynamics in the coastal area, as well as in the straits connecting the North Sea and Baltic Sea. The model response to atmospheric forcing in different frequency intervals is quantified. The results demonstrate that the effects of the two mechanical drivers, tides and wind, are not additive, yet non-linear interactions play an important role. There is a tendency for tidally and wind-driven circulations to be coupled, in particular in the coastal areas and straits. High-frequency atmospheric variability tends to amplify the mean circulation and modify the exchange between the North and the Baltic Sea. The ocean response to different frequency ranges in the wind forcing is area-selective depending on specific local dynamics. The work done by wind on the oceanic circulation depends strongly upon whether the regional circulation is tidally or predominantly wind-driven. It has been demonstrated that the atmospheric variability affects the spring-neap variability very strongly.  相似文献   

16.
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.  相似文献   

17.
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.  相似文献   

18.
Analyses of independent laboratory- and field-scale measurements from two sites on Sapelo Island, Georgia reveal heterogeneity in hydraulic parameters across the upland–estuary interface. Regardless of the method used (short-duration pumping tests, amplitude attenuation of tidal pumping data, sediment grain size distributions, and falling head permeameter tests), we obtain hydraulic conductivity of 10−4 m s−1 for the fine-grained, well-sorted, clean sands that make up the upland areas. Proximal to the upland–estuary boundary, the tidal pumping analyses and permeameter tests suggest that hydraulic conductivities decrease by more than two orders of magnitude, a result consistent with the presence of a clogging layer. Such a clogging layer may arise due to a variety of physical, chemical, or biological processes. The extent and orientation of the layers of reduced hydraulic conductivity near the upland–estuary boundary influence the nature of the aquifer's response to tidal forcing. Where the lower conductivity layer forms a relatively flat creek bank, tidal pumping produces a primarily mechanical response in the adjacent aquifer. Where the creek bank is nearly vertical, there is a more direct hydraulic connection between the tidal creek and the adjacent aquifer. The clogging layer likely contributes to the development of complicated flow pathways across the upland–estuary boundary. Effective flow paths calculated from tidal pumping data terminate within the marsh, beyond the boundary of the upland aquifer, suggesting a diffuse regime of groundwater discharge in the marsh. We postulate that, in many settings, submarsh flow may be as important as seepage faces for groundwater discharge into the marsh–estuary complex.  相似文献   

19.
A three-dimensional primitive-equation model is used to simulate the Long Island Sound (LIS) outflow for a 1-year (2001) period. The model domain includes LIS and New York Bight (NYB). Tidal and wind forcing are included, and seasonal salinity and temperature variations are assimilated. The model results are validated with the HF radar, moored acoustic Doppler current profiler (ADCP), and ferry-based ADCP observations. The agreement between simulated and observed flow patterns generally is very good. The difference in seasonal mean currents between the model and moored ADCP is about 0.01 m/s; the correlation of dominant velocity fluctuations between the model and HF radar is 0.83; and the difference in mean LIS transport between the model and shipboard ADCP is about 5%. However, the model predicts a prominent tidally generated headland eddy not supported by the HF radar observation. The model sensitivity study indicates that the tides, winds, and ambient coastal front all have important impact on the buoyant outflow. The tides and winds cause stronger vertical mixing, which reduces the surface plume strength. The ambient coastal front, on the other hand, tends to enhance the plume.  相似文献   

20.
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  相似文献   

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

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