UC Santa Barbara

Joint estimation of transmissivity (T) and storativity (S) in a confined aquifer is done via maximum likelihood (ML). The differential equation of groundwater flow is discretized by the finite-element method, leading to equation ψφ{symbol}t+Γxt=ut. Elements of matrices ψ and Γ, as well as estimated covariance matrix of noise term ut, are functions of T and S. By minimizing the negative loglikelihood function corresponding to discretized groundwater flow equation with respect to T and S, ML estimators are obtained. The ML approach is found to yield accurate estimates of T and S (within 9 and 10% of their actual values, respectively) and showed quadratic convergence in Newton's search technique. Prediction of aquifer response, using ML estimators, results in estimated piezometric heads accurate to ±0.5 m from their actual, exact values. Statistical properties of ML estimators are derived and some basic results for statistical inference are given. © 1986 Plenum Publishing Corporation.