Electric-hydraulic conductivity correlation in fractured crystalline bedrock: Central Landfill, Rhode Island, USA

Remote sensing and geoelectrical methods were used to find water-bearing fractures in the Scituate granite under the Central Landfill of Rhode Island. These studies were necessary to evaluate the integrity of the sanitary landfill and for planning safe landfill extensions. The most useful results were obtained with fracture trace analysis using Landsat and SLAR imagery in combination with ground-based resistivity measurements using Schlumberger vertical electrical soundings based on the assumption of horizontally layered strata. Test borings and packer tests confirmed, in the presence of a lineament and low bedrock resistivity, the probable existence of high bedrock fracture density and high average hydraulic conductivity. However, not every lineament was found to be associated with high fracture density and high hydraulic conductivity. Lineaments alone are not a reliable basis for characterising a landfill site as being affected by fractured bedrock. Horizontal fractures were found in borings located away from lineaments. High values of hydraulic conductivity were correlated with low bedrock resistivities. Bedrock resistivities between 60 and 700 Ω m were associated with average hydraulic conductivities between 4 and 60 cm/day. In some cases very low resistivities were confined to the upper part of the bedrock where the hydraulic conductivity was very large. These types of fractures apparently become narrower in aperture with depth. Bedrock zones having resistivities greater than 1000 Ω m showed, without exception, no flow to the test wells. Plots of bedrock resistivity versus the average hydraulic conductivity indicate that the resistivity decreases with increasing hydraulic conductivity. This relationship is inverse to that found in most unconsolidated sediments and is useful for estimating the hydraulic conductivity in groundwater surveys in fractured bedrock. In appropriate settings such as the Central Landfill site in New England, this electric-hydraulic correlation relationship, supplemented by lineament trace analysis, can be used effectively to estimate the hydraulic conductivity in bedrock from only a limited number of resistivity depth soundings and test wells.     

Time scales of DNAPL migration in sandy aquifers examined via numerical simulation

The time required for dense nonaqueous phase liquid (DNAPL) to cease migrating following release to the subsurface is a valuable component of a site conceptual model. This study uses numerical simulation to investigate the migration of six different DNAPLs in sandy aquifers. The most influential parameters governing migration cessation time are the density and viscosity of the DNAPL and the mean hydraulic conductivity of the aquifer. Releases of between 1 and 40 drums of chlorinated solvent DNAPLs, characterized by relatively high density and low viscosity, require on the order of months to a few years to cease migrating in a heterogeneous medium sand aquifer having an average hydraulic conductivity of 7.4 x 10(-3) cm/s. In contrast to this, the release of 20 drums of coal tar (rho(D)= 1061 kg/m(3), micro(D)= 0.161 Pa.s) requires more than 100 years to cease migrating in the same aquifer. Altering the mean hydraulic conductivity of the aquifer results in a proportional change in cessation times. Parameters that exhibit relatively little influence on migration time scales are the DNAPL-water interfacial tension, release volume, source capillary pressure, mean aquifer porosity, and ambient ground water hydraulic gradient. This study also demonstrates that low-density DNAPLs (e.g., coal tar) give rise to greater amounts of lateral spreading and greater amounts of pooling on capillary barriers than high-density DNAPLs such as trichloroethylene or tetrachloroethylene.     

Comparison of Measured and Modelled Hydraulic Conductivities of Fractured Sandstone Cores

— A new method for characterising the detailed fracture geometry in sandstone cores is presented. This method is based on the impregnation of samples with coloured resin, without significant disturbance of the fractures. The fractures are made clearly visible by the resin, thus allowing the fracture geometry to be examined digitally. In order to model the bulk hydraulic conductivity, the samples are sectioned serially perpendicular to the flow direction. The hydraulic conductivity of individual sections is estimated by summing the contribution of the matrix and each fracture from the digital data. Finally, the hydraulic conductivity of the bulk sample is estimated by a harmonic average in series along the flow path. Results of this geometrical method are compared with actual physical conductivity values measured from fluid experiments carried out prior to sectioning. The predicted conductivity from the fracture geometry parameters (e.g., fracture aperture, fracture width, fracture length and fracture relative roughness all measured using an optical method) is in good agreement with the independent physical measurements, thereby validating the approach.     

