Time-lapse Ground-penetrating Radar Research Articles

Monitoring the evolution and migration of a methane gas plume in an unconfined sandy aquifer using time-lapse GPR and ERT

Fugitive methane (CH4) leakage associated with conventional and unconventional petroleum development (e.g., shale gas) may pose significant risks to shallow groundwater. While the potential threat of stray (CH4) gas in aquifers has been acknowledged, few studies have examined the nature of its migration and fate in a shallow groundwater flow system. This study examines the geophysical responses observed from surface during a 72day field-scale simulated CH4 leak in an unconfined sandy aquifer at Canadian Forces Base Borden, Canada, to better understand the transient behaviour of fugitive CH4 gas in the subsurface. Time-lapse ground-penetrating radar (GPR) and electrical resistivity tomography (ERT) were used to monitor the distribution and migration of the gas-phase and assess any impacts to groundwater hydrochemistry. Geophysical measurements captured the transient formation of a CH4 gas plume emanating from the injector, which was accompanied by an increase in total dissolved gas pressure (PTDG). Subsequent reductions in PTDG were accompanied by reduced bulk resistivity around the injector along with an increase in the GPR reflectivity along horizontal bedding reflectors farther downgradient. Repeat temporal GPR reflection profiling identified three events with major peaks in reflectivity, interpreted to represent episodic lateral CH4 gas release events into the aquifer. Here, a gradual increase in PTDG near the injector caused a sudden lateral breakthrough of gas in the direction of groundwater flow, causing free-phase CH4 to migrate much farther than anticipated based on groundwater advection. CH4 accumulated along subtle permeability boundaries demarcated by grain-scale bedding within the aquifer characteristic of numerous Borden-aquifer multi-phase flow experiments. Diminishing reflectivity over a period of days to weeks suggests buoyancy-driven migration to the vadose zone and/or CH4 dissolution into groundwater. Lateral and vertical CH4 migration was primarily governed by subtle, yet measurable heterogeneity and anisotropy in the aquifer.

Form and function in hillslope hydrology: in situ imaging and characterization of flow-relevant structures

Abstract. The study deals with the identification and characterization of rapid subsurface flow structures through pedo- and geo-physical measurements and irrigation experiments at the point, plot and hillslope scale. Our investigation of flow-relevant structures and hydrological responses refers to the general interplay of form and function, respectively. To obtain a holistic picture of the subsurface, a large set of different laboratory, exploratory and experimental methods was used at the different scales. For exploration these methods included drilled soil core profiles, in situ measurements of infiltration capacity and saturated hydraulic conductivity, and laboratory analyses of soil water retention and saturated hydraulic conductivity. The irrigation experiments at the plot scale were monitored through a combination of dye tracer, salt tracer, soil moisture dynamics, and 3-D time-lapse ground penetrating radar (GPR) methods. At the hillslope scale the subsurface was explored by a 3-D GPR survey. A natural storm event and an irrigation experiment were monitored by a dense network of soil moisture observations and a cascade of 2-D time-lapse GPR trenches. We show that the shift between activated and non-activated state of the flow paths is needed to distinguish structures from overall heterogeneity. Pedo-physical analyses of point-scale samples are the basis for sub-scale structure inference. At the plot and hillslope scale 3-D and 2-D time-lapse GPR applications are successfully employed as non-invasive means to image subsurface response patterns and to identify flow-relevant paths. Tracer recovery and soil water responses from irrigation experiments deliver a consistent estimate of response velocities. The combined observation of form and function under active conditions provides the means to localize and characterize the structures (this study) and the hydrological processes (companion study Angermann et al., 2017, this issue).

Form and function in hillslope hydrology: characterization of subsurface flow based on response observations

Abstract. The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches, allowed us to link processes and put them in a larger context. Transfer to other scales beyond observational scale and generalizations, however, rely on the knowledge of structures (form) and remain speculative. The complementary approach with a methodological focus on form (i.e., structure exploration) is presented and discussed in the companion paper by Jackisch et al.(2017).

Ground penetrating radar for fracture mapping in underground hazardous waste disposal sites: A case study from an underground research tunnel, South Korea

