Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Joint interpretation of hydrological and geophysical data: electrical resistivity tomography results from a process hydrological research site in the Black Forest Mountains,Germany

The use of electrical resistivity tomography (ERT; non‐intrusive geophysical technique) was assessed to identify the hydrogeological conditions at a surface water/groundwater test site in the southern Black Forest, Germany. A total of 111 ERT transects were measured, which adopted electrode spacings from 0·5 to 5 m as well as using either Wenner or dipole‐dipole electrode arrays. The resulting two‐dimensional (2D) electrical resistivity distributions are related to the structure and water content of the subsurface. The images were interpreted with respect to previous classical hillslope hydrological investigations within the same research basin using both tracer methods and groundwater level observations. A raster‐grid survey provided a quasi 3D resistivity pattern of the floodplain. Strong structural heterogeneity of the subsurface could be demonstrated, and (non)connectivities between surface and subsurface bodies were mapped. Through the spatial identification of likely flow pathways and source areas of runoff, the deep groundwater within the steeper valley slope seems to be much more connected to runoff generation processes within the valley floodplain than commonly credited in such environmental circumstances. Further, there appears to be no direct link between subsurface water‐bodies adjacent to the stream channel. Deep groundwater sources are also able to contribute towards streamflow from exfiltration at the edge of the floodplain as well as through the saturated areas overlying the floodplain itself. Such exfiltrated water then moves towards the stream as channelized surface flow. These findings support previous tracer investigations which showed that groundwater largely dominates the storm hydrograph of the stream, but the source areas of this component were unclear without geophysical measurements. The work highlighted the importance of using information from previous, complementary hydrochemical and hydrometric research campaigns to better interpret the ERT measurements. On the other hand, the ERT can provide a better spatial understanding of existing hydrochemical and hydrometric data. Copyright © 2009 John Wiley & Sons, Ltd.     

Seasonal variations in groundwater upwelling zones in a Danish lowland stream analyzed using Distributed Temperature Sensing (DTS)

The distribution of groundwater inflows in a stream reach plays a major role in controlling the stream temperature, a vital component shaping the riverine ecosystem. In this study, the Distributed Temperature Sensing (DTS) system was installed in a small Danish lowland stream, Elverdamsåen, to assess the seasonal dynamics of groundwater inflow zones using high spatial (1 m) and temporal (3 minutes) resolution of water temperature measurements. Four simple criteria consisting of 30 min average temperature at 16:00, mean and standard deviation of diurnal temperatures, and the day–night temperature difference were applied to three DTS datasets representing stream temperature responses to the variable meteorological and hydrological conditions prevailing in summer, winter and spring. The standard deviation criterion was useful to identify groundwater discharge zones in summer and spring conditions, while the mean temperature criterion was better for the winter conditions. In total, 20 interactions were identified from the DTS datasets representing summer, 16 in winter and 19 in spring, albeit with only two interactions contributing in all three seasons. Higher baseflow to streamflow ratio, antecedent precipitation and presence of fractured clayey till in the stream reach were deemed as the vital factors causing apparent seasonal variation in the locations of upwelling zones, prompting use of DTS not only in preconceived scenarios of large diurnal temperature change but rather a long‐term deployment covering variable meteorological and hydrological scenarios. Copyright © 2012 John Wiley & Sons, Ltd.     

Does hillslope trenching enhance groundwater recharge and baseflow in the Peruvian Andes?

