Groundwater–surface water interactions in a river estuary and the importance of geomorphology: Insights from hydraulic,thermal and geophysical observations

This study presents the groundwater flow and salinity dynamics along a river estuary, the Werribee River in Victoria, Australia, at local and regional scales. Along a single reach, salinity across a transverse section of the channel (~80 m long) with a point bar was monitored using time-lapse electrical resistivity (ER) through a tidal cycle. Groundwater fluxes were concurrently estimated by monitoring groundwater levels and temperature profiles. Regional porewater salinity distribution was mapped using 6-km long longitudinal ER surveys during summer and winter. The time-lapse ER across the channel revealed a static electrically resistive zone on the side of the channel with a pronounced cut bank. Upward groundwater flux and steep vertical temperature gradients with colder temperatures deeper within the sediment suggested a stable zone of fresh groundwater discharge along this cut bank area. Generally, less resistive zones were observed at the shallow portion of the inner meander bank and at the channel center. Subsurface temperatures close to surface water values, vertical head gradients indicating both upward and downward groundwater flux, and higher porewater salinity closer to that of estuary water suggest strong hyporheic circulation in these zones. The longitudinal surveys revealed higher ER values along deep and sinuous segments and low ER values in shallow and straighter reaches in both summer and winter; these patterns are consistent with the local channel-scale observations. This study highlights the interacting effects of channel morphology, broad groundwater–surface water interaction, and hyporheic exchange on porewater salinity dynamics underneath and adjacent to a river estuary.     

Imaging hyporheic zone solute transport using electrical resistivity

Traditional characterization of hyporheic processes relies upon modelling observed in‐stream and subsurface breakthrough curves to estimate hyporheic zone size and infer exchange rates. Solute data integrate upstream behaviour and lack spatial coverage, limiting our ability to accurately quantify spatially heterogeneous exchange dynamics. Here, we demonstrate the application of near‐surface electrical resistivity imaging (ERI) methods, coupled with experiments using an electrically conductive stream tracer (dissolved NaCl), to provide in situ imaging of spatial and temporal dynamics of hyporheic exchange. Tracer‐labelled water in the stream enters the hyporheic zone, reducing electrical resistivity in the subsurface (to which subsurface ERI is sensitive). Comparison of background measurements with those recording tracer presence provides distributed characterization of hyporheic area (in this application, ∼0·5 m²). Results demonstrate the first application of ERI for two‐dimensional imaging of stream‐aquifer exchange and hyporheic extent. Future application of this technique will greatly enhance our ability to quantify processes controlling solute transport and fate in hyporheic zones, and provide data necessary to inform more complete numerical models. Copyright © 2010 John Wiley & Sons, Ltd.     

Form-based river restoration decreases wetland hyporheic exchange: Lessons learned from the Upper Colorado River

Restoration of river–wetland systems to recover lost ecosystem services and restore consistent flood regimes is commonly directed at modifying in-channel storage and hyporheic exchange. Here, we monitored the hydrologic response to channel realignment in a montane river–wetland system by comparing pre- and post-restoration measurements. In 2015, an earthen berm and 190 m segment of the Upper Colorado River were constructed to consolidate flow from multiple channels into the historic thalweg. We injected a sodium chloride tracer during baseflow and used mass-balance calculations and electrical resistivity imaging to assess changes in near-channel hyporheic exchange. Results indicate a decrease in hyporheic exchange within the wetland due to lost complexity along the consolidated flow path. Subsurface complexity appears to control hyporheic exchange more than surface complexity. Flow consolidation increased the area-adjusted wetland water yield by 231 mm, indicating a loss of wetland water storage capacity. One year of post-restoration monitoring suggests that the form-based channel restoration directed at consolidating flow into a single thread adversely affected the hyporheic exchange functioning in the pre-restoration system. Results from this case study are applicable to restoration planners as they consider the effects of form-based projects on water storage capacity in similar systems. © 2018 John Wiley & Sons, Ltd.     

Electrical resistivity imaging of hydrologic flow through surface coal mine valley fills with comparison to other landforms

