Surface Ground Penetrating Radar Research Articles

Application of Subset Full-Waveform Inversion to Multioffset Surface Ground-Penetrating Radar Data

Application of Subset Full-Waveform Inversion to Multioffset Surface Ground-Penetrating Radar Data

Review of Ground Penetrating Radar Applications for Water Dynamics Studies in Unsaturated Zone

For water dynamics investigation in unsaturated (vadose) zones, ground penetrating radar is a popular hydro-geophysical method because it is non-invasive for soil, has high resolution and the results have a direct link with water content. Soil water content and soil hydraulic properties are two key factors for describing the water dynamics in vadose zones. There has been tremendous progress in soil water content and soil hydraulic properties estimation with ground penetrating radar. The purpose of this paper is to provide an overview of the application of ground penetrating radar for soil water dynamics studies. This paper first summarizes various methods for the determination of soil water content. including traditional methods in the surveys of surface ground penetrating radar, borehole ground penetrating radar, and off-ground ground penetrating radar, as well as relatively new methods, such as full waveform inversion, the average envelope amplitude method, and the frequency shift method. This paper further provides a review for estimating soil hydraulic properties with GPR according to the types of ground penetrating radar data. We hope that this review can provide a reference for the application of ground penetrating radar in soil water dynamics studies in the future.

Identification of buried historical mineshaft using ground-penetrating radar

Buried historical mineshafts constitute significant sinkhole hazards and sources of geotechnical risks. The locations of buried mineshafts on old mining maps are prone to inaccuracies. An efficient way to limit these inaccuracies is to use geophysical methods. However, their effectiveness strongly depends on the local survey conditions.The aim of the study was to identify the buried 14th century Wielki mineshaft in Bochnia Salt Mine, a UNESCO Worldwide Heritage Site, Poland. The need to locate this shaft was due to the following reasons. Firstly, the Bochnia Salt Mine had difficulties with water flowing into the workings from the supposed location of the Wielki mineshaft. The shaft area was one of the most waterlogged in the salt mine's terrain. An urgent task was to determine the shaft's exact location to implement a waterproof barrier to protect the historic mine workings. Secondly, the modernisation of the market square in Bochnia required clarification of the location of the Wielki mineshaft and checking as to whether it posed a sinkhole hazard to the terrain surface. Another important reason for locating the medieval shaft was its archaeological value.A ground surface and borehole ground-penetrating radar (GPR) study and geological-engineering recognition using boreholes were employed to indicate the location of the Wielki shaft. The anthropogenic layer of heterogeneous material and the urbanised site resulted in challenges with regard to detecting evident georadar anomalies. Several dozen intersecting GPR profiles were gathered in the research area using 250 MHz and 100 MHz antennae. The initial surface GPR results, indicative of the GPR anomaly, were consolidated by cored borehole investigation and a borehole GPR study. In particular, the presence of a void, loosened zones, and wooden fragments of the shaft lining was confirmed in the inclined control borehole. Following our classification of typical causes of the occurrence of sinkholes over buried mineshafts, the consolidation of backfill material in a mineshaft was determined as the cause of the sinkhole originated above the Wielki shaft (type B sinkhole). Destruction of the shaft's wooden cap may result in the occurrence of a sinkhole on the ground surface. Therefore, the void was backfilled and a hydraulic barrier was consequently made to stop the seepage of water into the Bochnia Salt Mine's historic workings.This case study demonstrates that surface and borehole GPR surveying integrated with borehole data offers the possibility of obtaining the high-resolution data needed to accurately locate an old buried mineshaft in an urbanised area. The GPR anomaly was indirectly detected on the basis of the deflection of the boundary between the natural ground and an anthropogenic layer with sufficient contrast in their electromagnetic properties. This approach can be useful wherever the georadar penetration range is limited by high attenuation and dispersiveness of the anthropogenic layer.

Estimating infiltration front depth using time-lapse multioffset gathers obtained from ground-penetrating-radar antenna array

