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.     

Effective groundwater–surface water exchange at watershed scales

Coupled groundwater–surface water (GW–SW) models are capable of simulating complex hydrological systems when used at fine resolutions. However, properly characterizing bulk GW–SW fluxes for either coarsely resolved integrated models or basin‐discretized surface water models remains a challenge. Loss of subgrid detail, while beneficially decreasing computational cost, leads to a decrease in model accuracy as scale effects become important. Ideally, coarse low‐resolution models should be informed by expected subgrid behaviour, reducing the impact of scale effects. Determining how to best represent these fine‐scale details in lower‐resolution models is important for improving the accuracy and appropriateness of these models. To investigate some of these scale effects, we here explore the relationships between area‐averaged hydraulic head and bulk GW–SW exchange fluxes (e.g. evapotranspiration and discharge), all of which are presumed to be controlled predominantly by subgrid topographic effects. These relationships may be useful for simply upscaling models without the complete loss of crucial fine‐resolution subgrid details. Using finely resolved simulation output from Modflow for a fine‐resolution simulation and post‐processed results generated to represent coarser resolutions, upscaled flux relationships (UFRs) are generated for multiple terrains; these UFRs define the relationships that exist between average hydraulic head and average fluxes in unconfined aquifer systems. It is found that, for steady‐flow regimes, similar one‐to‐one power law relationships consistently exist between area‐averaged hydraulic heads, exchange fluxes and saturated area for a variety of terrains. Additionally, when the averaged values are properly normalized, the generated steady‐state UFRs for a single terrain are independent of hydraulic conductivity and potential evapotranspiration rates and apparently insensitive to the presence of mild heterogeneity. While some hysteresis is apparent in the relationships under transient conditions, transient artefacts are shown to be minor under some circumstances, indicating that UFRs may be applied to both steady‐state and transient scenarios. Simpler tests performed under saturated and variably saturated conditions in a cross‐sectional model show similar trends, suggesting that the UFR representation is extendable to systems where the vadose zone plays a significant role. It is suggested that relatively simple UFRs such as these may find use as an alternative to direct point upscaling or multi‐resolution models for estimating GW–SW exchange fluxes in coarse‐scale models. They also appear to justify the functional form of some classical models of baseflow and evapotranspiration used in conceptual surface water models. Copyright © 2015 John Wiley & Sons, Ltd.     

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

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

The impact of groundwater–surface water interactions on the water balance of a mesoscale lowland river catchment in northeastern Germany

The glacially formed northeastern German lowlands are characterized by extensive floodplains, often interrupted by relatively steep moraine hills. The hydrological cycle of this area is governed by the tight interaction of surface water dynamics and the corresponding directly connected shallow groundwater aquifer. Runoff generation processes, as well as the extent and spatial distribution of the interaction between surface water and groundwater, are controlled by floodplain topography and by surface water dynamics. A modelling approach based on extensive experimental analyses is presented that describes the specific water balance of lowland areas, including the interactions of groundwater and surface water, as well as reflecting the important role of time‐variable shallow groundwater stages for runoff generation in floodplains. In the first part, experimental investigations of floodplain hydrological characteristics lead to a qualitative understanding of the water balance processes and to the development of a conceptual model of the water balance and groundwater dynamics of the study area. Thereby model requirements which allow for an adequate simulation of the floodplain hydrology, considering also interactions between groundwater and surface water have been characterized. Based on these analyses, the Integrated Modelling of Water Balance and Nutrient Dynamics (IWAN) approach has been developed. This consists of coupling the surface runoff generation and soil water routines of the deterministic, spatially distributed hydrological model WASIM‐ETH‐I with the three‐dimensional finite‐difference‐based numerical groundwater model MODFLOW and Processing MODFLOW. The model was applied successfully to a mesoscale subcatchment of the Havel River in northeast Germany. It was calibrated for two small catchments (1·4 and 25 km²), where the importance of the interaction processes between groundwater and surface waters and the sensitivity of several controlling parameters could be quantified. Validation results are satisfying for different years for the entire 198 km² catchment. The model approach was further successfully tested for specific events. The experimental area is a typical example of a floodplain‐dominated landscape. It was demonstrated that the lateral flow processes and the interactions between groundwater and surface water have a major importance for the water balance and periodically superimposed on the vertical runoff generation. Copyright © 2006 John Wiley & Sons, Ltd.     