Surface coal mining has altered land cover, near‐surface geologic structure, and hydrologic processes of large areas in central Appalachia, USA. These alterations are associated with changes in water quality such as elevated total‐dissolved solids, which is usually measured via its surrogate, specific conductance (SC). The SC of valley fill effluent streams is a function of fill construction methods, materials, and age; yet hydrologic studies that relate these variables to water quality are sparse due to the difficulty of conducting traditional hydrologic studies in mined landscapes. We used electrical resistivity imaging (ERI) to visualize the subsurface geologic structure and hydrologic flow paths within a valley fill. ERI is a noninvasive geophysical technique that maps spatiotemporal changes in resistivity of the subsurface. We paired ERI with artificial rainfall experiments to track infiltrated water as it moved through the valley fill. Results indicate that ERI can be used to identify subsurface geologic structure and track advancing wetting fronts or preferential flow paths. Our results suggest that the upper portion of the fill contains significant fines, whereas the deeper profile is primarily large rocks and void spaces. Water tended to pond on the surface of compacted areas until it reached preferential flow paths, where it appeared to infiltrate quickly down to >15 m depth in 75 min. ERI applications can improve understanding of how fill construction techniques influence subsurface water movement, and in turn may aid in the development of valley fill construction methods to reduce water quality effects.     

A Comparison of Hyporheic Transport at a Cross‐Vane Structure and Natural Riffle

While restoring hyporheic flowpaths has been cited as a benefit to stream restoration structures, little documentation exists confirming that constructed restoration structures induce comparable hyporheic exchange to natural stream features. This study compares a stream restoration structure (cross‐vane) to a natural feature (riffle) concurrently in the same stream reach using time‐lapsed electrical resistivity (ER) tomography. Using this hydrogeophysical approach, we were able to quantify hyporheic extent and transport beneath the cross‐vane structure and the riffle. We interpret from the geophysical data that the cross‐vane and the natural riffle induced spatially and temporally unique hyporheic extent and transport, and the cross‐vane created both spatially larger and temporally longer hyporheic flowpaths than the natural riffle. Tracer from the 4.67‐h injection was detected along flowpaths for 4.6 h at the cross‐vane and 4.2 h at the riffle. The spatial extent of the hyporheic zone at the cross‐vane was 12% larger than that at the riffle. We compare ER results of this study to vertical fluxes calculated from temperature profiles and conclude significant differences in the interpretation of hyporheic transport from these different field techniques. Results of this study demonstrate a high degree of heterogeneity in transport metrics at both the cross‐vane and the riffle and differences between the hyporheic flowpath networks at the two different features. Our results suggest that restoration structures may be capable of creating sufficient exchange flux and timescales of transport to achieve the same ecological functions as natural features, but engineering of the physical and biogeochemical environment may be necessary to realize these benefits.     

Application of ERT,Saline Tracer and Numerical Studies to Delineate Preferential Paths in Fractured Granites

Accurate quantification of in situ heterogeneity and flow processes through fractured geologic media remains elusive for hydrogeologists due to the complexity in fracture characterization and its multiscale behavior. In this research, we demonstrated the efficacy of tracer-electrical resistivity tomography (ERT) experiments combined with numerical simulations to characterize heterogeneity and delineate preferential flow paths in a fractured granite aquifer. A series of natural gradient saline tracer experiments were conducted from a depth window of 18 to 22 m in an injection well (IW) located inside the Indian Institute of Technology Hyderabad campus. Tracer migration was monitored in a time-lapse mode using two cross-sectional surface ERT profiles placed in the direction of flow gradient. ERT data quality was improved by considering stacking, reciprocal measurements, resolution indicators, and geophysical logs. Dynamic changes in subsurface electrical properties inferred via resistivity anomalies were used to highlight preferential flow paths of the study area. Temporal changes in electrical resistivity and tracer concentration were monitored along the vertical in an observation well located at 48 m to the east of the IW. ERT-derived tracer breakthrough curves were in agreement with geochemical sample measurements. Fracture geometry and hydraulic properties derived from ERT and pumping tests were further used to evaluate two mathematical conceptualizations that are relevant to fractured aquifers. Results of numerical analysis conclude that dual continuum model that combines matrix and fracture systems through a flow exchange term has outperformed equivalent continuum model in reproducing tracer concentrations at the monitoring wells (evident by a decrease in RMSE from 199 to 65 mg/L). A sensitivity analysis on model simulations conclude that spatial variability in hydraulic conductivity, local-scale dispersion, and flow exchange at fracture-matrix interface have a profound effect on model simulations.     

