Temporal and spatial response of hyporheic zone geochemistry to a storm event

Although there has been recent focus on understanding spatial variability in hyporheic zone geochemistry across different morphological units under baseflow conditions, less attention has been paid to temporal responses of hyporheic zone geochemistry to non‐steady‐state conditions. We documented spatial and temporal variability of hyporheic zone geochemistry in response to a large‐scale storm event, Tropical Storm Irene (August 2011), across a pool–riffle–pool sequence along Chittenango Creek in Chittenango, NY, USA. We sampled stream water as well as pore water at 15 cm depth in the streambed at 14 locations across a 30 m reach. Sampling occurred seven times at daily intervals: once during baseflow conditions, once during the rising limb of the storm hydrograph, and five times during the receding limb. Principal component analysis was used to interpret temporal and spatial changes and dominant drivers in stream and pore water geochemistry (n = 111). Results show the majority of spatial variance in hyporheic geochemistry (62%) is driven by differential mixing of stream and ground water in the hyporheic zone. The second largest driver (17%) of hyporheic geochemistry was temporal dilution and enrichment of infiltrating stream water during the storm. Hyporheic sites minimally influenced by discharging groundwater (‘connected’ sites) showed temporal changes in water chemistry in response to the storm event. Connected sites within and upstream of the riffle reflected stream geochemistry throughout the storm, whereas downstream sites showed temporally lagged responses in some conservative and biogeochemically reactive solutes. This suggests temporal changes in hyporheic geochemistry at these locations reflect a combination of changes in infiltrating stream chemistry and hyporheic flowpath length and residence time. The portion of the study area strongly influenced by groundwater discharge increased in size throughout the storm, producing elevated Ca²⁺ and concentrations in the streambed, suggesting zones of localized groundwater inputs expand in response to storms. Copyright © 2013 John Wiley & Sons, Ltd.     

Impact of debris dams on hyporheic interaction along a semi‐arid stream

Hyporheic exchange increases the potential for solute retention in streams by slowing downstream transport and increasing solute contact with the substrate. Hyporheic exchange may be a major mechanism to remove nutrients in semi‐arid watersheds, where livestock have damaged stream riparian zones and contributed nutrients to stream channels. Debris dams, such as beaver dams and anthropogenic log dams, may increase hyporheic interactions by slowing stream water velocity, increasing flow complexity and diverting water to the subsurface. Here, we report the results of chloride tracer injection experiments done to evaluate hyporheic interaction along a 320 m reach of Red Canyon Creek, a second order stream in the semi‐arid Wind River Range of Wyoming. The study site is part of a rangeland watershed managed by The Nature Conservancy of Wyoming, and used as a hydrologic field site by the University of Missouri Branson Geologic Field Station. The creek reach we investigated has debris dams and tight meanders that hypothetically should enhance hyporheic interaction. Breakthrough curves of chloride measured during the field experiment were modelled with OTIS‐P, a one‐dimensional, surface‐water, solute‐transport model from which we extracted the storage exchange rate α and cross‐sectional area of the storage zone A_s for hyporheic exchange. Along gaining reaches of the stream reach, short‐term hyporheic interactions associated with debris dams were comparable to those associated with severe meanders. In contrast, along the non‐gaining reach, stream water was diverted to the subsurface by debris dams and captured by large‐scale near‐stream flow paths. Overall, hyporheic exchange rates along Red Canyon Creek during snowmelt recession equal or exceed exchange rates observed during baseflow at other streams. Copyright © 2005 John Wiley & Sons, Ltd.     

Effect of morphology and discharge on hyporheic exchange flows in two small streams in the Cascade Mountains of Oregon,USA

