Published in Soil Sci. Soc. Am. J. 68:66-73 (2004).
© 2004 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
DIVISION S-1SOIL PHYSICS
A Simplified Falling-Head Technique for Rapid Determination of Field-Saturated Hydraulic Conductivity
V. Bagarello*,a,
M. Iovinoa and
D. Elrickb
a Dipartimento di Ingegneria e Tecnologie Agro-Forestali, Università degli Studi di Palermo, Viale delle Scienze, 90128, Palermo, Italy
b Dep. of Land Resource Science, Univ. of Guelph, Guelph, ON, Canada N1G 2W1
* Corresponding author (bagav{at}unipa.it).
 |
ABSTRACT
|
|---|
Simplified measurements of the field-saturated hydraulic conductivity, Kfs, require short duration experiments, small water volumes, and easily transportable equipment. A simplified falling-head (SFH) technique for the rapid determination of Kfs has been developed and tested. The technique consists in applying a small volume of water on a soil surface, confined by a ring inserted a short distance into the soil, and then measuring the time from the application of water to the instant at which the surface area is no longer covered by water. A measurement of the initial and field-saturated soil water contents, and an estimate of the
* parameter of the Gardner's exponential model are then used to calculate Kfs using a simple solution that includes gravity. The Kfs of both repacked and undisturbed soil cores was determined in the laboratory by the SFH and the early time constant-head (ECH) techniques. The SFH and (constant-head) pressure infiltrometer (PI) techniques were then compared in the field. The maximum discrepancy between the mean Kfs results obtained within an experiment was of a factor of approximately two. This difference is negligible in most practical applications and it was concluded that the SFH technique compared favorably with the ECH technique in the laboratory and to the PI technique in the field. The SFH technique appears promising for determining Kfs in a relatively short period of time without the need for extensive instrumentation or analytical methodology, and therefore it appears suitable for detailed field measurements over large areas.
Abbreviations: CV, coefficient of variation ECFH, sequential early time constant head/early time falling head ECH, early time constant head EFH, early time falling head GA, Green and Ampt Kfs, field-saturated hydraulic conductivity OPD, one-ponding depth PI, pressure infiltrometer SFH, simplified falling head TPD, two-ponding depth
 |
INTRODUCTION
|
|---|
THE HYDRAULIC CONDUCTIVITY of saturated soil is one of the most important soil properties controlling water infiltration and surface runoff, leaching of pesticides from agricultural lands, and migration of pollutants from contaminated sites to the ground water.
Saturated hydraulic conductivity depends strongly on soil texture and structure, and therefore can vary widely in space. Since hydraulic conductivity is determined essentially at points on a field scale, a large number of determinations is required to assess the magnitude and structure of the variation within the selected area (Logsdon and Jaynes, 1996). Spatially distributed determinations of hydraulic conductivity have to be repeated at different times, particularly in soils where structure varies over time because of natural or anthropogenic factors (Prieksat et al., 1994). For structured soils in particular, saturated hydraulic conductivity has to be measured directly in the field to minimize disturbance of the sampled soil volume and to maintain its functional connection with the surrounding soil (Bouma, 1982). Using small volumes of water, easily transportable equipment, and conducting short-duration experiments is desirable to obtain Kfs data at a great number of locations over a large area and with the realistic use of resources in terms of time and costs.
Most field techniques, such as the traditional well permeameter and ring infiltrometer constant-head techniques (Reynolds and Elrick, 1986, 1990), rely on the attainment of a steady-state flow rate of water into the soil (Elrick and Reynolds, 1992b; Reynolds, 1993), and this can limit their field use since short-duration steady-state experiments can be conducted only in relatively permeable soils (Elrick and Reynolds, 1992a). Other studies have shown, however, that short duration experiments can also be conducted in low permeability soils using early time transient flow measurements from a surface ring to obtain estimates of Kfs (Fallow et al., 1994; Elrick et al., 1995; Odell et al., 1998). The simplest technique that uses the initial period of the infiltration process to deduce Kfs is probably the ECH technique. This technique was specifically developed for short-time measurements of the saturated hydraulic conductivity of very low permeability soils (i.e., Kfs < 1 x 109 ms1) (Fallow et al., 1994; Odell et al., 1998). For both steady-state and ECH infiltration experiments, a constant head device has to be used, and transported throughout the area of measurement. Using a falling-head experiment to determine Kfs in situ is preferable since this type of experiment is simpler than a constant-head experiment (Philip, 1992). This approach has been used recently for determining Kfs from borehole measurements (Munoz-Carpena et al., 2002).
Using previous work conducted by Smith (1999) and Elrick et al. (2002), we developed a SFH technique for measuring Kfs. Field application of this technique requires only a small ring that is inserted to a short distance into the soil, a small volume of water, and a stopwatch. In this paper, the simplified method for measuring Kfs is presented and a laboratory and a field evaluation of its performance is performed.
 |
A Review of Some Transient Flow Equations and Their Applications
|
|---|
The initial flow of water from a surface ring into initially unsaturated soil is controlled almost solely by the pressure head and capillary sorption and it can be considered to be one-dimensional and free from gravitational effects so that the cumulative infiltration of water, I (L), can be described by (Philip, 1957, 1969):
 | [1] |
where SH (L T1/2) is the soil sorptivity for a ponded head, H (L), and t (T) is the time.
During constant-head conditions, SH is invariant and it can be expressed by the following approximate analytical relationship, that was deduced using previous work by White and Sully (1987) (Fallow et al., 1994):
 | [2] |
