Determining long time‐scale hyporheic zone flow paths in Antarctic streams

In the McMurdo Dry Valleys of Antarctica, glaciers are the source of meltwater during the austral summer, and the streams and adjacent hyporheic zones constitute the entire physical watershed; there are no hillslope processes in these systems. Hyporheic zones can extend several metres from each side of the stream, and are up to 70 cm deep, corresponding to a lateral cross‐section as large as 12 m², and water resides in the subsurface year around. In this study, we differentiate between the near‐stream hyporheic zone, which can be characterized with stream tracer experiments, and the extended hyporheic zone, which has a longer time‐scale of exchange. We sampled stream water from Green Creek and from the adjacent saturated alluvium for stable isotopes of D and ¹⁸O to assess the significance and extent of stream‐water exchange between the streams and extended hyporheic zones over long time‐scales (days to weeks). Our results show that water residing in the extended hyporheic zone is much more isotopically enriched (up to 11‰ D and 2·2‰ ¹⁸O) than stream water. This result suggests a long residence time within the extended hyporheic zone, during which fractionation has occurred owing to summer evaporation and winter sublimation of hyporheic water. We found less enriched water in the extended hyporheic zone later in the flow season, suggesting that stream water may be exchanged into and out of this zone, on the time‐scale of weeks to months. The transient storage model OTIS was used to characterize the exchange of stream water with the extended hyporheic zone. Model results yield exchange rates (α) generally an order magnitude lower (10^?5 s^?1) than those determined using stream‐tracer techniques on the same stream. In light of previous studies in these streams, these results suggest that the hyporheic zones in Antarctic streams have near‐stream zones of rapid stream‐water exchange, where ‘fast’ biogeochemical reactions may influence water chemistry, and extended hyporheic zones, in which slower biogeochemical reaction rates may affect stream‐water chemistry at longer time‐scales. Copyright © 2003 John Wiley & Sons, Ltd.     

Macropores as preferential flow paths in meander bends

Macropores are subsurface connected void spaces caused by processes such as fracture of soils, micro‐erosion, and fauna burrows. They are common near streams (e.g. hyporheic and riparian zones) and may act as preferential flow paths between surface and groundwaters, affecting hydrologic and biogeochemical processes. We tested the hydrologic function of macropores by constructing an artificial macropore within the saturated zone of a meander bend (open macropore, ‘OM’) and later filling its upstream end (partially filled macropore, ‘PFM’). For each treatment, we injected saline tracer at an upgradient monitoring well within the meander and monitored downgradient hydraulics and tracer transport. Pressure transducers in monitoring wells indicated hydraulic gradients within the meander were 32% higher perpendicular to and 6% higher parallel to the macropore for the OM than for the PFM. Additionally, hydraulic conductivities measured via falling head tests were 29 to 550 times higher along the macropore than in nearby sediment. We used electrical conductivity probes in wells and electrical resistivity imaging to track solute transport. Transport velocities through the meander were on average 9 and 21% higher (per temporal moment analysis and observed tracer peak, respectively) for the OM than for the PFM. Furthermore, temporal moments of tracer breakthrough analysis indicated downgradient longitudinal dispersion and breakthrough tracer curve tailing were on average 234% and 182% higher for the OM, respectively. This suggests the OM enabled solute transport at overall shorter timescales than the matrix but also increased tailing. Our results demonstrate the importance of macropores to meander bend hydrology and solute transport. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of oscillating flow on hyporheic zone development

The hyporheic zone is an ecologically important ecotone that describes the extent to which nutrient-rich surface waters penetrate the shallow subsurface adjacent to a flowing surface water body. Although steady-state models satisfactorily explain the incursion of surface water into the subsurface as a function of head gradients developed across streambed riffles, they fail to account for the depth that surface water is observed to penetrate the subsurface or for the extent to which the hyporheic zone develops adjacent to the stream channel. To investigate these issues, transient flow modeling has been conducted at the riffle scale and supported by data for an instrumented site in northern Ontario where stream-stage fluctuations are strictly regulated. Model results show that daily stream-stage fluctuations between 0.6 and 4 m produce oscillating solute flow paths that typically reduce residence times of water and solutes in the hyporheic zone from 60 days or more under steady-state conditions to less than 1 day. Furthermore, similar stream-stage fluctuations increase the depth that solutes pervade the subsurface and banks lateral to the stream from around 1 m under steady-state conditions to as much as 2 and 10 m, respectively. The results demonstrate that the transient flow conditions triggered in the subsurface by variable stream stage can exert a strong influence on hyporheic zone development and have important implications for the hyporheos. The results are especially important for hyporheic communities that may survive gradual changes to their living conditions by migrating to more hospitable aquatic habitats, but are unable to respond to rapid changes provoked by more extreme hydrological events.     