Stream‐tracer injections were used to examine the effect of channel morphology and changing stream discharge on hyporheic exchange flows. Direct observations were made from well networks to follow tracer movement through the hyporheic zone. The reach‐integrated influence of hyporheic exchange was evaluated using the transient storage model (TSM) OTIS‐P. Transient storage modelling results were compared with direct observations to evaluate the reliability of the TSM. Results from the tracer injection in the bedrock reach supported the assumption that most transient storage in headwater mountain streams results from hyporheic exchange. Direct observations from the well networks in colluvial reaches showed that subsurface flow paths tended to parallel the valley axis. Cross‐valley gradients were weak except near steps, where vertical and cross‐valley hydraulic gradients indicated a strong potential for stream water to downwell into the hyporheic zone. The TSM parameters showed that both size and residence time of transient storage were greater in reaches with a few large log‐jam‐formed steps than in reaches with more frequent, but smaller steps. Direct observations showed that residence times in the unconstrained stream were longer than in the constrained stream and that little change occurred in the location and extent of the hyporheic zone between low‐ and high‐baseflow discharges in any of the colluvial reaches. The transient storage modelling results did not agree with these observations, suggesting that the TSM was insensitive to long residence‐time exchange flows and was very sensitive to changes in discharge. Disagreements between direct observations and the transient storage modelling results highlight fundamental problems with the TSM that confound comparisons between the transient storage modelling results for tracer injections conducted under differing flow conditions. Overall, the results showed that hyporheic exchange was little affected by stream discharge (at least over the range of baseflow discharges examined in this study). The results did show that channel morphology controlled development of the hyporheic zone in these steep mountain stream channels. Copyright © 2005 John Wiley & Sons, Ltd.     

Temporal Hyporheic Zone Response to Water Table Fluctuations

Expansion and contraction of the hyporheic zone due to temporal hydrologic changes between stream and riparian aquifer influence the biogeochemical cycling capacity of streams. Theoretical studies have quantified the control of groundwater discharge on the depth of the hyporheic zone; however, observations of temporal groundwater controls are limited. In this study, we develop the concept of groundwater‐dominated differential hyporheic zone expansion to explain the temporal control of groundwater discharge on the hyporheic zone in a third‐order stream reach flowing through glacially derived terrain typical of the Great Lakes region. We define groundwater‐dominated differential expansion of the hyporheic zone as: differing rates and magnitudes of hyporheic zone expansion in response to seasonal vs. storm‐related water table fluctuation. Specific conductance and vertical hydraulic gradient measurements were used to map changes in the hyporheic zone during seasonal water table decline and storm events. Planar and riffle beds were monitored in order to distinguish the cause of increasing hyporheic zone depth. Planar bed seasonal expansion of the hyporheic zone was of a greater magnitude and longer in duration (weeks to months) than storm event expansion (hours to days). In contrast, the hyporheic zone beneath the riffle bed exhibited minimal expansion in response to seasonal groundwater decline compared to storm related expansion. Results indicated that fluctuation in the riparian water table controlled seasonal expansion of the hyporheic zone along the planar bed. This groundwater induced hyporheic zone expansion could increase the potential for biogeochemical cycling and natural attenuation.     

The role of the hyporheic zone across stream networks

Many hyporheic papers state that the hyporheic zone is a critical component of stream ecosystems, and many of these papers focus on the biogeochemical effects of the hyporheic zone on stream solute loads. However, efforts to show such relationships have proven elusive, prompting several questions: Are the effects of the hyporheic zone on stream ecosystems so highly variable in place and time (or among streams) that a consistent relationship should not be expected? Or, is the hyporheic zone less important in stream ecosystems than is commonly expected? These questions were examined using data from existing groundwater modelling studies of hyporheic exchange flow at five sites in a fifth‐order, mountainous stream network. The size of exchange flows, relative to stream discharge (Q_HEF:Q), was large only in very small streams at low discharge (area ≈ 100 ha; Q < 10 l/s). At higher flows (flow exceedance probability > 0·7) and in all larger streams, Q_HEF:Q was small. These data show that biogeochemical processes in the hyporheic zone of small streams can substantially influence the stream's solute load, but these processes become hydrologically constrained at high discharge or in larger streams and rivers. The hyporheic zone may influence stream ecosystems in many ways, however, not just through biogeochemical processes that alter stream solute loads. For example, the hyporheic zone represents a unique habitat for some organisms, with patterns and amounts of upwelling and downwelling water determining the underlying physiochemical environment of the hyporheic zone. Similarly, hyporheic exchange creates distinct patches of downwelling and upwelling. Upwelling environments are of special interest, because upwelling water has the potential to be thermally or chemically distinct from stream water. Consequently, micro‐environmental patches created by hyporheic exchange flows are likely to be important to biological and ecosystem processes, even if their impact on stream solute loads is small. Published in 2011 by John Wiley & Sons, Ltd.     

