10.3 Groundwater exploration and hydrogeophysics

Electrical Resistivity and Electromagnetic Methods for Groundwater Exploration

Groundwater exploration depends on geophysical methods to map subsurface structures and locate water resources without drilling. Electrical resistivity, electromagnetic surveys, seismic techniques, and ground-penetrating radar each reveal different aquifer properties and groundwater dynamics. These same methods are also critical for tracking contaminant plumes and evaluating remediation efforts.

Principles of Electrical Resistivity Methods

Electrical resistivity methods work by injecting current into the ground through electrodes and measuring the resulting voltage differences. The measured resistance depends on lithology, porosity, and water content of subsurface materials.

The resistivity of rocks and sediments is primarily controlled by the amount and salinity of pore water. This relationship is described by Archie's Law, which relates bulk resistivity to porosity, water saturation, and pore fluid resistivity. In practice, this means:

• Saturated zones exhibit lower resistivity than unsaturated zones because pore water conducts current far better than rock matrix or air.
• Saline groundwater drives resistivity even lower, while fresh groundwater in clean sands may still show moderate resistivity.
• Clay-rich materials also reduce resistivity due to surface conduction along clay mineral surfaces, independent of pore water salinity.

Resistivity surveys use various electrode configurations, each with trade-offs:

• Wenner array: Good vertical resolution, straightforward field setup, but relatively low lateral resolution.
• Schlumberger array: Better depth sounding capability with fewer electrode moves; commonly used for vertical electrical sounding (VES).
• Dipole-dipole array: Strong lateral resolution, making it well-suited for 2D imaging, but more sensitive to noise at large separations.

Increasing the electrode spacing increases the depth of investigation, but at the cost of reduced resolution.

Principles of Electromagnetic Methods

Electromagnetic (EM) methods induce electromagnetic fields in the subsurface and measure the secondary response. Unlike resistivity methods, EM techniques don't require galvanic contact with the ground, which makes them faster to deploy over large areas.

Two main categories exist:

• Frequency-domain EM (FDEM): Transmits a continuous signal at one or more frequencies. Lower frequencies penetrate deeper. Well-suited for rapid mapping of lateral conductivity variations.
• Time-domain EM (TDEM): Transmits a pulse and measures the decay of the secondary field over time. Later time gates sample deeper. Particularly effective for detecting conductive targets at depth.

EM methods are especially sensitive to conductive materials. Clay-rich sediments and saline groundwater produce strong responses, making EM useful for mapping aquitards, saltwater intrusion zones, and the boundaries between fresh and saline water. The depth of investigation depends on transmitter-receiver separation (FDEM), measurement time window (TDEM), and the conductivity structure itself: highly conductive near-surface layers limit penetration into deeper targets.

Interpreting Geophysical Data for Groundwater Resources

Inverting Geophysical Data for Subsurface Models

Raw geophysical measurements (apparent resistivity, EM voltages) don't directly show you the subsurface structure. They need to be inverted to produce models of true resistivity or conductivity distribution. The inversion process works in these steps:

1. Collect field data along profiles or in grid patterns.
2. Define an initial subsurface model (a grid of cells, each with an assigned resistivity).
3. Use forward modeling to predict what the data would look like for that model.
4. Compare predicted data to observed data and calculate the misfit.
5. Iteratively adjust the model to minimize the misfit, subject to regularization constraints that keep the model geologically reasonable.
6. Evaluate the final model for resolution and reliability, particularly at depth where sensitivity decreases.

The resulting 2D or 3D resistivity models can then be interpreted in hydrogeological terms:

• Low-resistivity zones (typically < 10–50 Ω·m, depending on geology) may indicate saturated, permeable aquifers or clay-rich units.
• High-resistivity zones (hundreds to thousands of Ω·m) often represent unsaturated material, clean bedrock, or low-permeability formations.
• Conductive anomalies in EM data could mean clay-rich aquitards or saline groundwater. Distinguishing between these requires additional information.

This ambiguity is a key challenge: a low-resistivity anomaly could be a productive freshwater aquifer in saturated sand, a clay aquitard, or a saline zone. You can't tell from resistivity alone.

Integrating Geophysical Data with Other Information

Because of the non-uniqueness problem described above, geophysical data should always be integrated with:

• Borehole logs (lithological descriptions, geophysical well logs) to calibrate resistivity values against known geology.
• Hydraulic test data (pumping tests, slug tests) to relate geophysical models to actual hydraulic conductivity and transmissivity.
• Water quality data to distinguish salinity-driven conductivity from clay-driven conductivity.

