The Application of Combined Seismic and Electrical Measurements to the Determination of the Hydraulic Conductivity of a Lake Bed

The hydraulic properties of lake beds control the interactions between lakes and ground water systems, but these properties are normally difficult to measure directly. The authors'method combines seismic reflection and electrical measurements to map the relative hydraulic conductivity of lake bed sediments. A shipboard seismic profiling system provides sediment thickness, while a towed electrical array yields longitudinal conductance and electrical chargeability. The sediment's leakance (hydraulic conductivity/thickness) can be calculated from the longitudinal conductance data. Leakance may then be converted to relative hydraulic conductivity through the seismically derived sediment thicknesses. Simultaneously acquired electrical chargeability provides an independent measure of clay content. The seismic and electrical systems are computer automated and yield production rates of approximately five line-kilometers/hour or 300 electrical soundings/hour. The systems provide continuous hydraulic information along the ship track rather than the point information derived from coring.
The procedure and systems have been used to map the bed of Lake Michigan offshore from an area of heavy pumpage. This location has been chosen to test the method because lake water has intruded the aquifer in plumes largely controlled by lake bed hydraulics. Mapping these plumes onshore permits the inference of the spatial distribution of offshore hydraulic conductivities. Offshore seepage measurements and numerical, chemical transport modeling of this site have confirmed the reliability of the geophysically derived hydraulic conductivities and have also demonstrated the improvement in numerical results achieved through the availability of spatially determined hydraulic conductivities.     

River Seepage Conductance in Large‐Scale Regional Studies

Flow exchange between surface and groundwater is of great importance be it for beneficial allocation and use of water resources or for the proper exercise of water rights. In large‐scale regional studies, most numerical models use coarse grid sizes, which make it difficult to provide an accurate depiction of the phenomenon. In particular, a somewhat arbitrary leakance coefficient in a third type (i.e., Cauchy, General Head) boundary condition is used to calculate the seepage discharge as a function of the difference of head in the river and in the aquifer, whose value is often found by calibration. A different approach is presented to analytically estimate that leakance coefficient. It is shown that a simple equivalence can be deduced from the analytical solution for the empirical coefficient, so that it provides the accuracy of the analytical solution while the model maintains a very coarse grid, treating the water‐table aquifer as a single calculation layer. Relating the empirical leakance coefficient to the exact conductance, derived from physical principles, provides a physical basis for the leakance coefficient. Factors such as normalized wetted perimeter, degree of penetration of the river, presence of a clogging layer, and anisotropy can be included with little computational demand. In addition the river coefficient in models such as MODFLOW, for example, can be easily modified when grid size is changed without need for recalibration.     

Transient Leakance and Infiltration Characteristics during Lake Bank Filtration

Infiltration capacity of bank filtration systems depends on water extraction and hydraulic resistance of the bed sediments. Lakebed hydraulics may be especially affected by clogging, which is dependent on settlement of fine particles, redox potential, and other factors. In the field, most of these processes are difficult to quantify, and thus, when calculating response to pumping the water flux across the sediment surface is assumed to be linearly dependent on the hydraulic gradient. However, this assumption was not adequate to describe conditions at a bank filtration site located at Lake Tegel, Berlin, Germany. Hence, we first assumed the leakage coefficient (or leakance) is spatially distributed and also temporally variant. Furthermore, observations show that the leakance is considerably higher in shallow than in deeper areas; hence, leakance was assumed to be dependent on the existence and thickness of an unsaturated zone below the lake. The proposed explanation of spatial and temporal variability in leakance involves a hypothesis for redox dependent and reversible biogeochemical clogging, supported by geochemical observations in surface water and ground water. Four leakance approaches are implemented in the ground water flow code MODFLOW2000 and calibrated by inverse modeling using the parameter estimation software PEST. These concepts are evaluated by examining the fit to the hydraulic heads, to infiltration measurements, transport modeling results, and considering the degrees of freedom due to the number of calibration parameters. The leakage concept based on the assumption of the influence of an unsaturated zone on clogging processes best explains the field data.     

Assessing Lakebed Hydraulic Conductivity and Seepage Flux by Potentiomanometer