Diffusive equilibrium in thin films provides evidence of suppression of hyporheic exchange and large‐scale nitrate transformation in a groundwater‐fed river

The hyporheic zone of riverbed sediments has the potential to attenuate nitrate from upwelling, polluted groundwater. However, the coarse‐scale (5–10 cm) measurement of nitrogen biogeochemistry in the hyporheic zone can often mask fine‐scale (<1 cm) biogeochemical patterns, especially in near‐surface sediments, leading to incomplete or inaccurate representation of the capacity of the hyporheic zone to transform upwelling NO₃^?. In this study, we utilised diffusive equilibrium in thin‐films samplers to capture high resolution (cm‐scale) vertical concentration profiles of NO₃^?, SO₄^2?, Fe and Mn in the upper 15 cm of armoured and permeable riverbed sediments. The goal was to test whether nitrate attenuation was occurring in a sub‐reach characterised by strong vertical (upwelling) water fluxes. The vertical concentration profiles obtained from diffusive equilibrium in thin‐films samplers indicate considerable cm‐scale variability in NO₃^? (4.4 ± 2.9 mg N/L), SO₄^2? (9.9 ± 3.1 mg/l) and dissolved Fe (1.6 ± 2.1 mg/l) and Mn (0.2 ± 0.2 mg/l). However, the overall trend suggests the absence of substantial net chemical transformations and surface‐subsurface water mixing in the shallow sediments of our sub‐reach under baseflow conditions. The significance of this is that upwelling NO₃^?‐rich groundwater does not appear to be attenuated in the riverbed sediments at <15 cm depth as might occur where hyporheic exchange flows deliver organic matter to the sediments for metabolic processes. It would appear that the chemical patterns observed in the shallow sediments of our sub‐reach are not controlled exclusively by redox processes and/or hyporheic exchange flows. Deeper‐seated groundwater fluxes and hydro‐stratigraphy may be additional important drivers of chemical patterns in the shallow sediments of our study sub‐reach. © 2015 The Authors. Hydrological Processes Published by 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.     

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.     

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.     

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.     

Spatial patterns of groundwater‐lake exchange – implications for acid neutralization processes in an acid mine lake

Exchange of groundwater and lake water with typically quite different chemical composition is an important driver for biogeochemical processes at the groundwater‐lake interface, which can affect the water quality of lakes. This is of particular relevance in mine lakes where anoxic and slightly acidic groundwater mixes with oxic and acidic lake water (pH < 3). To identify links between groundwater‐lake exchange rates and acid neutralization processes in the sediments, exchange rates were quantified and related to pore‐water pH, sulfate and iron concentrations as well as sulfate reduction rates within the sediment. Seepage rates measured with seepage meters (?2.5 to 5.8 L m^‐2 d^‐1) were in reasonable agreement with rates inverted from modeled chloride profiles (?1.8 to 8.1 L m^‐2 d^‐1). Large‐scale exchange patterns were defined by the (hydro)geologic setting but superimposed by smaller scale variations caused by variability in sediment texture. Sites characterized by groundwater upwelling (flow into the lake) and sites where flow alternated between upwelling and downwelling were identified. Observed chloride profiles at the alternating sites reflected the transient flow regime. Seepage direction, as well as seepage rate, were found to influence pH, sulfate and iron profiles and the associated sulfate reduction rates. Under alternating conditions proton‐consuming processes, for example, sulfate reduction, were slowed. In the uppermost layer of the sediment (max. 5 cm), sulfate reduction rates were significantly higher at upwelling (>330 nmol g^‐1 d^‐1) compared to alternating sites (<220 nmol g^‐1 d^‐1). Although differences in sulfate reduction rates could not be explained solely by different flux rates, they were clearly related to the prevailing groundwater‐lake exchange patterns and the associated pH conditions. Our findings strongly suggest that groundwater‐lake exchange has significant effects on the biogeochemical processes that are coupled to sulfate reduction such as acidity retention and precipitation of iron sulfides. Copyright © 2012 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.     

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.     

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.     

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.     

Geomorphologic control on pollutant retardation at the groundwater–surface water interface

