Soil Science Society of America Journal 64:1211-1219 (2000)
© 2000 Soil Science Society of America
DIVISION S-1-SOIL PHYSICS
Hydraulic Gradient and Wetting Rate Effects on the Hydraulic Conductivity of Two Calcium Vertisols
M. Moutiera,
I. Shainbergb and
G.J. Levyb
a Université catholique de Louvain, Unité de Génie Rural, Place Croix du Sud 2, bte 2, B-1348 Louvain-la-Neuve, Belgium
b Inst. of Soils, Water, and Environmental Sciences, Agricultural Research Organization (ARO), The Volcani Center, P.O. Box 6, Bet Dagan 50-250, Israel
marylene.moutier{at}skynet.be
 |
ABSTRACT
|
|---|
Understanding the combined role of intrinsic (e.g., clay content) and extrinsic soil conditions in determining the saturated hydraulic conductivity (K) is a key factor in improving soil and irrigation management. Hydraulic gradient effects on K have been investigated on two Ca vertisols (Chromic Haploxerert) from Yizreel and Kedma, Israel. Samples were packed in columns and subjected to two prewetting rates (4.5 and 70 mm h-1) and two hydraulic gradients (2.8 and 11.6 for Yizreel; 3.1 and 13.6 for Kedma). Saturated K was determined during leaching with CaCl2 solutions having total electrolyte concentrations (TEC) of 0.5, 0.01 M Cl-, and deionized water (DW). The average hydraulic conductivity at the end of the leaching with the 0.5 M solution
of the two vertisols increased with a decrease in prewetting rate. The effect of prewetting rate was more pronounced in Yizreel, where the high clay content (70.3% clay) resulted in a more stable structure. Upon leaching with the 0.01 M solution,
r0.01 first decreased and then increased. This increase in
r0.01 was explained by cohesive bond formation, which increased with an increase in soil clay content and with increased proximity between the clay particles. High hydraulic gradient enhanced clay to clay contacts, and a steep increase in
r0.01. Following fast prewetting, the effect of the hydraulic gradient on cohesive bond formation was more pronounced in Yizreel than in Kedma (46.5% clay). Leaching with DW decreased the
rDW, regardless of the prewetting treatment or soil clay content, suggesting that swelling was the governing mechanism in the reduction of
rDW.
Abbreviations: DW, deionized water EC, electrical conductivity SAR, Na adsorption ratio TEC, total electrolyte concentration
 |
INTRODUCTION
|
|---|
THE SATURATED HYDRAULIC CONDUCTIVITY (K) of soils depends on soil permanent properties (e.g., soil texture, clay mineralogy), soil structure, presence of entrapped air (e.g., Christiansen, 1944; De Backer, 1967; Constantz et al., 1988; Faybishenko, 1995), type of exchangeable cations, and salt concentration of the leaching solution. Saturated hydraulic conductivity tends to decrease with increasing exchangeable Na percentage and decreasing TEC (Quirk and Schofield, 1955). Recently, aggregate stability, infiltration rate, and erosion studies demonstrated that soil structure and its hydraulic properties depended on, in addition to sodicity and water quality, the antecedent water content, aging, and prewetting rate (Kemper and Rosenau, 1984; Shainberg et al., 1996; Levy et al., 1997).
The prewetting rate influences aggregate stability (Le Bissonnais, 1988, 1989; Levy et al., 1997) and, hence, K. Whereas slow prewetting minimizes degradation of the original structure of the soil, fast prewetting disintegrates the aggregates and deteriorates the hydraulic properties of soils (Kemper and Koch, 1966). Aggregate breakdown (slaking) results from the development of internal pressures that cause them to explode (Kemper and Rosenau, 1984). These pressures can arise from a differential hydration and swelling of the clay fraction (Le Bissonnais, 1989) or the compression of occluded air in the capillary pores (Panabokke and Quirk, 1957; Boiffin, 1984). Murray and Quirk (1990) suggested that aggregate slaking was due to both the weakening of the cohesive interparticle binding forces by wetting and the appearance of repulsive interparticle forces. The presence of organic matter and inorganic cementing agents such as Fe oxides will contribute to the stability of aggregates. The extent of slaking and the relative importance of each process controlling it were found to depend on TEC, Na adsorption ratio (SAR) (Abu-Sharar et al., 1987), and clay content (Le Bissonnais, 1988).
Hydraulic properties are also affected by the presence of entrapped air, the amount of which depends on the method used for the initial saturation (Faybishenko, 1995). Following an upward initial saturation, which occurs in soils during the rise of the groundwater table, the volume of entrapped air is <5% and mobile air is almost absent. Upon wetting a dry soil from the top, the volume of entrapped air is large and it blocks a significant volume of the water-conducting pores. As entrapped air dissolves, K increases. Increases of K by a factor of 5 to 10 were reported by Constantz et al. (1988), whereas Christiansen (1944) determined that the relative permeability increased by 2 to 40 times. It should be noted that when soil columns are wetted from below, the volume of entrapped air is small and most of it is immobile; changes in K upon air dissolution are therefore insignificant (Faybishenko, 1995).
An increase of aggregate mechanical strength with time has been reported by several authors (e.g., Blake and Gilman, 1970; Singer et al., 1992; Shainberg et al., 1996). Blake and Gilman (1970) observed significant thixotropic changes in artificial aggregates (Webster clay loam [fine-loamy, mixed, superactive, mesic Typic Endoaquoll] with 34.9% clay) within 20 to 30 h. Significant development of cohesive forces in aggregates of a vertisol (46.5% clay) were reported during 24 h of aging (Shainberg et al., 1996). The development of cohesive bonds has been recently inferred from K measurements in a smectitic vertisol (Moutier et al., 1998).
The hydraulic conductivity of a disturbed vertisol (A horizon, Chromic Haploxerert) wetted by capillary rise was shown to be affected by aging and by the hydraulic gradient applied to the soil column (Moutier et al., 1998). The relative hydraulic conductivity (Kr) initially decreased before increasing as a function of leaching volume or time. A combination of three mechanisms was suggested to explain the observed changes in saturated Kr (Moutier et al., 1998):
 |
