Modeling of Methane Migration in Shallow Aquifers from Shale Gas Well Drilling

The vertical portion of a shale gas well, known as the “tophole” is often drilled using an air‐hammer bit that may introduce pressures as high as 2400 kPa (350 psi) into groundwater while penetrating shallow aquifers. A 3‐D TOUGH2 model was used to simulate the flow of groundwater under the high hydraulic heads that may be imposed by such trapped compressed air, based on an observed case in West Virginia (USA) in 2012. The model realizations show that high‐pressure air trapped in aquifers may cause groundwater to surge away from the drill site at observable velocities. If dissolved methane is present within the aquifer, the methane can be entrained and transported to a maximum distance of 10.6 m per day. Results from this study suggest that one cause of the reported increase in methane concentrations in groundwater near shale gas production wells may be the transport of pre‐existing methane via groundwater surges induced by air drilling, not necessarily direct natural gas leakage from the unconventional gas reservoir. The primary transport mechanisms are advective transport of dissolved methane with water flow, and diffusive transport of dissolved methane.     

Tracking Groundwater Discharge to a Large River using Tracers and Geophysics

Few studies have investigated large reaches of rivers in which multiple sources of groundwater are responsible for maintaining baseflow. This paper builds upon previous work undertaken along the Fitzroy River, one of the largest perennial river systems in north‐western Australia. Synoptic regional‐scale sampling of both river water and groundwater for a suite of environmental tracers (⁴He, ⁸⁷Sr/⁸⁶Sr, ²²²Rn and major ions), and subsequent modeling of tracer behavior in the river, has enabled definition and quantification of groundwater input from at least three different sources. We show unambiguous evidence of both shallow “local” groundwater, possibly recharged to alluvial aquifers beneath the adjacent floodplain during recent high‐flow events, and old “regional” groundwater introduced via artesian flow from deep confined aquifers. We also invoke hyporheic exchange and either bank return flow or parafluvial flow to account for background ²²²Rn activities and anomalous chloride trends along river reaches where there is no evidence of the local or regional groundwater inputs. Vertical conductivity sections acquired through an airborne electromagnetic (AEM) survey provide insights to the architecture of the aquifers associated with these sources and general groundwater quality characteristics. These data indicate fresh groundwater from about 300 m below ground preferentially discharging to the river, at locations consistent with those inferred from tracer data. The results demonstrate how sampling rivers for multiple environmental tracers of different types—including stable and radioactive isotopes, dissolved gases and major ions—can significantly improve conceptualization of groundwater—surface water interaction processes, particularly when coupled with geophysical techniques in complex hydrogeological settings.     

Net Inflow: An Important Target on the Path to Aquifer Sustainability

Aquifers supporting irrigated agriculture are a resource of global importance. Many of these systems, however, are experiencing significant pumping-induced stress that threatens their continued viability as a water source for irrigation. Reductions in pumping are often the only option to extend the lifespans of these aquifers and the agricultural production they support. The impact of reductions depends on a quantity known as “net inflow” or “capture.” We use data from a network of wells in the western Kansas portions of the High Plains aquifer in the central United States to demonstrate the importance of net inflow, how it can be estimated in the field, how it might vary in response to pumping reductions, and why use of “net inflow” may be preferred over “capture” in certain contexts. Net inflow has remained approximately constant over much of western Kansas for at least the last 15 to 25 years, thereby allowing it to serve as a target for sustainability efforts. The percent pumping reduction required to reach net inflow (i.e., stabilize water levels for the near term [years to a few decades]) can vary greatly over this region, which has important implications for groundwater management. However, the reduction does appear practically achievable (less than 30%) in many areas. The field-determined net inflow can play an important role in calibration of regional groundwater models; failure to reproduce its magnitude and temporal variations should prompt further calibration. Although net inflow is a universally applicable concept, the reliability of field estimates is greatest in seasonally pumped aquifers.     

Palaeosol Control of Arsenic Pollution: The Bengal Basin in West Bengal,India