Nondestructive and noninvasive visualization and quantification of dynamic subsurface hydrologic processes are needed. Using a ground-penetrating-radar (GPR) antenna array, time-lapse common-offset gather and common midpoint (CMP) data can be collected by fixing the antenna at a given location when scanning the subsurface. We have determined wetting front depths continuously during a field infiltration experiment by estimating electromagnetic (EM) wave velocities at given elapsed times using GPR antenna array data. A surface GPR antenna array system, consisting of 10 transmitters (Tx) and 11 receivers (Rx), that can scan each Tx-Rx combination in 10 s at a millisecond scale has been used to acquire all 110 Rx-Tx combinations in approximately 1.5 s. The field infiltration experiment has been conducted at an experimental field near the Tottori Sand Dunes in Japan. Using the estimated EM wave velocity from the CMP data, the depth to the wetting front is computed every minute. The estimated wetting front arrival time agrees with the time at which a sudden increase in the moisture sensor output is observed at a depth from 20 cm and below. We determine that time-lapsed CMP data collected with the GPR antenna array system could be used to estimate EM wave velocities continuously during the infiltration. The GPR antenna array is capable of accurate and quantitative tracking of the wetting front.

Evaluation of geophysical methods for characterizing industrial and municipal waste dumps

The main aim of this study was to assess the suitability of geophysical methods for investigating old waste dumps. Electrical Resistivity Tomography (ERT), Seismic Refraction Tomography (SRT), Multichannel Analysis of Surface Waves (MASW) and Ground Penetrating Radar (GPR) were the techniques used to characterize a waste dump in the town of Dąbrowa Górnicza (S Poland). GPR and SRT were the most difficult methods to use because of the dense vegetation, which attenuated the passage of electromagnetic and seismic waves to the ground. However, GPR did turn out to be an appropriate tool for characterizing variations in the surface cover. ERT, SRT and MASW clearly highlighted the transition between the waste deposit and the host sediments, and determined the approximate thickness of the waste deposits. With MASW, however, the waste type and the boundary between the waste layer and surface cover could not be delineated. In some places, the surface cover was identified using SRT. With both these methods, the problem with identification may be due to the small contrast in the S- and P-wave velocities through two kinds of waste (municipal and industrial), the thinness of the waste layer, and the considerable differentiation of the surface cover. The most accurate results regarding the waste deposits were obtained using ERT and different electrode spacings. ERT pinpointed the exact location of the stored waste, distinguished between the types of waste, and identified the soil cover. Data from shallow boreholes confirmed the interpretations of the methods.

Selection of geophysical methods based on matter-element analysis with analytic hierarchy process

The selection of the most appropriate method(s) for a specific engineering geophysical site characterisation can be challenging due to the complexity and potential ambiguity of the factors influencing the utility and cost-effectiveness of the methods. This paper, using matter-element analysis techniques with analytic hierarchy process (AHP), proposes a demonstration matter-element model for selecting the most appropriate geophysical methods, and establishes the general flow and calculation process for geophysical method selection. On the basis of the analysis of the factors influencing the utility and cost-effectiveness of geophysical methods, an evaluation criteria system is constructed. To demonstrate the application of this process, the proposed matter-element geophysical methods selection mode is applied to the selection of geophysical methods for an engineering geophysical site characterisation project that ultimately involved the acquisition of electrical resistivity tomography (ERT), multichannel analysis of surface wave data (MASW) and ground penetrating radar (GPR). The optimal results, based on the matter-element model, indicate the recommended geophysical methods are consistent with the field methods ultimately employed. The research proposes a new quantitative and objective computing approach to the selection of geophysical methods. The matter-element model can minimise subjective influences to a certain extent, and is conducive to enhancing the utility and cost-effectiveness of a geophysical survey.

Crosshole reflection imaging with ground-penetrating radar data: Applications in near-surface sedimentary settings

Crosshole ground-penetrating radar (GPR) is applied in areas that require a very detailed subsurface characterization. Analysis of such data typically relies on tomographic inversion approaches providing an image of subsurface parameters. We have developed an approach for processing the reflected energy in crosshole GPR data and applied it on GPR data acquired in different sedimentary settings. Our approach includes muting of the first arrivals, separating the up- and the downgoing wavefield components, and backpropagating the reflected energy by a generalized Kirchhoff migration scheme. We obtain a reflection image that contains information on the location of electromagnetic property contrasts, thus outlining subsurface architecture in the interborehole plane. In combination with velocity models derived from different tomographic approaches, these images allow for a more detailed interpretation of subsurface structures without the need to acquire additional field data. In particular, a combined interpretation of the reflection image and the tomographic velocity model improves the ability to locate layer boundaries and to distinguish different subsurface units. To support our interpretations of our field data examples, we compare our crosshole reflection results with independent information, including borehole logs and surface GPR data.