Time-lapse geophysical surveys add a temporal dimension. By repeating measurements over weeks, months, or seasons, you can monitor:

• Changes in groundwater levels (the water table boundary shifts in resistivity images).
• Seasonal or pumping-induced salinity changes.
• Aquifer storage variations, which inform sustainability assessments.

Time-lapse approaches are powerful because even if the absolute resistivity is ambiguous, changes in resistivity over time are directly linked to changes in saturation or water chemistry.

Seismic and Ground-Penetrating Radar in Hydrogeological Studies

Seismic Methods for Hydrogeological Characterization

Seismic methods measure how elastic waves propagate through the subsurface. Wave velocities depend on the density and elastic moduli of rocks and sediments, which in turn relate to lithology, porosity, and saturation.

Two main approaches are used in hydrogeological work:

• Seismic refraction: Measures the travel times of waves refracted along layer boundaries. Particularly useful for mapping the depth to bedrock or the water table, because seismic velocity increases sharply at these interfaces. P-wave velocity jumps from roughly 200–800 m/s in unsaturated sediments to 1,500+ m/s below the water table (since water is nearly incompressible).
• Seismic reflection: Records waves reflected from subsurface interfaces. Provides detailed images of stratigraphy, bedrock topography, and structural features like faults or buried channels that may control groundwater flow paths.

Seismic methods complement resistivity and EM surveys because they respond to mechanical properties rather than electrical properties. A clay layer and a sand aquifer might have similar seismic velocities but very different resistivities, so combining both datasets reduces interpretation ambiguity.

Ground-Penetrating Radar for Shallow Subsurface Imaging

Ground-penetrating radar (GPR) transmits high-frequency electromagnetic pulses (typically 10 MHz to 1 GHz) into the ground and records reflections from boundaries where the dielectric permittivity changes. The dielectric permittivity of a material is strongly influenced by water content, since water has a much higher permittivity ($\epsilon_r \approx 80$) than most minerals ($\epsilon_r \approx 4\text{–}8$) or air ($\epsilon_r = 1$).

GPR is effective for mapping:

• Shallow aquifer geometry and the water table surface.
• Bedrock surfaces and sedimentary structures (e.g., channel fills, gravel lenses).
• Subsurface utilities, cavities, or buried infrastructure that may affect groundwater flow.

The trade-off with GPR is between resolution and penetration depth. Higher antenna frequencies give finer resolution but attenuate faster. In conductive materials (clay-rich soils, saline groundwater), GPR signal is rapidly absorbed, limiting penetration to less than a meter in some cases. In resistive materials (dry sand, gravel, limestone), penetration can reach 10–30 m or more.

Geophysics in Contaminant Transport and Remediation

Characterizing Subsurface Properties for Contaminant Transport

Contaminant transport through the subsurface is governed by porosity, permeability, and hydraulic conductivity. Geophysical methods can characterize these properties indirectly:

• Electrical resistivity and EM surveys detect conductive contaminant plumes. Leachate from landfills, industrial effluent, or road salt runoff typically has elevated dissolved solids, making the plume more conductive than the surrounding groundwater. Repeated surveys can map the extent and migration direction of the plume over time.
• Seismic and GPR methods identify preferential flow pathways such as fracture zones, faults, or high-permeability sand and gravel channels. These features can cause contaminants to travel much faster and farther than predicted by simple uniform-flow models.

This information directly guides practical decisions: where to place monitoring wells, how to orient extraction wells, and where barriers or reactive zones should be installed.

Monitoring Remediation Efforts with Geophysical Methods

Time-lapse geophysical monitoring is one of the most valuable applications in remediation. Rather than relying solely on point measurements from monitoring wells, geophysical surveys provide spatially continuous information about what's happening between wells.

Specific applications include:

• Pump-and-treat systems: Resistivity monitoring can track whether the conductive plume is shrinking as contaminated water is extracted and treated.
• In-situ bioremediation: Microbial activity changes pore water chemistry (pH, dissolved gases, ionic strength), which can alter the bulk resistivity. Time-lapse surveys can map where bioremediation is active.
• Chemical oxidation/reduction: Injected amendments (permanganate, zero-valent iron, etc.) change the electrical properties of the subsurface. Geophysical monitoring can track the distribution of injected reagents and verify they've reached the target zone.

Integrating geophysical monitoring data with numerical hydrogeological models improves predictions of contaminant fate and transport. This feedback loop supports adaptive management: if the geophysical data show the plume migrating in an unexpected direction, the remediation strategy can be adjusted accordingly.