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


     


Published online 12 March 2007
Published in Soil Sci Soc Am J 71:500-506 (2007)
DOI: 10.2136/sssaj2006.0128
© 2007 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Moret, D.
Right arrow Articles by Arrúe, J. L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Moret, D.
Right arrow Articles by Arrúe, J. L.
Agricola
Right arrow Articles by Moret, D.
Right arrow Articles by Arrúe, J. L.
Related Collections
Right arrow Hydraulic Conductivity
Right arrow Soil Methods/Instrumentation
Right arrow Structure and Properties

SOIL PHYSICS

Characterizing Soil Water-Conducting Macro and Mesoporosity as Influenced by Tillage Using Tension Infiltrometry

D. Moret* and J. L. Arrúe

Departamento de Suelo y Agua, Estación Experimental de Aula Dei, Consejo Superior de Investigaciones, Científicas (CSIC), PO Box 202, 50080 Zaragoza, Spain

* Corresponding author (david{at}eead.csic.es).


    ABSTRACT
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Knowledge of soil water conductive porosity is of paramount importance for understanding the water and solute movement in soil. The objective of this study was to describe and evaluate an alternative pore index, i.e., the representative mean pore radius for two consecutive soil water tensions ({lambda}{Delta}{Psi}), for characterizing the water-transmitting macro- and mesoporosity of soil. The hydraulic conductivity (K) and related hydraulic parameters were measured using a tension disk infiltrometer at –14, –4, –1, and 0 cm water pressure heads at a soil depth of 2 cm on a long-term conservation tillage experiment (after 9 yr of trial) under conventional tillage (CT), reduced tillage (RT), and no-tillage (NT) treatments. The soil was loam (Xerollic Calciorthid). The measurements were performed within the first half of a 16- to 18-mo-long fallow period. Unlike the model based on the classical capillary rise theory, which assumes no differences between measurement sites (e.g., tillage treatments) for the corresponding equivalent pore radius (C0) for macro- and mesopores, significant differences in {lambda}{Delta}{Psi} values for macro- and mesopores were found between the three tillage systems tested. The number of water-transmitting macro- and mesopores per unit area estimated by means of the capillary rise approach was significantly greater than that calculated using the {lambda}{Delta}{Psi} pore index. Although significantly higher values of {lambda}{Delta}{Psi} for both macro- and mesopores were observed under NT, a significantly lower near-saturation K was found under NT than CT and RT treatments due to the lower density of water-transmitting macro- and mesopores.

Abbreviations: CT, Conventional tillage • NT, no-tillage • RT, reduced tillage


    INTRODUCTION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The measurement of soil hydraulic properties, which includes the quantification of macropores and preferential flow, is of paramount importance for many soil-related studies involving disciplines such as agriculture, forestry, and hydrology. The characterization of these properties is difficult because of the fragile and transient nature of the soil macropores and the lack of appropriate measurement techniques (Messing and Jarvis, 1993).

In recent decades, the tension disk infiltrometer (Perroux and White, 1988) has become a very popular and valuable device for studying saturated and near-saturated soil hydraulic properties. This is a relatively rapid technique that can be applied in situ. It operates in the near-zero soil water pressure head range, where the soil pores are highly hydraulically active in the transmission of water and solutes (Ankeny et al., 1991). Disk infiltrometer measurements make it possible to calculate the hydraulic conductivity, K, and sorptivity, S (Perroux and White, 1988), and estimate the flow-weighted mean pore radius, {lambda}{psi} (White and Sully, 1987), of undisturbed soils. Watson and Luxmoore (1986) used the disk infiltrometer to estimate the concentration of effective macro- and mesopores on the soil surface. In this case, the minimum equivalent pore radius (C0) for a given tension range was used to compute the maximum number of effective pores per square meter using Poiseuille's Law for flow in a capillary tube. This approach, however, involves an inconsistency because K is related to the range of pore sizes participating in water transmission, while C0 relates to the pore size of water storage for a given tension under static conditions (Reynolds et al., 1995). Accordingly, Reynolds et al. (1995) proposed using the representative mean pore radius, {lambda}{psi}, instead of C0 to estimate the number of water-conductive pores per unit area. More recently, Bodhinayake et al. (2004) have presented a new general equation for water-conducting porosity and derived analytical solutions for the water-conducting porosity based on ponded- and tension-infiltration measurements in conjunction with four commonly used hydraulic conductivity–pressure head functions.