Swelling
|
|---|
Swelling of the quasi-crystals (reduction in the intraaggregate porosity) was offered to explain the decrease of K (Blackmore and Miller, 1961) when a more dilute solution replaced a more concentrated one.
 |
Physical Compression of the Clay Fabric
|
|---|
Kemper et al. (1972) observed the physical compression of the clay fabric at the outflow side of a clay column. This compression of clay resulted in lower K values in the outflow side of the column, and thus, lower average K values for the whole column. A similar observation was made in our laboratory in a leaching experiment on Ca-saturated Yizreel soil samples packed in columns equipped with piezometers (I. Shainberg, 1997, unpublished data). It was noted that with leaching, the hydraulic gradient in the lower part of the column increased (i.e., a decrease of K). Because water flow in such a column is controlled by the layer with the lowest K (i.e., a bottleneck for the water flow), the average K decreased. It should be noted that the decrease in K was not accompanied by a measurable change in the soil column length. The degree of compression increased with increasing hydraulic gradient and decreasing SARs (Moutier et al., 1998), as greater repulsion forces develop with higher SAR values (Kemper et al., 1972).
 |
Development of Cohesive Bonds
|
|---|
The development of cohesive bonds between clay fabric structural units and microaggregates in close contact (Kemper and Rosenau, 1984; Bresson and Boiffin, 1990; Shainberg et al., 1996; Levy et al., 1997), unlike the former two mechanisms, leads to an increase in the average size of transmission pores, and hence an increase in Kr (Moutier et al., 1998). This concept of "reaggregation" is similar to the "coalescence" idea presented by Bresson and Boiffin (1990). Coalescence is regarded as the welding of initially loose aggregates into larger, stronger structural units by plastic deformation that is induced by the prewetting procedure. Such structural changes in seedbeds of a red-brown earth wetted by capillary rise were also observed by Bresson and Moran (1995).
The rate of cohesive bonds development, as inferred from the increase of Kr with time, was also shown to depend on the hydraulic gradient applied and water quality (Moutier et al., 1998). The higher the hydraulic gradient and the lower the SAR of the leaching solution, the steeper the increase of Kr with time of aging. This suggested that the proximity of the clay particles was essential for the formation of the cohesive forces at the aggregate level (Zubkova, 1998). The "proximity" consideration explains some conflicting results related to cohesive bond development under saturated conditions where water layers surrounding the clay particles may prevent close contacts (Blake and Gilman, 1970; Kemper and Rosenau, 1984; Shainberg et al., 1996).
The combined effect of intrinsic soil properties and extrinsic conditions in the measurement of saturated K should be studied in greater detail. We hypothesize that the hydraulic gradient effect on controlling the saturated hydraulic conductivity might be affected by the initial structural state of the soil sample. Thus, this study aims to investigate the effect of a hydrostatically induced pressure on saturated K, as measured with the constant head method, for two Ca-saturated vertisols of differing clay contents. The contrasting initial structure was achieved by varying the prewetting treatment prior to the K measurements.
 |