Potentiomanometers (PMs) are commonly used to determine flux directions across interfaces between surface waters and aquifers. We describe a complementary function: estimating small‐scale hydraulic conductivity (K) in a lakebed, using the constant‐head injection test (CHIT) by Cardenas and Zlotnik (2003) with the PM designed by Winter et al. (1988). A piezometer with a small screen is inserted into the lakebed. Local head potential is obtained by measuring the head difference between the test point and the aquifer interface. The piezometer is then used for water injection. This technique is illustrated by measurements taken from Alkali Lake in the Sand Hills, Nebraska, United States. Lakebed K and seepage fluxes ranged from 0.037 to 0.090 m/d and Darcy velocities ranged from 0.004 to 0.027 m/d. Results were consistent with the supplementary data gathered using a modified CHIT and a cone penetrometer. The compact size of the device and the small volumes used for injection enable this method to estimate lakebed K values as low as 0.01 to 0.1 m/d, a range seldom explored in lake‐aquifer interface systems.     

Tide‐induced groundwater head fluctuation in a coastal aquifer system with a submarine outcrop covered by a thin silt layer

We present an analytical solution of groundwater head response to tidal fluctuation in a coastal multilayered aquifer system consisting of an unconfined aquifer, a leaky confined aquifer and a semi‐permeable layer between them. The submarine outcrop of the confined aquifer is covered by a thin silt layer. A mathematical model and the analytical solution of this model are given. The silt layer reduces the amplitude of the hydraulic head fluctuation by a constant factor, and shifts the phase by a positive constant (time lag), both of which depend on the leakances of the silt layer and the semi‐permeable layer. The time lag is less than 1·5 h and 3·0 h for semi‐diurnal and diurnal sea tides respectively. When the leakance of the semi‐permeable layer or the silt layer assumes certain special values, the solution becomes the existing solutions derived by previous researchers. The amplitude of the hydraulic head fluctuation in the confined aquifer increases with the leakance of the silt layer and decreases with the leakance of the semi‐permeable layer, whereas the phase shift of the fluctuation decreases with both of them. A hypothetical example shows that neglecting the silt layer may result in significant parameter estimation discrepancy between the amplitude attenuation and the time‐lag fittings. Copyright © 2007 John Wiley & Sons, Ltd.     

Use of an electromagnetic seepage meter to investigate temporal variability in lake seepage

A commercially available electromagnetic flowmeter is attached to a seepage cylinder to create an electromagnetic seepage meter (ESM) for automating measurement of fluxes across the sediment/water interface between ground water and surface water. The ESM is evaluated through its application at two lakes in New England, one where water seeps into the lake and one where water seeps out of the lake. The electromagnetic flowmeter replaces the seepage-meter bag and provides a continuous series of measurements from which temporal seepage processes can be investigated. It provides flow measurements over a range of three orders of magnitude, and contains no protruding components or moving parts. The ESM was used to evaluate duration of seepage disturbance following meter installation and indicated natural seepage rates resumed approximately one hour following meter insertion in a sandy lakebed. Lakebed seepage also varied considerably in response to lakebed disturbances, near-shore waves, and rainfalls, indicating hydrologic processes are occurring in shallow lakebed settings at time scales that have largely gone unobserved.     

Simulation of Lake-Groundwater Interaction under Steady-State Flow

MODFLOW is one of the most popular groundwater simulation tools available; however, the development of lake modules that can be coupled with MODFLOW is lacking apart from the LAK3 package. This study proposes a new approach for simulating lake - groundwater interaction under steady-state flow, referred to as the sloping lakebed method (SLM). In this new approach, discretization of the lakebed in the vertical direction is independent of the spatial discretization of the aquifer system, which can potentially solve the problem that the lake and groundwater are usually simulated at different scales. The lakebed is generalized by a slant at the bottom of each lake grid cell, which can be classified as fully submerged, dry, and partly submerged. The SLM method accounts for all lake sources and sinks, establishing a governing equation that can be solved using Newton's method. A benchmarking case study was conducted using a modified model setup in the LAK3 user manual. It was found that when there is a sufficient number of layers at the top of the groundwater model, SLM simulates an almost identical groundwater head as the LAK3-based model; when the number of layers decreases, SLM is unaffected while LAK3 may be at a risk of giving unrealistic results. Additionally, the SLM can reflect the relationship between the simulated lake surface area and lake water depth more accurately. Therefore, the SLM method is a promising alternative to the LAK3 package when simulating lake - groundwater interaction.     

Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods

Abstract

