A comparison of peak and plate electrodes in electrical resistivity tomography: application to the chalky groundwater of the Beauvais aquifer (northern part of the Paris basin,France)

Electrical resistivity tomography was used in order to explore an experimental site of the LaSalle Beauvais Polytechnic Institute (France). The test was conducted along a profile line of 315 m length, using 64 electrodes deployed at an inter‐electrode spacing of 5 m, and the data were recorded using gradient, Wenner and pole–dipole arrays. The performance of plate electrodes (non‐conventional flat‐based) is compared with the performance of peak electrodes (conventional spike). The hydrogeophysical investigation of the chalk aquifer system of Beauvais shows that the performance of plate electrodes is satisfactory, leading to inversions of small root‐mean‐square errors. Peak and the plate electrodes were tested before and after injection of a salt tracer in the piezometer of the experimental site. The study demonstrates the usefulness of plate electrodes (efficient, less time consuming) and the possibility of aquifer characterization by a salt tracer. Copyright © 2010 John Wiley & Sons, Ltd.     

A simple constant-head injection test for streambed hydraulic conductivity estimation

A fast, efficient constant-head injection test (CHIT) for in situ estimation of hydraulic conductivity (K) of sandy streambeds is presented. This test uses constant-head hydraulic injection through a manually driven piezometer. Results from CHIT compare favorably to estimates from slug testing and grain-size analysis. The CHIT combines simplicity of field performance, data interpretation, and accuracy of K estimation in flowing streams.     

Depth‐dependent hydraulic conductivity distribution patterns of a streambed

Characterization of streambed hydraulic conductivity from the channel surface to a great depth below the channel surface can provide needed information for the determination of stream‐aquifer hydrologic connectedness, and it is also important to river restoration. However, knowledge on the streambed hydraulic conductivity for sediments 1 m below the channel surface is scarce. This study describes a method that was used to determine the distribution patterns of streambed hydraulic conductivity for sediments from channel surface to a depth of 15 m below. The method includes Geoprobe's direct‐push techniques and Permeameter tests. Direct‐push techniques were used to generate the electrical conductivity (EC) logs and to collect sequences of continuous sediment cores from river channels, as well as from the alluvial aquifer connected to the river. Permeameter tests on these sediment cores give the profiles of vertical hydraulic conductivity (K_v) of the channel sediments and the aquifer materials. This method was applied to produce K_v profiles for a streambed and an alluvial aquifer in the Platte River Valley of Nebraska, USA. Comparison and statistical analysis of the K_v profiles from the river channel and from the proximate alluvial aquifer indicates a special pattern of K_v in the channel sediments. This depth‐dependent pattern of K_v distribution for the channel sediments is considered to be produced by hyporheic processes. This K_v‐distribution pattern implied that the effect of hyporheic processes on streambed hydraulic conductivity can reach the sediments about 9 m below the channel surface. Copyright © 2010 John Wiley & Sons, Ltd.     

The usefulness of outcrop analogue air permeameter measurements for analyzing aquifer heterogeneity: quantifying outcrop hydraulic conductivity and its spatial variability

Saturated hydraulic conductivity (K) is one of the most important parameters determining groundwater flow and contaminant transport in both unsaturated and saturated porous media. Although several well‐established laboratory methods exist for determining K, in situ measurements of this parameter remain very complex and scale dependent. Often, the limited accessibility of subsurface sediments for sampling means an additional impediment to our ability to quantify subsurface K heterogeneity. One potential solution is the use of outcrops as analogues for subsurface sediments. This paper investigates the use of air permeameter measurements on outcrops of unconsolidated sediments to quantify K and its spatial heterogeneity on a broad range of sediment types. The Neogene aquifer in northern Belgium is used as a case study for this purpose. To characterize the variability in K, 511 small‐scale air permeability measurements were performed on outcrop sediments representative over five of the aquifer's lithostratigraphic units. From these measurements, outcrop‐scale equivalent K tensors were calculated using numerical upscaling techniques. Validation of the air permeameter‐based K values by comparison with laboratory constant head K measurements reveals a correlation of 0.93. Overall, the results indicate that hand‐held air permeameters are very efficient and accurate tools to characterize saturated K, as well as its small‐scale variability and anisotropy on a broad range of unconsolidated sediments. The studied outcrops further provided a qualitative understanding of aquifer hydrostratigraphy and quantitative estimates about K variability at the centimetre‐scale to metre‐scale. Copyright © 2013 John Wiley & Sons, Ltd.     