As Andean glaciers rapidly retreat due to climate change, the balance of groundwater and glacial meltwater contributions to stream discharge in tropical, proglacial watersheds will change, potentially increasing vulnerability of water resources. The Shullcas River Watershed, near Huancayo, Peru, is fed only partly by the rapidly receding Huaytapallana glaciers (<20% of dry season flow). To potentially increase recharge and therefore increase groundwater derived baseflow, the government and not‐for‐profit organizations have installed trenches along large swaths of hillslope in the Shullcas Watershed. Our study focuses on a nonglacierized subcatchment of the Shullcas River Watershed and has 2 objectives: (a) create a model of the Shullcas groundwater system and assess the controls on stream discharge and (b) investigate the impact of the infiltration trenches on recharge and baseflow. We first collected hydrologic data from the field including a year‐long hydrograph (2015–2016), meteorological data (2011–2016), and infiltration measurements. We use a recharge model to evaluate the impact of trenched hillslopes on infiltration and runoff processes. Finally, we use a 3‐dimensional groundwater model, calibrated to the measured dry season baseflow, to determine the impact of trenching on the catchment. Simulations show that trenched hillslopes receive approximately 3.5% more recharge, relative to precipitation, compared with unaltered hillslopes. The groundwater model indicates that because the groundwater flow system is fast and shallow, incorporating trenched hillslopes (~2% of study subcatchment area) only slightly increases baseflow in the dry season. Furthermore, the location of trenching is an important consideration: Trenching higher in the catchment (further from the river) and in flatter terrain provides more baseflow during the dry season. The results of this study may have important implications for Andean landscape management and water resources.     

Bidirectional stream–groundwater flow in response to ephemeral and intermittent streamflow and groundwater seasonality

It is often assumed that the net groundwater flow direction is towards the channel in headwater streams in humid climates, with magnitudes dependent on flow state. However, studies that characterize stream–groundwater interactions in ephemeral and intermittent streams in humid landscapes remain sparse. Here, we examined seasonally driven stream–groundwater interactions in response to temporary streamflow on the basis of field observations of streamflow and groundwater on an adjacent hillslope. The direction of hydraulic head gradients between the stream and groundwater shifted seasonally. The stream gained water (head gradients were towards the stream) when storage state was high. During this period, streamflow was persistent. The stream lost water to the groundwater system (head gradients were away from the stream) when storage state was low. During this period, streamflow only occurred in response to precipitation events, and head gradients remained predominantly away from the stream during events. This suggested that mechanisms other than deep groundwater contributions produced run‐off when storage was low, such as surface and perched subsurface flowpaths above the water table. Analysis of the annual water balance for the study period showed that the residual between precipitation inputs and streamflow and evapotranspiration outputs, which were attributed to the loss of water to the deeper, regional groundwater system, was similar in magnitude to streamflow. This, coupled with results that showed bidirectionality in stream–groundwater head gradients, indicated that headwaters composed of temporary (e.g., ephemeral and intermittent) streams can be important focal areas for regional groundwater recharge, and both contribute to and receive water, solutes, and materials from the groundwater system.     

Air/water/sediment temperature contrasts in small streams to identify groundwater seepage locations

The need to identify groundwater seepage locations is of great importance for managing both stream water quality and groundwater sourced ecosystems due to their dependency on groundwater‐borne nutrients and temperatures. Although several reconnaissance methods using temperature as tracer exist, these are subjected to limitations related to mainly the spatial and temporal resolution and/or mixing of groundwater and surface water leading to dilution of the temperature differences. Further, some methods, for example, thermal imagery and fiber optic distributed temperature sensing, although relative efficient in detecting temperature differences over larger distances, these are labor‐intensive and costly. Therefore, there is a need for additional cost‐effective methods identifying substantial groundwater seepage locations. We present a method expanding the linear regression of air and stream temperatures by measuring the temperatures in dual‐depth; in the stream column and at the streambed‐water interface (SWI). By doing so, we apply metrics from linear regression analysis of temperatures between air/stream and air/SWI (linear regression slope, intercept, and coefficient of determination), and the daily water temperature cycle (daily mean temperatures, temperature variance, and the mean diel temperature fluctuation). We show that using metrics from only single‐depth stream temperature measurements are insufficient to identify substantial groundwater seepage locations in a head‐water stream. Conversely, comparing the metrics from dual‐depth temperatures show significant differences; at groundwater seepage locations, temperatures at the SWI merely explain 43–75% of the variation opposed to ? 91% at the corresponding stream column temperatures. In general, at these locations at the SWI, the slopes ( < 0.25) and intercepts ( > 6.5 °C) are substantially lower and higher, respectively, while the mean diel temperature fluctuations ( < 0.98 °C) are decreased compared to remaining locations. The dual‐depth approach was applied in a post‐glacial fluvial setting, where metrics analyses overall corroborated with field measurements of groundwater fluxes and stream flow accretions. Thus, we propose a method reliably identifying groundwater seepage locations along streambeds in such settings.     

