SSSAJ Grow Your Career with SSSA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (5)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Naasz, R.
Right arrow Articles by Charpentier, S.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Naasz, R.
Right arrow Articles by Charpentier, S.
Agricola
Right arrow Articles by Naasz, R.
Right arrow Articles by Charpentier, S.
Related Collections
Right arrow Soil Physics
Right arrow Hysteresis
Right arrow Laboratory Column Studies
Published in Soil Sci. Soc. Am. J. 69:13-22 (2005).
© 2005 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA

Division S-1—Soil Physics

Measuring Hysteretic Hydraulic Properties of Peat and Pine Bark using a Transient Method

R. Naasz*, J.-C. Michel and S. Charpentier

Unité Mixte de Recherche A_462 ‘SAGAH’ INRA/INH/Université d‘Angers, Sciences Agronomiques Appliquées à l‘Horticulture, 2 rue Le Nôtre, 49045 Angers Cedex, France

* Corresponding author (remi.naasz{at}inh.fr).


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The precise and continuous measurement of the hydraulic properties of growing media is of vital importance for the effective management of irrigation and fertilization. The main purpose of this study was to use a transient method to characterize the hydraulic properties of two horticultural substrates (peat and composted pine bark) and, second, to test the applicability of well-known hydraulic models commonly used on mineral soils. Substrate water retention and hydraulic conductivity curves were determined in the laboratory by the Instantaneous Profile Method (IPM), during a drying–wetting cycle. The results showed hysteresis for peat in the {theta}({psi}) curves (30% vol.) and in the K({psi}) curves [one order of magnitude in K({psi})], whereas this phenomenon was very limited for pine bark [approximately 10% vol. in {theta}({psi}) and approximately a half-order of magnitude in K({psi})]. The K({theta}) curves for both substrates considerably decreased after the considerable variation in water content observed in {theta}({psi}), that is to say, after –5 and –3 kPa for peat and pine bark, respectively. Fitting the van Genuchten retention model to observed data resulted in a high correlation (0.96 ≤ R2 ≤ 0.99) and fitting unsaturated hydraulic conductivity models (VGM and BC) resulted in a lower correlation (0.43 ≤ R2 ≤ 0.64). Our results seemed to be in agreement with other hydraulic studies of substrates. These overall results implied that hydraulic properties of tested substrates widely fluctuate in a narrow range of water potentials that could therefore rapidly affect water and air availability to the roots. Knowledge of K({psi}) and K({theta}) curves, in addition to the {theta}({psi}) curve, could contribute to alleviating stress conditions.

Abbreviations: AC, air capacity • BC, Brooks and Corey • IPM, instantaneous profile method • PVDF, polyvinylidenfluoride • TDR, time domain reflectometry • VG, van Genuchten • VGM, van Genuchten-Mualem • WHC, water holding capacity


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
KNOWLEDGE OF WATER dynamics is essential for a better understanding of how the soil–plant system functions, particularly in terms of fertilization–irrigation management and of pollutant leaching as well. This is also and especially true in soilless cultures such as container media, due to the specificity of this production system with its extensive use of chemicals. Indeed, plants growing in containers have a limited volume of substrate in which water, gas, and solute availability highly fluctuate over a short period of time (a few hours) (Polak and Wallach, 2001), involving frequent cycles of watering (fertigation) and drying during growing management. These variations in water content could lead to a reorganization of the solid phase due to shrinkage and swelling (Heiskanen, 1995) and to a change in wettability (Valat et al., 1991; Michel et al., 2001). As a result, most organic substrates exhibit pronounced hysteresis phenomena (da Silva et al., 1993; Otten et al., 1999; Heinen and Raats, 1999), which can greatly influence the properties of water and air circulation.

