Analytical and Semi‐Analytical Tools for the Design of Oscillatory Pumping Tests

Oscillatory pumping tests—in which flow is varied in a periodic fashion—provide a method for understanding aquifer heterogeneity that is complementary to strategies such as slug testing and constant‐rate pumping tests. During oscillatory testing, pressure data collected at non‐pumping wells can be processed to extract metrics, such as signal amplitude and phase lag, from a time series. These metrics are robust against common sensor problems (including drift and noise) and have been shown to provide information about aquifer heterogeneity. Field implementations of oscillatory pumping tests for characterization, however, are not common and thus there are few guidelines for their design and implementation. Here, we use available analytical solutions from the literature to develop design guidelines for oscillatory pumping tests, while considering practical field constraints. We present two key analytical results for design and analysis of oscillatory pumping tests. First, we provide methods for choosing testing frequencies and flow rates which maximize the signal amplitude that can be expected at a distance from an oscillating pumping well, given design constraints such as maximum/minimum oscillator frequency and maximum volume cycled. Preliminary data from field testing helps to validate the methodology. Second, we develop a semi‐analytical method for computing the sensitivity of oscillatory signals to spatially distributed aquifer flow parameters. This method can be quickly applied to understand the “sensed” extent of an aquifer at a given testing frequency. Both results can be applied given only bulk aquifer parameter estimates, and can help to optimize design of oscillatory pumping test campaigns.     

Interpretation of Pumping Tests in Heterogeneous Aquifers with Constant Head Boundary

This paper investigates the impact of heterogeneity of the transmissivity field on the interpretation of steady-state pumping test data from aquifer systems delimited by constant head boundaries such as aquifers adjacent to lakes or rivers. Spatially variable transmissivity fields are randomly generated and used to simulate the drawdown due to a pumping well located at different distances from a constant head boundary. The steady-state drawdown simulated at different observation wells are then interpreted using the Hantush method (Hantush 1959). The numerical simulations show that, in contrast to the case of infinite aquifer domains, the interpreted transmissivity varies depending on well locations and the separation distance between pumping well and boundary relative to the correlation length. The ensemble-averaged estimated transmissivity varies between the geometric mean and the arithmetic mean, and can even exceed the arithmetic mean in a narrow domain adjacent to the boundary. It approaches the geometric mean of the underlying transmissivity field only if the distance between the pumping well is more than 20 times the characteristic length of the transmissivity field.     

Mathematical Modeling of Displacement Processes in Confined Flows in Coastal Aquifers

Mathematical models of fresh subsoil water flows in a confined aquifer to a sea (pit, pool, and the like) containing salt water are considered. To study these models mixed boundary value problems of the theory of analytical functions are formulated and solved with the use of the Polubarinova-Kochina method. The models were used to develop algorithms for calculating displacement in situations where subsoil water flows discharge into the sea at lateral inflow or inflow from below. The effect of the model physical parameters on the character and extent of displacement was analyzed with the use of obtained exact analytical relationships and numerical calculations. The hydrodynamic structure is described, and specific features of the flows being modeled are established.     

Impact of Sea-Level Rise on Sea Water Intrusion in Coastal Aquifers

Despite its purported importance, previous studies of the influence of sea-level rise on coastal aquifers have focused on specific sites, and a generalized systematic analysis of the general case of the sea water intrusion response to sea-level rise has not been reported. In this study, a simple conceptual framework is used to provide a first-order assessment of sea water intrusion changes in coastal unconfined aquifers in response to sea-level rise. Two conceptual models are tested: (1) flux-controlled systems, in which ground water discharge to the sea is persistent despite changes in sea level, and (2) head-controlled systems, whereby ground water abstractions or surface features maintain the head condition in the aquifer despite sea-level changes. The conceptualization assumes steady-state conditions, a sharp interface sea water-fresh water transition zone, homogeneous and isotropic aquifer properties, and constant recharge. In the case of constant flux conditions, the upper limit for sea water intrusion due to sea-level rise (up to 1.5 m is tested) is no greater than 50 m for typical values of recharge, hydraulic conductivity, and aquifer depth. This is in striking contrast to the constant head cases, in which the magnitude of salt water toe migration is on the order of hundreds of meters to several kilometers for the same sea-level rise. This study has highlighted the importance of inland boundary conditions on the sea-level rise impact. It identifies combinations of hydrogeologic parameters that control whether large or small salt water toe migration will occur for any given change in a hydrogeologic variable.     

Numerical Simulations of Seasonally Oscillated Groundwater Dynamics in Coastal Confined Aquifers

