|
|
||||||||
Department of Hydrology and Water Resources, University of Arizona. P.O. Box 210110, Tucson, AZ 85721
* Corresponding author (druck{at}hwr.arizona.edu)
Received 3 December 2002.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: BGPR, borehole ground penetrating radar EM, electromagnetic TDR, time domain reflectometry ZOP, zero-offset profiling
| INTRODUCTION |
|---|
|
|
|---|
![]() | [1] |
(cm3 cm-3), can be determined from the apparent dielectric permittivity using the linearized form of the empirical relationship presented by Topp et al. (1980) as given by Ferré et al. (1996):
![]() | [2] |
The preceding analysis can be used to determine the water content at any depth if the direct wave is associated with the first energy to arrive at the receiver at that depth. However, the direct route does not always have the shortest travel time. Rather, under some conditions, a critically refracted wave may arrive before the direct arrival (Fig. 1). Critical refraction occurs at any interface across which the velocity increases (Sheriff and Geldart, 1995). The contrast between the EM wave velocity in air (vair = 0.3 m ns-1) and in soil (0.17 > vsoil > 0.05 m ns-1) gives rise to critical refraction (Bohidar and Hermance, 2002) when the antennae are located below the ground surface. Critically refracted energy arrives at the receiver through secondary waves generated at the surface (Parkin et al., 2000), where a head wave is created. The receiver can intersect the head wave at any depth.
|
In this study, equations are derived to relate the refraction termination depth to the velocity of EM energy through the near-surface soil and to the antennae separation. This is similar to the development presented by Hammon et al. (2002). The analysis is extended here to determine the near-surface water content and to develop guidelines to predict the refraction termination depth.
It is important to note that these developments explicitly assume that the dielectric permittivity of the medium is both homogeneous and isotropic in the shallow subsurface. Using existing equipment to collect ZOP data, it is not possible to determine dielectric anisotropy. Similarly, ZOP measurements cannot be used to determine horizontal variations in dielectric permittivity between BGPR antennae. The precision with which vertical heterogeneity due to soil layering or water content distributions can be determined is a function of the sample depth interval. Presumably, there is some minimum thickness that can be resolved given the nature of signals generated by dipole antennae. In addition, we assumed that attenuation does not affect the ability to identify first arrivals on the BGPR trace. Under some conditions, refracted waves may experience greater attenuation than direct waves because of to their longer travel paths and transmission through different media. This could result in a termination depth that is controlled primarily by attenuation rather than relative velocities of adjacent layers.
| THEORY |
|---|
|
|
|---|
![]() | [3] |
For the case of BGPR operated just below the ground surface, v1 and v2 are the velocities in the soil immediately below the ground surface and in the air, respectively, and the normal is a line that is perpendicular to the ground surface. For a given velocity contrast, only one incident angle will give rise to a wave that has an angle of refraction of 90° from the normal and travels along the ground surface at the velocity of air. This angle of incidence is referred to as the critical angle, ic, and the refracted wave is known as the critical refraction. The critical angle can be defined from Eq. [3] as
![]() | [4] |
Considering BGPR antennae located close to the ground surface (Fig. 1), the travel time of the critically refracted wave (trefr) is equal to the sum of the travel times of the propagating EM wave from the transmitter to the airsoil interface (ttx-surface), along the surface (tair), and down from the interface to the receiving antennae (tsurface-rx):
![]() | [5] |
From Huygen's principle, the angle of refraction from the interface to the receiving antennae is also equal to the critical angle. As a result, the distance that the wave travels to or from the ground surface is
![]() | [6] |
The travel time over this distance is
![]() | [7] |
The distance that the wave travels in the air above the ground surface is
![]() | [8] |
The travel time along the ground surface is
![]() | [9] |
Substituting Eq. [7] and [9] into Eq. [5] yields
![]() | [10] |
The direct travel time between the antennae at the same depth is
![]() | [11] |
Note that while tdirect is independent of z, trefr increases linearly with z. This linear increase in trefr with z is a useful characteristic for identifying critical refractions in ZOP profiles.
Critically refracted waves can arrive before direct waves when the antennae are near the ground surface (Fig. 2). However, below the refraction termination depth, zrtd, the travel time of the direct wave is shorter than that of the critically refracted wave. The refraction termination depth can be determined by equating the travel times of the critically refracted and direct waves (Equations 10 and 11) and solving for zrtd:
![]() | [12] |
|
![]() | [13] |
Alternatively, if the travel time profile resembles that shown in Fig. 2, then zrtd can be obtained directly from the profile. Then vsoil can be found by rearranging Eq. [12] to
![]() | [14] |
From the definition of sinic, the value of cosic can be expressed as
![]() | [15] |
Substituting Eq. [15] into Eq. [14] and solving for vsoil gives
![]() | [16] |
It should be noted that Eq. [16] requires a precise measure of zrtd. In practice, the sampling depth interval (depth that the antennae are lowered between readings) is commonly 0.25 m or more. The actual zrtd for the example system is 1.061 m. Typical measurement resolution gives rise to uncertainty as great as one-half the sampling depth interval: 1.061 ± 0.125. Water content estimates, considering these uncertainties, range from 0.08 to 0.33 cm3 cm-3, with the correct water content being 0.17 cm3 cm-3. Moreover, if the shallow subsurface is layered, or water content varies with depth, the calculated zrtd may be deeper than the base of the near-surface layer. In this case, the inferred water content applies only above a depth associated with an abrupt change in measured velocity.
Alternatively, the slope of the travel time profile near the ground surface,
, can be used to estimate vsoil by taking the derivative of Eq. [10] with respect to depth:
![]() | [17] |
![]() | [18] |
Substituting Eq. [15] and [18] into Eq. [17] and solving for vsoil gives
![]() | [19] |
From Fig. 2, the slope of the travel time profile near the ground surface is 18.9 ns m-1 (20 ns/1.061 m), equating to a soil velocity of 0.10 m ns-1. Using Eq. [2] gives a water content of 0.17 cm3 cm-3. Actual measurements will give rise to some uncertainty of the slope. Choosing an arbitrary 5% underestimation (
= 18) or overestimation (
= 19.8), the calculated water contents range from 0.16 and 0.19, respectively. It is suggested that linear regression techniques be used to properly assess the slope error to give some measure of the accuracy of shallow water content measurements. Because zrtd may not be evident on a travel time profile plot, especially if the water content varies with depth in the shallow subsurface, we recommend using the slope analysis, which relies only on measurements made near the ground surface.
Because first-arriving critically refracted waves can complicate BGPR profiling applications, it is reasonable to ask whether a separation distance can be chosen that eliminates critically refracted waves as first arrivals. Equation [12] shows that zrtd increases with increasing antennae separation distance, x. Therefore, reducing the antennae separation will minimize the depth over which critically refracted waves arrive before direct waves. However, to allow for simple ray path analysis of travel times, the antennae separation distance must be greater than one-half of the wavelength of the center frequency in air. For example, a 100-MHz EM wave has a wavelength of 3 m in air, so the antennae separation distance when using 100-MHz antennae should be at least 1.5 m. Theoretically, a 900-MHz antenna has a minimum antennae separation of 0.18 m. However, this small separation is not practical for most field applications, given typical borehole annuli are on the order of tens of centimeters and vertical borehole deviations could be significant. The minimum zrtd is associated with a low water content, high velocity condition. Considering a medium with a volumetric water content of 0.05 cm3 cm-3, the minimum separation allowable for 100-MHz antennae results in a zrtd of 0.43 m. This result suggests that it is not possible to eliminate completely first-arriving critical refractions for measurements made very near the ground surface.
The preceding analysis only considered critical refraction at the airsoil interface. However, critical refraction can occur whenever EM waves traveling through a low-velocity medium contact a boundary with a higher-velocity medium. There are many examples of this below the ground surface. For example, the antennae may be located within a high water content, fine-textured layer that is in contact with a lower water content, coarse-textured layer. Or, the antennae may be located below the water table and the transmitted energy may refract critically at a boundary within unsaturated regions above the capillary fringe. Conclusions that are identical to those made for the impacts of refractions at the ground surface can be drawn regarding critical refractions at depth. That is, each layer will have a refraction termination depth below (or above) the layer boundary. Furthermore, if the critical refractions can be identified on the travel time profile, the slope of the travel time profile can be used to determine the water content in the high water content layers. However, given that the refraction termination depth increases with the contrast in velocities between adjacent layers, it is likely that the critical refractions will primarily affect measurements within layers that contrast sharply with adjacent layers.
| FIELD EXPERIMENT |
|---|
|
|
|---|
The infiltration experiment was conducted at the Western Campus Agricultural Center, an experimental farm along the bank of the Santa Cruz River in Tucson, AZ. During the experiment, water was applied to the ground surface at a constant rate of 0.0136 m h-1 for 66 h (Rucker and Ferré, 2002). For this study, BGPR measurements were made using two access tubes located within the infiltration gallery (Fig. 3). A Sensors and Software PulseEkko100 system with 100-MHz antennae (Sensors and Software, Mississauga, ON) was used for all BGPR measurements. Postprocessing of first-arrival picking was accomplished through an automated picking software, Picker, provided by Sensors and Software.
|
In addition to measurements with BGPR, two buried time domain reflectometry (TDR) probes were located within the infiltration gallery. One TDR probe was located outside of the gallery, and one multilevel TDR sensor was located within the gallery (Fig. 3). The buried TDR probes were placed perpendicular with the ground surface within hand-augered holes that were 10 cm in diameter and backfilled with native material. The multilevel TDR (Environmental Sensors Inc., Victoria, BC, Canada) was installed at the surface and extended to 1.2 m below ground surface. The multilevel TDR is discretized into five segments, with the segments having lengths of 0.15, 0.15, 0.3, 0.3, and 0.3 m from top to bottom. Before infiltration began, the TDR-measured volumetric water contents for these depth intervals were 0.163, 0.123, 0.125, 0.13, and 0.17 cm3 cm-3. These values give a length-weighted water content of 0.145 cm3 cm-3. Once infiltration began, the multilevel TDR malfunctioned, so measurements with this instrument were not available for comparison with later travel time profiles.
| RESULTS |
|---|
|
|
|---|
1 = 19.10 ± 0.09 ns m-1, where the uncertainty represents one standard deviation above and below the mean. This yields a mean velocity of 0.0989 m ns-1, with velocities ranging from 0.0985 to 0.0993 m ns-1 over one standard deviation of the value of
1. The mean volumetric water content determined from Eq. [2] using this slope was 0.170 cm3 cm-3, which is in good agreement with that measured using the multilevel TDR. The range of water contents based on the range of slopes measured was 0.173 to 0.176 cm3 cm-3. Similarly, high resolution was found for other measurements based on the uncertainty of the fitted slope. The velocity measured just before infiltration ceased (see "After Infiltration" in Fig. 4) was 0.056 m ns-1, which equates to a volumetric water content of 0.44 cm3 cm-3. The profile measured 1 wk after infiltration ceased yielded a volumetric water content of 0.25 cm3 cm-3. The value of zrtd calculated using Eq. [12] is reported in Table 1 for each travel time profile.
|
|
During the infiltration experiment, the volumetric water content was also measured using three buried TDR probes (Fig. 5). The measurement depths were 1.2, 1.6, and 1.67 m below ground surface. All of the TDR time series showed similar response as the wetting front moved through the soil profile. At early time, the water content was constant with time at each depth, although it differs among probes. The wetting front reached all of the TDR probes at approximately the same time. After this time, the water content increased to a maximum value and remained constant until infiltration ceased. Throughout the remaining measurement time, the water content decreased with time to a constant value. The two TDR probes located within the infiltration gallery, at depths of 1.6 m and greater, showed lower maximum water contents (0.28 cm3 cm-3) than the TDR probe buried just outside of the gallery at 1.2 m depth (0.37 cm3 cm-3).
|
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Lambot, J. Rhebergen, I. van den Bosch, E. C. Slob, and M. Vanclooster Measuring the Soil Water Content Profile of a Sandy Soil with an Off-Ground Monostatic Ground Penetrating Radar Vadose Zone J., November 1, 2004; 3(4): 1063 - 1071. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Lambot, M. Antoine, I. van den Bosch, E. C. Slob, and M. Vanclooster Electromagnetic Inversion of GPR Signals and Subsequent Hydrodynamic Inversion to Estimate Effective Vadose Zone Hydraulic Properties Vadose Zone J., November 1, 2004; 3(4): 1072 - 1081. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. F. Rucker and T. P. A. Ferre Parameter Estimation for Soil Hydraulic Properties Using Zero-Offset Borehole Radar: Analytical Method Soil Sci. Soc. Am. J., September 1, 2004; 68(5): 1560 - 1567. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. F. Rucker and T. P. A. Ferre Correcting Water Content Measurement Errors Associated with Critically Refracted First Arrivals on Zero Offset Profiling Borehole Ground Penetrating Radar Profiles Vadose Zone J., February 1, 2004; 3(1): 278 - 287. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Huisman, J. A. Huisman, S. S. Hubbard, J. D. Redman, and A. P. Annan Measuring Soil Water Content with Ground Penetrating Radar: A Review Vadose Zone J., November 1, 2003; 2(4): 476 - 491. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Serbin and D. Or Near-Surface Soil Water Content Measurements Using Horn Antenna Radar: Methodology and Overview Vadose Zone J., November 1, 2003; 2(4): 500 - 510. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||