Stream and bed temperature variability in a coastal headwater catchment: influences of surface‐subsurface interactions and partial‐retention forest harvesting

Stream temperature was recorded between 2002 and 2005 at four sites in a coastal headwater catchment in British Columbia, Canada. Shallow groundwater temperatures, along with bed temperature profiles at depths of 1 to 30 cm, were recorded at 10‐min intervals in two hydrologically distinct reaches beginning in 2003 or 2004, depending on the site. The lower reach had smaller discharge contributions via lateral inflow from the hillslopes and fewer areas with upwelling (UW) and/or neutral flow across the stream bed compared to the middle reach. Bed temperatures were greater than those of shallow groundwater during summer, with higher temperatures in areas of downwelling (DW) flow compared to areas of neutral and UW flow. A paired‐catchment analysis revealed that partial‐retention forest harvesting in autumn 2004 resulted in higher daily maximum stream and bed temperatures but smaller changes in daily minima. Changes in daily maximum stream temperature, averaged over July and August of the post‐harvest year, ranged from 1.6 to 3 °C at different locations within the cut block. Post‐harvest changes in bed temperature in the lower reach were smaller than the changes in stream temperature, greater at sites with DW flow, and decreased with depth at both UW and DW sites, dropping to about 1 °C at a depth of 30 cm. In the middle reach, changes in daily maximum bed temperature, averaged over July and August, were generally about 1 °C and did not vary significantly with depth. The pre‐harvest regression models for shallow groundwater were not suitable for applying the paired‐catchment analysis to estimate the effects of harvesting. However, shallow groundwater was warmer at the lower reach following harvesting, despite generally cooler weather compared to the pre‐harvest year. Copyright © 2012 John Wiley & Sons, Ltd.     

Floodplain connection buffers seasonal changes in urban stream water quality

Urban streams in the Northeastern United States have large road salt inputs during the winter, increased nonpoint sources of inorganic nitrogen and decreased short‐term and permanent storage of nutrients. Restoration activities that re‐establish connection between streams and riparian environments may be effective for improving urban stream water quality. Meadowbrook Creek, a first‐order stream in Syracuse, NY, provides a unique setting to explore impacts of stream–floodplain connection because it flows along a negative urbanization gradient, from channelized and armoured headwaters to a broad, vegetated floodplain with a riparian aquifer. In this study, we investigated how reconnection to groundwater and introduction of riparian vegetation impacted urban surface water chemistry by making biweekly longitudinal surveys of stream water chemistry in the creek from May 2012 until June 2013. We used multiple methods to measure groundwater discharge rates along the creek. Chloride concentrations in the upstream, disconnected reach were influenced by discharge of road salt during snow melt events and ranged from 161.2 to 1440 mg/l. Chloride concentrations in the downstream, connected reach had less temporal variation, ranging from 252.0 to 1049 mg/l, because of buffering by groundwater discharge, as groundwater chloride concentrations ranged from 84.0 to 655.4 mg/l. In the summer, there was little to no nitrate in the disconnected reach because of limited sources and high primary productivity, but concentrations reached over 1 mg N/l in the connected reach because of the presence of riparian vegetation. During the winter, when temperatures fell below freezing, nitrate concentrations in the disconnected reach increased to 0.58 mg N/l but were still lower than the connected reach, which averaged 0.88 mg N/l. Urban stream restoration projects that restore floodplain connection may impact water quality by storing high salinity road run‐off during winter overbank events and discharging that water year‐round, thereby attenuating seasonal fluctuations in chloride. Contrary to prior findings, we observed that floodplain connection and riparian vegetation may alter nitrate sources and sinks such that nitrate concentrations increase longitudinally in connected urban streams. Copyright © 2014 John Wiley & Sons, Ltd.     