The effect of transient storage on nitrate uptake lengths in streams: an inter‐site comparison

Surface water–groundwater interaction in the hyporheic zone may enhance biogeochemical cycling in streams, and it has been hypothesized that streams exchanging more water with the hyporheic zone should have more rapid nitrate utilization. We used simultaneous conservative solute and nitrate addition tracer tests to measure transient storage (which includes hyporheic exchange and in‐stream storage) and the rate of nitrate uptake along three reaches within the Red Canyon Creek watershed, Wyoming. We calibrated a one‐dimensional transport model, incorporating transient storage (OTIS‐P), to the conservative solute breakthrough curves and used the results to determine the degree of transient storage in each reach. The nitrate uptake length was quantified from the exponential decrease in nitrate concentration with distance during the tracer tests. Nitrate uptake along the most downstream reach of Red Canyon Creek was rapid (turnover time K^?1_c = 32 min), compared with nitrate uptake reported in other studies (K^?1_c = 12 to 551 min), but other sites within the watershed showed little nitrate retention or loss. The uptake length S_w‐NO^?₃ for the most downstream reach was 500 m and the mass transfer coefficient V_f‐NO^?₃ was 6·3 m min^?1. Results from 15 other nitrate‐addition tracer tests were used to create a regression model relating transient storage and measures of stream flow to nitrate uptake length. The model, which includes specific discharge and transient storage area, explains almost half the variability in nitrate uptake length (adjusted R² = 0·44) and is most effective for comparing sites with very different stream characteristics. Although large differences in specific discharge and storage zone area explain inter‐site differences in nitrate uptake, other unmeasured variables, such as available organic carbon and microbial community composition, are likely important for predicting differences in nitrate uptake between sites with similar specific discharge rates and storage zone areas, such as when making intra‐site comparisons. Copyright © 2007 John Wiley & Sons, Ltd.     

Identifying and linking source areas of flow and P transport in dairy‐grazed headwater catchments,North Island,New Zealand

We examined the applicability of the critical‐source area (CSA) concept to the dairy‐grazed 192‐ha Upper Toenepi catchment and its 8·7‐ha Kiwitahi sub‐catchment, New Zealand. We evaluated if phosphorus (P) transport from land into stream is dominated by saturation‐excess (SE) and infiltration‐excess (IE) runoff during stormflow and by sub‐surface (<1·5 m depth) flows during baseflow. We measured stream flow and shallow groundwater levels, collected monthly stream, tile drain (TDA) and groundwater samples, and flow‐proportional stream samples from the Kiwitahi sub‐catchment, and determined their dissolved reactive phosphorus (DRP) and total phosphorus (TP) concentrations. In the Kiwitahi sub‐catchment, during storm events, IE contributions were significant. Contributions from SE appeared significant in the Upper Toenepi catchment. However, in both catchments, sub‐surface contributions dominated stormflow and baseflow periods. Absence of water table at the surface and the water table gradient towards the stream indicated that P transport during events was not limited to surface runoff. The dynamics of the groundwater table and the occurrence of SE areas were influenced by proximity to the stream and hillslope positions. Baseflow accounted for 42% of the annual flow in the Kiwitahi sub‐catchment, and contributed 37 and 52% to the DRP and TP loads, respectively. The P transport during baseflow appeared equally important as P losses from CSAs during stormflow. The close resemblance in P levels between groundwater and stream samples during baseflow demonstrates the importance of shallow groundwater for stream flow. In the Upper Toenepi catchment, contributions from effluent ponds (EFFs) dominated P loads. Management strategies should focus on controlling P release from EFFs, and on decreasing Olsen P concentrations in soil to minimize leaching of P via sub‐surface flow to streams. Research is needed to quantify the role of sub‐surface flow as well as to expand management strategies to minimize P transfers during stormflow and baseflow conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

Implementing a nonlinear groundwater module in the soil and water assessment tool (SWAT)

