High-Resolution Integrated Transport Model for Studying Surface Water–Groundwater Interaction

Transport processes that lead to exchange of mass between surface water and groundwater play a significant role for the ecological functioning of aquatic systems, for hydrological processes and for biogeochemical transformations. In this study, we present a novel integral modeling approach for flow and transport at the sediment–water interface. The model allows us to simultaneously simulate turbulent surface and subsurface flow and transport with the same conceptual approach. For this purpose, a conservative transport equation was implemented to an existing approach that uses an extended version of the Navier–Stokes equations. Based on previous flume studies which investigated the spreading of a dye tracer under neutral, losing and gaining flow conditions the new solver is validated. Tracer distributions of the experiments are in close agreement with the simulations. The simulated flow paths are significantly affected by in- and outflowing groundwater flow. The highest velocities within the sediment are found for losing condition, which leads to shorter residence times compared to neutral and gaining conditions. The largest extent of the hyporheic exchange flow is observed under neutral condition. The new solver can be used for further examinations of cases that are not suitable for the conventional coupled models, for example, if Reynolds numbers are larger than 10. Moreover, results gained with the integral solver provide high-resolution information on pressure and velocity distributions at the rippled streambed, which can be used to improve flow predictions. This includes the extent of hyporheic exchange under varying ambient groundwater flow conditions.     

Hyporheic exchange flows induced by constructed riffles and steps in lowland streams in southern Ontario,Canada

Stream–subsurface water interaction induced by natural riffles and constructed riffles/steps was examined in lowland streams in southern Ontario, Canada. The penetration of stream water into the subsurface was analysed using hydrometric data, and the zone of > 10% stream water was calculated from a chemical mixing equation using tracer injection of bromide and background chloride concentrations. The constructed riffles studied induced more extensive hyporheic exchange than the natural riffles because of their steeper longitudinal hydraulic head gradients and coarser streambed sediments. The depth of > 10% stream water zone in a small and a large constructed riffle extended to > 0·2 m and > 1·4 m depths respectively. Flux and residence time distribution of hyporheic exchange were simulated in constructed riffles using MODFLOW, a finite‐difference groundwater flow model. Hyporheic flux and residence time distribution varied along the riffles, and the exchange occurring upstream from the riffle crest was small in flux and had a long residence time. In contrast, hyporheic exchange occurring downstream from the riffle crest had a relatively short residence time and accounted for 83% and 70% of total hyporheic exchange flow in a small and large riffle respectively. Although stream restoration projects have not considered the hyporheic zone, our data indicate that constructed riffles and steps can promote vertical hydrologic exchange and increase the groundwater–surface water linkage in degraded lowland streams. Copyright © 2006 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.     

A Spatially Periodic Solute Boundary for MT3DMS and PHT3D

The assumption of spatial repetition is commonly made when producing bedform scale models of the hyporheic zone. Two popular solute transport codes, MT3DMS and PHT3D, do not currently provide the necessary boundary condition required to simulate spatial periodicity in hyporheic zone transport problems. In this study, we develop a spatially periodic boundary (SPB) for solutes that is compatible with a SPB that was previously developed for MODFLOW to simulate the flow component of spatially periodic problems. The approach is ideal for simulating groundwater flow and transport patterns under repeating surface features, such as ripples or dunes on the bottom of a lake or stream. The appropriate block‐centered finite‐difference approach to implement the boundary is presented and the necessary source code modifications are discussed. The performance of the solute SPB, operating in conjunction with the groundwater flow SPB, is explored through comparison of a multi‐bedform hyporheic‐zone model with a single bedform variant. The new boundary conditions perform well in situations where both dispersive effects and lateral seepage flux in the underflow regime beneath the hyporheic zone are minimal.     

Focused Groundwater Controlled Feedbacks into the Hyporheic Zone During Baseflow Recession

Groundwater surface water interaction in the hyporheic zone remains an important challenge for water resources management and ecosystem restoration. In heterogeneous stratified glacial sediments, reach‐scale environments contain an uneven distribution of focused groundwater flow occurring simultaneously with diffusely discharging groundwater. This results in a variation of stream‐aquifer interactions, where focused flow systems are able to temporally dominate exchange processes. The research presented here investigates the direct and indirect influences focused groundwater discharge exerts on the hyporheic zone during baseflow recession. Field results demonstrate that as diffuse sources of groundwater deplete during baseflow recession, focused groundwater discharge remains constant. During baseflow recession the hyporheic zone is unable to expand, while the high nitrate concentration from focused discharge changes the chemistry of the stream. The final result is a higher concentration of nitrate in the hyporheic zone as this altered surface water infiltrates into the subsurface. This indirect coupling of focused groundwater discharge and the hyporheic zone is unaccounted for in hyporheic studies at this time. Results indicate important implications for the potential reduction of agricultural degradation of water quality.     

