ANNOUNCEMENT
The Foreign Branch office of PARSAN becomes Operational in Singapore, Kingdom of Saudi Arabia and Bahrain!

Blog

Ground Penetrating Radar- Principle & Applications in Vadose Zone study




Ground Penetrating Radar- Principle & Applications in Vadose Zone study

Dr. Sanjay Rana

Chief Geophysicist & Director

PARSAN Overseas Private Limited

sanjay@parsan.biz


Abstract :


Soil erosion, groundwater depletion, salinization, and pollution have been recognized for decades as major threats to ecosystems and human health. The top part of the earth between the surface and the water table is called the vadose zone. The vadose zone mediates many of the processes that govern water resources and quality, such as the partition of precipitation into infiltration and runoff, groundwater recharge, contaminant transport, plant growth, evaporation, and energy exchanges between the earths surface and its atmosphere. It also determines soil organic carbon sequestration and carbon-cycle feedbacks, which could substantially affect climate change.


Information about near surface soil water content is a vital component for vadose zone, agricultural and ecological studies, as well as for climate models that require input about processes at the air-soil interface. This paper focuses on investigating the applicability of a surface geophysical method, ground penetrating radar (GPR), for use as a water content estimation tool.


The paper presents review of methods to measure soil water content with ground penetrating radar (GPR). Four methodologies have been discussed: soil water content determined from reflected wave velocity, soil water content determined from ground wave velocity, soil water content determined from transmitted wave velocity between boreholes, and soil water content determined from the surface reflection coefficient. Basic methodologies, capabilities and limitations of each approach have been discussed in the paper.

Introduction:

Water at the land surface is a vital resource for both human needs and natural ecosystems. Society's fresh water needs for agriculture, sanitation, municipal, and industrial supply are ever increasing. At the same time, natural hazards involving water, such as floods, droughts and landslides are major natural threats to society in many countries. The vadose zone, which may be defined as the transition zone between the atmosphere and groundwater reservoirs, is important for water resource management, because it regulates the water availability for vegetation, including crops, while at the same time provides a protective buffer zone between land surface and groundwater against solutes and pollutants.

Hydrologists, soil scientists, ecologists, meteorologists, and agronomists all study the space and time variability of water in the vadose zone, hereafter referred to as soil water content, over a range of scales and for a variety of reasons. At the regional to continental scale, the exchange of moisture and energy between soil, vegetation, and the atmosphere has an impact on near-surface atmospheric moisture and temperature, which in turn define the regional climate. For example, soil water content determines to a large extent the relative magnitudes of sensible and latent heat fluxes and therefore determines the diurnal evolution of the atmospheric boundary layer. Currently, there is a need to establish and quantify the contribution of soil water content  regulated land atmosphere coupling to regional climate anomalies, such as continental droughts and large-scale precipitation events.

The vadose zone mediates many of the processes that govern water resources and quality, such as the partition of precipitation into infiltration and runoff, groundwater recharge, contaminant  transport, plant growth, evaporation, and energy exchanges between the earths surface and its atmosphere. It also determines soil organic carbon sequestration and carbon-cycle feedbacks, which could substantially affect climate change.

The vadose zones inherent spatial variability and inaccessibility make direct observation of the important below-ground (termed subsurface) processes difficult. Conventional soil sampling is destructive, laborious, expensive, and may not be representative of the actual variability over space and time. In a societal context where the development of sustainable and optimal environmental management strategies has become a priority, there is a strong prerequisite for the development of noninvasive characterization and monitoring techniques of the vadose zone.

In particular, approaches integrating water movement, geological, and physical principles (called hydro-geophysics) applied at relevant scales are required to appraise dynamic belowground phenomena and to develop optimal sustainability, exploitation, and remediation strategies.

Among existing geophysical techniques, ground-penetrating radar (GPR) technology is of particular interest for providing high-resolution subsurface images and specifically addressing water-related questions. GPR is based on the transmission and reception of electromagnetic waves into the ground, whose propagation velocity and signal strength is determined by the soil electromagnetic properties and spatial distribution. As the electric permittivity of water overwhelms the permittivity of other soil components, the presence of water in the soil principally governs GPR wave propagation. Therefore, GPR-derived dielectric permittivity is usually used as surrogate measure for soil water content.

In the areas of unsaturated zone hydrology and water resources, GPR has been used to identify soil layering, locate water tables, follow wetting front movement, estimate soil water content, assist in subsurface hydraulic parameter identification, assess soil salinity, and support the monitoring of contaminants.

GPR has known a rapid development over the last decade, notes Sbastien Lambot, a well known vadose zone expert. Yet, several challenges must still be overcome before we can benefit from the full potential of GPR. In particular, more exact GPR modeling procedures together with the integration of other sources of information, such as other sensors or process knowledge, are required to maximize quantitative and qualitative information retrieval capabilities of GPR. Once this is achieved, GPR will be established as a powerful tool to support the understanding of the vadose zone hydrological processes and the development of optimal environmental and agricultural management strategies for our soil and water resources.

Clearly there is a need for soil water content measurements across a range of spatial scales. High-frequency electromagnetic techniques are the most promising category of soil water content sensors to fulfill this need because this category contains a range of techniques that measure the same soil water content proxy, namely dielectric permittivity, at different spatial scales. Remote sensing with either passive microwave radiometry or active radar instruments is the most promising technique for measuring soil water content variations over large regions. The passive instruments have low spatial resolution and can either be airborne with pixel sizes of thousands of square meters or satellite-borne with footprints in the order of tens of square kilometers. In the near future, passive satellite remote sensing will provide global coverage of critical hydrological data, including soil water content. Active radar instruments have smaller pixel sizes ranging from 1 to several 1000 m2. Although remote sensing will surely play an important role in many future hydrological studies, currently there is still a need to establish transfer functions between remote sensing and the more familiar in situ soil water content measurements. Additionally, because remote sensing approaches for estimating water content estimate water content in the uppermost 0.05 m of the soil and require that the vegetation cover is minimal, remote sensing methods are not applicable in all types of vadose zone studies.
 