The Soil and Water Assessment Tool (SWAT) is widely used in modeling water quantity and quality. In the original SWAT, groundwater flow is calculated using a linear‐reservoir model, with outflow proportional to storage. However, observations show that this assumption is not always applicable; for example, macropores in Karst formations would seriously affect the groundwater behavior. A nonlinear groundwater algorithm was introduced in a new version of the SWAT model, called ISWAT. The Shenandoah Valley area in the Eastern U.S., which includes a number of geologic formations including Karst, was selected to test the modified ISWAT model. Parameter ESTimation (PEST) was coupled with ISWAT to auto‐calibrate the nonlinear parameter values. Ten years of record at 15 stream gauges were used to calibrate the model. The nonlinear ISWAT, statistically and visually, performed better in stream discharge estimation especially during baseflow recession and low‐flow periods. This indicated that the nonlinear algorithm can better represent groundwater behavior. The coupled ISWAT‐PEST approach can be used in future stream discharge simulation. Copyright © 2013 John Wiley & Sons, Ltd.     

Impact of redox and transport processes in a riparian wetland on stream water quality in the Fichtelgebirge region,southern Germany

Biologically mediated redox processes in the riparian zone, like denitrification, can have substantially beneficial impacts on stream water quality. The extent of these effects, however, depends greatly on the hydrological boundary conditions. The impact of hydrological processes on a wetland's nitrogen sink capacity was investigated in a forested riparian fen which is drained by a first‐order perennial stream. Here, we analysed the frequency distributions and time‐series of pH and nitrogen, silica, organic carbon and oxygen concentrations in throughfall, soil solution, groundwater and stream water, and the groundwater levels and stream discharges from a 3‐year period. During baseflow conditions, the stream was fed by discharging shallow, anoxic groundwater and by deep, oxic groundwater. Whereas the latter delivered considerable amounts of nitrogen (～0·37 mg l^?1) to the stream, the former was almost entirely depleted of nitrogen. During stormflow, near‐surface runoff in the upper 30 cm soil layer bypassed the denitrifying zone and added significant amounts to the nitrogen load of the stream. Nitrate‐nitrogen was close to 100% of deep groundwater and stream‐water nitrogen concentration. Stream‐water baseflow concentrations of nitrate, dissolved carbon and silica were about 1·6 mg l^?1, 4 mg l^?1 and 7·5 mg l^?1 respectively, and >3 mg l^?1, >10 mg l^?1 and <4 mg l^?1 respectively during discharge peaks. In addition to that macroscale bypassing effect, there was evidence for a corresponding microscale effect: Shallow groundwater sampled by soil suction cups indicated complete denitrification and lacked any seasonal signal of solute concentration, which was in contrast to piezometer samples from the same depth. Moreover, mean solute concentration in the piezometer samples resembled more that of suction‐cup samples from shallower depth than that of the same depth. We conclude that the soil solution cups sampled to a large extent the immobile soil‐water fraction. In contrast, the mobile fraction that was sampled by the piezometers exhibited substantially shorter residence time, thus being less exposed to denitrification, but predominating discharge of that layer to the stream. Consequently, assessing the nitrogen budget based on suction‐cup data tended to overestimate the nitrogen consumption in the riparian wetland. These effects are likely to become more important with the increased frequency and intensity of rainstorms that are expected due to climate change. Copyright © 2006 John Wiley & Sons, Ltd.     

Quantifying nitrogen cycling beneath a meander of a low gradient,N‐impacted,agricultural stream using tracers and numerical modelling