Soil water infiltration is directly affected by the soil tillage management. Tillage and the ensuing compaction can drastically modify the total "hydraulically active" porosity. The tension disk infiltrometer is a valuable tool for studying the effect of soil tillage practices on the soil surface hydraulic properties (Malone et al., 2003). Chan and Heenan (1993) and McGarry et al. (2000), among others, have reported higher infiltration rates under no-tillage than tilled treatments due to the greater number of macropores (Logsdon et al., 1990). Other researchers, by contrast, have found similar (Sauer et al., 1990) or lower (Miller et al., 1998; Evett et al., 1999) values of hydraulic conductivity under the no-tillage treatment. In general, the water flow of structured soils is mainly conducted by macropores even though these constitute only a very small fraction of the total porosity (Sauer et al., 1990; Reynolds et al., 1995; Angulo-Jaramillo et al., 1997; Cameira et al., 2003).

The objective of this study was to propose and evaluate an alternative pore index derived from {lambda}{psi}, the equivalent mean pore radius for two consecutive soil water tensions, {lambda}{Delta}{psi}, for characterizing and quantifying the water-conductive soil macro- and mesoporosity. To this end, measurements of soil hydraulic conductivity using the disk infiltrometry technique and the corresponding water-conductive pore analysis were performed under three different tillage systems on a long-term conservation tillage experiment in a dryland semiarid cereal-growing field in central Aragon.


    MATERIALS AND METHODS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Theory: Soil Water Conductive Porosity
A steady flow from a tension disk infiltrometer can be described by the simplified algebraic expression (Wooding, 1968)

Formula 1[1]
where Q{psi} (m3 s–1) is the steady-state infiltrating flux when a constant pressure head, {psi} (m of water), is applied to the tension infiltrometer membrane, R (m) is the radius of the infiltration soil surface, K{psi} (m s–1) is the hydraulic conductivity, and {phi}{psi} (m2 s–1) is the matric flux potential (Gardner, 1958) for the infiltration surface, which is defined by

Formula 2[2]
where {psi}i (m) is the reference soil water potential.

The maximum equivalent pore radius that remains full of water under a given applied pore water pressure {psi}, C0, is defined by the classical capillary rise theory as

Formula 3[3]
where {sigma} (g s–2) is the surface tension of water, {rho} (g cm–3) is the density of water, and g (cm s–2) is the acceleration due to gravity. In our study, we defined macropores as those pores that drain at {psi} > –4 cm (C0 > 0.375 mm; Clothier and White, 1981) and mesopores as those pores draining at {psi} between –4 and –14 cm (0.375 > C0 > 0.107 mm).

White and Sully (1987) defined the representative mean pore radius, {lambda}{psi}, as

Formula 4[4]

This pore index represents an effective "equivalent mean" pore radius that conducts water when infiltration occurs at a given value of {psi}. However, Reynolds and Elrick (2005) have argued that {lambda}{psi} cannot be an actual physical pore size, but a kind of index of "water conductiveness" that relates primarily to flow impedance and only secondarily to pore size or volume. The representative mean pore radius, {lambda}{psi}, can also be expressed as

Formula 5[5]
where {alpha}{psi}* = K{psi}/{phi}{psi} is the slope of the lnK vs.{psi} curve (Ankeny, 1992). The constant {alpha}{psi}* value between adjacent {psi} settings is one of the main assumptions used in steady-flow analysis using tension infiltrometry (Reynolds and Elrick, 1991; Ankeny et al., 1991). When {alpha}{psi}* is expressed for two consecutive soil water tensions, {Delta}{psi} = {psi}i{psi}i–1, the representative mean pore radius (Eq. [5]) defined for {Delta}{psi}, {lambda}{Delta}{psi}, can be expressed as