The Impact of the Degree of Aquifer Confinement and Anisotropy on Tidal Pulse Propagation

Oceanic tidal fluctuations which propagate long distances up coastal rivers can be exploited to constrain hydraulic properties of riverbank aquifers. These estimates, however, may be sensitive to degree of aquifer confinement and aquifer anisotropy. We analyzed the hydraulic properties of a tidally influenced aquifer along the Meghna River in Bangladesh using: (1) slug tests combined with drilling logs and surface resistivity to estimate Transmissivity (T); (2) a pumping test to estimate T and Storativity (S) and thus Aquifer Diffusivity (D_PT); and (3) the observed reduction in the amplitude and velocity of a tidal pulse to calculate D using the Jacob‐Ferris analytical solution. Average Hydraulic Conductivity (K) and T estimated with slug tests and borehole lithology were 27.3 m/d and 564 m²/d, respectively. Values of T and S determined from the pumping test ranged from 400 to 500 m²/d and 1 to 5 × 10^?4, respectively with D_PT ranging from 9 to 40 × 10⁵ m²/d. In contrast, D estimated from the Jacob‐Ferris model ranged from 0.5 to 9 × 10⁴ m²/d. We hypothesized this error resulted from deviations of the real aquifer conditions from those assumed by the Jacob‐Ferris model. Using a 2D numerical model tidal pulses were simulated across a range of conditions and D was calculated with the Jacob‐Ferris model. Moderately confined (K_top/K_aquifer < 0.01) or anisotropic aquifers (K_x/K_z > 10) yield D within a factor of 2 of the actual value. The order of magnitude difference in D between pumping test and Jacob‐Ferris model at our site argues for little confinement or anisotropy.     

Maximizing Net Extraction Using an Injection‐Extraction Well Pair in a Coastal Aquifer

In this study, we examine the maximum net extraction rate from the novel arrangement of an injection‐extraction well pair in a coastal aquifer, where fresh groundwater is reinjected through the injection well located between the interface toe and extraction well. Complex potential theory is employed to derive a new analytical solution for the maximum net extraction rate and corresponding stagnation‐point locations and recirculation ratio, assuming steady‐state, sharp‐interface conditions. The injection‐extraction well‐pair system outperforms a traditional single extraction well in terms of net extraction rate for a broad range of well placement and pumping rates, which is up to 50% higher for an aquifer with a thickness of 20 m, hydraulic conductivity of 10 m/d, and fresh water influx of 0.24 m²/d. Sensitivity analyses show that for a given fresh water discharge from an inland aquifer, a larger maximum net extraction is expected in cases with a smaller hydraulic conductivity or a smaller aquifer thickness, notwithstanding physical limits to drawdown at the pumping well that are not considered here. For an extraction well with a fixed location, the optimal net extraction rate linearly increases with the distance between the injection well and the sea, and the corresponding injection rate and recirculation ratio also increase. The analytical analysis in this study provides initial guidance for the design of well‐pair systems in coastal aquifers, and is therefore an extension beyond previous applications of analytical solutions of coastal pumping that apply only to extraction or injection wells.     

Sorption Behaviour of Phenols on Natural Sandy Aquifer Material during Flow‐through Column Experiments: The Effect of pH

