Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.

Architectural characterization of a delta-front reservoir analogue combining ground penetrating radar and electrical resistivity tomography: Roda Sandstone (Lower Eocene, Graus-Tremp basin, Spain)

Three-dimensional reconstruction of reservoir analogues can be improved combining data from different geophysical methods. Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT) data are valuable tools, since they provide subsurface information from internal architecture and facies distribution of sedimentary rock bodies, enabling the upgrading of depositional models and heterogeneity reconstruction. The Lower Eocene Roda Sandstone is a well-known deltaic complex widely studied as a reservoir analogue that displays a series of sandstone wedges with a general NE to SW progradational trend. To provide a better understanding of internal heterogeneity of a 10m-thick progradational delta-front sandstone unit, 3D GPR data were acquired. In addition, common midpoints (CMP) to measure the sandstone subsoil velocity, test profiles with different frequency antennas (25, 50 and 100MHz) and topographic data for subsequent correction in the geophysical data were also obtained. Three ERT profiles were also acquired to further constrain GPR analysis. These geophysical results illustrate the geometry of reservoir analogue heterogeneities both depositional and diagenetic in nature, improving and complementing previous outcrop-derived data. GPR interpretation using radar stratigraphy principles and attributes analysis provided: 1)tridimensional geometry of major stratigraphic surfaces that define four units in the GPR Prism, 2) image the internal architecture of the units and their statistical study of azimuth and dips, useful for a quick determination of paleocurrent directions. These results were used to define the depositional architecture of the progradational sandbody that shows an arrangement in very-high-frequency sequences characterized by clockwise paleocurrent variations and decrease of the sedimentary flow, similar to those observed at a greater scale in the same system. This high-frequency sequential arrangement has been attributed to the autocyclic dynamics of a supply-dominated delta- front where fluvial and tidal currents are in competition. The resistivity models enhanced the viewing of reservoir quality associated with cement distribution caused by depositional and early diagenetic processes related to the development of transgressive and regressive systems tracts in igh-frequency sequences.

Use of ground-penetrating radar to map subsurface features at the Lapa do Santo archaeological site (Brazil)

In this article we present some results of ground-penetrating radar (GPR) studies carried out at the Lapa do Santo archaeological site. This cave is within the Lagoa Santa karstic region, Minas Gerais State, Brazil. Results from 44 GPR profiles obtained with 400 MHz shielded antennas indicated anomalous hyperbolic reflections and areas with high sub-horizontal reflection amplitude suggesting archaeological and geological potential targets, respectively. These results were encouraging and were used to guide excavations at this site. Excavation of test units (metre by metre) allowed identifying an anthropogenic feature, e.g., a fire hearth structure and natural features, such as a stalagmite and top of bedrock. Results also indicated the importance of the GPR survey as a tool for orienting archaeological researches, increasing the probability of finding archaeological interest targets in an excavation program in an area of environmental protection.

The use of Ground Penetrating Radar and alternative geophysical techniques for assessing embankments and dykes safety

The research is part of a survey for the detection of the hydraulic and geotechnical conditions of river embankments funded by the Reno River Basin Regional Technical Service of the Region Emilia-Romagna. The hydraulic safety of the Reno River, one of the main rivers in North-Eastern Italy, is indeed of primary importance to the Emilia-Romagna regional administration. The large longitudinal extent of the banks (several hundreds of kilometres) has placed great interest in non-destructive geophysical methods, which, compared to other methods such as drilling, allow for the faster and often less expensive acquisition of high-resolution data. The present work aims to experience the Ground Penetrating Radar (GPR) for the detection of local non-homogeneities (mainly stratigraphic contacts, cavities and conduits) inside the Reno River and its tributaries embankments, taking into account supplementary data collected with traditional destructive tests (boreholes, cone penetration tests etc.). A comparison with non-destructive methodologies likewise electric resistivity tomography (ERT), Multi-channels Analysis of Surface Waves (MASW), FDEM induction, was also carried out in order to verify the usability of GPR and to provide integration of various geophysical methods in the process of regular maintenance and check of the embankments condition. The first part of this thesis is dedicated to the explanation of the state of art concerning the geographic, geomorphologic and geotechnical characteristics of Reno River and its tributaries embankments, as well as the description of some geophysical applications provided on embankments belonging to European and North-American Rivers, which were used as bibliographic basis for this thesis realisation. The second part is an overview of the geophysical methods that were employed for this research, (with a particular attention to the GPR), reporting also their theoretical basis and a deepening of some techniques of the geophysical data analysis and representation, when applied to river embankments. The successive chapters, following the main scope of this research that is to highlight advantages and drawbacks in the use of Ground Penetrating Radar applied to Reno River and its tributaries embankments, show the results obtained analyzing different cases that could yield the formation of weakness zones, which successively lead to the embankment failure. As advantages, a considerable velocity of acquisition and a spatial resolution of the obtained data, incomparable with respect to other methodologies, were recorded. With regard to the drawbacks, some factors, related to the attenuation losses of wave propagation, due to different content in clay, silt, and sand, as well as surface effects have significantly limited the correlation between GPR profiles and geotechnical information and therefore compromised the embankment safety assessment. Recapitulating, the Ground Penetrating Radar could represent a suitable tool for checking up river dike conditions, but its use has significantly limited by geometric and geotechnical characteristics of the Reno River and its tributaries levees. As a matter of facts, only the shallower part of the embankment was investigate, achieving also information just related to changes in electrical properties, without any numerical measurement. Furthermore, GPR application is ineffective for a preliminary assessment of embankment safety conditions, while for detailed campaigns at shallow depth, which aims to achieve immediate results with optimal precision, its usage is totally recommended. The cases where multidisciplinary approach was tested, reveal an optimal interconnection of the various geophysical methodologies employed, producing qualitative results concerning the preliminary phase (FDEM), assuring quantitative and high confidential description of the subsoil (ERT) and finally, providing fast and highly detailed analysis (GPR). Trying to furnish some recommendations for future researches, the simultaneous exploitation of many geophysical devices to assess safety conditions of river embankments is absolutely suggested, especially to face reliable flood event, when the entire extension of the embankments themselves must be investigated.