The resolution of the freshwater and saline water aquifers in a coastal terrain (Mahanadi Basin, India) is updated. We analysed electrical borehole log data at four sites and compared the water resistivity regime of the freshwater and saline water zones obtained from electrical borehole logging, with the resistivity regime obtained by interpreting vertical electrical sounding (VES) data. The multilayer VES data interpretation is modified to a simple model, containing only the freshwater zone and the saline water zone. The composite geophysical parameters of the freshwater and saline water zones, in particular the resistivity and longitudinal unit conductance regime, are identified. The resolution obtained from the composite geophysical data analyses is very clear and convincing. The composite longitudinal unit conductance regime of the saline water zones is very high compared to that of the freshwater zones. This makes the identification of the two aquifers easy and increases its reliability. A technique which enables analysis of composite geophysical data of freshwater and saline water zones at VES sites in the vicinity of the borehole log sites is proposed. The significance of longitudinal unit conductance in resolving the freshwater and saline water aquifers is illustrated graphically. The proposed technique is validated by correlating the longitudinal unit conductance and resistivity with the total dissolved solids. The efficiency of the technique is validated by carrying out discriminant function analysis.

Citation Hodlur, G. K., Dhakate, R., Sirisha, T. & Panaskar, D. B. (2010) Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods. Hydrol. Sci. J. 55(3), 414–434.     

Simulation of multi-support depth-varying earthquake ground motions within heterogeneous onshore and offshore sites

This paper presents a novel approach to model and simulate the multi-support depth-varying seismic motions (MDSMs) within heterogeneous offshore and onshore sites. Based on 1D wave propagation theory, the three-dimensional ground motion transfer functions on the surface or within an offshore or onshore site are derived by considering the effects of seawater and porous soils on the propagation of seismic P waves. Moreover, the depth-varying and spatial variation properties of seismic ground motions are considered in the ground motion simulation. Using the obtained transfer functions at any locations within a site, the offshore or onshore depth-varying seismic motions are stochastically simulated based on the spectral representation method (SRM). The traditional approaches for simulating spatially varying ground motions are improved and extended to generate MDSMs within multiple offshore and onshore sites. The simulation results show that the PSD functions and coherency losses of the generated MDSMs are compatible with respective target values, which fully validates the effectiveness of the proposed simulation method. The synthesized MDSMs can provide strong support for the precise seismic response prediction and performance-based design of both offshore and onshore large-span engineering structures.     

Drawdown and stream depletion produced by pumping in the vicinity of a partially penetrating stream

Commonly used analytical approaches for estimation of pumping-induced drawdown and stream depletion are based on a series of idealistic assumptions about the stream-aquifer system. A new solution has been developed for estimation of drawdown and stream depletion under conditions that are more representative of those in natural systems (finite width stream of shallow penetration adjoining an aquifer of limited lateral extent). This solution shows that the conventional assumption of a fully penetrating stream will lead to a significant overestimation of stream depletion (> 100%) in many practical applications. The degree of overestimation will depend on the value of the stream leakance parameter and the distance from the pumping well to the stream. Although leakance will increase with stream width, a very wide stream will not necessarily be well represented by a model of a fully penetrating stream. The impact of lateral boundaries depends upon the distance from the pumping well to the stream and the stream leakance parameter. In most cases, aquifer width must be on the order of hundreds of stream widths before the assumption of a laterally infinite aquifer is appropriate for stream-depletion calculations. An important assumption underlying this solution is that stream-channel penetration is negligible relative to aquifer thickness. However, an approximate extension to the case of nonnegligible penetration provides reasonable results for the range of relative penetrations found in most natural systems (up to 85%). Since this solution allows consideration of a much wider range of conditions than existing analytical approaches, it could prove to be a valuable new tool for water management design and water rights adjudication purposes.     

Geophysical and Hydrologic Studies of Lake Seepage Variability

Variations in lake seepage were studied along a 130 m shoreline of Mirror Lake NH. Seepage was downward from the lake to groundwater; rates measured from 28 seepage meters varied from 0 to ?282 cm/d. Causes of this variation were investigated using electrical resistivity surveys and lakebed sediment characterization. Two‐dimensional (2D) resistivity surveys showed a transition in lakebed sediments from outwash to till that correlated with high‐ and low‐seepage zones, respectively. However, the 2D survey was not able to predict smaller scale variations within these facies. In the outwash, fast seepage was associated with permeability variations in a thin (2 cm) layer of sediments at the top of the lakebed. In the till, where seepage was slower than that in the outwash, a three‐dimensional resistivity survey mapped a point of high seepage associated with heterogeneity (lower resistivity and likely higher permeability). Points of focused flow across the sediment–water interface are difficult to detect and can transmit a large percentage of total exchange. Using a series of electrical resistivity geophysical methods in combination with hydrologic data to locate heterogeneities that affect seepage rates can help guide seepage meter placement. Improving our understanding of the causes and types of heterogeneity in lake seepage will provide better data for lake budgets and prediction of mass transfer of solutes or contaminants between lakes and groundwater.     