Modeling of Plume Capture by Continuous, Low-Permeability Barriers

Modeling was performed to simulate ground water flow through reactive barriers of lower hydraulic conductivity than the surrounding aquifer to determine the plume capture widths. As a plume approaches such a barrier, it spreads laterally. Therefore, to intercept an entire plume, the barrier must be wider than the upgradient width of the undisturbed plume. The results indicate that, for practical values of barrier thickness and plume width, hydraulic conductivities ten-fold less than that of the aquifer can be accommodated by making the width of the barrier approximately 20% greater than the upgradient width of the plume. Barrier hydraulic conductivities one-hundred-fold less than that of the aquifer may require barrier widths up to twice the width of the upgradient plume for plumes 100 feet wide (33 m) and as little as 1.1 times for plumes 1000 feet wide (325 m). The results presented here lend support to the view that novel emplacement methods that create zones of slightly lower hydraulic conductivity than the native aquifer may be viable alternatives to the excavation-and-backfill approaches which have thus far been used for installing permeable reactive barriers.     

Transport of waste leakage in stratified formations

Contaminant transport in a strongly heterogeneous stratified formation whose log hydraulic conductivity distribution has a variance greater than unity is investigated. Four kinds of waste leakage scenario are studied. They are: (1) continuous waste leakage from landfills; (2) temporal waste leakage from landfills; (3) continuous deep-well injection wastes; and (4) temporal deep-well injection wastes. Ensemble average concentrations and variances of concentration distributions are calculated for the four scenarios. The results in this paper show that when heterogeneity of a formation increases, transport in this formation differs significantly from the linear solutions which assume that the variances of log hydraulic conductivity are less than unity.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Hydraulic conductivity in upland blanket peat: measurement and variability

A key parameter used in wetland hydrological and landform development models is hydraulic conductivity. Head recovery tests are often used to measure hydraulic conductivity, but the calculation techniques are usually confined to rigid soil theory. This is despite reports demonstrating the misapplication of rigid soil theory to non‐rigid soils such as peats. Although values of hydraulic conductivity calculated using compressible techniques have been presented for fenland peats, these data have never, to the authors' knowledge, been compared with such calculations in other peat types. Head recovery tests (slug withdrawal) were performed on piezometers at depths ranging from 10 to 80 cm from the surface on north Pennines blanket peats. Results were obtained using both rigid and compressible soil theories, thus allowing comparison of the two techniques. Compressible soil theory gives values for hydraulic conductivity that are typically a factor of five times less than rigid soil calculations. Hydraulic conductivity is often assumed to decrease with depth in upland peats, but at the study site in the northern Pennines it was not found to vary significantly with depth within the range of peat depths sampled. The variance within depth categories was not significantly different to the variance between depth categories showing that individual peat layers did not have characteristic hydraulic conductivity values. Thus, large lateral and vertical differences in hydraulic conductivity over short distances create problems for modelling but may help account for the high frequency of preferential flow pathways within what is otherwise a low matrix hydraulic conductivity peat. Hydraulic conductivity was found to vary significantly between sampling sites, demonstrating that hillslope‐ or catchment‐scale variability may be more important than plot‐scale variability. Values for compressibility of the peats are also reported. These generally decline with depth, and they also vary significantly between sampling sites. There are implications for the way in which measurements of hydraulic conductivity and other properties of blanket peat are interpreted, as the effects of environmental change in one part of a peat catchment may be very different to those in another. Copyright © 2002 John Wiley & Sons, Ltd.     

Are Visible Fractures Accurate Predictors of Flow and Mass Transport in Fractured Till?

