Conditional first‐order second‐moment method and its application to the quantification of uncertainty in groundwater modeling

[1] Decision making in water resources management usually requires the quantification of uncertainties. Monte Carlo techniques are suited for this analysis but imply a huge computational effort. An alternative and computationally efficient approach is the first-order second-moment (FOSM) method which directly propagates parameter uncertainty into the result. We apply the FOSM method to both the groundwater flow and solute transport equations. It is shown how conditioning on the basis of measured heads and/or concentrations yields the “principle of interdependent uncertainty” that correlates the uncertainties of feasible hydraulic conductivities and recharge rates. The method is used to compute the uncertainty of steady state heads and of steady state solute concentrations. It is illustrated by an application to the Palla Road Aquifer in semiarid Botswana, for which the quantification of the uncertainty range of groundwater recharge is of prime interest. The uncertainty bounds obtained by the FOSM method correspond well with the results obtained by the Monte Carlo method. The FOSM method, however, is much more advantageous with respect to computational efficiency. It is shown that at the planned abstraction rate the probability of exceeding the natural replenishment of the Palla Road Aquifer by overpumping is 30%.

[1]  J. P. Delhomme,et al.  Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach , 1979 .

[2]  George F. Pinder,et al.  Simulation of groundwater flow and mass transport under uncertainty , 1977 .

[3]  Harald Kunstmann,et al.  Groundwater Resources Management UnderUncertainty , 1998 .

[4]  E. Eriksson,et al.  Chloride concentration in groundwater, recharge rate and rate of deposition of chloride in the Israel Coastal Plain , 1969 .

[5]  Gedeon Dagan,et al.  Analysis of flow through heterogeneous random aquifers: 2. Unsteady flow in confined formations , 1982 .

[6]  W. Kinzelbach,et al.  Computation of stochastic wellhead protection zones by combining the first-order second-moment method and Kolmogorov backward equation analysis , 2000 .

[7]  M. S. Bedinger Estimation of Natural Groundwater Recharge , 1989 .

[8]  W. Durner Groundwater recharge: A guide to understanding and estimating natural recharge , 1992 .

[9]  W. W. Wood,et al.  Chemical and Isotopic Methods for Quantifying Ground‐Water Recharge in a Regional, Semiarid Environment , 1995 .

[10]  A. Kiureghian,et al.  First‐order reliability approach to stochastic analysis of subsurface flow and contaminant transport , 1987 .

[11]  Allan L. Gutjahr,et al.  Cross‐correlated random field generation with the direct Fourier Transform Method , 1993 .

[12]  Dennis McLaughlin,et al.  Stochastic analysis of nonstationary subsurface solute transport: 1. Unconditional moments , 1989 .

[13]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[14]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

[15]  J. L. Devary,et al.  Pore velocity estimation uncertainties , 1982 .

[16]  N. Cox Statistical Models in Engineering , 1970 .

[17]  W. Yeh Review of Parameter Identification Procedures in Groundwater Hydrology: The Inverse Problem , 1986 .

[18]  W. Edmunds,et al.  A geochemical and isotopic approach to recharge evaluation in semi-arid zones; past and present , 1980 .

[19]  A. Scheidegger General Theory of Dispersion in Porous Media , 1961 .

[20]  E. Wood,et al.  A distributed parameter approach for evaluating the accuracy of groundwater model predictions: 1. Theory , 1988 .

[21]  L. Gelhar Stochastic Subsurface Hydrology , 1992 .

[22]  S. P. Neuman,et al.  Effects of kriging and inverse modeling on conditional simulation of the Avra Valley Aquifer in southern Arizona , 1982 .

[23]  Adrian E. Scheidegger,et al.  Statistical Hydrodynamics in Porous Media , 1954 .

[24]  Richard L. Cooley,et al.  A method of estimating parameters and assessing reliability for models of steady state Groundwater flow: 2. Application of statistical analysis , 1979 .

[25]  C. Tiedeman,et al.  Analysis of uncertainty in optimal groundwater contaminant capture design , 1993 .

[26]  Wilson H. Tang,et al.  Probability concepts in engineering planning and design , 1984 .

[27]  L. Townley,et al.  Computationally Efficient Algorithms for Parameter Estimation and Uncertainty Propagation in Numerical Models of Groundwater Flow , 1985 .

[28]  Curtis M. Oldenburg,et al.  Linear and Monte Carlo uncertainty analysis for subsurface contaminant transport simulation , 1997 .

[29]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[30]  Srikanta Mishra,et al.  Probabilistic Multiphase Flow Modeling Using the Limit‐State Method , 1997 .

[31]  M. Hughes,et al.  The use of environmental chloride and tritium to estimate total recharge to an unconfined aquifer , 1978 .

[32]  Joel P. Conte,et al.  Probabilistic Screening Tool for Ground-Water Contamination Assessment , 1995 .

[33]  J. Bear Dynamics of Fluids in Porous Media , 1975 .

[34]  Peter K. Kitanidis,et al.  Comparison of Gaussian Conditional Mean and Kriging Estimation in the Geostatistical Solution of the Inverse Problem , 1985 .

[35]  L. Connell An analysis of perturbation based methods for the treatment of parameter uncertainty in numerical groundwater models , 1995 .

[36]  Michael D. Dettinger,et al.  First order analysis of uncertainty in numerical models of groundwater flow part: 1. Mathematical development , 1981 .

[37]  M. Sharma,et al.  Groundwater recharge estimation using chloride, deuterium and oxygen-18 profiles in the deep coastal sands of Western Australia , 1985 .

[38]  R. Allan Freeze,et al.  Stochastic analysis of steady state groundwater flow in a bounded domain: 2. Two‐dimensional simulations , 1979 .

[39]  D. McLaughlin,et al.  Stochastic analysis of nonstationary subsurface solute transport: 2. Conditional moments , 1989 .

[40]  D. A. Barry,et al.  The first‐order reliability method of predicting cumulative mass flux in heterogeneous porous formations , 1997 .