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

VL - 14

SP - 252

EP - 272

JO - Advances in Water Resources

JF - Advances in Water Resources

SN - 0309-1708

IS - 5

ER -