In watersheds impacted by nitrate from agricultural fertilizers, nitrification and denitrification may be decoupled as denitrification in the hyporheic zone is not limited to naturally produced nitrate. While most hyporheic research focuses on the 1–2 m of sediment beneath the stream bed, there are a limited number of studies that quantify nitrogen (N) cycling at larger hyporheic scales (10s of metres to kms). We conducted an investigation to quantify N cycling through a single meander of a low gradient, meandering stream, draining an agricultural watershed. Chemistry (major ions and N species) and hydrologic data were collected from the stream and groundwater beneath the meander. Evidence indicates that nearly all the shallow groundwater flowing beneath the meander originates as stream water on the upgradient side of the meander, and returns to the stream on the downgradient side. We quantified the flux of water beneath the meander using a numerical model. The flux of N into and out of the meander was quantified by multiplying the concentration of the important N species (nitrate, ammonium, dissolved organic nitrogen (DON)) by the modelled water fluxes. The flux of N into the meander is dominated by nitrate, and the flux of N out of the meander is dominated by ammonium and DON. While stream nitrate varied seasonally, ammonium and DON beneath the meander were relatively constant throughout the year. When stream nitrate concentrations are high (>2 mg litre^?1), flow beneath the meander is a net sink for N as more N from nitrate in stream water is consumed than is produced as ammonium and DON. When stream nitrate concentrations are low (<2 mg litre^?1), the flux of N entering is less than exiting the meander. On an annual basis, the meander hyporheic flow serves as a net sink for N. Copyright © 2007 John Wiley & Sons, Ltd.     

Unsaturated zone fracture flow contributions to stream flow: evidence for the process in South Africa and its importance

Understanding hydrological processes has always been important to the development and successful application of conceptual hydrological models. It can also contribute to informed water resources management, particularly in the context of understanding the potential impacts of both land use and climate change. Improved conceptual and quantitative understanding of near‐surface hydrological processes emerged through field studies during the 1960s to1980s; however, there remains a degree of ambiguity about the processes that link surface water and groundwater. This is especially the case in South Africa where a great deal of confusion has arisen about the source of the ‘baseflow’ signal in stream flow observations. This paper suggests that fracture flow within the unsaturated zone could have a lateral component and therefore re‐emerge and contribute to stream flow in catchments with relatively steep topography. The implication is that ‘baseflows’ could be made up of groundwater contributions (caused by intersection of the water table with stream channels) as well as an unsaturated zone flow component. Evidence for the existence of the process is presented on the basis of small‐scale observations and interpretations of stream flow observations. The potential importance of the process relates to interpreting different methods of recharge estimation, assessing the impacts of groundwater abstraction on stream flow, as well as the application and interpretation of the results of hydrological models. The conclusions are that the process does exist, but that there is less than conclusive evidence for its importance. There is therefore a need for further studies that can quantify the scale of the process and therefore its importance. Only then will it be possible to develop a consistent understanding of the processes of surface water and groundwater interaction and therefore manage water resources in a truly integrated manner. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Collocation of hydrological and biological attenuation of nitrate in an urban stream

The hydrologic and biogeochemical processes that control nutrient export in urban streams are not well understood. Attenuation can occur by tributary dilution, groundwater discharge, and biological processing both in the water column and the hyporheic zone. A wastewater treatment plant on Pennypack Creek, an urban stream near Philadelphia, PA, provided high nitrate concentrations for analysis of downstream attenuation processes. Longitudinal sampling for an 8‐km reach revealed decreases in nitrate concentration of 2 mg l^?1 at high flow and 4.5 mg l^?1 during low flow. During high flow, δ¹⁵N‐NO₃ increased from 9.5 to 10.5‰ and during low flow increased from 10.1 to 11.1‰. Two reaches were sampled at fine spatial intervals (approximately 200 m) to better identify attenuation processes. Mixing analysis indicated that groundwater discharge and biological processing both control nitrate concentration and isotope signatures. However, fine‐scaled sampling did not reveal spatially discrete zones; instead, these processes were occurring simultaneously. While both processes attenuate nitrate, they have opposite isotope signatures, which may have muted changes in δ¹⁵N‐NO₃. At high flow, a decrease in Cl/NO₃ ratios helped distinguish groundwater discharge occurring along both finely sampled reaches. At low flow, biological processing seemed to be occurring more extensively, but the δ¹⁵N‐NO₃ signature was not consistent with either a single process or a sequential combination of groundwater dilution and biological nitrate attenuation. The collocation of processes makes it more difficult to assess biological processing hot spots and predict how urbanization and subsequent stream restoration influence nitrate attenuation. Copyright © 2016 John Wiley & Sons, Ltd.     

Riparian influence on hyporheic‐zone formation downstream of a small dam in the Blackland Prairie region of Texas