Sensitivity of Simulated Hyporheic Exchange to River Bathymetry: The Steinlach River Test Site

This study determines the aspects of river bathymetry that have the greatest influence on the predictive biases when simulating hyporheic exchange. To investigate this, we build a highly parameterized HydroGeoSphere model of the Steinlach River Test Site in southwest Germany as a reference. This model is then modified with simpler bathymetries, evaluating the changes to hyporheic exchange fluxes and transit time distributions. Results indicate that simulating hyporheic exchange with a high-resolution detailed bathymetry using a three-dimensional fully coupled model leads to nested multiscale hyporheic exchange systems. A poorly resolved bathymetry will underestimate the small-scale hyporheic exchange, biasing the simulated hyporheic exchange towards larger scales, thus leading to overestimates of hyporheic exchange residence times. This can lead to gross biases in the estimation of a catchment's capacity to attenuate pollutants when extrapolated to account for all meanders along an entire river within a watershed. The detailed river slope alone is not enough to accurately simulate the locations and magnitudes of losing and gaining river reaches. Thus, local bedforms in terms of bathymetric highs and lows within the river are required. Bathymetry surveying campaigns can be more effective by prioritizing bathymetry measurements along the thalweg and gegenweg of a meandering channel. We define the gegenweg as the line that connects the shallowest points in successive cross-sections along a river opposite to the thalweg under average flow conditions. Incorporating local bedforms will likely capture the nested nature of hyporheic exchange, leading to more physically meaningful simulations of hyporheic exchange fluxes and transit times.     

Surface water-groundwater exchange dynamics in buried-valley aquifer systems

Groundwater is a primary source of drinking water worldwide, but excess nutrients and emerging contaminants could compromise groundwater quality and limit its usage as a drinking water source. As such contaminants become increasingly prevalent in the biosphere, a fundamental understanding of their fate and transport in groundwater systems is necessary to implement successful remediation strategies. The dynamics of surface water-groundwater (hyporheic) exchange within a glacial, buried-valley aquifer system are examined in the context of their implications for the transport of nutrients and contaminants in riparian sediments. High conductivity facies act as preferential flow pathways which enhance nutrient and contaminant delivery, especially during storm events, but transport throughout the aquifer also depends on subsurface sedimentary architecture (e.g. interbedded high and low conductivity facies). Temperature and specific conductance measurements indicate extensive hyporheic mixing close to the river channel, but surface water influence was also observed far from the stream-aquifer interface. Measurements of river stage and hydraulic head indicate that significant flows during storms (i.e., hot moments) alter groundwater flow patterns, even between consecutive storm events, as riverbed conductivity and, more importantly, the hydraulic connectivity between the river and aquifer change. Given the similar mass transport characteristics among buried-valley aquifers, these findings are likely representative of glacial aquifer systems worldwide. Our results suggest that water resources management decisions based on average (base) flow conditions may inaccurately represent the system being evaluated, and could reduce the effectiveness of remediation strategies for nutrients and emerging contaminants.     

Spatial and temporal dynamics of nitrogen exchange in an upwelling reach of a groundwater-fed river and potential response to perturbations changing rainfall patterns under UK climate change scenarios