The results of research on the pollutant retardation potential of permeable riverbed sediments in catchments with significant groundwater–surface water (GW-SW) interaction are presented. The fraction of organic carbon and cation exchange capacity of fluvial sediments in various geomorphologic environments have been quantified. Sediments in selected reaches of the rivers Tern and Leith (UK), from the underlying Permian sandstone aquifers, and from along the length of the rivers Severn and Eden into which the Tern and Leith discharge have been investigated. Statistical analyses show significant variation in the geochemistry and pollutant retardation potential of sediments from different geomorphologic features, and between upland and lowland rivers. The sorption potential of fine-grained sediments deposited in pools was greater than sand in runs and coarser deposits in riffles. Similarly, sediments in lowland rivers were found to have a greater retardation potential than those in upland rivers. There was generally greater retardation potential in fluvial sediments of all types than in the underlying aquifers, and in lowland rivers the fluvial sediment retardation potential greatly dominated that of the aquifer. The findings demonstrate the potential for pollutant retardation processes in riverbed sediments of sandstone catchments, and suggest that consideration of retardation processes at the groundwater–surface water interface should be included into environmental risk-assessment studies, in order to better assess and manage the effects of contaminated groundwater discharges to rivers, particularly in lowland catchments. Copyright © 2008 John Wiley & Sons, Ltd.     

Groundwater–surface water interactions in a lowland watershed: source contribution to stream flow

The lower coastal plain of the Southeast USA is undergoing rapid urbanisation as a result of population growth. Land use change has been shown to affect watershed hydrology by altering stream flow and, ultimately, impairing water quality and ecologic health. However, because few long‐term studies have focused on groundwater–surface water interactions in lowland watersheds, it is difficult to establish what the effect of development might be in the coastal plain region. The objective of this study was to use an innovative improvement to end‐member mixing analysis (EMMA) to identify time sequences of hydrologic processes affecting storm flow. Hydrologic and major ion chemical data from groundwater, soil water, precipitation and stream sites were collected over a 2‐year period at a watershed located in USDA Forest Service's Santee Experimental Forest near Charleston, South Carolina, USA. Stream flow was ephemeral and highly dependent on evapotranspiration rates and rainfall amount and intensity. Hydrograph separation for a series of storm events using EMMA allowed us to identify precipitation, riparian groundwater and streambed groundwater as main sources to stream flow, although source contribution varied as a function of antecedent soil moisture condition. Precipitation, as runoff, dominated stream flow during all storm events while riparian and streambed groundwater contributions varied and were mainly dependent on antecedent soil moisture condition. Sensitivity analyses examined the influence of 10% and 50% increases in analyte concentration on EMMA calculations and found that contribution estimates were very sensitive to changes in chemistry. This study has implications on the type of methodology used in traditional forms of EMMA research, particularly in the recognition and use of median end‐member water chemistry in hydrograph separation techniques. Potential effects of urban development on important hydrologic processes (groundwater recharge, interflow, runoff, etc.) that influence stream flow in these lowland watersheds were qualitatively examined. Copyright © 2011 John Wiley & Sons, Ltd.     

High frequency stream bed mobility of a low‐gradient agricultural stream with implications on the hyporheic zone

Little Kickapoo Creek (LKC), a low‐gradient stream, mobilizes its streambed–fundamentally altering its near‐surface hyporheic zone–more frequently than do higher‐gradient mountain and karst streams. LKC streambed mobility was assessed through streambed surveys, sediment sampling, and theoretical calculations comparing basal shear stress (τ_b) with critical shear stress (τ_c). Baseflow τ_b is capable of entraining a d₅₀ particle; bankfull flow could entrain a 51·2 mm particle. No particle that large occurs in the top 30 cm of the substrate, suggesting that the top 30 cm of the substrate is mobilized and redistributed during bankfull events. Bankfull events occur on average every 7·6 months; flows capable of entraining d₅₀ and d₈₅ particles occur on average every 0·85 and 2·1 months, respectively. Streambed surveys verify streambed mobility at conditions below bankfull. While higher gradient streams have higher potential energy than LKC, they achieve streambed‐mobilization thresholds less frequently. Heterogeneous sediment redistribution creates an environment where substrate hydraulic conductivity (K) varies over four orders of magnitude. The frequency and magnitude of the substrate entrainment has implications on hyporheic zone function in fluid, solute and thermal transport models, interpretations of hyporheic zone stability, and understanding of LKC's aquatic ecosystem. Copyright © 2008 John Wiley & Sons, Ltd.