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 Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (2) Citing Articles via Google Scholar Google Scholar Articles by Tartakovsky, A. M. Articles by Tartakovsky, D. M. Search for Related Content PubMed Articles by Tartakovsky, A. M. Articles by Tartakovsky, D. M. GeoRef GeoRef Citation Agricola Articles by Tartakovsky, A. M. Articles by Tartakovsky, D. M. Related Collections Spatial Variability Analytical Solutions

SPECIAL SECTION: UNDERSTANDING SUBSURFACE FLOW AND TRANSPORT PROCESSES AT THE IDAHO NATIONAL ENGINEERING & ENVIRONMENTAL LABORATORY (INEEL) SITE

Transient Flow in a Heterogeneous Vadose Zone with Uncertain Parameters

^a Geosciences Research Group, Idaho National Engineering and Environmental Laboratory, P.O. Box 1625, MS 2025, Idaho Falls, ID 83415-2025
^b Department of Mathematics, University of Arizona, Tucson, AZ 85721-0089
^c Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545
^* Corresponding author (tartam{at}inel.gov).

Received 3 June 2003.

We consider transient flow in unsaturated heterogeneous porous media with uncertain hydraulic parameters. Our aim is to provide unbiased predictions (estimates) of system states, such as pressure head, water content, and fluxes, and to quantify the uncertainty associated with such predictions. We achieve this goal by treating hydraulic parameters as random fields and the corresponding flow equations as stochastic. Current stochastic analyses of transient flow in partially saturated soils require linearization of the constitutive relations, which may lead to significant inaccuracies when these relations are highly nonlinear. If relative conductivity and saturation vary exponentially with pressure and the corresponding scaling parameters are random variables, the transient Richards equation is mapped onto a linear equation by means of the Kirchhoff transformation. This allows us to develop deterministic differential equations for the first and second ensemble moments of pressure and saturation. We solve these equations analytically, for vertical infiltration, and compare them with direct Monte Carlo simulations.

Abbreviations: INEEL, Idaho National Engineering and Environmental Laboratory • MCS, Monte Carlo simulation • MDE, moment differential equation • PDF, probability density function

This article has been cited by other articles:

				A. M. Tartakovsky, D. Bolster, and D. M. Tartakovsky Hydrogeophysical Approach for Identification of Layered Structures of the Vadose Zone from Electrical Resistivity Data Vadose Zone J., November 26, 2008; 7(4): 1207 - 1214. [Abstract] [Full Text] [PDF]