Towed vehicle observations of thin layer structure and a low-salinity intrusion in Northern Monterey Bay,CA

Thin phytoplankton layers are common features in the coastal environment; however sampling these fine-scale optical features across broad horizontal scales remains a challenge. To investigate the horizontal spatial structure of thin phytoplankton layers, we performed an overnight survey in northern Monterey Bay, CA, USA using a SeaSciences Acrobat towed-vehicle. Physical and optical measurements were collected between the surface and near-bottom-depths along four parallel, across-shore transects. Three coherent chlorophyll features were observed: (1) a broad, sub-surface patch at the offshore end, (2) a near-surface patch at the nearshore end, and (3) a deep patch located between the nearshore and offshore patches. The offshore and nearshore patch were separated by a change in seafloor slope and a region of compressed, shoaling isopycnals. Both the offshore and nearshore features were located at the pycnocline, had similar optical properties, and were co-located with a low-salinity intrusion. The deep chlorophyll patch had associated physical and optical properties that were distinct from the patches at the pycnocline. The results from this study further underscore the heterogeneous horizontal spatial structure of thin layers and also add to the growing evidence suggesting that low-salinity intrusions may be strongly linked to the formation of thin phytoplankton layers over the northern shelf of Monterey Bay.     

Effect of seawater on incident plane P and SV waves at ocean bottom and engineering characteristics of offshore ground motion records off the coast of southern California,USA

The effect of seawater on vertical ground motions is studied via a theoretical method and then actual offshore ground motion records are analyzed using a statistical method. A theoretical analysis of the effect of seawater on incident plane P and SV waves at ocean bottom indicate that on one hand, the affected frequency range of vertical ground motions is prominent due to P wave resonance in the water layer if the impedance ratio between the seawater and the underlying medium is large, but it is greatly suppressed if the impedance ratio is small; on the other hand, for the ocean bottom interface model selected herein, vertical ground motions consisting of mostly P waves are more easily affected by seawater than those dominated by SV waves. The statistical analysis of engineering parameters of offshore ground motion records indicate that:(1) Under the infl uence of softer surface soil at the seafl oor, both horizontal and vertical spectral accelerations of offshore motions are exaggerated at long period components, which leads to the peak spectral values moving to a longer period.(2) The spectral ratios(V/H) of offshore ground motions are much smaller than onshore ground motions near the P wave resonant frequencies in the water layer; and as the period becomes larger, the effect of seawater becomes smaller, which leads to a similar V/H at intermediate periods(near 2 s). These results are consistent with the conclusions of Boore and Smith(1999), but the V/H of offshore motion may be smaller than the onshore ground motions at longer periods(more than 5 s).     

Simple scaling properties of Canadian annual average streamflow

The spatial scaling properties of Canadian annual average streamflow (abbreviated as AASF) are assessed using both the product moments (PMs) and the probability weighted moments (PWMs) of AASF across the entire country and in its sub-climatic regions. By the PMs, the log relationship between the kth moments of AASF and the drainage area can be almost represented by a perfect straight line across the entire country and in its sub-climatic regions, whose regression parameters are a linear function of the moment order. By the PWMs, the logarithm of the kth PWM is a linear function of the logarithm of drainage area for the entire country and its sub-climatic regions, where its slope (or scale exponent) in a region is constant and is independent of the order. These results indicate that Canadian AASF exhibits simple scaling and drainage area alone may describe most of the variability in the moments of AASF. The third approach, based on the log linearity between quantiles and drainage area, is applied to Region 2, also demonstrate simple scaling of AASF in that region, as concluded from using PMs and PWMs methods, which indicates that all three methods are consistent. The simple scaling results provide a basis for using the index flood method to conduct regional frequency analysis of AASF in Canada.     

Limitations of ROI Testing for Venting Design: Description of an Alternative Approach Based on Attainment of a Critical Pore-Gas Velocity in Contaminated Media