Small‐order streams have highly variable flows that can result in large temporal and spatial variation of the hyporheic zone. Dam construction along these intermittent headwater streams alters downstream flow and influences the hydrologic balance between stream water and the adjacent riparian zone. A 3‐year site study was conducted along an impounded second‐order stream to determine the water balance between stream, unsaturated zone, groundwater and riparian vegetation. The presence of the upstream impoundment provided near‐perennial water flow in the stream channel. The observed woody plant transpiration accounted for 71% of average annual water loss in the site. The overall contribution of stream water via the hyporheic zone to site water balance was 73 cm, or 44% of total inputs. This exceeded both rainfall and upland subsurface contribution to the site. A highly dynamic hyporheic zone was indicated by high water use from woody plants that fluctuated seasonally with stream water levels. We found leaf area development in the canopy layer to be closely coupled with stream and groundwater fluctuations, indicating its usefulness as a potential indicator of site water balance for small dam systems. The net result of upstream impoundment increased riparian vegetation productivity by influencing movement of stream water to storage in the groundwater system. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Mixing interfaces,fluxes, residence times and redox conditions of the hyporheic zones induced by dune-like bedforms and ambient groundwater flow

Recent studies highlighted the importance of the interface between streams and their surrounding sediment, known as the hyporheic zone, where stream waters flow through the alluvium. These pore water fluxes stem from the interaction among streambed morphology, stream hydraulics and surrounding groundwater flow. We analytically model the hyporheic hydraulics induced by a spatially uniform ambient groundwater flow made of a horizontal, underflow, and a vertical, basal, component, which mimics gaining and losing stream conditions. The proposed analytical solution allows to investigate the control of simple hydromorphological quantities on the extent, residence time and redox conditions of the hyporheic zone, and the thickness of the mixing interface between hyporheic and groundwater cells. Our analysis shows that the location of the mixing zone shallows or deepens in the sediment as a function of bedform geometry, surface hydraulic and groundwater flow. The point of stagnation, where hyporheic flow velocities vanish and where the separation surface passes through, is shallower than or coincides with the deepest point of the hyporheic zone only due to underflow. An increase of the ambient flow causes a reduction of the hyporheic zone volume similarly in both losing and gaining conditions. The hyporheic residence time is lognormally distributed under neutral, losing and gaining conditions, with the residence time moments depending on the same set of parameters describing dune morphology and stream flow.     

Control of river stage on the reactive chemistry of the hyporheic zone

We examined the influence of river stage on subsurface hydrology and pore water chemistry within the hyporheic zone of a groundwater‐fed river during the summer baseflow period of 2011. We found river stage and geomorphologic environment to control chemical patterns in the hyporheic zone. At a high river stage, the flux of upwelling water in the shallow sediments (>20 cm) decreased at sample sites in the upper section of our study reach and increased substantially at sites in the lower section. This differential response is attributed to the contrasting geomorphology of these subreaches that affects the rate of the rise and fall of a river stage relative to the subsurface head. At sites where streamward vertical flux decreased, concentration profiles of a conservative environmental tracer suggest surface water infiltration into the riverbed below depths recorded at a low river stage. An increase in vertical flux at sites in the lower subreach is attributed to the movement of lateral subsurface waters originating from the adjacent floodplain. This lateral‐moving water preserved or decreased the vertical extent of the hyporheic mixing zone observed at a low river stage. Downwelling surface water appeared to be responsible for elevated dissolved organic carbon (DOC) and manganese (Mn) concentrations in shallow sediments (0–20 cm); however, lateral subsurface flows were probably important for elevated concentrations of these solutes at deeper levels. Results suggest that DOC delivered to hyporheic sediments during a high river stage from surface water and lateral subsurface sources could enhance heterotrophic microbial activities. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessing the Transport of Pharmaceutical Compounds in a Layered Aquifer Discharging to a Stream