A well-established in situ electromagnetic technique for soil water content investigations is time domain reflectometry (TDR), which was introduced in vadose zone hydrology in the early 1980s. Time domain reflectometry has developed into a reliable method for soil water content determination that can easily be automated. Furthermore, TDR can simultaneously measure dielectric permittivity and bulk soil conductivity, which allows the study of water and solute transport within the same soil volume. Although TDR is highly suited for monitoring the development of soil water content at one location with a high temporal resolution, the small measurement volume makes it sensitive to small-scale soil water content variation (e.g., macropores, air gaps due to TDR insertion) within this volume. Furthermore, assessment of spatial soil water content variation with TDR is labor intensive because TDR sensors need to be installed at each measurement location.

Clearly, there is a scale gap between remote sensing and TDR measurements of soil water content. At intermediate spatial scales, such as agricultural land and small catchments, reliance on sparse TDR measurements or coarse remote sensing measurements might not provide the accurate soil water content information required at these scales (e.g., crop management, precision farming). Therefore, there is a need for soil water content measurement techniques that can provide dense and accurate measurements at an intermediate scale. Since the early days of electromagnetic measurement techniques, ground penetrating radar (GPR) has been conceived as the natural intermediate-scale counterpart of TDR for soil water content measurements. Although the number of TDR applications has increased immensely in the past 20 yr, the number of GPR applications for measuring soil water content has only recently increased. Probably, the most important reason behind this delay is the more complicated behavior of the unguided waves used in GPR as compared with waves guided by a TDR sensor. Furthermore, recent improvements in GPR technology allow more accurate travel time measurements, which are needed for soil water content determination with GPR. In this review, we focus on investigating the applicability of a surface geophysical method, ground penetrating radar (GPR), for use as a water content estimation tool.


Vadose Zone :

The vadose zone, also termed the unsaturated zone, is the portion of Earth between the land surface and the phreatic zone or zone of saturation ("vadose" is Latin for "shallow"). It extends from the top of the ground surface to the water table. Water in the vadose zone has a pressure head less than atmospheric pressure, and is retained by a combination of adhesion (funiculary groundwater), and capillary action (capillary groundwater). If the vadose zone envelops soil, the water contained therein is termed soil moisture.

Movement of water within the vadose zone is studied within soil physics and hydrology, particularly hydrogeology, and is of importance to agriculture, contaminant transport, and flood control. The Richards equation is    often used to mathematically describe the flow of water, which is based partially on Darcy's law. Groundwater recharge, which is an important process that refills aquifers, generally occurs through the vadose zone from precipitation.

In speleology, cave passages formed in the vadose zone tend to be canyon-like in shape, as the water dissolves bedrock on the floor of the passage. Passages created in completely water-filled conditions are called phreatic passages and tend to be circular in cross-section.


                                                                                
 

This zone also includes the capillary fringe above the water table, the height of which will vary according to the grain size of the sediments. In coarse-grained mediums the fringe may be flat at the top and thin, whereas in finer grained material it will tend to be higher and may be very irregular along the upper surface. The vadose zone varies widely in thickness, from being absent to many hundreds of feet, depending upon several factors. These include the environment and the type of earth material present. Water within this interval, which is moving downward under the influence of gravity, is called vadose water, or gravitational water.

The vadose zone mediates many of the processes that govern water resources and quality, such as the partition of precipitation into infiltration and runoff, groundwater recharge, contaminant transport, plant growth, evaporation, and energy exchanges between the earths surface and its atmosphere. It also determines soil organic carbon sequestration and carbon-cycle feedbacks, which could substantially affect climate change.


                                                                                        


                                                                                              
Principles of Electromagnetic Methods

Electrical properties come from the interaction between electrical fields and charged particles, particularly the electron. Electrical conduction (transport) is the result of charge motion and results in energy dissipation (energy loss or conversion to heat). Electrical polarization (dielectric permittivity) is the result of charge separation over a distance, storing energy. Magnetic polarization (permeability or susceptibility) is the result of electron spin and motion in atomic orbits, and also results in energy loss and storage. Electrical and magnetic processes are also coupled, so accelerating electrons generate electromagnetic radiation, moving charges (currents) generate a magnetic field, and time varying magnetic fields cause charges to move. The velocity of electromagnetic wave propagation (speed of light) is the reciprocal of the square-root of the product of permittivity times permeability. The velocity in low loss, non-magnetic materials is (to a good approximation) the speed of light in vacuum divided by the square root of the relative dielectric permittivity (relative to that in free space).
 


                                                                                           v =  c
                                                                                                 ------
                                                                                                 √e'
  

Maxwell's equations (Maxwell, 1864, 1991) describe the propagation of an electromagnetic field. It is a coupled process, propagating as a three-dimensional, polarized, vector wave field. At low frequencies and high losses, the equations reduce to the diffusion equation and are called electromagnetic induction. At the high frequencies of radar, the energy storage in dielectric and magnetic polarization creates wave propagation. In the ideal, lossless case (vacuum), the electric and magnetic fields are in phase, orthogonal polarized vector fields, propagating at the speed of light. In real materials, they are out of phase, not completely polarized, propagating with a velocity lower than the speed of light in vacuum, scattered by changes in electric and magnetic properties, and with all of the preceding varying as functions of frequency.

The polarized vector field propagates in a straight line (neglecting relativistic gravitational effects) until it encounters a change in electrical or magnetic properties. At the change, the wave is scattered (reflected, refracted or diffracted) with amplitudes determined by the Fresnel reflection coefficient, angles determined by Snell's Law, and a polarization change described by the Stokes-Mueller matrices. Amplitudes are also varying with direction from the antenna (controlled by the antenna pattern), and with distance from the antenna (geometric spreading losses and material property dissipation losses) as described by the radar equation. Part of what determines whether or not significant scattering occurs is the spatial scale over which the change in properties occurs. This is both a detectability and a resolution issue. Resolution is determined by the spatial geometry of change versus the size of the wavelength of the propagating field.

