Published online 27 October 2006
Published in Soil Sci Soc Am J 70:1998-2007 (2006)
DOI: 10.2136/sssaj2006.0046
© 2006 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
Soil Physics
Tillage Effects on Hydraulic Properties and Macroporosity in Silty and Sandy Soils
U. Buczkoa,*,
O. Bensb and
R. F. Hüttlb
a Univ. of Rostock, Institute for Land Use, Justus-von-Liebig-Weg 6, D-18059 Rostock, Germany
b Brandenburg Univ. of Technology, Cottbus, Chair of Soil Protection and Recultivation, P.O. Box 101344, D-03013 Cottbus, Germany
* Corresponding author (uwe.buczko{at}uni-rostock.de)
 |
ABSTRACT
|
|---|
In agricultural soils, macroporosity and hydraulic properties are influenced by tillage practices. The objective of this study was to characterize macroporosity and surface soil hydraulic properties in two soils of different texture (Lietzen: sandy loamHumic Dystrudept; Adenstedt: silt loam; Typic Hapludoll) under conventional (CT) and conservational (RT) tillage systems. Soil hydraulic conductivity was assessed in situ by ponded infiltration with single rings (n = 70) and tension infiltration by means of a "closed-top" hood infiltrometer (HIF; n = 48). Macroporosity (pore diameters >1 mm) was estimated from differences in infiltration at saturation and at 3 cm H2O soil matric potential. Mean saturated hydraulic conductivity (Ks) for Lietzen was 3.1 x 105 m s1 and for Adenstedt was 4.3 x 105 m s1. These values are by one order of magnitude higher than values estimated from soil texture. This implies that soil structure has a dominant influence on hydraulic conductivity. Mean values of macroporosity were 0.005% for Lietzen and 0.018% for Adenstedt (using the method of Watson and Luxmoore). The respective values were 0.0008 and 0.0013% when the method of Bodhinayake et al. was used. For Adenstedt, RT showed higher macroporosity than CT (not significant at P < 0.05 for n = 12). Such treatment-induced differences were less developed for Lietzen. The Ks values measured with the ponded ring infiltrometer (RIF) at the sandy Lietzen site were higher than the corresponding values measured with the tension infiltrometer. These differences may be caused by subcritical soil water repellency (i.e., contact angles of the soil-water-air interface below 90°), although further factors could also be important (e.g., air entrapment, differences in water saturation, geometry of infiltration devices).
Abbreviations: CT, conventional tillage HIF, hood infiltrometer RIF, ring infiltrometer RT, reduced (conservation) tillage TIF, tension infiltrometer WRC, water retention characteristics
 |
INTRODUCTION
|
|---|
HYDRAULIC CONDUCTIVITY and the infiltration capacity of soils are governed to a large extent by macropores, that is, the macropore volume fraction, diameter distribution, and the continuity and connectivity of the macropore network (Ehlers, 1975; Logsdon et al., 1990; Shipitalo et al., 2000). In agricultural soils, management and tillage practices have a strong impact on the physical properties of the topsoil and the characteristics of the macropore system (e.g., Frede et al., 1994; Tebrügge and Düring, 1999; Shipitalo et al., 2000). The generic term "conventional" tillage refers to soil management systems that involve plowing and harrowing, whereas "conservation tillage" is a practice with less soil disturbance, reduced penetration depth and without inversion of the topsoil. In no-till systems, the seeds are directly drilled into the soil without any further tillage practice.
For RT and no-till, many studies reported enhanced infiltration capacity and hydraulic conductivity compared with CT (e.g., Logsdon et al., 1993; Frede et al., 1994; Reynolds et al., 1995, 2000; Tebrügge and Düring, 1999; Cameira et al., 2003). In several cases, however, lower hydraulic conductivity and infiltration rates under RT practices were reported (e.g., Culley et al., 1987; Zachmann et al., 1987; Heard et al., 1988; Unger, 1992). In other studies, no significant differences of hydraulic properties under different tillage practices were observed (e.g., Starr, 1990; Ankeny et al., 1990), or the differences of infiltration properties changed with time (Dunn and Phillips, 1991). The greater numbers of macropores and higher connectivity of the macropore system in vertical direction, which is found commonly in RT compared with CT systems is mostly attributed to greater abundances of earthworms and less disturbance of the topsoil (Ehlers, 1975; Zachmann et al., 1987; Meek et al., 1990; Dunn and Phillips, 1991; Reynolds et al., 1995).
Soil hydraulic properties can be characterized by various methods. Widely used are in situ measurements of infiltration rates, either under saturated conditions with ponded infiltration rings and pressure infiltrometers (Bond and Collis-George, 1981; Reynolds and Elrick, 1990), or tension infiltrometers under slightly negative hydraulic potentials (Ankeny et al., 1988; Perroux and White, 1988). A special kind of tension infiltrometer is the "closed-top" (Dixon, 1975; Topp and Zebchuk, 1985) or hood infiltrometer (HIF). Compared with the standard tension infiltrometer, measurements with the HIF require much less effort to prepare the soil surface and consequently less disturbance. Moreover, the effect of the contact material on measured hydraulic parameters (Reynolds and Zebchuk, 1996) is eliminated.
It is generally accepted that, especially in structured macroporous soils, field methods to evaluate hydraulic properties are preferable to laboratory methods, because the volume of soil cores for laboratory measurements is often too small to be representative, and continuous macropores are often dissected by the core walls (Angulo-Jaramillo et al., 2000). The correlation between hydraulic parameters measured with field methods and values determined using soil cores in the laboratory is often poor (Reynolds et al., 2000).
Tension disc permeameters allow estimation of the proportion which different size ranges of macropores contribute to total soil water flow by adjusting the hydraulic potential of the water supply. From the differences in infiltration rates under different hydraulic supply potentials, numbers and volume fractions of hydraulically effective macropores have been derived in agricultural (Dunn and Phillips, 1991; Trojan and Linden, 1998; Bodhinayake and Si, 2004) and forest (Watson and Luxmoore, 1986, 1988; Buttle and McDonald, 2000) soils. Use of contact sand with traditional membrane-based tension infiltrometers can cause clogging of the macropores near the soil surface. Therefore, characterization of hydraulically active macroporosity could be a strong motivation for the use of HIFs, which do not require preparation of the soil surface with contact sand. Other methods to characterize macroporosity in soils include dye tracer experiments (Andreini and Steenhuis, 1990; Droogers et al., 1998), inventory of visible macropores in the field (e.g., Logsdon et al., 1990; Ela et al., 1992), image analysis (e.g., Edwards et al., 1988), resin impregnation techniques (Singh et al., 1991), X-ray tomography (Anderson et al., 1990), calculation from laboratory-measured water retention characteristics (Carter, 1988) or from soil water contents directly beneath tension infiltrometer measurement sites (Bodhinayake and Si, 2004). While there is no universally accepted threshold diameter for macropores, many studies use a minimum value of 1 mm. The different methods to characterize macroporosity seem to yield systematically different results for macroporosity: Whereas macroporosities derived from tension infiltrometry are typically within the range of 0.001 to 0.05%, the macroporosities obtained by other methods are generally higher. Direct comparisons (Bodhinayake and Si, 2004) revealed macroporosities derived from water retention characteristics, which are by a factor of about 1000 higher than macroporosities derived from tension infiltrometry. Infiltrometry-derived macroporosity can be interpreted as hydraulically active porosity, whereas macroporosity calculated with other methods represents the total, "static" macroporosity.
The objectives of this study are twofold:
- To characterize and compare the surface saturated soil hydraulic conductivity, infiltration properties and macroporosity of two differently textured agricultural soils which are subject to long-term CT and RT sytems;
- To compare hydraulic properties determined with ponded ring infiltration measurements with data derived by a new type of "closed-top," membrane-free tension infiltrometer device, the "hood infiltrometer."
The novelty of this study compared with similar previous work by other authors is threefold: (i) application of a relatively new measurement device (HIF); (ii) comparison of tension infiltration with falling head ponded infiltration measurements; (iii) utilization and comparison of two different methods to calculate hydraulically active macroporosity based on tension infiltration data (namely, the approaches of Watson and Luxmoore, 1986, and Bodhinayake and Si, 2004).
 |