Which fractures are imaged with Ground Penetrating Radar? Results from an experiment in the Äspö Hardrock Laboratory, Sweden

Identifying fractures in the subsurface is crucial for many geomechanical and hydrogeological applications. Here, we assess the ability of the Ground Penetrating Radar (GPR) method to image open fractures with sub-mm apertures in the context of future deep disposal of radioactive waste. GPR experiments were conducted in a tunnel located 410 m below sea level within the Äspö Hard Rock Laboratory (Sweden) using 3-D surface-based acquisitions (3.4 m × 19 m) with 160 MHz, 450 MHz and 750 MHz antennas. The nature of 17 identified GPR reflections was analyzed by means of three new boreholes (BH1-BH3; 9–9.5 m deep). Out of 21 injection and outflow tests in packed-off 1-m sections, only five provided responses above the detection threshold with the maximum transmissivity reaching 7.0 × 10−10 m2/s. Most GPR reflections are situated in these permeable regions and their characteristics agree well with core and Optical Televiewer data. A 3-D statistical fracture model deduced from fracture traces on neighboring tunnel walls show that the GPR data mainly identify fractures with dips between 0 and 25°. Since the GPR data are mostly sensitive to open fractures, we deduce that the surface GPR method can identify 80% of open sub-horizontal fractures. We also find that the scaling of GPR fractures in the range of 1–10 m2 agrees well with the statistical model distribution indicating that fracture lengths are preserved by the GPR imaging (no measurement bias). Our results suggests that surface-GPR carries the resolution needed to identify the most permeable sub-horizontal fractures even in very low-permeability formations, thereby, suggesting that surface-GPR could play an important role in geotechnical workflows, for instance, for industrial-scale siting of waste canisters below tunnel floors in nuclear waste repositories.

Development of push moraines in deeply frozen sediment adjacent to a cold‐based glacier in the McMurdo Dry Valleys, Antarctica

AbstractPrevious studies of push moraines have concluded that structural architecture is controlled by a combination of stress exerted by flowing ice and the rheology of the deforming sediment. However, the relationship between deformation processes and the thermal/hydrologic conditions within the sediment remains ambiguous. Using a combination of surface mapping, logging exposures and ground penetrating radar we examine the relationships between surface morphology, material properties and development of deformation structures in deeply frozen sediment that has been deformed by Joyce Glacier, a cold‐based glacier in the McMurdo Dry Valleys. The structural architecture of the push moraines is characterized by a combination of brittle and ductile deformation structures that have produced a structurally complex pattern of thrust faults, low angle listric thrusts and recumbent folds that extend 400 m beyond the glacier margin. Deformation is ductile where the ice concentration exceeds c. 65% volume and predominantly brittle where the ice concentration is less than c. 65%. The change in rheology reflects transition in behaviour from the material having a predominantly frictional character when the ice is limited to pore spaces to a non‐frictional character in which strength is primarily determined by the cohesive strength of the ice. This work shows that glaciotectonic deformation can occur in deeply frozen permafrost where there is no liquid pore water. We conclude that the presence of liquid porewater is not a necessary condition for the development of glaciotectonic deformation or for the formation of push moraines. © 2019 John Wiley & Sons, Ltd.

The application of non-destructive methods in the diagnostics of the approach pavement at the bridges

The article presents the possibility of using non-destructive methods of road pavement diagnostics as an alternative to traditional means to assess the reasons for premature cracks adjacent to bridge objects. Two scanning methods were used: laser scanning to measure geometric surface deformation and ground penetrating radar (GPR) inspection to assess the road pavement condition. With the use of a laser scanner, an effective tool for road deformation assessment several approach pavement surfaces next to the bridges were scanned. As the result, a point cloud was obtained including spatial information about the pavement deformation. The data accuracy was about 3 mm, the deformations were presented in the form of deviation maps between the reference surface and the actual surface. Moreover characteristic pavement surface cross-sections were presented. The in situ measurements of the GPR method were performed and analysed in order to detect non-homogeneity in the density of structural layers of the pavement. Due to the analysis of the permittivity of individual layers, it was possible to detect non-homogeneity areas. The performed GPR measurements were verified by standard invasive tests carried out by drilling boreholes and taking cores from the pavement and testing the compaction and air voids content in asphalt layers. As a result of the measurements made by both methods significant differences in layer compacting factor values were diagnosed. The factor was much smaller in the area directly next to the bridgehead and much larger in the zone located a few meters away. The research showed the occurrence of both design and erection errors as well as those related to the maintenance of engineering structures.