Ground penetrating radar early-time technique for soil electromagnetic parameters estimation

In recent years, thanks to the technological advances, electromagnetic methods for non-invasive shallow subsurface characterization have been increasingly used in many areas of environmental and geoscience applications. Among all the geophysical electromagnetic methods, the Ground Penetrating Radar (GPR) has received unprecedented attention over the last few decades due to its capability to obtain, spatially and temporally, high-resolution electromagnetic parameter information thanks to its versatility, its handling, its non-invasive nature, its high resolving power, and its fast implementation. The main focus of this thesis is to perform a dielectric site characterization in an efficient and accurate way studying in-depth a physical phenomenon behind a recent developed GPR approach, the so-called early-time technique, which infers the electrical properties of the soil in the proximity of the antennas. In particular, the early-time approach is based on the amplitude analysis of the early-time portion of the GPR waveform using a fixed-offset ground-coupled antenna configuration where the separation between the transmitting and receiving antenna is on the order of the dominant pulse-wavelength. Amplitude information can be extracted from the early-time signal through complex trace analysis, computing the instantaneous-amplitude attributes over a selected time-duration of the early-time signal. Basically, if the acquired GPR signals are considered to represent the real part of a complex trace, and the imaginary part is the quadrature component obtained by applying a Hilbert transform to the GPR trace, the amplitude envelope is the absolute value of the resulting complex trace (also known as the instantaneous-amplitude). Analysing laboratory information, numerical simulations and natural field conditions, and summarising the overall results embodied in this thesis, it is possible to suggest the early-time GPR technique as an effective method to estimate physical properties of the soil in a fast and non-invasive way.

Detecting Improvised Explosive Devices Via Forward Looking Ground Penetrating Radar

Forward-looking ground penetrating radar shows promise for detection of improvised explosive devices in active war zones. Because of certain insurmountable physical limitations, post-processing algorithm development is the most popular research topic in this field. One such investigative avenue explores the worthiness of frequency analysis during data post-processing. Using the finite difference time domain numerical method, simulations are run to test both mine and clutter frequency response. Mines are found to respond strongest at low frequencies and cause periodic changes in ground penetrating radar frequency results. These results are called into question, however, when clutter, a phenomenon generally known to be random, is also found to cause periodic frequency effects. Possible causes, including simulation inaccuracy, are considered. Although the clutter models used are found to be inadequately random, specular reflections of differing periodicity are found to return from both the mine and the ground. The presence of these specular reflections offers a potential alternative method of determining a mine’s presence.

Geophysical Investigation of Subsurface Three-Dimensional Architecture of Icy Debris Fans With Ground Penetrating Radar in the Wrangell Mountains, Alaska

