Evaluation of confidence intervals for a steady-state leaky aquifer model

Abstract The fact that dependent variables of groundwater models are generally nonlinear functions of model parameters is shown to be a potentially significant factor in calculating accurate confidence intervals for both model parameters and functions of the parameters, such as the values of dependent variables calculated by the model. The Lagrangian method of Vecchia and Cooley [Vecchia, A.V. & Cooley, R.L., Water Resources Research, 1987, 23(7), 1237–1250] was used to calculate nonlinear Scheffe-type confidence intervals for the parameters and the simulated heads of a steady-state groundwater flow model covering 450 km2 of a leaky aquifer. The nonlinear confidence intervals are compared to corresponding linear intervals. As suggested by the significant nonlinearity of the regression model, linear confidence intervals are often not accurate. The commonly made assumption that widths of linear confidence intervals always underestimate the actual (nonlinear) widths was not correct. Results show that nonlinear effects can cause the nonlinear intervals to be asymmetric and either larger or smaller than the linear approximations. Prior information on transmissivities helps reduce the size of the confidence intervals, with the most notable effects occurring for the parameters on which there is prior information and for head values in parameter zones for which there is prior information on the parameters.

[1]  Franklin A. Graybill,et al.  Introduction to the Theory of Statistics, 3rd ed. , 1974 .

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

[3]  D. G. Watts,et al.  Relative Curvature Measures of Nonlinearity , 1980 .

[4]  M. Hill.  A computer program (MODFLOWP) for estimating parameters of a transient, three-dimensional ground-water flow model using nonlinear regression , 1992 .

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

[6]  Richard L. Cooley,et al.  Simultaneous confidence and prediction intervals for nonlinear regression models with application to a groundwater flow model , 1987 .

[7]  Keith Beven,et al.  The future of distributed models: model calibration and uncertainty prediction. , 1992 .

[8]  Douglas P. Wiens,et al.  Correction factors for F ratios in nonlinear regression , 1987 .

[9]  J. Wolfowitz,et al.  Introduction to the Theory of Statistics. , 1951 .

[10]  Irwin Guttman,et al.  On Beale's Measures of Non-linearity , 1965 .

[11]  Steen Christensen,et al.  On the Strategy of Estimating Regional‐Scale Transmissivity Fields , 1997 .

[12]  G. P. Clarke,et al.  Approximate Confidence Limits for a Parameter Function in Nonlinear Regression , 1987 .

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

[14]  H. Linssen,et al.  Non-linearity measures : a case study , 1975 .

[15]  E. Beale,et al.  Confidence Regions in Non‐Linear Estimation , 1960 .

[16]  N. Draper,et al.  Applied Regression Analysis , 1966 .

[17]  Richard L. Cooley,et al.  Exact Scheffé-type confidence intervals for output from groundwater flow models. 1. Use of hydrogeologic information , 1993 .

[18]  S. Christensen.  Hydrological Model for the Tude Å Catchment , 1994 .

[19]  Roger J. Brooks,et al.  Determining the range of predictions of a groundwater model which arises from alternative calibrations , 1994 .

[20]  Norman R. Draper,et al.  Applied regression analysis (2. ed.) , 1981, Wiley series in probability and mathematical statistics.