Microwave Surface Reflection Method for Soil Moisture Determination Using frequency of 1.7-2.6 GHz

Soil moisture is one of the critical components to be investigated in civil, geological and agricultural works. This is because the parameter can affect the physical and electromagnetic characteristics of soil, such as density and permittivity and this can further restricts soil application. Unfortunately, finding a suitable non-destructive model for accurate soil moisture determination is challenging. In this article, the concept and development of soil moisture determination via ground penetrating radar (GPR) principle and surface reflection method is explained. The system is designed to be used with standard horn antenna with a frequency of 1.7-2.6 GHz along with vector network analyzer (VNA). The proposed system can measure soil moisture of three types of soil samples such as sand, loamy, and clay with high degree of accuracy. In this research, microwave surface reflection method is applied to analyze the effect of soil moisture with its electrical properties using our novel GPR principle. The result of the research is promising with high percentage of agreement with Topp theoretical value. The values are 31% to 61% for sand, 5% to 42% for clay, and 44% to 54% for loamy. For validation on the system, a new type of soil is used for measurement, and the result has an accuracy of 93%. By using the proposed developed models, soil moisture estimation can be easily determined with minimal data input through a novelty GPR surface reflection method.

Tomogram-based comparison of geostatistical models: Application to the Macrodispersion Experiment (MADE) site

We propose a new framework to compare alternative geostatistical descriptions of a given site. Multiple realizations of each of the considered geostatistical models and their corresponding tomograms (based on inversion of noise-contaminated simulated data) are used as a multivariate training image. The training image is scanned with a direct sampling algorithm to obtain conditional realizations of hydraulic conductivity that are not only in agreement with the geostatistical model, but also honor the spatially varying resolution of the site-specific tomogram. Model comparison is based on the quality of the simulated geophysical data from the ensemble of conditional realizations. The tomogram in this study is obtained by inversion of cross-hole ground-penetrating radar (GPR) first-arrival travel time data acquired at the MAcro-Dispersion Experiment (MADE) site in Mississippi (USA). Various heterogeneity descriptions ranging from multi-Gaussian fields to fields with complex multiple-point statistics inferred from outcrops are considered. Under the assumption that the relationship between porosity and hydraulic conductivity inferred from local measurements is valid, we find that conditioned multi-Gaussian realizations and derivatives thereof can explain the crosshole geophysical data. A training image based on an aquifer analog from Germany was found to be in better agreement with the geophysical data than the one based on the local outcrop, which appears to under-represent high hydraulic conductivity zones. These findings are only based on the information content in a single resolution-limited tomogram and extending the analysis to tracer or higher resolution surface GPR data might lead to different conclusions (e.g., that discrete facies boundaries are necessary).

Combination of H-Alpha Decomposition and Migration for Enhancing Subsurface Target Classification of GPR

Polarimetric technology has been one of the most important advances in microwave remote sensing during recent decades. H-alpha decomposition, which is a type of polarimetric analysis technique, has been common for terrain and land-use classification in polarimetric synthetic aperture radar. However, the technique has been less common in the ground penetrating radar (GPR) community. In this paper, we apply the H-alpha decomposition to analyze the surface GPR data to obtain polarimetric attributes for subsurface target classification. Also, by combining H-alpha decomposition and migration, we can obtain a subsurface H-alpha color-coded reconstructed target image, from which we can use both the polarimetric attributes and the geometrical features of the subsurface targets to enhance the ability of subsurface target classification of surface GPR. A 3-D full polarimetric GPR data set was acquired in a laboratory experiment, in which four targets, a scatterer with many branches, a ball, a plate, and a dihedral scatter, were buried in dry sand under flat ground surface, and used to test these techniques. As results, we obtained the subsurface H-alpha distribution and classified the subsurface targets. Also, we derived a subsurface H-alpha color-coded reconstructed target image and identified all four targets in the laboratory experiment.

