TY - JOUR
T1 - Monte Carlo experiments with robust estimation of aquifer model parameters
AU - Wittmeyer, G.
AU - Neuman, S. P.
PY - 1991/10
Y1 - 1991/10
N2 - Automatic calibration procedures for fitting groundwater models to measured data have almost exclusively used L2 error norms for estimation criteria. Quadratic estimation criteria result from applying the method of maximum likelihood under the assumption of jointly normal errors in the measured head data which, for the case of uncorrelated errors with constant variance, results in the ordinary least squares criterion. Quadratic estimation criteria are extremely sensitive to errors which violate the Gaussian assumption and pose an additional source of instability for the inverse problem. We consider the use of statistically robust estimation procedures to reduce these destabilizing effects. In particular, we test the robust M-estimator of Huber with a series of Monte Carlo experiments, compute the sample mean squared error, squared bias and variance of the estimated log-hydraulic conductivities for realizations of head measurements contaminated with mixed normal noise and compare these measures to their counterparts for the OLS estimator. We show that Huber's M-estimator achieves the same estimation variance as the OLS estimator with reduced need for prior information to regularize the problem.
AB - Automatic calibration procedures for fitting groundwater models to measured data have almost exclusively used L2 error norms for estimation criteria. Quadratic estimation criteria result from applying the method of maximum likelihood under the assumption of jointly normal errors in the measured head data which, for the case of uncorrelated errors with constant variance, results in the ordinary least squares criterion. Quadratic estimation criteria are extremely sensitive to errors which violate the Gaussian assumption and pose an additional source of instability for the inverse problem. We consider the use of statistically robust estimation procedures to reduce these destabilizing effects. In particular, we test the robust M-estimator of Huber with a series of Monte Carlo experiments, compute the sample mean squared error, squared bias and variance of the estimated log-hydraulic conductivities for realizations of head measurements contaminated with mixed normal noise and compare these measures to their counterparts for the OLS estimator. We show that Huber's M-estimator achieves the same estimation variance as the OLS estimator with reduced need for prior information to regularize the problem.
UR - http://www.scopus.com/inward/record.url?scp=0026360458&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=0026360458&partnerID=8YFLogxK
U2 - 10.1016/0309-1708(91)90038-P
DO - 10.1016/0309-1708(91)90038-P
M3 - Article
AN - SCOPUS:0026360458
SN - 0309-1708
VL - 14
SP - 252
EP - 272
JO - Advances in Water Resources
JF - Advances in Water Resources
IS - 5
ER -