HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Free Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Lu, Z. Articles by Zhang, D. Search for Related Content PubMed Articles by Lu, Z. Articles by Zhang, D. GeoRef GeoRef Citation Agricola Articles by Lu, Z. Articles by Zhang, D. Related Collections Flow Other Models Vadose Zone Processes and Chemical Transport

Stochastic Analysis of Transient Flow in Heterogeneous, Variably Saturated Porous Media

The van Genuchten–Mualem Constitutive Model

Hydrology, Geochemistry, and Geology Group (EES-6), MS T003, Los Alamos National Laboratory, Los Alamos, NM 87545
* Corresponding author (donzhang{at}lanl.gov)

Received 31 January 2002.

In this study, on the basis of the van Genuchten–Mualem constitutive relationship, we develop a general nonstationary stochastic model for transient, variably saturated flow in randomly heterogeneous media with the method of moment equations. We first derive partial differential equations governing the statistical moments of the flow quantities by perturbation expansions and then implement these equations under general conditions with the method of finite differences. The nonstationary stochastic model developed is applicable to the entire domain of bounded, multidimensional vadose zones or integrated unsaturated–saturated systems in the presence of random or deterministic recharge and sink–source and in the presence of multiscale, nonstationary medium features. We demonstrate the model with some two-dimensional examples of unsaturated and integrated unsaturated–saturated flows. The validity of the developed stochastic model is confirmed through high-resolution Monte Carlo simulations. We also investigate the relative contributions of the soil variabilities (K_S, , and n) as well as the variability in recharge Q to the pressure head variance. It is found that the pressure head variance is sensitive to these variabilities, in the order of n, , K_S, and Q. Though the variability of and n is usually smaller than that of K_S and Q, their effect on the pressure head variance should not be ignored.

Abbreviations: MC, Monte Carlo approach • ME, moment equation–based approach