where 
is the difference between the field-saturated water content,
fs (L3 L3), and the initial water content,
i (L3 L3),
m (L2 T1) is the field-saturated matric flux potential, b is a constant equal to 0.55 (White and Sully, 1987), and Kfs (L T1) is the field-saturated hydraulic conductivity. The first term in Eq. [2] gives the sorptivity, S0, for H = 0 and the second term gives the increase in sorptivity due to the positive (ponded) head, H. The Green and Ampt (1911) (GA) equation can also be used to derive another approximate expression of SH for constant-head conditions:
 | [3] |
where
f (L) is the soil water pressure head at the wetting front, and
f is negative. Since, for rectangular hydraulic conductivity function like GA,
m and
f are related by (Elrick et al., 2002):
 | [4] |
the only difference between Eq. [2] and [3] is in the estimation of b that is equal to 0.5 for GA soils. Solving Eq. [2] for Kfs gives:
 | [5] |
and solving Eq. [3] for Kfs gives:
 | [6] |
where the
* parameter (L1) is given by (Elrick and Reynolds, 1992a):
 | [7] |
Equations [5] and [6] can be used for determining Kfs from early time flow data under a constant head (Fallow et al., 1994; Odell et al., 1998; Elrick et al., 1995, 2002). The procedure, the ECH technique by Fallow et al. (1994), evaluates SH from the slope of the linear portion of the I vs. t1/2 graph, and needs the measurement of 
and the estimation of the
* parameter (Elrick and Reynolds, 1992b).
During early time falling head (EFH) infiltration, SH is variable and therefore does not exhibit a linear behavior on a graph of I vs. t1/2 (Odell et al., 1998). For falling-head conditions, the cumulative infiltration is a function of H as follows (Elrick et al., 1995):
 | [8] |
where a (L2) is the cross-sectional area of the falling-head tube, A (L2) is the cross-sectional area of the infiltrating surface and H0 (L) is the height of the ponded head at t = 0. Fallow et al. (1994) used Eq. [2] to approximate a falling-head infiltration process based on a succession of constant-head conditions:
 | [9] |
where the impact of the falling head is given by the term H(t). Substitution of Eq. [8] into Eq. [9] yields (Fallow et al., 1994; Elrick et al., 1995):
 | [10] |
Equation [10] was fitted to EFH data by a nonlinear, least-squares procedure to determine simultaneously Kfs and
m (EFH technique) (Fallow et al., 1994). Since H0 was difficult to obtain, the EFH technique was later extended by Elrick et al. (1995) by the addition of an initial period of constant-head infiltration at H = H0 (Sequential Early-Time Constant-Head/Falling-Head [ECFH] technique). Recently, Elrick et al. (2002) used the GA approach to deduce a new implicit equation in I for falling-head conditions that neglects gravity:
 | [11] |
where:
 | [12] |
The GA expression for a falling head that includes gravity was derived by Philip (1992):
 | [13] |
where:
 | [14] |
Equations [11] and [13] are valid until the falling head drops to zero.
To evaluate the applicability of early time techniques, the accuracy of the Kfs determination has to be established by also distinguishing between early time techniques requiring an independent estimation of
* (ECH technique) and early time techniques allowing simultaneous determination of Kfs and
m or
* (EFH and ECFH techniques). See Elrick and Reynolds (1992b) for a discussion of the two procedures referred to above as applied to the constant-head permeameter. Only a few investigations have been conducted to evaluate the reliability of the estimates of Kfs and
m (or
*) obtained from early time ponded flow experiments and most investigations have used the solutions deduced from Eq. [2] (Vauclin et al., 1994; Elrick et al., 1995; Odell et al., 1998).
It has been well established that field-saturated hydraulic conductivity tends to be obtained more reliably than
* from steady, ponded flow out of a ring since the hydrostatic pressure and gravity components of flow are maximized rather than the capillary component (Reynolds and Elrick, 1990; Vauclin et al., 1994; Elrick et al., 1995; Bagarello et al., 2000; Mertens et al., 2002). Also,
* is obtained from Kfs/
m and thus contains the variability in the calculation of both Kfs and
m. The estimates of Kfs appear more stable than those of
* (Elrick et al., 2002) and not greatly affected by the early time technique used (Odell et al., 1998), suggesting that early time ponded flow experiments can be used to estimate Kfs. In this sense, a simplified approach with a priori estimate of the
* parameter appears more attractive than a more complicated approach involving additional loads of an experimental and computational nature.
 |
The Simplified Falling-Head Technique
|
|---|
The SFH technique consists in applying quickly a small volume of water, V (L3), on the soil surface confined by a ring (a/A = 1) inserted to some distance into the soil and in measuring the time, ta (T), from the application of water to the instant at which the surface area is no longer covered by water. At t = ta, I(ta) = H0 = D, where D = V/A (L) is the depth of water corresponding to V. Equation [13], with Eq. [7] inserted, then becomes:
 | [15] |
