Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 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.     

Scales and magnitude of hyporheic,river–aquifer and bank storage exchange fluxes

Many studies have investigated the exchange processes that occur between rivers and groundwater systems and have successfully quantified the water fluxes involved. Specifically, these exchange processes include hyporheic exchange, river–aquifer exchange (groundwater discharge and river loss) and bank storage exchange. Remarkably, there are relatively few examples of field studies where more than one exchange process is quantified, and as a consequence, the relationships between them are not well understood. To compare the relative magnitudes of these common exchange processes, we have collected data from 54 studies that have quantified one or more of these exchange flux types. Each flux value is plotted against river discharge at the time of measurement to allow the different exchange flux types to be compared. We show that there are positive relationships between the magnitude of each exchange flux type and increasing river discharge across the different studies. For every one order of magnitude increase in river discharge, the hyporheic, river–aquifer and bank storage exchange fluxes increase by factors of 2.7, 2.9 and 2.5, respectively. On average, hyporheic exchange fluxes are almost an order of magnitude greater than river–aquifer exchange fluxes, which are, in turn, approximately four times greater than bank storage exchange fluxes for the same river discharge. Unless measurement approaches that can distinguish between different types of exchange flux are used, there is potential for hyporheic exchange fluxes to be misinterpreted as river–aquifer exchange fluxes, with possible implications for water resource management decisions. Copyright © 2015 John Wiley & Sons, Ltd.     

Surface water–groundwater interaction at restored streams and associated reference reaches

Stream restoration goals include improving habitat and water quality through reconstruction of morphological features found at analogous, pristine stream reaches. Enhancing hyporheic exchange may facilitate achieving these goals. Although hyporheic exchange at restoration sites has been explored in a few previous studies, comparative studies of restored versus reference or control streams are largely absent. We hypothesized that restoration cross‐vanes enhance hyporheic exchange, resulting in biogeochemical alteration of stream water chemistry in the streambed. Two streams restored using cross‐vanes to control erosion and improve habitat were compared with their associated reference reaches, which provided the basis for the restoration design. Thirteen temperature profile rods with vertically stacked sensors were installed at each site for 2 weeks. Heat tracing was used to quantify vertical flux in the streambed from the diurnal temperature fluctuations in the subsurface. Stream water and bed pore waters from mini‐piezometers were analysed for ion and nutrient chemistry. In general, mean vertical flux rates through the streambed were small throughout reference sites (?0.3 to 0.3 m/day) and at most locations at restored sites. Immediately adjacent to cross‐vanes, vertical flux rates were larger (up to 3.5 m/day). Geochemistry of pore waters shows distinct differences in the sources for the reference and restored sites. Strong downwelling zones adjacent to cross‐vanes showed high dissolved oxygen (10.75 mg/l) and geochemistry in the streambed similar to surface water. Reference sites had lower dissolved oxygen in the streambed (0.66–5.14 mg/l), and geochemical patterns suggest a mixture of discharging groundwater and surface water in the hyporheic zone. Restored sites also clearly show sulfate and nitrate reduction occurring in the streambed, which is not observed at the reference sites. The stream restoration sites studied here enhance rapid hyporheic exchange, but upwelling of groundwater has a stronger influence on streambed geochemistry at reference sites. Copyright © 2012 John Wiley & Sons, Ltd.     

Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns

This study investigates spatial patterns and temporal dynamics of aquifer–river exchange flow at a reach of the River Leith, UK. Observations of sub‐channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up‐welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0–15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the ‘window of detection’, i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0·91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non‐monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater‐surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer‐river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright © 2011 John Wiley & Sons, Ltd.     

Vertical surface water–groundwater exchange processes within a headwater floodplain induced by experimental floods