Flow‐through column experiments were carried out to investigate the influence of pH on the sorption of three phenols (2‐methyl‐4, 6‐dinitrophenol, 2, 4, 6‐trichlorophenol, pentachlorophenol) onto a natural sandy aquifer material collected from a bank filtration site of River Elbe, Germany. For the phenols investigated, an increase in sorption (retardation) with decreasing pH is observed indicating a stronger sorption of the neutral species in comparison to that of the anions formed by dissociation. The anions of 2‐methyl‐4, 6‐dinitrophenol and 2, 4, 6‐trichlorophenol do not show significant sorption. On the contrary, pentachlorophenol showed sorption not only in neutral form but also in ionic form significantly which should be taken into account while assessing the fate and transport of such compound. A linear model based on the degree of protonation (calculated from pH and pK_a) can be used to resolve the apparent (observed) sorption coefficient (K_{d, app}) into its neutral (K_{d, n}) and ionised (K_{d, i}) components. Knowing pK_a, K_{d, n}, and K_{d, i} the apparent sorption coefficient for pH values other than experimentally investigated can be predicted.     

TABLE DES MATIÈRES DU BULLETIN DE 1956 A 1960 / CONTENTS OF THE BULLETIN 1956–1960

Abstract

Piezometers and wells installed for water quality monitoring are frequently used to assess the saturated hydraulic conductivity (K) in the surrounding formation. A series of recovery tests was conducted to evaluate how purging, required to obtain representative water quality samples, affected measured values of hydraulic conductivity in 15 newly installed and undeveloped piezometers placed to between 2 and 15 m depth (in oxidized and unoxidized material) in a loamy glacial till (K range from 10^?6 to 10^?9 m s^?1). Piezometers were purged between 9 and 11 times for sampling over a period of five months. The effect of the purgings on piezometer development was evaluated by changes in slope of the water level recovery curves which were used to calculate hydraulic conductivity. The first five purgings following piezometer installation increased K in the 15 piezometers by an average of 34%. The average increase in a value of K after 10 purgings was 44%. Values measured for hydraulic conductivity in a 75 mm diameter auger hole appeared stable after four purgings but piezometers installed in larger diameter boreholes (100 mm to 280 mm) snowed increases in K with up to 10 purgings. The hydraulic conductivity determined for piezometers installed at a 30° angle to the vertical showed greater variability than was observed in the adjacent vertically installed piezometers at the same depth.     

Evaluating the quality of hydraulic conductivity estimates from piezometer slug tests in peat

Although widely used in wetland hydrological studies, hydraulic conductivity (K) estimates from piezometer slug tests are often of questionable validity. Frequently, this is because insufficient attention is paid to the details of the test procedure. Further, in a potentially heterogeneous and anisotropic medium such as peat, the use of slug tests is prone to error. In this paper we address some of the methodological issues surrounding piezometer slug tests in peat. We compare slug test data with laboratory determinations of vertical and horizontal K obtained using a new method. Piezometers were installed at three depths in a floodplain fen peat in Norfolk, UK. Slug tests were initiated by both slug insertion and slug withdrawal, and repeat tests were conducted to examine the robustness of our K estimates. Most of the tests displayed departures from the log‐linear model of Hvorslev, the form of departure being consistent with compressible soil behaviour. The results suggest that insertion tests gave similar results to those initiated by withdrawal. Repeat testing showed that withdrawal data, in particular, gave highly reproducible normalized responses that were independent of the initial head. Values for K estimated using the slug tests were in the range 1 × 10⁻⁴ to 1·6 × 10⁻³ cm s⁻¹, which is towards the upper end of the range reported for peats generally. Laboratory tests yielded similar values of K to those obtained from the slug tests. Although the laboratory tests showed that the peat was anisotropic, the K values generated by slug testing proved relatively good estimates of both vertical and horizontal K. Copyright © 2005 John Wiley & Sons, Ltd.     

Hydraulic Tomography: Continuity and Discontinuity of High‐K and Low‐K Zones