Groundwater in the Bengal Basin is badly polluted by arsenic (As) which adversely affects human health. To provide low‐As groundwater for As mitigation, it was sought across 235 km² of central West Bengal, in the western part of the basin. By drilling 76 boreholes and chemical analysis of 535 water wells, groundwater with <10 µg/L As in shallow aquifers was found under one‐third of a study area. The groundwater is in late Pleistocene palaeo‐interfluvial aquifers of weathered brown sand that are capped by a palaeosol of red clay. The aquifers form two N‐S trending lineaments that are bounded on the east by an As‐polluted deep palaeo‐channel aquifer and separated by a shallower palaeo‐channel aquifer. The depth to the top of the palaeo‐interfluvial aquifers is mostly between 35 and 38 m below ground level (mbgl). The palaeo‐interfluvial aquifers are overlain by shallow palaeo‐channel aquifers of gray sand in which groundwater is usually As‐polluted. The palaeosol now protects the palaeo‐interfluvial aquifers from downward migration of As‐polluted groundwater in overlying shallow palaeo‐channel aquifers. The depth to the palaeo‐interfluvial aquifers of 35 to 38 mbgl makes the cost of their exploitation affordable to most of the rural poor of West Bengal, who can install a well cheaply to depths up to 60 mbgl. The protection against pollution afforded by the palaeosol means that the palaeo‐interfluvial aquifers will provide a long‐term source of low‐As groundwater to mitigate As pollution of groundwater in the shallower, heavily used, palaeo‐channel aquifers. This option for mitigation is cheap to employ and instantly available.     

Managing Groundwater to Ensure Ecosystem Function

Groundwater is a critical resource not only for human communities but also for many terrestrial, riparian, and aquatic ecosystems and species. Yet groundwater planning and management decisions frequently ignore or inadequately address the needs of these natural systems. As a consequence, ecosystems dependent on groundwater have been threatened, degraded, or eliminated, especially in arid regions. There is growing acknowledgment that governmental protections for these ecological resources are necessary, but current legal, regulatory and voluntary provisions are often inadequate. Groundwater management premised on “safe yield,” which aims to balance human withdrawals with natural recharge rates, typically provides little to no consideration for water needed by ecosystems. Alternatively, the “sustainable yield” concept aims to integrate social, economic and environmental needs for groundwater, but the complexity of groundwater systems creates substantial uncertainty about the impact that current or future groundwater withdrawals will have on ecosystems. Regardless of the legal or regulatory framework, guidance is needed to help ensure environmental water needs will be met, especially in the face of pressure to increase human uses of groundwater resources. In this paper, we describe minimum provisions for planning, managing, and monitoring groundwater that collectively can lower the risk of harm to groundwater-dependent ecosystems and species, with a special emphasis on arid systems, where ecosystems and species may be especially reliant upon and sensitive to groundwater dynamics.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ¹⁸O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.     

A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

Quantifying River‐Groundwater Interactions of New Zealand's Gravel‐Bed Rivers: The Wairau Plain

New Zealand's gravel‐bed rivers have deposited coarse, highly conductive gravel aquifers that are predominantly fed by river water. Managing their groundwater resources is challenging because the recharge mechanisms in these rivers are poorly understood and recharge rates are difficult to predict, particularly under a more variable future climate. To understand the river‐groundwater exchange processes in gravel‐bed rivers, we investigate the Wairau Plain Aquifer using a three‐dimensional groundwater flow model which was calibrated using targeted field observations, “soft” information from experts of the local water authority, parameter regularization techniques, and the model‐independent parameter estimation software PEST. The uncertainty of simulated river‐aquifer exchange flows, groundwater heads, spring flows, and mean transit times were evaluated using Null‐space Monte‐Carlo methods. Our analysis suggests that the river is hydraulically perched (losing) above the regional water table in its upper reaches and is gaining downstream where marine sediments overlay unconfined gravels. River recharge rates are on average 7.3 m³/s, but are highly dynamic in time and variable in space. Although the river discharge regularly hits 1000 m³/s, the net exchange flow rarely exceeds 12 m³/s and seems to be limited by the physical constraints of unit‐gradient flux under disconnected rivers. An important finding for the management of the aquifer is that changes in aquifer storage are mainly affected by the frequency and duration of low‐flow periods in the river. We hypothesize that the new insights into the river‐groundwater exchange mechanisms of the presented case study are transferable to other rivers with similar characteristics.     

Low‐Resolution Modeling of Dense Drainage Networks in Confining Layers

Groundwater‐surface water (GW‐SW) interaction in numerical groundwater flow models is generally simulated using a Cauchy boundary condition, which relates the flow between the surface water and the groundwater to the product of the head difference between the node and the surface water level, and a coefficient, often referred to as the “conductance.” Previous studies have shown that in models with a low grid resolution, the resistance to GW‐SW interaction below the surface water bed should often be accounted for in the parameterization of the conductance, in addition to the resistance across the surface water bed. Three conductance expressions that take this resistance into account were investigated: two that were presented by Mehl and Hill (2010) and the one that was presented by De Lange (1999). Their accuracy in low‐resolution models regarding salt and water fluxes to a dense drainage network in a confined aquifer system was determined. For a wide range of hydrogeological conditions, the influence of (1) variable groundwater density; (2) vertical grid discretization; and (3) simulation of both ditches and tile drains in a single model cell was investigated. The results indicate that the conductance expression of De Lange (1999) should be used in similar hydrogeological conditions as considered in this paper, as it is better taking into account the resistance to flow below the surface water bed. For the cases that were considered, the influence of variable groundwater density and vertical grid discretization on the accuracy of the conductance expression of De Lange (1999) is small.     