Residence time of bedform-driven hyporheic exchange

Biogeochemical and ecological transformations in hyporheic zones are dependent on the timing of hyporheic exchange. We show through linked modeling of open channel turbulent flow, groundwater flow, and solute transport that the residence time distributions of solutes advected by hyporheic flow induced by current–bedform interaction follow power-laws. This tailing behavior of solutes exiting the sediments is explained by the presence of multiple path lengths coupled with very large variability in Darcy flow velocity, both occurring without heterogeneity in sediment permeability. Hyporheic exchange through bedforms will result in short-time fractal scaling of stream water chemistry.     

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.     

Stream water infiltration, bank storage, and storage zone changes due to stream-stage fluctuations

During a flood period, stream-stage increases induce infiltration of stream water into an aquifer; subsequent declines in stream stage cause a reverse motion of the infiltrated water. This paper presents the results of the water exchange rate between a stream and aquifer, the storage volume of the infiltrated stream water in the surrounding aquifer (bank storage), and the storage zone. The storage zone is the part of aquifer where groundwater is replaced by stream water during the flood. MODFLOW was used to simulate stream–aquifer interactions and to quantify rates of stream infiltration and return flow. MODPATH was used to trace the pathlines of the infiltrated stream water and to determine the size of the storage zone. Simulations were focused on the analyses of the effects of the stream-stage fluctuation, aquifer properties, the hydraulic conductivity of streambed sediments, regional hydraulic gradients, and recharge and evapotranspiration (ET) rates on stream–aquifer interactions. Generally, for a given stream–aquifer system, larger flow rates result from larger stream-stage fluctuations; larger storage volumes and storage zones are produced by larger and longer-lasting fluctuations. For a given stream-stage hydrograph, a lower-permeable streambed, an aquitard, or an anisotropic aquifer of low vertical hydraulic conductivity can significantly reduce the rate of infiltration and limit the size of the storage zone. The bank storage solely caused by the stage fluctuation differs slightly between gaining and losing streams. Short-term rainfall recharge and ET loss in the shallow groundwater slightly influence on the flow rate, but their effects on bank storage in a larger area for a longer period can be considerable.     

The influence of local hydrogeologic forcings on near‐stream event water recharge and retention (Upper San Pedro River,Arizona)

The rise in stream stage during high flow events (floods) can induce losing stream conditions, even along stream reaches that are gaining during baseflow conditions. The aquifer response to flood events can affect the geochemical composition of both near‐stream groundwater and post‐event streamflow, but the amount and persistence of recharged floodwater may differ as a function of local hydrogeologic forcings. As a result, this study focuses on how vertical flood recharge varies under different hydrogeologic forcings and the significance that recharge processes can have on groundwater and streamflow composition after floods. River and shallow groundwater samples were collected along three reaches of the Upper San Pedro River (Arizona, USA) before, during and after the 2009 and 2010 summer monsoon seasons. Tracer data from these samples indicate that subsurface floodwater propagation and residence times are strongly controlled by the direction and magnitude of the dominant stream–aquifer gradient. A reach that is typically strongly gaining shows minimal floodwater retention shortly after large events, whereas the moderately gaining and losing reaches can retain recharged floodwater from smaller events for longer periods. The moderately gaining reach likely returned flood recharge to the river as flow declined. These results indicate that reach‐scale differences in hydrogeologic forcing can control (i) the amount of local flood recharge during events and (ii) the duration of its subsurface retention and possible return to the stream during low‐flow periods. Our observations also suggest that the presence of floodwater in year‐round baseflow is not due to long‐term storage beneath the streambed along predominantly gaining reaches, so three alternative mechanisms are suggested: (i) repeated flooding that drives lateral redistribution of previously recharged floodwater, (ii) vertical recharge on the floodplain during overbank flow events and (iii) temporal variability in the stream–aquifer gradient due to seasonally varying water demands of riparian vegetation. Copyright © 2011 John Wiley & Sons, Ltd.     

From streambed temperature measurements to spatial‐temporal flux quantification: using the LPML method to study groundwater–surface water interaction