We report the complex spatial and temporal dynamics of hyporheic exchange flows (HEFs) and nitrogen exchange in an upwelling reach of a 200 m groundwater-fed river. We show how research combining hydrological measurement, geophysics and isotopes, together with nutrient speciation techniques provides insight on nitrogen pathways and transformations that could not have been captured otherwise, including a zone of vertical preferential discharge of nitrate from deeper groundwater, and a zone of rapid denitrification linking the floodplain with the riverbed. Nitrate attenuation in the reach is dominated by denitrification but is spatially highly variable. This variability is driven by groundwater flow pathways and landscape setting, which influences hyporheic flow, residence time and nitrate removal. We observed the spatial connectivity of the river to the riparian zone is important because zones of horizontal preferential discharge supply organic matter from the floodplain and create anoxic riverbed conditions with overlapping zones of nitrification potential and denitrification activity that peaked 10–20 cm below the riverbed. Our data also show that temporal variability in water pathways in the reach is driven by changes in stage of the order of tens of centimetres and by strength of water flux, which may influence the depth of delivery of dissolved organic carbon. The temporal variability is sensitive to changes to river flows under UK climate projections that anticipate a 14%–15% increase in regional median winter rainfall and a 14%–19% reduction in summer rainfall. Superimposed on seasonal projections is more intensive storm activity that will likely lead to a more dynamic and inherently complex (hydrologically and biogeochemically) hyporheic zone. We recorded direct evidence of suppression of upwelling groundwater (flow reversal) during rainfall events. Such flow reversal may fuel riverbed sediments whereby delivery of organic carbon to depth, and higher denitrification rates in HEFs might act in concert to make nitrate removal in the riverbed more efficient.     

Observing lingering hyporheic storage using electrical resistivity: variations around stream restoration structures,Crabby Creek,PA

Time‐lapse geophysical surveys can map lingering hyporheic storage by detecting changes in response to saline tracer. Tracer tests were conducted in Crabby Creek, an urban stream outside Philadelphia, to examine the influence of stream restoration structures and variable sediment thickness. We compared electrical resistivity surveys with extensive well sampling (57 wells) in two 13.5‐m‐long reaches, each with a step drop created by a J‐hook. The two step drops varied in tracer behaviour, based on both the well data and the geophysical data. The well data showed more variation in arrival time where the streambed sediment was thick and was more uniform where sediment was thin. The resistivity in the reach with thin sediment showed lingering tracer in the hyporheic zone both upstream and downstream from the J‐hook. In the second reach where the sediment was thicker, the lingering tracer in the hyporheic zone was more extensive downstream from the J‐hook. The contrasting results between the two reaches from both methods suggested that sediments influenced hyporheic exchange more than the step at this location. Resistivity inversion differed from well data in both reaches in that it showed evidence for tracer after well samples had returned to background, mapping lingering tracer either upstream or downstream of a step. We conclude that resistivity surveys may become an important tool for hyporheic zone characterization because they provide information on the extent of slow moving fluids in the hyporheic zone, which have the potential to enhance chemical reactions. Copyright © 2012 John Wiley & Sons, Ltd.     

Geomorphology,hydrology, and aquatic vegetation drive seasonal hyporheic flow patterns across a gravel‐dominated floodplain

Across 1·7 km² of the Umatilla River floodplain (Oregon, USA), we investigated the influences of an ephemeral tributary and perennial ‘spring channel’ (fed only by upwelling groundwater) on hyporheic hydrology. We derived maps of winter and summer water‐table elevations from data collected at 46 monitoring wells and 19 stage gauges and used resulting maps to infer groundwater flow direction. Groundwater flow direction varied seasonally across the floodplain and was influenced by main channel stage, flooding, the tributary creek, and the location and direction of hyporheic exchange in the spring channel. Hyporheic exchange in the spring channel was evaluated with a geochemical mixing model, which confirmed patterns of floodplain groundwater movement inferred from water‐table maps and showed that the spring channel was fed predominantly by hyporheic water from the floodplain aquifer (87% during winter, 80% during summer), with its remaining flow supplied by upslope groundwater from the adjacent catchment aquifer. Summertime growth of aquatic macrophytes in the spring channel also influenced patterns of hyporheic exchange and groundwater flow direction in the alluvial aquifer by increasing flow resistance in the spring channel, locally raising surface water stage and adjacent water‐table elevation, and thereby altering the slope of the water‐table in the hyporheic zone. The Umatilla River floodplain is larger than most sites where hyporheic hydrology has been investigated in detail. Yet, our results corroborate other research that has identified off‐channel geomorphic features as important drivers of hyporheic hydrology, including previously published modeling efforts from a similar river and field observations from smaller streams. Copyright © 2007 John Wiley & Sons, Ltd.     

Application of electrical resistivity imaging (ERI) for the assessment of peat properties: a case study of the Całowanie Fen,Central Poland

