Abstract
Comprehensive information about the spatial distribution of the subsurface hydraulic properties is crucial to model groundwater flow, to predict solute transport in aquifers and to design remediation actions. In this work, a Bayesian Geostatistical approach, as implemented in bgaPEST, was adopted to estimate the hydraulic properties of a well field located at the Campus of Science and Technology of the University of Parma (Northern Italy), in a contest of a highly parameterized inversion. Head data, collected by means of multi frequency oscillatory pumping tests, were used to both estimate the hydraulic parameters and validate the results. The groundwater flow processes were modelled by means of MODFLOW 2005 and an adjoint-state formulation of the same software was used to efficiently calculate the sensitivity matrix, required by the inverse procedure. The Bayesian Geostatistical approach estimated the hydraulic conductivity and specific storage fields, handling a large number of parameters. The results of the inversion are consistent with the alluvial nature of the investigated aquifer and the preliminary traditional pumping tests carried out at the site.
Similar content being viewed by others
References
Alberti L, Lombi S, Zanini A (2011) Identifying sources of chlorinated aliphatic hydrocarbons in a residential area in Italy using the integral pumping test method. Hydrogeol J 19(6):1253–1267. https://doi.org/10.1007/s10040-011-0742-1
Butera I, Tanda MG, Zanini A (2013) Simultaneous identification of the pollutant release history and the source location in groundwater by means of a geostatistical approach. Stoch Environ Res Risk Assess 27(5):1269–1280. https://doi.org/10.1007/s00477-012-0662-1
Cardiff M, Bakhos T, Kitanidis PK, Barrash W (2013) Aquifer heterogeneity characterization with oscillatory pumping: sensitivity analysis and imaging potential. Water Resour Res 49:5395–5410. https://doi.org/10.1002/wrcr.20356
Cardiff M, Barrash W (2011) 3-D transient hydraulic tomography in unconfined aquifers with fast drainage response. Water Resour Res 47:W12518. https://doi.org/10.1029/2010WR010367
Cardiff M, Barrash W, Kitanidis PK (2012) A field proof-of-concept of aquifer imaging using 3-D transient hydraulic tomography with modular, temporarily emplaced equipment. Water Resour Res 48:W05531. https://doi.org/10.1029/2011WR011704
Carrera J, Alcolea A, Medina A, Hidalgo J, Slootenet LJ (2005) Inverse problem in hydrogeology. Hydrogeol J 13:206. https://doi.org/10.1007/s10040-004-0404-7
Chen Y, Oliver DS (2012) Ensemble randomized maximum likelihood method as an iterative ensemble smoother. Math Geosci 44:1–26. https://doi.org/10.1007/s11004-011-9376-z
Clemo T (2007) MODFLOW-2005 ground water model—user guide to the adjoint state based sensitivity process (ADJ). Center for the Geophysical Investigation of the Shallow Subsurface, Boise State University, Boise
Cupola F, Tanda MG, Zanini A (2015) Contaminant release history identification in 2-D heterogeneous aquifers through a minimum relative entropy approach. SpringerPlus 4:656. https://doi.org/10.1186/s40064-015-1465-x
Doherty JE (2010) PEST, model-independent parameter estimation—user manual, 5th edn. Watermark Numerical Computing, Brisbane (with slight additions)
Doherty JE, Hunt RJ (2010) Approaches to highly parameterized inversion—a guide to using PEST for groundwater-model calibration: U.S. Geological Survey Scientific Investigations Report 5169. https://pubs.usgs.gov/sir/2010/5169
D’Oria M, Fienen MN (2012) MODFLOW-style parameters in underdetermined parameter estimation. Ground Water 50:149–153. https://doi.org/10.1111/j.1745-6584.2011.00803.x
D’Oria M, Mignosa P, Tanda MG (2014) Bayesian estimation of inflow hydrographs in ungauged sites of multiple reach systems. Adv Water Resour 63:143–151. https://doi.org/10.1016/j.advwatres.2013.11.007
D’Oria M, Mignosa P, Tanda MG (2015) An inverse method to estimate the flow through a levee breach. Adv Water Resour 82:166–175. https://doi.org/10.1016/j.advwatres.2015.05.002
D’Oria M, Tanda MG (2012) Reverse flow routing in open channels: a Bayesian geostatistical approach. J Hydrol 460:130–135. https://doi.org/10.1016/j.jhydrol.2012.06.055
Emerick AA, Reynolds AC (2013) Ensemble smoother with multiple data assimilation. Comput Geosci 55:3–15. https://doi.org/10.1016/j.cageo.2012.03.011
Evensen G (2009) Data assimilation: the ensemble Kalman filter, 2nd edn. Springer, Berlin. https://doi.org/10.1007/978-3-642-03711-5
Fienen MN, Clemo TM, Kitanidis PK (2008) An interactive Bayesian geostatistical inverse protocol for hydraulic tomography. Water Resour Res 44:W00B01. https://doi.org/10.1029/2007WR006730
Fienen M, Hunt R, Krabbenhoft D, Clemo T (2009) Obtaining parsimonious hydraulic conductivity fields using head and transport observations: a Bayesian geostatistical parameter estimation approach. Water Resour Res 45:W08405. https://doi.org/10.1029/2008wr007431
Fienen MN, D’Oria M, Doherty JE, Hunt RJ (2013) Approaches in highly parameterized inversion: bgaPEST, a Bayesian geostatistical approach implementation with PEST-Documentation and instructions. In: U.S. Geological Survey Techniques and Methods 7-C9. https://pubs.usgs.gov/tm/07/c09
Grana D, Fjeldstad T, Omre H (2017) Bayesian Gaussian mixture linear inversion for geophysical inverse problems. Math Geosci 49:493. https://doi.org/10.1007/s11004-016-9671-9
Hansen TM, Journel AG, Tarantola A, Mosegaard K (2006) Linear inverse Gaussian theory and geostatistics. Geophysics 71(6):R101–R111. https://doi.org/10.1190/1.2345195
Hantush MS (1956) Analysis of data from pumping tests in leaky aquifers. Trans Am Geophys Union 37(6):702–714
Harbaugh AW (2005) MODFLOW-2005, The U.S. Geological Survey modular ground-water model—the Ground-Water Flow Process. In: U.S. Geological Survey Techniques and Methods 6-A16. https://pubs.usgs.gov/tm/2005/tm6A16
Hendricks Franssen HJ, Alcolea A, Riva M, Bakr M, van der Wiel N, Stauffer F, Guadagnini A (2009) A comparison of seven methods for the inverse modelling of groundwater flow. Application to the characterisation of well catchments. Adv Water Resour 32(6):851–872. https://doi.org/10.1016/j.advwatres.2009.02.011
Hill MC (2006) Using models to manage systems subject to sustainability indicators. In: Webb B, Hirata R, Kruse E, Vrba J (eds) Sustainability of groundwater resources and its indicators, IAHS publication 302. IAHS Press, Wallingford
Hoeksema RJ, Kitanidis PK (1984) An application of the geostatistical approach to the inverse problem in two-dimensional groundwater modeling. Water Resour Res 20(7):1003–1020. https://doi.org/10.1029/WR020i007p01003
Kitanidis PK, Vomvoris EG (1983) A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one-dimensional simulations. Water Resour Res 19(3):677–690. https://doi.org/10.1029/WR019i003p00677
Kitanidis PK (1993) Generalized covariance functions in estimation. Math Geol 25(5):525–540. https://doi.org/10.1007/BF00890244
Kitanidis PK (1995) Quasi-linear geostatistical theory for inversing. Water Resour Res 31(10):2411–2419. https://doi.org/10.1029/95WR01945
Kitanidis PK (1997) Introduction to geostatistics, applications in hydrogeology. Cambridge University Press, Cambridge
Kitanidis PK, Lee J (2014) Principal component geostatistical approach for large-dimensional inverse problems. Water Resour Res 50:5428–5443. https://doi.org/10.1002/2013WR014630
Leonhardt G, D’Oria M, Kleidorfer M, Rauch W (2014) Estimating inflow to a combined sewer overflow structure with storage tank in real time: evaluation of different approaches. Water Sci Technol 70(7):1143–1151. https://doi.org/10.2166/wst.2014.331
Li L, Srinivasan S, Zhou H, Gomez-Hernandez JJ (2015) Two-point or multiple point statistics? A comparison between the ensemble Kalman filtering and the ensemble pattern matching inverse methods. Adv Water Resour 86(B):297–310. https://doi.org/10.1016/j.advwatres.2015.05.014
Liu X, Lee J, Kitanidis PK, Parker J, Kim U (2012) Value of information as a context-specific measure of uncertainty in groundwater remediation. Water Resour Manage 26(6):1513–1535. https://doi.org/10.1007/s11269-011-9970-3
McLaughlin D, Townley LR (1996) A reassessment of the groundwater inverse problem. Water Resour Res 32(5):1131–1161. https://doi.org/10.1029/96WR00160
Michalak AM, Bruhwiler L, Tans PP (2004) A geostatistical approach to surface flux estimation of atmospheric trace gases. J Geophys Res 109:D14. https://doi.org/10.1029/2003jd004422
Nowak W, Cirpka OA (2004) A modified Levenberg–Marquardt algorithm for quasi-linear geostatistical inversing. Adv Water Resour 27(7):737–750. https://doi.org/10.1016/j.advwatres.2004.03.004
Rabinovich A, Barrash W, Cardiff M, Hochstetler DL, Bakhos T, Dagan G, Kitanidis PK (2015) Frequency dependent hydraulic properties estimated from oscillatory pumping tests in an unconfined aquifer. J Hydrol 531(1):2–16. https://doi.org/10.1016/j.jhydrol.2015.08.021
Regione Emilia-Romagna, Servizio Geologico Sismico e dei Suoli, ENI-AGIP (1998) Riserve idriche sotterranee della Regione Emilia-Romagna, p 119 (in Italian)
Rimstad K, Avseth P, Omre H (2012) Hierarchical Bayesian lithology/fluid prediction: a North Sea case study. Geophysics 77(2):B69–B85
Rubin Y, Chen X, Murakami H, Hahn M (2010) A Bayesian approach for inverse modeling, data assimilation, and conditional simulation of spatial random fields. Water Resour Res 46:W10523. https://doi.org/10.1029/2009WR008799
Snodgrass MF, Kitanidis PK (1997) A geostatistical approach to contaminant source identification. Water Resour Res 33(4):537–546. https://doi.org/10.1029/96WR03753
Soares A, Nunes R, Azevedo L (2017) Integration of uncertain data in geostatistical modelling. Math Geosci 49:253. https://doi.org/10.1007/s11004-016-9667-5
Tarantola A (2005) Inverse problem theory and methods for model parameter estimation. SIAM, Philadelphia
Theis CV (1935) The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using groundwater storage. Trans Am Geophys Union 2:519–524
Yeh WWG (1986) Review of parameter identification procedures in groundwater hydrology: the inverse problem. Water Resour Res 22(2):95–108. https://doi.org/10.1029/WR022i002p00095
Yeh TCJ, Liu S (2000) Hydraulic tomography: development of a new aquifer test method. Water Resour Res 36(8):2095–2105. https://doi.org/10.1029/2000WR900114
Veling EJM, Maas C (2010) Hantush well function revisited. J Hydrol 393(3):381–388. https://doi.org/10.1016/j.jhydrol.2010.08.033
Xu T, Gomez-Hernandez JJ (2016) Joint identification of contaminant source location, initial release time, and initial solute concentration in an aquifer via ensemble Kalman filtering. Water Resour Res 52:6587–6595. https://doi.org/10.1002/2016WR019111
Zanini A, Kitanidis PK (2009) Geostatistical inversing for large-contrast transmissivity fields. Stoch Environ Res Risk Assess 23:565. https://doi.org/10.1007/s00477-008-0241-7
Zanini A, Woodbury AD (2016) Contaminant source reconstruction by empirical Bayes and Akaike’s Bayesian information criterion. J Contam Hydrol 185–186:74–86. https://doi.org/10.1016/j.jconhyd.2016.01.006
Zanini A, Tanda MG, Woodbury AD (2017) Identification of transmissivity fields using a Bayesian strategy and perturbative approach. Adv Water Resour 108:69–82. https://doi.org/10.1016/j.advwatres.2017.07.022
Zhou H, Gómez-Hernández JJ, Liangping L (2014) Inverse methods in hydrogeology: evolution and recent trends. Adv Water Resour 63:22–37. https://doi.org/10.1016/j.advwatres.2013.10.014
Zhou H, Li L, Hendricks Franssen HJ, Gómez-Hernández JJ (2012) Pattern recognition in a bimodal aquifer using the normal-score ensemble Kalman filter. Math Geosci 44:169. https://doi.org/10.1007/s11004-011-9372
Zhou YQ, Lim D, Cupola F, Cardiff M (2016) Aquifer imaging with pressure waves—evaluation of low-impact characterization through sandbox experiments. Water Resour Res 52:2141–2156. https://doi.org/10.1002/2015WR017751
Zimmerman DA, de Marsily G, Gotway CA, Marietta MG, Axness CL, Beauheim RL, Bras RL, Carrera J, Dagan G, Davies PB, Gallegos DP, Galli A, Gómez-Hernández J, Grindrod P, Gutjahr AL, Kitanidis PK, Lavenue AM, McLaughlin D, Neuman SP, RamaRao BS, Ravenne C, Rubin Y (1998) A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow. Water Resour Res 34(6):1373–1413. https://doi.org/10.1029/98WR00003
Acknowledgements
The authors are grateful to the editor and the anonymous reviewers for their valuable comments and suggestions on this work.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Marco, D., Andrea, Z. & Fausto, C. Oscillatory Pumping Test to Estimate Aquifer Hydraulic Parameters in a Bayesian Geostatistical Framework. Math Geosci 50, 169–186 (2018). https://doi.org/10.1007/s11004-017-9717-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11004-017-9717-7