首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 546 毫秒
1.
Delivery of dissolved chemicals to bioremediate oil from the Exxon Valdez oil spill was investigated at Beach EL056C of Eleanor Island, Prince William Sound, Alaska. The delivery technique was high pressure injection (HPI) of an inert tracer, lithium, at the approximate depth of 1.0 m into the beach near the mid‐tide line. The results revealed that the maximum injection flow rate was 3.0 L/min and the associated pressure was around 196 kPa. Therefore, exceeding any of these values would probably cause failure of the injection system. The injected tracer was monitored at multiple depths of four surrounding observation wells, and the results showed that the tracer plume occupied an area of 12 m2 within 24 h. The tracer plume traveled at the average speeds of 10 m/d in the seaward direction and 1.7 m/d in other directions. The rapid movement under HPI and the large diameter of influence (3.0 m) indicated that bioremediation of the Exxon Valdez oil on this beach via injection of chemicals is logistically feasible.  相似文献   

2.
Measurements of velocity and density profiles were used to describe the tidal and mean flow structure across and along a sill in Refugio Channel, a fjord-like inlet in Southern Chile (43.9°S). These are the first oceanographic measurements of any kind effected in Refugio Channel. Current profiles were obtained with a 307.2-kHz acoustic Doppler current profiler during two semidiurnal cycles along a repeated triangular circuit. Two along-channel transects formed the sides of the triangle that crossed the sill and were identified as the western and eastern transects. One cross-channel transect, the base of the triangle, was located on the seaward side of the sill. Density profiles were obtained at the corners of the triangle. The longitudinal mean flow in the western transect showed a two-layer exchange structure over the landward side of the sill. The structure of net seaward flow at the surface and landward flow at depth was disrupted by the sill in such a way that over the seaward side of the sill, only seaward flow was observed throughout the water column. This likely resulted from the blocking of landward net flow by the sill. In the eastern transect, two-layer exchange dominated over most of the transect and was consistent with the observed density profiles. Over the seaward side of the sill, a surface layer, ∼10m deep, flowed landward as a third layer. This feature should have been caused by river input further seaward (to the north) and produced a surface convergence region over the sill. In terms of tidal flows, the greatest tidal current amplitudes were 40cm s−1 over the sill as the flow accelerated through the reduced cross-sectional area of the channel. Near-surface flow convergences were identified over both along-channel transects.  相似文献   

3.
The geomorphology of the southern Yellow Sea（SYS） is characterized by offshore radial sand ridges（RSR）.An offshore tidal channel（KSY Channel） is located perpendicular to the coast,comprised of a main and a tributary channel separated by a submarine sand ridge（KSY Sand Ridge） extending seaward.In order to investigate the interactions among water flow,sediment transport,and topography,current velocity and suspended sediment concentration（SSC） were observed at 11 anchor stations along KSY Channel in RSR during a spring tide cycle.High resolution bottom topography was also surveyed.Residual currents and tidally averaged suspended sediment fluxes were calculated and analyzed by using the decomposition method.Results suggested that the water currents became stronger landward but with asymmetrical current speed and temporal duration of flood and ebb tides.Residual currents showed landward water transport in the nearshore channel and a clockwise circulation around the KSY Sand Ridge.Tidally-averaged SSC also increased landward along the channel.The main mechanisms controlling SSC variations were resuspension and horizontal advection,with spatial and temporal variations in the channel,which also contributed to sediment redistribution between channels and sand ridges.Residual flow transport and the tidal pumping effect dominated the suspended sediment flux in the KSY Channel.The KSY Sand Ridge had a potential southward migration due to the interaction between water flow,sediment transport,and topography.  相似文献   

