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-1Soil 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é dAngers, Sciences Agronomiques Appliquées à lHorticulture, 2 rue Le Nôtre, 49045 Angers Cedex, France
* Corresponding author (remi.naasz{at}inh.fr).
 |
ABSTRACT
|
|---|
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 dryingwetting cycle. The results showed hysteresis for peat in the
(
) curves (30% vol.) and in the K(
) curves [one order of magnitude in K(
)], whereas this phenomenon was very limited for pine bark [approximately 10% vol. in
(
) and approximately a half-order of magnitude in K(
)]. The K(
) curves for both substrates considerably decreased after the considerable variation in water content observed in
(
), 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(
) and K(
) curves, in addition to the
(
) 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
|
|---|
KNOWLEDGE OF WATER dynamics is essential for a better understanding of how the soilplant system functions, particularly in terms of fertilizationirrigation 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
|
|---|
Water Dynamics Theory
Water movement in unsaturated porous media is generally described by Richards' equation (Richards, 1931):
 | [1] |
where t is time,
represents the pressure head (m),
is the volumetric water content (m3 m3) 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
(
) and the hydraulic conductivity characteristic K(
), are required. In this study, we used the van Genuchten (1980) function to describe the water retention characteristic:
 | [2] |
where
r and
s are the residual volumetric water content and the volumetric water content at saturation (m3 m3), respectively,
is a curve fitting parameter, sometimes interpreted as the inverse of the air entry value (1/
e) (m1), 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
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
and volumetric water content
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: 103 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).
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
/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
0,
1,
2, and
3 are fitting coefficients.
The third model (Eq. [8]) used is a partly deterministic model (the
-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
represents the porosity. The
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),
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 contentwater 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
(TDR) and
(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 dryingrewetting 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
and
ka (which is similar to Topp's equation). The dryingrewetting 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 (
r,
s,
, 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
) 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 (
r,
s,
, n, m, Ks, and
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
|
|---|
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 contentdielectric 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 (
-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
model (R2 = 0.89). Therefore, we used Topp's polynomial regression to estimate water content values for the two substrates.
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).
Curves obtained on the finest fractions of raw material (peat with a particle-size range of 05 mm and pine bark with a particle-size range of 010 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(
) 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(
) 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(
) 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(
) 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 GenuchtenMualem [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(
) was also represented (Fig. 6)
over a large potential range (0 to 26 kPa). For peat, the variation of K(
) 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 min1), where the drying curve decreases more rapidly than the wetting curve. As opposed to K(
) curves, this result suggest that hysteresis phenomena are observable during the evolution of peat hydraulic conductivity and appear before the rapid decrease in K(
) previously observed in Fig. 4 and 5. On the other hand, the general evolution in K(
) 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(
) 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(
) and K(
) 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 (
r,
s,
, n, and m) estimated in Eq. [2] and by simultaneously fitting Ks and
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.
The overall results concerning hydraulic conductivity [K(
) and K(
)] curves suggest different behaviors during a dryingwetting 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
|
|---|
Substrate water retention and hydraulic conductivity curves were determined in the laboratory by a transient procedure (IPM), during a dryingwetting 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 shrinkageswelling 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
(
) and the K(
) curves for peat, whereas this phenomenon was very limited for pine bark. Furthermore, the K(
) 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(
) and K(
) curves, in addition to the
(
) 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 watergas flow model to better understand the substrate-plant system and improve fertigation management.
Received for publication October 17, 2003.
 |
REFERENCES
|
|---|
- Allaire-Leung, S.E., J. Caron, and L.E. Parent. 1999. Changes in physical properties of peat substrates during plant growth. Can. J. Soil Sci. 79:137139.
- Brooks, R.H., and A.T. Corey. 1964. Hydraulic properties of porous media. Hydrol. Pap. No. 3. Colorado State Univ., Fort Collins.
- Bunt, A.C. 1961. Some physical properties of pot-plant composts and their effect on plant growth. I. Bulky physical conditioners. Plant Soil 13:322332.
- Caron, J., H.L. Xu, P.Y. Bernier, I. Duchesne, and P. Tardif. 1998. Water availability in three artificial substrates during Prunus x cistena growth: Variable threshold values. J. Am. Soc. Hortic. Sci. 123:931936.[Abstract/Free Full Text]
- da Silva, F.F., R. Wallach, and Y. Chen. 1993. Hydraulic properties of sphagnum peat moss and tuff (scoria) and their potential effects on water availability. Plant Soil 154:119126.
- De Boodt, M., and O. Verdonck. 1972. The physical properties of the substrates in horticulture. Acta Hort. 26:3744.
- Deléens, E., J.F. Morot-Gaudry, F. Martin, A. Thoreux, and A. Gojon. 1997. Méthodologie 15N. p. 265280. In J.F. Morot-Gaudry (ed.) Assimilation de l'azote chez les plantes supérieures. Aspects physiologique, biochimique et moléculaire. INRA, Paris.
- Dirksen, C., and S. Dasberg. 1993. Improved calibration of time domain reflectometry soil water content measurements. Soil Sci. Soc. Am. J. 57:660667.[Abstract/Free Full Text]
- Ferré, P.A., D.L. Rudolph, and R.G. Kachanoski. 1996. Spatial averaging of water content by time domain reflectometry: Implications for twin rod probes with and without dielectric coatings. Water Resour. Res. 32:271279.
- Galvin, L.F. 1976. Physical properties of Irish peats. Ir. J. Agric. Res. 15:207221.
- Green, R.E., L.R. Ahuja, and S.K. Chong. 1986. Hydraulic conductivity, diffusivity and sportivity of unsaturated soils: Fields methods. p. 771798. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Heinen, M., and P.A.C. Raats. 1999. Hydraulic properties of root zone substrates used in greenhouse horticulture. p. 467476. In M.Th. Van Genuchten et al. (ed.) Proc. Int. Workshop on the characterization and measurement of the hydraulic properties of unsaturated porous media. Univ. of California, Riverside.
- Heiskanen, J. 1995. Physical properties of two-component growth media based on sphagnum peat and their implications for plant-available water and aeration. Plant Soil 172:4554.
- Hudson, D.B., P.J. Wierenga, and R.G. Hills. 1996. Unsaturated hydraulic properties from upward flow into soil cores. Soil Sci. Soc. Am. J. 60:388396.[Abstract/Free Full Text]
- Laurén, A., and J. Heiskanen. 1997. Physical properties of the mor layer in a scots pine stand: I. Hydraulic conductivity. Can. J. Soil Sci. 77:627634.
- Lemaire, F., A. Dartigues, L.M. Riviere, S. Charpentier, and P. Morel. 2003. Cultures en pots et conteneurs. Principes agronomiques et applications. 2nd ed. INRA, Paris.
- Michel, J.C., L.M. Rivière, and M.N. Bellon-Fontaine. 2001. Measurement of the wettability of organic materials in relation to water content by capillary rise method. Eur. J. Soil Sci. 52:459467.
- Milks, R.R., W.C. Fonteno, and R.A. Larson. 1989. Hydrology of horticultural substrates: I. Mathematical models for moisture characteristics of horticultural container media. J. Am. Soc. Hortic. Sci. 114:4852.
- Mualem, Y. 1976. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 12:513522.
- NF EN 13037. 2000. Amendements du sol et supports de cultureDétermination du pH. Association Française de Normalisation, Paris, France.
- NF EN 13039. 2000. Amendements du sol et supports de cultureDétermination de la matière organique et des cendres. Association Française de Normalisation, Paris, France.
- NF EN 13041. 2000. Amendements du sol et supports de cultureDétermination des propriétés physiquesMasse volumique apparente sèche, volume d'air, volume d'eau, valeur de rétraction et porosité totale. Association Française de Normalisation, Paris, France.
- Otten, W., P.A.C. Raats, H. Challa, and P. Kabat. 1999. Saptial and temporal dynamics of water in the root environment of potted plants on a flooded bench fertigation system. Neth. J. Agric. Sci. 47:5165.
- Paige, G.B., and D. Hillel. 1993. Comparison of three methods for assessing soil hydraulic properties. Soil Sci. 155:175189.
- Paquet, J.M., J. Caron, and O. Banton. 1993. In situ determination of the water desorption characteristics of peat substrates. Can. J. Soil Sci. 73:329339.
- Pepin, S., A.P. Plamandon, and J. Stein. 1992. Peat water content measurement using time domain reflectometry. Can. J. For. Res. 22:534540.
- Polak, A., and R. Wallach. 2001. Measuring soil moisture dynamics in an irrigated orchard by time domain reflectometry method. Acta Hort. 562:3946.
- Richards, L.A. 1931. Capillary conduction of liquids through porous mediums. Physics 1:318333.
- Rivière, L.M., J.C. Foucard, and F. Lemaire. 1990. Irrigation of container crops according to the substrate. Sci. Hortic. (Amsterdam) 43:339349.
- Roth, K., R. Schulin, H. Flühler, and W. Attinger. 1990. Calibration of time domain reflectometry for water content measurement using a composite dielectric approach. Water Resour. Res. 26:22672273.
- Schlotzhauer, S.M., and J.S. Price. 1999. Soil water flow dynamics in a managed cutover peat field, Quebec: Field and laboratory investigations. Water Resour. Res. 35:36753683.
- Schwärzel, K., M. Renger, R. Sauerbrey, and G. Wessolek. 2002. Soil physical characteristics of peat soils. J. Plant Nutr. 165:479486.
- Simunek, J., and M.Th. Van Genuchten. 1997. Estimating unsaturated soil hydraulic properties from multiple tension disc infiltrometer data. Soil Sci. 162:383398.
- Simunek, J., O. Wendroth, N. Wypler, and M.Th. Van Genuchten. 2001. Non-equilibrium water flow characterized by means of upward infiltration experiments. Eur. J. Soil Sci. 52:1324.
- Stolte, J., J.I. Frijer, W. Bouten, C. Dirksen, J. Halbertsma, J.C. Van Dam, J.A. Van den Berg, G.J. Veerman, and J.H.M. Wösten. 1994. Comparison of six methods to determine unsaturated soil hydraulic conductivity. Soil Sci. Soc. Am. J. 58:15961603.[Abstract/Free Full Text]
- Tamari, S., L. Bruckler, J. Halbertsma, and J. Chadoeuf. 1993. A simple method for determining soil hydraulic properties in the laboratory. Soil Sci. Soc. Am. J. 57:642651.[Abstract/Free Full Text]
- Todoroff, P., and P. Langellier. 1998. Comparison of empirical and partly deterministic methods of time domain reflectometry calibration, based on a study of two tropical soils. Soil Tillage Res. 45:325340.
- Topp, G.C., J.L. Davis, and A.P. Annan. 1980. Electromagnetic determination of soil water content: Measurements in coaxial transmission lines. Water Resour. Res. 16:574582.
- Valat, B., C. Jouany, and L.M. Rivière. 1991. Characterization of the wetting properties of air-dried peats and composts. Soil Sci. 152:100107.
- van Genuchten, M.Th. 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44:892898.[Abstract/Free Full Text]
- van Genuchten, M.Th., F.J. Leij, and S.R. Yates. 1991. The RETC code for quantifying hydraulic functions of unsaturated soils, Version 1.0. EPA Report 600/291/065, USDA ARS, U.S. Salinity Laboratory, Riverside, CA.
- Wallach, R., F.F. da Silva, and Y. Chen. 1992. Unsaturated hydraulic characteristics of composted agricultural wastes, tuff and their mixtures. Soil Sci. 153:434441.
- Weiss, R., J. Alm, R. Laiho, and J. Laine. 1998. Modelling moisture retention in peat soils. Soil Sci. Soc. Am. J. 62:305313.[Abstract/Free Full Text]
- Wendroth, O., W. Ehlers, J.W. Hopmans, H. Kage, J. Halbertsma, and J.H.M. Wösten. 1993. Reevaluation of the evaporation method for determining hydraulic functions in unsaturated soils. Soil Sci. Soc. Am. J. 57:14361443.[Abstract/Free Full Text]
- Wessolek, G., R. Plagge, F.J. Leij, and M.Th. Van Genuchten. 1994. Analysing problems in describing field and laboratory measured soil hydraulic properties. Geoderma 64:93110.
- Wind, P. 1969. Capillary conductivity data estimated by a simple method. p. 181198. In P.E. Rijtema and H.Wassink (ed.) Water in the unsaturated zone. Proc.Wageningen Symp. June 1966. Vol. 1 IASAH, Genbrugge, Belgium.
- Young, M.H., A. Karagunduz, J. Simunek, and K.D. Pennell. 2002. A modified upward infiltration method for characterizing soil hydraulic properties. Soil Sci. Soc. Am. J. 66:5764.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
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]
|
 |
|