Stochastic analysis of saturated – unsaturated flow in heterogeneous media by combining Karhunen-Loeve expansion and perturbation method Jinzhong Yang a , Dongxiao Zhang b,c, * , Zhiming Lu b a National Key Laboratory of Water Resources and Hydropower Engineering Sciences, Wuhan University, Wuhan, China b Hydrology, Geochemistry, and Geology Group (EES-6), Los Alamos National Laboratory, Los Alamos, NM 87545, USA c Mewbourne School of Petroleum & Geological Eng, The University of Oklahoma, Norman, OK73019, USA Received 8 April 2003; revised 30 September 2003; accepted 20 October 2003 Abstract In this study, a stochastic model for transient saturated –unsaturated flow is developed based on the Karhunen-Loeve expansion of the input random soil properties combined with a perturbation method. The log-transformed saturated hydraulic conductivity f ðxÞ and the soil pore size distribution parameter aðxÞ are assumed to be normal random functions with known covariances. We decompose f ðxÞ and aðxÞ as infinite series in a set of orthogonal normal random variables by the Karhunen- Loeve expansion and expand the pressure head as polynomial chaos with the same set of orthogonal random variables. The perfectly correlated and uncorrelated cases between f ðxÞ and aðxÞ are studied. By using the Karhunen-Loeve expansion of the input random parameters, polynomial chaos decomposition of pressure head, and the perturbation method, the saturated – unsaturated flow equation and the corresponding initial and boundary conditions are represented by a series of partial differential equations in which the dependent variables are the deterministic coefficients of the polynomial chaos expansion. Once the partial differential equations are solved subsequently by a numerical method, the random representation of pressure head is obtained by combining the deterministic coefficients obtained and the random variables from the Karhunen-Loeve expansion of the input random functions. The moments of pressure head and water content are determined directly from the random representation of the pressure head. We demonstrated the applicability of the proposed KL-based stochastic model with some examples of unsaturated and saturated – unsaturated flow in two dimensions, and compared the results with those from the moment-based stochastic model. It is shown that the KL-based models are more computationally efficient than the conventional moment-based models. q 2004 Elsevier B.V. All rights reserved. Keywords: Stochastic modeling; Saturated– unsaturated flow; Karhunen-Loeve expansion; Perturbation method; Numerical analysis; Vadose zone 1. Introduction It has been recognized that medium heterogeneity significantly impacts fluid flow and solute transport in the subsurface. The vadose zone connects Journal of Hydrology 294 (2004) 18–38 www.elsevier.com/locate/jhydrol 0022-1694/$ - see front matter q 2004 Elsevier B.V. All rights reserved. doi:10.1016/j.jhydrol.2003.10.023 * Corresponding author. Tel.: þ1-405-325-5928; fax: þ 1-405- 325-7477. E-mail addresses: cshy@public.wh.hb.cn (J. Yang); donzhang@ou.edu (D. Zhang); zhiming@lanl.gov (Z. Lu).