4.
This paper presents results of investigations (1983–1992) into rates of change, morphology and processes occurring during the current erosional phase in a Morecambe Bay cyclic saltmarsh, in which it has narrowed from c. 1000 m (1975) to c. 150 m (1992). Monthly monitoring of marsh edge erosion and creek knickpoint retreat has revealed temporal and spatial variations. Highest frequency changes of low magnitude coincided with non-storm conditions and overmarsh tides above 5·80 m OD, which submerged the whole marsh. Less frequent changes of greater magnitude were associated with both overmarsh tides and strong onshore winds over 15 ms?1, which generated high energy waves. The lowest frequency change of greatest magnitude occurred during an extreme onshore storm event and surge. Morphologically, during the erosional phase, a low angled landward slope was generated as erosion of the c. 0·5 m high active seaward cliff coincided with vertical accretion of 0·07 ma?1 of relatively coarse sediment on the marsh surface immediately landward. Tidal hydrodynamics strongly influence the saltmarsh, which is confined to the upper 2·5 m of the macrotidal range (maximum c. 10·5 m). During overmarsh spring tides (maximum creek flood flow rate 0·13 ms?1, up to bankfull level), flooding begins over lower landward creek banks before submerging the higher marsh edge. During ebb tides, water trapped by this higher edge can escape seaward only via the creeks (maximum ebb velocities 2·07 ms?1 below bankfull level). Wave erosion also is limited to spring tides. Monthly mapping of the Kent Estuary channel pattern seaward of the saltmarsh showed that medium term higher erosion rates were related to the presence of a large channel, which lowered the adjacent creek base level and allowed larger waves to attack the marsh edge than when a sandbank flanked the marsh. Major River Kent channel shifts appear to initiate accretional or erosional phases of cyclic saltmarsh development.  相似文献   

5.
A comprehensive numerical study was undertaken to investigate transport of a variable-density, conservative solute plume in an unconfined coastal aquifer subject to high and low frequency oceanic forcing. The model combined variable-density saturated flow for groundwater and solute transport, and wave hydrodynamics from a 2D Navier–Stokes solver. A sinusoidal tidal signal was specified by implementing time-varying heads at the seaward boundary. The solute plume behavior was investigated under different oceanic forcing conditions: no forcing, waves, tide, and combined waves and tide. For each forcing condition, four different injected solute densities (freshwater, brackish water, seawater, brine) were used to investigate the effects of density on the transport of the injected plume beneath and across the beach face. The plume’s low-order spatial moments were computed, viz., mass, centroid, variance and aspect ratio. The results confirmed that both tide- and wave-forcing produce an upper saline plume beneath the beach face in addition to the classical saltwater wedge. For the no-forcing and tide-only cases (during rising tides), an additional small circulation cell below the beach face was observed. Oceanic forcing affects strongly the solute plume’s flow path, residence time and discharge rate across the beach face, as well as its spreading. For the same oceanic forcing, solute plumes with different densities follow different trajectories from the source to the discharge location (beach face). The residence time and plume spreading increased with plume density. It was concluded that simulations that neglect the effect of waves or tides cannot reproduce accurately solute plume dispersion and also, in the case of coasts with small waves or tides, the solute residence time in the aquifer.  相似文献   

6.
Bifurcations in tidally influenced deltas distribute river discharge over downstream channels, asserting a strong control over terrestrial runoff to the coastal ocean. Whereas the mechanics of river bifurcations is well-understood, junctions in tidal channels have received comparatively little attention in the literature. This paper aims to quantify the tidal impact on subtidal discharge distribution at the bifurcations in the Mahakam Delta, East Kalimantan, Indonesia. The Mahakam Delta is a regular fan-shaped delta, composed of a quasi-symmetric network of rectilinear distributaries and sinuous tidal channels. A depth-averaged version of the unstructured-mesh, finite-element model second-generation Louvain-la-Neuve Ice-ocean Model has been used to simulate the hydrodynamics driven by river discharge and tides in the delta channel network. The model was forced with tides at open sea boundaries and with measured and modeled river discharge at upstream locations. Calibration was performed with water level time series and flow measurements, both spanning a simulation period. Validation was performed by comparing the model results with discharge measurements at the two principal bifurcations in the delta. Results indicate that within 10 to 15 km from the delta apex, the tides alter the river discharge division by about 10% in all bifurcations. The tidal impact increases seaward, with a maximum value of the order of 30%. In general, the effect of tides is to hamper the discharge division that would occur in the case without tides.  相似文献   