Estimation of Kfs by Eq. [15] requires measurement of both ta and 
and the estimation of the
* parameter. No theoretical limitations exist on the area to be sampled, but there are practical limitations. Knowledge of 
and of the depth of insertion of the ring into the soil allows the investigator to determine the volume of voids within the soil volume confined by the ring. A volume of water less than or equal to the volume of voids can then be selected to assure one-dimensional flow during the experiment. Equation [15] is precisely Eq. [5] of Philip (1992) solved for Kfs. Since Eq. [15] includes gravity, the only time limitation will occur if the wetting front emerges from the bottom of the ring and three-dimensional flow commences.
For three Kfs
* combinations (Kfs = 1 x 104 m s1 and
* = 36 m1; Kfs = 1 x 106 m s1 and
* = 12 m1; Kfs = 1 x 108 m s1 and
* = 4 m1) and two values of 
(
= 0.2 m3m3 and 
= 0.4 m3m3), Figure 1
shows the relationship between ta and D according to Eq. [15]. For illustrative purposes, D was allowed to vary between 0.005 and 0.03 m (for 
= 0.2 m3 m3) or 0.06 m (for 
= 0.4 m3 m3). These maximum values of D correspond to the volume of water required to fill the initially air-filled voids in the soil volume confined by the ring. An insertion depth of 0.15 m was considered since it appears to be a reasonable compromise between the need of avoiding excessive disturbance of the sampled soil volume and that of assuring one-dimensional flow. For Kfs = 1 x 104 m s1, ta is generally <5 min while ta is >3.5 h for Kfs = 1 x 108 m s1. It appears, therefore, that falling-head experiments of reasonable duration can potentially be conducted on a broad range of soils. In particular, the SFH experiment appears more attractive than the traditional constant-head field techniques for determining the field-saturated conductivity of relatively low permeability soils. The traditional methods typically rely on the attainment of a steady-state flow rate that can take several days to occur (Odell et al., 1998). The same falling-head experiment can be also used to determine soil sorptivity if gravity effects are negligible (Smith, 1999).
The effect of minor changes in D, 
, ta, and
* on the predictions of Kfs obtained by Eq. [15] was tested for each of the 54 (Kfs,
*, D, 
, ta) combinations considered in Fig. 1. The relative error in Kfs, eK, for a given data set of
*, D, 
, and ta values (denoted by the variable ß) was calculated using the following relationship:
 | [16] |
where Kfs(ßp) is the Kfs value obtained by perturbing one of the input data while maintaining the others constant. A variation of ±10% of the true value was considered for D, 
, and ta. A variation of plus or minus one category among those proposed by Elrick and Reynolds (1992a) was considered for the
* parameter; for example, 12 m1 was replaced by 4 m1. The values of eK ranged from 9.0 to +11.0% for both 
and ta. Values of eK ranging from 19.0 to +21.0% and from 74.0 to +179.0% were obtained for D and
*, respectively.
Sensitivity of Kfs to minor changes in D, 
, and ta was small and probably negligible in practice, given that Kfs varies over several orders of magnitude in nature. An implication of this result was that in most cases the determination of Kfs should not be greatly affected by the assumption that the infiltration surface was smooth or by small errors (i.e., ±10% of the true value) in the measurement of the soil water content or the duration of the infiltration process. Moreover, the assumption that the initial depth of ponding was established instantaneously on the infiltration surface appeared to be generally reasonable. An improper evaluation of
* by one category yielded a Kfs value differing appreciably from the true value. Similar levels of error have been reported for steady-state experiments conducted with the Guelph permeameter and PI techniques (Reynolds and Elrick, 1990; Elrick and Reynolds, 1992a; Bosch, 1997) and they have been considered acceptable for many practical applications (Elrick and Reynolds, 1992a). Moreover, an erroneous choice of
* appears rather unlikely, given that only four categories were suggested by Elrick and Reynolds (1992a) for describing a very wide range of soils.
 |
MATERIALS AND METHODS
|
|---|
Laboratory Experiments
Experiments were conducted on both undisturbed soil cores and packed soil columns to compare the hydraulic conductivity values obtained by the ECH and SFH techniques. These techniques show similarities since both use the initial portion of the infiltration process, need an estimate of the
* parameter and sample small and relatively similar soil volumes.
Nine undisturbed soil cores having a silty clay loam texture (Soil Survey Staff, 1998) were collected in 0.085-m diam. x 0.11-m high stainless steel cylinders from the upper layer of a soil classified as Chromic Calcixerert (Soil Survey Staff, 1998), after eliminating the first 0.02 to 0.03 m of soil. Each cylinder was inserted vertically into the soil by hammering gently on the top of the cylinder with a rubber hammer and progressively removing the surrounding soil up to the established depth to reduce disturbance during sampling. The height of the soil cores varied between 0.055 and 0.092 m with a mean value of 0.072 m. In the laboratory, a nylon filter and a rigid wire screen were connected to the base of the core to prevent soil loss from the bottom of the sample and the cores were weighed to determine the initial soil water content. A small depth of water, D = 0.015 m, was then applied on the soil surface and the time, ta, from the application of water to the instant at which the surface area was no longer covered by water was measured. Depending on the soil core, ta varied between 141 and 9914 s, with a mean value of 1575 s. For all but two experiments (cores C4 and C7, having a volume of air-filled voids of 75.5 and 79.8 cm3, respectively), the infiltrated volume of water (V = 85.1 cm3) was less than the volume of air-filled voids at the time of the experiment. Soil cores were then exposed to air in the laboratory to promote loss of water. After about 10 d, the cores were weighed again and a constant head of water, H, of 0.014 m was established on the soil surface. Early time infiltration data were collected during the first 679 to 14413 s (mean = 2645 s) and they were used to obtain an estimate of SH. Starting from the collection time, separation of the soil from the cylinder wall was never observed for these undisturbed cores. The soil cores were then oven-dried and weighed. For each soil core, two estimates of the Kfs were obtained. In particular, a value of Kfs was obtained from the SFH experiment, using Eq. [15], and another value of Kfs was obtained from the ECH experiment, using Eq. [6]. To maintain simplicity of the tested early time methods, the field-saturated soil water content was assumed to be equal to total porosity that was calculated from dry bulk density measurements. A value of
* = 12 m1 was used to calculate Kfs (Elrick and Reynolds, 1992a, b).
Laboratory experiments were also conducted to determine the Kfs of a compacted artificial soil composed by equal weights of a mono-granular sand (particle diameter, dp = 73 µm) and a commercial kaolin powder (median diameter d50 = 3 µm). This mixture was chosen since preliminary tests showed its ability to provide values of Kfs that were reasonably low for the purposes of the investigation and from cores that did not swell appreciably on wetting. Forty-two 0.05-m diam. by 0.04- or 0.05-m high columns having a mean bulk density of 1.38 g cm3 (coefficient of variation [CV] = 1.6%) were prepared by using the air-dried mixture. A step-by-step procedure was applied to prepare the core. At each step, a fixed amount of artificial soil, corresponding to a height of the core of about 0.01 m, was put in the cylinder and then it was manually compacted. Efforts were made to give similar amounts of compacting energy to each layer. The SFH technique was used to determine the Kfs of 22 cores. A depth of water, D = 0.01 m, was used for these experiments. Depending on the core, ta varied between 442 and 2841 s, with a mean value of 1301 s. The ECH technique was applied to determine the Kfs of 20 cores. A constant depth of ponding, equal to 0.01 m, was maintained to determine SH from early time infiltration measurements that were collected for the initial 133 to 347 s (mean = 270.2 s), depending on the experiment. Equation [6] was used to calculate Kfs. For both techniques (SFH, ECH),
i was set equal to the mean of the measured values for the air-dried material (
i = 0.03 m3m3) and
fs was estimated from bulk density measurements. The mean depth of the wetted layer was calculated to be equal to 0.022 m for the SFH experiments and to 0.010 m for the ECH experiments. Published values of
* were not available for the artificial soil used in the repacked soil experiments. Therefore, to establish the value of
* to be used in the calculations, the wetting retention curve of three replicated soil columns was determined in the pressure head range 1200
h
10 mm. The Brooks and Corey (1964) equation was then fitted to the water retention data by using the RETC code (van Genuchten et al., 1991) and an estimate of
* was obtained from (Fallow et al., 1994):
 | [17] |