However, until now, most studies dealing with the physical properties of substrates have only attempted to characterize water and gas distribution related to water potential (Bunt, 1961; De Boodt and Verdonck, 1972; Rivière et al., 1990), disregarding their availability to the plants, estimated by hydraulic conductivity and gas diffusivity measurements. Assessment of these transient properties are difficult and tedious due to the wide variety of natural and artificial, organic and mineral substrates which could present considerably different structures and physical properties as well. Furthermore, compared with most soils, horticultural substrates are coarse porous materials affected by handling (Paquet et al., 1993) and root activity (Allaire-Leung et al., 1999). However, efforts have recently been made to assess these transient properties (Laurén and Heiskanen, 1997; Caron et al., 1998) and to apply standard soil mathematical functions (van Genuchten, 1980; Mualem, 1976) to horticultural substrates (Milks et al., 1989; Wallach et al., 1992; Otten et al., 1999) with the aim of modeling water flows in these unsaturated media. But more precise and reliable characterization of the hydraulic properties of substrates is still necessary for the effective management of fertigation and for subsequent plant growth in soilless culture production.

In this context, the first objective of this study was to use a transient experimental technique (IPM) generally used on soil, estimating the hydraulic properties of two organic materials commonly used as horticultural substrates (peat and composted pine bark) during evaporation as well as during infiltration experiments in the laboratory.

The second objective of our research was to test whether or not hydraulic models (van Genuchten, 1980; Mualem, 1976; Brooks and Corey, 1964) could be used to correctly describe the substrate properties during a desiccation/infiltration cycle.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Water Dynamics Theory
Water movement in unsaturated porous media is generally described by Richards' equation (Richards, 1931):

[1]
where t is time, {psi} represents the pressure head (m), {theta} is the volumetric water content (m3 m–3) and z is depth (m, positive upward). Richards' equation can be solved numerically but knowledge of the hydraulic properties of the soil, namely the retention curve {theta}({psi}) and the hydraulic conductivity characteristic K({psi}), are required. In this study, we used the van Genuchten (1980) function to describe the water retention characteristic:

[2]
where {theta}r and {theta}s are the residual volumetric water content and the volumetric water content at saturation (m3 m–3), respectively, {alpha} is a curve fitting parameter, sometimes interpreted as the inverse of the air entry value (1/{psi}e) (m–1), and n and m are fitting constants reflecting the steepness of the retention curve. The hydraulic conductivity characteristic is described using the pore-size distribution model of Mualem (1976), combined with the water retention characteristic (Eq. [2]) hydraulic conductivity can be expressed as a function of the pressure head (Eq. [3]) or water content (Eq. [4]):

[3]

[4]
where Ks is the hydraulic conductivity at saturation and {tau} is an empirical pore connectivity factor, used as a curve fitting parameter. Parameters n and m were used independently.

We also used the model of Brooks and Corey (1964) (BC) to describe the hydraulic conductivity characteristic:

[5]
where B is an empirical fitting parameter.

Materials and Sample Preparation
Weakly decomposed sphagnum peat and composted pine bark were chosen as substrates in our study for several reasons. Indeed, sphagnum peat is internationally considered as the main reference substrate in horticultural production (in terms of quality and quantity), and composted pine bark is one of the main peat additives (particularly in France) to improve the physical properties of growing media. Furthermore, both of these materials were chosen because they present considerably different structures and physical properties as well: sphagnum peat is a fine fibrous material while pine bark is large, coarse, and platy. Consequently, peat is generally considered as a retention media whereas pine bark is considered to be an aerated material.

Sieving was performed to obtain the finest fractions possible of these substrates in this study: peat with a particle-size range of 0 to 5 mm and pine bark with a particle-size range of 0 to 10 mm. The physical characteristics of the substrates (Table 1) were calculated, taking the average value of four replicates. Bulk density, organic matter, and pH were determined according to European procedures NF EN 13041 (2000), NF EN 13039 (2000), and NF EN 13037 (2000), respectively. Particle density was estimated by the pycnometer method using petroleum spirit (Galvin, 1976) and cation exchange capacity according to the procedure of Lemaire et al. (2003). Organic elemental analysis (C and N) was performed using the procedure of Dumas and Pregl (Deléens et al., 1997).


View this table:
[in this window]
[in a new window]
 
