TY  - JOUR
AU  - Schoniger, A.
AU  - Nowak, W.
AU  - Hendricks-Franssen, H.J.
TI  - Parameter estimation by ensemble Kalman filters with transformed data: Approach and application to hydraulic tomography
JO  - Water resources research
VL  - 48
SN  - 0043-1397
CY  - Washington, DC
PB  - AGU
M1  - PreJuSER-21348
SP  - W04502
PY  - 2012
N1  - The authors would like to thank the German Research Foundation (DFG) for financial support of the project within the Cluster of Excellence in Simulation Technology (EXC 310/1) and within the International Research Training Group "Nonlinearities and upscaling in porous media" (NUPUS, IRTG 1398) at the University of Stuttgart.
AB  - Ensemble Kalman filters (EnKFs) are a successful tool for estimating state variables in atmospheric and oceanic sciences. Recent research has prepared the EnKF for parameter estimation in groundwater applications. EnKFs are optimal in the sense of Bayesian updating only if all involved variables are multivariate Gaussian. Subsurface flow and transport state variables, however, generally do not show Gaussian dependence on hydraulic log conductivity and among each other, even if log conductivity is multi-Gaussian. To improve EnKFs in this context, we apply nonlinear, monotonic transformations to the observed states, rendering them Gaussian (Gaussian anamorphosis, GA). Similar ideas have recently been presented by Beal et al. (2010) in the context of state estimation. Our work transfers and adapts this methodology to parameter estimation. Additionally, we address the treatment of measurement errors in the transformation and provide several multivariate analysis tools to evaluate the expected usefulness of GA beforehand. For illustration, we present a first-time application of an EnKF to parameter estimation from 3-D hydraulic tomography in multi-Gaussian log conductivity fields. Results show that (1) GA achieves an implicit pseudolinearization of drawdown data as a function of log conductivity and (2) this makes both parameter identification and prediction of flow and transport more accurate. Combining EnKFs with GA yields a computationally efficient tool for nonlinear inversion of data with improved accuracy. This is an attractive benefit, given that linearization-free methods such as particle filters are computationally extremely demanding.
KW  - J (WoSType)
LB  - PUB:(DE-HGF)16
UR  - <Go to ISI:>//WOS:000302531800001
DO  - DOI:10.1029/2011WR010462
UR  - https://juser.fz-juelich.de/record/21348
ER  -