where h is the pressure head, hi is the initial (negative) pressure head, Se = (
r)/(
s
r) is the effective saturation,
s and
r are the saturated and the residual water content, respectively,
is a pore-connectivity parameter, set equal to 0.5 (Mualem, 1976), and
is the exponent of the Brooks and Corey (1964) equation. A value of
* equal to 2.0 m1 was obtained and used in the calculations.
Each set of Kfs values was summarized by calculating the mean
and the CV of the individual measurements. The statistical frequency distributions for the Kfs data obtained from the undisturbed soil core experiments were found to be log-normal (data not shown), which is common for undisturbed soil (e.g., Lee et al., 1985; Warrick, 1998). Geometric means and the associated coefficients of variation were therefore calculated (Lee et al., 1985). The statistical frequency distributions for the Kfs data obtained from the repacked soil cores were found to be normally distributed (data not shown). This result was expected since unstructured soil was used for these experiments. Arithmetic means and the associated coefficients of variation were therefore determined.
Field Experiment
Field experiments were conducted to compare the SFH technique with the PI (Reynolds and Elrick, 1990) technique for measuring Kfs. Reasons for choosing this last technique for the comparison include: (i) the PI technique is one of the simplest techniques for measuring Kfs in the field, (ii) it appears to yield Kfs results that are comparable with other techniques (e.g., tension infiltrometer, laboratory soil cores) for different soil type and land management combinations (Reynolds et al., 2000), and (iii) the same surface area can be sampled with the SFH and PI techniques and this should simplify the interpretation of the Kfs results.
The field site is located near Menfi (western Sicily, Italy), approximately 90 km south of Palermo, at an altitude of 20 m above sea level. The field is flat and supports a 2-yr-old espalier system vineyard (Vitis vinifera) with a distance between the rows of vines of 2.20 m. The soil is classified as Chromic Calcixerert (Soil Survey Staff, 1998). An area of approximately 40 by 4.40 m2 was chosen for the investigation and 30 sites were randomly selected within this area. Soil texture of the 0- to 0.05- and 0.05- to 0.10-m layers was determined on 30 soil samples collected in the vicinity of each site (Table 1). The soil was classified as loam (Soil Survey Staff, 1998) in most cases (87% of the total) and sandy clay loam or sandy loam in a few cases (8 and 5% of the total, respectively). The gravel content was negligible.
View this table:
[in this window]
[in a new window]
|
Table 1. Basic statistics of the clay, silt, and sand contents measured at two different depths for each of the 30 sites selected for the field experiment.
|
|
Starting on 25 July 2002, a rectangular area of approximately 0.60 by 0.40 m2 was bounded at each site by small ridges and 10 to 12 L of water were poured to facilitate ring insertion, given that the soil appeared quite hard during the preliminary survey. Two to 6 d after wetting, a ring having an inner diameter of 0.135 to 0.149 m (depending on the ring) and an height of 0.25 m was inserted vertically by hand to a depth of 0.12 m. Ring insertion was conducted very gently and inspection of the soil surface after ring insertion showed that there was no appreciable soil compression within the area bounded by the ring. This suggested that soil compression was not a major factor affecting the applicability of the SFH technique.
At each site, an experiment with the SFH technique was conducted both before and after the PI experiment to also test the effect of the sequence of the experiments on the Kfs results. Different authors have shown that repeating an infiltration experiment on the same soil core can have an appreciable influence on the Kfs results (Bagarello et al., 2000; McKenzie et al., 2001). The first SFH experiment was conducted immediately after ring insertion and the PI experiment was performed 10 to 37 d after the SFH experiment. Eight to 30 d after the PI application, 15 L of water were poured on each bounded area, including the surface area confined by the ring, to reduce differences between the water content of the soil volume confined by the ring and that of the soil surrounding the ring. The second SFH experiment was conducted on 1 Oct. 2002, 20 d after this last wetting. A total rain amount of 18.2 mm was registered during the experimental period. Most of this rain (i.e., 14.4 mm) occurred after the PI application. The daily temperature varied between 17.5 and 31.0°C, with a mean value of 23.7°C.
For the SFH experiments, a pre-established volume of water was applied on the soil surface in a few seconds (i.e., approximately 5 s) and the time ta from the application of water to the instant at which the surface area was no longer covered by water was measured. For comparative purposes, the time, t0.5a, at which half the surface area was no longer covered by water was also recorded. Two undisturbed soil cores (0.05-m diam. by 0.05-m high) were collected within the bounded area at a depth of 0 to 0.05 and 0.05 to 0.10 m, respectively, before each SFH experiment. These cores were used to determine the bulk density and the initial volumetric soil water content that were averaged over the two depths. Field-saturated water content was calculated using bulk density results. Before conducting each SFH experiment (as well as each PI experiment), the soil surface confined by the ring was gently moved with a spatula to remove any surface crust and to improve the contact between the soil and the inside wall of the ring. A site-specific determination of 
was not obtained before each experiment. A constant value of V (V = 429.4 cm3, 0.025
D
0.030 m, depending on the ring), giving one-dimensional flow for 
0.25 m3 m3, was used for all the SFH experiments. One-dimensional flow was however checked after the SFH experiment. The depth of the wetting front, dw (L), at the end of the SFH experiment was calculated by the following relationship:
 | [18] |