Tracer experiments conducted in the laboratory on undisturbed core samples (<7.3-cm-diameter) have been a standard method for estimating hydraulic and transport properties of fractured till since the 1980s. This study assesses the relationship between visible fractures on the top and bottom of core samples and the resulting hydraulic and mass transport properties of the core. We hypothesized that more visible fractures would indicate the presence of a well-connected fracture network, leading to greater hydraulic conductivity (K) values and earlier chemical breakthrough times. To test this hypothesis, water flow and bromide (Br-) tracer experiments were performed on 10, 16-cm diameter, 16-cm-tall samples of fractured Dows Formation till from central Iowa. Visually identifiable fractures were present on the top and bottom of every sample. Results indicate that the visual identification of fractures does not predict a connected fracture network, as some samples produced breakthrough curves showing rapid first arrival times and shapes characteristic of solute transport in a fractured medium, while others appeared similar to an unfractured medium. No correlation was found between the number of visible fractures and K (Pearson's r = 0.25), or Br- first arrival time (r = −0.33), but a strong negative correlation between K and first arrival time (r = −0.92). Results indicate that the sample volume was not large enough to reliably contain a connected fracture network. Thus, testing large volumes of till at the field scale coupled with fracture-flow modeling likely represents the best approach for estimating hydraulic and mass transport properties for fractured till.     

A constitutive model for water flow in unsaturated fractured rocks

A conceptual model for describing effective saturation in fractured hard rock is presented. The fracture network and the rock matrix are considered as an equivalent continuum medium where each fracture is conceptualized as a porous medium of granular structure and the rock matrix is assumed to be impermeable. The proposed model is based on the representation of a rough‐walled fracture by an equivalent porous medium, which is described using classical constitutive models. A simple closed‐form equation for the effective saturation is obtained when the van Genuchten model is used to describe saturation inside fractures and fractal laws are assumed for both aperture and number of fractures. The relative hydraulic conductivity for the fractured rock is predicted from a simple relation derived by Liu and Bodvarsson. The proposed constitutive model contains three independent parameters, which may be obtained by fitting the proposed effective saturation curve to experimental data. Two of the model parameters have physical meaning and can be identified with the reciprocal of the air entry pressure values in the fractures of minimum and maximum apertures. Effective saturation and relative hydraulic conductivity curves match fairly well the simulated constitutive relations obtained by Liu and Bodvarsson. Copyright © 2008 John Wiley & Sons, Ltd.     

On the Relation Between Saturated Conductivity and Capillary Retention Characteristics

Abstract. A simple closed-form expression relating saturated hydraulic conductivity to the van Genuchten capillary retention model parameters is derived. Application of this equation to an experimental data set shows reasonable agreement between measured and predicted saturated conductivity values. The proposed equation provides a consistent theoretical basis for estimating both saturated and unsaturated hydraulic conductivity from statistical pore structure models.     

Constraints on Upward Migration of Hydraulic Fracturing Fluid and Brine

Recent increases in the use of hydraulic fracturing (HF) to aid extraction of oil and gas from black shales have raised concerns regarding potential environmental effects associated with predictions of upward migration of HF fluid and brine. Some recent studies have suggested that such upward migration can be large and that timescales for migration can be as short as a few years. In this article, we discuss the physical constraints on upward fluid migration from black shales (e.g., the Marcellus, Bakken, and Eagle Ford) to shallow aquifers, taking into account the potential changes to the subsurface brought about by HF. Our review of the literature indicates that HF affects a very limited portion of the entire thickness of the overlying bedrock and therefore, is unable to create direct hydraulic communication between black shales and shallow aquifers via induced fractures. As a result, upward migration of HF fluid and brine is controlled by preexisting hydraulic gradients and bedrock permeability. We show that in cases where there is an upward gradient, permeability is low, upward flow rates are low, and mean travel times are long (often >10⁶ years). Consequently, the recently proposed rapid upward migration of brine and HF fluid, predicted to occur as a result of increased HF activity, does not appear to be physically plausible. Unrealistically high estimates of upward flow are the result of invalid assumptions about HF and the hydrogeology of sedimentary basins.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Site characterization in densely fractured dolomite: comparison of methods

One of the challenges in characterizing fractured-rock aquifers is determining whether the equivalent porous medium approximation is valid at the problem scale. Detailed hydrogeologic characterization completed at a small study site in a densely fractured dolomite has yielded an extensive data set that was used to evaluate the utility of the continuum and discrete-fracture approaches to aquifer characterization. There are two near-vertical sets of fractures at the site; near-horizontal bedding-plane partings constitute a third fracture set. Eighteen boreholes, including five coreholes, were drilled to a depth of approximately 10.6 m. Borehole geophysical logs revealed several laterally extensive horizontal fractures and dissolution zones. Flowmeter and short-interval packer testing identified which of these features were hydraulically important. A monitoring system, consisting of short-interval piezometers and multilevel samplers, was designed to monitor four horizontal fractures and two dissolution zones. The resulting network consisted of >70 sampling points and allowed detailed monitoring of head distributions in three dimensions. Comparison of distributions of hydraulic head and hydraulic conductivity determined by these two approaches suggests that even in a densely fractured-carbonate aquifer, a characterization approach using traditional long-interval monitoring wells is inadequate to characterize ground water movement for the purposes of regulatory monitoring or site remediation. In addition, traditional multiwell pumping tests yield an average or bulk hydraulic conductivity that is not adequate for predicting rapid ground water travel times through the fracture network, and the pumping test response does not appear to be an adequate tool for assessing whether the porous medium approximation is valid.     

Seepage mechanisms and permeability differences between loess and paleosols in the critical zone of the Loess Plateau

The soil in the Loess Plateau has special permeability characteristics due to the alternating distribution of loess and paleosol layers. Using an analysis of the physical properties, microstructure and thermogravimetric analysis of loess and paleosol, this paper examines strata seepage mechanisms in the loess tableland area and considers the applicability of a hydraulic conductivity empirical formula. The analysis shows that hydraulic conductivity attenuation with depth can be represented by a negative exponential model, while hydraulic conductivity values are not normally distributed. The best-fitting models of hydraulic conductivity in the horizontal (K_H) and vertical (K_V) directions are Gaussian models, and both have strong spatial correlations. This study of the difference in permeability between loess and paleosol found that the depositional environment was the dominant causal factor, making the average hydraulic conductivity of most loess layers greater than that of the underlying paleosol layers. Different microstructures between loess and paleosol also confirmed the microscopic explanation in permeability anisotropy and their permeability difference. Thermogravimetric analysis determined temperature ranges for different types of water lost by heat, and then calculated ratios of bound water mass to liquid limit, with an average of 0.768. A modified formula suitable for loess was obtained by integrating the consistency index method and effective porosity ratio model into the hydraulic conductivity empirical formula. Compared with the results of laboratory tests and uncorrected formulas, the modified formula provides a good estimate of strata hydraulic conductivity. Accurate understandings of seepage mechanisms and permeability differences in the loess area are important, promoting ecological restoration and providing scientific guidance for the sustainable development of the Loess Plateau.     

Imaging pre‐existing natural fractures using microseismic data

In hydraulic fracturing treatments, locating not only hydraulic fractures but also any pre‐existing natural fractures and faults in a subsurface reservoir is very important. Hydraulic fractures can be tracked by locating microseismic events, but to identify the locations of natural fractures, an additional technique is required. In this paper, we present a method to image pre‐existing fractures and faults near a borehole with virtual reverse vertical seismic profiling data or virtual single‐well profiling data (limited to seismic reflection data) created from microseismic monitoring using seismic interferometry. The virtual source data contain reflections from natural fractures and faults, and these features can be imaged by applying migration to the virtual source data. However, the imaging zone of fractures in the proposed method is strongly dependent on the geographic extent of the microseismic events and the location and direction of the fracture. To verify our method, we produced virtual reverse vertical seismic profiling and single‐well profiling data from synthetic microseismic data and compared them with data from real sources in the same relative position as the virtual sources. The results show that the reflection travel times from the fractures in the virtual source data agree well with travel times in the real‐source data. By applying pre‐stack depth migration to the virtual source data, images of the natural fractures were obtained with accurate locations. However, the migrated section of the single‐well profiling data with both real and virtual sources contained spurious fracture images on the opposite side of the borehole. In the case of virtual single‐well profiling data, we could produce correct migration images of fractures by adopting directional redatuming for which the occurrence region of microseismic events is divided into several subdivisions, and fractures located only on the opposite side of the borehole are imaged for each subdivision.     

Delineating alluvial aquifer heterogeneity using resistivity and GPR data

Conceptual geological models based on geophysical data can elucidate aquifer architecture and heterogeneity at meter and smaller scales, which can lead to better predictions of preferential flow pathways. The macrodispersion experiment (MADE) site, with >2000 measurements of hydraulic conductivity obtained and three tracer tests conducted, serves as an ideal natural laboratory for examining relationships between subsurface flow characteristics and geophysical attributes in fluvial aquifers. The spatial variation of hydraulic conductivity measurements indicates a large degree of site heterogeneity. To evaluate the usefulness of geophysical methods for better delineating fluvial aquifer heterogeneities and distribution of preferential flow paths, a surface grid of two-dimensional ground penetrating radar (GPR) and direct current (DC) resistivity data were collected. A geological model was developed from these data that delineate four stratigraphic units with distinct electrical and radar properties including (from top to bottom) (1) a meandering fluvial system (MFS); (2) a braided fluvial system (BFS); (3) fine-grained sands; and (4) a clay-rich interval. A paleochannel, inferred by other authors to affect flow, was mapped in the MFS with both DC resistivity and GPR data. The channel is 2 to 4 m deep and, based on resistivity values, is predominantly filled with clay and silt. Comparing previously collected hydraulic conductivity measurements and tracer-plume migration patterns to the geological model indicates that flow primarily occurs in the BFS and that the channel mapped in the MFS has no influence on plume migration patterns.     

Electrokinetic Migration of Nitrate Through Heterogeneous Granular Porous Media

This study investigates and quantifies the influence of physical heterogeneity in granular porous media, represented by materials with different hydraulic conductivity, on the migration of nitrate, used as an amendment to enhance bioremediation, under an electric field. Laboratory experiments were conducted in a bench‐scale test cell under a low applied direct current using glass bead and clay mixes and synthetic groundwater to represent ideal conditions. The experiments included bromide tracer tests in homogeneous settings to deduce controls on electrokinetic transport of inorganic solutes in the different materials, and comparison of nitrate migration under homogeneous and heterogeneous scenarios. The results indicate that physical heterogeneity of subsurface materials, represented by a contrast between a higher‐hydraulic conductivity and lower‐hydraulic conductivity material normal to the direction of the applied electric field exerts the following controls on nitrate migration: (1) a spatial change in nitrate migration rate due to changes in effective ionic mobility and subsequent accumulation of nitrate at the interface between these materials; and (2) a spatial change in the voltage gradient distribution across the hydraulic conductivity contrast, due to the inverse relationship with effective ionic mobility. These factors will contribute to higher mass transport of nitrate through low hydraulic conductivity zones in heterogeneous porous media, relative to homogeneous host materials. Overall electrokinetic migration of amendments such as nitrate can be increased in heterogeneous granular porous media to enhance the in situ bioremediation of organic contaminants present in low hydraulic conductivity zones.     

Landfill Performance Assessment at a Semiarid Site: Modeling and Validation

The HELP model (Version 2.05) was applied to simulate the long-term percolation from a proposed landfill in southern New Mexico. The model predicted percolation would be about 0.0012 in/yr (10^-10 cm/s). This result compared very favorably with independent estimates of recharge at the site which used the chloride mass balance method and hydrogeologic properties. The recharge estimates at this site are also quite similar to values obtained at other sites in New Mexico and west Texas. The long-term percolation through a closed landfill at this site is very small and would be nearly 1000-fold smaller than the saturated hydraulic conductivity of a typical clay liner.     

Effective hydraulic conductivity of nonstationary aquifers

Assuming that the ln hydraulic conductivity in an aquifer is mathematically approximated by a spatial deterministic surface, or trend, plus a stationary random noise, we treat the problem of finding what the effective hydraulic conductivity of that aquifer is. This problem is tackled by spectral methods applied to a type of diffusion equation of groundwater flow, together with suitable coordinate transformations. Analytical (exact) solutions in terms of elementary functions are presented for one- and three-dimensional finite and infinite domains. Stability criteria are obtained for the solutions, in terms of a critical parameter, that turns out to involve the product of correlation scale and trend gradient. For the case of finite and symmetrical domains, additional provisions to insure the stability of numerical calculations of effective hydraulic conductivity are provided. Effective hydraulic conductivity is an important property, with potential applications in the calibrations of groundwater and transport numerical models.