Estimating geostatistical parameters and spatially-variable hydraulic conductivity within a catchment system using an ensemble smoother

Groundwater flow models are important tools in assessing baseline conditions and investigating management alternatives in groundwater systems. The usefulness of these models, however, is often hindered by insufficient knowledge regarding the magnitude and spatial distribution of the spatially-distributed parameters, such as hydraulic conductivity ( K ), that govern the response of these models. Proposed parameter estimation methods frequently are demonstrated using simplified aquifer representations, when in reality the groundwater regime in a given watershed is influenced by strongly-coupled surface-subsurface processes. Furthermore, parameter estimation methodologies that rely on a geostatistical structure of K often assume the parameter values of the geostatistical model as known or estimate these values from limited data. In this study, we investigate the use of a data assimilation algorithm, the Ensemble Smoother, to provide enhanced estimates of K within a catchment system using the fully-coupled, surface-subsurface flow model CATHY. Both water table elevation and streamflow data are assimilated to condition the spatial distribution of K . An iterative procedure using the ES update routine, in which geostatistical parameter values defining the true spatial structure of K are identified, is also presented. In this procedure, parameter values are inferred from the updated ensemble of K fields and used in the subsequent iteration to generate the K ensemble, with the process proceeding until parameter values are converged upon. The parameter estimation scheme is demonstrated via a synthetic three-dimensional tilted v-shaped catchment system incorporating stream flow and variably-saturated subsurface flow, with spatio-temporal variability in forcing terms. Results indicate that the method is successful in providing improved estimates of the K field, and that the iterative scheme can be used to identify the geostatistical parameter values of the aquifer system. In general, water table data have a much greater ability than streamflow data to condition K . Future research includes applying the methodology to an actual regional study site.

[1]  D. R. Nielsen,et al.  On describing and predicting the hydraulic properties of unsaturated soils , 1985 .

[2]  G. Evensen,et al.  Data assimilation and inverse methods in terms of a probabilistic formulation , 1996 .

[3]  G. Marsily,et al.  Cokriged estimation of aquifer transmissivity as an indirect solution of the inverse problem: A practical approach , 1993 .

[4]  P. Kitanidis.  On the geostatistical approach to the inverse problem , 1996 .

[5]  Johan Valstar,et al.  A representer‐based inverse method for groundwater flow and transport applications , 2004 .

[6]  Dongxiao Zhang,et al.  Investigation of flow and transport processes at the MADE site using ensemble Kalman filter , 2008 .

[7]  Andrew F. Bennett,et al.  Inverse Methods in Physical Oceanography: Frontmatter , 1992 .

[8]  Mario Putti,et al.  Coupled surface runoff and subsurface flow model for catchment simulations , 1999 .

[9]  Yann Kerr,et al.  Deriving catchment-scale water and energy balance parameters using data assimilation based on extended Kalman filtering , 2002 .

[10]  Cokriging estimation of the conductivity field under variably saturated flow conditions , 1999 .

[11]  M. Rader,et al.  A computer model , 1996, IEEE Conference Record - Abstracts. 1996 IEEE International Conference on Plasma Science.

[12]  Stefano Orlandini,et al.  Physically-based distributed model for coupled surface runoff and subsurface flow simulation at the catchment scale , 2000 .

[13]  Renzo Rosso,et al.  Diffusion Wave Modeling of Distributed Catchment Dynamics , 1996 .

[14]  E. G. Vomvoris,et al.  A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one‐dimensional simulations , 1983 .

[15]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

[16]  Dara Entekhabi,et al.  An ensemble‐based reanalysis approach to land data assimilation , 2005 .

[17]  H. Ngodock,et al.  The Representer Method, the Ensemble Kalman Filter and the Ensemble Kalman Smoother: A Comparison Study Using a Nonlinear Reduced Gravity Ocean Model , 2006 .

[18]  W. Kinzelbach,et al.  Real‐time groundwater flow modeling with the Ensemble Kalman Filter: Joint estimation of states and parameters and the filter inbreeding problem , 2008 .

[19]  D. Aubert,et al.  Sequential assimilation of soil moisture and streamflow data in a conceptual rainfall-runoff model , 2003 .

[20]  G. Evensen,et al.  An ensemble Kalman smoother for nonlinear dynamics , 2000 .

[21]  G. Evensen.  Data Assimilation: The Ensemble Kalman Filter , 2006 .

[22]  Domenico A Baú,et al.  Optimal design of pump-and-treat systems under uncertain hydraulic conductivity and plume distribution. , 2008, Journal of contaminant hydrology.

[23]  Andrew F. Bennett,et al.  Inverse Methods in Physical Oceanography: Bibliography , 1992 .