MATERIALS AND METHODS
|
|---|
Experimental Sites
The investigations were performed on two agricultural sites in Northern Germany: Lietzen in Brandenburg and Adenstedt in Lower Saxony. The most pertinent characteristics of those sites are compiled in Table 1. On both sites, different tillage practices have been applied in parallel for several years. Whereas Lietzen is a typical site for the eastern German dry sandy soils, Adenstedt is an example for the fertile soils of the central German loess belt.
View this table:
[in this window]
[in a new window]
|
Table 1. Pertinent characteristics of the Lietzen and Adenstedt sites (CT, conventional tillage; RT, conservation ["reduced"] tillage).
|
|
Soil Sampling and Laboratory Measurements
Soil sampling for texture and water retention analysis was done in April (Adenstedt) and May (Lietzen) 2000. For each site and tillage treatment, undisturbed cylindrical cores of 100 cm3 were taken according to soil horizon boundaries from 160 cm deep profile pits. Three cores per horizon were taken at the midpoint of each soil horizon. Soil water retention characteristics were determined using ceramic suction plates (Romano et al., 2002) at matric tension heads (corresponding to negative matric potentials) of 32, 63, 100, 300, and 630 cm H2O, and by means of a pressure chamber (Eijkelkamp Agrisearch Equipment BV, Giesbeek, The Netherlands) for 1000 and 15000 cm H2O (Dane and Hopmans, 2002), while the water content at saturation was determined from the difference in weight between fully saturated samples and after oven-drying at 105°C. Water content at 3-cm tension was estimated from interpolation between water contents at saturation (0 cm suction) and 32-cm tension, assuming a logarithmic scale of matric tension heads. We note that such an interpolation could possibly induce some error. The water retention data for this study were taken "as is" from a different project and it was not possible to change the pressure steps applied a posteriori. Since these data are not the main topic of this study, but serve only complementary purposes, we decided to take these data and perform the interpolation to obtain the water contents at the pressure head of interest, bearing in mind, that the values have a low level of accuracy.
Infiltration Measurements
Surface soil hydraulic properties were measured in situ between April 2000 and July 2002. Details for measurements at the Lietzen site are compiled in Table 2; at Adenstedt, all measurements were done in April/May 2000. Infiltration at near-saturated conditions was measured using a HIF, whereas ponded infiltration was measured with a single RIF. For each plot, measurements were done within an area of approximately 100 m2 and at locations unaffected by wheel traffic.
View this table:
[in this window]
[in a new window]
|
Table 2. Ks values (x 105 m s1) determined with the hood (HIF) and the ring (RIF) infiltrometers for the different dates at the Lietzen site; for most dates the values represent several data.
|
|
The HIF (Umwelt-Geräte-Technik, Müncheberg, Germany, http://www.ugt-online.de) is a modified version of a closed-top infiltrometer (Dixon, 1975; Topp and Zebchuk, 1985): Water infiltrates from a hemispherical hood (base diameter of 16 cm) with prescribed negative hydraulic heads. The circular contact line between the hood rim and the soil is sealed with medium textured sand. Steady-state infiltration rates at hydraulic pressures between 0 cm H2O and the air entry pressure of the soil system (at about 10 cm H2O) were measured. Apparent steady-state was achieved in general after about 15 min. At each measurement location, infiltration was recorded consecutively for the three hydraulic heads 5,3, and 0 cm H2O in ascending order (i.e., from unsaturated to more saturated conditions). We are aware, that the sequence of applied pressure heads can have an effect on the estimated Ks values (Bagarello et al., 2000).
The exact pressure heads applied during the measurements varied slightly, because in many instances, it was not possible to adjust the supply pressure head precisely to the predetermined value (although actually applied values could be measured with a precision of 0.1 cm). In such cases, the supply pressure head was recorded, and corresponding infiltration rates for 0 and 3 cm H2O were obtained by interpolation and extrapolation, assuming logarithmic scales of both pressure heads and infiltration rates.
Pore radii r (L) are calculated from the hydraulic head
(L) using the equation of capillarity:
 | [1] |
Here,
is the surface tension of water (M T2) (= 71.3 mN m1 at 15°C),
the contact angle between the water-air interface and the solid surface,
the density of water (M L3) (= 1000 kg m3), and g the acceleration due to gravity (L T2) (= 9.81 m s2). Assuming the pore space consists of capillary tubes, and
= 0, a hydraulic head of 3 cm H2O corresponds to a pore diameter of 1 mm, which is often adopted as a threshold diameter for macropores (Watson and Luxmoore, 1986; Trojan and Linden, 1998; Buczko et al., 2003).
However, measurements with the capillary rise method (Wöllecke, 2006) yielded contact angles distinctly greater than 0° (Table 1). For contact angles of 60°, the pore diameter corresponding to 3 cm H2O would be only 0.5 mm. However, the exact contact angles in the soil are strongly dependent on water content and probably approach 0° after prolonged infiltration. Therefore, the macropore threshold of 3 cm H2O is adopted here, and consequently, infiltration at 3 cm H2O excludes macropores, and the difference between infiltration rates at 0 cm H2O and 3 cm H2O, qm (L T1), can be used to estimate the hydraulically active macroporosity of a soil. An upper bound for the number of macropores per unit area, Nm, results from the HagenPoiseuille equation for laminar flow through a capillary tube (Watson and Luxmoore, 1986; Dunn and Phillips, 1991):
 | [2] |
Here,
denotes the dynamic viscosity of water (M L1 T1) (taken here as 0.00115 Pa s for a temperature of 15°C), and rm (L) is the minimum radius for macropores (here: 0.5 mm). Implicitly assumed in Eq. [2] is a unit hydraulic gradient, that is, steady-state conditions during infiltration flow. The effective macroporosity,
m, is estimated as:
 | [3] |
Since rm is the minimum radius of macropores,
m as calculated with Eq. [2] and [3] is an estimate of the maximum value for hydraulically active macroporosity, because rm occurs in the denominator of Eq. [2].
Recently, Bodhinayake et al. (2004) proposed a modified method for calculating the hydraulically effective macroporosity, which takes into account the hydraulic conductivity, K(
), in the pressure head interval corresponding to the two pore radii a and b. The resulting contribution of hydraulically active macroporosity for this interval of pore radii is calculated as:
 | [4] |
For known hydraulic property functions, an exact analytical solution for Eq. [4] may be obtained by direct integration. Here, Gardner's exponential model was used, and the resulting analytical solution given in Bodhinayake et al. (2004; Eq. [21]). In the present study, macroporosities calculated with both methods, Watson and Luxmoore (1986) and Bodhinayake et al. (2004), are given.
Hydraulic conductivity, K(
), was calculated from HIF data with the double tension approach according to Ankeny et al. (1991): Wooding's (1968) analytical solution for steady-state flow from a circular disc with radius rd may be written for two different hydraulic heads
1 and
2 and corresponding water volume flux rates Q1 and Q2 (L3T1) as:
 | [5] |
 | [6] |
An exponential relationship is assumed for the hydraulic conductivity function K(
) (LT1), with Ks (LT1) denoting the saturated hydraulic conductivity. The parameter
of this relation is calculated by dividing Eq. [5] by Eq. [6]:
 | [7] |
The hydraulic conductivity functions are then calculated with:
 | [8] |
and
 | [9] |
Ponded infiltration with variable ponding depths was measured using single rings (RIF) of sheet iron of 20 cm diameter, which were pressed 8 cm into the soil. Measurements were initiated with a ponded depth of about 20 cm and infiltration rates as a function of ponding depth recorded with a time resolution of 30 s, until all the water in the ring had infiltrated. Immediately thereafter, the filling and measurement procedure was repeated.
The calculation of Ks from the ponded infiltration data follows the procedure outlined in Reynolds and Elrick (1990) and Elrick et al. (1995). According to Reynolds and Elrick (1990), the water flux Q [L3 T1] from a RIF with a constant ponding depth H (L) is given by
 | [10] |
Here,
m denotes the matric flux potential (L2 T1), and G is a dimensionless geometry factor, which for H > 5 cm can be calculated as (Reynolds and Elrick, 1990)
 | [11] |
with d = depth of ring insertion; for d = 8 cm, and r = 10 cm, a value of G = 0.4368 is calculated with Eq. [11]. Based on two flux rates Q1 and Q2 (with Q2 > Q1) measured at the corresponding ponding depths H1 and H2 (H2 > H1), Ks is calculated as:
 | [12] |
and the corresponding matric flux potential
 | [13] |
This calculation procedure was developed for stationary infiltration at constant ponding depths. Further assumptions are, that the soil is homogeneous, the soil matrix is rigid, and hydraulic properties are isotropic.
A modification of this procedure for falling head conditions was developed by Elrick et al. (1995). Equation [10] may be written as
 | [14] |
Here, A is the base area of the ring (L2) and H(t) the time-variable ponding depth. Separation of variables and integration results in a forward formulation for H(t):
 | [15] |
The unknown parameters Ks and
m in Eq. [15] were estimated by an inverse calculation. An objective function was minimized which incorporates measured (H*) and predicted (H) values of ponding depths and prior estimates of Ks and
m (similarly as in Zeleke and Si, 2005):
 | [16] |
Here,
2 denotes the variance of the measurement error of H,
K2 that of Ks, and 
2 that of
m; n is the number of H(t) data pairs for each measurement sequence. The superscripts "ss" and "fh" denote "steady state" and "falling head,", respectively. Prior estimates of Ksss and
mss were obtained with Eq. [12] and [13], assuming a quasi steady state infiltration flux. The calculation was based on regression lines of ponding depths vs. infiltration rates, from which infiltration rates for ponding depths of 20, 15, 10, and 5 cm were estimated. These were subsequently used to calculate Ks and
m with Eq. [12] and [13]. The unknown variances of measurement errors were simply assumed to be 10% for H, Ks, and
m (Zeleke and Si, 2005).
Statistical Analysis
One-way analysis of variance (ANOVA) was used to test for statistical significance of differences among means of Ks and macroporosity, using the program SPSS (version 11.0.0). Comparisons between individual means were made using Fisher's Least Significant Difference when the F-value in the ANOVA was statistically significant. Customarily, the null hypothesis (equal means) is rejected at a probability P < 0.05. A test for normal distribution revealed, that macroporosity was better approximated by a log-normal distribution. For Ks values, some datasets proved to be log-normally distributed, whereas others were better approximated by a normal distribution. In ANOVA, log-transformed values were used for macroporosity, whereas both untransformed and log-transformed data were used for Ks values.
 |
RESULTS AND DISCUSSION
|
|---|
Hydraulic Conductivity
For the Lietzen site, Ks values (Fig. 1
, Table 3) are slightly higher when determined with the RIF in comparison with the HIF values and higher for conventional tillage (CT) than for conservation (reduced) tillage (RT). None of these differences is, however, significant at P < 0.05 (Table 3). The overall geometric mean of Ks values for all measurements at the Lietzen site is significantly lower than the overall geometric mean of Ks values at the Adenstedt site (Table 3).

View larger version (12K):
[in this window]
[in a new window]
|
Fig. 1. Boxplots of Ks values for the conventional and conservation tillage plots of the Lietzen site, determined with the hood (HIF) and the ring (RIF) infiltrometer. The numbers of measurements are given in Table 3. Boxes denote the 25th and 75th percentiles positions, the continuous bar inside the box shows the median and the dash-dotted line the mean value; the whiskers denote the 10th and 90th percentiles; outliers are indicated as triangles. This notation applies also to subsequent box plots.
|
|
View this table:
[in this window]
[in a new window]
|
Table 3. Summary statistics of Ks values (x 105 m s1) determined with the hood (HIF) and the ring (RIF) infiltrometers for the Lietzen and Adenstedt sites (CT: conventional tillage; RT: conservation ["reduced"] tillage).
|
|
In comparison with Ks values calculated with pedotransfer functions based on soil texture, the mean Ks values for both the Lietzen and the Adenstedt sites seem rather high: an estimation with the neural network model "Rosetta" (Schaap et al., 2001), based solely on the soil textural class, revealed for the sandy loam soil at Lietzen a Ks value 4.4 x 106 m s1 and for the silt loam at Adenstedt 2.1 x 106 m s1, that is, values about one order of magnitude lower than actually measured. These differences are probably due to the influence of soil structure (earthworm burrows, root channels) and plant residues. This is not accounted for in the pedotransfer prediction model. The greater differences between measured and estimated Ks values for the Adenstedt compared with the Lietzen site probably reflect the more dominant and better developed soil structure at Adenstedt (Hangen et al., 2002; Buczko et al., 2003). Similar differences between Ks values measured in situ and predicted with pedotransfer functions are described by Bodhinayake and Si (2004) for a silt loam soil under grass. A possible further factor for the explanation of the differences between Ks values predicted with pedotransfer functions and the values actually measured in situ could be the requirement of steady-state flow conditions in the analysis. This might not be attained during the relatively short measurement times of 20 min, as used in the HIF measurements. Especially for finer textured soils, true steady-state flow conditions are attained often only after longer times, sometimes several hours (Elrick et al., 1990). A measurement time of 20 min was chosen for practical reasons. Numerical simulations (not shown here in detail) indicated, that the error when using apparent steady-state infiltration rates after 20 min is mostly relatively low (<20%, depending on hydraulic properties and initial water contents), and at any rate could not explain the observed large differences between predicted and measured Ks values.
In contrast to the Lietzen data, at the Adenstedt site, Ks values derived from RIF measurements are not consistently higher than Ks values derived from HIF measurements (mean values are lower for RIF-Ks than for HIF-Ks values; Table 3; Fig. 2
). For the Lietzen site, the geometric mean of RIF-Ks values is 1.6 times higher than the HIF-Ks data. The different relations between Ks values derived from HIF and RIF measurements at the two different sites may be caused by slight differences in soil water repellency, and concomitantly, different degrees of wetting during HIF measurements. For both sites, soil water repellency measured with the water drop penetration time (WDPT) and the ethanol percentage (EP) test was very low (WDPT values were consistently below 5 s, and EP values 0%). However, measurements with the capillary rise method revealed "subcritical" repellency levels (i.e., contact angles of the solid-liquid-gas interface below 90°) with contact angles of about 60° (Wöllecke, 2006; Table 1). For the RT plot at the Adenstedt site, the measured contact angle is clearly lower (about 30°). However, the contact angle data given in Table 1 are based on only few soil samples, and a statistical treatment was not possible. A further indication of a slightly higher repellency level at the Lietzen site is given by the coarser texture compared with the Adenstedt site and the dryer climatic conditions (Table 1). Both factors, dryer climate and coarser texture, favor the formation of soil water repellency (Doerr et al., 2000). It is well known that in severely water repellent soils (i.e., contact angles > 90°) infiltration with positive hydraulic heads is necessary to overcome the initial resistance to wetting and that the hydraulic head corresponds to the water entry pressure which is directly related to the repellency level (Wang et al., 2000).

View larger version (13K):
[in this window]
[in a new window]
|
Fig. 2. Boxplots of Ks values for the conventional and conservation tillage plots of the Adenstedt site, determined with the hood (HIF) and the ring (RIF) infiltrometer. The numbers of measurements are given in Table 3.
|
|
Similar infiltration measurements at the forest site Kahlenberg, which exhibits in parts severe water repellency (Buczko et al., 2002, 2005; Buczko and Bens, 2006), revealed distinctly higher differences between Ks values derived from RIF measurements and values calculated using HIF data (Buczko et al., 2006). Higher Ks values estimated with ponded infiltration measurements compared with tension infiltration data are also reported by other authors: Reynolds et al. (2000) reported in a systematic comparison of different methods mostly higher RIF-Ks compared with TIF-Ks values, especially in sandy soils and under forest vegetation. Bodhinayake and Si (2004) reported for silt loam soils under grass RIF-Ks values which are higher by a factor of about 2 than TIF-Ks values. On the other hand, Bagarello et al. (2000) found in a sandy loam soil under citrus no significant differences in Ks values between both measurement methods. In severely repellent soils with contact angles greater than 90° the differences might be explained by incomplete and spatially heterogeneous wetting. In the agricultural soils investigated here, contact angles were <90° and the soils were thoroughly wetted during infiltration measurements. However, it has been shown by various authors, that subcritical repellency levels persist also in wetted soils and that they can have a distinct influence on infiltration and soil water flow (e.g., Lamparter et al., 2006). An influence of subcritical soil water repellency with contact angles <90° on soil hydrological processes is also indicated by the capillary rise Eq. [1].
An exact correlation between RIF- and HIF- Ks values at the Lietzen and Adenstedt sites is problematical, because corresponding RIF and HIF measurements were done at neighboring (about 50 cm apart), but not absolutely identical positions. Between these neighboring datasets of RIF and HIF Ks values there is no correlation (not shown here). Similar poor correlations between Ks values determined with ponded and tension infiltration measurements are reported by Reynolds et al. (2000).
Concerning Ks values of different tillage treatments, only for the Adenstedt site values for the RT plot are significantly higher (by a factor of about 2) than values at the CT plot. However, these differences are not significant when only RIF-Ks values are considered (Table 3, Fig. 2). Similarly, higher Ks values under RT compared with CT was described for many sites worldwide: Cameira et al. (2003) reported for a silt loam soil under corn higher Ks values (measured with TIF) under RT than under CT and Reynolds et al. (2000) reported for soils with loamy and clay-loamy texture higher Ks values at no-till sites compared with CT. According to Logsdon et al. (1993), ponded infiltration rates for RT under corn in loam soil are significantly higher than infiltration rates in CT plots. Under trafficked treatment, Ankeny et al. (1990) described for a silty clay loam soil higher TIF infiltration rates for no-till compared with CT.
The lack of significant differences in hydraulic properties between CT and RT at the Lietzen site (Table 3) can possibly be attributed to the short duration of the tillage experiment (since 1996; Table 1) at the time of the studies, and in general to the lower earthworm activity in this sandy soil compared with the Adenstedt site (Hangen et al., 2002).
A factor not accounted for in the present analysis is temporal variability of hydraulic conductivity: for many agricultural soils, seasonal variations of Ks values and infiltration rates have been described (Angulo-Jaramillo et al., 1997; Heddadj and Gascuel-Odoux, 1999; Cameira et al., 2003; Bagarello and Sgroi, 2004). Temporal variations of Ks were also observed at the Lietzen site (see Table 2), but the observed temporal trends were inconclusive and too few measurements were available for each date for a proper statistical treatment of the data. Moreover, the amplitude of temporal variation seemed low enough to justify a bulk treatment of all the data measured throughout the course of the year.
Also, spatial variability of hydraulic properties is often observed at agricultural sites (Logsdon and Jaynes, 1996; Mohanty et al., 1996a). Since the exact measurement locations at the Lietzen and Adenstedt sites were slightly different at the different dates, spatial variability probably plays a role, but could not be quantified and was superimposed by temporal variability.
Macroporosity
Values of hydraulically active macroporosities (Table 4, Fig. 3
) range for both sites largely between 0.001 and 0.2% when using the calculation procedure of Watson and Luxmoore (1986). These values are similar to those determined in many other studies (see Table 5). When the procedure of Bodhinayake et al. (2004) is used, macroporosities are by a factor of 5 to 10 lower (Table 4). Similarly large differences for macroporosities estimated with the two methods are reported by Bodhinayake et al. (2004). For better comparability with results obtained by other authors in previous studies (Table 5), we focus here on the values obtained with the method of Watson and Luxmoore (1986), bearing in mind the limitations of this method (see Bodhinayake et al., 2004). For Lietzen, the geometric mean macroporosity is 0.005%, whereas it is significantly higher at Adenstedt (0.018%). This is in line with a more abundant soil fauna (mainly earthworms) observed and the higher macropore continuity determined with dye-tracer experiments at Adenstedt (Hangen et al., 2002). For both sites, hydraulically active macroporosities are higher for the RT treatment (although not significant at P < 0.05, Table 4). For Adenstedt, macroporosities estimated from HIF infiltration data are similar to values obtained from classical tension infiltrometer measurements (Buczko et al., 2003; Table 6).

View larger version (16K):
[in this window]
[in a new window]
|
Fig. 3. Boxplots of macroporosity values for the conventional and conservation tillage plots of the Adenstedt and Lietzen site, determined with the hood infiltrometer (HIF). The numbers of measurements are given in Table 4. For each study plot, the left side box plot denotes macroporosity values calculated with the approach of Watson and Luxmoore (1986), and the right side box plot denotes the values calculated with the method of Bodhinayake et al. (2004).
|
|
Water retention characteristics (WRC) (Fig. 4
) are similar for the different tillage treatments at each site but distinctly different between both sites. For all soil water matric tension heads, water storage is higher at Adenstedt. Macroporosities calculated from WRC (Table 6) are higher by a factor of about 150 to 350 than values estimated from HIF infiltration data. For Adenstedt, they are in the same range as macroporosities estimated from dye-tracer experiments and inventories of visible unstained pores (Buczko et al., 2003). Similarly large differences between macroporosities estimated with these two (or similar) methods have been described also by other authors (Messing and Jarvis, 1993; Buttle and McDonald, 2000; Bodhinayake and Si, 2004). Both these indicators of macroporosity have different physical significance: whereas macroporosities estimated from differences in water contents at two different tensions represent total static macroporosities, values estimated from unsaturated infiltration measurements at different soil water tensions can be interpreted as hydraulically active macroporosities. At the pore scale, the differences between both types of macroporosities may be caused by pore necking, tortuosity and dead-end pores (Messing and Jarvis, 1993). At the core scale, a factor may be hysteresis of the water retention curve: whereas the WRC measured on soil cores in the laboratory represents the draining curve, unsaturated infiltration measurements represent the wetting curve. The ratio of macroporosities derived from tension infiltrometer measurements and that derived from static water content differences can be used as an indicator of macropore connectivity (Bodhinayake and Si, 2004). For the Adenstedt site, the macropore connectivity is clearly higher for RT compared with CT, whereas for Lietzen macropore connectivity is more or less equal for both tillage treatments (Table 6). Similar values of macropore connectivity were reported by Bodhinayake and Si (2004) for silt loam soils under grass, whereas Bodhinayake and Si (2004) described lower macropore connectivities in cultivated soils.

View larger version (13K):
[in this window]
[in a new window]
|
Fig. 4. Measured water retention curves (0- to 10-cm soil depth) for the conventional and conservation tillage plots of the Adenstedt and Lietzen sites. The number of replications for each suction head and treatment/plot is n = 3. For better clearness of this graphic, error bars are not given.
|
|
 |
CONCLUSIONS
|
|---|
Saturated hydraulic conductivity (Ks), measured with ponded ring infiltration and hood infiltration measurements, revealed for the sandy Lietzen site an overall geometric mean value of 3.1 x 105 m s1 and for the Adenstedt silt loam 4.3 x 105 m s1. These Ks values are distinctly (by a factor of about 10) higher than values estimated with pedotransfer functions based on the textural class. Consequently, soil structure has presumably a strong influence on hydraulic conductivity in these soils, especially in the Adenstedt silt loam. Comparing the different tillage treatments, only at the finer textured Adenstedt site was Ks significantly higher for the RT than for the CT treatment, whereas at the sandy loam soil site Lietzen, differences in Ks with tillage treatment were not significant. This may be due to the higher biological activity and the better developed soil structure at Adenstedt compared with Lietzen. Using the method of Watson and Luxmoore (1986), geometric mean values of hydraulically active macroporosity were only 0.005% for Lietzen, but 0.018% for the Adenstedt site, whereas the respective values were 0.0008 and 0.0013% when the method of Bodhinayake et al. (2004) was used. For Adenstedt, the RT treatment showed higher hydraulically active macroporosity and higher macropore connectivity than the CT treatment. Such differences in macroporosity for different tillage treatments were far less developed for Lietzen. For the Lietzen site, Ks measured with ponded ring infiltration yielded higher values than measurements with tension infiltration. The differences in Ks between the infiltration methods may be caused by subcritical soil water repellency (i.e., contact angles of the soil-water-air interface below 90°). However, other factors may also be important (e.g., air entrapment, differences in water saturation, geometry of infiltration devices). In summary, the differences in hydraulic properties (Ks, macroporosity) induced by different tillage treatments were clearly visible only for the finer textured Adenstedt site, but hardly apparent for the sandy Lietzen site. This is presumably caused to a high degree by the higher activity of soil fauna and concomitant better developed soil structure at the Adenstedt site.
 |
ACKNOWLEDGMENTS
|
|---|
This investigation was financially supported by the German Ministry of Education and Research (BMBF) as subproject A 4.3 of the German Research Network Natural Disasters (DFNK) under grant number 01SF9971/4. Thanks to B. Wöllecke (Cottbus) for assistance with soil sampling and water retention data, and to Dr. N.A. Wahl (Univ. Aalborg, Denmark) for assistance with infiltration measurements.
Received for publication February 1, 2006.
 |
REFERENCES
|
|---|
- Anderson, S.H., R.L. Peyton, and C.J. Gantzer. 1990. Evaluation of constructed and natural soil macropores using x-ray computed tomography. Geoderma 46:1329.[CrossRef][ISI]
- Andreini, M.S., and T.S. Steenhuis. 1990. Preferential paths of flow under conventional and conservational tillage. Geoderma 46:85102.[CrossRef][ISI]
- Angulo-Jaramillo, R., F. Moreno, B.E. Clothier, J.L. Thony, G. Vachaud, E. Fernandez-Boy, and J.A. Cayuela. 1997. Seasonal variation of hydraulic properties of soils measured using tension disk infiltrometer. Soil Sci. Soc. Am. J. 61:2732.[ISI]
- Angulo-Jaramillo, R., J.-P. Vandervaere, S. Roulier, J.L. Thony, J.-P. Gaudet, and M. Vauclin. 2000. Field measurement of soil surface hydraulic properties by disc and ring infiltrometersA review and recent developments. Soil Tillage Res. 55:129.[CrossRef]
- Ankeny, M.D., T.C. Kaspar, and R. Horton. 1988. Design for an automated tension infiltrometer. Soil Sci. Soc. Am. J. 52:893896.[Abstract/Free Full Text]
- Ankeny, M.D., T.C. Kaspar, and R. Horton. 1990. Characterization of tillage and traffic effects on unconfined infiltration measurements. Soil Sci. Soc. Am. J. 54:837840.[Abstract/Free Full Text]
- Ankeny, M.D., T.C. Kaspar, and R. Horton. 1991. Simple field method for determining unsaturated hydraulic conductivity. Soil Sci. Soc. Am. J. 55:467470.[ISI]
- 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]
- Bagarello, V., and A. Sgroi. 2004. Using the single-ring infiltrometer method to detect temporal changes in surface soil field-saturated hydraulic conductivity. Soil Tillage Res. 76:1324.[CrossRef]
- Bodhinayake, W., and B.C. Si. 2004. Near-saturated surface soil hydraulic properties under different land uses in the St Denis National Wildlife Area, Saskatchewan, Canada. Hydrol. Processes 18:28352850.[CrossRef]
- Bodhinayake, W., B.C. Si, and C. Xiao. 2004. New method for determining water-conducting macro- and mesoporosity from tension infiltrometer. Soil Sci. Soc. Am. J. 68:760769.[Abstract/Free Full Text]
- Bond, W.J., and N. Collis-George. 1981. Ponded infiltration into simple soil systems: 3. The behavior of infiltration rate with time. Soil Sci. 131:327333.
- Buczko, U., O. Bens, H. Fischer, and R.F. Hüttl. 2002. Water repellency in sandy luvisols under different forest transformation stages in Northeast-Germany. Geoderma 109:118.[CrossRef][ISI]
- Buczko, U., O. Bens, E. Hangen, J. Brunotte, and R.F. Hüttl. 2003. Infiltration and macroporosity of a silt loam soil under two contrasting tillage systems. Landbauforschung Völkenrode 53:181190.
- Buczko, U., O. Bens, and R.F. Hüttl. 2005. Variability of soil water repellency in sandy forest soils with different stand structure under Scots pine (Pinus sylvestris) and beech (Fagus sylvatica). Geoderma 126:317336.[CrossRef][ISI]
- Buczko, U., O. Bens, and R.F. Hüttl. 2006. Water infiltration and hydrophobicity in forest soils of a pine-beech transformation chronosequence. J. Hydrol. In press.
- Buczko, U., and O. Bens. 2006. Assessing soil hydrophobicity and its variability through the soil profile using two different methods. Soil Sci. Soc. Am. J. 70:718727.[Abstract/Free Full Text]
- Buttle, J.M., and D.J. McDonald. 2000. Soil macroporosity and infiltration characteristics of a forest podzol. Hydrol. Processes 14:831848.[CrossRef]
- Cameira, M.R., R.M. Fernando, and L.S. Pereira. 2003. Soil macropore dynamics affected by tillage and irrigation for a silty loam alluvial soil in southern Portugal. Soil Tillage Res. 70:131140.[CrossRef]
- Carter, M.R. 1988. Temporal variability of soil macroporosity in a fine sandy loam under mouldboard ploughing and direct drilling. Soil Tillage Res. 12:3751.
- Culley, J.L.B., W.E. Larson, and G.W. Randall. 1987. Physical properties of a typic Haplaquoll under conventional and no-tillage. Soil Sci. Soc. Am. J. 51:15871593.[Abstract/Free Full Text]
- Dane, J.H., and J.W. Hopmans. 2002. Pressure cell. P. 684688. In J.H. Dane and G.C. Topp, (ed.). Methods of soil analysis. Part 4. SSSA Book Series No. 5, SSSA, Madison, WI.
- Dixon, R.M. 1975. Design and use of closed-top infiltrometers. Soil Sci. Soc. Am. Proc. 39:755763.
- Doerr, S.H., R.A. Shakesby, and R.P.D. Walsh. 2000. Soil water repellency: Its causes, characteristics and hydrogeomorphological significance. Earth-Science Rev. 51:3365.
- Droogers, P., A. Stein, J. Bouma, and G. de Boer. 1998. Parameters for describing soil macroporosity derived from staining patterns. Geoderma 83:293308.[CrossRef][ISI]
- Dunn, G.H., and R.E. Phillips. 1991. Macroporosity of a well-drained soil under no-till and conventional tillage. Soil Sci. Soc. Am. J. 55:817823.[Abstract/Free Full Text]
- Edwards, W.M., L.D. Norton, and C.E. Redmond. 1988. Characterizing macropores that affect infiltration into nontilled soil. Soil Sci. Soc. Am. J. 52:482487.
- Ehlers, W. 1975. Observations of earthworm channels and infiltration on tilled and untilled loess soil. Soil Sci. 119:242249.
- Ela, S.D., S.C. Gupta, and W.J. Rawls. 1992. Macropore and surface seal interactions affecting water infiltration into soil. Soil Sci. Soc. Am. J. 56:714721.[ISI]
- Elrick, D.E., W.D. Reynolds, H.R. Geering, and K.A. Tan. 1990. Estimating steady infiltration rate times for infiltrometers and permeameters. Water Resour. Res. 26:759769.[CrossRef]
- 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 falling head conditions. Water Resour. Res. 31:18831893.[CrossRef]
- Frede, H.G., R. Beisecker, and S. Gäth. 1994. Long-term impacts of tillage on the soil ecosystem. J. Plant Nutr. Soil Sci. 157:197203.
- Hangen, E., U. Buczko, O. Bens, J. Brunotte, and R.F. Hüttl. 2002. Infiltration patterns into two soils under conventional and conservation tillage: Influence of plant structure and soil animal activity. Soil Tillage Res. 63:181186.[CrossRef]
- Heard, J.R., E.J. Kladivko, and J.V. Mannering. 1988. Soil macroporosity, hydraulic conductivity and air permeability of silty soils under long term conservation tillage in Indiana. Soil Tillage Res. 11:118.
- Heddadj, D., and C. Gascuel-Odoux. 1999. Topographic and seasonal variations of unsaturated hydraulic conductivity as measured by tension disc infiltrometers at the field scale. Eur. J. Soil Sci. 50:275283.
- Lamparter, A., M. Deurer, J. Bachmann, and W.H.M. Duijnisveld. 2006. Effect of subcritical hydrophobicity in a sandy soil on water infiltration and mobile water content. J. Plant Nutr. Soil Sci. 169:3846.[CrossRef]
- 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. 60:703709.[Abstract/Free Full Text]
- Logsdon, S.D., R.R. Allmaras, L. Wu, J.B. Swan, and G.W. Randall. 1990. Macroporosity and its relation to saturated hydraulic conductivity under different tillage practices. Soil Sci. Soc. Am. J. 54:10961101.[Abstract/Free Full Text]
- Logsdon, S.D., J.L. Jordahl, and D.L. Karlan. 1993. Tillage and crop effects on ponded and tension infiltration rates. Soil Tillage Res. 28:179189.[CrossRef]
- Meek, B.D., W.R. DeTar, D. Rolph, E.R. Rechel, and L.M. Carter. 1990. Infiltration rate as affected by an alfalfa and no-till cotton cropping system. Soil Sci. Soc. Am. J. 54:505508.[Abstract/Free Full Text]
- Messing, I., and N.J. Jarvis. 1993. Temporal variation in the hydraulic conductivity of a tilled clay soil as measured by tension infiltrometers. J. Soil Sci. 44:1124.
- Mohanty, B.P., M.D. Ankeny, R. Horton, and R. Kanwar. 1996a. Spatial analysis of hydraulic conductivity measured using disk infiltrometers. Water Resour. Res. 30:24892498.
- Mohanty, B.P., R. Horton, and M.D. Ankeny. 1996b. Infiltration and macroporosity under row crop agricultural field in a glacial till soil. Soil Sci. 161:205213.[CrossRef]
- Perroux, K.M., and I. White. 1988. Designs for disc permeameters. Soil Sci. Soc. Am. J. 52:12051215.[Abstract/Free Full Text]
- 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.[ISI]
- Reynolds, W.D., E.G. Gregorich, and W.E. Curnoe. 1995. Characterisation of water transmission properties in tilled and untilled soils using tension infiltrometers. Soil Tillage Res. 33:117131.[CrossRef]
- 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]
- Reynolds, W.D., and W.D. Zebchuk. 1996. Use of contact material in tension infiltrometer measurements. Soil Technol. 9:141159.
- Romano, N., J.W. Hopmans, and H. Dane. 2002. Suction table. p. 392398. In J.H. Dane, G.C. Topp, (ed.). Methods of soil analysis. Part 4. SSSA Book Series No.5. SSSA, Madison, WI.
- Schaap, M.G., F.J. Leij, and M.Th. van Genuchten. 2001. Rosetta: A computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions. J. Hydrol. 251:163176.[CrossRef]
- Shipitalo, M.J., W.A. Dick, and W.M. Edwards. 2000. Conservation tillage and macropore factors that affect water movement and the fate of chemicals. Soil Tillage Res. 53:167183.[CrossRef]
- Singh, P., R.S. Kanwar, and M.L. Thompson. 1991. Macropore characterization for two tillage systems using resin-impregnation techniques. Soil Sci. Soc. Am. J. 55:16741679.[Abstract/Free Full Text]
- Starr, J.L. 1990. Spatial and temporal variation of ponded infiltration. Soil Sci. Soc. Am. J. 54:629636.[Abstract/Free Full Text]
- Tebrügge, F., and R.A. Düring. 1999. Reducing tillage intensity: A review of results from a long-term study in Germany. Soil Tillage Res. 53:1528.
- Topp, G.C., and W.D. Zebchuk. 1985. A closed adjustable head infiltrometer. Canad. Agric. Eng. 27:99104.
- Trojan, M.D., and D.R. Linden. 1998. Macroporosity and hydraulic properties of earthworm-affected soils as influenced by tillage and residue management. Soil Sci. Soc. Am. J. 62:16871692.[Abstract/Free Full Text]
- Unger, P.W. 1992. Infiltration of simulated rainfall: Tillage system and crop residue effects. Soil Sci. Soc. Am. J. 56:283289.[Abstract/Free Full Text]
- Wang, Z., L. Wu, and Q.J. Wu. 2000. Water-entry value as an alternative indicator of soil water-repellency and wettability. J. Hydrol. 231232:7683.[CrossRef]
- Wang, D., B. Lowery, J.M. Norman, and K. McSweeney. 1996. Ant burrow effects on water flow and soil hydraulic properties of sparta sand. Soil Tillage Res. 37:8393.[CrossRef]
- Watson, K.W., and R.J. Luxmoore. 1986. Estimating macroporosity in a forest watershed by use of tension infiltrometer. Soil Sci. Soc. Am. J. 50:578582.[Abstract/Free Full Text]
- Wilson, G.V., and R.J. Luxmoore. 1988. Infiltration, macroporosity, and mesoporosity distributions on two forested watersheds. Soil Sci. Soc. Am. J. 52:329335.[Abstract/Free Full Text]
- Wöllecke, B. 2006. Einfluß von Acker- und Waldnutzung auf Wasseraufnahme und Wasserspeicherung von Löss- und Sandböden. Ph.D. Diss. Univ. Cottbus.
- Wooding, R.A. 1968. Steady infiltration from a shallow circular pond. Water Resour. Res. 4:12591273.
- Zachmann, J.E., D.R. Linden, and C.E. Clapp. 1987. Macroporous infiltration and redistribution as affected by earthworms, tillage and residue. Soil Sci. Soc. Am. J. 51:15801586.[Abstract/Free Full Text]
- Zeleke, T.B., and B.C. Si. 2005. Parameter estimation using the falling head infiltration model: Simulation and field experiment. Water Resour. Res. 41:W02027 doi:10.1029/2004WR003407.[CrossRef]
This article has been cited by other articles:

|
 |

|
 |
 
H. Blanco-Canqui and R. Lal
Impacts of Long-Term Wheat Straw Management on Soil Hydraulic Properties under No-Tillage
Soil Sci. Soc. Am. J.,
June 8, 2007;
71(4):
1166 - 1173.
[Abstract]
[Full Text]
[PDF]
|
 |
|