7.
A natural gradient tracer test using perdeuterated MTBE was conducted in an anaerobic aquifer to determine the relative importance of dispersion and degradation in reducing MTBE concentrations in ground water. Preliminary ground water chemistry and hydraulic conductivity data were used to place the tracer within an existing dissolved MTBE plume at Port Hueneme, California. Following one year of transport, the tracer plume was characterized in detail.
Longitudinal dispersion was identified as the dominant mechanism for lowering the perdeuterated MTBE concentrations. The method of moments was used to determine the longitudinal and lateral dispersion coefficients (0.85 m2/day and 0.08 m2/day, respectively). A mass-balance analysis, carried out after one year of transport, accounted for 110% of the injected mass and indicated that no significant mass loss occurred. The plume structure created by zones of higher and lower hydraulic conductivity at the site was complex, consisting of several localized areas of high tracer concentration in a lower concentration plume. This is important because the aquifer has generally been characterized as exhibiting fairly minor heterogeneity. In addition, the tracer plume followed a curved flowpath that deviated from the more macroscopic direction of ground water flow inferred from local ground water elevation measurements and the behavior of the existing plume. Understanding the mass balance, plume structure, curvature of the tracer plume, and consequently natural attenuation behavior required the detailed sampling approach employed in this study. These data imply that a detailed understanding of site hydrogeology and an extensive sampling network may be critical for the correct interpretation of monitored natural attenuation of MTBE.  相似文献   

8.
The increase of salt intrusion in recent years in the Modaomen Estuary, one of the estuaries of the Pearl River Delta in China, has threatened the freshwater supply in the surrounding regions, especially the cities of Zhongshan, Zhuhai in Guangdong Province and Macau. A numerical modeling system using nested grids was developed to investigate the salt transport mechanisms and the response of salt intrusion to changes in river discharge and tidal mixing. The steady shear transport induced by estuarine circulation reaches maximum and minimum, respectively, during neap and spring tides, while the tidal oscillatory transport shows an opposite pattern. The net transport is landward during neap tides and seaward during spring tides. The salt intrusion length responding to constant river discharges generally follows a power law of ?0.49. The dependence of salt intrusion on tidal velocity is less than that predicted by theoretical models for exchange flow dominated estuaries. The response of salt intrusion to change in tidal velocity depends largely on river discharge. When river flow increases, the impact of tidal velocity increases and the phase lag of response time decreases. The asymmetries of salt intrusion responding to increasing and decreasing river discharge (tidal velocity) are observed in the estuary.  相似文献   

9.
Resource extraction and transportation activities in subarctic Canada can result in the unintentional release of contaminants into the surrounding peatlands. In the event of a release, a thorough understanding of solute transport within the saturated zone is necessary to predict plume fate and the potential impacts on peatland ecosystems. To better characterize contaminant transport in these systems, approximately 13,000 L/day of sodium chloride tracer (200 mg/L) was released into a bog in the James Bay Lowland. The tracer was pumped into a fully penetrating well (1.5 m) between July 5 and August 18, 2015. Horizontal and vertical plume development was measured via in situ specific conductance and water table depth from an adaptive monitoring network. Over the spill period, the bulk of the plume travelled a lateral distance of 100 m in the direction of the slight regional groundwater and topographical slope. The plume shape was irregular and followed the hollows, indicating preferential flow paths due to the site microtopography. Saturated transport of the tracer occurred primarily at ~25 cm below ground surface (bgs), and at a discontinuous high hydraulic conductivity layer ~125 cm bgs due to a complex and heterogeneous vertical hydraulic conductivity profile. Plume measurement was confounded by a large amount of precipitation (233 mm over the study period) that temporarily diluted the tracer in the highly conductive upper peat layer. Longitudinal solute advection can be approximated using local water table information (i.e., depth and gradient); microtopography; and meteorological conditions. Vertical distribution of solute within the peat profile is far more complex due to the heterogeneous subsurface; characterization would be aided by a detailed understanding of the site‐specific peat profile; the degree of decomposition; and the type of contaminant (e.g., reactive/nonreactive). The results of this research highlight the difficulty of tracking a contaminant spill in bogs and provide a benchmark for the characterization of the short‐term fate of a plume in these complex systems.  相似文献   

