VZJ
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 20 November 2007
Published in Vadose Zone J 6:925-934 (2007)
DOI: 10.2136/vzj2006.0141
© 2007 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Full Text Free
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zhang, Z. F.
Right arrow Articles by Ward, A.L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Zhang, Z. F.
Right arrow Articles by Ward, A.L.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Zhang, Z. F.
Right arrow Articles by Ward, A.L.
Related Collections
Right arrow Multi-Fluid Systems

Saturation-Dependent Hydraulic Conductivity Anisotropy for Multifluid Systems in Porous Media

Z. F. Zhang, M. Oostrom* and A.L. Ward

Hydrology Group, Environmental Technology Division, Pacific Northwest National Lab., P.O. Box 999, Richland, WA 99352

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

 
FIG. 1. Schematic of the water and nonaqueous phase liquid (NAPL) pressure head variation with time at the top and bottom boundaries of the simulation domain. For each simulation set, the water pressure head was kept constant (–2.0, –1.75, –1.5, and –1.25 m) while the NAPL pressure was lowered in steps from zero to a value slightly larger than water pressure head.

 

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

 
FIG. 2. Nonaqueous phase liquid (NAPL) retention curves of anisotropic soils (R = 50) with different heterogeneity levels (Soils 3, 7, 11, and 15 in Table 1). R, ratio of correlation length at the direction parallel to soil strata and that normal to soil strata; {alpha}, inverse macroscopic capillary length; n, pore-size distribution parameter; Y = ln(Ks), with Ks being the saturated hydraulic conductivity; {sigma}Y2, variance of Y.

 

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

 
FIG. 3. Relationship of (a) {alpha} vs. {sigma}Y2 and (b) n vs. {sigma}Y2 for the total liquid retention curves for an air–water (data from Zhang et al., 2003) and air–NAPL–water system. {alpha}, inverse macroscopic capillary length; n, pore-size distribution parameter; Y = ln(Ks) with Ks being the saturated hydraulic conductivity; {sigma}Y2, variance of Y.

 

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

 
FIG. 4. The ratios of (a) parameter {alpha} and (b) parameter n of an air–NAPL–water system to those of an air-water system (data from Zhang et al., 2003). The solid line in each plot represents the average. {alpha}aw, {alpha} parameter for an air–water system; {alpha}ao, {alpha} parameter for an air–NAPL–water system; Y = ln(Ks) with Ks being the saturated hydraulic conductivity; {sigma}Y2, variance of Y.

 

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

 
FIG. 5. Nonaqueous phase liquid (NAPL) permeability as a function of total fluid saturation for anisotropic heterogeneous soils 9–12 in Table 1 and for hw = –2.0 m. Lp, connectivity–tortuosity coefficient in the direction parallel to strata; Ln, connectivity–tortuosity coefficient in the direction normal to strata; R, ratio of correlation lengths of the direction parallel and normal to strata; Y = ln(Ks) with Ks being the saturated hydraulic conductivity; {sigma}Y2, variance of Y.

 

Figure 6
View larger version (27K):
[in this window]
[in a new window]

 
FIG. 6. Nonaqueous phase liquid (NAPL) anisotropy coefficients as a function of NAPL saturation of soils 9–12 in Table 1 for imposed water pressure head of –2.0 m. Lines: computed relations using Eq. [16] with permeabilities listed in Table 2 and L factors values in plots. Circles: ratio of NAPL permeability values for simulations with flow direction parallel to those normal to soil strata.

 

Figure 7
View larger version (23K):
[in this window]
[in a new window]

 
FIG. 7. Relationships between L vs. ln(R*) of the soils with different levels of heterogeneity. L, tortuosity coefficient; R*, ratio of correlation length in the direction parallel and normal to flow; {sigma}Y2, variance of Y, where Y = ln(Ks). with Ks being the saturated hydraulic conductivity.

 

Figure 8
View larger version (10K):
[in this window]
[in a new window]

 
FIG. 8. Comparisons of the (a) intercepts and (b) slopes of the relationship between L vs. ln(R*) for an air–water (Zhang et al., 2003) and air–NAPL–water system. {sigma}Y2, variance of Y, where Y = ln(Ks), with Ks being the saturated hydraulic conductivity.

 





HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Journal of Natural Resources
and Life Sciences Education
Soil Science Society of America Journal
Journal of Plant Registrations Journal of
Environmental Quality
The Plant Genome
Copyright © 2007 by the Soil Science Society of America.