The real part of the permittivity of water within the megahertz to gigahertz bandwidth is approximately 80, whereas the permittivity of air is 1 and of most other common soil constituents is about 3 to 10. This large contrast in permittivity explains the success of soil water content measurements with electromagnetic techniques working within this frequency bandwidth.

The most commonly used relationship between apparent permittivity,  , and volumetric soil water content,  (m3 m-3), was proposed by Topp et al. (1980):


                                                                       



This equation was determined empirically for mineral soils having various textures. It has an accuracy of 0.022 m3 m-3 determined in an independent validation on mineral soils. The term apparent is used because the permittivity used in this equation is determined from the measured electromagnetic propagation velocity in the soil.

It is important to realize that most available calibration equations between permittivity and water content were derived using TDR, which mainly operates in the frequency range from 500 to 1000 MHz. However, it has long been recognized that high clay content soils exhibit significant permittivity dispersion at low frequencies. Recently, West et al. (2003) presented frequency-dependent permittivity measurements of fine-grained sandstone samples containing up to 5% clay and soil samples containing Ottawa sand and varying amounts of montmorillonite clay. Their result showed that both the sandstone and soil samples showed significant frequency dispersion below 350 MHz. This implies that site-specific calibration may be required for those applications that require accurate water content measurements with lower antenna frequencies, such as the commonly used 100-MHz antenna. However, even when using published petrophysical relationships derived with TDR (such as above equation) with permittivity values obtained from GPR data, reasonable information about water content variation and spatial patterns can be obtained.
 

Principles of Ground Penetrating Radar :

Ground Penetrating Radar, also known as GPR, Georadar, Subsurface Interface Radar, Geoprobing Radar, is a totally non-destructive technique to produce a cross section profile of subsurface without any drilling, trenching or ground disturbances. Ground penetrating radar (GPR) profiles are used for evaluating the location and depth of buried objects and to investigate the presence and continuity of natural subsurface conditions and features, including determination of soil water content.

Basic Principle:

The GPR operates by transmitting electromagnetic impulses into the ground through transmitter antenna. The transmitted energy is reflected from various buried objects or distinct contacts between different earth materials, across which there is a contrast in dielectric constant. The antenna then receives the reflected waves and displays them in real time on screen. Data is also saved in appropriate memory for later processing and interpretation.

Ground penetrating radar waves can reach depths upto 60 meters in low conductivity materials such as dry sand or granite. Clays, shale
and other high conductivity materials may attenuate or absorb GPR signals, greatly decreasing the depth of penetration.

The depth of penetrating in also determined by the GPR antenna used. Antennas with low frequency obtain reflections from deeper depths but have low resolution. These low frequency antennas are used for investigating the geology of a site, such as for locating sinkholes or fractures, and to locate large, deep buried objects.

                                                                                                                                       


Antennas with higher frequencies (300 to 2000 MHz) obtain reflections from shallow depths (0 to 10 meters) and have a high resolution. These high frequency antennas are used to investigate surface soils and to locate small or large shallow buried objects, pipes, cables and rebar in concrete.


GPR can detect objects of any material, metallic or non-metallic.

Application Areas :

Geological and hydro-geological investigations including mapping of bedrock topography, water levels, solution features, glacial structures, soils and aggregates.

Engineering investigations to evaluate dams, sea walls, tunnels, pavements, roadbeds, railway embankments, piles, bridge decks, river scour, buildings and monuments.

Location and evaluation of buried structures including utilities, foundations, reinforcing bars, cavities, tombs, archaeological artifacts, and animal burrows.
 
Site investigations: location of buried engineering structures and underground storage tanks.

Subsurface mapping for cables, pipes and other buried structures prior to trench-less operations.

Advantages:

Rapid ground coverage- Antenna towed either by hand or from a vehicle.

High-resolution coverage of the survey area, detecting even small objects.

On-site interpretation possible due to instant graphic display.

Limitations:

Data acquisition may be slow over difficult terrain.

Depth of penetration is limited in materials with high electrical conductivities, clays.

Energy may be reflected and recorded from aboveground features, walls, canopies, unless antennae are well shielded.

Artifacts in the near surface (reinforcing bars, boulders, components of made ground) may scatter the transmitted energy and complicate the received signal and/or reduce depth of penetration.


Ground Penetrating Radar Theory :

GPR includes a radio transmitter and receiver, connected to an antenna coupled to the ground. The  transmitted signal penetrates a short distance into the ground and some of it reflects off any object with  different electrical properties than the soil.  Since plastic pipe and voids are different than soil, these  are some of the objects  which  reflect the signal.  Receiver  records  the  time and amplitude of reflected signal.


                                                                                


Multiple pulses sent  into  ground, and signal displayed on screen. The radar is moved along the ground and each new echo is plotted alongside the
previous ones.  When enough of these signals are plotted side-by-side, the operator can see a pattern which he can interpret as an object.
 
                                                                             


When the object is ahead of the radar, it takes more time for the echo to bounce back to the antenna. As it passes over, the time grows shorter, and then longer again as it goes past the object. This effect causes the image to take the shape of a curve called a "hyperbola". Experienced users recognize that a hyperbola is actually the image of a smaller object  (like a pipe) located at the center. Other patterns are produced by different structures. For example, a buried tank might have a flat image with curves down from either end.


                                                                              



It is not easy to determine the depth of an object without some knowledge about the dielectric constant in the local soil. The GPR can measure time very accurately, but the speed of the radar signal varies considerably with the soil type. In air, radar waves travel at 186,000 miles per second (or about 1 foot in one nanosecond). If you know the dielectric constant (from experience) or if you can look at the image of an object with a known depth, the GPR can be adjusted to read depth accurately. Software today allows construction of synthetic Hyperbola  for depth calculations.



                                                                                

                                                                                              Note: Photographs used above: Courtesy Georadar Inc.

The figure as under presents possible propagation paths for surface GPR energy. Principally, all these waves can be used to measure soil water content. In the following section, we focus on soil water content estimation using reflected and ground wave travel time data. In addition, we also discuss the estimation of water content using borehole GPR travel time data and using ground surface reflection amplitude data.

                                                                      



Measuring Soil Water Content with Reflected Waves:

Two classes of methods to estimate soil water content from reflected wave travel time data can be distinguished. The first class contains the methods that use a single antenna separation for soil water content estimation (e.g., soil water content estimation from scattering objects and traditional GPR sections). The second class contains the methods that require multiple measurements with different antenna separations.

Single or common offset reflection methods:

The energy that GPR transmits into the soil will be (partly) reflected when contrasts in soil permittivity are encountered. Figure hereunder shows an idealized GPR section measured with surface radar and a fixed antenna separation (single offset) over an anomaly (e.g., a water-filled pipe) having a different permittivity than the host material.



                                                                             


As GPR emits waves in all directions, reflected energy is measured before the GPR is directly over it. The reflected events in the radar section trace out a hyperbola because the reflected energy of the GPR measurement directly above the anomaly has the shortest travel distance (time) and all other waves will have a larger distance to travel. The average wave velocity in the soil determines the convexity of the reflection hyperbola; i.e., it determines how much longer the waves need to travel the extra distance. The average velocity between the ground surface and the anomaly, vsoil, can be determined from a GPR transect by fitting the following hyperbola to measured arrival times at several positions x


                                                                                              


where x is the position relative to the position of the scattering object (apex of the hyperbola), d is
the depth of the scattering object, and trw,x is the arrival time of the reflected wave at position x that has been zero time corrected. If the GPR section is measured with a significant antenna   separation, a, this should also be included in the velocity determination as follows :


                                                                             


Most common GPR analysis software provides routines where the velocity can be determined interactively by manually fitting hyperbola to the limbs of the reflections of the scattering object. The velocity can then be used to calculate soil permittivity and soil water content.
 

                                                              



Although it is simple and straightforward to determine velocity from scattering objects, it is a method that has not been used often for soil water content determination. The main drawback of this method is that it can only be used in soils where scattering objects can be observed in the GPR section. Even when scattering objects are present, this method only provides the average soil water content to the depth of the reflector; that is, the user has no control over the depth resolution of the soil water content measurements.

The accuracy of the single offset GPR reflection method for estimating water content under natural conditions is not yet well established. Some researchers have tested the concept of using GPR reflections under natural conditions to estimate water content using depth measurements to the reflectors obtained at discrete locations from, for example, noting lithologic transitions during drilling or water table observations. The use of single offset GPR reflection data for estimating spatially variable water content under naturally heterogeneous conditions at the field scale is a topic of active research.

Multi offset Reflection Methods :

Single offset measurements cannot be used to determine water content from reflecting soil layers if no information about the depth of the reflector is available. In that case, one can use a multi-offset GPR acquisition geometry to determine soil water content from radar reflections. Two commonly used multi-offset GPR acquisition geometries are called Common-MidPoint (CMP) and Wide Angle Reflection and Refraction (WARR) measurements as shown below:


                                                                           
 



In CMP acquisition, the distance between the antennas is increased stepwise while keeping a common midpoint. In WARR acquisition, the distance between the antennas is increased stepwise with the transmitter at a fixed position. A schematic outcome of a multi-offset GPR measurement is given as under:


                                                                            



Most of the commercial software provides solutions for multi offset recording of data and determination of permittivity of various layers encountered. Following is a typical profile shot using Common Depth Point technique.




                                                                                           


Three signals are clearly visible in the given profile, namely:

1. Air wave path signal. Always appears as an inclined straight line.

2. Path of signal reflected from the first layer interface.

3. Path of signal refracted on the first layer interface and reflected from the second layer interface.

To calculate characteristics of a medium (wave velocity or permittivity) by hodographs of received signals, software provides easy input windows:



                                                                                  


Initial values are input are inserted by the user and then the software performs the calculations providing permittivity and thickness of various layers.

The numerical calculated results will be displayed for each layer, namely: layer thickness in meters, electromagnetic wave velocity in cm/ns, permittivity, and mean root square value of approximation error. The less is the last value, the more accurate is the result.

Although multi-offset measurements are widely used in GPR data processing for determining velocity profiles with depth, there are some distinct disadvantages to soil water content determination with this method:

1. As with soil water content determination from single offset measurements, there is no control over the measurement depth resolution.

2. Multi-offset measurements are cumbersome to make and do not allow reconnaissance studies of soil water content variation.

3. For heterogeneous media, soil water content determined from multi-offset measurements is biased toward the common midpoint in the case of CMP acquisition geometry and toward the position of the fixed antenna in the case of WARR acquisition geometry.


Soil Water Content Measurements with the Ground Wave

The ground wave is the part of the radiated energy that travels between the transmitter and receiver through the top of the soil. The ground wave is detected by the GPR receiver, even in the absence of clearly reflecting soil layers. The evanescent character of the ground wave measured by the GPR receiver antenna at the soil surface requires that both the transmitter and the receiver be placed close to the soil surface.
 
The ground wave can easily be recognized on data collected using a multi-offset GPR acquisition geometry, by the observed linear relationship between antenna separation and ground wave travel times, which starts at the origin of the multi-offset measurement set. The slope of the ground wave in a multi-offset measurement is directly related to the ground wave velocity and can, therefore, be used for soil water content determination. The accuracy of this soil water content measurement technique was determined by researchers by analyzing a set of 24 multi-offset measurements collected using 225-MHz antennas and independent gravimetric samples, and was found to be fairly accurate.


                                                                      

Wide angle reflection and refraction (WARR) measurement recorded on loamy sand with the 225-MHz antennas. The velocity of the ground wave is v and of the air wave is c, both in meters per nanosecond.

Estimation of soil water content using multi-offset GPR measurements is cumbersome and time-consuming, as was mentioned before. The ground wave velocity can also be determined from a single offset GPR measurement, provided that the approximate arrival time of the ground wave is known from a multi-offset GPR measurement. Following procedure for soil water content mapping with the ground wave of GPR has therefore been suggested:

1. Identify an approximate ground wave arrival time for different antenna separations in a multi-offset GPR measurement,

2. Choose an antenna separation where the ground wave is clearly separated from the air and reflected waves and

3. Use this antenna separation for single offset GPR measurements and relate the changes in ground wave arrival time to changes in soil permittivity.