Formula 6[6]
where n is the number of measurements performed in a sequence. Note that {lambda}{Delta}{psi} defines a kind of pore index of "water conduciveness" that relates to the flow impedance for a soil tension range corresponding to a specific "pore size."

Using Poiseuille's Law for flow in a capillary tube, Watson and Luxmoore (1986) defined the maximum number of effective macro- and mesopores per unit area, NC0, as

Formula 7[7]
where C0 is the minimum pore radius for the given soil tension range, and µ (g cm–1 s–1) is the dynamic viscosity of water. They likewise defined the effective porosity, {theta}C0, as

Formula 8[8]
However, as pointed out by Reynolds et al. (1995), this approach introduces a contradiction because K in Eq. [7] is related to the range of pore sizes participating in water transmission, while C0 relates to the minimum pore sizes of water storage for a given tension range under static conditions. They proposed using the representative mean pore radius, {lambda}{psi}, instead of C0, to estimate the number of conductive pores per unit area, N{psi}, required to produce the measured K{psi} (Reynolds et al., 1995) as

Formula 9[9]
Although N{psi} gives a general description of the soil pore concentration with the flow-weighted mean size {lambda}{psi}, (Eq. [9]), this index does not allow us to distinguish the relative contributions of specific pore sizes to the soil water flow. Accordingly, on the basis of the models of Watson and Luxmoore (1986) (Eq. [7]) and Reynolds et al. (1995) (Eq. [9]), we propose the use of {lambda}{Delta}{psi} (Eq. [6]) to compute the number of effective water-conductive pores per unit of area for two consecutive soil water tensions, N{Delta}{psi}, as

Formula 10[10]
The effective water-conductive porosity for two consecutive soil water tensions, {theta}{Delta}{psi}, can be expressed as

Formula 11[11]
Note that {theta}{Delta}{psi}, which is based on the {lambda}{Delta}{psi} parameter, is more likely to be dominated by flow impedances (pore tortuosity, roughness, air entrapment, and obstructions) than the pore size or volume effects. Since the physical interpretation of {theta}{Delta}{psi} values to describe the total effective macro- and mesoporosity is problematic and perhaps even misleading, the use of {theta}{Delta}{psi} should be taken with caution.

The contribution of both macropores and mesopores to the total saturated water flux, K{Delta}i (%), can be calculated (Cameira et al., 2003) according to

Formula 12[12]
where n is the number of measurements performed in a sequence, Ki and Ki–1 the hydraulic conductivity for two consecutive tensions, and K0 the saturated hydraulic conductivity.

Experimental Site
The site is located at the dryland research farm of the Estación Experimental de Aula Dei (CSIC) in the province of Zaragoza (41°44'N, 0°46'W, altitude 270 m). The climate is semiarid with an average annual precipitation of 390 mm and an average annual air temperature of 14.5°C. Soil at the research site is a loam (fine-loamy, mixed, thermic Xerollic Calciorthid) according to the USDA soil classification (Soil Survey Staff, 1975). Selected physical and chemical properties of the soil were given in López et al. (1996).

The study was conducted on a large block of plots, which were set up on a nearly level area (slope 0–2%) of land in 1991 within a long-term conservation tillage experiment. The field was in a winter barley (Hordeum vulgare L.)–fallow rotation. The study was conducted in the 1999 to 2000 fallow season when the field was in the 16- to 18-mo-long fallow phase of this rotation, which extends from harvest (June–July) to sowing (November–December) the following year.

Three different fallow management treatments were compared: CT, RT, and NT. The CT treatment consisted of moldboard plowing of fallow plots to a depth of 30 to 40 cm in late winter or early spring, followed by secondary tillage with a sweep cultivator to a depth of 10 to 15 cm in late spring. In the RT treatment, the primary tillage was chisel plowing to a depth of 25 to 30 cm (noninverting action), followed as in CT by a pass of the sweep cultivator in late spring. The NT treatment used exclusively herbicides (glyphosate [N-(phosphonomethyl)glycine]) for weed control throughout the fallow season.