Knowledge on groundwater–surface water interaction and especially on exchange fluxes between streams and aquifers is an important prerequisite for the study of transport and fate of contaminants and nutrients in the hyporheic zone. One possibility to quantify groundwater–surface water exchange fluxes is by using heat as an environmlental tracer. Modern field equipment including multilevel temperature sticks and the novel open‐source analysis tool LPML make this technique ever more attractive. The recently developed LPML method solves the one‐dimensional fluid flow and heat transport equation by combining a local polynomial method with a maximum likelihood estimator. In this study, we apply the LPML method on field data to quantify the spatial and temporal variability of vertical fluxes and their uncertainties from temperature–time series measured in a Belgian lowland stream. Over several months, temperature data were collected with multilevel temperature sticks at the streambed top and at six depths for a small stream section. Long‐term estimates show a range from gaining fluxes of ?291 mm day^?1 to loosing fluxes of 12 mm day^?1; average seasonal fluxes ranged from ?138 mm day^?1 in winter to ?16 mm day^?1 in summer. With our analyses, we could determine a high spatial and temporal variability of vertical exchange fluxes for the investigated stream section. Such spatial and temporal variability should be taken into account in biogeochemical cycling of carbon, nutrients and metals and in fate analysis of contaminant plumes. In general, the stream section was gaining during most of the observation period. Two short‐term high stream stage events, seemingly caused by blockage of the stream outlet, led to a change in flow direction from gaining to losing conditions. We also found more discharge occurring at the outer stream bank than at the inner one indicating a local flow‐through system. With the conducted analyses, we were able to advance our understanding of the regional groundwater flow system. Copyright © 2015 John Wiley & Sons, Ltd.     

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

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

A temporal stable isotopic (δ18O, δD,d‐excess) comparison in glacier meltwater streams,Taylor Valley,Antarctica

In this paper, we describe the importance of hyporheic dynamics within Andersen Creek and Von Guerard Stream, Taylor Valley, Antarctica, from the 2010–2011 melt season using natural tracers. Water collection started at flow onset and continued, with weekly hyporheic‐zone sampling. The water δ¹⁸O and δD values were isotopically lighter in the beginning and heavier later in the season. D‐excess measurements were used as an indicator of mixing because an evaporative signature was evident and distinguishable between 2 primary end‐members (glacier meltwater and hyporheic zone). Hyporheic‐zone influence on the channel water was variable with a strong control on streamwater chemistry, except at highest discharges. This work supports previous research indicating that Von Guerard Stream has a large, widespread hyporheic zone that varies in size with time and discharge. Andersen Creek, with a smaller hyporheic zone, displayed hyporheic‐zone solute interaction through the influence from subsurface hypersaline flow. Overall, the evolution of Taylor Valley hyporheic‐zone hydrology is described seasonally. In mid‐December, the hyporheic zone is a dynamic system exchanging with the glacier meltwater in the channel, and with diminishing flow in January, the hyporheic zone drains back into the channel flow also impacting stream chemistry. This work adds new information on the role of hyporheic zone–stream interaction in these glacier meltwater streams.     

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.     

Combination of CFCs and stable isotopes to characterize the mechanism of groundwater–surface water interactions in a headwater basin of the North China Plain

Mountainous areas are characterized by steep slopes and rocky landforms, with hydrological conditions varying rapidly from upstream to downstream, creating variable interactions between groundwater and surface water. In this study, mechanisms of groundwater–surface water interactions within a headwater catchment of the North China Plain were assessed along the stream length and during different seasons, using hydrochemical and stable isotope data, and groundwater residence times estimated using chlorofluorocarbons. These tracers indicate that the river is gaining, due to groundwater discharge in the headwater catchment both in the dry and rainy seasons. Residence time estimation of groundwater using chlorofluorocarbons data reveals that groundwater flow in the shallow sedimentary aquifer is dominated by the binary mixing of water approximating a piston flow model along 2 flow paths: old water, carried by a regional flow system along the direction of river flow, along with young water, which enters the river through local flow systems from hilly areas adjacent to the river valley (particularly during the rainy season). The larger mixing ratio of young water from lateral groundwater recharge and return flow of irrigation during the rainy season result in higher ion concentrations in groundwater than in the dry season. The binary mixing model showed that the ratio of young water versus total groundwater ranged from 0.88 to 0.22 and 1.0 to 0.74 in the upper and lower reaches, respectively. In the middle reach, meandering stream morphology allows some loss of river water back into the aquifer, leading to increasing estimates of the ratio of young water (from 0.22 to 1). This is also explained by declining groundwater levels near the river, due to groundwater extraction for agricultural irrigation. The switch from a greater predominance of regional flow in the dry season, to more localized groundwater flow paths in the wet season is an important groundwater–surface water interactions mechanism, with important catchment management implications.     

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.     

Does streambed heterogeneity matter for hyporheic residence time distribution in sand-bedded streams?