Hydraulic tomography is an emerging field and modeling method that provides a continuous hydraulic conductivity (K) distribution for an investigated region. Characterization approaches that rely on interpolation between one‐dimensional (1D) profiles have limited ability to accurately identify high‐K channels, juxtapositions of lenses with high K contrast, and breaches in layers or channels between such profiles. However, locating these features is especially important for groundwater flow and transport modeling, and for design and operation of in situ remediation in complex hydrogeologic environments. We use transient hydraulic tomography to estimate 3D K in a volume of 15‐m diameter by 20‐m saturated thickness in a highly heterogeneous unconfined alluvial (clay to sand‐and‐gravel) aquifer with a K range of approximately seven orders of magnitude at an active industrial site in Assemini, Sardinia, Italy. A modified Levenberg‐Marquardt algorithm was used for geostatistical inversion to deal with the nonlinear nature of the highly heterogeneous system. The imaging results are validated with pumping tests not used in the tomographic inversion. These tests were conducted from three of five clusters of continuous multichannel tubing (CMTs) installed for observation in the tomographic testing. Locations of high‐K continuity and discontinuity, juxtaposition of very high‐K and very low‐K lenses, and low‐K “plugs” are evident in regions of the investigated volume where they likely would not have been identified with interpolation from 1D profiles at the positions of the pumping well and five CMT clusters. Quality assessment methods identified a suspect high‐K feature between the tested volume and a lateral boundary of the model.     

A Correction on Coastal Heads for Groundwater Flow Models

We introduce a simple correction to coastal heads for constant‐density groundwater flow models that contain a coastal boundary, based on previous analytical solutions for interface flow. The results demonstrate that accurate discharge to the sea in confined aquifers can be obtained by direct application of Darcy's law (for constant‐density flow) if the coastal heads are corrected to ((α + 1)/α)h_s ? B/2α, in which h_s is the mean sea level above the aquifer base, B is the aquifer thickness, and α is the density factor. For unconfined aquifers, the coastal head should be assigned the value . The accuracy of using these corrections is demonstrated by consistency between constant‐density Darcy's solution and variable‐density flow numerical simulations. The errors introduced by adopting two previous approaches (i.e., no correction and using the equivalent fresh water head at the middle position of the aquifer to represent the hydraulic head at the coastal boundary) are evaluated. Sensitivity analysis shows that errors in discharge to the sea could be larger than 100% for typical coastal aquifer parameter ranges. The location of observation wells relative to the toe is a key factor controlling the estimation error, as it determines the relative aquifer length of constant‐density flow relative to variable‐density flow. The coastal head correction method introduced in this study facilitates the rapid and accurate estimation of the fresh water flux from a given hydraulic head measurement and allows for an improved representation of the coastal boundary condition in regional constant‐density groundwater flow models.     

Characteristics of preferential flow and groundwater discharge to Shingobee Lake,Minnesota, USA

Small‐scale heterogeneities and large changes in hydraulic gradient over short distances can create preferential groundwater flow paths that discharge to lakes. A 170 m² grid within an area of springs and seeps along the shore of Shingobee Lake, Minnesota, was intensively instrumented to characterize groundwater‐lake interaction within underlying organic‐rich soil and sandy glacial sediments. Seepage meters in the lake and piezometer nests, installed at depths of 0·5 and 1·0 m below the ground surface and lakebed, were used to estimate groundwater flow. Statistical analysis of hydraulic conductivity estimated from slug tests indicated a range from 21 to 4·8 × 10^?3 m day^?1 and small spatial correlation. Although hydraulic gradients are overall upward and toward the lake, surface water that flows onto an area about 2 m onshore results in downward flow and localized recharge. Most flow occurred within 3 m of the shore through more permeable pathways. Seepage meter and Darcy law estimates of groundwater discharge agreed well within error limits. In the small area examined, discharge decreases irregularly with distance into the lake, indicating that sediment heterogeneity plays an important role in the distribution of groundwater discharge. Temperature gradients showed some relationship to discharge, but neither temperature profiles nor specific electrical conductance could provide a more convenient method to map groundwater–lake interaction. These results suggest that site‐specific data may be needed to evaluate local water budget and to protect the water quality and quantity of discharge‐dominated lakes. Copyright © 2002 John Wiley & Sons, Ltd.     

Cost Comparisons of Aquifer Heterogeneity Characterization Methods

