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

被引:45
作者
Meyerhoff, Steven B. [1 ,4 ]
Maxwell, Reed M. [1 ]
Revil, Andre [2 ,3 ]
Martin, Jonathan B. [5 ]
Karaoulis, Marios [2 ]
Graham, Wendy D. [6 ]
机构
[1] Colorado Sch Mines, Dept Geol & Geol Engn, Hydrol Sci & Engn Program, Golden, CO 80401 USA
[2] Colorado Sch Mines, Dept Geophys, Golden, CO 80401 USA
[3] Univ Savoie, CNRS, UMR 5275, ISTerre, Le Bourget Du Lac, France
[4] Itasca Denver Inc, Lakewood, CO USA
[5] Univ Florida, Dept Geol Sci, Gainesville, FL USA
[6] Univ Florida, Water Inst, Gainesville, FL USA
基金
美国国家科学基金会;
关键词
groundwater and surface water interactions; geophysics; electrical resistivity tomography; synthetic modeling; BOREHOLE RADAR; DC-RESISTIVITY; LARGE-SCALE; FLOW; PARALLEL; MATRIX; MODEL; ERT; CONDUCTIVITY; INVERSION;
D O I
10.1002/2013WR013991
中图分类号
X [环境科学、安全科学];
学科分类号
08 ; 0830 ;
摘要
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 m(3)/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.
引用
收藏
页码:2566 / 2585
页数:20
相关论文
共 75 条
[1]   A parallel multigrid preconditioned conjugate gradient algorithm for groundwater flow simulations [J].
Ashby, SF ;
Falgout, RD .
NUCLEAR SCIENCE AND ENGINEERING, 1996, 124 (01) :145-159
[2]   Time variant cross correlation to assess residence time of water and implication for hydraulics of a sink-rise karst system [J].
Bailly-Comte, V. ;
Martin, Jonathan B. ;
Screaton, E. J. .
WATER RESOURCES RESEARCH, 2011, 47
[3]   Water exchange and pressure transfer between conduits and matrix and their influence on hydrodynamics of two karst aquifers with sinking streams [J].
Bailly-Comte, Vincent ;
Martin, Jonathan B. ;
Jourde, Herve ;
Screaton, Elizabeth J. ;
Pistre, Severin ;
Langston, Abigail .
JOURNAL OF HYDROLOGY, 2010, 386 (1-4) :55-66
[4]   Karst groundwater: a challenge for new resources [J].
Bakalowicz, M .
HYDROGEOLOGY JOURNAL, 2005, 13 (01) :148-160
[5]   Seasonal variation of moisture content in unsaturated sandstone inferred from borehole radar and resistivity profiles [J].
Binley, A ;
Winship, P ;
West, LJ ;
Pokar, M ;
Middleton, R .
JOURNAL OF HYDROLOGY, 2002, 267 (3-4) :160-172
[6]   Vadose zone flow model parameterisation using cross-borehole radar and resistivity imaging [J].
Binley, A ;
Cassiani, G ;
Middleton, R ;
Winship, P .
JOURNAL OF HYDROLOGY, 2002, 267 (3-4) :147-159
[7]   Quantifying contributions to storm runoff through end-member mixing analysis and hydrologic measurements at the Panola Mountain Research Watershed (Georgia, USA) [J].
Burns, DA ;
McDonnell, JJ ;
Hooper, RP ;
Peters, NE ;
Freer, JE ;
Kendall, C ;
Beven, K .
HYDROLOGICAL PROCESSES, 2001, 15 (10) :1903-1924
[8]  
Chalikakis K, 2011, HYDROGEOL J, V19, P1169, DOI 10.1007/s10040-011-0746-x
[9]   Analytical 1D dual-porosity equivalent solutions to 3D discrete single-continuum models. Application to karstic spring hydrograph modelling [J].
Cornaton, F ;
Perrochet, P .
JOURNAL OF HYDROLOGY, 2002, 262 (1-4) :165-176
[10]   A particle-tracking scheme for simulating pathlines in coupled surface-subsurface flows [J].
de Rooij, Rob ;
Graham, Wendy ;
Maxwell, Reed M. .
ADVANCES IN WATER RESOURCES, 2013, 52 :7-18