|
|
||||||||
a Department of Earth and Environmental Sciences, California State University, Fresno, CA 93740
b Department of Environmental Science, University of California, Riverside, CA 92521
* Corresponding author (zwang{at}csufresno.edu)
Received 19 November 2002.
ABSTRACT
Two- and three-dimensional laboratory experiments were conducted in a large 1 by 1 m2 HeleShaw cell and a 10-cm cylinder column to study the stability of the redistribution process in a uniformly packed porous medium of coarse sand. Our results demonstrate that fingers form and propagate rapidly as soon as infiltration stops when the soil is initially dry, but form more slowly and are larger when the soil is wet. Finger widths ranged from about 4.5 cm when the soil was dry to 17 cm when the soil was very wet, and the fraction of the cross-sectional area occupied by fingers also increased (from 36 to 66%). Finger velocities declined with time in all studies, with averages ranging from 3.6 to 9.6 cm min-1. We demonstrated that ponded infiltration was stable and uniform in our system, so that the fingering we observed was not due to soil heterogeneity. The porous medium retained a memory of the fingers formed in the first experiment, so fingers formed in subsequent redistribution cycles followed the old finger paths, even after 28 d had elapsed. Fingers provide channels for rapid drainage of previously infiltrated water, especially when the soil ahead of the front is dry. Our findings clearly contradict predictions made by the Richards equation, which calculates stable flow during redistribution in homogeneous soil.
ALTHOUGH UNSTABLE FLUID FLOW has been acknowledged to exist in field soils, it has generally been associated with specific conditions that interfere with the normal surface tension and viscous flow that stabilizes water as it moves through a porous medium. Aside from these extreme conditions, preferential flow of water and chemicals has generally been attributed to soil heterogeneity and not to instability within the fluid phase (White, 1985). However, it has long been known in fluid mechanics and petroleum engineering (e.g., Taylor, 1950; Saffman and Taylor, 1958, Chuoke et al., 1959) that under certain conditions fluid flow can become unstable in a perfectly uniform porous medium or even when one fluid flows through another in a medium without a porous solid framework. These concepts have been extended to hydrologic applications (e.g., Raats, 1973; Philip, 1975; Parlange and Hill, 1976) to predict the conditions under which unstable flow might occur in porous media. Once the flow becomes unstable, fingering will bypass a significant portion of the matrix and transport the invading fluid much faster and deeper than would occur under stable flow conditions. A thorough review of experimental and theoretical studies of unstable flow in the vadose zone is given in de Rooij (2000) and Nieber (2001).
Previous research on unstable flow in vadose zone hydrology has identified some of the important, necessary conditions for the onset of fluid instability. The most widely studied and documented occurrences of unstable flow in the vadose zone are of two types: vertical flow from a fine-textured layer into a coarse one (Hill and Parlange, 1972; Starr et al., 1978, 1986; Glass et al., 1988, 1989a,b; Hillel and Baker, 1988; Baker and Hillel, 1990, Nieber, 1996; Wang et al., 1998b) and infiltration into a hydrophobic medium (Van Ommen et al., 1988; Hendrickx et al., 1993; Ritsema et al., 1993; Burcar et al., 1994; Carrillo et al., 2000; Wang et al., 1998b, 2000b). Because so much of the attention has been focused on these cases, it is commonly believed that unstable flow is a result either of heterogeneity or soil water repellence. However, infiltration experiments in perfectly uniform porous media have also shown unstable flow induced by air compression ahead of the wetting front (White et al., 1976, 1977; Wang et al., 1997, 1998c) or unsaturated infiltration under low rates (Selker et al., 1992; Wang et al., 1998b; Geiger and Durnford, 2000). Another condition predicted by Raats (1973) and Philip (1975) to cause unstable flow in uniform porous media is redistribution following the cessation of infiltration. It has received relatively little attention in the literature, presumably because the redistribution water front is not considered to be very different from an infiltration front. However, because redistribution is one of the most common natural hydrologic processes, it is important to understand the extent to which the process may be unstable, and what factors influence its occurrence.
Although Nicholl et al. (1994) observed fingering during redistribution in initially dry fractures, Diment and Watson (1985) and Tamai et al. (1987) seem to be the only researchers who have demonstrated that redistribution causes unstable flow in extremely coarse-textured uniform materials. However, these authors concluded that instability only occurs when the medium ahead of the draining front is extremely dry. Diment and Watson's experiments, which were performed in a small HeleShaw cell (80 cm high, 40 cm wide, and 1.5 cm thick), showed that redistribution "stabilized" when the initial water content was increased to only a few percent of saturation. In addition, they performed a stability analysis on the Richards flow equation and showed that instability was extremely rare during redistribution (Diment and Watson 1985). It is possible, however, that their experimental apparatus might have been too small to observe unstable flow under a range of conditions. Fluid stability analyses (Saffman and Taylor, 1958; Chuoke et al., 1959; Philip, 1975) have shown that unstable flow fingers have a characteristic width and spatial frequency or wavelength, and, therefore, that their development might be suppressed if the lateral extent of the experimental apparatus is not sufficiently large. This might explain why studies of redistribution conducted in small columns have not produced unstable flow, because the cross-sectional area may have been too small to display several finger wavelengths. In addition, since fingers are predicted to widen when the initial soil water content increases (Jury et al., 2003), the spatial scale required to observe fingering in moist soil will also increase. Our recent study of redistribution in nonstructured field soils (Wang et al., 2003) showed that instabilities developed after infiltration ceased and that the fingers we observed were quite large (
13-cm diam.) and had a mean wavelength of nearly 30 cm. Therefore, it appears to be necessary to use a large observation frame to study unstable flow during redistribution.
The objectives of the study reported here were (i) to conduct experiments in a large HeleShaw cell to observe the flow patterns during redistribution following infiltration, (ii) to observe the appearance of fingering during subsequent cycles of redistribution if flow was unstable, (iii) to study the effects of high initial water content on the development of fingers, and (iv) to compare the results of two-dimensional experiments with three-dimensional observations. We also report digital pictures of fingering during redistribution. Detailed measurements of the processes and theoretical analyses of the mechanisms of fingering are presented in a companion paper (Jury et al., 2003).
MATERIALS AND METHODS
Laboratory Setup and Materials
A HeleShaw cell with transparent walls was constructed using two sheets of 100 by 100 cm2 Plexiglas bolted together with 1-cm spacing. Gridlines creating 10 by 10 cm2 squares were drawn on the front and back panels for delineating the wetting patterns. The bottom of the cell was open to the atmosphere through nine 3-mm holes to prevent air compression. The cell was filled with commercial silica sand that was first sieved to obtain a uniform 0.5- to 0.8-mm particle size. Before packing, the slab box was laid horizontal with the front (top) panel removed, and dry sand was poured directly on top of the bottom panel using a flat pan. After filling the sand to about 1-cm thickness over 99 cm of the panel, the front panel was mounted and bolted, and the slab box was set vertically upright. The box was then gently tapped using a rubber hammer until the sand was consolidated to 95 cm. The nonuniform surface layer of loose sand was removed using a vacuum, and the remaining sand was leveled to its final height of 91 cm.
The saturated hydraulic conductivity of the sand, estimated using a constant head permeameter, was Ks = 682 cm h-1. The water-entry value of the sand was about hwe = 3 cm, measured using a tension-pressure permeameter (Wang et al., 2000a). We measured the drying loop of the moisture release curve of the sand used in our experiments with a hanging water column (Jury et al., 1991), as shown in Fig. 1. A V-shaped funnel flume was constructed to hold and release a maximum of 2000 mL of water through a slotted tube attached to the bottom of the flume, which produced a uniform application of ponded water to the sand surface in 2 to 5 s.
|
|
RESULTS
Figures 2 and 3 show photographs of the wetting front taken at different times during each of the experiments. The redistribution time is indexed from t = 0, defined as the moment when ponded water first disappeared from the sand surface. By comparing these photographs and relating the observations to the experimental conditions, we can draw a number of inferences about the nature of the flow process during infiltration and redistribution.
|
|
The finger widths we observed are consistent with sizes predicted by the equation
![]() | [1] |
The great depths reached by the fingers following such a small water application (5 cm) are also understandable from the soil properties. As shown in Jury et al. (2003), maximum finger depth can be estimated approximately by assuming that the final profile is in hydrostatic equilibrium. The initial water infiltration I produces a water application W to the fingers that increases to I/
, where
is the fraction of the area covered by fingers. In our study,
0.36, so I/
= W = 8.3 cm. At equilibrium, the finger will have a water content distribution at heights z above the tip which is given by
eq[h* - z], where
eq[h] is the equilibrium drying loop of the water characteristic function, and h* is the matric potential of the drying loop corresponding to the water-entry value of the wetting loop. Thus, the finger will reach a depth zmax, which is given by
![]() | [2] |
The drying loop of the water characteristic function of our coarse sand is shown in Fig. 1, along with the best fit (
s = 0.371,
r = 0.039,
= 0.07, N = 7.5) of the van Genuchten (1980) parameters
![]() | [3] |
Using our best estimate h*
12 cm, we calculated that 5 cm of water would penetrate to a depth zmax approximately 115 cm below the depth where the fingers first form. Thus, the small amount of water storage behind the tip of the finger is responsible for its deep penetration. Fingers formed in a finer-textured soil might penetrate only a few centimeters following a 5-cm water application, both because of the higher water storage in the finger and because the fingers would be wider (Jury et al., 2003).
Finger Recurrence
The second experiment, shown in Fig. 2g through 2l, was designed to see how the residual traces of the fingers formed in the previous irrigation cycle would affect the wetting front during a subsequent irrigation. As soon as infiltration ended at t = 0, water started to funnel into the old finger paths. The wetting front of the region between the fingers moved only a small distance below the point reached by the first cycle, and once again halted its advance as soon as the fingers were reached. Two additional fingers were created in this cycle, one nearest the left side of the box and the other at the third position from the right side. Also, the area fraction occupied by fingers rose to about 46%. The finger size increased slightly to a mean of 5.1 ± 2.2 cm, and the finger speed decreased by nearly 50%, as shown in Table 1. More water drained out of the box through the fingers in this study than in the first one.
Porous Medium Memory of Fingers
The third study (Fig. 3) was virtually a replication of the second, conducted 4 d after Exp. 2 ended. The development, location, and speed of the fingers were all quite similar to the processes observed in Exp. 2, except for an increase (to 6.8 ± 2.2 cm) in the mean size of the fingers, and a growth (to 54%) in the finger area fraction (Table 1). The fourth experiment was conducted to determine if the porous medium retained a memory of the fingered flow paths for an extensive period of time. During the 28 d that elapsed between the end of Exp. 3 and the beginning of Exp. 4, water in the finger domains gradually diffused into the surrounding profile, erasing the visual memory of the flow paths, as shown in Fig. 3f. Despite the apparent homogeneity of the soil ahead of the front, however, the 5-cm irrigation added in Exp. 4 followed virtually the same flow paths that had been produced by the fingers in the previous cycles. Figure 4 shows a closeup view of the infiltration pattern from Fig. 3d at 16 min of the third cycle, and from Fig. 3j at 18 min into the cycle conducted 28 d later. The fingers from the second cycle were slightly larger (9.2 ± 2.3 cm) and occupied about 66% of the cross-sectional area. Consequently, the finger speed (see Table 1) was reduced significantly compared with earlier experiments. Thus, the effect of the long redistribution time was to widen the finger pathways but not to eliminate them.
|
|
The behavior in experiments 4 and 6 directly contradicts the results observed by Diment and Watson (1985) in their smaller HeleShaw cell; they noted that an increase to only a few percent water content eliminated the instability they had observed in oven-dry soil. According to their stability analysis, our results also contradict the predictions of the Richards equation.
Finger Speed
The dynamic changes of the wetting front velocity in all the experiments are shown in Fig. 6. The finger velocity, even in the earliest stages, is significantly less than the velocity of the ponded wetting front, indicating that the fingers are unsaturated. Finger velocity declined with time, as predicted by the model of Jury et al. (2003), and indicates that the rate of supply from the matrix is controlling the advance of the fingers. Velocity decreased as initial water content increased (Exp. 4 and 6), which may reflect the increase in area fraction associated with these studies. During redistribution, the downward flow in the matrix region between the fingers essentially stopped, which indicates that the pressure at the wetting front has dropped below the water-entry value required to continue its advance. Once the fingers fully develop, water behind the matrix front between the fingers experiences a lateral pressure gradient and reinforces the finger flow, thereby lowering the matrix water pressure below the threshold for downward movement (Jury et al., 2003). When the soil was initially dry, the wetted matrix region between the fingers did not move beyond 10 cm below the surface, and the fingers moved to the bottom of the cell, producing significant drainage (see Table 1). However, for the initially wet sand in Exp. 6, the wetting front remained stable for a longer time and moved to 42 cm below the surface before splitting.
|
|
CONCLUSION
Our laboratory experiments showed that the redistribution process following the cessation of infiltration is unstable in homogeneous soil, producing a series of fingers that move far ahead of the original front. We demonstrated that ponded infiltration is stable and uniform in our system, so that the fingering observed is not due to soil heterogeneity. The porous medium retained a memory of the fingers formed in the first experiment, so that fingers formed in subsequent redistribution cycles followed the old finger paths, even after 28 d had elapsed. Fingers provide channels for rapid drainage of previously infiltrated water, especially when the soil ahead of the front is dry.
Contrary to a previous observation by Diment and Watson (1985) that the wetting front is stabilized during redistribution when the initial moisture content was raised above 1%, our experiment in very wet sand that had been saturated 24 h earlier still resulted in unstable flow. The fingers observed in this case were significantly larger than when the soil was dry. In the three-dimensional experiments, we also witnessed unstable flow during redistribution after the cessation of saturated infiltration, resulting in fingers propagating in the center or along the wall of a small column.
Our findings clearly contradict predictions made by Richards equation, which calculates stable flow during redistribution in homogeneous soil. We speculate (Jury et al., 2003) that the reason for this failure is the absence of a water-entry matric potential in the continuum description of flow.
ACKNOWLEDGMENTS
The research was supported by Research Grant No. IS-2859-97 from BARD, the U.S.-Israel Binational Agricultural Research and Development Fund.
REFERENCES
This article has been cited by other articles:
![]() |
T. Annaka and S. Hanayama Pressure Head Profile within Growing Fingers in Initially Dry Glass Beads Soil Sci. Soc. Am. J., May 16, 2007; 71(3): 901 - 908. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Nicholl and R. J. Glass Infiltration into an Analog Fracture: Experimental Observations of Gravity-Driven Fingering Vadose Zone J., November 11, 2005; 4(4): 1123 - 1151. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Cho, G. H. de Rooij, and M. Inoue The Pressure Head Regime in the Induction Zone During Unstable Nonponding Infiltration: Theory and Experiments Vadose Zone J., September 12, 2005; 4(4): 908 - 914. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Wang, W. A. Jury, A. Tuli, and D.-J. Kim Unstable Flow during Redistribution: Controlling Factors and Practical Implications Vadose Zone J., May 1, 2004; 3(2): 549 - 559. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. A. Jury, Z. Wang, and A. Tuli A Conceptual Model of Unstable Flow in Unsaturated Soil during Redistribution Vadose Zone J., February 1, 2003; 2(1): 61 - 67. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Plant Registrations | Soil Science Society of America Journal | ||||
| Journal of Natural Resources and Life Sciences Education |
Journal of Environmental Quality |
||||