Table 1. Main characteristics for the two substrates studied. Means of four replicate experiments with standard deviations (numbers in parentheses).

 
Since the physical properties of organic substrates are largely influenced by preparation and, more precisely, by the packing of materials, the two substrates were prepared according to the European standardized procedure NF EN 13041 (2000): two PVC (polyvinylchloride) cylinders (14 cm in diameter and 14 cm high) were manually filled but without packing with the two different media, slowly wetted (30 min) from the bottom, saturated for 24 h and then allowed to equilibrate for 48 h to a water potential of –50 cm. The cylinders were then emptied, substrates homogenized and PVDF (polyvinylidenfluoride) cylinders (9.9 cm in diameter and 12 cm high; V = 931 cm3) were filled with each substrate without packing and slowly rewetted from the bottom for 24 h.

Laboratory Procedures
Experimental Design
For each core PVDF cylinder, the water potential {psi} and volumetric water content {theta} were determined by tensiometers and time domain reflectometry (TDR) probes, respectively, installed at two levels, h1 and h2, from the bottom of the cylinder (h1 = 9 cm and h2 = 3 cm, Fig. 1) . The mini-tensiometers (2.2 mm in diameter by 20 mm in length, ceramic cell, SDEC 220, Société Développement Et Commercialisation, Reignac, France) were connected to pressure transducers (differential pressure sensors, precision: ±0.03%, response time: 10–3 s) monitored by a real-time multitasking computer to control the measurement and collect data. The TDR system consisted of TDR mini-probes (three wires, 80 mm long, with 4-mm uncoated stainless steel rods and a spacing of 10 mm). The probes were connected to a Tektronix 1502C system (Tektronix, Beaverton, OR) via a multiplexer run by WINTDR Software (Time Domain Reflectometry Soil Sample Analysis Program V. 5.1, Utah State University, Logan, UT).



View larger version (33K):
[in this window]
[in a new window]
 
Fig. 1. Schematic representation of the experimental design.

 
Time Domain Reflectometry Calibration
The calibration procedure consisted of simultaneously measuring the weight (balance) and the dielectric constant Ka (single TDR probe inserted vertically) of a PVDF cylinder filled with substrates and initially prepared as previously described (See Materials and Sample Preparation section), during drying. The volumetric water content was determined by measuring the difference of the weight during the experiment and after oven-drying at 105°C. In the present study, three replicate calibration experiments were made for each substrate.

As a result of the particularities of materials used, three different models were tested to fit the results of the measured {theta}/Ka relationship. Two of them are based on empirical polynomial functions: a third-order polynomial model (Eq. [6]) developed by Topp et al. (1980), which can be considered as the most common one used for many soils, and a second-order polynomial model (Eq. [7]) defined by Pepin et al. (1992), already used on organic substrates:

[6]

[7]
where {alpha}0, {alpha}1, {alpha}2, and {alpha}3 are fitting coefficients.

The third model (Eq. [8]) used is a partly deterministic model (the {alpha}-model) based on the theory of mixtures, which assumes that the porous media is a ternary system of water, air, and the solid fraction. The model version tested here is the one reported by Roth et al. (1990) and Todoroff and Langellier (1998):

[8]
where Ka is the dielectric constant of the medium, Kwater, Ksoil, and Kair are dielectric constants of water, the solid fraction and air, respectively, while {phi} represents the porosity. The {alpha} coefficient accounts for the effects of the geometrical arrangement of the medium components and the consequences on its dielectric performance. According to Dirksen and Dasberg (1993), {alpha} and Ksoil were considered to be fitting coefficients and Kwater and Kair were defined as 79.63 and 1, respectively, at 22°C.

Evaporation and Infiltration Experiments
After the sample preparation procedure, tensiometers and TDR were installed, and the drying and rewetting experiments could begin when readings indicated that the cores were in hydrostatic equilibrium. The experiment lasted for approximately one month and represented approximately 8000 sets of water content–water potential data.

The evaporation experiment was based on a concept similar to the evaporation method of Wind (1969), revaluated and modified by Wendroth et al. (1993) and Tamari et al. (1993). We adapted this method using an instantaneous profile procedure (Paige and Hillel, 1993; Stolte et al., 1994; Wessolek et al., 1994), which provides a detailed description of {theta} (TDR) and {psi} (tensiometers) as a function of time and space. During the evaporation experiment, the bottom of the column was sealed to prevent water loss (zero-flux bottom boundary condition) and the top of the sample was subjected to evaporation using small ventilators. The experiment was terminated when the uppermost tensiometer in the substrate core reached a suction level of approximately –20 kPa and infiltration could then begin.