Materials and methods
|
|---|
Disturbed samples of two smectitic soils (Chromic Haploxerert) were taken from the A horizon (0250 mm) of two cultivated fields from the Yizreel Valley and the Pleshet Plains (Kedma), Israel. Soil samples were physically and chemically characterized by standard methods (Klute, 1986; Page et al., 1986) (Table 1) . The clay fraction of these soils contained
64% smectite, 13% kaolinite, 4% illite, 12% calcite, and 7% sesquioxides (Banin and Amiel, 1970).
Soil columns were prepared by packing 120 g of sieved (<2 mm) and air-dried soil into small cylinders (5.4 cm in diameter by
4 to 5 cm long with a metal screen covered with sand at the bottom). The dry bulk densities of Yizreel and Kedma soils were 1.33 ± 0.02 (one standard deviation) and 1.46 ± 0.04 Mg m-3, respectively. A filter paper covered the surface of the soil to minimize soil disturbance when replacing solutions. After packing, the soil columns were prewetted from below, either quickly (20 min, fast prewetting treatment) or slowly (8 h, slow prewetting treatment), with a CaCl2 solution (0.5 M Cl-). Calcium systems were studied because the effects of hydraulic gradient and aging were reported to be more pronounced in Ca-saturated soil samples (Moutier et al., 1998). The prewetting treatments were chosen to simulate likely field conditions. The fast prewetting was achieved by using a Mariotte bottle with zero water pressure head set at the bottom of the dry soil column. Once the water level reached the top of the sample, the pressure head was increased to expose the samples to an hydraulic gradient of approximately 1.2. Prewetting rate was quantified by weighting the Mariotte bottle during capillary rise. This prewetting resulted in an average rate of 160 mL h-1 (equivalent to 70 mm h-1). The slow prewetting was achieved by saturating the soil column from below using a peristaltic pump at an average flow rate of 10.2 mL h-1 (equivalent to 4.5 mm h-1) until a ponded water level of 1 cm above the soil surface was attained. During prewetting and saturation of the soil columns, the Yizreel and Kedma soils swelled from an initial dry soil column length of 3.94 ± 0.07 and 3.58 ± 0.10 cm, respectively, to 5.26 ± 0.16 and 4.44 ± 0.06 cm at the wet state. Swelling was more pronounced in Yizreel soil because of its higher clay content. The flow was then reversed and the columns were leached with a constant head device. Two hydraulic gradients (2.8 and 11.6 for Yizreel; 3.1 and 13.6 for Kedma) were applied, corresponding with hydrostatically induced pressures of 145 and 625 kg m-2 for the low and high hydraulic gradient treatment, respectively. In order to ensure saturation of the soil samples with exchangeable Ca, the soil columns were leached with 10 to 30 pore volumes of a CaCl2 solution (0.5 M Cl-), which took 1 to 4 h, depending on the prewetting treatment and the soil. However, the conditioning stage (leaching with the 0.5 M Cl- solution) was much longer for Yizreel soil samples exposed to a fast prewetting and low hydraulic gradient since leaching with 10 to 20 pore volumes took
22 h (Moutier et al., 1998). The hydraulic conductivity measured at the end of the conditioning stage was taken as reference hydraulic conductivity (K0.5).
Following leaching with the 0.5 M Cl- solution, the soil columns were leached with
80 pore volumes (
5 L) of dilute CaCl2 solution (0.01 M Cl-) or DW. The 0.01 M Cl- solution was chosen in order to induce some swelling (but no clay dispersion), whereas DW was used to simulate rainwater infiltration. As development of cohesive forces were reported to be time dependent (e.g., Blake and Gilman, 1970; Shainberg et al., 1996), leaching duration was limited to 20 h. However, the treatment was much shorter for soil columns exposed to a high hydraulic gradient and slow prewetting rate as it only took 200 min to infiltrate 80 pore volumes through these columns. In the few cases where steady-state hydraulic conductivities were not maintained, the experimental results were extrapolated to steady-state values, which are presented in Table 2
.
The hydraulic gradient was maintained constant during leaching for a given solution (0.5 M Cl-, 0.01 M Cl-, or DW) except for short periods (510 min) when flow was stopped to change the leaching solution. The 5-L Mariotte bottle with the original solution was first disconnected and the Mariotte bottle with the new solution then reconnected and the hydraulic gradient reestablished. This procedure minimized mixing of the two leaching solutions (Moutier et al., 1998).
The column leachates were collected in tubes with a fraction collector and electrical conductivity (EC), pH, and volume of leachate were measured. Results for both vertisols were compared in terms of relative hydraulic conductivity, Kr, defined as the ratio of the saturated hydraulic conductivity calculated for a given solution (K0.01 or KDW) to the reference hydraulic conductivity calculated at the end of the conditioning stage (K0.5). Up to five replicates were carried out and the results were analyzed by one-way analysis of variance using paired Student t test (comparison of two means with unknown but equal variance). A significance level of 0.05 was chosen for the tests.
 |