Improved estimation of soil clay content by the fusion of remote hyperspectral and proximal geophysical sensing

Planning sustainable soil exploitation and land resource evaluation require up-to-date and accurate maps of soil properties. In that respect, geophysical techniques present particular interests given their non-invasiveness and their fast data acquisition capacity, which permit to characterize large areas with fine spatial and/or temporal resolutions. We investigated the relevancy of combining data from airborne hyperspectral (Hs), electromagnetic induction (EMI) and far-field ground-penetrating radar (GPR) for mapping soil properties, in particular soil clay content, at the field scale. Data from the three techniques were acquired at a test site in Mugello (Italy) characterized by relatively strong spatial variations of soil texture. Soil samples were collected for determining ground truth clay content. For the frequencies used in this study (200–650MHz), the GPR surface reflection is mainly determined by soil dielectric permittivity, itself primarily influenced by soil moisture. In contrast, EMI is mostly sensitive to soil electrical conductivity, which integrates several soil properties including in particular soil moisture and clay content. Taking advantage of the complementary information provided by the two instruments, the GPR and EMI data were combined and correlated to local ground-truth clay content data to provide high-resolution clay content maps over the entire field area. Besides, a relationship was also observed between Hs data and clay content measurements, which permitted to produce a Hs-derived clay content map. EMI–GPR and Hs maps showed close spatial patterns and a relatively high correlation was observed between both clay content estimates, as well as between clay content estimates and ground-truth clay content measurements. Moreover, data fusion allowed constraining the EMI–GPR and Hs information and reduced the uncertainty of mapped clay content estimates. These results demonstrated great promise for integrated, digital soil mapping applications.

Detecting Moisture Damage in Archaeology and Cultural Heritage: a Brief Introduction

Moisture damage is the most important issue in the preservation and integrity of cultural heritage. This paper discusses the ability of geophysical instruments to detect this problem. Non-destructive techniques (NDTs), such as Ground Penetrating Radar (GPR), use electromagnetic (EM) impulses to investigate archaeological sites and building structures that are affected by moisture and can be used to locate and estimate the extent of damage and to develop restoration plans before permanent damage occurs. The main objective of this paper is to introduce the capacity of surface GPR to rapidly and non-invasively estimate physical soil properties, develop novel processing strategies and provide valuable information about the investigated material in archaeological and cultural heritage sites. This new approach analyzes the amplitude attributes of the GPR pulse obtained from conventional single-offset surface-coupled profiling. To achieve the objective of this study, the technique is examined in two different experimental test settings to show that GPR analyses clearly highlight dampness as ringing anomalies with a very low EM signal amplitudes that are caused by high attenuation, poor antenna coupling, and temporal stretching. These indicators are important for diagnosing cultural heritage sites by allowing for the correct and precise visualization of radargrams and time-slices of the moisture anomalies.

Characterization of Cavities Using the GPR, LIDAR and GNSS Techniques

The study of the many types of natural and manmade cavities in different parts of the world is important to the fields of geology, geophysics, engineering, architectures, agriculture, heritages and landscape. Ground-penetrating radar (GPR) is a noninvasive geodetection and geolocation technique suitable for accurately determining buried structures. This technique requires knowing the propagation velocity of electromagnetic waves (EM velocity) in the medium. We propose a method for calibrating the EM velocity using the integration of laser imaging detection and ranging (LIDAR) and GPR techniques using the Global Navigation Satellite System (GNSS) as support for geolocation. Once the EM velocity is known and the GPR profiles have been properly processed and migrated, they will also show the hidden cavities and the old hidden structures from the cellar. In this article, we present a complete study of the joint use of the GPR, LIDAR and GNSS techniques in the characterization of cavities. We apply this methodology to study underground cavities in a group of wine cellars located in Atauta (Soria, Spain). The results serve to identify construction elements that form the cavity and group of cavities or cellars. The described methodology could be applied to other shallow underground structures with surface connection, where LIDAR and GPR profiles could be joined, as, for example, in archaeological cavities, sewerage systems, drainpipes, etc.

Lidar snow cover studies on glaciers in the Ötztal Alps (Austria): comparison with snow depths calculated from GPR measurements