Although the results with ground wave data generally are promising, there are still some uncertainties associated with this method. An important, but unresolved issue is the effective measurement volume over which the ground wave averages. It is suggested that the influence depth is approximately one-half of the wavelength [  = c/(f )1/2], which would, for example, mean that for a central frequency of 225 MHz, the depth of influence could vary from 0.50 m (  = 4.0) to 0.22 m (  = 20.0).

Other drawbacks of using the GPR ground wave to estimate soil water content include the following:
 
1. It might be difficult to separate the ground wave arrival in the clutter of critically refracted and reflected waves.

2. It might be difficult to choose an antenna separation for which the arrival times of the air and ground wave can consistently be picked despite moving the antennas across a field with varying soil water content.

3. The ground wave is attenuated more quickly than other waves, which limits the range of antenna separation at which the ground wave can be observed.

Soil Water Content Measurements with Borehole GPR

For borehole GPR applications, the transmitting and receiving antenna are lowered into a pair of vertical access tubes. In the zero offset profile (ZOP) mode, the antennas are lowered such that their midpoints are always at the same depth. The arrival time of the direct wave between the boreholes and the known borehole separation is used to calculate the velocity and soil permittivity. The ZOP mode is an attractive approach for measuring the soil water content profile of the vadose zone with a high spatial resolution and a large sampling volume. Furthermore, each borehole GPR measurement requires only a couple of seconds and, therefore, the ZOP method is potentially capable of measuring transient processes within the unsaturated zone.

                                                                          



Soil water content can also be determined from a multi-offset profile (MOP). The first arrival times of all multi-offset measurements can be used to reconstruct a (tomographic) two-dimensional image of the soil water content distribution between the boreholes. To extract high-resolution, quantitative information from radar tomographic data, it is important to process the data as accurately as possible. The procedure for inverting radar tomographic data has been evolved, including methods for recognizing and correcting for errors caused by incorrect station geometry, incorrect zero time and zero time drift, geometric spreading, transmitter radiation pattern, transmitter amplitude, and high angle ray paths. The two-dimensional tomogram is obtained by discretizing the area between the boreholes in rectangular cells of constant velocity and determining the velocity of each cell by minimizing the difference between measured arrival times and arrival times calculated for raypaths passing through these cells. When necessary, three-dimensional tomograms can also be reconstructed. The drawback of acquiring the two-dimensional tomogram is the much longer time (typically several hours) required to obtain all the required measurements. Therefore, the MOP mode is best suited for steady-state water content conditions.
 


                                                                            

A comparison between estimates of soil water content (SWC) obtained from simple transformations of 200-  MHz zero offset profiling (ZOP) and multi-offset profiling (MOP) borehole ground penetrating radar (GPR) as well as cone penetrometer (CPT) data at the DOE Hanford Site in Washington (modified from Majer et al., 2002).

Recent publications show that borehole GPR is quickly becoming an important site-specific investigation tool in hydrogeological studies. An obvious but important difference from the other GPR methods presented in this review is that borehole GPR is not suitable as a reconnaissance tool. Despite the strong increase of borehole GPR applications, there are several points that require attention when using this technique:

1. To obtain quantitative information from borehole radar data, it is important to recognize and correct for errors caused by acquisition procedures and transmission characteristics, and to recognize the accuracy limitations.

2. It is important to consider the impact of refracted waves, especially those traveling in the soil, on the accuracy of soil water content measurements with borehole GPR.

3. The length of the borehole GPR antenna, the borehole separation distance, and the antenna frequencies have an impact on the maximum spatial resolution that can be achieved;

4. Soil heterogeneity affects the sampling volume of borehole GPR.
 

Soil Water Content Measurements with Surface Reflections

The measurement principle of soil water content measurements with air-launched surface reflections is illustrated as under.


                                                                  


The GPR antennas are operated at some distance above the ground by mounting them on a vehicle or a low-flying air platform. The soil property being measured is the reflection coefficient of the air soil interface, R, which is related to the soil permittivity,  soil, by

                                                                                                            




The reflection coefficient is determined from the measured amplitude, Ar, relative to the amplitude of a "perfect" reflector, Am, such as a metal plate larger than the footprint of the radar:

                                                                                                          


The determination of soil water content from the reflection coefficient requires accurate amplitude measurements. Clearly, the impact of surface roughness and soil water content profiles on the surface reflection coefficient are two key issues that need to be addressed when applying this technique. It is interesting to note that these problems are similar to those in active remote sensing, and it is possible that this soil water content measurement technique can profit from developments in the highly sponsored field of remote sensing.

Following picture shows a sample of data acquired with the 500-MHz GPR surface reflection measurement set up. Generally, the soil water content measured with GPR is similar to the soil water content measured with 0.20-m-long TDR sensors. However, there seems to be quite a large amount of very short distance variation in soil water content. Three likely explanations for the observed large variation of 0.10 m3 m-3 in soil water content include (i) the impact of the soil water content profile with depth on the reflection coefficient, (ii) the impact of surface roughness on the reflection coefficient, and (iii) reliability and accuracy of amplitude measurements.
 


                                                              

Water content profiles measured with air-launched 500-MHz ground penetrating radar (GPR) surface reflectivity method compared with measurements obtained using 20-cm-long time domain reflectometry (TDR) probes.



Conclusions :

The paper has discussed in details capabilities and limitations of using Ground Penetrating Radar for determining soil water content. Various examples from work carried out by researchers across the world clearly demonstrate that GPR can be used by experienced researchers to provide reliable estimates of subsurface water content. The exception to this is the GPR surface reflectivity method that, although promising, does not yet provide the accuracy and robustness of the other methods. In our opinion, the complexity of data acquisition and processing currently limits the applicability of GPR for water content estimation to the research community. However, continued development of specific and easy-to-use GPR instrumentation, coupled with increased experience and application of GPR for water content estimation, should help to bridge the gap between the advances made in the research community and the practical need for straightforward tools to accurately measure water content in high resolution and over large areas. Most likely, the borehole GPR method will continue to be used as an important research tool for studying infiltration processes and contaminant transport in the saturated and unsaturated zones. There may even be potential applications outside the research community for monitoring subsurface water content as part of groundwater protection strategies.
 