Results and discussion
|
|---|
Effects of Prewetting Rate
Mean reference hydraulic conductivities
, mean relative hydraulic conductivities
measured at the end of the corresponding leaching period, and their standard deviations are summarized in Table 2. Hydraulic conductivities
for the fast prewetting treatment of the Yizreel vertisol were not affected significantly by the hydraulic gradient. As expected, the slow prewetting treatment minimized structure degradation and led to significantly higher values when compared with the fast prewetting treatment. Relative hydraulic conductivities for the slow prewetting treatment for the Yizreel soil samples exposed to a high hydraulic gradient were significantly higher
0.5, than those exposed to a low hydraulic gradient for both leaching solutions.
A comparison of
r0.01 as a function of pore volume (Fig. 1a)
and cumulative time (Fig. 1b) is given for Yizreel soils samples exposed to fast and slow prewetting and two different hydraulic gradients. When the 0.5 M Cl- solution was replaced by the 0.01 M Cl- solution
r0.01 decreased, and attained a minimum before increasing as a function of pore volume or time. However, the amplitude and the rate of increase of
r0.01 differed among the treatments. Results for Yizreel soil samples exposed to fast prewetting were discussed elsewhere (Moutier et al., 1998) and only the main conclusions are presented. The decrease in
r0.01 was attributed to the combination of two mechanisms, namely swelling upon dilution of the leaching solution and physical compression of the clay fabric at the outflow end of the soil column (Kemper et al., 1972; Moutier et al., 1998). Physical compression increases with an increase in the hydrostatically induced pressure, and hence, with an increase in hydraulic gradient. The initial decrease was followed by an increase in
r0.01. Increases of the hydraulic conductivities were reported by several authors (Christiansen, 1944; Constantz et al., 1988; Faybishenko, 1995). Faybishenko (1995) attributed his observations to air entrapment and its subsequent dissolution. Our soil samples were prewetted by capillary rise, and only immobile entrapped air remained in the column, representing probably <5% of the total porosity (Faybishenko, 1995). This was also confirmed by flushing soil columns with CO2 prior to the prewetting treatment (data not presented). The small amount of entrapped CO2 dissolved readily during the soil conditioning phase when the soil was leached with the 0.5 M solution. The trends exhibited by
r0.01 (Fig. 1) were similar for both air and CO2 flushing prior to soil saturation, with no significant difference between the two treatments. The experiment with and without CO2 flushing confirmed (i) that the effect of entrapped air following the fast prewetting procedure applied to the soil columns was minimal and (ii) that dissolution of entrapped air was not responsible for the increase in
r0.01. Thus, the development of cohesive bonds between clay fabric structural units and microaggregates in close contact has been proposed to explain the observed increase in
r0.01 (Kemper and Rosenau, 1984; Shainberg et al., 1996; Moutier et al., 1998).

View larger version (25K):
[in this window]
[in a new window]
|
Fig. 1 Relative hydraulic conductivities as a function of (a) pore volume and (b) accumulated time for Yizreel leached with 0.01 M Cl- solutions (data for the fast prewetting treatment were taken from Moutier et al., 1998)
|
|
The hydraulic conductivity of the slow prewetted Yizreel soil columns also depended on the hydraulic gradient (Fig. 1a). The higher the hydraulic gradient, the greater the increase in
r0.01. In the low hydraulic gradient treatment, the effect of aging time on the increase of
r0.01 was less pronounced in slow prewetted columns than in fast prewetted ones. When aggregate breakdown was reduced by slow prewetting, the efficacy of cohesive bond formation and reaggregation was also limited. In the high hydraulic gradient treatment, fast prewetting led to higher
r0.01 values than slow prewetting (Table 2), although absolute values were much lower with the former treatment. K0.01 equaled 2.90 and 16.05 cm h-1 for the fast prewetting, high hydraulic gradient and slow prewetting, high hydraulic gradient treatments, respectively.
The rate of development of interparticle bonding forces, as expressed by the rate of increase in
r0.01, was evaluated by plotting
r0.01 as a function of cumulative time (Fig. 1b). The assumption that prolonged leaching was needed for cohesive forces to develop is not necessarily true (Moutier et al., 1998). A sharp increase in
r0.01 was observed for soil columns exposed to slow prewetting and high hydraulic gradient treatment in spite of the fact that leaching with the 0.01 M Cl- solution lasted only 200 min (Fig. 1b). It is suggested therefore that the rate of development of bonds between clay particles in close proximity (e.g., when high hydraulic gradient is used) is very fast. Conversely, the formation of bonds in soils exposed to low hydraulic gradient is slow because of the relatively large distance between the clay particles.
The observed interaction between the rate of bond formation, as expressed by the increase in
r0.01, and the proximity of clay particles, as determined by the hydrostatically induced pressure, was verified by studying the effect of aging at zero hydraulic gradient. The experiment at zero gradient mimics experiments at zero compaction such as those recently published by Ben-Hur et al. (1998). Following leaching with a 0.01 M Cl- solution (five pore volumes), the slow prewetting treatment for Yizreel soil samples were exposed to a zero gradient and allowed to age for 16 h (960 min). The flow was then reestablished and a high hydraulic gradient of 11.8 was applied (50 min). The rate of increase in
r0.01 in Yizreel soil samples exposed to a zero hydraulic gradient for 16 h was much lower than the subsequent rate of increase when high hydraulic gradient was applied (Fig. 2)
.
r0.01 increased from 1.14 to 1.6 within 16 h (equivalent to a rate of 0.03 h-1) and from 1.6 to 2.2 within the subsequent 50 min (equivalent to a rate of 0.72 h-1) when the high hydraulic gradient was applied. The latter was close to the rate of increase of
r0.01 observed in the slow prewetting, high hydraulic gradient treatment (0.54 h-1) (Fig. 1b). This experiment supported the hypothesis that the rate of hydraulic conductivity increase and the formation of cohesive bonds strongly depended on the mutual proximity of the clay particles. This is also clearly demonstrated in Fig. 1b. For a given prewetting treatment, the rate of hydraulic conductivity increase was higher in columns exposed to high hydraulic gradient than those exposed to low hydraulic gradient. This is in agreement with previously published data (Moutier et al., 1998) and confirms the fact that the rate of particle bond formation depends on the hydraulic gradient; the higher the hydraulic gradient, the closer the clay particles and the higher the rate of cohesive forces development. This confirms recently published findings that the "closeness" of clay particles is a prerequisite for the manifestation of these cohesive forces at the aggregate level (Zubkova, 1998). It also explains some conflicting evidence that cohesive interparticle bonds might develop under saturated conditions where water layers surrounding the clay particles impede close contacts (Blake and Gilman, 1970; Kemper and Rosenau, 1984; Shainberg et al., 1996).