and it was <0.12 m in all but two cases (dw = 0.121 and 0.138 m, respectively). The mean value of dw was equal to 0.088 m (standard deviation,
= 0.018 m) for the first set of SFH experiments and 0.080 m (
= 0.016 m) for the second set of SFH experiments. Equation [15] was used to calculate Kfs for each SFH experiment. According to Elrick and Reynolds (1992a), a value of the
* parameter equal to 12 m1 was used for the calculations.
The PI experiments were conducted by using two infiltrometers consisting of a Mariotte reservoir 1.0-m high with an inner diameter of 0.11 m. Two pressure heads were established on the soil in succession and in ascending order (H1 = 0.05 m, H2 = 0.107 m). After some preliminary tests, the duration of the experiment was set to 4 h for the first pressure head and to 3 h for the second pressure head. This duration allowed us to obtain apparent steady-state conditions for both steps of each experiment. For one site only, a unique ponding head (H = 0.05 m) was used because the infiltration process was particularly slow and the duration of the experiment exceeded 11 h. The rate of fall of the water level in the infiltrometer reservoir was monitored at 2- to 5-min time intervals. In many cases, high infiltration rates were observed and the water reservoir emptied before concluding the experiment. Refilling of the reservoir was conducted by maintaining ponded conditions on the soil surface confined by the ring. Apparent steady state was detected on a flow rate versus time plot. Flow at the end of the first step of the experiment was three-dimensional because the volume of water discharged from the PI reservoir was greater than the volume of the initially air-filled voids in the soil volume confined by the ring (a dry initial condition, i.e., 
= 0.4 m3m3, was assumed for these calculations) plus the volume of water required to maintain H1 = 0.05 m above the soil surface. The three-dimensional flow is taken into account by the PI equation. Both the two-ponding depth (TPD) and the one-ponding depth (OPD) approaches were applied to calculate Kfs (Reynolds and Elrick, 1990). A constant value of the
*-parameter, equal to 12 m1, was used for the OPD approach.
The geometric mean and the associated coefficient of variation were calculated for the Kfs results obtained in the field with both the SFH and PI experiments because they were found to be log-normally distributed (data not shown).
 |