Restoring hydrologic connectivity between channels and floodplains is common practice in stream and river restoration. Floodplain hydrology and hydrogeology impact stream hydraulics, ecology, biogeochemical processing, and pollutant removal, yet rigorous field evaluations of surface water–groundwater exchange within floodplains during overbank floods are rare. We conducted five sets of experimental floods to mimic floodplain reconnection by pumping stream water onto an existing floodplain swale. Floods were conducted throughout the year to capture seasonal variation and each involved two replicate floods on successive days to test the effect of varying antecedent moisture. Water levels and specific conductance were measured in surface water, soil, and groundwater within the floodplain, along with surface flow into and out of the floodplain. Vegetation density varied seasonally and controlled the volume of surface water storage on the floodplain. By contrast, antecedent moisture conditions controlled storage of water in floodplain soils, with drier antecedent moisture conditions leading to increased subsurface storage and slower flood wave propagation across the floodplain surface. The site experienced spatial heterogeneity in vertical connectivity between surface water and groundwater across the floodplain surface, where propagation of hydrostatic pressure, preferential flow, and bulk Darcy flow were all mechanisms that may have occurred during the five floods. Vertical connectivity also increased with time, suggesting higher frequency of floodplain inundation may increase surface water–groundwater exchange across the floodplain surface. Understanding the variability of floodplain impacts on water quality noted in the literature likely requires better accounting for seasonal variations in floodplain vegetation and antecedent moisture as well as heterogeneous exchange flow mechanisms. Copyright © 2016 John Wiley & Sons, Ltd.     

Heterogeneity in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream

Spatial and temporal variability in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream was investigated. Four locations in a 150‐m reach of the stream were studied using hydrometric and hydrochemical tracing techniques. A high degree of hydrological connectivity between the riparian hillslope and the stream channel was indicated at two locations, where hydrochemical changes and hydraulic gradients indicated that the hyporheic zone was dominated by upwelling ground water. The chemistry of ground water reflected relatively long residence times and reducing conditions with high levels of alkalinity and conductivity, low dissolved oxygen (DO) and nitrate. At the other locations, connectivity was less evident and, at most times, the hyporheic zone was dominated by downwelling stream water characterized by high DO, low alkalinity and conductivity. Substantial variability in hyporheic chemistry was evident at fine (<10 m) spatial scales and changed rapidly over the course of hydrological events. The nature of the hydrochemical response varied among locations depending on the strength of local ground water influence. It is suggested that greater emphasis on spatial and temporal heterogeneity in ground water–surface water interactions in the hyporheic zone is necessary for a consideration of hydrochemical effects on many aspects of stream ecology. For example, the survival of salmonid eggs in hyporheic gravels varied considerably among the locations studied and was shown to be associated with variation in interstitial chemistry. River restoration schemes and watershed management strategies based only on the surface expression of catchment characteristics risk excluding consideration of potentially critical subsurface processes. Copyright © 2002 John Wiley & Sons, Ltd.     

Field assessment of surface water–groundwater connectivity in a semi‐arid river basin (Murray–Darling,Australia)

In semi‐arid and arid river basins, understanding the connectivity between rivers and alluvial aquifers is one of the key challenges for the management of groundwater resources. The type of connection present (gaining, losing‐connected, transitional and losing‐disconnected) was assessed at 12 sites along six Murray–Darling Basin river reaches. The assessments were made by measuring the hydraulic head in the riparian zone near the rivers to evaluate if the water tables intersected the riverbeds and by measuring fluid pressure (ψ) in the riverbeds. The rationale for the latter was that ψ will always be greater than or equal to zero under connected conditions (either losing or gaining) and always lesser than or equal to zero under losing‐disconnected conditions. A mixture of losing‐disconnected, losing‐connected and gaining conditions was found among the 12 sites. The losing‐disconnected sites all had a riverbed with a lower hydraulic conductivity than the underlying aquifer, usually in the form of a silty clay or clay unit 0.5–2 m in thickness. The riparian water tables were 6 to 25 m below riverbed level at the losing‐disconnected sites but never lower than 1 m below riverbed level at the losing‐connected ones. The contrast in water table depth between connected and disconnected sites was attributed to the conditions at the time of the study, when a severe regional drought had generated a widespread decline in regional water tables. This decline was apparently compensated near losing‐connected rivers by increased infiltration rates, while the decline could not be compensated at the losing‐disconnected rivers because the infiltration rates were already maximal there. Copyright © 2012 John Wiley & Sons, Ltd.     

Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

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

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

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.     