Flow directions of stream-groundwater exchange in a headwater catchment during the hydrologic year

Understanding near-stream groundwater dynamics and flow directions is important for predicting hillslope-stream connectivity, streamflow generation, and hydrologic controls of streamwater quality. To determine the drivers of groundwater flow in the stream corridor (i.e., the stream channel and the adjacent groundwater in footslopes and riparian areas), we observed the water levels of 36 wells and 7 piezometers along a headwater stream section over a period of 18 months. Groundwater dynamics during events were controlled by the initial position of the groundwater table relative to the subsurface structure. The near-stream groundwater table displayed a fast and pronounced response to precipitation events when lying in fractured bedrock with low storage capacity, and responded less frequently and in a less pronounced way when lying in upper layers with high storage capacity. Precipitation depth, intensity, regolith thickness above the fractured bedrock, and proximity to and elevation above the stream channel also had an effect on the groundwater dynamics, which varied with hydrologic conditions. Our high-frequency and spatially dense measurements highlight the competing influence of groundwater inflow from upslope locations, streamwater level and bedrock properties on the spatiotemporal dynamics of flowpaths in the stream corridor. Near-stream groundwater pointed uniformly towards the stream channel when the stream corridor was hydrologically connected to upslope groundwater. However, local interruptions of the water inflow from upslope locations caused flow reversals towards the footslopes. The direction of near-stream groundwater followed the local fractured bedrock topography during dry hydrologic conditions on a few occasions after events. The outcomes of this research contribute to a better understanding of the drivers controlling spatiotemporal changes in near-stream groundwater dynamics and flow directions in multiple wetness states of the stream corridor.     

Quantifying groundwater discharge to a small perennial stream in southern Ontario, Canada

A study of the interaction between groundwater and surface water was undertaken within a small agricultural watershed in southern Ontario, Canada. Groundwater contributions to streamflow were measured along a section of stream during baseflow conditions and during rainfall events. Four techniques were used to estimate the contribution of groundwater to the stream along a 450 m reach (three during baseflow and one during stormflow conditions). Under baseflow conditions, streamflow measurements using the velocity–area technique indicated that the net groundwater flux to the stream during the summer months was 10 ml s⁻¹ m⁻¹. Hydrometric measurements (i.e. hydraulic gradient and hydraulic conductivity) taken using mini-piezometers installed in the sediments beneath the stream resulted in net groundwater flux estimates that were four to five times lower. Seepage meters failed to provide any measurements of water flux into or out of the stream. Therefore, based on these results, the velocity–area technique gives the best estimate of groundwater discharge. Hydrograph separations were conducted using isotopic ratios and electrical conductivity on two large rainfall events with different antecedent moisture conditions in the catchment. Both events showed that pre-event water (generally considered groundwater) dominated streamflow and tile drain flow with 64%–80% of the total discharge contributed by pre-event water. High water table conditions within the catchment resulted in greater stream discharge and a greater contribution of event water in the streamflow than that observed under low water table conditions for similar intensity storm events. The results also showed that differences in riparian zone width, vegetation and surface saturation conditions between the upper and lower catchment can influence the relative magnitude of streamflow response from the two catchment areas.     

Discharge estimation of submarine springs in the Dead Sea based on velocity or density measurements in proximity to the water surface

The Dead Sea is a closed lake, the water level of which is lowering at an alarming rate of about 1 m/year. Factors difficult to determine in its water balance are evaporation and groundwater inflow, some of which emanate as submarine groundwater discharge. A vertical buoyant jet generated by the difference in densities between the groundwater and the Dead Sea brine forms at submarine spring outlets. To characterize this flow field and to determine its volumetric discharge, a system was developed to measure the velocity and density of the ascending submarine groundwater across the center of the stream along several horizontal sections and equidistant depths while divers sampled the spring. This was also undertaken on an artificial submarine spring with a known discharge to determine the quality of the measurements and the accuracy of the method. The underwater widening of the flow is linear and independent of the volumetric spring discharge. The temperature of the Dead Sea brine at lower layers primarily determines the temperature of the surface of the upwelling, produced above the jet flow, as the origin of the main mass of water in the submarine jet flow is Dead Sea brine. Based on the measurements, a model is presented to evaluate the distribution of velocity and solute density in the flow field of an emanating buoyant jet. This model allows the calculation of the volumetric submarine discharge, merely requiring either the maximum flow velocity or the minimal density at a given depth.     