View larger version (15K):
[in this window]
[in a new window]
|
Fig. 2 Relative hydraulic conductivities as a function of accumulated time for Yizreel exposed to slow prewetting and high hydraulic gradient treatment (50 min) following a zero hydraulic gradient treatment (16 h)
|
|
Replacing the 0.01 M Cl- solution with DW resulted in a decrease in
rDW of Yizreel columns, irrespective of the prewetting rate or the hydraulic gradient applied to the soil columns (Fig. 3a and 3b)
. Also presented in Fig. 3b are the quasi-steady state EC values measured at the end of leaching with DW. In the calcareous Yizreel, EC values were affected by CaCO3 dissolution that, in turn, depended on flow velocity. When a high hydraulic gradient was applied or when a high
rDW was maintained because of slow prewetting, a high water flux through the column was maintained. Consequently there was less time for CaCO3 to dissolve, and EC was low (0.035 dS m-1). When the water flow was slowed by fast prewetting and a low hydraulic gradient, there was more time for CaCO3 to dissolve and the EC was higher (0.07 dS m-1). However, in both cases, the total electrolyte concentrations, inferred from the ECs of the leachates1
(Fig. 3b) were higher than the critical flocculation value of Ca-saturated smectite (0.25 mmolc L-1) (Oster et al., 1980), suggesting that swelling was the governing mechanism responsible for the decrease in
rDW. The long duration needed for the decrease in
rDW (Fig. 3b), coupled with the fact that swelling is a slow process compared with clay dispersion (Keren and Singer, 1988), further supported the swelling mechanism.

View larger version (28K):
[in this window]
[in a new window]
|
Fig. 3 Relative hydraulic conductivities as a function of (a) pore volume and (b) accumulated time for Yizreel leached with deionized water (data for the fast prewetting treatment were taken from Moutier et al., 1998)
|
|
After 1200 min of leaching, the Yizreel soil with deionized water,
rDW was approaching a quasi-steady state value of
0.6, in three of the four treatments (Fig. 3b). The near steady state
rDW indicated that the swelling in these columns was almost complete. Only in the soil columns subjected to a slow prewetting and high hydraulic gradient treatment, where the initial hydraulic conductivity was high (K0.01 = 16.05 cm h-1), the leaching period with DW was short (200 min) and
rDW was maintained at 2.2. Slow prewetting prevented aggregate degradation and maintained large conducting pores between the aggregates, whereas a high hydraulic gradient resulted in a short leaching time. Evidently, that peculiar combination prevented swelling from reaching its full extent and resulted in a stable
rDW. The deleterious effect of DW on
rDW in these slowly prewetted, highly smectitic soil samples was probably prevented by the development of inter- and intraaggregate cohesion forces resulting from the high hydraulic gradient, and limited swelling.
Effects of Clay Content
Mean reference hydraulic conductivities
, mean relative hydraulic conductivities
measured at the end of leaching and their standard deviations for the vertisol with a lower clay content (Kedma soil with 46.5% clay) are also summarized in Table 2. The reference values
of the Kedma exposed to fast prewetting were 1.07 and 2.11 cm h-1 under high and low hydraulic gradient treatments, respectively, and were similar to those for Yizreel despite the difference between their clay contents. As expected, slow prewetting resulted in higher
0.5. However, this increase in
0.5 was lower than the corresponding increase for Yizreel (Table 2). The more pronounced effect of prewetting rate on
0.5 of Yizreel was attributed to the more stable structure that its higher clay content imparted to this soil (Kemper and Koch, 1966); the higher the clay content, the more stable the aggregates. Slow prewetting did not cause aggregate breakdown in the Yizreel soil; therefore
0.5 was high. Due to its lower clay content, Kedma soil aggregates were less stable and some slaking probably occurred during the slow prewetting treatment, resulting in the lower
0.5 values. When fast prewetting was applied, aggregate disintegration occurred in both soils, resulting in comparable
0.5 values.
As cohesive forces are more prone to develop in soils with higher clay contents (Singer et al., 1992), the impact of hydraulic gradient on hydraulic conductivity was expected to be less pronounced in Kedma soil samples. All treatments (Fig. 4a)
, with the exception of the slow prewetting, low hydraulic gradient treatment, resulted in trends similar to those of Yizreel (Fig. 1a) but were less pronounced. For the Kedma vertisol, both prewetting treatments resulted in final
r0.01 values that were not significantly affected by the hydraulic gradient.