Stream water residence times within streambed sediments are key values to quantify hyporheic processes including sediment thermal regime, solute transient storage, dilution rates and biogeochemical transformations, such as those controlling degassing nitrous oxide. Heterogeneity of the streambed sediment hydraulic properties has been shown to be potentially an important factor to characterize hyporheic processes. Here, we quantify the importance of streambed heterogeneity on residence times of dune-like bedform induced hyporheic fluxes at the bedform and reach scales. We show that heterogeneity has a net effect of compression of the hyporheic zone (HZ) toward the streambed, changing HZ volume from the homogenous case and thus inducing remarkable differences in the flow field with respect to the homogeneous case. We unravel the physical conditions for which the commonly used homogenous field assumption is applicable for quantifying hyporheic processes thus explaining why predictive measures based on a characteristic residence time, like the Damköhler number, are robust in heterogeneous sand bedded streams.     

Bedform-induced hyporheic exchange with unsteady flows

The interaction between surface and subsurface water has a crucial influence on the biochemistry of stream environments. Even though the river discharge and the flow conditions can seldom be considered to be steady, the influence of this unsteadiness on the hyporheic exchange has often been neglected. In this work, a model for the study of hyporheic exchange during unsteady conditions has been developed. The model provides a sound analytical framework for the analysis of the effects of a varying stream discharge on the exchange between a stream and the hyporheic zone. The effects of the unsteadiness on the water exchange flux, the residence time of the solutes in the bed, and the stored mass are quantified. A synthetic example shows the substantial influence of a flood on the hyporheic exchange, and that the application of a steady model can lead to an underestimation of the exchanged mass, even after the flood has ended.     

From Grain to Floodplain: Evaluating heterogeneity of floodplain hydrostatigraphy using sedimentology,geophysics, and remote sensing

Stratigraphy is a fundamental component of floodplain heterogeneity and hydraulic conductivity and connectivity of alluvial aquifers, which affect hydrologic processes such as groundwater flow and hyporheic exchange. Watershed-scale hydrological models commonly simplify the sedimentology and stratigraphy of floodplains, neglecting natural floodplain heterogeneity and anisotropy. This study, conducted in the upper reach of the East River in the East River Basin, Colorado, USA, combines point-, meander-, and floodplain-scale data to determine key features of alluvial aquifers important for estimating hydrologic processes. We compare stratigraphy of two meanders with disparate geometries to explore floodplain heterogeneity and connectivity controls on flow and transport. Meander shape, orientation, and internal stratigraphy affected residence time estimates of laterally exchanged hyporheic water. Although the two meanders share a sediment source, vegetation, and climate, their divergent river migration histories resulted in contrasting meander hydrofacies. In turn, the extent and orientation of these elements controlled the effective hydraulic conductivity and, ultimately, estimates of groundwater transport and hyporheic residence times. Additionally, the meanders’ orientation relative to the valley gradient impacted the hydraulic gradient across the meanders—a key control of groundwater velocity. Lastly, we combine our field data with remotely sensed data and introduce a potential approach to estimate key hydrostratigraphic packages across floodplains. Prospective applications include contaminant transport studies, hyporheic models, and watershed models. © 2019 John Wiley & Sons, Ltd.     

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.     

Melt water driven stream and groundwater stage fluctuations on a glacier forefield (Dammagletscher,Switzerland)

In many mountain regions, large land areas with heterogeneous soils have become ice‐free with the ongoing glacier retreat. On these recently formed proglacial fields, the melt of the remaining glaciers typically drives pronounced diurnal stream level fluctuations that propagate into the riparian zone. This behaviour was measured on the Damma glacier forefield in central Switzerland with stage recorders in the stream and groundwater monitoring wells along four transects. In spite of the large groundwater stage variations, radon measurements in the near‐stream riparian zone indicate that there is little mixing between stream water and groundwater on daily time scales. At all four transects, including both losing and gaining reaches, the groundwater level fluctuations lagged the stream stage variations and were often damped with distance from the stream. Similar behaviours have been modelled using the diffusion equation in coastal regions influenced by tidal sea level variations. We thus tested the ability of such a model to predict groundwater level fluctuations in proglacial fields. The model reproduced several key features of the observed fluctuations at three of four locations, although discrepancies also arise due to non representative input data and model simplifications. Nevertheless, calibration of the model for the individual transects yielded realistic estimates of hydraulic diffusivities between the stream and groundwater monitoring wells. We conclude that studying diurnal groundwater fluctuations can provide important information about the subsurface hydrology of alpine watersheds dominated by glacier melt. Copyright © 2012 John Wiley & Sons, Ltd.