The effect of density-driven flow on the transport of high-concentration solutes in the hyporheic zone

In this study, both laboratory experiments and numerical simulations were conducted to investigate the effect of density-driven flow on the transport of high-concentration pollutants in the hyporheic zone. The results show that the density gradient can change the flow of pore water and the strong density-driven flow can lead to an unstable flow, which increases the effect of preferential flow and thus causes the appearance of solute fingers in the hyporheic zone. Notably, these solute fingers become more obvious with the increase of depth. The appearance of solute fingers depends on the relative strength of the pumping exchange and density gradient, which are represented by the dimensionless number M* and N* respectively. Finger flows appear near the interface when M^* is less than 0.5 N^*. This study may contribute to better understanding the transport and destination of solutes and thus may provide some insights into the assessment on pollution incidents.     

Effect of rapidly changing river stage on uranium flux through the hyporheic zone

Measurement of ground water/surface water interaction within the hyporheic zone is increasingly recognized as an important aspect of subsurface contaminant fate and transport. Understanding the interaction between ground water and surface water is critical in developing a complete conceptual model of contaminant transport through the hyporheic zone. At the Hanford Site near Richland, Washington, ground water contaminated with uranium discharges to the Columbia River through the hyporheic zone. Ground water flux varies according to changes in hydraulic gradient caused by fluctuating river stage, which changes in response to operation of dams on the Columbia River. Piezometers and continuous water quality monitoring probes were installed in the hyporheic zone to provide long-term, high-frequency measurement of hydraulic gradient and estimated uranium concentrations. Subsequently, the flux of water and uranium was calculated for each half-hour time period over a 15-month study period. In addition, measurement of water levels in the near-shore unconfined aquifer enhanced the understanding of the relationship between river stage, aquifer elevation, and uranium flux. Changing river stage resulted in fluctuating hydraulic gradient within the hyporheic zone. Further, influx of river water caused lower uranium concentrations as a result of dilution. The methods employed in this study provide a better understanding of the interaction between surface and ground water in a situation with a dynamically varying vertical hydraulic gradient and illustrate how the combination of relatively standard methods can be used to derive an accurate estimation of water and contaminant flux through the hyporheic zone.     

Colloid transport and distribution in the hyporheic zone

Colloids moving from the stream into the hyporheic zone may have a negative impact on aquatic ecosystems as they are potential contaminants or carriers of contaminants. Moreover, retained colloids in the hyporheic zone could not only reduce the exchange flux between the stream and streambed but also change the conditions of the bed, affecting the habitats for aquatic organisms. Previous studies focused on the exchange flux across the sediment–water interface, but the colloid transport processes and distribution of retained colloids in the streambed have received little attention. We conducted experiments within a laboratory flume to examine these processes in a streambed driven by bedform‐induced hyporheic flow. Retained colloids measured in the bed at the end of the experiments revealed colloid retention mainly in the shallow layer of hyporheic zone (0–5 cm below the interface). The results demonstrated significant effects of particle trapping and settling on the colloid transport and distribution in the streambed. Retention leads to the formation of a colloid‐filled shallow layer in the bed. Particle paths based on model simulations showed that colloid settling in pore water modifies the direction of colloid transport and allows the colloid particles to move more deeply in the bed.     

Hillslope water tracks in the High Arctic: Seasonal flow dynamics with changing water sources in preferential flow paths