References:

Alumbaugh, D., P. Chang, L. Paprocki, J. Brainard, R.J. Glass, and C.A. Rautman. 2002. Estimating moisture contents in the vadose zone using cross-borehole ground penetrating radar: A study of accuracy and repeatability. Water Resour. Res. 38:1309.

Annan, A.P. 1973. Radar interferometry depth sounding: Part I. Theoretical discussion. Geophysics 38:557580.

Annan, A.P. 1996. Transmission, dispersion and GPR. J. Environ. Eng. Geophys. 0:125136.

Berktold, A., K.G. Wollny, and H. Alstetter. 1998. Subsurface moisture determination with the ground wave of GPR. p. 675680. In Proc. Int. Conf. on Ground-Penetrating Radar, 7th, Lawrence, KS. May 1998. Univ. of Kansas, Lawrence.

Binley, A., P. Winship, R. Middleton, M. Pokar, and J. West. 2001. High-resolution characterization of vadose zone dynamics using cross-borehole radar. Water Resour. Res. 37:26392652.

Binley, A., P. Winship, L.J. West, M. Pokar, and R. Middleton. 2002. Seasonal variation of moisture content in unsaturated sandstone inferred from borehole radar and resistivity profiles. J. Hydrol. (Amsterdam) 267:160172.

Bohidar, R.N., and J.F. Hermance. 2002. The GPR refraction method. Geophysics 67:14741485.

Buursink, M.L., J.W. Lane, Jr., W.P. Clement, and M.D. Knoll. 2002. Use of vertical-radar profiling to estimate porosity at two New England sites and comparison with neutron log porosity. In Proc. of SAGEEP '02, Las Vegas, NV. 1014 Feb. 2002. Environmental and Engineering Geophysical Society, Denver, CO.

Callies, U., A. Rhodin, and D.P. Eppel. 1998. A case study on variational soil moisture analysis from atmospheric observations. J. Hydrol. (Amsterdam) 212213:95108.

Cerveny, V., and J.E.P. Soares. 1992. Fresnel volume ray tracing. Geophysics 57:902915.

Chanzy, A., A. Tarussov, A. Judge, and F. Bonn. 1996. Soil water content determination using a digital ground-penetrating radar. Soil Sci. Soc. Am. J. 60:13181326

Dalton, F.N., W.N. Herkelrath, D.S. Rawlins, and J.D. Rhoades. 1984. Time-domain reflectometry: Simultaneous measurement of soil water content and electrical conductivity with a single probe. Science 224:989990.

Dannowski, G., and U. Yaramanci. 1999. Estimation of water content and porosity using combined radar and geoelectrical measurements. Eur. J. Environ. Eng. Geophys. 4:7185.

Davis, J.L., and A.P. Annan. 1989. Ground-penetrating radar for high resolution mapping of soil and rock stratigraphy. Geophys. Prospect. 37:531551.

Davis, J.L., and A.P. Annan. 2002. Ground penetrating radar to measure soil water content. p. 446463. In J.H. Dane and G.C. Topp (ed.) Methods of soil analysis. Part 4. SSSA Book Ser. 5. SSSA, Madison WI.

Day-Lewis, F.D., J.M. Harris, and S. Gorelick. 2002. Time-lapse inversion of crosswell radar data. Geophysics 67:17401752.

Debye, P. 1929. Polar molecules. Dover Publ., Mineola, NY.

Dix, C.H. 1955. Seismic velocities from surface measurements. Geophysics 20:73.

Dobson, M.C., F.T. Ulaby, M.T. Hallikainen, and M.A. El-Rayes. 1985. Microwave dielectric behaviour of wet soil. Part II. Dielectric mixing models. IEEE Trans. Geosci. Remote Sens. 23:3546.

Du, S. 1996. Determination of water content in the subsurface with the ground wave of ground penetrating radar. Ph.D. thesis. Ludwig-Maximilians-Universitt, Munich, Germany.

Endres, A.L., W.P. Clement, and D.L. Rudolph. 2000. Ground penetrating radar imaging of an aquifer during a pumping test. Ground Water 38:566576.

Entekhabi, D., G.R. Asrar, A.K. Betts, K.J. Beven, R.L. Bras, C.J. Duffy, T. Dunne, R.D. Koster, D.P. Lettenmaier, D.B. McLaughlin, W.J. Shuttleworth, M.Th. van Genuchten, M.Y. Wei, and E.F. Wood. 1999. An agenda for land surface hydrology research and a call for the second international hydrological decade. Bull. Am. Meteorol. Soc. 80:20432058.

Eppstein, M.J., and D.E. Dougherty. 1998. Efficient three-dimensional data inversion: Soil characterization and moisture monitoring from cross-well ground-penetrating radar at a Vermont test site. Wat. Resour. Res. 34:18891900.

Famiglietti, J.S., J.A. Deveraux, C.A. Laymon, T. Tsegaye, P.R. Houser, T.J. Jackson, S.T. Graham, M. Rodell, and P.J. van Oevelen. 1999. Ground-based investigation of soil moisture variability within remote sensing footprints during the Southern Great Plains 1997 (SGP97) hydrology experiment. Wat. Resour. Res. 35:18391851.

Ferr, P.A., D.L. Rudolph, and R.G. Kachanoski. 1996. Spatial averaging of water content by time domain reflectometry: Implications for twin rod probes with and without dielectric coatings. Wat. Resour. Res. 32:271279.

Ferr, P.A., G. von Glinski, and L.A. Ferr. 2003. Monitoring the maximum depth of drainage in response to pumping using borehole ground penetrating radar. Available at www.vadosezonejournal.org. Vadose Zone J. 2:511518

Friedman, S.P. 1998. A saturation degree-dependent composite spheres model for describing the effective dielectric constant of unsaturated porous media. Wat. Resour. Res. 34:29492961.