Secure disposal or storage of nuclear waste within stable geologic environments hinges on the effectiveness of artificial and natural radiation barriers. Fractures in the bedrock are viewed as the most likely passage for the transport of radioactive waste away from a disposal site. We utilize ground penetrating radar (GPR) to map fractures in the tunnel walls of an underground research tunnel at the Korea Atomic Energy Research Institute (KAERI). GPR experiments within the KAERI Underground Research Tunnel (KURT) were carried out by using 200MHz, 500MHz, and 1000MHz antennas. By using the high-frequency antennas, we were able to identify small-scale fractures, which were previously unidentified during the tunnel excavation process. Then, through 3-D visualization of the grid survey data, we reconstructed the spatial distribution and interconnectivity of the multi-scale fractures within the wall. We found that a multi-frequency GPR approach provided more details of the complex fracture network, including deep structures. Furthermore, temporal changes in reflection polarity between the GPR surveys enabled us to infer the hydraulic characteristics of the discrete fracture network developed behind the surveyed wall. We hypothesized that the fractures exhibiting polarity change may be due to a combination of air-filled and mineralogical boundaries. Simulated GPR scans for the considered case were consistent with the observed GPR data. If our assumption is correct, the groundwater flow into these near-surface fractures may form the water-filled fractures along the existing air-filled ones and hence cause the changes in reflection polarity over the given time interval (i.e., 7days). Our results show that the GPR survey is an efficient tool to determine fractures at various scales. Time-lapse GPR data may be essential to characterize the hydraulic behavior of discrete fracture networks in underground disposal facilities.

Attribute-based analysis of time-lapse ground-penetrating radar data

Analysis of time-lapse ground-penetrating radar (GPR) data can provide information regarding subsurface hydrological processes, such as preferential flow. However, the analysis of time-lapse data is often limited by data quality; for example, for noisy input data, the interpretation of difference images is often difficult. Motivated by modern image-processing tools, we have developed two robust GPR attributes, which allow us to distinguish amplitude (contrast similarity) and time-shift (structural similarity) variations related to differences between individual time-lapse GPR data sets. We tested and evaluated our attributes using synthetic data of different complexity. Afterward, we applied them to a field data example, in which subsurface flow was induced by an artificial rainfall event. For all examples, we identified our structural similarity attribute to be a robust measure for highlighting time-lapse changes also in data with low signal-to-noise ratios. We determined that our new attribute-based workflow is a promising tool to analyze time-lapse GPR data, especially for imaging subsurface hydrological processes.

Cross-polarized GPR imaging of fracture flow channeling

Ground penetrating radar (GPR) can be used to image fractures and monitor fluid flow in the subsurface. Conventional GPR imaging uses single-polarization, co-polarized acquisition. We examine the use of cross-polarized GPR signals for imaging flow channeling in a discrete horizontal fracture. Numerical modeling (FDTD) demonstrates that when the fracture channel is oriented at an oblique angle to the survey line, depolarization of the GPR signal results in scattered energy along the cross-polarized components. When the channel is oriented parallel or orthogonal to the survey line, all scattered energy is captured by the co-polarized components and no signal is present in the cross-polarized orientation. Multipolarization, time-lapse 3D GPR field data were acquired at the Altona Flat Rock test site in New York State. The GPR surveys were conducted during background fresh fracture water conditions and during a natural gradient saline tracer test which was used to highlight flow channels along a sub-horizontal fracture. Amplitude analysis of the cross-polarized data reveals flow channeling that is in agreement with the co-polarized GPR images and with independent hydraulic tests. This investigation demonstrates that cross-polarized components of GPR signals can be used to enhance imaging of flow channels in fractured media.

Temporal Monitoring of the Soil Freeze-Thaw Cycles over a Snow-Covered Surface by Using Air-Launched Ground-Penetrating Radar

We tested an off-ground ground-penetrating radar (GPR) system at a fixed location over a bare agricultural field to monitor the soil freeze-thaw cycles over a snow-covered surface. The GPR system consisted of a monostatic horn antenna combined with a vector network analyzer, providing an ultra-wideband stepped-frequency continuous-wave radar. An antenna calibration experiment was performed to filter antenna and back scattered effects from the raw GPR data. Near the GPR setup, sensors were installed in the soil to monitor the dynamics of soil temperature and dielectric permittivity at different depths. The soil permittivity was retrieved via inversion of time domain GPR data focused on the surface reflection. Significant effects of soil dynamics were observed in the time-lapse GPR, temperature and dielectric permittivity measurements. In particular, five freeze and thaw events were clearly detectable, indicating that the GPR signals respond to the contrast between the dielectric permittivity of frozen and thawed soil. The GPR-derived permittivity was in good agreement with sensor observations. Overall, the off-ground nature of the GPR system permits non-invasive time-lapse observation of the soil freeze-thaw dynamics without disturbing the structure of the snow cover. The proposed method shows promise for the real-time mapping and monitoring of the shallow frozen layer at the field scale.

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.