RESULTS AND DISCUSSION
|
|---|
Laboratory Experiments
On average, the ECH Kfs estimate of the undisturbed soil cores was 1.82 times less than the value determined by the SFH technique and the two techniques produced values of CV differing by a factor of 1.04 (Table 2). Correlation between the log-transformed values of Kfs obtained with the two techniques was statistically significant (P = 0.05, coefficient of determination, r2 = 0.65). However, the observed discrepancies are negligible from a practical point of view, given that many authors have suggested that an error of the estimate of Kfs by a factor of two or three can be considered acceptable for many practical purposes (Elrick and Reynolds, 1992a; Reynolds and Zebchuk, 1996; Elrick et al., 2002). A tendency of the SFH technique to produce higher values of Kfs as compared with the ECH technique was however recognized. We felt two possible reasons for the observed discrepancies. First, a short-term surface soil structure degradation has been recognized to determine a reduction of Kfs measurements conducted repeatedly on a given soil core (Bagarello et al., 2000). The ECH experiments were conducted after the SFH experiments that, likely, were effective in promoting some form of degradation of the upper layer of the core. Second, the initial soil water content was higher for the ECH experiments (mean 
= 0.191 m3 m3, CV = 40.1%) than for the SFH experiments (mean 
= 0.248 m3 m3, CV = 23.8%). Probably closure of macropores occurred as the soil water content increased and this determined a reduction of hydraulic conductivity.
View this table:
[in this window]
[in a new window]
|
Table 2. Basic statistics of the field-saturated hydraulic conductivity values, Kfs (m s1), measured on nine undisturbed soil cores by the early time constant-head (ECH) and simplified falling-head (SFH) techniques.
|
|
For the packed soil columns, the arithmetic mean value of Kfs obtained by the ECH technique was 1.23 times higher than that determined by the SFH technique and the associated coefficients of variation differed by a factor of 1.05 (Table 3). Therefore, the two techniques produced practically identical estimates of both the mean and the coefficient of variation of Kfs.
View this table:
[in this window]
[in a new window]
|
Table 3. Basic statistics of the field-saturated hydraulic conductivity values, Kfs (m s1), measured on the packed soil columns by the early time constant-head (ECH) and simplified falling-head (SFH) techniques.
|
|
Field Experiment
Using t0.5a instead of ta for calculating Kfs by Eq. [15] yielded
values and associated CVs differing by 1.21 and 1.08 times, respectively, for the first SFH experiment and by 1.19 and 0.96 times, respectively, for the second SFH experiment. These differences were small, confirming that minor variations in the estimate of the duration of the infiltration process did not have a great effect on the Kfs results. The Kfs values determined by using the ta measurements were used for method comparison since ta is less subjectively evaluated than t0.5a.
The TPD approach produced positive values of Kfs and
* for 21 PI experiments (success rate = 72%,
= 6.29 x 106 m s1, CV = 198.6%). For these experiments, the geometric mean value of
* was equal to 9.5 m1 and this result supported the choice of using
* = 12 m1 for the OPD calculations. This level of support was considered to be enough for this investigation because the sensitivity of
* to measurement errors using the TPD approach is high and reasonably accurate estimates of Kfs can be obtained using
* values that are largely erroneous (Mertens et al., 2002). The Kfs values obtained by the OPD approach were used in the method comparison to include data from all sites.
The two sets of SFH experiments produced
values that differed by a factor of 1.93 and that were 1.14 (second set) to 2.20 (first set) times higher than the
value obtained with the PI technique (Table 4). The CVs of the Kfs data obtained with the transient approach were similar, differing by a factor of 1.11, and they were 1.73 to 1.92 times lower than the CV value calculated for the steady-state results. The first and the second application of the SFH technique yielded higher Kfs results than the PI technique for 24 (80% of the total) and 14 (47% of the total) experiments, respectively. The maximum values of Kfs obtained with the first set of SFH experiments (Kfs = 4.05 x 105 m s1) and the PI technique (Kfs = 3.58 x 105 m s1) were similar and they were approximately two times higher than the maximum Kfs measured with the second application of the SFH technique (Kfs = 1.79 x 105 m s1). In all cases, a statistically significant (P = 0.05) correlation was detected between the log-transformed Kfs data, with values of r2 ranging from 0.51 (first SFH experiment vs. PI experiment) to 0.77 (first SFH experiment vs. second SFH experiment).
View this table:
[in this window]
[in a new window]
|
Table 4. Basic statistics of the field-saturated hydraulic conductivity, Kfs (m s1) values obtained at the 30 field sites with the two simplified-falling head (SFH) and the pressure infiltrometer (PI) experiments.
|
|
In summary, the transient and steady-state techniques produced mean values of Kfs differing by a maximum factor of approximately two and the transient approach gave Kfs predictions that were less variable than the steady-state approach. Therefore, the two techniques produced estimates of
that were practically equivalent. This is particularly encouraging as the transient technique took much less time and equipment. As a matter of fact, 30 SFH experiments, including undisturbed soil core collection at two depths, were concluded within a working day. Working for approximately 15 h (also using artificial light) and using two instruments, a maximum of four PI experiments per day were performed.
An attempt to interpret the Kfs results was made notwithstanding that this can be an imprecise enterprise due to the lack of independent Kfs datum or benchmark on which evaluations and judgments can be made (Reynolds et al., 2000).
The sampled soil volumes were wetted only once before the first SFH experiment and on several occasions before the second SFH experiment and this probably contributed to the observed discrepancies between the two sets of SFH data. Maintaining moist conditions for an extended period likely promoted some form of soil structure weakening that reduced macroporosity (Beven and German, 1982; Hillel, 1998) and increased biological activity (McKenzie et al., 2001). These factors probably led to a general reduction in the Kfs of the upper soil layer sampled by the SFH technique.
An implication of this interpretation is that the results of the first set of SFH experiments should be considered as a more meaningful comparison between the SFH and PI techniques. Structure modification induced by wetting processes between the first SFH experiment and the PI experiment was negligible since in practice the soil was only wetted by the water supplied with the SFH experiment. Natural soil structure modifications were also negligible because the PI experiment was conducted a relatively short time after the first SFH experiment (with a few exceptions, <25 d).
The method comparison showed that the SFH technique tended to produce slightly higher and less variable Kfs results than the PI technique. Probably, the different duration of the two types of experiment contributed to this discrepancy. During a ponded infiltration experiment, the ability of the largest pores to transport flow likely decreases as the duration of the experiment and the sampled soil volume increases. This depends on the fact that both the tortuosity of the flow pathways (Hillel, 1998) and the potential for the occurrence of infiltration processes into the soil matrix through the large pore walls (Beven and German, 1982) tend to increase. Soil swelling can be also promoted by a long-duration infiltration experiment (Bagarello et al., 1999). Although the clay type was not determined, some soil swelling is expected to occur at the field site given that small cracks (i.e., <0.01 m large) were observed under dry conditions in the undisturbed areas of the field. All these factors reduce the rate of water entry into the soil through the infiltration surface. Consequently, a long-duration experiment can yield a value of Kfs equal to or lower than a short-duration experiment. Soil properties controlling flow transport (e.g., size, shape, and connectivity of the pores) vary greatly in space (Warrick, 1998) and therefore the differences in Kfs obtained by the two approaches are expected to be site-dependent. A higher frequency of low Kfs is however expected and this can produce a higher variability of the individual measurements.
 |