Using new mass balance methods to estimate gross surface water and groundwater exchange with naturally occurring tracer 222Rn in data poor regions: a case study in northwest China

Given that the concentration of ²²²Rn in groundwater is much higher than that in surface water and that its radioactive half‐life (3.83 d) is short, ²²²Rn is an effective tracer of groundwater–surface water interactions. In this study, a new mass balance method is presented, which can be used to estimate specific groundwater–surface water interactions within a river reach. Three possible situations of interaction between groundwater and surface water are considered, and equations based on the mass conservation of ²²²Rn are formulated for judging specific groundwater–surface water interaction processes and for calculating water flux. A case study was conducted for the Nalenggele River, Northwest China, to demonstrate the usefulness of this method. Samples of river water and groundwater containing ²²²Rn were collected from the study area to estimate the interactions between groundwater and surface water. The amount of water exchanged during these interactions was estimated and the results show that transformations between groundwater and surface water are frequent along the stream. The ²²²Rn mass balance method is highly sensitive for studying such interactions, even in areas for which conventional hydrologic data are sparse. Copyright © 2014 John Wiley & Sons, Ltd.     

A modelling study of hyporheic exchange pattern and the sequence,size, and spacing of stream bedforms in mountain stream networks,Oregon, USA

Studies of hyporheic exchange flows have identified physical features of channels that control exchange flow at the channel unit scale, namely slope breaks in the longitudinal profile of streams that generate subsurface head distributions. We recently completed a field study that suggested channel unit spacing in stream longitudinal profiles can be used to predict the spacing between zones of upwelling (flux of hyporheic water into the stream) and downwelling (flux of stream water into the hyporheic zone) in the beds of mountain streams. Here, we use two‐dimensional groundwater flow and particle tracking models to simulate vertical and longitudinal hyporheic exchange along the longitudinal axis of stream flow in second‐, third‐, and fourth‐order mountain stream reaches. Modelling allowed us to (1) represent visually the effect that the shape of the longitudinal profile has on the flow net beneath streambeds; (2) isolate channel unit sequence and spacing as individual factors controlling the depth that stream water penetrates the hyporheic zone and the length of upwelling and downwelling zones; (3) evaluate the degree to which the effects of regular patterns in bedform size and sequence are masked by irregularities in real streams. We simulated hyporheic exchange in two sets of idealized stream reaches and one set of observed stream reaches. Idealized profiles were constructed using regression equations relating channel form to basin area. The size and length of channel units (step size, pool length, etc.) increased with increasing stream order. Simulations of hyporheic exchange flows in these reaches suggested that upwelling lengths increased (from 2·7 m to 7·6 m), and downwelling lengths increased (from 2·9 m to 6·0 m) with increase in stream order from second to fourth order. Step spacing in the idealized reaches increased from 5·3 m to 13·7 m as stream size increased from second to fourth order. Simulated downwelling lengths increased from 4·3 m in second‐order streams to 9·7 m in fourth‐order streams with a POOL–RIFFLE–STEP channel unit sequence, and increased from 2·5 m to 6·1 m from second‐ to fourth‐order streams with a POOL–STEP–RIFFLE channel unit sequence. Upwelling lengths also increased with stream order in these idealized channels. Our results suggest that channel unit spacing, size, and sequence are all important in determining hyporheic exchange patterns of upwelling and downwelling. Though irregularities in the size and spacing of bedforms caused flow nets to be much more complex in surveyed stream reaches than in idealized stream reaches, similar trends emerged relating the average geomorphic wavelength to the average hyporheic wavelength in both surveyed and idealized reaches. This article replaces a previously published version (Hydrological Processes, 19 (17), 2915–2929 (2005) [ DOI:10.1002/hyp.5790 ]. See also retraction notice DOI:10.1002/hyp.6350 Copyright © 2006 John Wiley & Sons, Ltd.     

Retracted and replaced: A modelling study of hyporheic exchange pattern and the sequence,size, and spacing of stream bedforms in mountain stream networks,Oregon, USA

This article has been retracted and replaced. See Retraction and Replacement Notice DOI: 10.1002/hyp.6350 Studies of hyporheic exchange flows have identified physical features of channels that control exchange flow at the channel unit scale, namely slope breaks in the longitudinal profile of streams that generate subsurface head distributions. We recently completed a field study that suggested channel unit spacing in stream longitudinal profiles can be used to predict the spacing between zones of upwelling (flux of hyporheic water into the stream) and downwelling (flux of stream water into the hyporheic zone) in the beds of mountain streams. Here, we use two‐dimensional groundwater flow and particle tracking models to simulate vertical and longitudinal hyporheic exchange along the longitudinal axis of stream flow in second‐, third‐, and fourth‐order mountain stream reaches. Modelling allowed us to (1) represent visually the effect that the shape of the longitudinal profile has on the flow net beneath streambeds; (2) isolate channel unit sequence and spacing as individual factors controlling the depth that stream water penetrates the hyporheic zone and the length of upwelling and downwelling zones; (3) evaluate the degree to which the effects of regular patterns in bedform size and sequence are masked by irregularities in real streams. We simulated hyporheic exchange in two sets of idealized stream reaches and one set of observed stream reaches. Idealized profiles were constructed using regression equations relating channel form to basin area. The size and length of channel units (step size, pool length, etc.) increased with increasing stream order. Simulations of hyporheic exchange flows in these reaches suggested that upwelling lengths increased (from 2·7 m to 7·6 m), and downwelling lengths increased (from 2·9 m to 6·0 m) with increase in stream order from second to fourth order. Step spacing in the idealized reaches increased from 5·3 m to 13·7 m as stream size increased from second to fourth order. Simulated upwelling lengths increased from 4·3 m in second‐order streams to 9·7 m in fourth‐order streams with a POOL–RIFFLE–STEP channel unit sequence, and increased from 2·5 m to 6·1 m from second‐ to fourth‐order streams with a POOL–STEP–RIFFLE channel unit sequence. Downwelling lengths also increased with stream order in these idealized channels. Our results suggest that channel unit spacing, size, and sequence are all important in determining hyporheic exchange patterns of upwelling and downwelling. Though irregularities in the size and spacing of bedforms caused flow nets to be much more complex in surveyed stream reaches than in idealized stream reaches, similar trends emerged relating the average geomorphic wavelength to the average hyporheic wavelength in both surveyed and idealized reaches. Copyright © 2005 John Wiley & Sons, Ltd.     

The role of recharge and evapotranspiration as hydraulic drivers of ion concentrations in shallow groundwater on Everglades tree islands,Florida (USA)

Recently, evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year‐round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007 to 2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration‐to‐recharge rates were elevated, while low evapotranspiration‐to‐recharge rates resulted in a substantial increase of the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion, the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands. Copyright © 2012 John Wiley & Sons, Ltd.     

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).     