[24]  T.-C. Jim Yeh,et al.  A Geostatistical Inverse Method for Variably Saturated Flow in the Vadose Zone , 1995 .

[25]  Anthony J. Jakeman,et al.  An application of the Kalman filtering technique for Streamflow forecasting in the Upper Murray Basin , 2001 .

[26]  Christian L. Keppenne,et al.  Data Assimilation into a Primitive-Equation Model with a Parallel Ensemble Kalman Filter , 2000 .

[27]  Y. Gautier.  Geostatistical Parameters Estimation Using Well Test Data , 2004 .

[28]  T.-C. Jim Yeh,et al.  An iterative geostatistical inverse method for steady flow in the vadose zone , 1996 .

[29]  J. Gómez-Hernández,et al.  Uncertainty assessment and data worth in groundwater flow and mass transport modeling using a blocking Markov chain Monte Carlo method. , 2009 .

[30]  Andrés Sahuquillo,et al.  Joint simulation of transmissivity and storativity fields conditional to steady-state and transient hydraulic head data , 1999 .

[31]  Jianlin Fu,et al.  Uncertainty assessment and data worth in groundwater ∞ow and mass transport modeling using a blocking Markov chain , 2009 .

[32]  Xianhong Xie,et al.  Data assimilation for distributed hydrological catchment modeling via ensemble Kalman filter , 2010 .

[33]  Steven A. Margulis,et al.  A Bayesian approach to snow water equivalent reconstruction , 2008 .

[34]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[35]  M. G. Marietta,et al.  Pilot Point Methodology for Automated Calibration of an Ensemble of Conditionally Simulated Transmissivity Fields: 2. Application , 1995 .

[36]  Andrés Sahuquillo,et al.  Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data , 2003 .

[37]  Allan D. Woodbury,et al.  Simultaneous inversion of hydrogeologic and thermal data: 2. Incorporation of thermal data , 1988 .

[38]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .

[39]  Y. Gautier,et al.  Determination of Geostatistical Parameters Using Well Test Data , 2004 .

[40]  Yan Chen,et al.  Data assimilation for transient flow in geologic formations via ensemble Kalman filter , 2006 .

[41]  Mario Putti,et al.  Ensemble Kalman filter data assimilation for a process‐based catchment scale model of surface and subsurface flow , 2009 .

[42]  J. Gómez-Hernández,et al.  Stochastic conditional inverse modeling of subsurface mass transport: A brief review and the self-calibrating method , 2003 .

[43]  C. Paniconi,et al.  Surface‐subsurface flow modeling with path‐based runoff routing, boundary condition‐based coupling, and assimilation of multisource observation data , 2010 .

[44]  Steven M. Gorelick,et al.  Coupled process parameter estimation and prediction uncertainty using hydraulic head and concentration data , 1991 .

[45]  G. Evensen.  Using the Extended Kalman Filter with a Multilayer Quasi-Geostrophic Ocean Model , 1992 .

[46]  Mohammadali Tarrahi,et al.  Assessing the performance of the ensemble Kalman filter for subsurface flow data integration under variogram uncertainty , 2011 .

[47]  R. Ibbitt,et al.  Hydrological data assimilation with the ensemble Kalman filter: Use of streamflow observations to update states in a distributed hydrological model , 2007 .

[48]  Domenico Baù,et al.  Ensemble smoother assimilation of hydraulic head and return flow data to estimate hydraulic conductivity distribution , 2010 .

[49]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

[50]  Henrik Madsen,et al.  Calibration framework for a Kalman filter applied to a groundwater model , 2006 .

[51]  A. S. Cullick,et al.  Construction of geostatistical aquifer models integrating dynamic flow and tracer data using inverse technique , 2002 .

[52]  Mario Putti,et al.  Time step and stability control for a coupled model of surface and subsurface flow , 2004 .

[53]  D. Baù,et al.  Estimating spatially-variable first-order rate constants in groundwater reactive transport systems. , 2011, Journal of contaminant hydrology.

[54]  Miguel A. Mariño,et al.  Estimation of Spatially Variable Aquifer Hydraulic Properties Using Kalman Filtering , 1997 .

[55]  D. McLaughlin,et al.  Hydrologic Data Assimilation with the Ensemble Kalman Filter , 2002 .

[56]  G. Evensen.  Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[57]  Peter K. Kitanidis,et al.  Analysis of the Spatial Structure of Properties of Selected Aquifers , 1985 .

[58]  P. Kitanidis,et al.  An Application of the Geostatistical Approach to the Inverse Problem in Two-Dimensional Groundwater Modeling , 1984 .

[59]  Harrie-Jan Hendricks Franssen,et al.  Ensemble Kalman filtering versus sequential self-calibration for inverse modelling of dynamic groundwater flow systems , 2009 .