Characterization of groundwater and surface water mixing in a semiconfined karst aquifer using time‐lapse electrical resistivity tomography

Groundwater flow in karst includes exchange of water between large fractures, conduits, and the surrounding porous matrix, which impacts both water quality and quantity. Electrical resistivity tomography combined with end-member mixing analysis (EMMA) and numerical flow and transport modeling was used to study mixing of karst conduit and matrix waters to understand spatial and temporal patterns of mixing during high flow and base flow conditions. To our knowledge, this is the first time EMMA and synthetic geophysical simulations have been combined. Here we interpret an 8 week time-lapse electrical resistivity data set to assess groundwater-surface mixing. We simulate flow between the karst conduits and the porous matrix to determine fractions of water recharged to conduits that has mixed with groundwater stored in the pore space of the matrix using a flow and transport model in a synthetic time-lapse resistivity inversion. Comparing the field and synthetic inversions, our results enable us to estimate exchange dynamics, spatial mixing, and flow conditions. Results showed that mixing occurred at a volumetric flux of 56 m3/d with a dispersivity around 1.69 m during the geophysical experiment. For these conditions, it was determined that conduit water composition ranged from 75% groundwater during base flow conditions to less than 50% groundwater in high flow conditions. Though subject to some uncertainties, the time-lapse inversion process provides a means to predict changing hydrologic conditions, leading to mixing of surface water and ground water and thus changes to water quantity and quality, as well as potential for water-rock reactions, in a semiconfined, sink-rise system.

[1]  E. Screaton,et al.  Interactions of diffuse and focused allogenic recharge in an eogenetic karst aquifer (Florida, USA) , 2012, Hydrogeology Journal.

[2]  P. Jeannin,et al.  A Global Experimental System Approach of Karst Springs1 Hydrographs and Chemographs , 2002, Ground water.

[3]  J. McDonnell,et al.  Quantifying contributions to storm runoff through end‐member mixing analysis and hydrologic measurements at the Panola Mountain Research Watershed (Georgia, USA) , 2001 .

[4]  M. Schoor.  Detection of sinkholes using 2D electrical resistivity imaging , 2002 .

[5]  A. Binley,et al.  Vadose zone flow model parameterisation using cross-borehole radar and resistivity imaging , 2001 .

[6]  Barbara J. Denahan,et al.  Electrical resistivity investigation of potential cavities underlying a proposed ash disposal area , 1984 .

[7]  M. Stewart,et al.  Geophysical signatures of fracture traces in a karst aquifer (Florida, U.S.A.) , 1983 .

[8]  E. Screaton,et al.  Monitoring well responses to karst conduit head fluctuations: Implications for fluid exchange and matrix transmissivity in the Floridan aquifer , 2006 .

[9]  A. Sato,et al.  Three‐dimensional spatial and temporal monitoring of soil water content using electrical resistivity tomography , 2001 .

[10]  James P. McNamara,et al.  Application of time-lapse ERT imaging to watershed characterization , 2008 .

[11]  J. Girard,et al.  Locating water-filled karst caverns and estimating their volume using magnetic resonance soundings , 2008 .

[12]  F. Šumanovac,et al.  Evaluation of resistivity and seismic methods for hydrogeological mapping in karst terrains , 2001 .

[13]  R. Maxwell,et al.  Stochastic environmental risk analysis: an integrated methodology for predicting cancer risk from contaminated groundwater , 1999 .

[14]  S. Ashby,et al.  A parallel multigrid preconditioned conjugate gradient algorithm for groundwater flow simulations , 1996 .

[15]  Nicolas Florsch,et al.  Least squares inversion of self‐potential (SP) data and application to the shallow flow of ground water in sinkholes , 2006 .

[16]  R. Maxwell,et al.  Integrated surface-groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model , 2006 .

[17]  L. Cathles,et al.  Electrical conductivity in shaly sands with geophysical applications , 1998 .

[18]  The study of karstic aquifers by geodetic measurements in bus de la genziana station – cansiglio plateau (northeastern Italy) , 2011 .

[19]  Jennifer Martin.  QUANTIFICATION OF THE MATRIX HYDRAULIC CONDUCTIVITY IN THE SANTA FE RIVER SINK/RISE SYSTEM WITH IMPLICATIONS FOR THE EXCHANGE OF WATER BETWEEN THE MATRIX AND CONDUITS , 2003 .

[20]  Claire Welty,et al.  Revisiting the Cape Cod bacteria injection experiment using a stochastic modeling approach. , 2005, Environmental science & technology.

[21]  R. Maxwell.  A terrain-following grid transform and preconditioner for parallel, large-scale, integrated hydrologic modeling , 2013 .

[22]  Roelof Versteeg,et al.  Improved hydrogeophysical characterization and monitoring through parallel modeling and inversion of time-domain resistivity andinduced-polarization data , 2010 .

[23]  Jonathan B. Martin,et al.  Exchange of water between conduits and matrix in the Floridan aquifer , 2001 .

[24]  G. Kaufmann.  A model comparison of karst aquifer evolution for different matrix-flow formulations , 2003 .

[25]  R. Maxwell,et al.  Visualization of conduit‐matrix conductivity differences in a karst aquifer using time‐lapse electrical resistivity , 2012 .

[26]  Wendy D. Graham,et al.  A particle-tracking scheme for simulating pathlines in coupled surface-subsurface flows , 2013 .

[27]  E. Screaton,et al.  Conduit enlargement in an eogenetic karst aquifer , 2010 .

[28]  D. Ford,et al.  Karst Hydrogeology and Geomorphology: Ford/Karst Hydrogeology and Geomorphology , 2007 .

[29]  E. Screaton,et al.  Conduit Properties and Karstification in the Unconfined Floridan Aquifer , 2004, Ground water.