The tillage treatments were arranged in an incomplete block design based on geostatistical concepts, with three replications for the RT and NT treatments and four for the CT treatment (López and Arrúe, 1995). In this design, each pair of treatment plots (i.e., CT–RT, RT–NT and CT–NT) forms an incomplete block in three locations. The adverse effects of soil spatial variability are reduced by making short-distance treatment comparisons, and, by keeping this distance constant, it is ensured that all contrasts are made with equal precision. More details about this design and its efficiency can be found in López and Arrúe (1995). In this study, a 7- by 7-m region was delimited within each incomplete block for field measurements at two sampling points, located in front of each other in adjacent plots (i.e., one sampling point per treatment) and separated by a distance of 5 m. In accordance with this sampling scheme, a total of 18 measurements (six per treatment) was made. The size of the basic plot was 33.5 by 10 m, with a separation of 1 m between plots.

To contrast the effects of the tillage treatments on the soil hydrophysical properties, ANOVA for the incomplete block design was used (López and Arrúe, 1995). Duncan's multiple range test was used to compare treatment means. On the other hand, a paired t-test was performed to compare the C0 vs. {lambda}{Delta}{psi}, NC0 vs. N{Delta}{psi}, and {theta}C0 vs. {theta}{Delta}{psi} pore indexes for each tillage treatment.

Experimental Measurements
To analyze the soil bulk hydrophysical properties, excluding the soil surface crust effect, field measurements of soil bulk density and hydraulic conductivity were made on the 2- to 10-cm depth soil layer before the primary tillage was implemented in March.

The soil dry bulk density ({rho}b) was determined by the core method using commercial stainless steel cylinders (50-mm diameter by 50-mm height). Core samples were taken at the 2-cm depth near the locations for the hydraulic property measurements. This sampling was made on the same day as infiltration measurements to determine the antecedent soil bulk density and moisture. The volumetric soil water content ({theta}) was determined from the values for the gravimetric soil moisture (oven drying the soil at 105°C) and the bulk density.

The soil hydraulic properties were characterized at each sampling point using a modified Perroux and White (1988) tension disk infiltrometer with a base radius of 125 mm (Moret et al., 2004). To consider only the effects of tillage on soil water infiltration, the surface crust, clods, and crop residues were removed from the soil surface. To ensure good hydraulic contact between the disk and the soil, a thin layer of commercial sand (80–160-µm grain size) was also poured onto the soil surface. The initial volumetric water content of the contact sand was <0.05 m3 m–3. The base of the disk was covered with a nylon cloth of 20-µm mesh. Infiltration runs were performed at four values of soil tension, {psi} (namely, –14, –4, –1, and 0 cm, applied in this order and at the same place). Flow monitoring continued until steady-state flow from the disk was attained. The flow readings were automatically recorded every 30 s from the drop in the water level of the water supply reservoir of the infiltrometer, using the time domain reflectometry (TDR) application and procedure developed by Moret et al. (2004). To monitor {theta} during the water infiltration, a three-wire TDR probe (diameter 2 mm; length 100 mm; spacing between the outermost rods 25 mm) was horizontally inserted into the soil at 4-cm depth beneath the infiltrometer disk before the infiltration test. The TDR measurements for both cumulative infiltration and {theta} were synchronously and automatically recorded using the software WinTDR 98 (Or et al., 1998). To this end, the two TDR probes were connected to a multiplexer (Model SDMX50, Campbell Scientific, Logan, UT) (Moret et al., 2004). Steady-state values of {theta} beneath the infiltrometer disk at –14 cm of pressure head were not consistent enough. This probably was because the wetting front under the disk infiltrometer did not completely surround the area of influence of the TDR probe when steady state was attained.

The actual water pressure head on the soil surface, {psi}s, was recalculated according to (Reynolds and Zebchuk, 1996):