Galagedara, L.W., G.W. Parkin, J.D. Redman, and A.L. Endres. 2002. Temporal and spatial variation of soil water content measured by borehole GPR under irrigation and drainage. Proc. of the Ninth Conf. on Ground Penetrating Radar. Proc. SPIE 4758:180185.

Garambois, S., P. Senechal, and H. Perroud. 2002. On the use of combined geophysical methods to assess water content and water conductivity of near-surface formations. J. Hydrol. (Amsterdam) 259:3248.

Gilson, E.W., J.D. Redman, J. Pilon, and A.P. Annan. 1996. Near surface applications of borehole radar. p. 545553. In Proc. of the Symp. on the Application of Geophysics to Engineering and Environmental Problems, Keystone, CO. Environ. Eng. and Geophys. Soc., Denver, CO.

Greaves, R.J., D.P. Lesmes, J.M. Lee, and M.N. Toksoz. 1996. Velocity variations and water content estimated from multi-offset, ground-penetrating radar. Geophysics 61:683695.

Grote, K., S.S. Hubbard, and Y. Rubin. 2002. GPR monitoring of volumetric water content in soils applied to highway construction and maintenance. Leading Edge Explor. 21:482485.

Grote, K., S.S. Hubbard, and Y. Rubin. 2003. Field-scale estimation of volumetric water content using GPR ground wave techniques. Wat. Resour. Res.

Hammon, W.S. III, X. Zeng, R.M. Corbeanu, and G.A. McMechan. 2003. Estimation of the spatial distribution of fluid permeability from surface and tomographic GPR data and core, with a 2-D example from the Ferron Sandstone, Utah. Geophysics 67:15051515.

Hasted, J.B. 1973. Aqueous dielectrics. Chapman and Hall, London.

Heimovaara, T.J., and W. Bouten. 1990. A computer-controlled 36-channel time domain reflectometry system for monitoring soil water contents. Water Resour. Res. 26:23112316.

Herkelrath, W.N., S.P. Hamburg, and F. Murphy. 1991. Automatic, real-time monitoring of soil moisture in a remote field area with time domain reflectometry. Water Resour. Res. 27:857864.

Hubbard, S.S., K. Grote, and Y. Rubin. 2002. Mapping the volumetric soil water content of a California vineyard using high-frequency GPR ground wave data. Leading Edge Explor. 21:552559.
 
Hubbard, S.S., J.E. Peterson, Jr., E.L. Majer, P.T. Zawislanski, K.H. Williams, J. Roberts, and F. Wobber. 1997. Estimation of permeable pathways and water content using tomographic radar data. Leading Edge Explor. 16:16231630.

Huisman, J.A., and W. Bouten. 2003. Accuracy and reproducibility of measuring soil water content with the ground wave of ground penetrating radar. J. Environ. Eng. Geophys. 8:6573.

Huisman, J.A., J.J.J.C. Snepvangers, W. Bouten, and G.B.M. Heuvelink. 2002. Mapping spatial variation in surface soil water content: Comparison of ground-penetrating radar and time domain reflectometry. J. Hydrol. (Amsterdam) 269:194207.

Huisman, J.A., J.J.J.C. Snepvangers, W. Bouten, and G.B.M. Heuvelink. 2003. Monitoring temporal development of spatial soil water content variation: Comparison of ground penetrating radar and time domain reflectometry. Available at www.vadosezonejournal.org. Vadose Zone J. 2:519�529.

Huisman, J.A., S. S. Hubbard, J. D. Redman and A. P. Annan. 2003. Measuring soil water content with Ground Penetrating Radar. Available at www.vadosezonejournal.org. Vadose Zone J. 2:476-491.

Huisman, J.A., C. Sperl, W. Bouten, and J.M. Verstraten. 2001. Soil water content measurements at different scales: Accuracy of time domain reflectometry and ground-penetrating radar. J. Hydrol. (Amsterdam) 245:4858.

Jackson, T.J., J. Schmugge, and E.T. Engman. 1996. Remote sensing applications to hydrology: Soil moisture. Hydrol. Sci. J. 41:517530.

Jacobsen, O.H., and P. Schjnning. 1994. Comparison of TDR calibration functions for soil water determination. p. 923. In Proc. Time Domain Reflectometry, Applications in Soil Science, Research Centre Foulum, Denmark. 16 Sept. 1994. SP-Report 25-33. Danish Inst. of Plant and Soil Sci., Tjele, Denmark.

Jones, S.B., and S.P. Friedman. 2000. Particle shape effect on the effective permittivity of anisotropic or isotropic media consisting of aligned or randomly oriented ellipsoidal particles. Water Resour. Res. 36:28212833.

Knoll, M.D., and W.P. Clement. 1999. Vertical radar profiling to determine dielectric constant, water content and porosity values at well locations. p. 821830. In Proc. of SAGEEP'99, Oakland, CA. Environ. Eng. and Geophys. Soc., Denver, CO.

Ledieu, J., P. De Ridder, P. De Clerck, and S. Dautrebande. 1986. A method of measuring soil moisture by time domain reflectometry. J. Hydrol. (Amsterdam) 88:319328.

Lehmann, F., and A.G. Green. 1999. Semi-automated georadar data acquisition in three dimensions. Geophysics 64:719731

Lesmes, D.P., R. Herbstzuber, and D. Wertz. 1999. Terrain permittivity mapping: GPR measurements of near-surface soil moisture. p. 575582. In Proc. SAGEEP'99, Oakland, CA. Environ. Eng. and Geophys. Soc., Denver, CO.

Majer, E.L., K.H. Williams, J.E. Peterson, and T.M. Daley. 2002. High resolution imaging of vadose zone transport using crosswell radar and seismic methods. Rep. LBNL 49022. Lawrence Berkeley National Laboratory, Berkeley, Ca.

Merz, B., and A. Bardossy. 1998. Effects of spatial variability on the rainfall runoff process in a small loess catchment. J. Hydrol. (Amsterdam) 212213:304317.

Merz, B., and E.J. Plate. 1997. An analysis of the effects of spatial variability of soil and soil moisture on runoff. Water Resour. Res. 33:29092922.