[30]  Andreas Kemna,et al.  Time-lapse three-dimensional inversion of complex conductivity data using an active time constrained (ATC) approach , 2011 .

[31]  E. Screaton,et al.  Time variant cross correlation to assess residence time of water and implication for hydraulics of a sink‐rise karst system , 2011 .

[32]  M. Roth,et al.  Evaluation of site characterization methods for sinkholes in Pennsylvania and New Jersey , 1999 .

[33]  M. Bakalowicz.  Karst groundwater: a challenge for new resources , 2005 .

[34]  James A. Kuiper,et al.  Modeling complex flow in a karst aquifer , 2006 .

[35]  Analytical 1D dual-porosity equivalent solutions to 3D discrete single-continuum models. Application to karstic spring hydrograph modelling , 2002, 1108.5484.

[36]  C. Braitenberg,et al.  Karst deformations due to environmental factors: evidences from the horizontal pendulums of Grotta Gigante, Italy , 2012 .

[37]  Gregory S. Springer.  A pipe-based, first approach to modeling closed conduit flow in caves , 2004 .

[38]  Jim E. Jones,et al.  Newton–Krylov-multigrid solvers for large-scale, highly heterogeneous, variably saturated flow problems , 2001 .

[39]  R. Azzam,et al.  Surface geophysics and borehole inspection as an aid to characterizing karst voids and vadose ventilation patterns (Nerja research site, S. Spain) , 2012 .

[40]  Gerard T. Schuster,et al.  Migration methods for imaging different‐order multiples , 2007 .

[41]  L. Oudin,et al.  Tilt and strain deformation induced by hydrologically active natural fractures: Application to the tiltmeters installed in Sainte-Croix-aux-Mines observatory (France) , 2009 .

[42]  A. Binley,et al.  Cross-hole electrical imaging of a controlled saline tracer injection , 2000 .

[43]  Junfeng Zhu,et al.  Challenges of using electrical resistivity method to locate karst conduits—A field case in the Inner Bluegrass Region, Kentucky , 2011 .

[44]  H. Jourde,et al.  Water exchange and pressure transfer between conduits and matrix and their influence on hydrodynamics of two karst aquifers with sinking streams , 2010 .

[45]  R. Reeves,et al.  A case study of resistivity and self‐potential signatures of hydrothermal instabilities, Inferno Crater Lake, Waimangu, New Zealand , 2009 .

[46]  A. Binley,et al.  Snowmelt infiltration: monitoring temporal and spatial variability using time-lapse electrical resistivity , 2004 .

[47]  André Revil,et al.  4D time-lapse ERT inversion: introducing combined time and space constraints , 2014 .

[48]  E. Screaton,et al.  River reversals into karst springs: A model for cave enlargement in eogenetic karst aquifers , 2011 .

[49]  A. Binley,et al.  Seasonal variation of moisture content in unsaturated sandstone inferred from borehole radar and resistivity profiles. , 2002 .

[50]  Frédéric Nguyen,et al.  Assessing the contribution of electrical resistivity tomography (ERT) and self-potential (SP) methods for a water well drilling program in fractured/karstified limestones , 2011 .

[51]  J. B. Martin,et al.  Relative importance and chemical effects of diffuse and focused recharge in an eogenetic karst aquifer: an example from the unconfined upper Floridan aquifer, USA , 2009 .

[52]  M. M. Mohan Kumar,et al.  Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding , 2008 .

[53]  André Revil,et al.  IP4DI: A software for time-lapse 2D/3D DC-resistivity and induced polarization tomography , 2013, Comput. Geosci..

[54]  J. Urbanc,et al.  Quantification of karst aquifer discharge components during storm events through end-member mixing analysis using natural chemistry and stable isotopes as tracers , 2006 .

[55]  F. Nobilis.  Dry spells in the Alpine country Austria , 1986 .

[56]  P. Styles,et al.  Integrated high-resolution geophysical investigations as potential tools for water resource investigations in karst terrain , 2002 .

[57]  Pitambar Gautam,et al.  Mapping of subsurface karst structure with gamma ray and electrical resistivity profiles: a case study from Pokhara valley, central Nepal , 2000 .

[58]  Jan Vanderborght,et al.  Proof of concept of regional scale hydrologic simulations at hydrologic resolution utilizing massively parallel computer resources , 2010 .

[59]  M. Boucher,et al.  Geophysical characterisation of karstic networks – Application to the Ouysse system (Poumeyssen, France) , 2009 .

[60]  Y. Albouy,et al.  Localization of Saturated Karst Aquifer with Magnetic Resonance Sounding and Resistivity Imagery , 2003, Ground water.

[61]  Georg Teutsch,et al.  Simulation of the development of karst aquifers using a coupled continuum pipe flow model , 2003 .

[62]  H. Vereecken,et al.  Imaging and characterisation of subsurface solute transport using electrical resistivity tomography (ERT) and equivalent transport models , 2002 .

[63]  W. White,et al.  Storm response of the karstic carbonate aquifer of southcentral Kentucky , 1988 .

[64]  Roger Guerin,et al.  Contribution of geophysical methods to karst-system exploration: an overview , 2011 .

[65]  Y. Benderitter,et al.  Shallow karst exploration using MT-VLF and DC resistivity methods1 , 1995 .

[66]  Panagiotis Tsourlos,et al.  4D active time constrained resistivity inversion , 2011 .

[67]  A. Binley,et al.  Monitoring snowmelt induced unsaturated flow and transport using electrical resistivity tomography , 2002 .

[68]  A. Revil,et al.  Detection of preferential infiltration pathways in sinkholes using joint inversion of self‐potential and EM‐34 conductivity data , 2007 .