Formula 13[13]
where Q({psi}0) (m3 s–1) is the steady flow rate out of the tension infiltrometer at a set pressure head, {psi}0 (m), r (m) is the radius of the contact sand layer, Tcs (m) is the thickness of the contact sand layer, and Kcs (m s–1) is the saturated hydraulic conductivity of the contact sand layer. In our study, Tcs and Kcs were 0.0015 m and 1.05 10–4 m s–1, respectively. The soil hydraulic conductivity, K, at the different soil water pressure heads, {psi}s (i.e., K14, K4, K1, and K0) and the matric potential, {phi}, were thus calculated from the cumulative water infiltration using the multiple-head method (Ankeny et al., 1991). The contribution of both macro- and mesopores to the total saturated water flux, K{Delta}i (%), was calculated from K14, K4, and K0 according to Eq. [12].


    RESULTS AND DISCUSSION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Soil Bulk Density
Table 1 shows the soil bulk density ({rho}b) and corresponding {theta} values measured in the 2- to 7-cm-depth soil layer for the different tillage treatments. The {rho}b values after 9 yr of continuous NT were greater than those observed for the CT and RT treatments. Greater soil compaction under NT has also been observed in other long-term experiments (Logsdon et al., 1990; Evett et al., 1999; Lampurlanés and Cantero-Martínez, 2003). This is commonly associated with the gradual consolidation of the soil matrix with time owing to rainfall and the absence of annual tillage-induced loosening. The lower values of {rho}b under RT than under CT can be attributed to the greater persistence of the soil loosening after chisel plowing in comparison with moldboard plowing (López et al., 1996).


View this table:
[in this window]
[in a new window]

 
Table 1. Average dry bulk density ({rho}b) and volumetric water content ({theta}) of the surface soil (2–7 cm) for conventional tillage (CT), reduced tillage (RT), and no-tillage (NT).

 
Soil Hydraulic Conductivity
The average increase of the soil surface pressure head, {psi}s {psi}0, due to the effect of the thickness of the contact sand layer on {psi}0 was <1 mm for all the soil tensions. Figure 1 summarizes the values of the soil hydraulic conductivity (K) measured for the different {psi} and tillage treatments. As suggested by Sauer et al. (1990) and Reynolds et al. (1995), the large increase in K over the small increase in {psi} indicates that extensive macropore networks exist in the three tillage treatments near saturation. For structured consolidated soil conditions, the NT soil presented K values significantly lower than those observed in CT and RT soils for the entire range of applied soil water pressure heads. No differences in K were found, however, between the CT and RT treatments. These results are similar to those found in other studies (Sauer et al., 1990; Moreno et al., 1997; Miller et al., 1998; Evett et al., 1999), but deviate from the general long-term tendency for NT soil to increase its near-saturated K compared with CT practices (Green et al., 2003).