Nakashima, Y., H. Zhou, and M. Sato. 2001. Estimation of groundwater level by GPR in an area with multiple ambiguous reflections. J. Appl. Geophys. 47:241249.

Nelson, S.O. 1994. Measurement of microwave dielectric properties of particulate materials. J. Food Eng. 21:365384.

Olhoeft, G.R. 1987. Electrical properties from 10-3 to 10+9 Hz: Physics and chemistry. Proc. Phys. Chem. Porous Media II 154:281298.

Ollson, O., L. Falk, O. Forslund, L. Lundmark, and E. Sandberg. 1992. Borehole radar applied to the characterization of hydraulically conductive fracture zones in crystalline rock. Geophys. Prospect. 40:109142.
 
Or, D., and J.M. Wraith. 1999. Temperature effects on soil bulk dielectric permittivity measured by time domain reflectometry: A physical model. Water Resour. Res. 35:371383.

Parkin, G., D. Redman, P. von Bertoldi, and Z. Zhang. 2000. Measurement of soil water content below a wastewater trench using ground-penetrating radar. Water Resour. Res. 36:21472154.

Pauwels, V.R.N., R. Hoeben, N.E.C. Verhoest, and F.P. De Troch. 2001. The importance of the spatial patterns of remotely sensed soil moisture in the improvement of discharge predictions for small-scale basins through data assimilation. J. Hydrol. (Amsterdam) 251:88102.

Peterson, J.E., Jr. 2001. Pre-inversion corrections and analysis of radar tomographic data. J. Environ. Eng. Geophys. 6:118.

Redman, J.D., J.L. Davis, L.W. Galagedara, and G.W. Parkin. 2002. Field studies of GPR air launched surface reflectivity measurements of soil water content. Proc. of the Ninth Conf. on Ground-Penetrating Radar. Proc. SPIE 4758:156161.

Robinson, D.A., S.B. Jones, J.M. Wraith, D. Or, and S.P. Friedman. 2003. A review of advances in dielelectric and electric conductivity measurements using time domain reflectometry. Available at www.vadosezonejournal. Vadose Zone J. 2:444475

Roth, K., R. Schulin, H. Flhler, and W. Attinger. 1990. Calibration of time domain reflectometry for water content measurement using a composite dielectric approach. Water Resour. Res. 26:22672273.

Rubin, Y. 2003. Applied stochastic hydrogeology. Oxford University Press, Oxford, UK.

Rucker, D.F., and P.A. Ferr. 2003. Near-surface water content estimation with borehole ground penetrating radar using critically refracted waves. Available at www.vadosezonejournal.org. Vadose Zone J. 2:247252.

Shinn, J.D., D.A. Timian, R.M. Morey, G. Mitchell, C.L. Antle, and R. Hull. 1998. Development of a CPT deployed probe for in situ measurement of volumetric soil moisture content and electrical resistivity. Field Anal. Chem. Tech. 2:103109.

Sperl, C. 1999. Determination of spatial and temporal variation of the soil water content in an agro-ecosystem with ground-penetrating radar. (In German.) Ph.D. thesis. Technische Universitat Mnchen, Munich, Germany.

Stoffregen, H., T. Zenker, and G. Wessolek. 2002. Accuracy of soil water content measurements using ground penetrating radar: Comparison of ground penetrating radar and lysimeter data. J. Hydrol. (Amsterdam) 267:201206.

Tillard, S., and J.-C. Dubois. 1995. Analysis of GPR data: Wave propagation velocity determination. J. Appl. Geophys. 33:7791.

Topp, G.C., J.L. Davis, and A.P. Annan. 1980. Electromagnetic determination of soil water content: Measurements in coaxial transmission lines. Water Resour. Res. 16:574582.

Topp, G.C., M. Yanuka, W.D. Zebchuk, and S.J. Zegelin. 1988. Determination of electrical conductivity using time domain reflectometry: Soil and water experiments in coaxial lines. Water Resour. Res. 24:945952.

Ulaby, F.T., P.C. Dubois, and J. van Zyl. 1996. Radar mapping of surface soil moisture. J. Hydrol. (Amsterdam) 184:5784.

van Dam, R.L., W. Schlager, M.J. Dekkers, and J.A. Huisman. 2002. Iron oxides as a cause of GPR reflections. Geophysics 67:536545.

van Oevelen, P.J. 2000. Estimation of areal soil water content through microwave remote sensing. Ph.D. thesis. Wageningen University, The Netherlands.

van Overmeeren, R.A., S.V. Sariowan, and J.C. Gehrels. 1997. Ground penetrating radar for determining volumetric soil water content; results of comparative measurements at two sites. J. Hydrol. (Amsterdam) 197:316338.

Vasco, D.W., J.E. Peterson, and K.H. Lee. 1997. Ground-penetrating radar velocity tomography in heterogeneous and anisotropic media. Geophysics 62:17581773.

Vellidis, G., M.C. Smith, D.L. Thomas, and L.E. Asmussen. 1990. Detecting wetting front movement in a sandy soil with ground-penetrating radar. Trans. ASAE 33:18671874.
 
Weiler, K.W., T.S. Steenhuis, J. Boll, and K.-J.S. Kung. 1998. Comparison of ground penetrating radar and time domain reflectometry as soil water sensors. Soil Sci. Soc. of Am. J. 62:12371239.

West, L.J., K. Handley, Y. Huang, and M. Pokar. 2003. Radar frequency dielectric dispersion in sandstone: Implications for determination of moisture and clay content. Water Resour. Res. 39. DOI:10.1029/2001WR000923.

Wyseure, G.C.L., M.A. Mojid, and M.A. Malik. 1997. Measurement of volumetric water content by TDR in saline soils. Eur. J. Soil Sci. 48:347354.

Yilmaz, O. 1987. Seismic data processing. SEG Investigations in Geophysics 2. Society of Exploration Geophysicists, Tulsa, OK.

Zhou, C., L. Liu, and J.W. Lane. 2001. Nonlinear inversion of borehole-radar tomography data to reconstruct velocity and attenuation distribution in earth materials. J. Appl. Geophys. 47:271284.