Data assimilation and parameter estimation via ensemble Kalman filter coupled with stochastic moment equations of transient groundwater flow

M. Panzeri, M. Riva, A. Guadagnini, Shlomo P Neuman

Research output: Contribution to journalArticle

30 Citations (Scopus)

Abstract

The ensemble Kalman filter (EnKF) is a powerful tool for assimilating data in earth system models. The approach allows real time Bayesian updating of system states and parameters as new data become available. This paper focuses on EnKF data assimilation in models of groundwater flow through complex geologic media. It has become common to treat the hydraulic conductivity of such media as correlated random fields conditioned on measured conductivity (medium property) and/or hydraulic head (system state) values. This renders the conductivity nonstationary and the corresponding conditional flow equations stochastic. Solving these equations and coupling them with EnKF generally entails computationally intensive Monte Carlo (MC) simulation. We propose to circumvent the need for MC through a direct solution of approximate nonlocal (integrodifferential) equations that govern the space-time evolution of conditional ensemble means (statistical expectations) and covariances of hydraulic heads and fluxes. We illustrate and explore our approach on synthetic two-dimensional examples in which a well pumps water from a randomly heterogeneous aquifer subject to prescribed head and flux boundary conditions. Embedding the solution in EnKF provides sequential updates of conductivity and head estimates throughout the space-time domain of interest. We demonstrate the computational feasibility and accuracy of our methodology, showing that hydraulic conductivity estimates are more sensitive to early than to later head values and improve with increasing assimilation frequency at early time.

Original languageEnglish (US)
Pages (from-to)1334-1344
Number of pages11
JournalWater Resources Research
Volume49
Issue number3
DOIs
StatePublished - 2013

Fingerprint

transient flow
Kalman filter
data assimilation
groundwater flow
conductivity
hydraulic head
hydraulic conductivity
boundary condition
aquifer
methodology
parameter estimation
simulation

Keywords

  • conditioning
  • data assimilation
  • ensemble Kalman filter
  • moment equations
  • stochastic models
  • transient groundwater flow

ASJC Scopus subject areas

  • Water Science and Technology

Cite this

Data assimilation and parameter estimation via ensemble Kalman filter coupled with stochastic moment equations of transient groundwater flow. / Panzeri, M.; Riva, M.; Guadagnini, A.; Neuman, Shlomo P.

In: Water Resources Research, Vol. 49, No. 3, 2013, p. 1334-1344.

Research output: Contribution to journalArticle

@article{0e0d9b03f76440aab6c9b2e559347960,
title = "Data assimilation and parameter estimation via ensemble Kalman filter coupled with stochastic moment equations of transient groundwater flow",
abstract = "The ensemble Kalman filter (EnKF) is a powerful tool for assimilating data in earth system models. The approach allows real time Bayesian updating of system states and parameters as new data become available. This paper focuses on EnKF data assimilation in models of groundwater flow through complex geologic media. It has become common to treat the hydraulic conductivity of such media as correlated random fields conditioned on measured conductivity (medium property) and/or hydraulic head (system state) values. This renders the conductivity nonstationary and the corresponding conditional flow equations stochastic. Solving these equations and coupling them with EnKF generally entails computationally intensive Monte Carlo (MC) simulation. We propose to circumvent the need for MC through a direct solution of approximate nonlocal (integrodifferential) equations that govern the space-time evolution of conditional ensemble means (statistical expectations) and covariances of hydraulic heads and fluxes. We illustrate and explore our approach on synthetic two-dimensional examples in which a well pumps water from a randomly heterogeneous aquifer subject to prescribed head and flux boundary conditions. Embedding the solution in EnKF provides sequential updates of conductivity and head estimates throughout the space-time domain of interest. We demonstrate the computational feasibility and accuracy of our methodology, showing that hydraulic conductivity estimates are more sensitive to early than to later head values and improve with increasing assimilation frequency at early time.",
keywords = "conditioning, data assimilation, ensemble Kalman filter, moment equations, stochastic models, transient groundwater flow",
author = "M. Panzeri and M. Riva and A. Guadagnini and Neuman, {Shlomo P}",
year = "2013",
doi = "10.1002/wrcr.20113",
language = "English (US)",
volume = "49",
pages = "1334--1344",
journal = "Water Resources Research",
issn = "0043-1397",
publisher = "American Geophysical Union",
number = "3",

}

TY - JOUR

T1 - Data assimilation and parameter estimation via ensemble Kalman filter coupled with stochastic moment equations of transient groundwater flow

AU - Panzeri, M.

AU - Riva, M.

AU - Guadagnini, A.

AU - Neuman, Shlomo P

PY - 2013

Y1 - 2013

N2 - The ensemble Kalman filter (EnKF) is a powerful tool for assimilating data in earth system models. The approach allows real time Bayesian updating of system states and parameters as new data become available. This paper focuses on EnKF data assimilation in models of groundwater flow through complex geologic media. It has become common to treat the hydraulic conductivity of such media as correlated random fields conditioned on measured conductivity (medium property) and/or hydraulic head (system state) values. This renders the conductivity nonstationary and the corresponding conditional flow equations stochastic. Solving these equations and coupling them with EnKF generally entails computationally intensive Monte Carlo (MC) simulation. We propose to circumvent the need for MC through a direct solution of approximate nonlocal (integrodifferential) equations that govern the space-time evolution of conditional ensemble means (statistical expectations) and covariances of hydraulic heads and fluxes. We illustrate and explore our approach on synthetic two-dimensional examples in which a well pumps water from a randomly heterogeneous aquifer subject to prescribed head and flux boundary conditions. Embedding the solution in EnKF provides sequential updates of conductivity and head estimates throughout the space-time domain of interest. We demonstrate the computational feasibility and accuracy of our methodology, showing that hydraulic conductivity estimates are more sensitive to early than to later head values and improve with increasing assimilation frequency at early time.

AB - The ensemble Kalman filter (EnKF) is a powerful tool for assimilating data in earth system models. The approach allows real time Bayesian updating of system states and parameters as new data become available. This paper focuses on EnKF data assimilation in models of groundwater flow through complex geologic media. It has become common to treat the hydraulic conductivity of such media as correlated random fields conditioned on measured conductivity (medium property) and/or hydraulic head (system state) values. This renders the conductivity nonstationary and the corresponding conditional flow equations stochastic. Solving these equations and coupling them with EnKF generally entails computationally intensive Monte Carlo (MC) simulation. We propose to circumvent the need for MC through a direct solution of approximate nonlocal (integrodifferential) equations that govern the space-time evolution of conditional ensemble means (statistical expectations) and covariances of hydraulic heads and fluxes. We illustrate and explore our approach on synthetic two-dimensional examples in which a well pumps water from a randomly heterogeneous aquifer subject to prescribed head and flux boundary conditions. Embedding the solution in EnKF provides sequential updates of conductivity and head estimates throughout the space-time domain of interest. We demonstrate the computational feasibility and accuracy of our methodology, showing that hydraulic conductivity estimates are more sensitive to early than to later head values and improve with increasing assimilation frequency at early time.

KW - conditioning

KW - data assimilation

KW - ensemble Kalman filter

KW - moment equations

KW - stochastic models

KW - transient groundwater flow

UR - http://www.scopus.com/inward/record.url?scp=84876720538&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84876720538&partnerID=8YFLogxK

U2 - 10.1002/wrcr.20113

DO - 10.1002/wrcr.20113

M3 - Article

AN - SCOPUS:84876720538

VL - 49

SP - 1334

EP - 1344

JO - Water Resources Research

JF - Water Resources Research

SN - 0043-1397

IS - 3

ER -