CONCLUSIONS
|
|---|
Measuring Kfs at a great number of locations over a large area with a realistic use of resources in terms of time and costs needs a short duration experiment, small volumes of water and easily transportable equipment.
A SFH technique for the rapid determination of Kfs has been proposed and verified here. The technique consists of applying quickly a small volume of water on the soil surface confined by a ring inserted a short depth into the soil and then measuring the time from the application of water to the instant at which the surface area is no longer covered by water. This time, together with a measurement of initial and field-saturated soil water content and an estimate of the
*-parameter, is used to calculate Kfs by a simple solution based on the GA assumptions. This solution includes gravity. The duration of the experiment ranges from <5 min for permeable soils (i.e., Kfs
1 x 104 m s1) to more than 3.5 h for soils with a relatively low permeability (i.e., Kfs
1 x 108 m s1). This theoretical analysis suggested that the SFH technique could be potentially applied to a broad range of soils, working faster and easier in permeable and moderately permeable soils. Minor errors (i.e., ±10% of the true value) in the measurement of the initial water depth, the duration of the infiltration process or the soil water content do not affect appreciably the determination of Kfs. The SFH technique can also be used to measure Kfs of detached soil cores.
Laboratory tests were conducted to ascertain the accuracy to which this SFH technique could estimate Kfs, compared with the ECH technique. This last technique is one of the simplest techniques that use the initial period of the infiltration process to deduce Kfs. An ECH experiment was conducted a few days after an SFH experiment on nine undisturbed soil cores. The mean value of Kfs obtained by the ECH technique was 1.82 times less than the corresponding value determined by the SFH technique and the associated coefficients of variation differed by a factor of 1.04. The observed discrepancies were small from a practical point of view and they were attributed to some form of surface soil degradation and closure of macropores promoted by the SFH experiment. Forty-two identical soil columns were also prepared and each technique was applied to about one half of the soil columns. Very similar estimates of both mean (i.e., differing by a factor of 1.23) and CV (i.e., differing by a factor of 1.05) of Kfs were obtained. We feel these results represent a good agreement between the ECH and SFH techniques.
A field test was conducted in a loamy soil to compare the SFH technique to the classical PI technique. The values of Kfs measured at 30 sites by the SFH technique were 2.0 times higher and 1.92 times less variable than the corresponding values obtained by the PI technique, suggesting that the two techniques yielded a similar estimate of the mean value of Kfs. The discrepancies between the SFH and the PI results were attributed to a decrease in the ability of the largest pores to conduct flow as the duration of the experiments and the sampled soil volume increased.
In all cases, the SFH technique appeared very promising since it yielded estimates of
differing by a maximum factor of two from those obtained by other techniques, both in the laboratory and in the field. In addition, the SFH technique appears to be simpler than the ECH and PI techniques used in this investigation for comparative purposes. In fact, the falling head technique is based on a simple and generally rapid experiment in the field (especially at the soil surface) and it requires standard and easily usable laboratory equipment, such as an oven and a balance. A large ring (0.200.30 m diam.) and a small volume of water (<1 L) could be used so that a larger representative elementary volume could be sampled. Vertical variations of Kfs could be assessed with a high resolution by carrying out infiltration tests at closely spaced depths, given that a small depth of soil is sampled. Especially in moderately permeable soils, many experiments could be conducted simultaneously by one person. The proposed technique appears promising for determining Kfs in a relatively short period of time without the need for extensive instrumentation or analytical methodology and therefore appears suitable for detailed sampling over large areas.
 |
ACKNOWLEDGMENTS
|
|---|
V. Bagarello and M. Iovino thank Giovanna Costa for her competent technical assistance in field data collection and Barbara Calaciura who conducted most of the laboratory experiments used in this study for her Degree Thesis. This study was supported by grants from the Italian Ministero dell'Istruzione, dell'Università e della Ricerca, Program PRIN 2000 "Impiego di modelli di simulazione nella gestione dell'irrigazione".
Received for publication February 13, 2003.
 |