Complex studies were carried out to recognize the fen structure and peat properties in the Ca?owanie Fen area, belonging to the Natura 2000 network. The studies were conducted in two study areas that differ significantly in terms of peat thickness. Electrical resistivity imaging (ERI) was used to identify the properties of the peat and its substrate, such as thickness and electrical resistivity. Comparison of the field studies with the laboratory tests has shown that the ash content rises electrical resistivity in peat. In addition, the study has shown that the application of non-invasive geophysical methods in protected areas is justified. The fen, as a medium containing mostly water, was a proper test area for the ERI measurements.     

Effect of streamflow stochasticity on bedform-driven hyporheic exchange

The interactions between the stream and the geomorphologic units that compose the stream channel result in an exchange of water, heat, and chemicals that is an important component of the flows of energy and nutrients in the river ecosystem. This exchange is characterized by complex spatial and temporal dynamics that depend on the characteristics of the stream flow and morphology. At present, many studies have addressed the development of spatial patterns of hyporheic exchange that are induced by many geomorphological factors at different scales. However, much less is known about the temporal evolution of the surface–subsurface exchange in response to the dynamics of the stream discharge. In order to investigate this problem, the present work analyzes the influence of streamflow variability on the hyporheic exchange induced by fluvial bedforms. A stochastic approach is employed to generate streamflow series whose statistical properties are representative of streams with different hydrological regimes. The resulting exchange fluxes and travel times are then computed, and the relationships between the streamflow regime and the dynamics of the exchange flux and travel times are investigated. The results show that the mean stream discharge can be used to estimate the average features of the temporal dynamics of hyporheic exchange. Moreover, exchange fluxes and residence times distributions exhibit significant fluctuations, which are tightly related to the coefficient of variation of the streamflow hydrograph.     

Feasibility Assessment of Long-Term Electrical Resistivity Monitoring of a Nitrate Plume

Long-term monitoring solutions at contaminated sites are necessary to track plume migration and evaluate the performance of remediation efforts. Electrical resistivity imaging (ERI) can potentially provide information about plume dynamics; however, the feasibility and likelihood of success are seldom evaluated before conducting a field study. Coupling flow and transport models with geoelectrical models provide a powerful way to assess the potential effectiveness of an actual ERI field campaign. We present a coupled approach for evaluating the feasibility of monitoring nitrate migration and remediation using 4D time-lapse ERI at a legacy nuclear waste facility. This kilometer-scale study focuses on depths below the water table (∼70 m). A flow and transport model is developed to perform simulations of nitrate migration and removal via a hypothetical pump-and-treat system. A tracer injection is also simulated at the leading edge of the nitrate plume to enhance the conductivity contrast between the native subsurface and the groundwater fluids. Images of absolute bulk conductivity provide limited information concerning plume migration while time-lapse difference images, which remove the static effects of geology, provide more useful information concerning plume dynamics over time. A spatial moment analysis performed on flow and transport and ERI models matches well during the tracer injection; however, inversion regularization smoothing otherwise limits the value in terms of locating the center of mass. We find that the addition of a tracer enables ERI to characterize plume dynamics during pump-and-treat operations, and late-time ERI monitoring provides a conservative estimate of nitrate plume boundaries in this synthetic study.     

Wetland-Scale Mapping of Preferential Fresh Groundwater Discharge to the Colorado River

Quantitative evaluation of groundwater/surface water exchange dynamics is universally challenging in large river systems, because existing methodology often does not yield spatially-distributed data and is difficult to apply in deeper water. Here we apply a combined near-surface geophysical and direct groundwater chemical toolkit to refine fresh groundwater discharge estimates to the Colorado River through a 4-km² wetland that borders the town of Moab, Utah, USA. Preliminary characterization of raw electromagnetic imaging (EMI) data, collected by kayak and by walking, was used to guide additional direct-contact electrical measurements and installation of new monitoring wells. Chemical data from the wells strongly supported the EMI spatial characterization of preferential fresh groundwater discharge embedded in natural brine groundwaters and weighted to the southern wetland section. Inversion of the EMI data revealed sub-meter scale detail regarding bulk electrical conductivity zonation across approximately 15.5 km of transects, collected in only 3 days. This electrical detail indicates processes such as salinization of the unsaturated zone and direct discharge through the Colorado River sediments and a tributary creek bed. Overall, the study contributed to a substantial reduction in fresh groundwater discharge estimates previously made using sparse existing well data and a simplified assumption of diffuse fresh groundwater discharge below the entire wetland. EMI will likely become a widely used tool in systems with natural electrical contrast as groundwater/surface water hydrogeologists continue to recognize the prevalence of preferential groundwater discharge processes.     

