Surface-based Ground-penetrating Radar Research Articles

GPR-inferred fracture aperture widening in response to a high-pressure tracer injection test at the Äspö Hard Rock Laboratory, Sweden

We assess the performance of the Ground Penetrating Radar (GPR) method in fractured rock formations of very low transmissivity (e.g. T ≈ 10−9–10−10 m2/s for sub-mm apertures) and, more specifically, to image fracture widening induced by high-pressure injections. A field-scale experiment was conducted at the Äspö Hard Rock Laboratory (Sweden) in a tunnel situated at 410 m depth. The tracer test was performed within the most transmissive sections of two boreholes separated by 4.2 m. The electrically resistive tracer solution composed of deionized water and Uranine was expected to lead to decreasing GPR reflections with respect to the saline in situ formation water. The injection pressure was 5000 kPa leading to an injection rate of 8.6 mL/min (at steady state) that was maintained during 25 h, which resulted in a total injected volume of 13 L. To evaluate the fracture pathways between the boreholes, we conducted 3-D surface-based GPR surveys before and at the end of the tracer tests, using 160 MHz and 450 MHz antennas. Difference GPR data between the two acquisitions highlight an increasing fracture reflectivity in-between the boreholes at depths corresponding to the injection interval. GPR-based modeling suggests that the observed increasing reflectivity is not due to the tracer solution, but rather to a 50% widening of the fracture. Considering prevailing uncertainties in material properties, a hydromechanical analysis suggests that such a degree of widening is feasible. This research demonstrates that field-scale in situ GPR experiments may provide constraints on fracture widening by high-pressure injection and could help to constrain field-scale elastic parameters in fractured rock.

Ground-penetrating radar monitoring of fast subsurface processes

Earth and environmental sciences rely on detailed information about subsurface processes. Whereas geophysical techniques typically provide highly resolved spatial images, monitoring subsurface processes is often associated with enormous effort and, therefore, is usually limited to point information in time or space. Thus, the development of spatial and temporal continuous field monitoring methods is a major challenge for the understanding of subsurface processes. We have developed a novel method for ground-penetrating-radar (GPR) reflection monitoring of subsurface flow processes under unsaturated conditions and applied it to a hydrological infiltration experiment performed across a periglacial slope deposit in northwest Luxembourg. Our approach relies on a spatial and temporal quasicontinuous data recording and processing, followed by an attribute analysis based on analyzing differences between individual time steps. The results demonstrate the ability of time-lapse GPR monitoring to visualize the spatial and temporal dynamics of preferential flow processes with a spatial resolution in the order of a few decimeters and temporal resolution in the order of a few minutes. We observe excellent agreement with water table information originating from different boreholes. This demonstrates the potential of surface-based GPR reflection monitoring to observe the spatiotemporal dynamics of water movements in the subsurface. It provides valuable, and so far not accessible, information for example in the field of hydrology and pedology that allows studying the actual subsurface processes rather than deducing them from point information.

Quantitative imaging for civil engineering by joint full waveform inversion of surface-based GPR and shallow seismic reflection data

Full waveform inversion (FWI) of ground penetrating radar (GPR) and seismic data is a promising imaging tool for the detailed characterization of latent cavities and tunnels. In this study, surface-based GPR FWI and seismic FWI are combined to quantitatively image the geometry and geophysical attributes of the subsurface. These two FWI methods are integrated into a joint procedure and linked by cross-gradient functions which can provide structural constraints between the parameters. The elastic and electromagnetic properties of the medium (P-wave velocity, permittivity, and conductivity) are linked and simultaneously inverted in the joint FWI approach. This joint FWI approach is tested using three numerical examples. In the first two examples, target objects with various shapes are contained in both homogeneous and inhomogeneous media. In the third example, two sediment-filled cavities and a fractured area are analyzed in the calcarenite layers. The results show that the joint FWI approach can integrate the common structure of models and correct the model structures that are inaccurately reconstructed in the individual FWI methods, providing a comprehensive structural interpretation of the study areas.

Surface-based GPR underestimates below-stump root biomass