The characterization of heterogeneity in hydraulic conductivity (K) is a major challenge for subsurface remediation projects. There are a number of field studies that compare the K estimates obtained using various techniques, but to our knowledge, no field‐based studies exists that compare the performance of estimated K heterogeneity fields or the associated characterization costs. In this paper, we compare the costs of characterizing the three‐dimensional K heterogeneity and its uncertainty estimates of a glaciofluvial aquifer‐aquitard sequence at a 15 m × 15 m × 18 m field site situated on the University of Waterloo campus. We compare geostatistical analysis of high resolution permeameter K data obtained from repacked core samples in five boreholes and hydraulic tomography analysis of four pumping tests consisting of up to 41 monitoring points per test. Aside from the comparison of costs, we also assess the performance of each method by predicting several pumping tests. Our analysis reveals that hydraulic tomography is somewhat more costly than the geostatistical analysis of high resolution permeameter K data due to the higher capital costs associated with the method. However, the equipment may be reused at other sites; hence these costs may be recovered over the life of the equipment. More significantly, hydraulic tomography is able to capture the most important features of the aquifer‐aquitard sequence leading to more accurate predictions of independent pumping tests. This suggests that more robust remediation systems may be designed if site characterization is performed with hydraulic tomography.     

Hydrogeologic Controls on Induced Seismicity in Crystalline Basement Rocks Due to Fluid Injection into Basal Reservoirs

A series of M_b 3.8–5.5 induced seismic events in the midcontinent region, United States, resulted from injection of fluid either into a basal sedimentary reservoir with no underlying confining unit or directly into the underlying crystalline basement complex. The earthquakes probably occurred along faults that were likely critically stressed within the crystalline basement. These faults were located at a considerable distance (up to 10 km) from the injection wells and head increases at the hypocenters were likely relatively small (～70–150 m). We present a suite of simulations that use a simple hydrogeologic‐geomechanical model to assess what hydrogeologic conditions promote or deter induced seismic events within the crystalline basement across the midcontinent. The presence of a confining unit beneath the injection reservoir horizon had the single largest effect in preventing induced seismicity within the underlying crystalline basement. For a crystalline basement having a permeability of 2 × 10^?17 m² and specific storage coefficient of 10^?7/m, injection at a rate of 5455 m³/d into the basal aquifer with no underlying basal seal over 10 years resulted in probable brittle failure to depths of about 0.6 km below the injection reservoir. Including a permeable (k_z = 10^?13 m²) Precambrian normal fault, located 20 m from the injection well, increased the depth of the failure region below the reservoir to 3 km. For a large permeability contrast between a Precambrian thrust fault (10^?12 m²) and the surrounding crystalline basement (10^?18 m²), the failure region can extend laterally 10 km away from the injection well.     

Simulation Assessment of Direct Push Injection Logging for High-Resolution Aquifer Characterization

Direct push injection logging (DPIL) has become one of the most widely used approaches for obtaining vertical profiles of hydraulic conductivity (K) in environmental site investigations. Despite its widespread use, however, there has been no rigorous analysis of the underlying physical processes that take place during DPIL or how the approach would perform under different hydrogeological and operating conditions. We address these issues through a series of numerical simulations. Results show that the ratio of DPIL injection rate over pressure can be used for direct determination of K when K is >10⁻⁶ m/s. When K is <10⁻⁶ m/s and specific storage (Ss) is >10⁻³/m, the ratio becomes increasingly sensitive to Ss; in that case, additional information on Ss is needed for reliable K estimation. For unconsolidated formations of moderate K or higher, the ratio of injection rate over pressure should provide a reasonable K estimate when Ss is <10⁻³/m. Although water injection at previous depths during continuous DPIL has only a small impact on the pressure response measured at the current injection depth, probe advancement can have a significant impact when K and Ss are small. Consequently, in fine-grained materials, the advancement-generated pore water pressure increase can comprise a large portion of the measured pressure response. To diminish the impact of probe advancement in such materials, advancement speed should be kept as low as possible (e.g., 0.5 cm/s).     

A Tube Seepage Meter for In Situ Measurement of Seepage Rate and Groundwater Sampling