We utilize data from a Superfund site where radius of influence (ROI) testing was conducted in support of a venting design to describe limitations of ROI evaluation in more detail than has been done previously, and to propose an alternative method of design based on specification and attainment of a critical pore-gas velocity in contaminated subsurface media. Since accurate gas permeability estimation is critical to pore-gas velocity computation, we assess the usefulness of ROI testing data on estimation of radial permeability, vertical permeability, and leakance. We apply information from published studies on rate-limited vapor transport to provide the basis for selection of a critical design pore-gas velocity for soils at this site. Using single-well gas flow simulations, we evaluate whether this critical pore-gas velocity was achieved at measured ROIs. We then conduct a series of multi-well gas flow simulations to assess how variation in anisotropy and leakance affect three-dimensional vacuum and pore-gas velocity profiles and determination of an ROI. Finally, when attempting to achieve a critical design pore-gas velocity we evaluate whether it is more efficient to install additional wells or pump existing wells at a higher flow rate.     

基于混合波场地震动输入技术的近海场地地震反应分析方法

基于人工边界子结构模型,提出一种利用混合波场实现近海场地中地震P波和SV波垂直输入的方法.该方法中用于波动输入的混合波场由计算模型两侧截断边界的自由波场和底面边界的入射波场构成,避免了不规则近海场地的自由波场求解.同时采用基于声流体单元的流固耦合算法模拟场地-海水动力相互作用,利用流体介质人工边界和黏弹性人工边界单元模...     

A semi‐analytical solution for groundwater responses to stream‐stage variations and tidal fluctuations in a coastal aquifer

A two‐dimensional semi‐analytical solution to analyse stream–aquifer interactions in a coastal aquifer where groundwater level responds to tidal effects is presented. The conceptual model considered is a two‐dimensional subsurface system with stream and coastline boundaries at right angles. The dimensional and non‐dimensional boundary value problems were solved for water level in the aquifer by successive application of Laplace and Fourier transform techniques, and the results were obtained by numerical inversion of the transformed solution. The solution was then verified by reducing the solutions to one‐dimensional known problems and comparing the results with those from previous studies. Hypothetical examples were used to examine the characteristics of water‐level variations due to the variations in stream stage and the fluctuations in tide level. Sensitivity analysis indicated that streambed leakance has no influence over the amplitude of groundwater fluctuations, but that the effect of stream stage increases with increasing leakance. Little difference was observed in the water level for different aquifer penetration ratios with narrow stream width. Increases in streambed leakance caused increases in the effect of aquifer penetration by the stream on the water level. An increased specific yield value resulted in decreased amplitude of water fluctuations and mean water level, and showed that water‐level variations due to stream and tidal boundaries are sensitive to specific yield. Copyright © 2006 John Wiley & Sons, Ltd.     

Impact of erosion and accretion on the distribution of enterococci in beach sands

Bacterial pathogens in coastal sediments may pose a health risk to users of beaches. Although recent work shows that beach sands harbor both indicator bacteria and potential pathogens, it is neither known how deep within beach sands the organisms may persist nor if they may be exposed during natural physical processes. In this study, sand cores of approximately 100 cm depth were collected at three sites across the beach face in Kitty Hawk, North Carolina, before, during, and after large waves from an offshore hurricane. The presence of DNA from the fecal indicator bacterium Enterococci was detected in subsamples at different depths within the cores by PCR amplification. Erosion and accretion of beach sand at the three sites were also determined for each sampling day. The results indicate that ocean beach sands with persisting enterococci signals could be exposed and redistributed when wind, waves, and currents cause beach erosion or accretion.     

DIRECT APPLICATIONS OF THE DAR ZARROUK PARAMETERS IN GROUND WATER SURVEYS*

The combination of layer resistivity and thickness in the so-called Dar Zarrouk parameters S and T may be of direct use in aquifer protection studies and for the evaluation of hydrologic properties of aquifers. The protective capacity of a clayey aquifer overburden is proportional to its longitudinal unit conductance S which, in terms of aquifer protection, gets a dimension of time (e.g. infiltration time). Aquifer storage in fissured reservoirs may be determined from differential conductance measurements (ΔS). Combination of the expression for ΔS with an empirical expression for electric conduction in fissured media yields a simple formula for water content per unit surface area. Both principles and possible developments are illustrated for a set of carboniferous limestone basins.