|
|
||||||||
a EMBRAPA, Agricultural Instrumentation, P.O. Box 741, Sao Carlos-SP, 13560-970, Brazil
b Department of Land, Air and Water Resources, Hydrology Program, 123 Veihmeyer Hall, University of California, Davis, CA 95616
Corresponding author (jwhopmans{at}ucdavis.edu)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: PR, penetration resistance PVC, polyvinyl chloride TDR, time domain reflectometry
| INTRODUCTION |
|---|
|
|
|---|
Soil scientists have related changes in PR as caused by tillage, traffic, or soil genetic pans to root growth, crop yields, and soil physical properties. For example, correlation between PR and crop root growth and water and nutrient exploration have been obtained (Stelluti et al., 1998), and cone penetrometers have been used extensively in soil science studies to identify natural and induced compacted layers (Henderson, 1989) or to predict related soil properties (Ayers and Bowen, 1987).
Among the soil parameters that affect PR, soil water content and bulk density are the most significant (Vazquez et al., 1991). For example, Stitt et al. (1982) conducted a comprehensive study of factors affecting PR in coarse-textured soils in the Atlantic Coastal Plain, and used stepwise regression to relate mechanical impedance to various measured soil properties. The highest correlation coefficients were found for a regression model that included soil water content, soil particle roughness and bulk soil density. Shaw et al. (1942) concluded that soil moisture is the dominant factor influencing the force required to push a penetrometer into the soil, with PR increasing as the moisture content decreased. In an experimental study by Henderson et al. (1988) it was found that PR was only slightly affected with a decrease of soil water content to
70% of field capacity. However, the PR increased exponentially with a further reduction of the water content of the sandy soil. This study showed that PR increased with an increase of bulk density across the whole measured water content range. However, because soil moisture varies both spatially and temporally and is only one of the soil variables related to PR, the utility of using PR to determine compaction effects is marginal. Moreover, interpretation of penetrometer data is difficult because water content or density measurements can generally not be taken at the exact same spatial location as the penetration resistance measurement.
Considering the strong dependence of PR on soil water content within a soil profile and across a field, it would be beneficial if both soil water content and soil resistance could be measured simultaneously at the same location and depth with a single probe. Among available techniques for soil water content measurements, TDR is the most attractive. Advantages of TDR over other soil water content measurement techniques include (i) potential for variable measurement volume size using different probe sizes and geometry, (ii) the use of the same probe for both laboratory and field measurements, (iii) small influence of dissolved solutes on the TDR moisture measurement within a low salinity range, (iv) its potential for automatic data acquisition and multiplexing, and (v) that it does not pose a radiation hazard.
Most commercial TDR equipment uses standard waveguides or probes with a usual length of 15 to 30 cm. The soil water content value is obtained from calibration curves using the travel time of an electromagnetic wave along the waveguide to determine the bulk dielectric constant of the soil (Topp et al., 1980). A minimum probe length is controlled by the rise time of the electromagnetic square wave reflected at the beginning and end of the TDR probe (Nissen et al., 1998). Petersen at al. (1995) examined the importance of probe length and diameter, distance between wave guides, and horizontal installation depth. They obtained excellent waveforms using a 5-cm probe in a coarse sandy soil with a water content of 0.21 cm3 cm-3. Kelly et al. (1995) obtained an accuracy of 0.035 cm3 cm-3 using TDR probes as short as 2.5 cm and a high-band width TDR system of 20 MHz. Amato and Ritchie (1995) experimented with short probes ranging in length from 1 to 15 cm. They concluded that at water content values of 0.07 cm3 cm-3 with travel times larger than 100 ps, the error in the water content was less than three volume percent. However, for water content measurements in drier soils with shorter travel times, errors were larger than 4 to 5%. Malicki et al. (1992) and Sri Ranjan and Domytrak (1997) described successfully the use of TDR miniprobes (5 cm long) for a clay loam soil.
Selker et al. (1993) introduced a serpentine type surface probe (10 by 10 cm) by imbedding the conductor and ground wires of the TDR probe within an acrylic plate, enabling miniaturization of TDR probes for high spatial resolution measurements. For the coiled probe developed by Nissen et al. (1998) the conductor wire was guided around a cylindrical PVC rod with four straight ground wires along the PVC rod. Their TDR probe allowed a reduction in probe length of a factor of five without a loss in sensitivity. To avoid short-circuiting, the conductor and ground wires were lacquer-coated. In both designs, the increased conductor wire lengths ensured long enough travel times for accurate water content measurements despite the smaller measured bulk soil volume. Both designs (serpentine and coil) are innovative compared with the traditional two, three, or four rod probes and bring many new TDR applications.
The concept of a combined measurement of penetration resistance and water content was previously presented (Ward, 1994; Young et al., 1998; Adams et al., 1998; Newman and Hummel, 1999; Vaz et al., 1999), but to date details regarding construction and calibration for different soils has been limited. For that reason, the objective of this work was to design, construct, and evaluate a coiled TDR probe to be used in combination with a cone penetrometer to determine water content and penetration resistance along a soil profile in a field setting. After analysis of the testing in the laboratory, the combined penetrometerTDR soil moisture probe measurement results for the field are presented.
| MATERIALS AND METHODS |
|---|
|
|
|---|
coaxial cable is passed through the hollow steel shaft of the penetrometer probe and connected to a cable tester (Tektronix 1502C, Tektronix, Beaverton, OR). Both copper wires were soldered to the corresponding conductor and ground of the coaxial cable in two opposing 2-mm access holes, right above the coil, as shown in Fig. 1a. The spaces between the wires of the coil and the two access holes were filled with an epoxy resin (2-Ton crystal clear epoxy, Devcon, Riviera Beach, FL) and smoothed to avoid the creation of air spaces between the wires during soil insertion. However, probesoil contact is also largely affected by the probe operator as straight vertical insertion is required. Figure 2
shows the details of the combined TDRcone penetrometer probe. Cone and probing rod sizes satisfy the American Society of Agricultural Engineers Standards (American Society of Agricultural Engineers, 1994).
|
|
coaxial cable mounted in an epoxy resin base as shown in Fig. 1c.
Laboratory Calibration
The waveform or trace is transferred from the cable tester to a personal computer through the RS232 serial port and analyzed. The trace (Fig. 3a and 3b)
is a visualization of the amplitude of a reflected pulsed electromagnetic wave as a function of propagation or travel time along the TDR probe. The trace can be regarded as a signature of the physical status of the soil, and it can be shown that knowledge of the travel time is sufficient to determine the bulk material dielectric constant of the soil (Topp et al., 1980). Travel times and bulk dielectric constant are determined by identification of the first and second reflection at the beginning and end of both TDR probe types. The procedural steps used to identify these reflection points were (i) smoothing of the waveform using a moving average approach and (ii) computation of the first and second time-derivatives of the smoothed waveform. Electric shorting of the probe at its beginning and end was used to confirm correct identification of the beginning and end reflection points of the TDR probes (Hook et al., 1992; Kelly et al., 1995; Nissen et al., 1998). Since the dielectric constant measured by the combined probe is a weighted average dielectric constant of the soil and the PVC and epoxy material between the TDR wires, a conventional TDR probe with two parallel wave guides of 5-cm length was used (Fig. 1c) to estimate the bulk soil dielectric constant of the investigated soil samples (
soil). Using the mixing model approach of Nissen et al. (1998), the dielectric constant measured by the coiled TDR probe (
coil) can be related to the soil dielectric constant as determined by the conventional probe (
soil)
![]() | (1) |
|
probe) and the bulk soil (
soil). An optimal design of the coiled probe would minimize the contribution of the probe to the dielectric measurement (or minimize the value of w), thereby maximizing the sensitivity of the coiled probe measurement to bulk soil water content. The parameter n defines the probe's geometry, and
probe is the dielectric constant of the PVC and epoxy material in which the wire coils are imbedded. The dielectric constant of the soil (
soil) as determined by the conventional two-rod probe is written in terms of the fractional bulk volume of each of the three soil phases (1 -
,
-
, and
, for the solid, gas, and water phase, respectively), according to Dobson et al. (1985)
![]() | (2) |
(cm3 cm-3) and
(cm3 cm-3) denote the soil porosity and volumetric water content, respectively, and
s,
a, and
w are the dielectric constant of the soil solid material, air, and water, respectively, with assumed values of
, and
(Dasberg and Hopmans, 1992). The
s value varies slightly with mineralogical composition of the soil solid material (Yu et al., 1999). For instance, the dielectric constant of quartz can vary between 3.75 and 4.1 (Lide, 1996), whereas an Al silicate has a dielectric constant of 4.8 (Fink, 1978). Also, the presence of organic matter increases the dielectric constant of organic soils to values as high as 5.0. For the mineral soils studied here, an
s value of 3.9 appears to be a good estimation for the investigated mineral soils. The exponent
depends on the geometry of the soil solid phase and the soil's orientation with respect to the applied electric field and must be -1 <
< +1 (Roth et al., 1990).
After substitution of Eq. [2] into [1], the dielectric constant as measured with the coiled TDR probe (
coil) can be written as
![]() | (3) |
The presented mixing model approach is preferred to allow for a meaningful physical interpretation of the calibration results (Roth et al., 1990), rather than the model fitting of an arbitrary empirical functional relationship. Moreover, the application of Eq. [3] inherently corrects for the influence of bulk soil density on the bulk soil dieletric constant. Alternatively, one can simply use a polynomial to substitute for Eq. [3], writing
coil as a function of water content, and fit the data to estimate the regression coefficients as was done in Topp et al. (1980), and later for the field calibration results.
Three different soil materials were used to test the coiled TDR probe in the laboratory. These soils were a Yolo silt clay loam (Eching et al., 1994), a Columbia fine sandy loam (Liu et al., 1998), and a washed sand (SRI supreme sand-30, Silica Resources, Marysville, CA). Samples with different water content values (range of 0.00.35 cm3 cm-3) were obtained after mixing a known amount of water to a fixed amount of dry soil. Wetted soils were packed in brass cores (8.25-cm-i.d., 9 cm high) at approximately constant dry bulk densities (Table 1) by packing a predetermined mass of dry soil into the known core volume. Subsequently, samples were covered to prevent water loss by evaporation and put aside for at least 1 d to allow for water distribution before TDR measurements were taken.
|
Calibration curves of
soil vs.
gravimetric were obtained with the mixing model described by Eq. [1], [2], and [3] in two steps using the fitting-model software (Wraith and Or, 1998). First, from TDR measurements of the conventional (
soil) and the coiled (
coil) probe for all three soils together and across the whole water content range, values for w, n, and
probe were fitted to Eq. [1] using the fitting-model software. Subsequently, Eq. [2] was fitted for each soil type to estimate soil-specific
values using independent values for porosity (as estimated from the soil core density) and volumetric water content values for each soil core.
Fitted n, w, and
probe and specific
values for Columbia soil, Yolo soil, and sand were used in Eq. [3] to produce the soil-specific calibration curves using average values of bulk density and porosity (considering a soil particle density of 2.6 g cm-3).
Theory of Dynamic Penetrometer Resistance
The cone penetrometer as used in this study (Fig. 2) is classified as an impact-loading or hammer penetrometer, yielding dynamic penetration characteristics (Bradford, 1986). This type of penetrometer was selected because of its simplicity and ease of construction. Moreover, because soil penetration occurs through several impacts, there is time to measure the water content by the TDR between impacts (
1 min).
During the impact-loading test of the cone penetrometer, the energy stored in the weight at a known elevation is used to drive the penetrometer rod into the soil. The depth of penetration achieved by application of the constant amount of energy is used as a measure of soil PR. The PR can be determined considering that the potential energy of the impact body is converted into work of cone penetration. The total potential energy of the system after impact is equal to the energy of the impact body at height h (m) plus the potential energy of an additional penetration distance x (m). After consideration of the loss of energy due to the impact and the inelastic collision of the weight with the stopper (Fig. 2), we can write (Terzaghi and Peck, 1978; Stolf, 1991)
![]() | (4) |
![]() | (5) |
![]() | (6) |
![]() | (7) |
The PR is obtained from combination of Eq. [7] and [4], which yields after division of the penetration force, F, by the base area of the cone, A (m2),
![]() | (8) |
The characteristics of the penetrometer used in this study (Fig. 2) are
, and
. Substituting these data in Eq. [8], we obtain the equation used to determine the PR (MPa) related to the penetration distance x (m).
![]() | (9) |
Standard penetrometers have a rod diameter slightly smaller than the cone base diameter, to avoid friction between the penetrometer rod and the soil material. However, the combined penetrometermoisture probe has a rod diameter equal to the cone diameter, to avoid air gaps and to ensure good contact between the coiled probe and soil. Therefore, using Eq. [9] may overestimate penetration resistance, when compared with standard cone penetrometer measurements, since we neglect friction losses by the coiled TDR section of the penetrometer. Hence, additional work will be needed to quantify the friction effect on PR measurements with the combined probe.
Field Measurements
In addition to the laboratory tests, the coiled penetrometer probe was tested for the Yolo soil at the Campbell Tract experimental field of the University of California at Davis. The soil is a Yolo silty clay loam with an approximate clay content of 21%, with its texture approximately uniform within the top 60 cm. Measurements of PR and water content were carried out to the 60-cm soil depth. After each impact, penetration depth and water content, as calculated by the WinTDR98 software (Soil Physics Group, Utah State University, 1998) were recorded. Time between impacts was
1 min. After completion of the coilpenetrometer probe measurements, core samples were taken in 5-cm increments from the soil surface to the 60-cm depth using aluminum rings (5-cm diam. and 5-cm height) for subsequent bulk soil density and gravimetric water content determinations in the laboratory. To improve data interpretation, we will present the penetration resistance and the water content profile data combined in a single graph, using depth-average values along 5-cm depth intervals.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
coil) and the conventional probe (
soil), using the data of all three soils is shown in Fig. 4
. The apparent outliers for the sandy material for
coil values larger than 4.5 may be caused by inadequate probesoil contact of the coiled probe for that soil. Using model-fitting software (Wraith and Or, 1998), data were fitted to Eq. [1] yielding parameter values of
. A w value of 0.655 indicates the large influence of the probe material on the dielectric measurement for the coiled TDR probe. Possibly, the geometry of the coiled probe can be changed (wire thickness and spacing and epoxy thickness) to reduce this w value, thereby increasing the sensitivity of the coiled TDR probe. The fitted
probe value of 2.703 is sufficiently close to handbook values of PVC or epoxy of 3.3 and 3.6, respectively (Weast, 1982), thereby validating the parameter fitting approach to some extend. For their specific TDR probe, Nissen et al. (1998) determined a w value of 0.52 and an n value of -0.13, and their fitted relationship is included in Fig. 4 for comparative purposes. Differences in w and n parameters are caused by specific probe characteristics. As expected, there is a nonlinear relationship between the coiled and the conventional probe caused by the constant contribution of the probe material (see Eq. [1]) to the coiled probe (Nissen et al., 1998).
|
was calculated by fitting Eq. [2] to the experimental data obtained with the conventional probe (Fig. 5)
using the model-fitting software and core specific values of porosity and water content. Fitted
values obtained were 0.538, 0.554, and 0.320 for the Columbia soil, the Yolo soil, and the sand, respectively (Table 1). The
values found for the Columbia and Yolo soil were relatively close to reported values of about 0.5 for various soils (Dobson et al., 1995; Dasberg and Hopmans, 1992; Roth et al., 1992; Panizovsky et al., 1999). However, the low
value found for the sand can be attributed to inadequate soilprobe contact of the two-rod conventional probe as well, thereby causing deviations from the generally accepted Topp equation (Walley, 1993), which is also presented in Fig. 5.
|
|
(for each soil) in Eq. [3] and using average values of bulk density and porosity (considering soil particle density equal 2.6 g cm-3) presented in Table 1. Alternatively, as indicated earlier, rather than using the physically based approach presented, the data of Fig. 6 can be fitted directly to an empirical model, thereby circumventing the need for bulk soil dielectric measurements by the conventional two-prong TDR probe. The results in Fig. 5 and 6 show that the dielectric constant values for all soils are close in the water content range of 0.0 to 0.10 cm3 cm-3 for both conventional and coiled probes. However, as the water content further increases, the bulk dielectric constant of the sand is increasingly lower than for the other two soils, despite its higher bulk density compared with the Yolo soil. In a glance, quite the opposite would be expected, since a larger soil density (Table 1) will increase the bulk dielectric constant (Dirksen and Dasberg, 1993), whereas the presence of possible bound water in the finer-textured soil (Columbia and Yolo) may decrease the bulk soil dielectric constant in those soils (Dasberg and Hopmans, 1992). There are two aspects for consideration to better understand the calibration results of both probes (Fig. 5 and 6). The first aspect is the displacement of the sandy soil material around the probe during probe insertion. It was visually observed that the sandy material was displaced at the soil sample surface. The soil displacement created air spaces between the probe and the surrounding sand material, thereby resulting in low dielectric values (Knight et al., 1997) for both probes (Fig. 5 and 6), relative to the other two soil types. Displacement of the sand was more apparent for the coiled than the conventional TDR probe. The second aspect to be considered is the apparent soil compaction in the radial direction by probe insertion for the Columbia and Yolo soil. The compaction probably caused a slight increase in soil bulk density and water content in the immediate vicinity of the TDR probes. Compaction is expected to be more significant for the coiled probe because of its larger volume. According to Roth et al. (1997), who quantified the compaction effect using x-ray computed tomography, soil compaction caused by probe installation can increase the soil bulk dielectric constant, depending on soil type, rod diameter, and soil water content at the moment of probe installation. Although soil density changes by compaction appear significant only within a relatively small soil volume around the probe, it is expected that the travel time of the electromagnetic wave will be affected, since measurement volumes of the TDR signal are small as well (Nissen et al., 1998). Also, we conducted laboratory experiments to directly measure soil density changes in soil cores, from soil sampling around the TDR probe. However, although soil compaction was observed, the large measurement error of soil density due to small sample volumes deemed these results to be uncertain.
Hence, we postulate that differences in
values and calibration curves (Fig. 6) between soils are caused by compaction (Columbia and Yolo) and displacement (sand) of soil material in the immediate vicinity of the TDR probe during probe insertion. Possible compaction of the Columbia and Yolo soil near the probe will cause an increment in the soil bulk dielectric as measured by TDR, whereas displacement of the sandy material by the probe creates air spaces between the TDR probe and surrounding soil, thereby decreasing the soil bulk dielectric constant (Ferre and Rudolph, 1996).
Field Measurements
Volumetric water content and PR measurements with the combined probe were conducted in a bare field research plot with the Yolo soil. Using the combined penetrometerTDR coiled probe (Fig. 2) in the field required fast data acquisition and processing in order to accurately measure the water content and the penetration resistance simultaneously. For that reason we used the WinTDR software for water content determination. However, the software could not include the mixing model as the calibration curve. Instead, we used a second-order polynomial equation
to fit the field calibration data of the Yolo soil (see also Fig. 6).
Figure 7 shows the results of the field measurements using the combined probe, including both the water content and PR data as a function of soil depth in a single graph. The PR and water content values were determined after each hammer impact and were averaged for each 5-cm depth increment. For example, the high PR for the 15- to 20-cm depth increment required five hammer drops, and hence the water content data were averaged for each 5-cm depth increment, enabling comparison with gravimetric water content data (Fig. 8) .
|
|
, indicating that the laboratory calibration can be used satisfactorily for the field measurements.
|
| CONCLUSIONS |
|---|
|
|
|---|
Although it is shown that the concept of the combined penetrometercoiled TDR probe is valid, it is recommended that alternative TDR designs are considered to increase TDR probe sensitivity to bulk soil water content, while simultaneously conducting additional tests to better define the size of the measurement volume of the coiled TDR probe. Moreover, additional field testing may be needed to determine whether soil-specific calibration is needed or that a single calibration may be used for a range of field soils. Finally, we conclude that the combined penetrometerTDR moisture probe can be an excellent tool to investigate the water content dependence of soil resistance in field soils.
| ACKNOWLEDGMENTS |
|---|
Received for publication September 24, 1999.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Minasny and A. B. McBratney Comments on "Simultaneous Measurement of Soil Penetration Resistance and Water Content with a Combined Penetrometer-TDR Moisture Probe" and "A Dynamic Cone Penetrometer for Measuring Soil Penetration Resistance" Soil Sci. Soc. Am. J., May 6, 2005; 69(3): 925 - 926. [Full Text] [PDF] |
||||
![]() |
R. Stolf, K. Reichardt, and C. M. P. Vaz Response to "Comments on 'Simultaneous Measurement of Soil Penetration Resistance and Water Content with a Combined Penetrometer-TDR Moisture Probe' and 'A Dynamic Cone Penetrometer for Measuring Soil Penetration Resistance'" Soil Sci. Soc. Am. J., May 6, 2005; 69(3): 927 - 929. [Full Text] [PDF] |
||||
![]() |
Z. Lu, C. J. Hickey, and J. M. Sabatier Effects of Compaction on the Acoustic Velocity in Soils Soil Sci. Soc. Am. J., January 1, 2004; 68(1): 7 - 16. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Robinson, D. A. Robinson, S. B. Jones, J. M. Wraith, D. Or, and S. P. Friedman A Review of Advances in Dielectric and Electrical Conductivity Measurement in Soils Using Time Domain Reflectometry Vadose Zone J., November 1, 2003; 2(4): 444 - 475. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. C. Topp, D. R. Lapen, M. J. Edwards, and G. D. Young Laboratory Calibration, In-Field Validation and Use of a Soil Penetrometer Measuring Cone Resistance and Water Content Vadose Zone J., November 1, 2003; 2(4): 633 - 641. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Persson and S. Haridy Estimating Water Content from Electrical Conductivity Measurements with Short Time-Domain Reflectometry Probes Soil Sci. Soc. Am. J., March 1, 2003; 67(2): 478 - 482. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Persson and J. M. Wraith Shaft-Mounted Time Domain Reflectometry Probe for Water Content and Electrical Conductivity Measurements Vadose Zone J., November 1, 2002; 1(2): 316 - 319. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. M. P. Vaz, J. W. Hopmans, A. Macedo, L. H. Bassoi, and D. Wildenschild Soil Water Retention Measurements Using a Combined Tensiometer-Coiled Time Domain Reflectometry Probe Soil Sci. Soc. Am. J., November 1, 2002; 66(6): 1752 - 1759. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals |