Contaminant source reconstruction by empirical Bayes and Akaike's Bayesian Information Criterion.

The objective of the paper is to present an empirical Bayesian method combined with Akaike's Bayesian Information Criterion (ABIC) to estimate the contaminant release history of a source in groundwater starting from few concentration measurements in space and/or in time. From the Bayesian point of view, the ABIC considers prior information on the unknown function, such as the prior distribution (assumed Gaussian) and the covariance function. The unknown statistical quantities, such as the noise variance and the covariance function parameters, are computed through the process; moreover the method quantifies also the estimation error through the confidence intervals. The methodology was successfully tested on three test cases: the classic Skaggs and Kabala release function, three sharp releases (both cases regard the transport in a one-dimensional homogenous medium) and data collected from laboratory equipment that consists of a two-dimensional homogeneous unconfined aquifer. The performances of the method were tested with two different covariance functions (Gaussian and exponential) and also with large measurement error. The obtained results were discussed and compared to the geostatistical approach of Kitanidis (1995).

[1]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[2]  Ritei Shibata.  Information Criteria for Statistical Model Selection , 2002 .

[3]  T. Ulrych,et al.  Minimum Relative Entropy Inversion: Theory and Application to Recovering the Release History of a Groundwater Contaminant , 1996 .

[4]  A. Zanini,et al.  Laboratory sandbox validation of pollutant source location methods , 2014, Stochastic Environmental Research and Risk Assessment.

[5]  Y. Mitsuhata.  Adjustment of regularization in ill‐posed linear inverse problems by the empirical Bayes approach , 2004 .

[6]  P. Kitanidis.  Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[7]  Alexander Y. Sun,et al.  A robust geostatistical approach to contaminant source identification , 2007 .

[8]  I. Butera,et al.  Use of numerical modelling to identify the transfer function and application to the geostatistical procedure in the solution of inverse problems in groundwater , 2006 .

[9]  D. Cox,et al.  An Analysis of Transformations , 1964 .

[10]  I. Butera,et al.  Use of numerical modelling to identify the transfer function and application to the geostatistical procedure in the solution of inverse problems in groundwater , 2006 .

[11]  I. Butera,et al.  A geostatistical approach to recover the release history of groundwater pollutants , 2003 .

[12]  Tad J. Ulrych,et al.  Information theory measures with application to model identification , 1986, IEEE Trans. Acoust. Speech Signal Process..

[13]  I. Butera,et al.  Simultaneous identification of the pollutant release history and the source location in groundwater by means of a geostatistical approach , 2013, Stochastic Environmental Research and Risk Assessment.

[14]  Kelin Wang,et al.  Minimum relative entropy: Theory and application to surface temperature reconstruction from borehole temperature measurements , 2000 .

[15]  M. Hendry,et al.  Clay Aquitards as Archives of Holocene Paleoclimate: δ18O and Thermal Profiling , 2007 .

[16]  Yefang Jiang,et al.  Transmissivity estimation for highly heterogeneous aquifers: comparison of three methods applied to the Edwards Aquifer, Texas, USA , 2007 .

[17]  Tadeusz J. Ulrych,et al.  A Bayes tour of inversion: A tutorial , 2001 .

[18]  H. Beltrami,et al.  A comparison of five different analyses in the interpretation of five borehole temperature data sets , 1992 .

[19]  A. Sun,et al.  A constrained robust least squares approach for contaminant release history identification , 2006 .

[20]  M. Hendry,et al.  Clay aquitards as archives of holocene paleoclimate: delta18O and thermal profiling. , 2007, Ground water.

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

[22]  A. Zanini,et al.  Evaluation of dispersivity coefficients by means of a laboratory image analysis. , 2015, Journal of contaminant hydrology.

[23]  P. Kitanidis,et al.  Geostatistical inversing for large-contrast transmissivity fields , 2009 .

[24]  Tadeusz J. Ulrych,et al.  Information-Based Inversion and Processing with Applications , 2011 .

[25]  T. Ulrych,et al.  A full‐Bayesian approach to the groundwater inverse problem for steady state flow , 2000 .

[26]  J. Bear,et al.  Modeling groundwater flow and pollution , 1987 .

[27]  A. Bagtzoglou,et al.  State of the Art Report on Mathematical Methods for Groundwater Pollution Source Identification , 2001 .

[28]  Allan D. Woodbury,et al.  Minimum relative entropy, Bayes and Kapur , 2011 .

[29]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[30]  G. Ferguson,et al.  Ground surface paleotemperature reconstruction using information measures and empirical Bayes , 2006 .

[31]  P. Smart,et al.  An evaluation of some fluorescent dyes for water tracing , 1977 .

[32]  Hirotugu Akaike,et al.  Likelihood and the Bayes procedure , 1980 .

[33]  H. Akaike.  A new look at the statistical model identification , 1974 .

[34]  P. Kitanidis,et al.  Estimation of historical groundwater contaminant distribution using the adjoint state method applied to geostatistical inverse modeling , 2004 .

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

[36]  R. Parker.  Understanding Inverse Theory , 1977 .

[37]  T. Skaggs,et al.  Recovering the release history of a groundwater contaminant , 1994 .

[38]  T. Ulrych,et al.  Minimum relative entropy and probabilistic inversion in groundwater hydrology , 1998 .

[39]  P. Kitanidis,et al.  A geostatistical approach to contaminant source identification , 1997 .