10.
Halifax Harbour is located on the Atlantic coast of Nova Scotia, Canada. It is one of the world’s largest, ice-free natural harbours and of great economic importance to the region. A good understanding of the physical processes controlling tides, flooding, transport and dispersion, and hydrographic variability is required for pollution control and sustainable development of the Harbour. For the first time, a multi-nested, finite difference coastal ocean circulation model is used to reconstruct the three-dimensional circulation and hydrography of the Harbour and its variability on timescales of hours to months for 2006. The model is driven by tides, wind and sea level pressure, air-sea fluxes of heat, and terrestrial buoyancy fluxes associated with river and sewage discharge. The predictive skill of the model is assessed by comparing the model simulations with independent observations of sea level from coastal tide gauges and currents from moored instruments. The simulated hydrography is also compared against a new monthly climatology created from all available temperature and salinity observations made in the Harbour over the last century. It is shown that the model can reproduce accurately the main features of the observed tides and storm surge, seasonal mean circulation and hydrography, and wind driven variations. The model is next used to examine the main physical processes controlling the circulation and hydrography of the Harbour. It is shown that non-linear interaction between tidal currents and complex topography occurs over the Narrows. The overall circulation can be characterized as a two-layer estuarine circulation with seaward flow in the thin upper layer and landward flow in the broad lower layer. An important component of this estuarine circulation is a relatively strong, vertically sheared jet situated over a narrow sill connecting the inner Harbour to the deep and relatively quiescent Bedford Basin. Local wind driven variability is strongest in winter as expected but it is also shown that a significant part of the temperature and salinity variability is driven by physical processes occurring on the adjacent inner continental shelf, especially during storm and coastal upwelling events.  相似文献   

11.
The main regularities in the interaction of tides and storm surges at river mouths are discussed. A study of the Elbe River mouth area is used to describe the processes of interaction of the eustatic sea level rise, tides, surges, and river flow and special features of formation of maximum water levels. As shown, the intensification of cyclonic activity over the Northern Atlantic in the second half of the XX century resulted in more frequent extremely high storm surges at the Elbe River mouth. An assessment is given for possible changes in the regime of tides and surges at the Elbe River mouth in the XXI century, which may be caused by the acceleration of the eustatic sea level rise. The impact of local hydraulic engineering works (diking, dredging, and channel straightening) on maximum water levels within the town of Hamburg is analyzed.  相似文献   

12.
Field observations of instantaneous water surface slopes in the swash zone are presented. For free-surface flows with a hydrostatic pressure distribution the surface slope is equivalent to the horizontal pressure gradient. Observations were made using a novel technique which in its simplest form consists of a horizontal stringline extending seaward from the beach face. Visual observation, still photography or video photography is then sufficient to determine the surface slope where the free-surface cuts the line or between reference points in the image. The method resolves the mean surface gradient over a cross-shore distance of 5 m or more to within ±0.001, or 1/20th −1/100th of typical beach gradients. In addition, at selected points and at any instant in time during the swash cycle, the water surface slope can be determined exactly to be dipping either seaward or landward. Close to the location of bore collapse landward dipping water surface slopes of order 0.05–0.1 occur over a very small region (order 0.5 m) at the blunt or convex leading edge of the swash. In the middle and upper swash the water surface slope at this leading edge is usually very close to horizontal or slightly seaward. Behind the leading edge, the water surface slope was observed to be very close to horizontal or dipping seaward at all times throughout the swash uprush. During the backwash the water surface slope was observed to be always dipping seaward, approaching the beach slope, and remained seaward until a new uprush edge or incident bore passed any particular cross-shore location of interest. The observations strongly suggest that the swash boundary layer is subject to an adverse pressure gradient during uprush and a favourable pressure gradient during the backwash. Furthermore, assuming Euler's equations are a good approximation in the swash, the observations also show that the total fluid acceleration is negative (offshore) for almost the whole of the uprush and for the entire backwash. The observations are contrary to recent work suggesting significant shoreward directed accelerations and pressure gradients occur in the swash (i.e., ∂u/∂t>0∼∂p/∂x<0), but consistent with analytical and numerical solutions for swash uprush and backwash. The results have important implications for sediment transport modelling in the swash zone.  相似文献   

