SSSAJ Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Nichol, C.
Right arrow Articles by Smith, L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Nichol, C.
Right arrow Articles by Smith, L.
Agricola
Right arrow Articles by Nichol, C.
Right arrow Articles by Smith, L.
Related Collections
Right arrow Soil Physics
Right arrow Vadose Zone Processes and Chemical Transport
Right arrow Water Content
Soil Science Society of America Journal 66:1454-1465 (2002)
© 2002 Soil Science Society of America

DIVISION S-1—SOIL PHYSICS

Evaluation of Uncoated and Coated Time Domain Reflectometry Probes for High Electrical Conductivity Systems

Craig Nichol*, Roger Beckie and Leslie Smith

Dep. of Earth and Ocean Sciences, Univ. of British Columbia, 6339 Stores Road, Vancouver, BC, Canada V6T 1Z4

* Corresponding author (cfnichol{at}hotmail.com)


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
High sample electrical conductivity reduces the quality of a time domain reflectrometry (TDR) waveform by the loss of signal amplitude. Two strategies are examined to obtain higher signal/noise waveforms: (i) waveform differencing by remote diode shorting and (ii) covering probe conductors with resistive coatings. Experiments using electrically conductive water solutions (0–5 dS m-1) and three-rod Zegelin type probes show conventional dual-tangent waveform analysis and waveform differencing using manual short circuits can accurately determine travel time but the remote diode shorting method can be systematically biased by the electrical properties of the diodes used. Three-rod Zegelin-type probes with a high resistance coating on the central rod permit collection of analyzable waveforms for solutions with electrical conductivities at least as high as 70 dS m-1. Dual-tangent analysis of the raw waveform is found to be more accurate than the remote diode shorting method within water solutions and within silica sand saturated with an electrically conductive water solution. The probe coating creates a nonlinear relationship between the apparent dielectric permittivity estimated using a coated probe and the actual sample apparent dielectric permittivity. Experimental measurements of this relationship can be fitted using an equation of the form for a coaxial cell. A three-rod coated probe with a single diode at the probe head is a practical means to collect interpretable waveforms in media with high electrical conductivity. However, measurements of travel time alone may not be sufficient to determine water content in soils with high concentrations of dissolved ions in the soil water solution.

Abbreviations: DC, direct current • MTDR, MoisturePoint TDR instrument • PIN, positive-intrinsic-negative diodes • TDR, time domain reflectometry • TEM, transverse electric and magnetic mode • TTDR, Tektronix TDR instrument


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
SOIL WATER CONTENT can be measured in the field using gravimetric methods, neutron scattering, or techniques based upon the thermal or electrical properties of soil-air-water mixtures. Electrical methods are advantageous because they are easily automated to conduct unattended measurements with high sampling rates at multiple locations. Time domain reflectrometry is a commonly used electrical technique. With knowledge of the electrical properties of soil-air-water mixtures, and appropriate experimental calibrations, the bulk soil dielectric permittivity can be related to the volumetric water content and the bulk soil electrical conductivity can be related to the soil water electrical conductance.

While TDR methods have gained wide acceptance and usage in soils, relatively less attention has been given to TDR's application to high electrical conductivity materials such as mining waste. In these materials, conventional TDR probe designs and signal analysis techniques fail to provide sufficiently accurate estimates of the sample dielectric properties. This paper describes an approach to design an automated TDR system for the measurement of water content in unsaturated mine waste rock. The TDR system is used in both a laboratory column experiment and a field-scale experiment (Nichol et al., 2000). In our application, a high soil water electrical conductivity (5–20 dS m-1) caused by oxidizing sulfide minerals in the mine waste leads to poor signal quality using conventional TDR techniques. Sufficient signal amplitude could not be reliably obtained using a short probe of conventional design. Two techniques for obtaining improved signals were selected: (i) waveform differencing using remote diode shorting (Hook and Livingstone, 1992), and (ii) the use of a resistive probe coating (Kelly et al., 1995; Ferré et al., 1996; Mojid et al., 1998). Our aim is to derive a method to reliably collect waveforms that can be interpreted by automated waveform fitting to derive an estimate of travel time and hence apparent dielectric permittivity.

In this paper, we test both techniques for the collection of interpretable waveforms in high electrical conductivity media. We first provide a brief summary of the TDR method and discuss how it is affected by bulk soil electrical conductivity. We then describe the relevant theory of waveform differencing by remote diode shorting and the use of probes covered with a high-resistance coating. A series of experiments are then presented. The first examines the performance of a short (160 mm) uncoated probe using water solutions of variable electrical conductivity. The second experiments analyse the performance of the waveform differencing method using the remote diode shorting method and manual short circuits under the same conditions as the first experiments. We then examine the performance of a coated probe under the same conditions for both conventional TDR waveform collection and analysis, and using the remote-diode shorting method. We next present the performance of a coated probe in silica sand saturated with water solutions of varying electrical conductivity. A last experiment demonstrates a practical means of compensating for the effect of a probe coating on the estimation of sample dielectric permittivity.


    THEORY
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Time domain reflectometry can be used to determine the relative apparent bulk dielectric permittivity ({epsilon}app), the bulk electrical conductivity ({sigma}DC), and in certain circumstances, the frequency-dependent real and imaginary parts of the complex sample dielectric permittivity (Hoekstra and Delaney, 1974; Dalton et al., 1984; Topp et al., 1988; Heimovaara et al., 1996; Friel and Or, 1999). Time domain reflectometry measures the propagation of a fast rise time, step voltage pulse through a coaxial cable to a waveguide (probe) in contact with the sample. Part of the incident pulse energy reflects back to the TDR instrument from each impedance change in the transmission line and probe, with a full reflection of the remaining pulse energy at the end of the probe. The sum of the incident voltage and reflected pulse voltage measured at the TDR instrument is plotted in time and presented as a waveform. The arrival of reflected energy at the TDR instrument causes a change in voltage over time, and thus a deflection of the waveform. The apparent velocity ({nu}app) of the pulse along the probe is determined from the probe length (L) and the arrival times of the reflections from the head and end of the probe (to, t1). The relative dielectric permittivity of the medium surrounding the probe is related to the velocity of a electromagnetic wave propagating in the transverse electric and magnetic mode (TEM) along a waveguide by:

[1]
where c is the velocity of light in a vacuum. It is assumed that the magnetic permeability of soils and soil-air-water solutions (µ) equals that of free space (µo) and therefore o/µ) is unity (Topp et al., 1980). In this paper, all dielectric permittivities are expressed relative to the permittivity of free space, and the term relative is taken as understood.

Successful TDR measurement of the sample dielectric permittivity requires collecting a waveform in which the pulse travel time at the probe head and probe end can be accurately determined. In conventional waveform analysis, the timing of a reflection is determined using two lines fitted to the deflection of the waveform caused by the arrival of reflected pulse energy. The first is a tangent line fit to the steepest slope of the rising or falling limb of the deflection. The second line is either a horizontal line fit to the predeflection waveform (flat-tangent method), or a tangent fit to the slope of the predeflection waveform (dual-tangent method) (Hoekstra and Delaney, 1974; Topp et al., 1982). The time of the intercept of the pre and postdeflection lines is the arrival time of the reflected pulse energy. The dual-tangent method is more accurate for high electrical conductivity solutions (Wraith and Or, 1999) and is used as the reference method for the alternate methods we examined.

The dielectric permittivity estimated by tangent-line analysis of a TDR waveform is termed the apparent dielectric permittivity ({epsilon}app) because the estimate of dielectric permittivity is not obtained for a single frequency. Frequency domain analyses of TDR waveforms using Fourier transform techniques (Hoekstra and Delaney, 1974; Heimovaara et al., 1996; Friel and Or, 1999) demonstrate that the incident TDR pulse generated by common TDR instruments contains energy in a range of frequencies up to 1.5 GHz. The deflection of the waveform is the integration of the arrivals of all of the frequencies reflected. Tangent-line methods of waveform analysis are biased towards the higher frequencies as these create the sharpest waveform deflections. The frequencies which dominate travel time estimates from tangent line waveform analysis are in the range of 700 MHz to 1.0 GHz (Hoekstra and Delaney, 1974; Heimovaara et al., 1996; Friel and Or, 1999).

Effect of Solution Electrical Conductivity
Ions in solution affect both the waveform quality and the pulse velocity. Direct current (DC) conductance and ohmic losses because of current between the probe conductors dissipate the pulse energy within the waveguide. The energy reflecting from the end of the probe is reduced, the waveform deflection becomes less distinct, contributing to greater uncertainty in tangent-line fitting and travel time determination. In samples with high electrical conductivity, the amplitude of the signal reflected from the probe end may be fully dissipated within the exposed probe conductors and no travel time can be determined. The effect of electrical conductivity on the signal amplitude decreases with a shorter length of probe. The practical lower limit for probe length with common TDR instruments is 0.1 to 0.15 m (Kelly et al., 1995).

Dissolved ions directly decrease the dielectric permittivity of a water solution by reducing the number of highly polarizable water molecules per unit volume and indirectly by the interaction of ions with the electrostatic bonding structure of liquid phase water molecules (Hasted, 1973). In contrast, the DC conductance resulting from ions in solution contributes to the imaginary component of the complex valued dielectric permittivity, leading to an increase in dielectric permittivity. However, at the frequencies which contribute most to tangent-line analysis of a TDR waveform, the DC conductance of water has negligible effect on the pulse velocity (Or and Wraith, 1999). Tabulated experimental measurements of {epsilon}({omega})solution for common electrolytes are available in Hasted (1973) and Robinson and Stokes (1959).

The dielectric properties of the water phase in soil differs from that of water in a water solution. Both polar water molecules and dissolved ions interact with mineral surfaces. The dielectric and conductive properties of these diffuse, bonded layers are different from the rest of the water solution (Or and Wraith, 1999). Experiments in saturated and unsaturated soils have shown that increasing the concentration of dissolved ions in soil water increases the measured travel time when soil water content is constant (Wyseure et al., 1997; Sun et al., 2000; Topp et al., 2000). The measured increases in travel time are larger than can be explained by the direct effect of the increased bulk electrical conductivity on the real and imaginary parts of the complex dielectric permittivity (Topp et al., 2000). This indicates that increasing the concentration of dissolved ions in the soil water solution affects other components of the complex-valued dielectric permittivity of the soil/air/water/bound water/dissolved ions mixture in a manner that is not fully understood.

Waveform Differencing and Remote Diode Shorting
Hook et al. (1992) describe a method to improve TDR signal quality using waveform differencing. A raw waveform is first collected from the probe. A short circuit then is created between the signal and ground conductors at the probe head. A second waveform is then collected, termed the shorted waveform. The incident pulse energy will leave the signal conductor and return to the TDR instrument along the ground conductor. This causes a sharp downwards deflection in the shorted waveform. The difference waveform is the difference between the raw waveform and shorted waveform. The timing of the deflection in the difference waveform is determined by tangent-line fitting to derive the travel time for the pulse through the short circuit and back to the instrument. The same process is then repeated with a short circuit created at the probe end.

Positive-intrinsic-negative diodes (PIN) can be used to create the short circuits remotely, and thus allow the waveform differencing method to be automated. PIN diodes have a low resistance under a forward-bias DC current and a high resistance under no current or a reverse-bias DC current. Hook et al. (1992) do not provide details or specifications for the diodes used. We used MPN3404 diodes as suggested by the manufacturers of the MoisturePoint TDR system (Environmental Sensors, Victoria, BC, Canada) which incorporates the remote diode method (Young,1998a). The impedance of these diodes, measured at 100 MHz, switches from >5 M{Omega} to <0.8 {Omega} under a forward bias current (On Semiconductor, 2000).

Our experimental waste rock pile allowed for direct burial of the TDR probes on an open soil surface. This removed the requirement for an open-ended probe design that can be inserted into in situ soil materials. The probe construction is based upon a three-rod Zegelin-type probe design (Young, 1998a), with the inclusion of a single PIN diode at the probe head, and two PIN diodes at the probe base. A schematic of the probes is shown in Fig. 1 . The velocity of the TDR pulse within the sample is determined from the difference between the time for the TDR pulse to pass through the probe-head diode, and the probe-end diode. Calibration measurements in materials of known dielectric permittivity allow the measured time to be corrected for the distance between the diodes located within the probe head and end termination structures and the start of the exposed conductors (Nichol et al., 2002).



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 1. Schematic of TDR probe design: (1) F660 BEF Commscope 75 ohm RG6 coaxial cable; (2) female F-connector; (3) male F-Connector to plug; (4) polyethylene terminal block; (5) 3.2-mm diam. 316 stainless steel rod; (6) 3.2-mm diameter 316 stainless steel rod, optionally coated with 0.4-mm thick polyolefin heat shrink; (7) On Semiconductor MPN3404 PIN diode; and (8) terminal block sealed with a silicon sealant.

 
Resistive Probe Coatings
The ohmic loss of signal voltage between the signal conductor and the ground conductors can be reduced by introducing a high resistance coating (Kelly et al., 1995). In this work, we utilize a polyolefin heat shrink resistive coating on the center rod (Livingstone, personal communication, 1997; Mojid et al., 1998). The energy of the travelling TEM wave will still extend outside the coating. The inclusion of an additional dielectric material complicates the relationship between the measured electrical properties and the electrical properties of the sample (Annan, 1977; Knight et al., 1997). The polyolefin heat shrink used has an apparent dielectric permittivity of ~3, lower than most soils materials, and therefore the estimate of apparent dielectric permittivity using a coated probe ({epsilon}cp) will always be less than the apparent dielectric permittivity that would have been estimated using an uncoated probe ({epsilon}app).

Three-rod or multi-rod probes are designed to emulate a coaxial cell (Zegelin et al., 1989). A simple analytical solution for the effect of a dielectric material of uniform thickness around the central conductor of a coaxial cell can be derived from Annan (1977).

[2]
where Kcp, equals the apparent dielectric determined from travel time in coated probe, Kc is the dielectric constant of coating, Ks represents the dielectric constant of sample; rc equals the outer radius of coating, ri is the outer radius of the inner conductor, and ro represents the radius to the outer conductor.

Equations have been derived which exactly describe the effects of coatings applied to three-rod or multi-rod probes, but these equations do not have analytical solutions and are solved numerically (Zegelin et al., 1989; Knight et al., 1994, 1997). We therefore test the simpler analytical solution.

All probe designs with a fixed dielectric material within the sampling zone have a nonlinear relationship between the apparent dielectric permittivity estimated directly from the measured travel time and the actual sample dielectric permittivity. This nonlinear relationship means that the sample dielectric permittivity must be uniform along the probe length to properly determine the sample dielectric permittivity (Ferré et al., 1996).


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
The primary objective of this research is to obtain reliable TDR signals and therefore estimates of the apparent dielectric permittivity in samples with high electrical conductivity. We conducted a series of experiments to compare the results of traditional methods of TDR signal collection and analysis with alternate methods that have been proposed for increasing signal quality. The probes used for the experiments (Fig. 1) are unbalanced, Zegelin-type three-conductor probes constructed using 3.2-mm diam. 316 stainless steel rods of varying length. Coated probes include a polyolefin heat shrink, carefully applied to the center conductor to ensure no air gaps between rod and coating. The average thickness was 0.40 mm. Probe designs are quoted in the remainder of the paper as three numbers a/b/c where a is the length of exposed rods in millimeters, b is the rod diameter in millimeters, and c is the center rod to outer rod spacing in millimeters.

Reference raw waveforms with 5000 points in time were collected using a Tektronix 1502C instrument (TTDR) connected to a PC running WAT TDR Version 3.11 (Waterloo Center for Groundwater Research, Waterloo, ON, Canada). All waveform analysis of TTDR waveforms was performed manually using the dual-tangent analysis method. All waveforms are presented as the reflection coefficient displayed against the two-way travel time. This approach represents the currently accepted practice of TDR measurement against which we compare the alternate methods.

Raw waveforms and diode-shorted waveforms with 256 points in time were collected using a MoisturePoint TDR (MTDR) instrument (Environmental Sensors, Victoria, BC, Canada) with firmware Version 1.27. The MTDR instrument was purchased with proprietary software and hardware specifically designed for conducting automated, diode-shorted TDR measurements. Travel times were measured using the automated waveform fitting routines which are included in the instrument firmware. Additional waveforms for presentation and for detailed waveform analysis were collected manually using a PC connected to the MTDR and running the instrument control software ViewPoint Version 1.34. For all measurements, all instrument settings and waveform fitting parameters in the instrument firmware and software were optimized for maximum accuracy and stability of tangent-line fitting. The MTDR instrument warm up time was set to 15 s and a minimum of five measurements were discarded prior to all recorded measurements to ensure the instrument returned a stable value of travel time.

The first experiments determine the effect of solution electrical conductivity on conventional travel time analysis and on the performance of the remote diode shorting method. Waveform and travel time analysis was carried out for a 160/3.2/25 uncoated probe immersed in isothermal water of variable solution electrical conductivity. The probe was connected to the TTDR and MTDR using 6 m of Commscope F660BEF 75 {Omega} RG6 coaxial cable, a male F-connection termination, and an F-to-BNC adapter. Solution electrical conductivity was altered from 0 to 5 dS m-1 using KCl. Raw waveforms and probe head and end diode-shorted waveforms were manually collected for each solution.

The effectiveness of the MPN3404 diodes as short circuits was assessed by taking waveform measurements using manually created short circuits. For distilled water and for a solution electrical conductivity of 2 dS m-1, 6 gauge Cu wire was used to create a short circuit at the start of the exposed rods and then at the end of the exposed rods. Shorted waveforms of each short circuit were collected using both TTDR and MTDR.

The second experiments address the effect of a resistive coating. Waveform and travel time analysis was carried out for a 281/3.2/25 coated probe using the same experimental setup as the uncoated probe. Electrical conductivity was altered from 1 to 70 dS m-1 using NaCl. The measurements for 30 to 70.8 dS m-1 solutions have some conductive amplitude loss of signal, and different overall probe impedance because of failure of the probe-head seal and invasion of electrically conductive water into the probe head. Conclusions can still be drawn from these results and they are included in the following analysis. No manually created short circuits were tested in these probes because of the presence of the probe coating.

The performance of a 281/3.2/25 coated probe and remote diode method was then assessed for silica sand saturated with electrically conductive water solutions. Water electrical conductivity was altered from 0 to 20 dS m-1 using NaCl. Raw and diode shorted waveforms were collected using the MTDR. Travel time was determined using automated remote diode shorting at the probe head and end. Water solution and bulk soil electrical conductivity were measured independently using a four electrode conductivity cell and conductivity meter.

The effect of the resistive coating on the estimation of sample dielectric permittivity was determined in the final set of experiments. Travel time measurements were made using a coated probe and materials of known dielectric permittivity with low electrical conductivity, the latter determined using a calibrated 295/3.2/25 uncoated probe. Measurements were conducted using the MTDR instrument and averaging a minimum of 15 measurements. For the purposes of these experiments, the best dielectric materials are liquids, which allow simple probe immersion, uniformity of sample and a controlled range of dielectric properties. Organic liquids can be used to obtain a range of dielectric properties. However, the 281/3.2/25 coated probe used in this study was required for other calibration work and the compatibility of the probe construction materials with various organic liquids was not known. To prevent damage to this probe, a more limited range of inert materials was used. Air, oven-dried silica sand, wetted silica sand, methanol, and water were used as dielectric media. Liquids were placed in a 15.2-cm diam. polyvinyl chloride (PVC) container. Granular materials were packed within a 12 by 12 by 63 cm watertight container.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Uncoated Probes
The conventional method of TDR raw waveform analysis is first demonstrated for an uncoated probe in water. Raw waveforms are presented in Fig. 2 for solution electrical conductivities from 0 to 5 dS m-1. The reflection from the probe head is virtually identical for all solution electrical conductivities. The negative slope on the probe segment of the waveform decreases and the probe end reflection shows a decrease in the final amplitude as solution electrical conductivity increases. Manually determined tangent-lines to the end reflections are presented for 0 and 2 dS m-1. The signal quality is significantly degraded at 5.0 dS m-1 making manual waveform fitting challenging. In the range of 0 to 5.0 dS m-1, no changes in travel time greater than the method uncertainty were detected. By definition, the transition time of a reflection is half the time for the signal to rise from the prereflection tangent to the maximum amplitude and is related to the highest frequencies contained within the reflected signal (Hook and Livingstone, 1992). The transition time of the probe end reflection was consistently 1.2 ns from 0 to 5 dS m-1. These water solution electrical conductivities do not contribute to either increased transition time or decreased amplitude of the high frequency components of the TDR pulse.



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 2. Raw time domain reflectrometry waveforms with manually determined tangent lines to the probe end reflection for a 160-mm length uncoated probe.

 
Signals obtained using remote diode shorting are shown in Fig. 3 . The gain- and offset-corrected MTDR raw waveforms, diode-shorted waveforms, and diode-difference waveforms are shown for distilled water (Fig. 3A,B) and a solution electrical conductivity of 2 dS m-1 (Fig. 3C,D). Waveforms are separated for when the diode is located at the probe head (Fig. 3A,C) and at the probe base (Fig. 3B,D). The probe head and probe end diode-shorted waveforms for distilled water both drop rapidly because of the short circuit created by the diodes. The resulting diode-difference waveforms have distinct signals, with sharp waveform rises and short transition times, making tangent-line fitting both simple and accurate. Note that the vertical position of the difference waveforms in this figure, and following figures, has been shifted vertically for greater clarity of presentation. However, at 2 dS m-1, for both the probe head (Fig. 3C) and probe end (Fig. 3D), the diode-shorted waveforms more closely follow the raw waveform. As a result, the diode-difference waveforms demonstrate reduced total amplitude, increased time to maximum amplitude, and a less distinct shape than those in distilled water. The 2 dS m-1 probe head and probe end diode-difference waveforms are more difficult to analyse by tangent-line fitting than the raw waveforms presented in Fig. 2.



View larger version (33K):
[in this window]
[in a new window]
 
Fig. 3. Remote diode shorting method waveforms for diodes located at the probe head (A, B) and at the probe base (C, D) for a 160-mm length uncoated probe in water solutions of varying electrical conductivity.

 
The performance of the MPN3404 diodes can be understood by examining the data for near-perfect short circuits created with Cu wire. Figure 4 presents the waveforms obtained using manually created short circuits in distilled water (Fig. 4A,B) and in a 2 dS m-1 water solution (Fig. 4C,D) for the probe head (Fig. 4A,C) and probe end (Fig. 4B,D). The shorted waveforms at the probe head and end are similar for both solution electrical conductivities. The difference waveforms have different amplitudes because of the difference in amplitude of the raw waveform, but similar transition times. The poor signal quality of the difference waveforms in Fig. 3 is therefore not the result of the waveform differencing method itself, but results from the MPN3404 diodes not creating adequate short circuits.



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 4. Manual probe shorting method waveforms for short circuits created using copper wire located at the probe-head end (A, B) of the exposed conductors and at the probe base end (C, D) of the exposed conductors for a 160 mm length uncoated probe in water solutions of varying electrical conductivity.

 
At the probe head, the multifrequency incident TDR pulse encounters two possible transmission paths when a short circuit is present. Wave energy will be reflected and transmitted along the waveguide or into the short circuit in proportions determined by the frequency dependent complex impedance of each path. The impedance of the manual short circuit is near zero, and thus all wave energy diverts from the signal conductor to the ground conductor. After 85 ns, the manually shorted waveforms in Fig. 4A and C are both are flat, illustrating that there are no reflections present in the shorted waveform from any energy that passed the short circuit and entered the probe.

In contrast, the pathway created by the remote diode short has a low, but finite impedance. Pulse energy does divide between the diode pathway and the probe pathway when it encounters the two impedances in parallel. The diode-shorted waveform collected at 2 dS m-1 at the probe head closely follows the raw waveform (Fig. 3C), indicating that significant signal energy continued past the location of the diode in the probe head, entered the probe and returned reflections from impedance changes within the probe. The transition time of the probe end reflections in the diode-shorted and the raw waveforms indicates that the frequency spectra of the energy pulse transmitted into the probe was similar. The difference between the diode-shorted waveform and the raw waveform represents the energy that returned via the diode. The gradual departure of the diode-shorted waveform from the raw waveform indicates the energy passing through the diode was dominated by lower frequencies. The probe impedance at low frequency (Z{infty}) was calculated from the reflection coefficient at long times using the method of Mallants et al. (1996). The probe impedance is 4000 {Omega} in distilled water and 20 {Omega} in 2 dS m-1 water. The MPN3404 PIN diodes used have a rated impedance of 0.8 {Omega} at 100 MHz. It is clear from Fig. 3 that the MPN3404 diode and our probe impedances are more similar at higher frequencies, and thus remote diode shorting becomes less effective.

At the probe end, a similar division of pulse energy occurs as at the probe head. Connecting diodes across the end of the probe conductors requires the presence of a probe termination structure, similar to the probe head, to house the diodes. Energy is reflected and transmitted at both the end of the exposed waveguides, and again at the diodes, which connect the end of the center conductor to ground within the probe end termination structure. At 2 dSm-1 (Fig. 3D), the higher final amplitude of the diode-shorted waveform than in distilled water (Fig. 3B) indicates a greater proportion of the energy was reflected back along the signal conductor, instead of through the diode short.

In high conductivity systems, the combination of the effects of imperfect diode shorting at the probe head and end leads to large errors in travel-time estimation, and hence water content. The automated instrument-determined travel time for varying solution conductivity is shown in Fig. 5 , with the theoretical travel time calculated from Eq. [1], and the travel time estimated from the manually created shorts in Fig. 4. The travel time estimated using manually created shorts does not include travel time within the probe head and end termination structures from the diodes to the ends of the exposed conductors. The automated determination of travel time increased by 6 ns. Remote-diode shorting hinders travel time estimation under the conditions examined. In contrast, manually fitted tangent-line analysis of the manual short circuit data presented in Fig. 4 determined no change in travel time greater than the uncertainty in the analysis, similar to the results of manual tangent-line fitting to the raw waveforms in Fig. 2.



View larger version (16K):
[in this window]
[in a new window]
 
Fig. 5. Automated remote diode shorted method measured two-way travel time, travel time measured using manual probe shorting and waveform differencing and travel time calculated from Hasted (1973) data for a 160-mm length uncoated probe in water solutions of varying electrical conductivity.

 
The method of waveform differencing does not significantly extend the range of waveform analysis to media with higher conductivities. At the probe end, waveform differencing by either manual short circuits, or diode short circuits relies upon the collection of a raw waveform with nonzero reflected signal amplitude. This method will therefore only be successful up to similar solution conductivities as conventional raw waveform analysis, or ~5 dS m-1 for a 0.16-m probe.

Coated Probes in Water
The performance of coated probes is first assessed using conventional waveform analysis. Raw waveforms for solution conductivities between 1 to 70.8 dS m-1 collected using the TTDR are shown in Fig. 6 . The raw waveforms from 1 to 20 dS m-1 indicate little signal amplitude loss and a consistent trend in probe impedance. The transition time of the probe head reflections is similar for 1 to 70.8 dS m-1. No attempt was made to precisely match the impedances of the cable, connectors, and probe head. Multiple reflections are seen between 67 and 73 ns caused by these impedance changes and the start of the exposed conductors. The accidental penetration of high electrical conductance water into the probe head causes the sharp drop in probe-head impedance observed in the lowest three waveforms. The multiple internal reflections originating from within the probe base are smoothed out at higher electrical conductivities.



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 6. Raw time domain reflectrometry waveforms with manually determined tangent lines to the probe end reflection for a 281-mm length coated probe in water solutions of varying electrical conductivity.

 
The dual-tangent method of waveform analysis was assessed by manual tangent-line fitting. The better preservation of probe head internal reflections in coated probes does not affect manual tangent-line fitting, but would make automated tangent-line analysis of the probe head reflection more difficult. The signal from the end reflection of a coated probe has lower transition times than the signal from the uncoated probe, more signal amplitude and therefore tangent lines are easier to fit to the waveform. Manually determined tangent lines to the probe end reflection are presented on Fig. 6. Between 1 and 70.8 dS m-1 the travel time determined by the tangent-line method decreases by 0.4 ns. This decrease in travel time will be discussed in final section. The signal attenuation encountered in uncoated probes, caused by high sample conductivity, is reduced in coated probes and raw waveforms are interpretable to >70 dS m-1.

The waveform differencing method was assessed for coated probes using only the remote diode shorting method as no manual short circuits could be created because of the presence of the probe coating. The diode waveforms collected at 1 (Fig. 7A,B) and 20 dS m-1 (Fig. 7C,D) are presented in Fig. 7, with the probe head waveforms on the left (Fig. 7A,C) and the probe end waveforms on the right (Fig. 7B,D). In contrast to the uncoated probes, the overall probe impedance calculated from the final reflection coefficient remains high throughout the range of solution electrical conductivities presented. The probe impedance at low frequency changes from 6000 to 1000 {Omega} from distilled water to 20 dS m-1. The impedance contrast between the shorted diode and the probe remains high and therefore the diode acts as a better short circuit than in uncoated probes. The probe head diode-shorted waveforms indicate negligible changes in the slope and intercept of the diode waveforms between 1 and 20 dS m-1. The probe head diode-difference waveforms eliminate internal reflections within the probe head, and are smoother and more easily interpreted using tangent-line analysis than the raw waveforms in Fig. 6.



View larger version (36K):
[in this window]
[in a new window]
 
Fig. 7. Remote diode shorting method waveforms for diodes located at the probe head (A, B) and at the probe base (C, D) for a 281-mm length coated probe in water solutions of varying electrical conductivity.

 
The diode-shorted waveform for the probe end diode in a 1-dS m-1 solution shows a very small reflection from the conductor to probe end transition; none is visible at 20 dS m-1. The probe end diode-shorted waveform at 20 dS m-1 has a 4-ns longer transition time than the 1 dS m-1 waveform because of frequency-dependent transmission by the diode and frequency filtering by the solution. The automated measurements of travel time using remote diode shorting from 1 to 70 dS m-1 are presented in Fig. 8 . The travel time calculated using the tangent-line method remains constant up to 5 dS m-1. Between 5 and 20 dS m-1, the travel time rises by 0.5 ns, beyond which it decreases. The combined effects of the longer diode waveform transition times and the raw waveform smoothing contribute to the inconsistent estimation of travel time. The decrease after 20 dS m-1 may be the result of the invasion of water into the probe head. The probe coating therefore allows collection of a waveform interpretable by travel time analysis in samples with electrical conductivities that are much higher than uncoated probes. Waveform differencing improves accuracy of travel time determination at the probe head, even with the MPN3404 diodes. We can infer from the success of the waveform differencing method with manual short circuits in uncoated probes that either manual short circuits or perfect remote diode short circuits created would be similarly successful in a coated probe, but the well preserved end reflection removes the necessity for waveform differencing at the probe end.



View larger version (11K):
[in this window]
[in a new window]
 
Fig. 8. Automated remote diode shorted method measured two-way travel time for a 281-mm length coated probe in water solutions of varying electrical conductivity

 
Coated Probes in Saturated Sand
The first experiments presented used electrically conductive water solutions to isolate the effects of DC conductance on signal quality from any other effects arising from the interaction of water and soil particles. We now examine data collected in a saturated silica sand to determine the quality of the waveforms collected by a coated probe using a sample material closer to our expected field conditions. Figure 9 presents waveforms collected in silica sand saturated with distilled water (Fig. 9A,B), and saturated with a water solution with a electrical conductivity of 19.6 dS m-1 (Fig. 9C,D). Waveforms are displayed for the probe head (Fig. 9A,C) and probe end (Fig. 9B,D). When the saturating solution was 19.6 dS m-1, the bulk electrical conductivity of the saturated silica sand measured with a four electrode conductivity probe was 8.3 dS m-1. The diode-difference waveforms at the probe head are good quality, and again, probe-head travel time is more easily estimated from the diode-difference waveform than the raw waveform. The diode-difference waveforms have the same rise time and waveform fitting derives the same travel time. At the probe end, the waveforms in the silica sand saturated with 19.6 dS m-1 water (Fig. 9, Panel D) are significantly more rounded, with longer transition times. The transition time of the probe end reflection is the same as in 20 dS m-1 water (Fig. 7D), but the final amplitude is actually lower.



View larger version (39K):
[in this window]
[in a new window]
 
Fig. 9. Remote diode shorting method waveforms for diodes located at the probe head (A, B) and at the probe base (C, D) for a 281-mm length coated probe in silica sand saturated with water solutions of varying electrical conductivity

 
Figure 10 presents the two-way travel times estimated from the waveforms in Fig. 9 using automated remote diode shorting at both the probe head and end (Curve A, closed circles) and manual dual-tangent fitting to the raw waveforms (Curve B, open circles). The measured travel times for both the remote diode method, and the manual method both increase, but the remote diode method overestimates the increase. Curve C (closed squares) presents the travel time measured in water of the same bulk electrical conductivity, taken from Fig. 8. This indicates that the increase in travel time determined by Curves A and B is not directly the result of increased electrical conductivity of the water solution, but arises from the interaction of the electrically conductive water solution with the mineral particles.



View larger version (16K):
[in this window]
[in a new window]
 
Fig. 10. Measured two-way travel time for 281-mm length coated probe in silica sand saturated with water solutions of varying electrical conductivity (A, B) and water of varying electrical conductivity (C). (A, closed circles) Automated remote diode shorting at both the probe head and end. (B, open circles) Manual dual-tangent fitting to the raw waveforms. (C, closed squares) Measurements in water solutions using automated remote diode shorting at both the probe head and end.

 
The measured apparent dielectric permittivity increases in both the manual and remote diode waveform analysis methods and we can therefore infer this is a real increase in apparent dielectric permittivity, similar to observations by Wyseure et al. (1997), Sun et al. (2000), and Topp et al. (2000). As discussed in the theory section, the apparent dielectric permittivity is not an actual physical property of the sample media, but is an average measurement over a nonspecific range of frequencies defined by the tangent-line method of waveform analysis. An increase in apparent dielectric permittivity may therefore be the result of either changes in the frequency range used to define the apparent dielectric permittivity or actual changes in the frequency-dependent dielectric properties of the media. It is clear from the waveform smoothing in Fig. 7 and 9 that changes in the frequency content of the TDR pulse reflections are occurring with increasing electrical conductivity. Changes in the actual dielectric permittivity of the soil-water mixture may also be occurring.

Detailed investigation of changes in frequency content of the TDR pulse, and changes in the dielectric permittivity of soil/air/water/bound water/dissolved ions mixtures are best determined using frequency domain techniques where it is possible to determine both the real and imaginary parts of the complex dielectric permittivity across all the frequencies in the TDR pulse (Hoekstra and Delaney, 1974; Heimovaara et al., 1996; Friel and Or, 1999). This type of analysis, using both uncoated and coated probes, would increase our understanding of the effects of high concentrations of dissolved ions on dielectric permittivity, TDR signal frequency content and on the apparent dielectric permittivity derived from tangent-line analysis.

These methods are not currently suited to automation in a field situation, and therefore we have focussed our analysis on conventional travel time measurements in this paper. We are primarily concerned with signal quality, and obtaining a waveform interpretable by automated waveform fitting. The results obtained in saturated silica sand indicate that a coated probe is sufficient to determine a measurement of travel time and hence apparent dielectric permittivity, but that apparent dielectric permittivity alone is not sufficient to determine the water content of a soil containing a high electrical conductivity water solution. Calibration of the apparent dielectric permittivity to water content for mine waste material is presented further in Nichol et al., 2002.

Effect of Probe Coating
The apparent dielectric permittivity estimated using a coated probe must be related to the apparent dielectric permittivity of the soil to accurately estimate water contents. The third set of experiments determined the measured travel time for both a coated probe and an uncoated probe in materials with different dielectric permittivity. Apparent dielectric permittivity for coated and uncoated probes were estimated from these measured travel times using Eq. [1]. The apparent dielectric permittivity ({epsilon}cp) derived from a coated probe is plotted against the apparent dielectric permittivity ({epsilon}app) determined with an uncoated probe in Fig. 11 .



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 11. Relationship of coated probe measured apparent dielectric permittivity to apparent dielectric permittivity: linear correction (A); Eq. [2] using known probe dimensions (B): ro = 12.5 mm, Kc = 2.8 (B); Eq. [2] using ro as fitting parameter (C): ro = 53.5 mm., Kc =2.8.

 
The dielectric permittivity of the probe coating can be estimated from where the apparent dielectric permittivity from the coated probe and the actual apparent dielectric permittivity coincide. This is estimated to be at a dielectric permittivity of 2.8. Three fitted curves are shown in Fig. 11. Curve A represents a two-point linear correction for the presence of a coating material using the air and water measurements as suggested by Young (1998b). The maximum overestimation of the apparent dielectric permittivity using Curve A occurs at a dielectric permittivity of 30, which corresponds to the apparent dielectric permittivity of saturated silica sand. The estimate of apparent soil dielectric permittivity using Curve A would be double the true value. Curve B is derived from Eq. [2] using the known probe dimensions and the estimated probe coating dielectric permittivity. Experiments by Friel and Or (1999) demonstrate that three-rod probes may emulate coaxial cells up to frequencies of 0.7 GHz. Equation [2] may be appropriate up to that frequency, but based upon the poor fit between the measured data and Curve B, it is not applicable at the 0.7 to 1 GHz that is estimated to dominate the travel time analysis of TDR waveforms. Curve C is the best-fit equation to the measured data obtained by modifying Eq. [2]. The rod and coating dimensions are retained as the known values, but ro is used as a fitting parameter. Using the form of Eq. [2], the correct curve shape is captured while using a relatively simple fitting equation.

It is now possible to assess the accuracy of a coated probe for the determination of the apparent dielectric permittivity of a high electrical conductivity sample. A 0.4 ns change in the measured travel time was noted in the previous discussion of the performance of coated probes in solutions between 1 and 70 dS m-1 (Fig. 6). The highest solution electrical conductance, 70.8 dS m-1, corresponds to a NaCl concentration of 0.77 mol L-1 (Mojid et al., 1998). Using the data tabulated in Hasted (1973), the dielectric constant would be expected to decrease from 80.7 to 69. Using Eq. [1], the two-way travel time in an uncoated probe of the same length as the coated probe should decrease by 1.3 ns. The actual decrease in measured travel time for the coated probe is affected by the probe coating. Using the best fit curve (Fig. 9C), this calculated change of dielectric permittivity would correspond to a change in measured dielectric permittivity from 29.3 to 27.3 and a travel time difference of 0.38 ns. This corresponds to the measured 0.4 ns change in travel time between 1 and 70 dS m-1. The coated probes are therefore sensitive enough to determine small changes in the dielectric permittivity.

Recommendations
The method of waveform differencing for determination of the travel time is accurate for high electrical conductivity samples (Fig. 4) when the waveguides are perfectly shorted. As this method requires a raw waveform with an end reflection, it has the same upper limit of solution electrical conductivity as does the interpretation of raw waveforms. The remote diode method of creating short circuits using the MPN3404 diodes recommended by Young (1998) is inaccurate for media with high electrical conductivities. The MPN3404 diodes have frequency-dependent impedances that are incompatible with our probe design. The results obtained for the remote diode method may be improved by the use of diodes other than the MPN3404 that have lower impedance characteristics at frequencies of 700 MHz to 1 GHz. We have not investigated alternate diodes, but can recommend that any diodes selected can be easily checked for performance by measuring diode waveforms of an uncoated probe in an electrically conductive water solution. The presence of significant signal amplitude in the diode-shorted waveform is an indication of unsuitable diode characteristics.

Measurement of apparent dielectric permittivity in high electrical conductivity systems requires the use of a probe coating to obtain an interpretable waveform. Our three-rod Zeglin-type probe with a coating on the center rod obtained interpretable waveforms in both water solutions and saturated silica sand. The accuracy of the probe head measurements for coated probes (Fig. 7) demonstrates that probe head travel time for coated probes is more easily determined using remote diode shorting, even with the MPN3404 diode. The accuracy of locating the probe-end reflection using conventional waveform analysis (Fig. 6) indicates that probe base travel time can be determined from the raw waveform using a dual-tangent method. No probe-end diodes or terminal block (Fig. 1) are therefore required. An empirical correction (Fig. 9) for the probe coating should be derived by TDR measurements in liquids of known dielectric properties. The calibration media should have dielectric permittivities close to the range of expected dielectric permittivities in the soil and calibration of individual probes is required. The use of a resistive coating permits the determination of travel time and hence apparent dielectric permittivity in high electrical conductivity systems. Further work using frequency domain techniques is required to investigate the changes in apparent dielectric permittivity with both soil water content and the concentration of dissolved ions in the soil water solution.


    ACKNOWLEDGMENTS
 
This project is part of the Waste Rock Hydrology Research Program, a joint research program initiated between the University of British Columbia, the University of Saskatchewan, Cogema Resources Incorporated, Cameco Corporation, and the Natural Sciences and Engineering Research Council of Canada.

Received for publication May 10, 2001.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 THEORY
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 




This article has been cited by other articles:


Home page
Vadose Zone JHome page
M. Persson, D. Bendz, and P. Flyhammar
Time-Domain Reflectometry Probe for Water Content and Electrical Conductivity Measurements in Saline Porous Media
Vadose Zone J., November 1, 2004; 3(4): 1146 - 1151.
[Abstract] [Full Text] [PDF]


Home page
Vadose Zone JHome page
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]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Nichol, C.
Right arrow Articles by Smith, L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Nichol, C.
Right arrow Articles by Smith, L.
Agricola
Right arrow Articles by Nichol, C.
Right arrow Articles by Smith, L.
Related Collections
Right arrow Soil Physics
Right arrow Vadose Zone Processes and Chemical Transport
Right arrow Water Content