|
|
||||||||
USDA-ARS Soil and Water Management Unit, 1991 Upper Buford Circle, St. Paul, MN 55108
* Corresponding author (dcobos{at}soils.umn.edu).
Received 12 March 2003.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The most common method of measuring soil heat flux is the combination method (Fuchs and Tanner, 1968; Fuchs, 1986) where soil heat flux plates are used to measure the heat flux at a known depth, and the heat storage in the layer above the plates is derived by calorimetry, where storage is the product of the temperature change during a given time period and the specific heat of the soil above the plate. The standard design of a soil heat flux plate is a thermopile sensor embedded in a resin disk of known thermal conductivity. The heat flux by conduction through the soil is described by Fourier's Law
![]() | [1] |
soil is the thermal conductivity of the soil (W m-1 K-1), and dT/dz is the vertical temperature gradient (K m-1). Similarly, heat flux through the soil heat flux plate is the product of the thermal conductivity of the resin material (
plate) and the temperature gradient across the disk. However, unless
plate =
soil, the heat flux field will either converge or diverge at the plate, so that the heat flux through the disk is not an accurate measure of the soil heat flux. This phenomenon, known as deflection error, was described by Philip (1961), who showed that deflection error could be minimized by reducing
plate and decreasing the thickness of the plate in proportion to the surface area of the horizontal plate surface. Philip (1961) also developed an analytical correction for deflection error based on the geometry of the plate, its thermal conductivity, and that of the soil. Since this correction requires additional measurements, it is generally ignored. Recently, a new type of heat flux plate has become commercially available. These plates have a heating element integrated into one side of the sensor. By generating a heat pulse of known power, and measuring the proportion of it that travels through the plate relative to that traveling into the soil, the deflection error can be corrected. However, these heat flux plates are quite large (diam. = 8 cm). Because the heat flux plates are impervious to the flow of fluids in the soil, the transport of both liquid water and water vapor is prevented. This would be expected to cause differences between the thermal properties of the soil above and below the sensor after rainfall or irrigation. Additionally, since the heat flux plates physically block latent heat transport through the soil, condensation on either heat flux plate surface would result in erroneous estimates of soil heat flux by the plate (Fuchs, 1986; Buchan, 1991).
All heat flux plates are also plagued by problems associated with installation. The general installation procedure is to expose a vertical soil face, scrape out a horizontal groove into the face, and insert the heat flux plate into the groove. It is preferable to have the groove slightly thinner than the thickness of the heat flux plate so that the heat flux plate fits tightly into the groove, and is therefore in good thermal contact with the soil on both the upper and lower faces. This method requires some skill if good soilplate contact is to be achieved without undue perturbation of the surrounding soil. Even when great care is taken in installation, good soilplate contact may not be achieved until the soil has been thoroughly wetted and has reconsolidated.
An alternative to the use of soil heat flux plates is direct measurement of soil heat flux with the gradient technique, where the heat flux is computed directly from Fourier's Law (Eq. [1]) as the product of the soil thermal conductivity and the gradient of temperature with depth. This method has been proposed a number of times (e.g., Staley and Gerhardt, 1957; Tanner, 1963; Fuchs, 1986, Sauer, 2002) and used in the laboratory to calibrate heat flux plates (van Loon et al., 1998), but the difficulty of determining
soil accurately has limited the use of this method in the field. Kimball et al. (1976) used the gradient approach to calculate soil heat fluxes with the DeVries (1963) method to compute
soil from the weighted averaging of the thermal conductivities of the soil constituents. Because of uncertainty in the DeVries (1963) approach to calculate
soil, Kimball et al. (1976) recommended that the gradient technique be used at a depth of 20 cm where large uncertainties in
soil have minimal impact on the determination of surface soil heat flux.
Recent advances in the use of transient-heated needle techniques have allowed the direct, accurate, in situ quantification of soil thermal conductivity (Bristow et al., 1994). Use of heated-needle techniques in conjunction with a concurrent, colocated measure of the soil temperature gradient permits a direct calculation of soil heat flux by the gradient technique. With the improved accuracy in the determination of
soil, the gradient technique can be employed at 5- to 10-cm depths where soil heat flux plates are commonly deployed with the combination method. Since the needles used are small in diameter (generally <2 mm), this method is much less intrusive than heat flux plates, and heat and fluid flow in soil are unaltered, preventing both deflection errors and any errors arising from the blockage of liquid water and water vapor. Furthermore, the needle-type sensors can be inserted horizontally into a vertical soil face without excavation of the soil volume directly involved in the measurement, thus minimizing soil disturbance, and promoting good soilsensor contact. We describe here both laboratory and field evaluations of the applicability of this method for measuring soil heat flux.
| MATERIALS AND METHODS |
|---|
|
|
|---|
![]() | [1] |
soil is the thermal conductivity of the soil (W m-1 K-1), and dT/dz is the temperature gradient (K m-1). The temperature gradient can be approximated by measuring the soil temperature at two depths, but
soil is more difficult to quantify. There are two main methods for direct in situ measurement of soil thermal conductivity by transient heating of needle probes; the single- and dual-probe methods. Both methods use an electrically heated needle, which approximates an infinite line source of heat. The rate at which the heat is conducted away from the line source depends on the thermal properties of the soil, so the thermal conductivity can be determined from the temperature response across time, either adjacent to the heat source (single-needle method; e.g., De Vries, 1952; Shiozawa and Campbell, 1990), or at some finite distance from the heat source (dual-needle method; e.g., Lubimova et al., 1961; Campbell et al., 1991). We have chosen here to employ the single-needle method.
With the single-needle transient heating method, a heater and temperature sensor are mounted together in a single needle. The temperature response of the probe is measured during heating, and the thermal conductivity is derived from a simplified analytical solution of the heat conduction equation. The complete solution for the heat conduction equation as it relates to thermal conductivity measurement by transient heating probe methods is well known. For a thorough discussion of the theory behind soil thermal property measurement by both the single and dual needle methods, see Bristow et al. (1994). Briefly, the soil thermal conductivity can be closely approximated with the single-needle method by
![]() | [2] |
Sensors
The sensors used here were custom fabricated by Thermal Logic Devices (Pullman, WA). They consist of three parallel stainless steel needles extending within a plane from a cylindrical epoxy body (Fig. 1). The needles have a diameter of 1.25 mm, and are 60 mm in length. All three needles contain a chromel-constantan thermocouple centered lengthwise in the needle. The center needle contains an additional 70-ohm heating element that extends the length of the needle. All needles are filled with an epoxy of high thermal conductivity to promote good thermal contact between the thermocouple, heating element, and outlying needle wall. Two of the sensors had an interneedle spacing of 6 mm, and the third sensor had an interneedle spacing of 12 mm to evaluate possible effects of needle spacing on measured values of soil heat flux. The sensors are installed so that the plane containing the three needles is normal to the soil surface. The center needle is used to measure
soil with the single-probe heat pulse technique described above, and the outer two needles are used to determine dT/dz.
|
plate =
soil), exactly 50% of this heat flux would travel through the heat flux plate, and the other half would conduct into the soil above the plate. However, because this is not generally the case, a correction is calculated as
![]() | [3] |
Lab Study
A comparison between two HFP01SC plates run with the self-calibration routine and three tri-needle probes was conducted under controlled conditions in the laboratory. We used a heat flux plate calibration chamber designed and constructed at the University of Wisconsin, and described in detail by Watts et al. (1990). This apparatus consists of a 45.7- by 50.8-cm rectangular heater plate at the bottom of a 10.2-cm deep chamber that is well insulated on three sides. The upper surface of the chamber is an aluminum heat sink, facilitating heat dissipation through the top of the chamber, thus creating a one-dimensional heat flux through the medium placed in the chamber. The total energy input to the heater plate can be easily determined, and heat losses through the insulation have been previously determined in detail, so the absolute magnitude of the heat flux through the chamber is well known (C.B. Tanner, unpublished data). The walls of the calibration chamber contain four sets of thermocouples in good thermal contact with the media in the chamber. These thermocouple pairs are oriented vertically, 2.86 cm apart, and are used to measure the dT/dz in the chamber. Since dT/dz and S are known, the thermal conductivity of the media in the calibration chamber can be calculated.
We filled the calibration apparatus with commercial quartz sandbox sand. This sand had a medium texture with particle sizes ranging from 0.4 to 0.5 mm. To avoid possible effects from flux divergence through the sides of the calibration chamber, the five sensors were arranged approximately equidistant (1112 cm) from the edges of the chamber, and were alternated in a circular pattern. This arrangement ensured that each sensor was at least 10 cm from adjacent sensors and chamber walls. The heat flux plates and the center needle of the tri-needle sensors were centered vertically in the chamber.
Soil thermal conductivity was determined hourly. A constant current module (CE8, Campbell Scientific, Logan, UT) was used to supply a 50-mA current to the heating elements in the center needle of each tri-needle probe, resulting in power dissipation (q from Eq. 2) of 2.917 W m-1 from each probe. Current was supplied to the heating elements for a total of 25 s. The temperature response was measured by the thermocouple in each center needle, and recorded by a datalogger (CR21X, Campbell Scientific, Logan, UT). Preliminary analysis indicated that 5 s of excluded data was sufficient to limit the effects of probe thermal properties, and achieve excellent linearity between ln t and T in sand (R2 > 0.999). However, we used a 10-s data rejection period to ensure that the effects of possible contact resistance in the field would have minimal effect on the calculation of
soil. The datalogger was programmed to calculate the slope of the least-squares linear regression of ln t vs. T (s in Eq. [2]). Temperature data from thermocouples in the outer two needles of each probe were recorded every 30 s and averaged across 0.5-h periods. Temperature gradient data were excluded for the first 5 min after the execution of the thermal conductivity measurement routine to minimize the effects of soil heat field perturbations caused by the probe heating. Soil heat flux was calculated every 30 min with the 30-min averaged value of dT/dz, and the most recent hourly measurement of
soil.
The HFP01SC correction routine was performed every two hours, and heat flux data were recorded every 30 s by a separate datalogger. These values were averaged across 0.5-h periods. The chamber thermocouples were read every 5 s, and averaged half-hourly.
Field Study
A comparative study was conducted between the tri-needle probes and two HFP01SC plates during the 2002 growing season at the University of Minnesota's Rosemount Research and Outreach Center, located 24 km south of St. Paul, MN (44°44' N, 93°05' W). The soil at this site is predominantly a Waukegan silt loam (fine-silty over sandy or sandy-skeletal, mixed, superactive, mesic Typic Hapludolls) with an average bulk density of 1.25 g cm-3. This field is farmed in a conventionally tilled corn (Zea mays L.) and soybean [Glycine max (L.) Merr.] rotation, and was planted in soybeans during this study. Three tri-needle probes (two with 6-mm needle spacing and one with 12-mm spacing) and two heat flux plates were installed on Day 128 (8 May 2002), and measurements were taken until Day 253 (10 Sept. 2002). All heat flux sensors were buried at a depth of 10 cm from the soil surface. The heat flux plates were installed by careful insertion into horizontal grooves that had been cut into a vertical soil face. The tri-needle probes were inserted into the vertical soil face with a guide to ensure that the needles remained parallel, and therefore that the spacing between the thermocouples in the outer needles remained uniform, and oriented so that the plane containing the needles was normal to the soil surface. Upon careful removal at the end of the field study, no changes in sensor spacing were apparent. However, such changes are a point of concern with this type of sensor.
The tri-needle sensors were operated in a manner identical to that described in the laboratory section. The HFP01SC correction routine was performed every 3 h, and heat flux data were recorded every 30 s. These values were averaged across 0.5-h periods. The heat storage in the soil layer above each sensor was calculated, but we have chosen to focus on sensor-based measurements, so these data have not been included in the following discussion. All sensors were within 1 m of each other, and after removal for planting, all sensors were placed between rows.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
soil was evaluated under controlled conditions in the laboratory. When undisturbed,
soil measurements by individual sensors varied <3% in all cases, and generally varied <2% during the course of several days in the medium-textured sand. In the laboratory, a uniform, homogenous media (sand) with a uniform thermal conductivity of 0.273 W m-1 K-1 was used. This allowed an evaluation of the precision of the thermal conductivity, thermal gradient, and soil heat flux measurements among the three tri-needle sensors. Table 1 shows intersensor variability between the three tri-needle sensors in the medium-textured sand during a soil warming period (Fig. 2) where S was altered from 0 to 126 W m-2. The performance of the tri-needle sensors was somewhat disappointing, with an average intersensor difference in soil heat flux of 5%. However, on further inspection it was determined that the majority of the error occurred as a result of differing response times of individual sensors during the period of rapid heating.
|
|
soil by comparing these quantities measured by the tri-needle sensors with independent measurements by the calibration chamber. Interestingly, it was found that the tri-needle probes overreported
soil by 3%, and underreported dT/dz by 5%. The excellent overall agreement in reference and measured S, and the fact that the errors in
soil and dT/dz are opposite in sign suggest that edge effects and nonhomogeneity of the sand near the built-in chamber thermocouples could result in flawed reference values for dT/dz (and consequently, reference
soil). However, even if the errors in
soil and dT/dz had been additive, the magnitude of the absolute error of 8% is still far less than what was observed in the heat flux plates. The large overestimation of S by the Hukseflux plates may be the result of the very low thermal conductivity of the dry sand (0.273 W m-1 K-1). The greater thermal conductivity of the heat flux plate material (0.8 W m-1 K-1) will cause flux convergence through the plate, and thus artificially high heat fluxes. The self-calibration routine of the Hukseflux plates should account for this, but may not be effective in media with very low thermal conductivities.
|
Other researchers have noted errors in the measurement of thermal conductivity arising during times of rapid soil heating (K. Bristow, personal communication, 2002). Under these conditions, the slope of the T vs. ln t plot (s in Eq. [2]) can be affected by changes in ambient soil temperature. The temperature change in the chamber immediately after the onset of heating averaged 4.25°C hr-1, which is similar to the maximum rate of soil temperature change measured in an agricultural field in central Minnesota at a depth of 5 cm. The thermal conductivity measurements in the dry sand changed by an average of 0.22%, indicating that errors in the determination of thermal conductivity due to rapid changes in ambient soil temperature are not significant under normal field conditions at depths where soil heat flux sensors are normally deployed. These results, and the good accuracy that has been achieved in the determination of soil thermal conductivity indicate that these sensors are suitable substitutes for soil heat flux plates at depths where they are normally installed in the combination method (510 cm). This is a significant improvement in the application of the gradient method over that described by Kimball et al. (1976) where lack of accuracy in the calculation of soil thermal conductivity prevented application of the gradient method at a depth of <20 cm.
Field Study
The three tri-needle probes were deployed in a field setting during the growing season of 2002 alongside two HFP01SC plates, with the intention of evaluating the overall performance of the needle sensors, isolating errors associated with fluid flow blockage by the HFP01SC plates. In the field, spatial variability in the thermal properties of the soil and the surface energy balance, especially under a partially closed plant canopy, confounds attempts to intercompare soil heat flux sensors. However, some general conclusions about sources of error can be drawn from an extended field data set.
The overall agreement among all five sensors was generally good (Fig. 3). However, the lack of a reliable reference, the effects of spatial variability on the thermal properties of the soil, and differences in heat storage above the sensors make extended analysis of agreement among the sensors impossible. Soil thermal conductivity during the season ranged from 0.672 to 2.24 W m-1 K-1 as water content changed, with a mean of 1.26 W m-1 K-1. Agreement between the heat flux plates and tri-needle sensors was not correlated with soil thermal conductivity, indicating that the Hukseflux self-calibration routine was effective in correcting for flux divergence in the field setting, in contrast with the measurements in the calibration chamber with dry sand.
|
| CONCLUSIONS |
|---|
|
|
|---|
A field intercomparison among the tri-needle sensors and commercially available heat flux plates yielded good agreement, with relative differences that were generally small. Expected errors associated with heat flux plates acting as barriers to fluid flow in soils were not observed in more than 100 d of field use.
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. E. Ochsner and J. M. Baker In Situ Monitoring of Soil Thermal Properties and Heat Flux during Freezing and Thawing Soil Sci. Soc. Am. J., July 1, 2008; 72(4): 1025 - 1032. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. Sauer, T. E. Ochsner, and R. Horton Soil Heat Flux Plates: Heat Flow Distortion and Thermal Contact Resistance Agron. J., January 1, 2007; 99(1): 304 - 310. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. E. Ochsner, T. J. Sauer, and R. Horton Field Tests of the Soil Heat Flux Plate Method and Some Alternatives Agron. J., June 5, 2006; 98(4): 1005 - 1014. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Mori, J. W. Hopmans, A. P. Mortensen, and G. J. Kluitenberg Estimation of Vadose Zone Water Flux from Multi-Functional Heat Pulse Probe Measurements Soil Sci. Soc. Am. J., April 11, 2005; 69(3): 599 - 606. [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 | |||