Preferential subsurface flow paths known as water tracks are often the principal hydrological pathways of headwater catchments in permafrost areas, exerting an influence on slope physical and biogeochemical processes. In polar deserts, where water resources depend on snow redistribution, water tracks are mostly found in hydrologically active areas downslope from snowdrifts. Here, we measured the flow through seeping water track networks and at the front of a perennial snowdrift, at Ward Hunt Island in the Canadian High Arctic. We also used stable isotope analysis to determine the origin of this water, which ultimately discharges into Ward Hunt Lake. These measurements of water track hydrology indicated a glacio‐nival run‐off regime, with flow production mechanisms that included saturation overland flow (return flow) in a low sloping area, throughflow or pipe‐like flow in most seepage locations, and infiltration excess overland flow at the front of the snowdrift. Each mechanism delivered varying proportions of snowmelt and ground water, and isotopic compositions evolved during the melting season. Unaltered snowmelt water contributed to >90% of total flow from water track networks early in the season, and these values fell to <5% towards the end of the melting season. In contrast, infiltration excess overland flow from snowdrift consisted of a steady percentage of snowmelt water in July (mean of 69%) and August (71%). The water seeping at locations where no snow was left in August 2015 was isotopically enriched, indicating a contribution of the upper, ice‐rich layer of permafrost to late summer discharge during warmer years. Air temperature was the main driver of snowmelt, but the effect of slope aspect on solar radiation best explained the diurnal discharge variation at all sites. The water tracks in this polar desert are part of a patterned ground network, which increases connectivity between the principal water sources (snowdrifts) and the bottom of the slope. This would reduce soil–water interactions and solute release, thereby favouring the low nutrient status of the lake.     

Spatial variability in the hyporheic zone refugium of temporary streams

A key ecological role hypothesized for the hyporheic zone is as a refugium that promotes survival of benthic invertebrates during adverse conditions in the surface stream. Many studies have investigated use of the hyporheic refugium during hydrological extremes (spates and streambed drying), and recent research has linked an increase in the abundance of benthic invertebrates within hyporheic sediments to increasing biotic interactions during flow recession in a temporary stream. This study examined spatial variability in the refugial capacity of the hyporheic zone in two groundwater-dominated streams in which flow permanence varied over small areas. Two non-insect taxa, Gammarus pulex and Polycelis spp. were common to both streams and were investigated in detail. Hydrological conditions in both streams comprised a four-month period of flow recession and low flows, accompanied by reductions in water depth and wetted width. Consequent declines in submerged benthic habitat availability were associated with increases in population densities of mobile benthic taxa, in particular G. pulex. The reduction in the spatial extent of the hyporheic zone was minimal, and this habitat was therefore a potential refugium from increasing biotic interactions in the benthic sediments. Concurrent increases in the hyporheic abundance and hyporheic proportion of a taxon’s total (benthic + hyporheic) population were considered as evidence of active refugium use. Such evidence was species-specific and site-specific, with refugium use being observed only for G. pulex and at sites dominated by downwelling water. A conceptual model of spatial variability in the refugial capacity of the hyporheic zone during habitat contraction is presented, which highlights the potential importance of the direction of hydrologic exchange.     

Heterogeneous oxygenation of wetland soils with increasing inundation: Redox potential,water depth,and preferential flow paths

Saturated wetland soils are typically depleted in oxygen due to rapid consumption by biological processes. Fine-grained sediments limit the vertical flux of oxygenated surface water. In forested wetlands potential exists for advective transport along preferential flow pathways created by decomposing roots or tree remains, resulting in isolated zones of oxygenated pore water. This potential is enhanced in wetlands that are perched above the regional water table, where increases in surface-water depth create significant downward hydraulic gradients. Along the Mississippi River floodplain, many wetlands are now perched as a result of groundwater mining. This study was conducted at Sky Lake, MS, through a growing season, with high temporal-resolution monitoring of redox potential as a proxy for the delivery of oxygenated surface water at 30 and 60 cm depths at six locations. Results showed reducing conditions at all locations under moderate inundation, but redox potential rose significantly at a subset of locations when surface-water exceeded one meter in depth. Rising redox potential is argued to represent localised advective transport of oxygen in excess of the rate of biological consumption. Wetland plants growing in perennially saturated soils may benefit from periodic increases in water depth and subsequent increase in the flux of oxygen along preferential flow paths.     

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.     

Hydrologic connectivity at the hillslope scale through intra-snowpack flow paths during snowmelt