High-resolution space–time quantification of soil moisture along a hillslope using joint analysis of ground penetrating radar and frequency domain reflectometry data

We combined ground-penetrating radar (GPR) and frequency domain reflectometry (FDR) to assess the space–time variability of soil moisture along a hillslope. Time-lapse GPR and FDR measurements were conducted weekly during the period 23/03–08/06/2011 along a cultivated hillslope in the Belgian loam belt. A full-wave GPR model, a soil dielectric mixing model and the Debye equation were combined to directly estimate soil moisture from GPR measurements. Measured GPR data were well reproduced by the full-wave GPR model, resulting in a relatively good agreement between the GPR and FDR-derived soil moisture. Subsequently, we merged the soil moisture obtained from both techniques in a data fusion framework and we investigated its spatial and temporal variability. The results indicate that there was a high correlation between the spatial variability of soil moisture and topography as well as between its temporal variability and rainfall. A temporal stability analysis showed that soil moisture at the footslope is higher and more stable than that at the summits and backslopes. The proposed approach appears to be promising for assessing soil moisture at the hillslope scale with a relatively high space–time resolution.

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.

Coupled hydrogeophysical inversion of time-lapse surface GPR data to estimate hydraulic properties of a layered subsurface

[1] A major challenge in vadose zone hydrology is to obtain accurate information on the temporal changes of the vertical soil water distribution and its feedback with the atmosphere and groundwater. A state of the art coupled hydrogeophysical inversion scheme is applied to evaluate soil hydraulic properties of a synthetic model and a field soil in southern Ontario based on time-lapse monitoring of soil dynamics with surface ground-penetrating radar (GPR). Film flow was included in the hydrological model to account for noncapillary water flow in a sandy medium during dry conditions. The synthetic study illustrated that GPR data contain sufficient information to accurately constrain soil hydraulic parameters within a coupled inversion framework and led to an accurate estimation of the soil hydraulic properties. When film flow was not accounted for within the inversion, an equally good fit could still be achieved. In this case, errors introduced by neglecting film flow were compensated by different hydraulic parameters. For the field data, the coupled inversion reduced the overall misfit compared to an uncalibrated model using hydraulic parameters obtained from laboratory data. Although the data fit improved significantly for water content in the deeper soil layers, accounting for film flow in the uppermost subsurface layer did not lead to a better fit of the GPR data. Further research is needed to describe the processes controlling water content in the dry range, in particular coupled heat and vapor transport. This study illustrates the suitability of surface GPR measurements combined with coupled inversion for near-surface characterization of soil hydraulic parameters.

Improving soil moisture profile reconstruction from ground-penetrating radar data: a maximum likelihood ensemble filter approach

Abstract. The vertical profile of shallow unsaturated zone soil moisture plays a key role in many hydro-meteorological and agricultural applications. We propose a closed-loop data assimilation procedure based on the maximum likelihood ensemble filter algorithm to update the vertical soil moisture profile from time-lapse ground-penetrating radar (GPR) data. A hydrodynamic model is used to propagate the system state in time and a radar electromagnetic model and petrophysical relationships to link the state variable with the observation data, which enables us to directly assimilate the GPR data. Instead of using the surface soil moisture only, the approach allows to use the information of the whole soil moisture profile for the assimilation. We validated our approach through a synthetic study. We constructed a synthetic soil column with a depth of 80 cm and analyzed the effects of the soil type on the data assimilation by considering 3 soil types, namely, loamy sand, silt and clay. The assimilation of GPR data was performed to solve the problem of unknown initial conditions. The numerical soil moisture profiles generated by the Hydrus-1D model were used by the GPR model to produce the "observed" GPR data. The results show that the soil moisture profile obtained by assimilating the GPR data is much better than that of an open-loop forecast. Compared to the loamy sand and silt, the updated soil moisture profile of the clay soil converges to the true state much more slowly. Decreasing the update interval from 60 down to 10 h only slightly improves the effectiveness of the GPR data assimilation for the loamy sand but significantly for the clay soil. The proposed approach appears to be promising to improve real-time prediction of the soil moisture profiles as well as to provide effective estimates of the unsaturated hydraulic properties at the field scale from time-lapse GPR measurements.