Contribution of baseflow nitrate export to non-point source pollution

As a common pollutant of nitrogen in groundwater, nitrate contamination has become a major concern worldwide. Baseflow, one of the dominant hydrological pathways for nitrate migration to streamflow, has been confirmed as a leading nitrate source for stream water where groundwater or subsurface flow contaminated heavily by nitrate. That is, sufficient improvements of water quality may not be attained without proper management for baseflow, even if non-point sources (NPS) pollutants discharged through surface runoff are being well managed. This article reviews the primary nitrate sources, the main factors affecting its transport, and the methodologies for baseflow nitrate estimation, to give some recommendations for future works, including: (1) giving sufficient consideration for the effects of climatological, morphological, and geological factors on baseflow recessions to obtain more reliable and accurate baseflow separation; (2) trying to solve calibration and validation problems for baseflow loads determining in storm flow period; (3) developing a simple and convenient algorithm with certain physics that can be used to separate baseflow NPS pollution from the total directly in different regions, for a reliable estimation of baseflow NPS pollution at larger scale (e.g., national scale); (4) improving groundwater quality simulation module of existing NPS pollution models to have a better simulation for biogeochemical processes in shallow aquifers; (5) taking integrated measures of “source control”, “process interception” and “end remediation” to prevent and control NPS nitrate pollution effectively, not just only the strict control of nutrients loss from surface runoff.     

Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Simulating Stable Isotope Ratios in Plumes of Groundwater Pollutants with BIOSCREEN‐AT‐ISO

BIOSCREEN is a well‐known simple tool for evaluating the transport of dissolved contaminants in groundwater, ideal for rapid screening and teaching. This work extends the BIOSCREEN model for the calculation of stable isotope ratios in contaminants. A three‐dimensional exact solution of the reactive transport from a patch source, accounting for fractionation by first‐order decay and/or sorption, is used. The results match those from a previously published isotope model but are much simpler to obtain. Two different isotopes may be computed, and dual isotope plots can be viewed. The dual isotope assessment is a rapidly emerging new approach for identifying process mechanisms in aquifers. Furthermore, deviations of isotope ratios at specific reactive positions with respect to “bulk” ratios in the whole compound can be simulated. This model is named BIOSCREEN‐AT‐ISO and will be downloadable from the journal homepage.     

Boundary depletion rate and drawdown in leaky wedge‐shaped aquifers

This study presents analytical solutions of the three‐dimensional groundwater flow to a well in leaky confined and leaky water table wedge‐shaped aquifers. Leaky wedge‐shaped aquifers with and without storage in the aquitard are considered, and both transient and steady‐state drawdown solutions are derived. Unlike the previous solutions of the wedge‐shaped aquifers, the leakages from aquitard are considered in these solutions and unlike similar previous work for leaky aquifers, leakage from aquitards and from the water table are treated as the lower and upper boundary conditions. A special form of finite Fourier transforms is used to transform the z‐coordinate in deriving the solutions. The leakage induced by a partially penetrating pumping well in a wedge‐shaped aquifer depends on aquitard hydraulic parameters, the wedge‐shaped aquifer parameters, as well as the pumping well parameters. We calculate lateral boundary dimensionless flux at a representative line and investigate its sensitivity to the aquitard hydraulic parameters. We also investigate the effects of wedge angle, partial penetration, screen location and piezometer location on the steady‐state dimensionless drawdown for different leakage parameters. Results of our study are presented in the form of dimensionless flux‐dimensionless time and dimensionless drawdown‐leakage parameter type curves. The results are useful for evaluating the relative role of lateral wedge boundaries and leakage source on flow in wedge‐shaped aquifers. This is very useful for water management problems and for assessing groundwater pollution. The presented analytical solutions can also be used in parameter identification and in calculating stream depletion rate and volume. Copyright © 2011 John Wiley & Sons, Ltd.     

Predicting Aquifer Response Time for Application in Catchment Modeling