Groundwater dynamics in a till hillslope: flow directions,gradients and delay

Knowledge of groundwater dynamics is important for the understanding of hydrological controls on chemical processes along the water flow pathways. To increase our knowledge of groundwater dynamics in areas with shallow groundwater, the groundwater dynamics along a hillslope were studied in a boreal catchment in Southern Sweden. The forested hillslope had a 1‐ to 2‐m deep layer of sandy till above bedrock. The groundwater flow direction and slope were calculated under the assumption that the flow followed the slope of the groundwater table, which was computed for different triangles, each defined by three groundwater wells. The flow direction showed considerable variations over time, with a maximum variation of 75°. During periods of high groundwater levels the flow was almost perpendicular to the stream, but as the groundwater level fell, the flow direction became gradually more parallel to the stream, directed in the downstream direction. These findings are of importance for the interpretation of results from hillslope transects, where the flow direction usually is assumed to be invariable and always in the direction of the hillslope. The variations in the groundwater flow direction may also cause an apparent dispersion for groundwater‐based transport. In contrast to findings in several other studies, the groundwater level was most responsive to rainfall and snowmelt in the upper part of the hillslope, while the lower parts of the slope reached their highest groundwater level up to 40 h after the upper parts. This can be explained by the topography with a wetter hollow area in the upper part. Copyright © 2011 John Wiley & Sons, Ltd.     

Stream temperature dynamics in two hydrogeomorphically distinct reaches

The objective of this study was to analyse stream temperature variability during summer in relation to both surface heat exchanges and reach‐scale hydrology for two hydrogeomorphically distinct reaches. The study focused on a 1·5‐km wildfire‐disturbed reach of Fishtrap Creek located north of Kamloops, British Columbia. Streamflow measurements and longitudinal surveys of electrical conductivity and water chemistry indicated that the upper 750 m of the study reach was dominated by flow losses. A spring discharged into the stream at 750 m below the upper reach boundary. Below the spring, the stream was neutral to losing on three measurement days, but gained flow on a fourth day that followed a rain event. Continuous stream temperature measurements typically revealed a downstream warming along the upper 750 m of the study reach on summer days, followed by a pronounced cooling associated with the spring, with little downstream change below the spring. Modelled surface energy exchanges were similar over the upper and lower sub‐reaches, and thus cannot explain the differences in longitudinal temperature patterns. Application of a Lagrangian stream temperature model provided reasonably accurate predictions for the upper sub‐reach. For the lower sub‐reach, accurate prediction required specification of concurrent flow losses and gains as a hydrological boundary condition. These findings are consistent with differences in the hydrogeomorphology of the upper and lower sub‐reaches. The modelling exercise indicated that substantial errors in predicted stream temperature can occur by representing stream‐surface exchange as a reach‐averaged one‐directional flux computed from differences in streamflow between the upper and lower reach boundaries. Further research should focus on reliable methods for quantifying spatial variations in reach‐scale hydrology. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Mapping of road‐salt‐contaminated groundwater discharge and estimation of chloride load to a small stream in southern New Hampshire,USA