View larger version (25K):
[in this window]
[in a new window]
|
Fig. 4 Relative hydraulic conductivities as a function of (a) pore volume and (b) accumulated time for Kedma leached with 0.01 M Cl- solutions
|
|
When a high hydraulic gradient was applied following fast prewetting,
r0.01 first decreased. However, the duration of the hydraulic conductivity decrease was shorter in the Kedma than in the Yizreel (10 pore volumes for Kedma compared with 30 pore volumes for Yizreel). This difference was probably due to the more limited swelling in the Kedma with its low clay content. Having a lower clay content (number of clay to clay contacts), the rate of cohesive bond formation was lower for Kedma than for Yizreel. As a result, the initial decrease in
r0.01 was followed by a more gradual increase in the Kedma vertisol (Fig. 4) than the corresponding increase in the Yizreel vertisol (Fig. 1).
Data in Fig. 4a support our hypothesis that the initial structural state of Kedma affects the impact of hydraulic gradient on relative hydraulic conductivity. The increase in
r0.01 for Kedma soil columns exposed to a low hydraulic gradient was more pronounced following a slow prewetting than following a fast prewetting treatment (Fig. 4a). Due to the lower clay content in Kedma soil samples, fast prewetting resulted in a more extensive breakdown of soil aggregates, and the effects of aging and the increase in
r0.01 were less pronounced. Conversely, under slow prewetting, the increase in
r0.01 was more pronounced because aging was more extensive in partially broken aggregates. These results (Fig. 4a) supported a previously published assumption (Levy et al., 1997) that aging is more effective in partially disintegrated aggregates than in aggregates where extensive slaking has taken place. Thus, the increase of
r0.01 was greatest when Kedma columns were exposed to slow prewetting and a low hydraulic gradient. Unlike the other three treatments, this treatment seemed to be more effective in improving
r0.01 in the Kedma than in the Yizreel soil.
The change in
r0.01 with time for Kedma soil is shown in Fig. 4b. The results confirmed the strong dependence between the rate of development of interparticle binding forces, as inferred from the rate of increase in
r0.01, and the mutual proximity of clay particles within the clay fabric, as determined by the hydraulic gradient applied to the soil column for a given treatment. The rate of bond formation increased with increasing hydraulic gradient. Because of the more stable structure resulting from slow prewetting, a similar hydraulic gradient treatment resulted in a higher rate of cohesive bond development than that which followed fast prewetting.
The effect of clay content on the rate of increase in
r0.01 is demonstrated by comparing the results shown in Fig. 1b and 4b. Following slow prewetting, the rate of development of interparticle binding forces increased in the following order: low hydraulic gradient (Yizreel) < low hydraulic gradient (Kedma) < high hydraulic gradient (Kedma) < high hydraulic gradient (Yizreel). The lowest and the highest rates were both ascribed to the high clay content of Yizreel. High clay content contributed to two opposing tendencies: (i) a limited effect of aging because of the more stable clay fabric and (ii) an enhanced number of potential clay to clay contacts, which increases with increasing hydraulic gradient. Because of its lower clay content, Kedma exhibited intermediate rates of bond formation for similar hydraulic gradient treatments. Following fast prewetting, the effect of hydraulic gradient on the rate of cohesive bond formation was more pronounced in Yizreel (Fig. 1b) than in Kedma because fast prewetting resulted in extensive slaking, which affected the Kedma to a greater extent, and the rate of cohesive bond formation was lower in the latter (Fig. 4b).
When the 0.01 M Cl- solution was replaced with DW,
rDW of the Kedma decreased irrespective of the prewetting rate or the hydraulic gradient applied to the soil columns (Fig. 5)
. The quasi-steady state ECs of the effluent are presented in Fig. 5b. As in the case of the Yizreel soil, the high EC in the Kedma was due to CaCO3 dissolution and depended on flow rates. Despite similarities in the ECs of the effluents when the two soils were leached with deionized water, the following should be noted:
- The fast prewetting and high hydraulic gradient treatment resulted in comparable relative hydraulic conductivity values in the two soils (Table 2). Evidently, swelling affected both soils in a similar way. Because conditions for cohesive forces development were not favorable when Kedma was subjected to the fast prewetting, low hydraulic gradient treatment,
rDW decreased to 0.14.
- In the slow prewetting, high hydraulic gradient treatment, the hydraulic conductivity of the Yizreel soil at the end of the leaching with DW was similar to that at the end of the leaching with the 0.01 M Cl- solution, and more than twice the initial value (Table 2). In the Kedma,
rDW decreased to 0.85. Following the slow prewetting treatment, aggregate breakdown in Yizreel was much less pronounced than in Kedma. During leaching with DW the adverse effect of swelling in this soil was probably offset by the development of cohesive forces, and the net effect was that the size of the water-conducting pores was not affected. Conversely, in the Kedma, in which some aggregate slaking took place during the slow prewetting, swelling of the microaggregates probably reduced the size of the conducting pores, and
rDW decreased.
- In the slow prewetting, low hydraulic gradient treatment, the
rDW of the Yizreel soil decreased to 0.61 (Table 2). Under a low hydraulic gradient, formation of cohesive bonds between clay particles is not as effective; therefore, because of swelling,
rDW decreased. Similar conditions prevailed in the Kedma, and the final
rDW in the two soils reached similar values.