Beaver dams along an agricultural stream in southern Ontario,Canada: their impact on riparian zone hydrology and nitrogen chemistry

The hydrology and nitrogen biogeochemistry of a riparian zone were compared before and after the construction of beaver dams along an agricultural stream in southern Ontario, Canada. The beaver dams increased surface flooding and raised the riparian water table by up to 1·0 m. Increased hydraulic gradients inland from the stream limited the entry of oxic nitrate‐rich subsurface water from adjacent cropland. Permeable riparian sediments overlying dense till remained saturated during the summer and autumn months, whereas before dam construction a large area of the riparian zone was unsaturated in these seasons each year. Beaver dam construction produced significant changes in riparian groundwater chemistry. Median dissolved oxygen concentrations were lower in riparian groundwater after dam construction (0·9–2·1 mg L^?1) than in the pre‐dam period (2·3–3·9 mg L^?1). Median NO₃‐N concentrations in autumn and spring were also lower in the post‐dam (0·03–0·07 mg L^?1) versus the pre‐dam period (0·1–0·3 mg L^?1). In contrast, median NH₄‐N concentrations in autumn and spring months were higher after dam construction (0·3–0·4 mg L^?1) than before construction (0·13–0·14 mg L^?1). Results suggest that beaver dams can increase stream inflow to riparian areas that limit water table declines and increase depths of saturated riparian soils which become more anaerobic. These changes in subsurface hydrology and chemistry have the potential to affect the transport and transformation of nitrate fluxes from adjacent cropland in agricultural landscapes. Copyright © 2009 John Wiley & Sons, Ltd.     

