TY - JOUR
T1 - Groundwater flow and mixing in a wetland-stream system: Field study and numerical modeling
AU - Karan, Sachin
AU - Engesgaard, Peter
AU - Looms, Majken C.
AU - Laier, Troels
AU - Kazmierczak, Jolanta
N1 - Funding Information:
We would like to thank Louise A. Andreasen, Mie Andreasen, Astrid M. Bisgaard, Carlos Duque, Valentina Marconi, Sascha Müller, Helene H.M. Nielsen, Eva Sebok and the students who have been helping over the years with collecting field data. We also thank an anonymous reviewer and the associate editor. This work is part of the Hydrology Observatory, HOBE ( http://www.hobecenter.dk ) funded by the Villum Foundation.
PY - 2013/4/30
Y1 - 2013/4/30
N2 - We combined electrical resistivity tomography (ERT) on land and in a stream with zone-based hydraulic conductivities (from multi-level slug testing) to investigate the local geological heterogeneity of the deposits in a wetland-stream system. The detailed geology was incorporated into a numerical steadystate groundwater model that was calibrated against average head observations. The model results were tested against groundwater fluxes determined from streambed temperature measurements. Discharge varied up to one order of magnitude across the stream and the model was successful in capturing this variability. Water quality analyses from multi-level sampling underneath the streambed and in the wetland showed a stratification in groundwater composition with an aerobic shallow zone with oxygen and nitrate (top ~3 m) overlying a reduced, anoxic zone. While NO
-
3 concentrations up to 58 mg L
-1 were found in the top of the aquifer and immediately underneath the streambed no NO
-
3 was detected deeper within the aquifer. An inverse relationship between NO
-
3 and SO
2-
4 suggests that pyrite oxidation takes place in the deeper parts of the aquifer. Simulated flow path lines showed very different origins for deeper groundwater samples. No nitrate reduction is believed to occur in the shallow zone, where oxygen is present, and the residence time is on the order of 1 year. Nitrate reduction can, however, occur in the deeper parts, which are oxygen-free, and where the residence time is on the order of 7 years. A simulation with a homogeneous model did not match the observations nearly as well as a heterogeneous model based on ERT and a spatially distributed hydraulic conductivity. Furthermore, the origin of the sampled groundwater could not have been predicted from groundwater hydraulic head and the groundwater chemistry alone. The presented approach of integrating such methods in groundwater-surface water exchange studies, proved efficient to obtain information of the controlling factors.
AB - We combined electrical resistivity tomography (ERT) on land and in a stream with zone-based hydraulic conductivities (from multi-level slug testing) to investigate the local geological heterogeneity of the deposits in a wetland-stream system. The detailed geology was incorporated into a numerical steadystate groundwater model that was calibrated against average head observations. The model results were tested against groundwater fluxes determined from streambed temperature measurements. Discharge varied up to one order of magnitude across the stream and the model was successful in capturing this variability. Water quality analyses from multi-level sampling underneath the streambed and in the wetland showed a stratification in groundwater composition with an aerobic shallow zone with oxygen and nitrate (top ~3 m) overlying a reduced, anoxic zone. While NO
-
3 concentrations up to 58 mg L
-1 were found in the top of the aquifer and immediately underneath the streambed no NO
-
3 was detected deeper within the aquifer. An inverse relationship between NO
-
3 and SO
2-
4 suggests that pyrite oxidation takes place in the deeper parts of the aquifer. Simulated flow path lines showed very different origins for deeper groundwater samples. No nitrate reduction is believed to occur in the shallow zone, where oxygen is present, and the residence time is on the order of 1 year. Nitrate reduction can, however, occur in the deeper parts, which are oxygen-free, and where the residence time is on the order of 7 years. A simulation with a homogeneous model did not match the observations nearly as well as a heterogeneous model based on ERT and a spatially distributed hydraulic conductivity. Furthermore, the origin of the sampled groundwater could not have been predicted from groundwater hydraulic head and the groundwater chemistry alone. The presented approach of integrating such methods in groundwater-surface water exchange studies, proved efficient to obtain information of the controlling factors.
KW - Flow paths
KW - Groundwater-stream exchange
KW - Numerical modeling
KW - Nutrients
KW - Riparian zone
UR - http://www.scopus.com/inward/record.url?scp=84886101714&partnerID=8YFLogxK
U2 - 10.1016/j.jhydrol.2013.02.030
DO - 10.1016/j.jhydrol.2013.02.030
M3 - Article
SN - 0022-1694
VL - 488
SP - 73
EP - 83
JO - Journal of Hydrology
JF - Journal of Hydrology
ER -