Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information
- 1 February 1986
- journal article
- Published by American Geophysical Union (AGU) in Water Resources Research
- Vol. 22 (2) , 199-210
- https://doi.org/10.1029/wr022i002p00199
Abstract
In this series of three papers a method is presented to estimate the parameters of groundwater flow models under steady and nonsteady state conditions. The parameters include values and directions of principal hydraulic conductivities (or transmissivities) in anisotropic media, specific storage (or storativity), interior and boundary recharge or leakage rates, coefficients of head‐dependent interior and boundary sources, and boundary heads. In transient situations, the initial head distribution can also be estimated if the system is originally at a steady state. Paper 1 of the series discusses some of the advantage in treating the inverse problem statistically and in regularizing its solution by means of penalty criteria based on prior estimates of the parameters. The inverse problem is posed in the framework of maximum likelihood theory cast in a manner that accounts for prior information about the parameters. Since not all the factors which contribute to the prior errors can be quantified statistically at the outset, the covariance matrices of these errors are expressed in terms of several parameters which, if unknown, can be estimated jointly with the hydraulic parameters by a stagewise optimization process. When transient head data are separated by a fixed time interval, the temporal structure of these data is approximated by a lag‐one autoregressive model with a correlation coefficient that can be treated as another unknown parameter. Estimation errors are analyzed by examining the lower bound of their covariance matrix in the eigenspace. Paper 1 concludes by suggesting that certain model identification criteria developed in the time series literature, all of which are based on the maximum likelihood concept, might be useful for selecting the best groundwater model (or the best method of parameterizing a particular model) among a number of given alternatives.This publication has 52 references indexed in Scilit:
- Analysis of nonintrinsic spatial variability by residual kriging with application to regional groundwater levelsMathematical Geology, 1984
- A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one‐dimensional simulationsWater Resources Research, 1983
- Incorporation of prior information on parameters into nonlinear regression groundwater flow models: 2. ApplicationsWater Resources Research, 1983
- A statistical approach to the inverse problem of aquifer hydrology: 1. TheoryWater Resources Research, 1979
- A method of estimating parameters and assessing reliability for models of steady state Groundwater flow: 2. Application of statistical analysisWater Resources Research, 1979
- Estimating the Dimension of a ModelThe Annals of Statistics, 1978
- A method of estimating parameters and assessing reliability for models of steady state groundwater flow: 1. Theory and numerical propertiesWater Resources Research, 1977
- A new look at the statistical model identificationIEEE Transactions on Automatic Control, 1974
- Calibration of distributed parameter groundwater flow models viewed as a multiple‐objective decision process under uncertaintyWater Resources Research, 1973
- Use of Linear Programing for Estimating Geohydrologic Parameters of Groundwater BasinsWater Resources Research, 1971