We designed and evaluated a “tube seepage meter” for point measurements of vertical seepage rates (q), collecting groundwater samples, and estimating vertical hydraulic conductivity (K) in streambeds. Laboratory testing in artificial streambeds show that seepage rates from the tube seepage meter agreed well with expected values. Results of field testing of the tube seepage meter in a sandy‐bottom stream with a mean seepage rate of about 0.5 m/day agreed well with Darcian estimates (vertical hydraulic conductivity times head gradient) when averaged over multiple measurements. The uncertainties in q and K were evaluated with a Monte Carlo method and are typically 20% and 60%, respectively, for field data, and depend on the magnitude of the hydraulic gradient and the uncertainty in head measurements. The primary advantages of the tube seepage meter are its small footprint, concurrent and colocated assessments of q and K, and that it can also be configured as a self‐purging groundwater‐sampling device.     

Modeling Transient Streaming Potentials in Falling‐Head Permeameter Tests

We present transient streaming potential data collected during falling‐head permeameter tests performed on samples of two sands with different physical and chemical properties. The objective of the work is to estimate hydraulic conductivity (K) and the electrokinetic coupling coefficient (C_l) of the sand samples. A semi‐empirical model based on the falling‐head permeameter flow model and electrokinetic coupling is used to analyze the streaming potential data and to estimate K and C_l. The values of K estimated from head data are used to validate the streaming potential method. Estimates of K from streaming potential data closely match those obtained from the associated head data, with less than 10% deviation. The electrokinetic coupling coefficient was estimated from streaming potential vs. (1) time and (2) head data for both sands. The results indicate that, within limits of experimental error, the values of C_l estimated by the two methods are essentially the same. The results of this work demonstrate that a temporal record of the streaming potential response in falling‐head permeameter tests can be used to estimate both K and C_l. They further indicate the potential for using transient streaming potential data as a proxy for hydraulic head in hydrogeology applications.     

Applicability of spectral analysis to determine hydraulic diffusivity

This study is to evaluate the applicability of estimating the one-dimensional horizontal hydraulic diffusivity of an unconfined aquifer with time-dependent fluctuation of lateral head and vertical recharge boundaries using observed water level spectra. Different models of boundary condition are imposed to evaluate the statistical significance between the calculated hydraulic diffusivity (ξ^′) with the given hydraulic diffusivity (ξ). The auto-spectra of the water level in observation wells tapping the same aquifer are closely related to those at the disturbed boundaries. For an aquifer with a constant hydraulic diffusivity, the water level fluctuation in the monitoring wells is linearly related to the water level spectra observed at the boundaries. The spectral density function of aquifer hydraulic head varies inversely with specific yield (S _y) and directly with recharge. Given small variation in water level spectra at the disturbed boundaries, the water level fluctuation in the aquifer is affected by the recharge condition and the aquifer spectral density function is sensitive to S _y. Using an iterative technique to estimate ξ from 1400 sets of given parameters, 99% of the ξ^′/ξ values deviated within only one order of magnitude with the model length (L) being equal to 1 km and 10 km. For L equal to 100 m, approximately 82% of the ξ^′/ξ population falls within two orders of magnitude. Therefore, spectral analysis of aquifer hydraulic head response can be used to estimate the hydraulic diffusivity of an unconfined aquifer which is affected by periodic variations in recharge and head at boundaries.     

Stochastic analysis of groundwater flow in aquifers with leakage

 Stochastic analysis of one- and two-dimensional flow through a shallow semi-confined aquifer with spatially variable hydraulic conductivity K represented by a stationary (statistically homogeneous) random process is carried out by using the spectral technique. The hydraulic head covariance functions for flows in a semi-confined aquifer bounded by a leaky layer above and an impervious stratum below are derived by assuming that the randomness forcing the head variation to originate from the hydraulic conductivity field of the aquifer. The head covariance functions are studied using two convenient forms of the logarithmic hydraulic conductivity process. The results demonstrate the significant reduction in the head variances and covariances due to the presence of a leaky layer. The hydraulic head correlation distance is also reduced greatly due to the presence of the leaky layer.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 John Wiley & Sons, Ltd.