View larger version (30K):
[in this window]
[in a new window]
|
Fig. 5 Relative hydraulic conductivities as a function of (a) pore volume and (b) accumulated time for Kedma leached with deionized water
|
|
 |
Summary and conclusions
|
|---|
We investigated the effect of hydrostatical pressure on K for two Ca-saturated vertisols of contrasting initial structure. Under conditions that led to limited aggregate breakdown (i.e., low clay content and slow prewetting, or high clay content and fast prewetting) development of cohesive forces between structural units was suggested to explain the increase in hydraulic conductivity. The rate of increase in K also depended on the hydraulic gradient. Under greater hydraulic gradient, the structural units of the soils moved closer together, leading to greater number of clay to clay contacts and to a faster increase in the hydraulic conductivity. Conversely, when a stable soil structure was maintained (i.e., high clay content and slow prewetting), the development of cohesive bonds was of minor importance in increasing K. Understanding the combined role of intrinsic (e.g., clay content) and extrinsic (e.g., prewetting rate and hydraulic gradient) soil conditions on K will allow development of better soil and irrigation strategies for sustaining soil hydraulic properties. Maintaining a stable structure is essential for ensuring high K in vertisols. This can be obtained by prevention of soil structure degradation or enhancing structural stability. Structure degradation can be reduced by using a combination of zero-tillage and irrigation methods that cause slow prewetting (e.g., microdrippers, microsprinklers). Structural stability can be augmented by enabling conditions favorable for the development of cohesive forces (zero tillage). Extrinsic factors, such as hydrostatic or mechanical pressure or partial drying, decrease the distance among soil particles and enhance the development of cohesive bonds that stabilize soil structure.
 |
ACKNOWLEDGMENTS
|
|---|
The first author is grateful to the Agricultural Research Organization, Bet Dagan, Israel; the Catholic University of Louvain, Louvain-la-Neuve, Belgium; and the Rotary Club of Wavre-Europe, Belgium, in providing her with grants enabling the realization of this work. This study was supported by grant no. TA-MOU-96-CA-16-016 U.S.-Israel Cooperative Development Research Program, Office of the Science Advisor, U.S. Agency for International Development (AID).
 |
NOTES
|
|---|
Contribution from the Agricultural Research Organization, The Volcani Center, P.O. Box 6, Bet Dagan 50-250, Israel. no 604/99 series.
1 The total electrolyte concentration (TEC) can be inferred from the EC using the following approximation: TEC (mmolc L-1) = 10 EC (dS m-1). 
Received for publication December 8, 1998.
 |