Examining the information content of time-lapse crosshole GPR data collected under different infiltration conditions to estimate unsaturated soil hydraulic properties

Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten–Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.

Field GPR monitoring of biostimulation in saturated porous media

Abstract Time-lapse borehole radar was used to monitor a sandy aquifer with dielectric property changes hypothesized to result from biomass growth and biodegradation of hydrocarbons. A blend of gasoline and ethanol was released below the water table at Canadian Forces Base (CFB) Borden, Canada, and a source of oxygen was used for cyclic stimulation of microbial activity over a period of two years. A dense grid of fourteen borehole ground-penetrating radar (GPR) pairs monitoring the bioactive region showed 200 MHz radar wave velocity changes of ± 4 % and signal attenuation changes of ± 25 % during three cycles of variable biostimulation. GPR signal changes correlated spatially and temporally to independent measurements of groundwater velocity flow changes and geochemical variations that occurred in response to microbial activity. Greater relative changes in radar wave velocity of propagation were observed in the region of enhanced bacterial stimulation, adjacent to the oxygen release wells, and along the path of the hydrocarbon plume where biomass growth was the greatest. In the zone of maximum biological activity, radar wave velocity decreased during two cycles of enhanced biostimulation, whereas it increased during the intervening period of low-level biostimulation. Monitored regions at the fringe of the plume and farther down-gradient of the oxygen source, where biostimulation was less or absent, exhibited GPR signal changes of lesser magnitude and consistency. Spatial and temporal variabilities of GPR observations indicate a high degree of transient heterogeneity within the relatively homogeneous Borden aquifer as a result of hydrocarbon biodegradation. Mechanisms postulated to cause the dielectric property changes observed in the radar data are biogenic gas formation, mineral dissolution/precipitation, and direct biomass growth. Geochemical data, in-situ groundwater velocity measurements and biomass quantification suggest biomass formation in the pore space as a plausible mechanism for the changes observed in the aquifer. Therefore, at 200 MHz frequency, biomass formation is observed to cause a decrease of EM wave velocity of propagation suggesting an increase of the bulk dielectric constant of the water saturated medium. This finding is in agreement with published laboratory-scale GPR monitoring of microbial growth in porous media. This study shows that time-lapse GPR along with supporting biogeochemical observations can be used for monitoring of biological activity at the field scale.

Efficient Electromagnetic Imaging of an Artificial Infiltration Process in the Vadose Zone Using Cross-Borehole Radar

Cross-borehole ground-penetrating radar (GPR) has been widely used to characterize the shallow subsurface and to monitor hydrogeologic processes. To investigate an infiltration process in the vadose zone, an artificial groundwater infiltration test was conducted in Nagaoka, Japan. Time-lapse cross-borehole GPR data were collected using zero-offset profiling (ZOP) mode. The infiltration process was observed as a variation of GPR traveltimes, which can be transformed into a dielectric constant, and further converted to volumetric water content. A standard ZOP analysis, for which all first arrivals are assumed to be direct waves, results in an underestimation of the dielectric constant because of the existence of critically refracted waves. This letter presents an efficient algorithm using the maximum first-cycle amplitude to approximately determine the traveltime of direct arrival, deriving a dielectric constant model more accurately than the standard ZOP analysis from ZOP data. Tests on synthetic and real field data show that the proposed approach is effective in building accurate water content profile without iterative calculations as in the standard ZOP analysis.

Time-lapse borehole radar monitoring of an infiltration experiment in the vadose zone

Ground-penetrating radar (GPR) is an effective tool for imaging the spatial distribution of water content. An artificial groundwater recharge test was conducted in Nagaoka City in Japan, and time-lapse crosshole GPR data were collected to monitor the infiltration process in the vadose zone. Since electromagnetic wave velocities in the vadose zone are largely controlled by variations in water content, an increase in traveltime is interpreted as an increase in saturation. In the test zone, the infiltrated water penetrated downward with an average velocity of about 2.7 m/h. A finite-difference time-domain method using two-dimensional cylindrical coordinates is applied to simulate radargrams associated with the advancing wetting front and to quantify the effects of critical refraction. Standard zero-offset profiling for which all first-arrivals are assumed to be direct waves results in an underestimation of water content in the transition zone above the wetting front. As a result, correct velocity analysis requires identification of first-arriving critically refracted waves from the traveltime profile to accurately determine a water content profile.