13.
The flow caused by the discharge of freshwater underneath a glacier into an idealized fjord is simulated with a 2D non-hydrostatic model. As the freshwater leaves horizontally the subglacial opening into a fjord of uniformly denser water it spreads along the bottom as a jet, until buoyancy forces it to rise. During the initial rising phase, the plume meanders into complex flow patterns while mixing with the surrounding fluid until it reaches the surface and then spreads horizontally as a surface seaward flowing plume of brackish water. The process induces an estuarine-like circulation. Once steady-state is reached, the flow consists of an almost undiluted buoyant plume rising straight along the face of the glacier that turns into a horizontal surface layer thickening as it flows seaward. Over the range of parameters examined, the estuarine circulation is dynamically unstable with gradient Richardson number at the sheared interface having values of <1/4. The surface velocity and dilution factors are strongly and non-linearly related to the Froude number. It is the buoyancy flux that primarily controls the resulting circulation with the momentum flux playing a secondary role.  相似文献   

14.
As an initial evaluation of the potential of digital elevation models (DEMs) and geographic information systems (GISs) for geomorphic characterization of rocky shorelines, airborne laser scan (ALS) data have been used to characterize shore platforms around Shag Point, southeastern New Zealand. The platforms have been characterized using field‐based techniques in previously published research, and therefore offer an ideal site for evaluation purposes. The main challenge involved the delineation of the shore platform area in terms of landward and seaward extents. The cliff top and landward edge of the shore platform was readily mapped, whereas the seaward edge of platforms was mapped with lesser precision due to difficulties associated with tidal inundation and the interference of wave action and surface water. In the central region of the study area (～0·1 km2) higher platform elevations and dense point cloud data enabled the generation of a high‐resolution (1 m) DEM. In analysing the DEM, ALS offered an advantage over the previous field survey in respect of the ability to assess continuous topography in plan‐view. The extent and form of two distinctive erosional surfaces is clearly apparent and was revealed through classifications based on slope and elevation. The spatial continuity of the upper surface implies that, in addition to the role of rock structure described in previous work, sea level and wave exposure may have been important factors in the generation and preservation of platform morphology at Shag Point. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

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

16.
A three-dimensional, hydrostatic, primitive equation numerical model with modern turbulence closures is used to explore lateral circulation and the associated transport of sediments in idealized, moderately to highly stratified estuaries. The model results suggest that boundary mixing on a sloping bottom can drive a significant amount of lateral circulation. This mechanism has received little attention to date in the estuarine literature. Good agreement with an analytical solution and similar vertical structures of lateral flows to observations from the Hudson River estuary support the importance of the boundary mixing mechanism. Boundary mixing is at least as important as differential advection for the modeled scenarios, when the two mechanisms are evaluated using the salt balance equation for model runs without rotation. Linearly superposing analytical solutions for lagged boundary mixing lateral flow and Ekman-forced lateral flow yields a good representation of the near-bottom lateral flow from the model with rotation. The 2 h lag required for the boundary mixing solution is roughly equal to the vertical diffusion time scale, indicating that lateral flow adjustment depends on development of a bottom mixed layer. Sediment dynamics at cross sections seaward and landward of the salt intrusion are very different. Seaward of the salt intrusion, sediments are eroded in the channel and preferentially deposited on the right slope (looking seaward), mainly due to the combination of high sediment concentration in the channel during flood with strong up-slope transport on that side (tidal pumping). Lateral sediment re-distribution landward of the salt intrusion is negligible due to weak residual lateral circulation.  相似文献   