The infiltration experiment was inspired by the Upward Infiltration Method (UIM) developed by Hudson et al. (1996). This study used constant flux bottom boundary conditions, which reduced the need for using the flux as an optimization parameter since the flux is independent of the soil properties (Simunek and van Genuchten, 1997). We preferred using constant pressure head bottom boundary conditions with a Mariotte system, recently tested on soil by Simunek et al. (2001) and Young et al. (2002). Three pressure levels were used during each infiltration experiment (stepwise increase): –15, –5, and 0 kPa.

Experimental Constraints
The introduction of TDR mini-probes and mini-tensiometers in the substrate column could possibly lead to a small change in the structure of materials and, as a result, in the hydraulic properties. The volume occupied by all sensors in the column was calculated and only represented <0.5% of the whole substrate core volume.

Furthermore, measurement of the water content by TDR technique during the calibration experiment was performed in a different way than during the drying–rewetting experiments. For calibration, a TDR mini-probe was installed vertically and was directly in contact with the column of substrate. A vertical water content gradient in the samples close to saturation could not be avoided but these water content variations were partially taken into account by the measurement itself: according to Ferré et al. (1996), time domain reflectometry measurements give a weighted average value of water content layers when the calibration function shows a linear relationship between {theta} and {surd}ka (which is similar to Topp's equation). The drying–rewetting experiments are based on the use of two mini-probes installed in a horizontal position and going through the entire thickness of the PVDF column, over a distance of 5.3 mm, influencing ka and water content measurements as a result. Since we knew the dielectric constant of PVDF (ka = 6.7), we approximated and verified the effective ka of the substrate column.

Hydraulic Property Measurements
Water Retention Curve
The water retention characteristic was directly measured by means of TDR and tensiometer sensors during evaporation and infiltration experiments. Data were then fitted with the van Genuchten water retention model (Eq. [2]), considering all parameters ({theta}r, {theta}s, {alpha}, n, and m) as fitting coefficients.

Unsaturated Hydraulic Conductivity
Using the IPM, the unsaturated hydraulic conductivity was calculated by direct measurement from water content and water pressure data obtained during evaporation and infiltration experiments with the Darcy equation (Green et al., 1986). The water flow through the column q is calculated from the temporal changes in water storage at two depths (h1 and h2) as follows:

[9]
and the unsaturated hydraulic conductivity is then obtained by dividing the fluxes calculated above with the hydraulic gradient (d{psi}) at the same positions (dz = h1 h2) and times:

[10]

Hydraulic conductivity models (van Genuchten-Mualem [VGM] and BC functions) and the van Genuchten retention model were simultaneously fitted to conductivity and retention data using the nonlinear least squares procedure RETC (van Genuchten et al., 1991). In the program, all hydraulic parameters ({theta}r, {theta}s, {alpha}, n, m, Ks, and {tau} or B) contained in Eq. [2] and [3] and in Eq. [2] and [5] were considered as seven independant coefficients. Their values were estimated for the best fit between retention and conductivity models to the observed data during program optimization.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Time Domain Reflectometry Calibration
Substrate moisture is plotted against the dielectric constant for peat and pine bark in Fig. 2 . Fitting the peat data simultaneously to Eq. [6] to [8] showed that the water content–dielectric constant relationship could also be represented by empirical equations (third-order polynomial regression of Topp et al. (1980) and second-order polynomial equation of Pepin et al., 1992) as a partly deterministic function ({alpha}-model from Todoroff and Langellier, 1998), where all three R2 values > 0.99 (Table 2). On the contrary, the Topp equation gave better results for pine bark (R2 = 0.94) than the Pepin empirical regression (R2 = 0.901) and {alpha} model (R2 = 0.89). Therefore, we used Topp's polynomial regression to estimate water content values for the two substrates.



View larger version (23K):
[in this window]
[in a new window]
 
Fig. 2. A comparison between observed data (symbols) and curves calculated with the equations of Topp et al. (1980), Pepin et al. (1992), and the {alpha}-model (Roth et al., 1990). Each curve was simultaneously fitted to the data from the three replicate experiments for (a) peat and for (b) pine bark.

 

View this table:
[in this window]
[in a new window]
 
Table 2. Results of fitting Eq. [6] to [8] observed data. All parameters are fitted except Kwater and Kair.

 
Retention Curves
During desiccation, the two substrates differed significantly in their water retention properties (Fig. 3) . Peat exhibited high water content (>85%) in the 0 to –1 kPa range, a very low air capacity (1%; AC = air content at –1 kPa), and a high water holding capacity (45%; WHC = water loss between –1 and –10 kPa). On the contrary, pine bark presented a lower water content (>75%) between 0 and –1 kPa, higher AC (5%), but lower WHC (25%) than peat. More than 20% of the water volume was lost in the –1- to –3-kPa range.



View larger version (39K):
[in this window]
[in a new window]
 
Fig. 3. Measured (symbols) and fitted (lines, VG) water retention curves during drying and wetting for (a) peat and (b) pine bark. Squares, diamonds, and triangles correspond to Replicates 1 to 3, respectively.

 
The physical behavior of these two substrates during rewetting also indicated significant differences. Peat exhibited pronounced hysteresis phenomena (da Silva et al., 1993; Otten et al., 1999) in water retention for a range of water potentials varying between –1 and –10 kPa. In this range of water potential, the difference was as high as 30% volume in water retention between drying and rewetting curves. On the other hand, these phenomena were not observed near saturation (between saturation and –1 kPa) and for desiccation greater than –10 kPa. This implies a decrease of WHC (35% between –1 and –10 kPa), a higher AC at –1 kPa (10%) and similar water content at saturation (87%). Hysteresis phenomena were less pronounced in the case of pine bark than in peat: the difference in water content between drying and wetting reached only 10% maximum for the same water potential value, and the AC (10%), the WHC (22%), and the water content at saturation (76%) were little modified during wetting as compared with drying.

Comparison of fitted curves with observed data (Fig. 3) showed the flexibility of the van Genuchten model (VG) describing substrate retention properties, as previously shown by Milks et al. (1989) and Weiss et al. (1998). Regardless of the substrate, a high correlation between measured and fitted data was obtained for both the drying and wetting curves and throughout the tested potential range (Table 3). Parameters resulting from the VG are presented and compared with other substrates studies in Table 2. Our parameters are in good agreement with the steady-state measurements of Milks et al. (1989) but in poorer agreement with those of da Silva et al. (1993) or the transient results of Otten et al. (1999) with this kind of material. This could be the result of the fitting procedure, where only the parameter n was fitted in the studies of da Silva et al. (1993) and Otten et al. (1999), while in our study n and m parameters were fitted as in the study of Milks et al. (1989).


View this table:
[in this window]
[in a new window]
 
Table 3. Hydraulic model parameters of Eq. [2]: fitted values of {theta}s, {theta}r, {alpha}, n, and m.

 
Curves obtained on the finest fractions of raw material (peat with a particle-size range of 0–5 mm and pine bark with a particle-size range of 0–10 mm) were in agreement with data obtained on peat soils during transient procedures (Laurén and Heiskanen, 1997; Schlotzhauer and Price, 1999; Schwärzel et al., 2002).

Hydraulic Conductivity
Calculated and fitted K({theta}) curves (Fig. 4 and 5) for both substrates were established over a wide range of water content. The hydraulic conductivity of peat decreased by about four orders of magnitude for water contents varying from 85 to 35%. However, K({theta}) varied slowly between 85 and 50% (more than one order of magnitude) and then decreased about three orders of magnitude between 50 and 35%. For the pine bark, higher K({theta}) decreased (five orders of magnitude) for a smaller water content range varying between 80 and 50%. The decrease was approximately one order of magnitude, between 80 and 55%), and was then very rapid (four orders of magnitude) around 50 to 55%. These results showed that the hydraulic conductivity K({theta}) of both substrates considerably decreased after the considerable variation in water content previously observed (Fig. 3), that is to say, after –5 and –3 kPa for peat and pine bark, respectively.