It is well established that changes in catchment land use can lead to significant impacts on water resources. Where land‐use changes increase evapotranspiration there is a resultant decrease in groundwater recharge, which in turn decreases groundwater discharge to streams. The response time of changes in groundwater discharge to a change in recharge is a key aspect of predicting impacts of land‐use change on catchment water yield. Predicting these impacts across the large catchments relevant to water resource planning can require the estimation of groundwater response times from hundreds of aquifers. At this scale, detailed site‐specific measured data are often absent, and available spatial data are limited. While numerical models can be applied, there is little advantage if there are no detailed data to parameterize them. Simple analytical methods are useful in this situation, as they allow the variability in groundwater response to be incorporated into catchment hydrological models, with minimal modeling overhead. This paper describes an analytical model which has been developed to capture some of the features of real, sloping aquifer systems. The derived groundwater response timescale can be used to parameterize a groundwater discharge function, allowing groundwater response to be predicted in relation to different broad catchment characteristics at a level of complexity which matches the available data. The results from the analytical model are compared to published field data and numerical model results, and provide an approach with broad application to inform water resource planning in other large, data‐scarce catchments.     

Aquifer Vulnerability Assessment Based on Sequence Stratigraphic and 39Ar Transport Modeling

A large‐scale groundwater flow and transport model is developed for a deep‐seated (100 to 300 m below ground surface) sedimentary aquifer system. The model is based on a three‐dimensional (3D) hydrostratigraphic model, building on a sequence stratigraphic approach. The flow model is calibrated against observations of hydraulic head and stream discharge while the credibility of the transport model is evaluated against measurements of ³⁹Ar from deep wells using alternative parameterizations of dispersivity and effective porosity. The directly simulated 3D mean age distributions and vertical fluxes are used to visualize the two‐dimensional (2D)/3D age and flux distribution along transects and at the top plane of individual aquifers. The simulation results are used to assess the vulnerability of the aquifer system that generally has been assumed to be protected by thick overlaying clayey units and therefore proposed as future reservoirs for drinking water supply. The results indicate that on a regional scale these deep‐seated aquifers are not as protected from modern surface water contamination as expected because significant leakage to the deeper aquifers occurs. The complex distribution of local and intermediate groundwater flow systems controlled by the distribution of the river network as well as the topographical variation (Tóth 1963) provides the possibility for modern water to be found in even the deepest aquifers.     

Investigating Groundwater Flow Between Edwards and Trinity Aquifers in Central Texas

Understanding the nature of communication between aquifers can be challenging when using traditional physical and geochemical groundwater sampling approaches. This study uses two multiport wells completed within Edwards and Trinity aquifers in central Texas to determine the degree of groundwater inter‐flow between adjacent aquifers. Potentiometric surfaces, hydraulic conductivities, and groundwater major ion concentrations and Sr isotope values were measured from multiple zones within three hydrostratigraphic units (Edwards and Upper and Middle Trinity aquifers). Physical and geochemical data from the multiport wells were combined with historical measurements of groundwater levels and geochemical compositions from the region to characterize groundwater flow and identify controls on the geochemical compositions of the Edwards and Trinity aquifers. Our results suggest that vertical groundwater flow between Edwards and Middle Trinity aquifers is likely limited by low permeability, evaporite‐rich units within the Upper and Middle Trinity. Potentiometric surface levels in both aquifers vary with changes in wet vs. dry conditions, indicating that recharge to both aquifers occurs through distinct recharge areas. Geochemical compositions in the Edwards, Upper, and Middle Trinity aquifers are distinct and likely reflect groundwater interaction with different lithologies (e.g., carbonates, evaporites, and siliceous sediments) as opposed to mixing of groundwater between the aquifers. These results have implications for the management of these aquifers as they indicate that, under current conditions, pumping of either aquifer will likely not induce vertical cross‐formational flow between the aquifers. Inter‐flow between the Trinity and the Edwards aquifers, however, should be reevaluated as pumping patterns and hydrogeologic conditions change.     

Influence of Seasonal Variations in Sea Level on the Salinity Regime of a Coastal Groundwater–Fed Wetland

Seasonal variations in sea level are often neglected in studies of coastal aquifers; however, they may have important controls on processes such as submarine groundwater discharge, sea water intrusion, and groundwater discharge to coastal springs and wetlands. We investigated seasonal variations in salinity in a groundwater‐fed coastal wetland (the RAMSAR listed Piccaninnie Ponds in South Australia) and found that salinity peaked during winter, coincident with seasonal sea level peaks. Closer examination of salinity variations revealed a relationship between changes in sea level and changes in salinity, indicating that sea level–driven movement of the fresh water‐sea water interface influences the salinity of discharging groundwater in the wetland. Moreover, the seasonal control of sea level on wetland salinity seems to override the influence of seasonal recharge. A two‐dimensional variable density model helped validate this conceptual model of coastal groundwater discharge by showing that fluctuations in groundwater salinity in a coastal aquifer can be driven by a seasonal coastal boundary condition in spite of seasonal recharge/discharge dynamics. Because seasonal variations in sea level and coastal wetlands are ubiquitous throughout the world, these findings have important implications for monitoring and management of coastal groundwater–dependent ecosystems.