17.
The broad (～500 km) southeastern Bering Sea continental shelf contains three fronts; outer (shelf break, ～170-m depth), middle (～100-m depth), and inner (～50-m depth). The shelf break and inner fronts appear to be analogous to similar fronts reported from other mid-latitude continental shelves; extensively studied examples are from the mid-Atlantic bight, off Nova Scotia, and around the British Isles. The middle front may have counterparts on the broad North Sea and East China Sea shelves.One-month current and temperature records from either side of the middle front, ～150 km landward from the shelf break, showed convergence in the layers deeper than 30 m in both the cross-shelf flow field and heat flux. The convergence was ～3 cm s?1, so an average upwelling at ～1 × 10?3cm s?1 and divergence in the surface layer were required to maintain continuity. Variations in the degree of convergence arose primarily from 1 to 5-day fluctuations in sub-tidal flow across the outer shelf domain seaward of the front.Diffusive landward heat flux was dominated by tidal scales. Horizontal eddy conductivities describing the flux were ～1 ? 106 on the landward side and ～5 × 106cm2 s?1 on the seaward side, and were less in the layers above the bottom layer. Advective flux by the mean flow was the same order as diffusive flux, but landward in the bottom layer and seaward in the mid-water column layers, in agreement with deductions from water mass analyses. Frontal effects reduced the net cross-shelf heat flux beneath ～30 m by about 50%. The observation of a flow convergence in the middle of a broad, flat continental shelf poses an important question of dynamics.  相似文献   

18.
The Adriatic coast of Punta Marina (Ravenna) is protected by 3-km long low crested breakwater structures (LCSs). Through a 3-years long multidisciplinar study, we assessed the impact of such defensive structures on environmental and biological condition. LCSs create pools where conditions are very different from the surrounding nearshore system. Mechanical disturbance by currents and waves varied greatly in intensity and frequency between seaward and landward sides of the structures. Sedimentary budget was positive at the landward side, but it was due to a gain on the seafloor and not on the emerged beach. The budget at seaward was negative. LCSs determine differences in benthic assemblages, alter the seasonal pattern of communities, and modify seasonal fluctuations of animal assemblages. Landward sheltered areas can be seen as “lagoonal island” surrounded by a “sea of marine habitat”. Differences in ecological quality status, obtained through M-AMBI, are due to the sum of these factors.  相似文献   

19.
The most important natural function of exposed sandy beaches is the filtration of large volumes of sea water which are flushed through the interstitial pore system by the actions of tides and waves. In this process organic material is mineralized by the interstitial fauna and nutrients returned to the sea. This paper examines the physical effects of oil in blocking or reducing this filtration. Reduced filtration of sea water may be caused by prevention or retardation of interstitial water flow and/or decrease in pore space. The magnitude of the effects depends on the volume of oil, its state of weathering, its location on the beach and its degree of admixture with sand. It appears that under all but the heaviest conditions of pollution these effects are relatively small and only become significant in situ when discrete layers of weathered oil occur near the surface in the upper intertidal zone where maximum water input normally occurs.  相似文献   

20.
The problem of seawater intrusion in a confined coastal aquifer is investigated. The aquifer is overlain by a leaky aquitard and both units extend continuously out under the sea. The problem is formulated in terms of the two governing differential equations, the fluid continuity equation written conveniently in terms of equivalent freshwater head, and the mass continuity equation. Use of linear rectangular finite elements, with direct integration and an iterative solution technique lead to an efficient numerical scheme that is capable of handling long simulation periods. The results, for a 300 m thick aquifer overlain by a 100 m thick aquitard, show that the aquitard has a controlling influence on the salt distribution. A zone of mixing in the aquifer is found to extend for several kilometres in the seaward as well as the landward direction. The time period required by the system to attain a state of dynamic equilibrium after a perturbation is applied may be of the order of centuries. The aquitard is found to cause a downward and seaward displacement of the average salt front.  相似文献   

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

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