View larger version (46K):
[in this window]
[in a new window]
 
Fig. 4. Calculated (symbols) and predicted (lines) hydraulic conductivity as a function of volumetric water content for (a) peat and (b) pine bark. The drying and wetting hydraulic conductivity curves (van Genuchten–Mualem [VGM]) were predicted (Eq. [4]) from the fitted drying and wetting water retention curves (Eq. [2]).

 


View larger version (45K):
[in this window]
[in a new window]
 
Fig. 5. Calculated (symbols) and predicted (lines) hydraulic conductivity as a function of volumetric water content for (a) peat and (b) pine bark. The drying and wetting hydraulic conductivity curves (BC) were predicted (Eq. [5]) from the fitted drying and wetting water retention curves (Eq. [2]).

 
Calculated K({psi}) was also represented (Fig. 6) over a large potential range (0 to –26 kPa). For peat, the variation of K({psi}) related to the water potential is around four orders of magnitude during drying and three orders of magnitude during wetting. Differences between drying and wetting appeared after –3 kPa in the desiccation direction (K = 0.1 cm min–1), where the drying curve decreases more rapidly than the wetting curve. As opposed to K({theta}) curves, this result suggest that hysteresis phenomena are observable during the evolution of peat hydraulic conductivity and appear before the rapid decrease in K({theta}) previously observed in Fig. 4 and 5. On the other hand, the general evolution in K({psi}) for pine bark was different from that of peat. A higher general decrease (three orders of magnitude during drying and wetting curves) was observed. The K({psi}) values for the drying curve were always lower than those of the wetting curve, and differences between the two branches seemed to appear from saturation but do not exceed a half-order of magnitude. This last result seemed to indicate poor hysteresis phenomena for pine bark as compared with peat, which is in line with retention curve observations (Fig. 3).