The seasonal snowmelt period is a critical component of the hydrologic cycle for many mountainous areas. Changes in the timing and rate of snowmelt as a result of physical hydrologic flow paths, such as longitudinal intra-snowpack flow paths, can have strong implications on the partitioning of meltwater amongst streamflow, groundwater recharge, and soil moisture storage. However, intra-snowpack flow paths are highly spatially and temporally variable and thus difficult to observe. This study utilizes new methods to non-destructively observe spatio-temporal changes in the liquid water content of snow in combination with plot experiments to address the research question: What is the scale of influence that intra-snowpack flow paths have on the downslope movement of liquid water during snowmelt across an elevational gradient? This research took place in northern Colorado with study plots spanning from the rain-snow transition zone up to the high alpine. Results indicate an increasing scale of influence from intra-snowpack flow paths with elevation, showing higher hillslope connectivity producing larger intra-snowpack contributing areas for meltwater accumulation, quantified as the upslope contributing area required to produce observed changes in liquid water content from melt rate estimates. The total effective intra-snowpack contributing area of accumulating liquid water was found to be 17, 6, and 0 m² for the above tree line, near tree line, and below tree line plots, respectively. Dye tracer experiments show capillary and permeability barriers result in increased number and thickness of intra-snowpack flow paths at higher elevations. We additionally utilized aerial photogrammetry in combination with ground penetrating radar surveys to investigate the role of this hydrologic process at the small watershed scale. Results here indicate that intra-snowpack flow paths have influence beyond the plot scale, impacting the storage and transmission of liquid water within the snowpack at the small watershed scale.     

Distribution and species composition of macroinvertebrates in the hyporheic zone of bed sediment

The hyporheic zone is a layer of substrate on a river bed where benthic animals normally live,grow,feed,reproduce,and exist for any portion of their life cycle.The hyporheic zone was studied by samplin...     

Methane in the hyporheic zone of a small lowland stream (Sitka, Czech Republic)

Methane concentrations and selected chemical parameters in interstitial water were examined along subsurface flowpath in two subsystems (hyporheic and parafluvial sediments) in the Sitka stream, Czech Republic. Interstitial methane concentrations exhibited a distinct spatial pattern. In the hyporheic downwelling zone where the sediments are relatively well oxygenated due to high hydrologic exchange with the surface water, low interstitial methane concentrations, averaging 9.3 μg CH₄/l, were found. In contrast, upwelling sediments and parafluvial sediments (active channel sediments lateral to the wetted channel) had significantly higher methane concentrations (p < 0.05, and p < 0.01, respectively), averaging 43.2 μg CH₄/l and 160.5 μg CH₄/l, respectively. Dissolved oxygen was the highest where surface water entered hyporheic/parafluvial sediments and decreased with water residence time in the sediments (p < 0.01). Nitrate concentrations decreased along the flowpath and were significantly lower at downstream end of the riffle (p < 0.001). Sulfate concentrations also show a slight decline with the water residence time, but differences were not significant. Effect of both nitrate and sulfate on methanogenesis is also discussed. The interstitial methane concentration significantly increased with surface water temperature (p < 0.001) and was negatively correlated with redox potential (p < 0.01) and dissolved oxygen (p < 0.05).     

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.     

Impacts of river bed gas on the hydraulic and thermal dynamics of the hyporheic zone