Concentrations of chloride in excess of State of New Hampshire water‐quality standards (230 mg/l) have been measured in watersheds adjacent to an interstate highway (I‐93) in southern New Hampshire. A proposed widening plan for I‐93 has raised concerns over further increases in chloride. As part of this effort, road‐salt‐contaminated groundwater discharge was mapped with terrain electrical conductivity (EC) electromagnetic (EM) methods in the fall of 2006 to identify potential sources of chloride during base‐flow conditions to a small stream, Policy Brook. Three different EM meters were used to measure different depths below the streambed (ranging from 0 to 3 m). Results from the three meters showed similar patterns and identified several reaches where high EC groundwater may have been discharging. Based on the delineation of high (up to 350 mmhos/m) apparent terrain EC, seven‐streambed piezometers were installed to sample shallow groundwater. Locations with high specific conductance in shallow groundwater (up to 2630 mmhos/m) generally matched locations with high streambed (shallow subsurface) terrain EC. A regression equation was used to convert the terrain EC of the streambed to an equivalent chloride concentration in shallow groundwater unique for this site. Utilizing the regression equation and estimates of one‐dimensional Darcian flow through the streambed, a maximum potential groundwater chloride load was estimated at 188 Mg of chloride per year. Changes in chloride concentration in stream water during streamflow recessions showed a linear response that indicates the dominant process affecting chloride is advective flow of chloride‐enriched groundwater discharge. Published in 2010 by John Wiley & Sons, Ltd.     

Impact of beaver dam analogues on hydrology in a semi-arid floodplain

Natural beaver ponds help connect the stream to the floodplain, maintain late summer low flows and reduce peak flow during high flow events by offering temporary surface water (SW) storage. When beavers are extirpated from the landscape, stream degradation often ensues. This study assesses the impact of beaver dam analogues (BDA) as a stream restoration technique to help maintain low flow water levels and enhance stream-floodplain interactions on a seasonal basis in Red Canyon Creek, Lander, WY. BDAs increased SW and groundwater (GW) levels, favoured the occurrence of flow reversals (i.e., stream-to-floodplain GW flow) during high flow events associated with mid-winter and early-spring thaw events, and reduced the groundwater-to-stream hydraulic gradient on an annual basis. Although GW temperatures varied seasonally, relatively cooler GW temperatures were observed in the BDA impacted reach compared to the control reach away from BDA influence. BDAs however did not significantly impact stream temperatures. Overall, results suggest that when installed in sequence, BDA complexes can successfully reconnect the stream to its floodplain, and ultimately increase SW-GW exchange at the floodplain scale by allowing flow reversals to occur and by reducing the GW to stream hydraulic gradient. Although BDAs built with fence posts, willow branches, sediments and small boulders are naturally porous and require regular maintenance, this study also highlights the viability of small BDAs as a restoration practice to enhance landscape resilience to drought and high flow events in deeply incised channels where beavers would not come back naturally.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.     

Influence of groundwater and topography on stream drying in semi-arid headwater streams

Non-perennial streams comprise over half of the global stream network and impact downstream water quality. Although aridity is a primary driver of stream drying globally, surface flow permanence varies spatially and temporally within many headwater streams, suggesting that these complex drying patterns may be driven by topographic and subsurface factors. Indeed, these factors affect shallow groundwater flows in perennial systems, but there has been only limited characterisation of shallow groundwater residence times and groundwater contributions to intermittent streams. Here, we asked how groundwater residence times, shallow groundwater contributions to streamflow, and topography interact to control stream drying in headwater streams. We evaluated this overarching question in eight semi-arid headwater catchments based on surface flow observations during the low-flow period, coupled with tracer-based groundwater residence times. For one headwater catchment, we analysed stream drying during the seasonal flow recession and rewetting period using a sensor network that was interspersed between groundwater monitoring locations, and linked drying patterns to groundwater inputs and topography. We found a poor relationship between groundwater residence times and flowing network extent (R² < 0.24). Although groundwater residence times indicated that old groundwater was present in all headwater streams, surface drying also occurred in each of them, suggesting old, deep flowpaths are insufficient to sustain surface flows. Indeed, the timing of stream drying at any given point typically coincided with a decrease in the contribution from near-surface sources and an increased relative contribution of groundwater to streamflow at that location, whereas the spatial pattern of drying within the stream network typically correlated with locations where groundwater inputs were most seasonally variable. Topographic metrics only explained ~30% of the variability in seasonal flow permanence, and surprisingly, we found no correlation with seasonal drying and down-valley subsurface storage area. Because we found complex spatial patterns, future studies should pair dense spatial observations of subsurface properties, such as hydraulic conductivity and transmissivity, to observations of seasonal flow permanence.