Figure 1
View larger version (17K):
[in this window]
[in a new window]

 
Fig. 1. Soil hydraulic conductivity (K) vs. pressure head ({Psi}) relationships for conventional tillage (CT), reduced tillage (RT), and no-tillage (NT). Bars represent the shortest significant difference (P < 0.05) for comparison among tillage treatments where significant differences were found (Duncan's test).

 
Analysis of the Pore Indexes
Representative Mean Pore Radius for Two Consecutive Soil Water Tensions
Figure 2 shows a comparison between the {lambda}{Delta}{psi} pore index (Eq. [6]) and the minimum equivalent pore radius (C0) used by Watson and Luxmoore (1986) for characterizing the water conductivity of the soil macro- and mesopores. In theory, {lambda}{Delta}{psi} for a given tension range should be found between the maximum and minimum C0 values for the corresponding range of {psi}. The fact that {lambda}{Delta}{psi} was larger than the minimum C0 for the macro- and mesopore tension range and smaller than the maximum C0 for the {psi} range between –4 and –14 cm (0.375 mm; Fig. 2) suggests the consistency of this {lambda}{Delta}{psi} pore index. Moreover, the fact that {lambda}{Delta}{psi} relates to pores participating in soil water transmission, and not to water-storage pores as in the case of the Watson and Luxmoore (1986) index, indicates that water transmission porosity would be better characterized when {lambda}{Delta}{psi} is used.


Figure 2
View larger version (19K):
[in this window]
[in a new window]

 
Fig. 2. Minimum equivalent pore radius, C0 (Eq. [3], dark points), and representative mean pore radius for two consecutive soil water tensions, {lambda}{Delta}{Psi}, (Eq. [6], white points-) estimated under conventional tillage (CT), reduced tillage (RT), and no-tillage (NT). Macropores are defined as those pores that drain at {Psi} > –4 cm (C0 > 0.375 mm) and mesopores as those pores draining at {Psi} between –4 and –14 cm (0.375 > C0 > 0.107 mm). Asterisks indicate significant difference between C0 and {lambda}{Delta}{Psi} (P < 0.05, t-test).

 
On average, {lambda}{Delta}{psi} for macropores was 86% greater than that measured for mesopores. Unlike the Watson and Luxmoore (1986) model, in which they assumed that the C0 for macro- and mesopores did not vary with measurement site, {lambda}{Delta}{psi} under NT was significantly greater than that calculated for CT and RT (Fig. 2). The greatest {lambda}{Delta}{psi} values under NT may be attributed to the increased fauna activity and the litter residues formed by accumulated organic matter (Logsdon and Kaspar, 1995) under this treatment. No differences in {lambda}{Delta}{psi} were found between CT and RT for either macro- or mesopores (Fig. 2).

Number of Effective Water-Conductive Pores per Unit Area for Two Consecutive Soil Water Tensions
The number of water-transmitting macro- and mesopores per unit area calculated according to the Watson and Luxmoore (1986) model, NC0 (Eq. [7]), was significantly greater than that estimated using the N{Delta}{psi} index (Eq. [10]; Fig. 3 ). This is because the pore radius, C0, used to calculate NC0 is significantly smaller than the {lambda}{Delta}{psi} values used for estimating N{Delta}{psi} (Fig. 2). We may thus conclude that, while NC0 estimates the theoretical maximum number of macro- and mesopores per unit area, N{Delta}{psi} characterizes the actual concentration of the soil's water-conductive macro- and mesoporosity. Differences between NC0 and N{Delta}{psi} are more evident under NT, where greater differences between {lambda}{Delta}{psi} and C0 are found (Fig. 2). The number of water-conductive mesopores per unit area, N{Delta}{psi}, was 99% greater than that calculated for macropores (Fig. 3). As regards the different tillage treatments, significantly lower values of N{Delta}{psi} for mesopores were observed under NT than under the CT and RT systems (Fig. 3).


Figure 3
View larger version (11K):
[in this window]
[in a new window]

 
Fig. 3. Maximum number of effective pores per unit area, NC0, (Eq. [7], black points) and number of effective water-transmitting pores per unit area, N{Delta}{Psi}, (Eq. [10], white points) for two consecutive water soil tensions measured under conventional tillage (CT), reduced tillage (RT), and no-tillage (NT). Macropores are defined as those pores that drain at {Psi} > –4 cm and mesopores as those pores draining at {Psi} between –4 and –14 cm. Asterisks indicate significant difference between NC0 and N{Delta}{Psi} (P < 0.05, t-test).

 
Effective Water-Conductive Porosity for Two Consecutive Soil Water Tensions
Figure 4 shows a comparison between the effective porosity estimated according to the Watson and Luxmoore (1986) model, {theta}C0 (Eq. [8]), and the volume of soil pores with the flow-weighted mean size for two consecutive soil water tensions, {theta}{Delta}{psi} (Eq. [11]). The results show that the Watson and Luxmoore (1986) model tends to overestimate the effective porosity for both water-conductive macro- and mesoporosity (Fig. 3). Similar results were obtained by Bodhinayake et al. (2004), who calculated the total conducting macro- and mesoporosity by using the analytical solutions derived from a general equation for water-conducting porosity based on tension infiltrometer measurements. In our study, the total macroporosity, calculated as the difference between the volumetric water content measured under the disk infiltrometer at 0 and –4 cm of pressure head (average value of 6%), was significantly greater than {theta}{Delta}{psi} (Fig. 4). This fact indicates that only 0.06% of total macroporosity contributes to fast water flow in soil.


Figure 4
View larger version (12K):
[in this window]
[in a new window]

 
Fig. 4. Effective porosity,{theta}C0, calculated using the minimum equivalent pore radius (C0) for a given tension range (Eq. [8], dark points) and effective porosity, {theta}{Delta}{Psi}, computed from the representative mean pore radius for two consecutive soil water tensions, {lambda}{Delta}{Psi}, (Eq. [11], white points) measured under conventional tillage (CT), reduced tillage (RT), and no-tillage (NT). Macropores are defined as those pores that drain at {Psi} > –4 cm and mesopores as those pores draining at {Psi} between –4 and –14 cm. Asterisks indicate significant difference between {theta}C0 and {theta}{Delta}{Psi} (P < 0.05, t-test).

 
The {theta}{Delta}{psi} index was also affected by the tillage treatments. Compared with CT and RT, the NT treatment shows significantly lower values of {theta}{Delta}{psi} for both macro- and mesopores (Fig. 4). By contrast, whereas no differences in {theta}{Delta}{psi} for macropores were detected between CT and RT, {theta}{Delta}{psi} for mesopores was greater under CT than under RT. This can be related to differences in soil loosening between chisel and moldboard plowing. In agreement with Miller et al. (1998), it can be concluded that the lower K values under NT are associated with the smaller number of large pores (Fig. 3), even though NT has greater values of {lambda}{Delta}{psi} (Fig. 2).

Analysis of Representative Mean Pore Radius and Number of Effective Water-Conducting Pores
The {lambda}{psi} and N{psi} vs. K relationships for structured soil conditions (Fig. 5 ) were similar to those reported by Reynolds et al. (1995). The effective equivalent mean pore radius, {lambda}{psi}, was relatively constant at its minimum value of about 0.05 mm for low K values, but then increased as K rose to higher levels. This behavior could be the result of constrictions (e.g., pore necks, entrapped air bubbles) and discontinuous pores within the {psi} range where {lambda}{psi} is constant (Reynolds et al., 1995). On the other hand, N{psi}, which is inversely related to {lambda}{psi}, increases when K decreases (Fig. 5). For all three tillage treatments, and throughout the entire range of {psi}, {lambda}{psi} varied between about 0.05 and 0.25 mm, as reported by other researchers (Sauer et al., 1990; Reynolds et al., 1995; Angulo-Jaramillo et al., 1997). Likewise, N{psi} values for the three tillage systems varied between 1000 and 17000 pores m–2 for {psi} = 0 and –14 cm, respectively. Both {lambda}{psi} and N{psi} were affected by the different tillage treatments. For K14 and K4, {lambda}{psi} was smaller and N{psi} was greater under NT than under CT and RT. These results are consistent with the higher soil bulk density found under NT (Table 1) as a consequence of the gradual consolidation of the soil matrix with time due to rainfall and the absence of annual tillage-induced loosening (Reynolds et al., 1995).


Figure 5
View larger version (21K):
[in this window]
[in a new window]

 
Fig. 5. Number of effective water-transmitting pores per unit area (N{Psi}) and representative mean pore radius ({lambda}{Psi}) vs. soil hydraulic conductivity (K) under conventional tillage (CT), reduced tillage (RT), and no-tillage (NT).

 
Compared with {lambda}{Delta}{psi} (Fig. 2), {lambda}{psi} (Fig. 5) shows a smaller pore radius for the corresponding macro- and mesopore soil water tension. These differences arise because {lambda}{psi} describes the representative average pore size for the total soil tension range, including the soil microporosity, whereas {lambda}{Delta}{psi} specifies the water-conducting pore size for a specific fraction of the soil water tension head. This would explain why the number of pores per unit area estimated using the Reynolds et al. (1995) model, N{psi}, (Eq. [9]), is higher than that calculated by N{Delta}{psi} (Eq. [10]). We may thus conclude that the soil water-conducting macro- and mesopores are better characterized using the {lambda}{Delta}{psi} and N{Delta}{psi} instead of the {lambda}{psi} and N{psi} indexes.

Contribution of Macropores and Mesopores to the Total Saturated Water Flux
On average and regardless of the tillage system, 75% of the total saturated water flux (Eq. [12]) was controlled by the macropores, and the remaining 16% by the mesopores. As Messing and Jarvis (1993) and Cameira et al. (2003) observed, these results indicate that the soil macropores have a larger influence on water flow than mesopores, even though they are a much smaller fraction of total soil water-conducting porosity. The contribution of macropores to the total saturated water flux under NT (79%) was significantly greater than that observed for the CT (69%) and RT (77%) treatments. In contrast, the contribution of mesopores to the total saturated water flux under CT (21%) was significantly higher than that estimated for RT (15%) and NT (13%).


    CONCLUSIONS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
This study evaluated an alternative pore index, {lambda}{Delta}{psi}, to characterize the water-conducting soil macro- and mesopores. The index has been tested under three different soil tillage management conditions in a long-term tillage experiment in central Aragon. The results indicate that {lambda}{Delta}{psi} is physically consistent. Unlike the Watson and Luxmoore (1986) model, which calculates the maximum concentration of water-conducting macro- and mesopores (C0) by assuming no differences in pore sizes between tillage treatments, significant differences in {lambda}{Delta}{psi} were found among the three tillage systems. The concentration of water-transmitting macro- and mesopores estimated by the capillary rise approach (Watson and Luxmoore, 1986) was too high if compared with the number of pores per unit area calculated using the {lambda}{Delta}{psi} pore index. The soil water-conducting macro- and mesoporosity are more properly characterized using {lambda}{Delta}{psi} and N{Delta}{psi} instead of {lambda}{psi} and N{psi}.

Although the largest values of {lambda}{Delta}{psi} for macropores were observed under NT, the soil hydraulic conductivity near saturation under NT was significantly lower than under CT and RT, due to the lower density of water-transmitting macro- and mesopores (N{Delta}{psi}). Regardless of the tillage system, the soil water flow at the soil surface was mainly regulated by macropores, even though this pore size represents a very small fraction of total soil porosity.


    ACKNOWLEDGMENTS
 
This research was supported by the Comisión Interministerial de Ciencia y Tecnología of Spain (Grants no. AGF98-0261-C02-02 and AGL 2001-2238-CO2-01 and PNFPI predoctoral fellowship and research contract awarded to D. Moret) and the European Union (FEDER funds). We acknowledge M.V. López, R. Gracia, and M.J. Salvador for their help in various aspects of this study.


    NOTES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
All rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher. Permission for printing and for reprinting the material contained herein has been obtained by the publisher.

Received for publication March 17, 2006.


    REFERENCES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 




This article has been cited by other articles:


Home page
Soil Sci.Home page
F. Zvomuya, F. J. Larney, P. R. DeMaere, and A. F. Olson
Hydraulic Properties of a Sandy Loam Soil following Spent Drilling Mud Application on Native Prairie
Soil Sci. Soc. Am. J., May 13, 2009; 73(4): 1108 - 1112.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Moret, D.
Right arrow Articles by Arrúe, J. L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Moret, D.
Right arrow Articles by Arrúe, J. L.
Agricola
Right arrow Articles by Moret, D.
Right arrow Articles by Arrúe, J. L.
Related Collections
Right arrow Hydraulic Conductivity
Right arrow Soil Methods/Instrumentation
Right arrow Structure and Properties


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Journal of Natural Resources
and Life Sciences Education
Vadose Zone Journal
Journal of Plant Registrations Journal of
Environmental Quality
The Plant Genome