Impact of a low-permeability lens on dune-induced hyporheic exchange

Hyporheic exchange induced by dunes is a key process controlling water fluxes and biogeochemical processes in river networks. Owing to the limitations of instrumental detection at small spatial scales, previous studies have focused mainly on dune-induced hyporheic exchange in homogeneous systems. A low-permeability lens is a natural, widespread heterogeneity in stream beds, and probably affects the processes of water flow and contaminant transportation significantly. To quantitatively analyse the response mechanism of hyporheic exchange to a low-permeability lens, a two-dimensional dune-generated hyporheic exchange model was developed using the VS2DH model. The results indicate a lens in a stream bed can hinder or enhance hyporheic exchange processes, depending on its relative spatial location to dunes. Both the increase in length and thickness of the lens could strengthen its impacts on hyporheic exchange. Regional groundwater–surface water interactions of higher intensity suppress the flow of hyporheic exchange in a stream bed with a low-permeability lens.     

Transient streambed hydraulic conductivity in channel and bar environments,Loup River,Nebraska

Streambed hydraulic conductivity (K) and vertical K (K _v) are key controls on groundwater and surface water exchange and biogeochemical fluxes through the hyporheic zone, but drivers of transient hydraulic properties in different fluvial environments are poorly understood. This study combines hydrogeology, geophysics, and sedimentology to reveal mechanisms of K and K _v transience in the upper 0.5 m of a sandy streambed during low discharge. Hydraulic tests (44 slug tests, 130 falling-head permeameter tests) and 130 grain-size analyses were repeated three times over 8 weeks on a 1,200 m² grid spanning: (a) a channel with continuously flowing water and mobile bed load; (b) an adjacent mid-channel bar that was stationary and infrequently submerged. Aerial photographs and ground-penetrating radar show scour and complete reworking of fluvial sediments in the channel. Bar sediments below the water table remained immobile, but infrequent flows of moderate discharge reworked the uppermost few centimetres of the bar top. Despite differences in sediment mobility and stream flow characteristics across environments, K and K _v exhibited order-of-magnitude differences in spatial heterogeneity and temporal variability in both the channel and bar. Mean K and K _v values in the channel were comparatively stable over time. In the immobile bar, mean K declined 20% and K _v declined 26% after increased discharge temporarily inundated the bar. Grain-size distributions were steady across both environments over time, but repeat geophysical surveys of the bar show a decrease in electrical conductivity, likely from porosity reduction. These findings suggest that sediment dynamics and stream flow characteristics in different streambed environments are important drivers of K transience during low discharge conditions. Specifically, pore clogging can be an important mechanism of transience over short durations (weeks to months) in immobile sediments subject to infrequent flows and minor reworking.     

The effects of hydropeaking on hyporheic interactions based on field experiments

Hydropeaking power production has the potential to pose serious challenges towards hydrology, water quality and ecology in the downstream water bodies. The effects of such abrupt changes of flow in hyporheic exchange have been explored in a few cases in the literature. This paper extends previous works with a study of finer time resolution in a river of a smaller size and with different climatic characteristics, adding to the current knowledge of peaking‐hyporheic interactions. A high‐frequency logging field experiment measuring hyporheic flow and temperature was conducted on a ~30 × 20 m gravel bar frequently exposed to dry conditions because of fast and abrupt flow changes. This study demonstrates that hyporheic processes are sensitive to hydropeaking with respect to rates of change, durations and temperature. Differences BETWEEN individual events, seasons, watering and dewatering processes and positions in the river bed that can be potentially relevant to ecology were investigated. Understanding the complexity of those processes at the fine scale from the physical point of view is both important for the judgement of potential ecological impacts and for the future management of such regulated systems. Copyright © 2014 John Wiley & Sons, Ltd.