Exact Scheffé-type confidence intervals for output from groundwater flow models. 1. Use of hydrogeologic information

A new method is developed to efficiently compute exact Scheffe-type confidence intervals for output (or other function of parameters) g([beta]) derived from a groundwater flow model. The method is general in that parameter uncertainty can be specified by any statistical distribution having a log probability density function (log pdf) that can be expanded in a Taylor series. However, for this study parameter uncertainty is specified by a statistical multivariate beta distribution that incorporates hydrogeologic information in the form of the investigator's best estimates of parameters and a grouping of random variables representing possible parameter values so that each group is defined by maximum and minimum bounds and an ordering according to increasing value. The new method forms the confidence intervals from maximum and minimum limits of g([beta]) on a contour of a linear combination of (1) the quadratic form for the parameters used by Cooley and Vecchia (1987) and (2) the log pdf for the multivariate beta distribution. Three example problems are used to compare characteristics of the confidence intervals for hydraulic head obtained using different weights for the linear combination. Different weights generally produced similar confidence intervals, whereas the method of Cooley and Vecchia (1987) often produced much largermore » confidence intervals.« less

[1]  R. L. Cooley.  Exact Scheffé-type confidence intervals for output from groundwater flow models: 2. Combined use of hydrogeologic information and calibration data , 1993 .

[2]  David Mautner Himmelblau,et al.  Applied Nonlinear Programming , 1972 .

[3]  Lynn W. Gelhar,et al.  Stochastic subsurface hydrology from theory to applications , 1986 .

[4]  Mary C. Hill,et al.  Analysis of accuracy of approximate, simultaneous, nonlinear confidence intervals on hydraulic heads in analytical and numerical test cases , 1989 .

[5]  Franklin A. Graybill,et al.  Theory and Application of the Linear Model , 1976 .

[6]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

[7]  Richard L. Cooley,et al.  A Comparison of Several Methods of Solving Nonlinear Regression Groundwater Flow Problems , 1985 .

[8]  Richard L. Cooley,et al.  CALCULATION OF NONLINEAR CONFIDENCE AND PREDICTION INTERVALS FOR GROUND-WATER FLOW MODELS. , 1987 .

[9]  L. Townley.  Second Order Effects of Uncertain Transmissivities on Predictions of Piezometric Heads , 1984 .

[10]  Lynn W. Gelhar,et al.  Stochastic Analysis of Flow in Heterogeneous Porous Media , 1984 .

[11]  Dan Rosbjerg,et al.  A Comparison of Four Inverse Approaches to Groundwater Flow and Transport Parameter Identification , 1991 .

[12]  Calyampudi Radhakrishna Rao,et al.  Linear Statistical Inference and its Applications , 1967 .

[13]  R. W. Andrews,et al.  Sensitivity Analysis for Steady State Groundwater Flow Using Adjoint Operators , 1985 .

[14]  Bryson C. Bates,et al.  Improved methodology for parameter inference in nonlinear, hydrologic regression models , 1992 .

[15]  Richard L. Cooley,et al.  Regression modeling of ground-water flow , 1990 .

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

[17]  Leonard F. Konikow,et al.  Predictive Accuracy of a Ground–Water Model — Lessons from a Postaudit , 1986 .

[18]  G. Dagan.  Statistical Theory of Groundwater Flow and Transport: Pore to Laboratory, Laboratory to Formation, and Formation to Regional Scale , 1986 .