View larger version (50K):
[in this window]
[in a new window]
 
Fig. 6. Calculated (symbols) and predicted (lines) hydraulic conductivity as a function of water potential for (a) peat and (b) pine bark. The drying and wetting hydraulic conductivity curves (VGM) were predicted (Eq. [3]) from the fitted drying and wetting water retention curves (Eq. [2]).

 
Comparison between calculated and fitted K({theta}) and K({psi}) curves for both substrates could be observed from Fig. 4 to 6, knowing that the fitted curves (VGM and BC) were obtained from all parameters ({theta}r, {theta}s, {alpha}, n, and m) estimated in Eq. [2] and by simultaneously fitting Ks and {tau} in Eq. [4] and Ks and B in Eq. [5]. All parameters can be found in Table 4, with the correlation coefficient R2. For the two substrates, a relatively low correlation between calculated and fitted hydraulic conductivity data was obtained (R2 varying between 0.43 and 0.64). The correlation is always better (except for the peat drying curve) when fitting data with the BC model in comparison with the VGM model. The m parameter in the VGM hydraulic conductivity model did not seem to provide additional accuracy in the fitting procedure. A comparison of our results with the study of Otten et al. (1999) (Table 4) and the results of Wallach et al. (1992), da Silva et al. (1993), and Caron et al. (1998) revealed similar hydraulic conductivity curves but with a greater decrease in their hydraulic conductivity curves. This was probably due to our fine materials with less hydraulic conductivity under wet conditions but decreasing slowly and, consequently, with more hydraulic conductivity during drying.


View this table:
[in this window]
[in a new window]
 