Icy debris fans have are newly-described landforms (Kochel and Trop, 2008 and 2012) as landforms developed immediately after deglaciation on Earth and similar features have been observed on Mars. Subsurface characteristics of Icy debris fans have not been previously investigated. Ground penetrating radar (GPR) was used to non-invasively investigate the subsurface characteristics of icy debris fans near McCarthy, Alaska, USA. The three fans investigated in Alaska are the East, West, and Middle fans (Kochel and Trop, 2008 and 2012) which below the Nabesna ice cap and on top of the McCarthy Creek Glacier. Icy debris fans in general are a largely unexplored suite of paraglacial landforms and processes in alpine regions. Recent field studies focused on direct observations and depositional processes. Their results showed that the fan's composition is primarily influenced by the type and frequency of depositional processes that supply the fan. Photographic studies show that the East Fan receives far more ice and snow avalanches whereas the Middle and West Fans receive fewer mass wasting events but more clastic debris is deposited on the Middle and West fan from rock falls and icy debris flows. GPR profiles and Wide-angle reflection and refraction (WARR) surveys consisting of both, common mid-point (CMP), and common shot-point (CSP) surveys investigated the subsurface geometry of the fans and the McCarthy Creek Glacier. All GPR surveys were collected in July of 2013 with 100MHz bi-static antennas. Four axial profiles and three cross-fan profiles were done on the West and Middle fans as well as the McCarthy Creek Glacier in order to investigate the relationship between the three features. GPR profiles yielded reflectors that were continuous for 10+ m and hyperbolic reflections in the subsurface. The depth to these reflections in the subsurface requires knowledge of the velocity of the subsurface. To find the velocity of the subsurface eight WARR surveys collected on the fans and on the McCarthy Creek glacier to provide information on variability of subsurface velocities. The profiles of the Middle and West fan have more reflections in their profiles compared to profiles done on the McCarthy Creek Glacier. Based on the WARR surveys, we interpret the lower energy return in the glacier to be caused by two reasons. 1) The increased attenuation due to wet ice versus drier ice and on the fan with GPR velocities >0.15m/ns. 2) Lack of interfaces in the glacier compared to those in the fans which are inferred to be produced by the alternating layers of stratified ice and lithic-rich layers. The GPR profiles on the West and Middle Fans show the shallow subsurface being dominated by lenticular reflections interpreted to be consistent with the shape of surficial deposits. The West Fan is distinguished from the Middle Fan by the nature of its reflections patterns and thicknesses of reflection packages that clearly shows the Middle fan with a greater thickness. The changes in subsurface reflections between the Middle and West Fans as well as the McCarthy Creek Glacier are thought to reflect the type and frequency of depositional processes and surrounding bedrock and talus slopes.

In Situ Monitoring of Free-Phase Gas Accumulation and Release in Peatlands Using Ground Penetrating Radar (GPR)

We tested a set of surface common mid-point (CMP) ground penetrating radar (GPR) surveys combined with elevation rods ( to monitor surface deformation) and gas flux measurements to investigate in-situ biogenic gas dynamics and ebullition events in a northern peatland ( raised bog). The main findings are: ( 1) changes in the two-way travel time from the surface to prominent reflectors allow estimation of average gas contents and evolution of free-phase gas (FPG); ( 2) peat surface deformation and gas flux measurements are strongly consistent with GPR estimated changes in FPG content over time; ( 3) rapid decreases in atmospheric pressure are associated with increased gas flux; and ( 4) single ebullition events can induce releases of methane much larger ( up to 192 g/m(2)) than fluxes reported by others. These results indicate that GPR is a useful tool for assessing the spatial distribution, temporal variation, and volume of biogenic gas deposits in peatlands.

Spatial Variability in Biogenic Gas Accumulations in Peat Soils Is Revealed By Ground Penetrating Radar (GPR)

We performed surface and borehole ground penetrating radar (GPR) tests, together with moisture probe measurements and direct gas sampling to detect areas of biogenic gas accumulation in a northern peatland. The main findings are: (1) shadow zones (signal scattering) observed in surface GPR correlate with areas of elevated CH4 and CO2 concentration; (2) high velocities in zero offset profiles and lower water content inferred from moisture probes correlate with surface GPR shadow zones; (3) zero offset profiles depict depth variable gas accumulation from 0-10% by volume; (4) strong reflectors may represent confining layers restricting upward gas migration. Our results have implications for defining the spatial distribution, volume and movement of biogenic gas in peatlands at multiple scales.

Radar investigations of Apollinaris Mons on Mars: Exploring the origin of the fan deposits

Apollinaris Mons is an isolated volcano on Mars straddling the boundary between the southern highlands and the northern plains. One of its most distinctive features is its massive fan-shaped deposit that extends from a breach on its summit to distances of more than 150 km and drapes its entire southern flank. The composition and formation mechanism of these deposits remains controversial. We investigate the radar properties of the fan deposits (FD) of Apollinaris Mons using low-frequency sounding radar data in combination with high-resolution images and crater-size frequency analysis to constrain their inner shape and bulk composition. Our analysis indicates that the FD attains an irregular thickness and is gradually thinner towards their lateral margins. The crater-size frequency analysis shows that they may have undergone repeated resurfacing, which is suggestive of long-term evolution. Our analysis of Shallow Radar (SHARAD) radargrams traversing different sections of the FD reveals multiple and different subsurface interfaces among the radargrams crossing the thinnest part, which suggests a layered and complex inner shape. Our estimates for the bulk real part of the dielectric constant of the FD ranges from 3 to 5, which is consistent with an icy-silicate mixture or pyroclastic composition. Therefore, we conclude that lahars or pyroclastic flows are the most likely mechanism that created the FD, yet we cannot rule out additional contributions from lava flows. A combination of multiple processes is also possible since the deposits appear to have been modified by fluvial processes at a later stage of their formation.