Effect of temporal resolution of water level and temperature inputs on numerical simulation of groundwater–surface water flux exchange in a heavily modified urban river

Groundwater interacts with surface water features nearly in all types of landscapes. Understanding these interactions has practical consequences on the quantity and quality of water in either system, because the depletion or contamination of one of the systems will eventually affect the other one. Many studies have shown that the use of heat as natural tracer in conjunction with water level measurements is an effective method for estimating water flow (fluxes) between groundwater and surface water. A number of studies have explored the effects of spatial and temporal variability of groundwater–surface water flux exchanges using temperature and water level measurements; however, the effect of temporal resolution of water level and temperature data on estimating flux remains unexplored. Therefore, this study investigated the effect of temporal resolution of input data on temporal variation of groundwater–surface water flux exchanges. To this end, we calibrated a variably saturated two‐dimensional groundwater flow and heat transport model (VS2DH) at hourly and daily time scales using temperatures measured at multiple depths below the riverbed of the Zenne River, located at a well‐known Belgian brownfield site. Results of the study showed that the computed water flux through the streambed ranged between ?32 mm/day and +25 mm/day using the hourly model and from ?10 mm/day to ?37 mm/day using the daily model. The hourly model resulted in detecting reversal of flow direction inducing short‐term surface water flow into the streambed. However, such events were not captured if daily temperature and water level measurements were used as input. These findings have important implications for understanding contaminant mass flux and their attenuation in the mixing zone of groundwater and surface water. Copyright © 2012 John Wiley & Sons, Ltd.     

Dissolution in a variably confined carbonate platform: effects of allogenic runoff,hydraulic damming of groundwater inputs,and surface–groundwater exchange at the basin scale

In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SI_CAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.     

The importance of incorporating diurnally fluctuating stream discharge in stream temperature energy balance models

Although stream temperature energy balance models are useful to predict temperature through time and space, a major unresolved question is whether fluctuations in stream discharge reduce model accuracy when not exactly represented. However, high‐frequency (e.g., subdaily) discharge observations are often unavailable for such simulations, and therefore, diurnal streamflow fluctuations are not typically represented in energy balance models. These fluctuations are common due to evapotranspiration, snow pack or glacial melt, tidal influences within estuaries, and regulated river flows. In this work, we show when to account for diurnally fluctuating streamflow. To investigate how diurnal streamflow fluctuations affect predicted stream temperatures, we used a deterministic stream temperature model to simulate stream temperature along a reach in the Quilcayhuanca Valley, Peru, where discharge varies diurnally due to glacial melt. Diurnally fluctuating streamflow was varied alongside groundwater contributions via a series of computational experiments to assess how uncertainty in reach hydrology may impact simulated stream temperature. Results indicated that stream temperatures were more sensitive to the rate of groundwater inflow to the reach compared with the timing and amplitude of diurnal fluctuations in streamflow. Although incorporating observed diurnal fluctuations in discharge resulted in a small improvement in model RMSE, we also assessed other diurnal discharge signals and found that high amplitude signals were more influential on modelled stream temperatures when the discharge peaked at specific times. Results also showed that regardless of the diurnal discharge signal, the estimated groundwater flux to the reach only varied from 1.7% to 11.7% of the upstream discharge. However, diurnal discharge fluctuations likely have a stronger influence over longer reaches and in streams where the daily range in discharge is larger, indicating that diurnal fluctuations in stream discharge should be considered in certain settings.