Despite the presence of gas in river beds being a well known phenomenon, its potential feedbacks on the hydraulic and thermal dynamics of the hyporheic zone has not been widely studied. This paper explores hypotheses that the presence of accumulated gas impacts the hydraulic and thermal dynamics of a river bed due to changes in specific storage, hydraulic conductivity, effective porosity, and thermal diffusivity. The hypotheses are tested using data analysis and modelling for a study site on the urban River Tame, Birmingham, UK. Gas, predominantly attributed to microbial denitrification, was observed in the river bed up to around 14% by volume, and to at least 0.8 m depth below river bed. Numerical modelling indicates that, by altering the relative hydraulic conductivity distribution, the gas in the river bed leads to an increase of groundwater discharge from the river banks (relative to river bed) by a factor of approximately 2 during river low flow periods. The increased compressible storage of the gas phase in the river bed leads to an increase in the simulated volume of river water invading the river bed within the centre of the channel during storm events. The exchange volume can be more than 30% greater in comparison to that for water saturated conditions. Furthermore, the presence of gas also reduces the water-filled porosity, and so the possible depth of such invading flows may also increase markedly, by more than a factor of 2 in the observed case. Observed diurnal temperature variations within the gaseous river bed at 0.1 and 0.5 m depth are, respectively, around 1.5 and 6 times larger than those predicted for saturated sediments. Annual temperature fluctuations are seen to be enhanced by around 4 to 20% compared to literature values for saturated sediments. The presence of gas may thus alter the bulk thermal properties to such a degree that the use of heat tracer techniques becomes subject to a much greater degree of uncertainty. Although the likely magnitude of thermal and hydraulic changes due to the presence of gas for this site have been demonstrated, further research is needed into the origins of the gas and its spatial and temporal variability to enable quantification of the significance of these changes for chemical attenuation and hyporheic zone biology.     

Salinity‐induced increase of the hydraulic conductivity in the hyporheic zone of coastal wetlands

In coastal zones globally, salinization is rapidly taking place due to the combined effects of sea level rise, land subsidence, altered hydrology, and climate change. Although increased salinity levels are known to have a great impact on both biogeochemical and hydrological processes in aquatic sediments, only few studies have included both types of processes and their potential interactions. In the present paper, we used a controlled 3‐year experimental mesocosm approach to test salinity induced interactions and discuss mechanisms explaining the observed hydrological changes. Surface water salinity was experimentally increased from 14 to 140 mmol Cl per L (0.9 and 9 PSU) by adding sea salt which increased pore water salinity but also increased sulfate reduction rates, leading to higher sulfide, and lower methane concentrations. By analyzing slug test data with different slug test analysis methods, we were able to show that hydraulic conductivity of the hyporheic zone increased 2.8 times by salinization. Based on our hydrological and biogeochemical measurements, we conclude that the combination of pore dilation and decreased methane production rates were major controls on the observed increase in hydraulic conductivity. The slug test analysis method comparison allowed to conclude that the adjusted Bouwer and Rice method results in the most reliable estimate of the hydraulic conductivity for hyporheic zones. Our work shows that both physical and biogeochemical processes are vital to explain and predict hydrological changes related to the salinization of hyporheic zones in coastal wetlands and provides a robust methodological approach for doing so.     

Residence time control on hot moments of net nitrate production and uptake in the hyporheic zone

The retention capacity for biologically available nitrogen within streams can be influenced by dynamic hyporheic zone exchange, a process that may act as either a net source or net sink of dissolved nitrogen. Over 5 weeks, nine vertical profiles of streambed chemistry (NO₃^? and NH₄⁺) were collected above two beaver dams along with continuous high‐resolution vertical hyporheic flux data. The results indicate a non‐linear relation of net NO₃^? production followed by net uptake in the hyporheic zone as a function of residence time. This Lagrangian‐based relation is consistent through time and across varied morphology (bars, pools, glides) above the dams, even though biogeochemical and environmental factors varied. The empirical continuum between net NO₃^? production and uptake and residence time is useful for identifying two crucial residence time thresholds: the transition to anaerobic respiration, which corresponds to the time of peak net nitrate production, and the net sink threshold, which is defined by a net uptake in NO₃^? relative to streamwater. Short‐term hyporheic residence time variability at specific locations creates hot moments of net production and uptake, enhancing NO₃^? production as residence times approach the anaerobic threshold, and changing zones of net NO₃^? production to uptake as residence times increase past the net sink threshold. The anaerobic and net sink thresholds for beaver‐influenced streambed morphology occur at much shorter residence times (1.3 h and 2.3 h, respectively) compared to other documented hyporheic systems, and the net sink threshold compares favorably to the lower boundary of the anaerobic threshold determined for this system with the new oxygen Damkohler number. The consistency of the residence time threshold values of NO₃^? cycling in this study, despite environmental variability and disparate morphology, indicates that NO₃^? hot moment dynamics are primarily driven by changes in physical hydrology and associated residence times. Copyright © 2013 John Wiley & Sons, Ltd.     

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.