A groundwater plume containing high concentrations of pharmaceutical compounds, mainly sulfonamides, barbiturates, and ethyl urethane, in addition to chlorinated ethenes and benzene was investigated. The contamination originating from a former pharmaceutical industry discharges into a multilayered aquifer system and a downgradient stream. In this study, geological and hydrogeological data were integrated into a numerical flow model to examine identified trends using statistical approaches, including principal component analysis and hierarchal cluster analysis. A joint interpretation of the groundwater flow paths and contaminant concentrations in the different compartments (i.e., groundwater and hyporheic zone) provided insight on the transport processes of the different contaminant plumes to the stream. The analysis of historical groundwater concentrations of pharmaceutical compounds at the site suggested these compounds are slowly degrading. The pharmaceutical compounds migrate in both a deep semiconfined aquifer, as well as in the shallow unconfined aquifer, and enter the stream along a 2-km stretch. This contrasted with the chlorinated ethenes, which mainly discharge to the stream as a focused plume from the unconfined aquifer. The integrated approach developed here, combining groundwater flow modeling and statistical analyses of the contaminant concentration data collected in groundwater and the hyporheic zone, lead to an improved understanding of the observed distribution of contaminants in the unconfined and semiconfined aquifers, and thus to their discharge to the stream. This approach is particularly relevant for large and long-lasting contaminant sources and plumes, such as abandoned landfills and industrial production sites, where field investigations may be very expensive.     

Hyporheic and total transient storage in small,sand‐bed streams

Key processes in stream ecosystems are linked to hydraulic retention, which is the departure of stream flow from ideal ‘plug flow’, and reflects fluid movement through surface and hyporheic storage zones. Most existing information about hyporheic exchange is based on flume studies or field measurements in relatively steep streams with beds coarser than sand. Stream tracer studies may be used to quantify overall hydraulic retention, but disaggregation of surface and hyporheic retention remains difficult. A stream tracer approach was used to compute the rates at which stream water is exchanged with water in storage zones (total storage) in short reaches of two small, sand‐bed streams under free and obstructed flow conditions. Tracer curves were fit to the one‐dimensional transport with inflow storage model OTIS‐P. Networks of piezometers were used to measure specific discharge between the stream and the groundwater. In the sand‐bed streams studied, parameters describing total retention were in the upper 50% of data compiled from the literature, most of which represented streams with beds coarser than sand. However, hyporheic storage was an insignificant component of total hydraulic retention, representing only 0·01–0·49% of total exchange, and this fraction did not increase after installation of flow obstructions. Total retention did not vary systematically with bed material size, but increased 50–100% following flow obstruction. Removal of roughness elements, such as large wood and debris dams, is detrimental to processes dependent upon transient storage in small, sand‐bed streams. Copyright © 2007 John Wiley & Sons, Ltd.     

Potential mechanistic causes of increased baseflow across northern Eurasia catchments underlain by permafrost

Warming in the Arctic is occurring at twice the rate of the global average, resulting in permafrost thaw and a restructuring of the Arctic hydrologic cycle as indicated by increased stream discharge during low-flow periods. In these cold regions, permafrost thaw is postulated to increase low-flow discharge, or baseflow, through either: (a) localized increases in groundwater storage and discharge to streams due to increased aquifer transmissivity from thickening of the freeze–thaw layer above permafrost known as the active layer or (b) long-term increases in regional groundwater circulation via enhancement of groundwater–surface water interactions due to extensive permafrost loss over decades. While increasing baseflow has been observed throughout northern Eurasia, the precise mechanistic causes remain elusive. In this study, we differentiate between where these two subsurface physical mechanisms of baseflow increase are occurring by performing a baseflow recession analysis using daily streamflow records from 1913 to 2003 for 139 stations in northern Eurasia underlain by varying permafrost areal extents. Results indicate that from 1913 to 2003, the majority of catchments underlain by continuous permafrost have an increasing trend in their recession flow intercepts, a proxy for increasing active layer thickness. Alternatively, the majority of catchments underlain by permafrost types that are less spatially extensive (e.g., discontinuous, sporadic, isolated, or no permafrost) have decreasing trends in their recession flow intercepts, indicating that a potential increase in active layer thickness is not the driving factor of baseflow variations in these catchments. This may indicate that in catchments underlain by continuous permafrost, active layer thickening correlates with increases in baseflow, whereas, in other catchments with less extensive permafrost, increases in baseflow may be caused by wholesale permafrost loss and vertical talik expansion that enhances regional groundwater circulation. The results of this work may inform our understanding of the subsurface mechanisms responsible for the changing Arctic hydrologic cycle.