Table 4. Hydraulic model parameters of Eq. [4] and [5]. All are fitted parameters.

 
The overall results concerning hydraulic conductivity [K({theta}) and K({psi})] curves suggest different behaviors during a drying–wetting cycle. First, there was a sharp decrease in the K curve arising after the considerable water variation in the retention curve (Fig. 3). If we start from the concept of available water (AW = defined as the difference between the water content of the substrate at –1 and –10 kPa) introduced by De Boodt and Verdonck (1972), our results (like those of da Silva et al., 1993), suggested that the sharp decrease in hydraulic conductivity could occur well before this empirical threshold of –10 kPa (–5 and –3 kPa for peat and pine bark, respectively). Moreover, very small changes in suction appeared to cause K to decrease considerably and, consequently, affect the availability of water to the root environment. According to Laurén and Heiskanen (1997), in samples with larger pores (pine bark), the hydraulic conductivity is considerable under wet conditions but deceases more rapidly during the drying process. On the other hand, material with smaller pores (peat) has lower hydraulic conductivity when wet but the K decreases at lower rates during drying.

The appearance of hysteretic phenomena could be related to the different structure of these two materials, fine and fibrous in the case of peat, and coarse and platy for pine bark. With regard to peat, these hysteresis phenomena were probably due to the change in the solid phase organization (swelling/shrinkage phenomena) and consequently in the pore interconnection, as well as to variations in wettability as was suggested by Michel et al. (2001). As was shown by Valat et al. (1991) and Michel et al. (2001), peat can acquire a hydrophobic character during extreme desiccation and could cause difficulties rewetting the substrate.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Substrate water retention and hydraulic conductivity curves were determined in the laboratory by a transient procedure (IPM), during a drying–wetting cycle. The relative simplicity and rapidity of this experimental procedure and its ability to directly estimate hydraulic characteristics make the method very appealing and seemed to be well adapted for determining these substrate characteristics. Nevertheless, methodological improvements should be implemented, particularly with regard to the shrinkage–swelling phenomena. Organic materials as peat showed substantial volume change (until 20%, Heiskanen, 1995) leading to a reorganization of the solid phase and could certainly affect hydraulic properties of substrate.

Results showed differences in the physical behavior of the two substrates studied. Hysteresis phenomena were evident in the {theta}({psi}) and the K({psi}) curves for peat, whereas this phenomenon was very limited for pine bark. Furthermore, the K({theta}) curve suggested that the hydraulic conductivity of both substrates considerably decreased after the high variation in water content observed in Fig. 3, that is to say, after –5 and –3 kPa for peat and pine bark, respectively.

The use of the VG (van Genuchten, 1980) retention model to describe the water retention characteristics of our two materials revealed a high correlation and seemed to be in agreement with other hydraulic studies of substrates and, more particularly, with results obtained on peat soils, due to the specific fine raw materials used. The use of unsaturated hydraulic conductivity models (VGM and BC) showed less agreement but they were also in agreement with results previously obtained on substrates.

These overall results implied that the hydraulic properties of tested substrates widely fluctuate within a narrow range of water potentials. Threshold values of water availability, obtained for the two substrates, seemed not to correspond to the common empirical value of –10 kPa (De Boodt and Verdonck, 1972) and could provide inaccurate predictions as to potential plant response. In practice, inaccurate predictions could rapidly affect water and air availability and lead to stress conditions in the root environment. Knowledge of K({psi}) and K({theta}) curves, in addition to the {theta}({psi}) curve, could contribute to alleviating water stress conditions and improve quality of growing media. Moreover, hydraulic properties (unsaturated hydraulic conductivity and water retention) should be associated with an aeration index of plant growth (e.g., gas relative diffusivity) and both combined in water–gas flow model to better understand the substrate-plant system and improve fertigation management.

Received for publication October 17, 2003.


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




This article has been cited by other articles:


Home page
Vadose Zone JHome page
K. Schwarzel, J. Simunek, H. Stoffregen, G. Wessolek, and M. Th. van Genuchten
Estimation of the Unsaturated Hydraulic Conductivity of Peat Soils: Laboratory versus Field Data
Vadose Zone J., May 26, 2006; 5(2): 628 - 640.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (5)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Naasz, R.
Right arrow Articles by Charpentier, S.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Naasz, R.
Right arrow Articles by Charpentier, S.
Agricola
Right arrow Articles by Naasz, R.
Right arrow Articles by Charpentier, S.
Related Collections
Right arrow Soil Physics
Right arrow Hysteresis
Right arrow Laboratory Column Studies


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