While lateral root mass is readily detectable with ground penetrating radar (GPR), the roots beneath a tree (below-stump) and overlapping lateral roots near large trees are problematic for surface-based antennas operated in reflection mode. We sought to determine if tree size (DBH) effects GPR root detection proximal to longleaf pine (Pinus palustris Mill) and if corrections for could be applied to stand-level estimates of root mass. GPR (1500 MHz) was used to estimate coarse root mass proximal to 33 longleaf pine trees and compared to the amount of biomass excavated from pits proportional in area to tree basal diameter. Lateral roots were excavated to a depth of 1 m and taproots were excavated in their entirety. GPR underestimated longleaf pine below-stump mass and the magnitude of the underestimation increased with tree DBH. Non-linear regressions between GPR estimated root mass/excavated root mass and tree diameter at breast height (DBH) were highly significant for both below-stump (lateral + taproot) root mass (p < 0.0001, R2 0.77) and lateral coarse root mass (p < 0.0001, R2 0.65). GPR underestimates root mass proximal to trees, and this needs to be accounted for to accurately estimate stand-level belowground biomass.

Resolving precipitation induced water content profiles by inversion of dispersive GPR data: A numerical study

Summary Surface-based ground-penetrating radar (GPR) measurements have significant potential for monitoring dynamic hydrologic processes at multiple scales in time and space. At early times during infiltration into a soil, the zone above the wetting front may act as a low-velocity waveguide that traps GPR waves, thereby causing dispersion and making interpretation of the data using standard methods difficult. In this work, we show that the dispersion is dependent upon the distribution of water within the waveguide, which is controlled by soil hydrologic properties. Simulations of infiltration were performed by varying the n-parameter of the Mualem–van Genuchten equation using HYDRUS-1D; the associated GPR data were simulated to evaluate the influence of dispersion. We observed a notable decrease in wave dispersion as the sharpness of the wetting front profile decreased. Given the sensitivity of the dispersion effect to the wetting front profile, we also evaluated whether the water content distribution can be determined through inversion of the dispersive GPR data. We found that a global grid search combined with the simplex algorithm was able to estimate the average water content when the wetted zone is divided into 2 layers. This approach was incapable, however, of representing the gradational nature of the water content distribution behind the wetting front. In contrast, the shuffled complex evolution algorithm was able to constrain a piece-wise linear function to closely match the shallow gradational water content profile. In both the layered and piece-wise linear case, the sensitivity of the dispersive data dropped sharply below the wetting front, which in this case was around 20 cm, i.e., twice the average wavelength, for a 900 MHz GPR survey. This study demonstrates that dispersive GPR data has significant potential for capturing the early-time dynamics of infiltration that cannot be obtained with standard GPR analysis approaches.

Quantitative high-resolution observations of soil water dynamics in a complicated architecture using time-lapse ground-penetrating radar

Abstract. High-resolution time-lapse ground-penetrating radar (GPR) observations of advancing and retreating water tables can yield a wealth of information about near-surface water content dynamics. In this study, we present and analyze a series of imbibition, drainage and infiltration experiments that have been carried out at our artificial ASSESS test site and observed with surface-based GPR. The test site features a complicated but known subsurface architecture constructed with three different kinds of sand. It allows the study of soil water dynamics with GPR under a wide range of different conditions. Here, we assess in particular (i) the feasibility of monitoring the dynamic shape of the capillary fringe reflection and (ii) the relative precision of monitoring soil water dynamics averaged over the whole vertical extent by evaluating the bottom reflection. The phenomenology of the GPR response of a dynamically changing capillary fringe is developed from a soil physical point of view. We then explain experimentally observed phenomena based on numerical simulations of both the water content dynamics and the expected GPR response.

Automated high-resolution GPR data collection for monitoring dynamic hydrologic processes in two and three dimensions

An automated GPR data-collection system was constructed to monitor dynamic hydrologic processes with high spatiotemporal resolution. The design of the system allows for fast acquisition of constant-offset profiles (COP) and common-midpoint surveys (CMP) to monitor unsaturated flow at multiple locations. A fast and accurate motion-control system enables the high-resolution monitoring of spatial patterns of reflectors through time while preserving important information about the normal moveout of reflectors for velocity analysis. In two experiments, infiltration was monitored using time-lapse ground-penetrating radar (GPR) measurements in two and three dimensions. The data from both experiments provide substantial qualitative insight about the dynamics of hydrologic events and a path forward for quantitative analysis of surface-based GPR monitoring data. Analysis shows (1) the advantages of collecting high-resolution time-lapse data, (2) the complexities of patterns associated with the wetting of the soil, and (3) evidence of nonuniform propagation of a wetting front through the soil column.