REFERENCES
|
|---|
- Abu-Sharar T.M., Bingham F.T., Rhoades J.D. Reduction in the hydraulic conductivity in relation to clay dispersion and disaggregation. Soil Sci. Soc. Am. J. 1987;51:342-346.[Abstract/Free Full Text]
- Banin and Amiel. 1970. A correlative study of the chemical and physical properties of a group of natural soils of Israel. Geoderma 3:185198.
- Ben-Dor E., Banin A. Determination of organic content in arid-zone soils using a simple "loss on ignition" method. Soil Sci. Plant Anal. 1989;20:1675-1695.
- Ben-Hur M., Agassi M., Keren R., Zhang J. Compaction, aging, and raindrop-impact effects on hydraulic properties of saline and sodic vertisols. Soil Sci. Soc. Am. J. 1998;62:1377-1383.[Abstract/Free Full Text]
- Blackmore A.V., Miller R.D. Tactoid size and osmotic swelling in calcium montmorillonite. Soil Sci. Soc. Am. J. 1961;25:169-173.[Abstract/Free Full Text]
- Blake G.R., Gilman R.D. Thixotropic changes with aging of synthetic soil aggregates. Soil Sci. Soc. Am. Proc. 1970;34:561-564.
- Boiffin J. La dégradation structurale des couches superficielles du sol sous l'action des pluies. Paris, Grignon: Thèse de Docteur-Ingénieur. Institut National Agronomique, 1984.
- Bresson L.M., Boiffin J. Morphological characterization of soil crust development stages on an experimental field. Geoderma 1990;47:301-325.
- Bresson L.M., Moran C.J. Structural change induced by wetting and drying in seedbeds of hardsetting soil with contrasting aggregate size distribution. Eur. J. Soil Sci. 1995;46:205-214.
- Christiansen J.O. Effect of entrapped air upon the permeability of soils. Soil Sci. 1944;58:355-365.
- Constantz J., Herkelrath W.N., Murphy F. Air encapsulation during infiltration. Soil Sci. Soc. Am. J. 1988;52:10-16.[Abstract/Free Full Text]
- De Backer L.W. The measurement of entrapped gas in the study of unsaturated flow phenomena. Water Resour. Res. 1967;3:245-249.
- Faybishenko B.A. Hydraulic behavior of quasi-saturated soils in the presence of entrapped air: laboratory experiments. Water Resour. Res. 1995;31:2421-2435.
- Kemper W.D., Koch E.J. Aggregate stability of soils from western U.S. and Canada. Washington, DC: U.S. Gov. Print. Office, 1966 USDA Tech. Bull. 1355..
- Kemper W.D., Shainberg I., Quirk J.P. Swelling pressures, electric potentials and ion concentrations: Their role in hydraulic and osmotic flow through clays. Soil Sci. Soc. Am. J. 1972;36:229-236.
- Kemper W.D., Rosenau R.C. Soil cohesion as affected by time and water content. Soil Sci. Soc. Am. J. 1984;48:1001-1006.[Abstract/Free Full Text]
- Keren R., Singer M.J. Effect of low electrolyte concentration on hydraulic conductivity of sodium/calcium montmorillonitesand system. Soil Sci. Soc. Am. J. 1988;52:368-373.[Abstract/Free Full Text]
- Klute, A. (ed.) 1986. Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Le Bissonnais, Y. 1988. Analyse des mécanismes de désagrégation et de la mobilisation des particules de terre sous l'action des pluies. Thèse de doctorat, Univ. d'Orleans, France.
- Le Bissonnais Y. Contribution à l'étude de la dégradation structurale superficielle: analyse des processus de microfissuration des agrégats par l'eau. Sci. Sol. 1989;27:187-199.
- Levy G.J., Levin J., Shainberg I. Prewetting rate and aging effects on seal formation and interrill soil erosion. Soil Sci. 1997;162:131-139.
- Moutier M., Shainberg I., Levy G.J. Hydraulic gradient, aging and water quality effects on the hydraulic conductivity of a vertisol. Soil Sci. Soc. Am. J. 1998;62:1488-1496.[Abstract/Free Full Text]
- Murray R.S., Quirk J.P. Interparticle forces in relation to the stability of soil aggregates. In: De Boodt M.F., et al. , ed. Soil colloids and their associations in aggregates. New York: Plenum Press, 1990:439-461.
- Oster J.D., Shainberg I., Wood J.D. Flocculation value and gel structure of sodium/calcium montmorillonite and illite suspensions. Soil Sci. Soc. Am. J. 1980;44:955-959.[Abstract/Free Full Text]
- Page, A.L., R.H. Miller, and D.R. Keeney (ed.) 1986. Methods of soil analysis. Part 2. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Panabokke C.R., Quirk J.P. Effect of initial water content on stability of soil aggregates in water. Soil Sci. 1957;83:185-195.
- Quirk J.P., Schofield R.K. The effect of electrolyte concentration on soil permeability. J. Soil Sci. 1955;6:163-178.
- Shainberg I., Goldstein D., Levy G.J. Rill erosion dependence on soil water content, aging and temperature. Soil Sci. Soc. Am. J. 1996;60:916-922.[Abstract/Free Full Text]
- Singer M.J., Southard R.J., Warrington D.N., Janitzky P. Stability of synthetic sandclay aggregates after wetting and drying cycles. Soil Sci. Soc. Am. J. 1992;56:1843-1848.[Abstract/Free Full Text]
- Zubkova T.A. The nature of mechanical strength of oven-dry aggregates of soil. Eurasian Soil Sci. 1998;31:254-262.
This article has been cited by other articles:

|
 |

|
 |
 
G. J. Levy, D. Goldstein, and A. I. Mamedov
Saturated Hydraulic Conductivity of Semiarid Soils: Combined Effects of Salinity, Sodicity, and Rate of Wetting
Soil Sci. Soc. Am. J.,
April 11, 2005;
69(3):
653 - 662.
[Abstract]
[Full Text]
[PDF]
|
 |
|