Abstract. The storage of water within the seasonal snow cover is a substantial source of runoff in high mountain catchments. Information about the spatial distribution of snow accumulation is necessary for calibration and validation of hydro-meteorological models. Generally, only a small number of precipitation measurements deliver precipitation input for modelling in mountain areas. The spatial interpolation and extrapolation of measurements of precipitation is still difficult. Multi-temporal application of lidar techniques from aircraft, so-called airborne laser scanning (ALS), provides surface elevations changes even in inaccessible terrain. These ALS surface elevation changes can be used to derive changes in snow depths of the mountain snow cover for seasonal or subseasonal time periods. However, since glacier surfaces are not static over time, ablation, densification of snow, densification of firn and ice flow contribute to surface elevation changes. ALS-derived surface elevation changes were compared to snow depths derived from 35.4 km of ground penetrating radar (GPR) profiles on four glaciers. With this combination of two different data acquisitions, it is possible to evaluate the effect of the summation of these processes on ALS-derived snow depth maps in the high alpine region of the Ötztal Alps (Austria). A Landsat 5 Thematic Mapper image was used to distinguish between snow covered area and bare ice areas of the glaciers at the end of the ablation season. In typical accumulation areas, ALS surface elevation changes differ from snow depths calculated from GPR measurements by −0.4 m on average with a mean standard deviation of 0.34 m. Differences between ALS surface elevation changes and GPR derived snow depths are small along the profiles conducted in areas of bare ice. In these areas, the mean absolute difference of ALS surface elevation changes and GPR snow depths is 0.004 m with a standard deviation of 0.27 m. This study presents a systematic approach to analyze deviations from ALS generated snow depth maps to ground truth measurements on four different glaciers. We could show that ALS can be an important and reliable data source for the spatial distribution of snow depths for most parts of the here investigated glaciers. However, within accumulation areas, just utilizing ALS data may lead to systematic underestimation of total snow depth distribution.

Exploring Shallow Subsurface of Mars and Introducing the GPR Technique for Planetary Sciences(Exploring Mars beyond the Surface Features)

Several geophysical techniques are available for the investigation and characterization of shallow subsurface. Most commonly used are electrical methods, seismic refraction and reflection methods, multichannel analysis of surface waves, gravity, magnetic, electromagnetic induction and ground penetrating radar. Mars is the one of the terrestrial planet, under discussion from several years due to interesting facts obtained from surface features including the evidences for the presence of water. Surface features on Mars include plans, mountains, impact crater, lava flows and dust cover. The possibility of life on the Mars is a question mark which is yet to be defined. In this article GPR technique is proposed for planetary missions to explore the subsurface of Mars. Ground Penetrating Radar (GPR) is a recent technology which provides high resolution picture for near surface structures. Due to the small size, light weight, high resolution, and simple operational system the GPR is supposed to be a best tool for the exploration of near surface features and hydrology of Mars. Launching a planetary mission at Mars with GPR could provide the information about the subsurface structures. This information could be used to define many question marks about red planet.

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.

Estimating water content in an active landfill with the aid of GPR

Landfill gas (LFG) receives a great deal of attention due to both negative and positive environmental impacts, global warming and a green energy source, respectively. However, predicting the quantity of LFG generated at a given landfill, whether active or closed is difficult due to the heterogeneities present in waste, and the lack of accurate in situ waste parameters like water content. Accordingly, ground penetrating radar (GPR) was evaluated as a tool for estimating in situ water content. Due to the large degree of subsurface heterogeneity and the electrically conductive clay cap covering landfills, both of which affect the transmission of the electromagnetic pulses, there is much scepticism concerning the use of GPR to quantify in situ water content within a municipal landfill.Two landfills were studied. The first landfill was used to develop the measurement protocols, while the second landfill provided a means of confirming these protocols. GPR measurements were initially completed using the surface GPR approach, but the lack of success led to the use of borehole (BH) GPR. Both zero offset profiling (ZOP) and multiple offset gathers (MOG) modes were tried, with the results indicating that BH GPR using the ZOP mode is the most simple and efficient method to measure in situ water content. The best results were obtained at a separation distance of 2m, where higher the water content, smaller the effective separation distance. However, an increase in water content did appear to increase the accuracy of the GPR measurements. For the effective separation distance of 2m at both landfills, the difference between GPR and lab measured water contents were reasonable at 33.9% for the drier landfill and 18.1% for the wetter landfill. Infiltration experiments also showed the potential to measure small increases in water content.