Vertical radar profiling: Combined analysis of traveltimes, amplitudes, and reflections

Vertical radar profiling (VRP) is a single-borehole geophysical technique, in which the receiver antenna is located within a borehole and the transmitter antenna is placed at one or various offsets from the borehole. Today, VRP surveying is primarily used to derive 1D velocity models by inverting the arrival times of direct waves. Using field data collected at a well-constrained test site in Germany, we evaluated a VRP workflow relying on the analysis of direct-arrival traveltimes and amplitudes as well as on imaging reflection events. To invert our VRP traveltime data, we used a global inversion strategy resulting in an ensemble of acceptable velocity models, and thus, it allowed us to appraise uncertainty issues in the estimated velocities as well as in porosity models derived via petrophysical translations. In addition to traveltime inversion, the analysis of direct-wave amplitudes and reflection events provided further valuable information regarding subsurface properties and architecture. The used VRP amplitude preprocessing and inversion procedures were adapted from ray-based crosshole ground-penetrating radar (GPR) attenuation tomography and resulted in an attenuation model, which can be used to estimate variations in electrical resistivity. Our VRP reflection imaging approach relied on corridor stacking, which is a well-established processing sequence in vertical seismic profiling. The resulting reflection image outlines bounding layers and can be directly compared to surface-based GPR reflection profiling. Our results of the combined analysis of VRP, traveltimes, amplitudes, and reflections were consistent with independent core and borehole logs as well as GPR reflection profiles, which enabled us to derive a detailed hydrostratigraphic model as needed, for example, to understand and model groundwater flow and transport.

High‐resolution water content estimation from surface‐based ground‐penetrating radar reflection data by impedance inversion

Mapping hydrological parameter distributions in high resolution is essential to understand and simulate groundwater flow and contaminant transport. Of particular interest is surface‐based ground‐penetrating radar (GPR) reflection imaging in electrically resistive sediments because of the expected close link between the subsurface water content and the dielectric permittivity, which controls GPR wave velocity and reflectivity. Conventional tools like common midpoint (CMP) velocity analysis provide physical parameter models of limited resolution only. We present a novel reflection amplitude inversion workflow for surface‐based GPR data capable of resolving the subsurface dielectric permittivity and related water content distribution with markedly improved resolution. Our scheme is an adaptation of a seismic reflection impedance inversion scheme to surface‐based GPR data. Key is relative‐amplitude‐preserving data preconditioning including GPR deconvolution, which results in traces with the source‐wavelet distortions and propagation effects largely removed. The subsequent inversion for the underlying dielectric permittivity and water content structure is constrained by in situ dielectric permittivity data obtained by direct‐push logging. After demonstrating the potential of our novel scheme on a realistic synthetic data set, we apply it to two 2‐D 100 MHz GPR profiles acquired over a shallow sedimentary aquifer resulting in water content images of the shallow (3–7 m depth) saturated zone having decimeter resolution.

Constraining 3-D electrical resistance tomography with GPR reflection data for improved aquifer characterization

Abstract Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240 × 40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with the underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (> 300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~ 200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.

Subsurface Water-filled Fracture Detection by Borehole Radar: A Case History

Operating radar equipment in boreholes offers the possibility of greater depth penetration and higher resolution than is achievable with surface-based ground-penetrating radar systems. We have acquired single-hole radar reflection data in a series of vertical boreholes situated within a granite body west of Beijing, China. Geological logs demonstrate that numerous fractures intersect the boreholes. After processing the single-hole reflection data, linear reflections from many of these fractures are observed. In addition, a large number of linear features that do not intersect the boreholes are interpreted as fracture reflections; they can be traced to distances of [Formula: see text] from the boreholes. Information from single-hole radar reflection data allows the minimum lengths, dips and distances of planar fractures from boreholes to be estimated, but not the azimuths. To determine the azimuths, information from two or more boreholes is required. For our survey site, this was achieved for three of the fractures.