Studies investigating the effects of inland recharge on coastal groundwater dynamics were carried out typically in unconfined aquifers, with few in confined aquifers. This study focused on the groundwater dynamics in confined aquifers with seasonally sinusoidally fluctuated inland groundwater head and constant sea level by numerical simulations. It is known that the mixing zone (MZ) of saltwater wedge in response to the seasonal oscillations of inland groundwater head swings around the steady-state MZ. However, our simulation results indicate that even the most landward freshwater-saltwater interface over a year is seaward from the steady-state location when the hydraulic conductivity K is ≤10⁻⁴ m/s under certain boundary conditions with given parameter values. That is, seasonal oscillations of inland groundwater head may reduce seawater intrusion in confined coastal aquifers when K ≤ 10⁻⁴ m/s. Sensitivity analysis indicates that for aquifers of K ≤ 10⁻⁴ m/s, the larger the inland head fluctuation amplitude is, the less the seawater intrudes. This is probably due to the reason that the seawater intrusion time decreases with the increase of fluctuation amplitude when K ≤ 10⁻⁴ m/s. Numerical simulations demonstrate that seasonal inland groundwater head oscillations promote the annual averaged recirculated seawater discharge across the seaward boundary.     

Electrical Resistivity Imaging and the Saline Water Interface in High-Quality Coastal Aquifers

Population growth and changing climate continue to impact on the availability of natural resources. Urbanization of vulnerable coastal margins can place serious demands on shallow groundwater. Here, groundwater management requires definition of coastal hydrogeology, particularly the seawater interface. Electrical resistivity imaging (ERI) appears to be ideally suited for this purpose. We investigate challenges and drivers for successful electrical resistivity imaging with field and synthetic experiments. Two decades of seawater intrusion monitoring provide a basis for creating a geo-electrical model suitable for demonstrating the significance of acquisition and inversion parameters on resistivity imaging outcomes. A key observation is that resistivity imaging with combinations of electrode arrays that include dipole–dipole quadrupoles can be configured to illuminate consequential elements of coastal hydrogeology. We extend our analysis of ERI to include a diverse set of hydrogeological settings along more than 100 km of the coastal margin passing the city of Perth, Western Australia. Of particular importance are settings with: (1) a classic seawater wedge in an unconfined aquifer, (2) a shallow unconfined aquifer over an impermeable substrate, and (3) a shallow multi-tiered aquifer system over a conductive impermeable substrate. We also demonstrate a systematic increase in the landward extent of the seawater wedge at sites located progressively closer to the highly urbanized center of Perth. Based on field and synthetic ERI experiments from a broad range of hydrogeological settings, we tabulate current challenges and future directions for this technology. Our research contributes to resolving the globally significant challenge of managing seawater intrusion at vulnerable coastal margins.     

A Method to Determine the Formation Constants of Leaky Aquifers, and Its Application to Pumping Test Data

Abstract. A method to calculate aquifer transmissivity, storage coefficient, and the leakage coefficient from pumping test data for a leaky aquifer is presented. The method is carried out by a computer program and is based on a minimization of the sum of squares of differences between drawdown in the observation well and the theoretical values from the Hantush and Jacob formula. No user defined starting points are necessary. Random error estimates for the parameters are given. Applications of the method are illustrated using data from pumping tests performed in leaky aquifers at the Cauca River Valley, Colombia.     

Analytical Solutions of Travel Time to a Pumping Well with Variable Evapotranspiration

Analytical solutions of groundwater travel time to a pumping well in an unconfined aquifer have been developed in previous studies, however, the change in evapotranspiration was not considered. Here, we develop a mathematical model of unconfined flow toward a discharge well with redistribution of groundwater evapotranspiration for travel time analysis. Dependency of groundwater evapotranspiration on the depth to water table is described using a linear formula with an extinction depth. Analytical solutions of groundwater level and travel time are obtained. For a typical hypothetical example, these solutions perfectly agree with the numerical simulation results based on MODFLOW and MODPATH. As indicated in a dimensionless framework, a lumped parameter which is proportional to the pumping rate controls the distributions of groundwater evapotranspiration rate and the travel time along the radial direction.     

Reliable Monitoring of the Transition Zone Between Fresh and Saline Waters in Coastal Aquifers

This study deals with the reliability of monitoring the transition zone between fresh and saline waters in coastal aquifers, considering the effect of tides in long‐perforated boreholes. Electric conductivity (EC) fluctuations in the coastal aquifer of Israel, as measured in long‐perforated borehole, were found to have the same periodicities as the sea tide, though some orders of magnitude larger than sea‐level or groundwater level fluctuations. Direct measurements in the aquifer through buried EC sensors demonstrate that EC measurements within the long‐perforated boreholes might be distorted due to vertical flow in the boreholes, whereas actual fluctuations of the transition zone within the aquifer are some orders of magnitude smaller. Considering these field data, we suggest that monitoring of the transition zone between fresh and saline water adjacent to the sea through long‐perforated boreholes is unreliable. EC fluctuations in short‐perforated boreholes (1 m perforation at the upper part of the transition zone) were somewhat larger than in the aquifer, but much smaller than those in the long‐perforated borehole. The short‐perforation diminishes the vertical flow and the distortion and therefore is more reliable for monitoring the transition zone in the shoreline vicinity.