REFERENCES
|
|---|
- Bagarello, V., M. Iovino, and W.D. Reynolds. 1999. Measuring hydraulic conductivity in a cracking clay soil using the Guelph permeameter. Trans. ASAE 42:957964.
- Bagarello, V., M. Iovino, and G. Tusa. 2000. Factors affecting measurement of the near saturated soil hydraulic conductivity. Soil Sci. Soc. Am. J. 64:12031210.[Abstract/Free Full Text]
- Beven, K., and P. German. 1982. Macropores and water flow in soils. Water Resour. Res. 18:13111325.
- Bosch, D.D. 1997. Constant head permeameter formula dependence on the alpha parameter. Trans. ASAE 40:13771379.
- Bouma, J. 1982. Measuring the hydraulic conductivity of soil horizons with continuous macropores. Soil Sci. Soc. Am. J. 46:438441.[Abstract/Free Full Text]
- Brooks, R.H., and A.T. Corey. 1964. Hydraulic properties of porous media. Hydrol. Paper No.3, Colorado State University, Ft. Collins.
- Elrick, D.E., and W.D. Reynolds. 1992a. Methods of analyzing constant-head well permeameter data. Soil Sci. Soc. Am. J. 56:320323.[Abstract/Free Full Text]
- Elrick, D.E., and W.D. Reynolds. 1992b. Infiltration from constant-head well permeameters and infiltrometers. p. 124. In G.C. Topp et al. (ed.) Advances in measurement of soil physical properties: Bringing theory into practice. SSSA Spec. Publ. 30. SSSA, Madison, WI.
- Elrick, D.E., G.W. Parkin, W.D. Reynolds, and D.J. Fallow. 1995. Analysis of early-time and steady-state single ring infiltration under constant head and falling head conditions. Water Resour. Res. 31:18831893.
- Elrick, D.E., R. Angulo-Jaramillo, D.J. Fallow, W.D. Reynolds, and G.W. Parkin. 2002. Analysis of infiltration under constant head and falling head conditions. p. 4753. In P.A.C. Raats et al. (ed.) Environmental mechanics: Water, mass and energy transfer in the biosphere. Geophysical Monograph Series, Vol.129, AGU, Washington, DC.
- Fallow, D.J., D.E. Elrick, W.D. Reynolds, N. Baumgartner, and G.W. Parkin. 1994. Field measurement of hydraulic conductivity in slowly permeable materials using early-time infiltration measurements in unsaturated media. p. 375389. In D.E. Daniel and S.J. Trautwein (ed.) Hydraulic conductivity and waste contaminant transport in soil. ASTM STP 1142, Am. Soc. for Testing and Materials, Philadelphia, PA.
- Green, W.H., and G.A. Ampt. 1911. Studies in soil physics. I. The flow of air and water through soils. J. Agric. Sci. 4:124.
- Hillel, D. 1998. Environmental soil physics. Academic Press, San Diego, CA.
- Lee, D.M., W.D. Reynolds, D.E. Elrick, and B.E. Clothier. 1985. A comparison of three field methods for measuring saturated hydraulic conductivity. Can. J. Soil Sci. 65:563573.
- Logsdon, S.D., and D.B. Jaynes. 1996. Spatial variability of hydraulic conductivity in a cultivated field at different times. Soil Sci. Soc. Am. J. 602:703709.
- McKenzie, N.J., H.P. Cresswell, H. Rath, and D. Jacquier. 2001. Measurement of unsaturated hydraulic conductivity using tension and drip infiltrometers. Aust. J. Soil Res. 39:823836.
- Mertens, J., D. Jacques, J. Vanderborght, and J. Feyen. 2002. Characterisation of the field-saturated hydraulic conductivity on a hillslope: In situ single ring pressure infiltrometer measurements. J. Hydrol. (Amsterdam) 263:217229.
- Mualem, Y. 1976. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 12:513522.
- Munoz-Carpena, R., C.M. Regalado, J. Alvarez-Benedi, and F. Bartoli. 2002. Field evaluation of the new Philip-Dunne permeameter for measuring saturated hydraulic conductivity. Soil Sci. 167:924.
- Odell, B.P., P.H. Groenevelt, and D.E. Elrick. 1998. Rapid determination of hydraulic conductivity in clay liners by early-time analysis. Soil Sci. Soc. Am. J. 62:5662.[Abstract/Free Full Text]
- Philip, J.R. 1957. The theory of infiltration: 4. Sorptivity and algebraic infiltration equations. Soil Sci. 84:257264.
- Philip, J.R. 1969. Theory of infiltration. Adv. Hydrosci. 5:215296.
- Philip, J.R. 1992. Falling head ponded infiltration. Water Resour. Res. 28:21472148.
- Prieksat, M.A., T.C. Kaspar, and M.D. Ankeny. 1994. Positional and temporal changes in ponded infiltration in a corn field. Soil Sci. Soc. Am. J. 58:181184.[Abstract/Free Full Text]
- Reynolds, W.D. 1993. Saturated hydraulic conductivity: Field measurement. p. 599613. In M.R. Carter (ed.) Soil sampling and methods of analysis. Canadian Soc. Soil Sci., Lewis Publishers, Boca Raton, FL.
- Reynolds, W.D., and D.E. Elrick. 1986. A method for simultaneous in situ measurement in the vadose zone of field-saturated hydraulic conductivity, sorptivity and the conductivity-pressure head relationship. Ground Water Monit. Rev. winter:8495.
- Reynolds, W.D., and D.E. Elrick. 1990. Ponded infiltration from a single ring: I. Analysis of steady flow. Soil Sci. Soc. Am. J. 54:12331241.[Abstract/Free Full Text]
- Reynolds, W.D., and W.D. Zebchuk. 1996. Hydraulic conductivity in a clay soil: Two measurement techniques and spatial characterization. Soil Sci. Soc. Am. J. 60:16791685.[Abstract/Free Full Text]
- Reynolds, W.D., B.T. Bowman, R.R. Brunke, C.F. Drury, and C.S. Tan. 2000. Comparison of tension infiltrometer, pressure infiltrometer, and soil core estimates of saturated hydraulic conductivity. Soil Sci. Soc. Am. J. 64:478484.[Abstract/Free Full Text]
- Smith, R.E. 1999. Technical note: Rapid measurement of soil sorptivity. Soil Sci. Soc. Am. J. 63:5557.[Abstract/Free Full Text]
- Soil Survey Staff. 1998. Keys to soil taxonomy. 8th ed. Pocahontas Press, Blacksburg, VA.
- van Genuchten, M.Th., F.J. Leij, and S.R. Yates. 1991. The RETC code for quantifying the hydraulic functions of unsaturated soils. U.S. Salinity Laboratory, Riverside, Ca.
- Vauclin, M., D.E. Elrick, J.L. Thony, G. Vachaud, Ph. Revol, and P. Ruelle. 1994. Hydraulic conductivity measurements of the spatial variability of a loamy soil. Soil Technol. 7:181195.
- Warrick, A.W. 1998. Spatial variability. p. 655675. In D. Hillel (ed.) Environmental Soil Physics. Academic Press, San Diego, CA.
- White, I., and M.J. Sully. 1987. Macroscopic and microscopic capillary length and time scales from field infiltration. Water Resour. Res. 23:15141522.
This article has been cited by other articles:

|
 |

|
 |
 
V. Bagarello and A. Sgroi
Testing Soil Encasing Materials for Measuring Hydraulic Conductivity of a Sandy-Loam Soil by the Cube Methods
Soil Sci. Soc. Am. J.,
June 18, 2008;
72(4):
1048 - 1057.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
W. D. Reynolds
Alternative Unsaturated Flow Analyses for the Falling-Head Ring Infiltrometer
Vadose Zone J.,
January 23, 2008;
7(1):
131 - 135.
[Abstract]
[Full Text]
[PDF]
|
 |
|