Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ¹⁸O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.     

Using geophysical surveys to test tracer‐based storage estimates in headwater catchments

Hydrogeophysical surveys were carried out in a 3.2 km² Scottish catchment where previous isotope studies inferred significant groundwater storage that makes important contributions to streamflow. We used electrical resistivity tomography (ERT) to characterize the architecture of glacial drifts and make an approximation of catchment‐scale storage. Four ERT lines (360–535 m in length) revealed extensive 5–10 m deep drift cover on steeper slopes, which extends up to 20–40 m in valley bottom areas. Assuming low clay fractions, we interpret variable resistivity as correlating with variations in porosity and water content. Using Archie's Law as a first approximation, we compute likely bounds for storage along the ERT transects. Areas of highest groundwater storage occur in valley bottom peat soils (up to 4 m deep) and underlying drift where up to 10 000 mm of precipitation equivalent may be stored. This is consistent with groundwater levels which indicate saturation to within 0.2 m of the surface. However, significant slow groundwater flow paths occur in the shallower drifts on steeper hillslopes, where point storage varies between ~1000 mm–5000 mm. These fluxes maintain saturated conditions in the valley bottom and are recharged from drift‐free areas on the catchment interfluves. The surveys indicate that catchment scale storage is >2000 mm which is consistent with tracer‐based estimates. Copyright © 2016 John Wiley & Sons, Ltd.     

Characterizing hydrological processes in a semiarid rangeland watershed: A hydrogeophysical approach

The complex ecohydrological processes of rangelands can be studied through the framework of ecological sites (ESs) or hillslope‐scale soil–vegetation complexes. High‐quality hydrologic field investigations are needed to quantitatively link ES characteristics to hydrologic function. Geophysical tools are useful in this context because they provide valuable information about the subsurface at appropriate spatial scales. We conducted 20 field experiments in which we deployed time‐lapse electrical resistivity tomography (ERT), variable intensity rainfall simulation, ground‐penetrating radar (GPR), and seismic refraction, on hillslope plots at five different ESs within the Upper Crow Creek Watershed in south‐east Wyoming. Surface runoff was measured using a precalibrated flume. Infiltration data from the rainfall simulations, coupled with site‐specific resistivity–water content relationships and ERT datasets, were used to spatially and temporally track the progression of the wetting front. First‐order constraints on subsurface structure were made at each ES using the geophysical methods. Sites ranged from infiltrating 100% of applied rainfall to infiltrating less than 60%. Analysis of covariance results indicated significant differences in the rate of wetting front progression, ranging from 0.346 m min^?1/2 for sites with a subsurface dominated by saprolitic material to 0.156 m min^?1/2 for sites with a well‐developed soil profile. There was broad agreement in subsurface structure between the geophysical methods with GPR typically providing the most detail. Joint interpretation of the geophysics showed that subsurface features such as soil layer thickness and the location of subsurface obstructions such as granite corestones and material boundaries had a large effect on the rate of infiltration and subsurface flow processes. These features identified through the geophysics varied significantly by ES. By linking surface hydrologic information from the rainfall simulations with subsurface information provided by the geophysics, we can characterize the ES‐specific hydrologic response. Both surface and subsurface flow processes differed among sites and are directly linked to measured characteristics.     

A case study in catchment hydrochemistry: Conflicting interpretations from hydrological and chemical observations

Soil water, stream water, groundwater and rain water were sampled through a storm event in a moorland catchment. Samples were analysed for major ions and deuterium. Chloride and deuterium are used as tracers to enable separation of the stream runoff hydrograph into three components: rain, soil and groundwater. The results indicate that rain water arrives in the stream quickly during the event and contributes a significant volume to the runoff peak. The chemical signal in the rain water is, however, significantly damped, apparently due to mixing with soil water held in the catchment before the event. This is further modified before reaching the stream, apparently through mixing with a deeper groundwater component. Interpretation of tracer, chemistry and hydrological data to present an integrated picture of catchment hydrochemical response is difficult due to problems in the chemical and conceptual definition of the flow components.     

INFERRING HYDROLOGICAL PROCESSES IN A TEMPERATE BASIN USING ISOTOPIC AND GEOCHEMICAL HYDROGRAPH SEPARATION: A RE-EVALUATION

Simultaneous monitoring of conservative and non-conservative tracers in streamflow offers a valuable means of obtaining information on the age and flow paths of water reaching the basin outlet. Previous studies of stormflow generation in a small forested basin on the Canadian Shield used isotopic (IHS) and geochemical hydrograph separations (GHS) to infer that some event water during snowmelt reaches the stream via subsurface pathways, and that surface water runoff is generated by direct precipitation on to saturated areas (DPSA) in the stream valley. These hypotheses were tested for rainfall inputs using simultaneous IHS (¹⁸O) and GHS (dissolved silica) of basin stormflow, supplemented by hydrochemical and hydrometric data from throughflow troughs installed on basin slopes. Comparison of pre-event and subsurface water hydrographs did not provide conclusive evidence for subsurface movement of event water to the stream, owing to the appreciable uncertainty associated with the hydrograph separations. However, IHSs of runoff at the soil–bedrock interface on basin slopes indicated that event water comprised 25–50% of total runoff from areas with deep soil cover, and that these contributions supplied event water flux from the basin in excess of that attributable to DPSA. The surface water component of stormflow estimated from the GHS was also largely the result of DPSA. GHS assumes that dissolved silica is rapidly and uniformly taken up by water infiltrating the soil and that water moving via surface pathways retains the low dissolved silica level of rainfall; however, neither assumption was supported by the hillslope results. Instead, results suggest that the observed depression of silica levels in basin stormflow previously attributed to dilution by DPSA was partly a function of transport of dilute event water to the channel via preferential pathways. Implications of these results for the general use of simultaneous IHS and GHS to infer hydrological processes are discussed. © 1997 by John Wiley & Sons, Ltd.     

Experimental evidence of fast groundwater responses in a hillslope/floodplain area in the Black Forest Mountains,Germany

Analysis of water flow pathways from hillslopes to streams is essential for the optimal protection of water resources as well as for ecohydrological studies. This study addresses runoff generation processes at a hillslope and near‐stream shallow groundwater system in the Black Forest Mountains, southwestern Germany. The changing spatial and temporal flow patterns during differing hydrological situations were examined using a combined hydraulic and hydrochemical approach. Groundwater levels at 10 wells, discharge at a near‐stream saturated area, and several natural tracers (deuterium, dissolved silica, and major anions and cations) were observed at different locations during high and low flows. The importance of the groundwater component during flood formation was clearly demonstrated: its contribution was about 80% during a double peak flood event at the saturated area. In addition, a rapid change of the shallow groundwater levels was observed along two transects of groundwater wells in the floodplain. This led to an enhanced groundwater discharge into the saturated area located at the end of one study transect. The amount of groundwater additionally activated during the event was about 30% of total discharge recorded at the outlet of the saturated area. Two alternative hypotheses are discussed to explain this phenomenon: the establishment of locally confined conditions and the development of a pressure wave (hypothesis A), or the significant change of the three‐dimensional groundwater flow lines that caused a large increase of the groundwater catchment at the saturated area during the investigated event (hypothesis B). Even if the exact flow paths and mechanisms could not be clearly identified, the importance of rapid responding hillslope groundwater was undoubtedly demonstrated by a combination of tracer and hydrometric methods. Copyright © 2004 John Wiley & Sons, Ltd.     

Paired‐basin comparison of hydrological response in harvested and undisturbed hardwood forests during snowmelt in central Ontario: I. Streamflow,groundwater and flowpath behaviour

Impacts of forest harvesting on groundwater properties, water flowpaths and streamflow response were examined 4 years after the harvest using a paired‐basin approach during the 2001 snowmelt in a northern hardwood landscape in central Ontario. The ability of two metrics of basin topography (Beven and Kirkby's ln(a/tan β) topographic index (TI) and distance to stream channel) to explain intra‐basin variations in groundwater dynamics was also evaluated. Significant relationships between TI and depth to potentiometric surface for shallow groundwater emerged, although the occurrence of these relationships during the melt differed between harvested and control basins, possibly as a result of interbasin differences in upslope area contributing to piezometers used to monitor groundwater behaviour. Transmissivity feedback (rapid streamflow increases as the water table approaches the soil surface) governed streamflow generation in both basins, and the mean threshold depths at which rapid streamflow increases corresponded to small rises in water level were similar for harvested (0·41 ± 0·05 m) and forested (0·38 ± 0·04 m) basins. However, topographic properties provided inconsistent explanations of spatial variations in the relationship between streamflow and depth to water at a given piezometer for both basins. Streamflow from the harvested basin exceeded that from the forested basin during the 2001 melt, and hydrometric and geochemical tracer results indicated greater runoff from the harvested basin via surface and near‐surface pathways. These differences are not solely attributable to harvesting, since the difference in spring runoff from the harvested basin relative to the forested control was not consistently larger than under pre‐harvest conditions. Nevertheless, greater melt rates following harvesting appear to have increased the proportion of water delivery to the stream channel via surface and near‐surface pathways. Copyright © 2005 John Wiley & Sons, Ltd.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Geoelectrical and hydrochemical investigations for characterizing the salt water intrusion in the Khanasser valley, northern Syria

An integrated approach of geoelectrical and hydrochemical investigation surveys was proposed for indicating contact regions between saline and fresh groundwater in the Khanasser valley region, northern Syria. The qualitative and quantitative interpretations of 34 vertical electrical soundings (VES) enable to characterize the salt water intrusion laterally and vertically. The established iso-apparent resistivity maps for different AB/2 spacings obviously indicate the presence of a lowresistivity (less than 4 Ohm·m) zone related to the salt water intrusion in the Quaternary and Paleogene deposits. The different hydrochemical and geophysical parameters, such as electrical resistivity, total dissolved solids (TDS) and major ions concentrations used to characterize the salt water intrusion gave almost similar results in locating and mapping the different boundaries of the groundwater salinity. The proposed approach is useful for mapping the interface between different groundwater qualities, and can be therefore used to successfully characterize the salt water intrusion phenomenon in other semi-arid regions. The application of such an approach is a powerful tool and can be used for water resource management in the water scarce areas.     

Time‐Lapse Electrical Geophysical Monitoring of Amendment‐Based Biostimulation

Biostimulation is increasingly used to accelerate microbial remediation of recalcitrant groundwater contaminants. Effective application of biostimulation requires successful emplacement of amendment in the contaminant target zone. Verification of remediation performance requires postemplacement assessment and contaminant monitoring. Sampling‐based approaches are expensive and provide low‐density spatial and temporal information. Time‐lapse electrical resistivity tomography (ERT) is an effective geophysical method for determining temporal changes in subsurface electrical conductivity. Because remedial amendments and biostimulation‐related biogeochemical processes often change subsurface electrical conductivity, ERT can complement and enhance sampling‐based approaches for assessing emplacement and monitoring biostimulation‐based remediation. Field studies demonstrating the ability of time‐lapse ERT to monitor amendment emplacement and behavior were performed during a biostimulation remediation effort conducted at the Department of Defense Reutilization and Marketing Office (DRMO) Yard, in Brandywine, Maryland, United States. Geochemical fluid sampling was used to calibrate a petrophysical relation in order to predict groundwater indicators of amendment distribution. The petrophysical relations were field validated by comparing predictions to sequestered fluid sample results, thus demonstrating the potential of electrical geophysics for quantitative assessment of amendment‐related geochemical properties. Crosshole radar zero‐offset profile and borehole geophysical logging were also performed to augment the data set and validate interpretation. In addition to delineating amendment transport in the first 10 months after emplacement, the time‐lapse ERT results show later changes in bulk electrical properties interpreted as mineral precipitation. Results support the use of more cost‐effective surface‐based ERT in conjunction with limited field sampling to improve spatial and temporal monitoring of amendment emplacement and remediation performance.     

Defining hydrochemical evolution of streamflow through flowpath dynamics in Kawakami headwater catchment,Central Japan

The hydrochemical behaviour of catchments is often investigated by inferring stream chemistry through identification of source areas involved in hydrograph separation analysis, yet its dynamic evolution of hydrologic pathways has received little attention. Intensive hydrometric and hydrochemical measurements were performed during two different storms on March 29, 2001 and August 21–22, 2001 to define hydrochemical evolution under the dynamic of flow pathways in a 5·2 ha first‐order drainage of the Kawakami experimental basin (KEB), Central Japan, a forested headwater catchment with various soil depths (1·8 to 5 m) overlying late Neogene of volcanic bedrocks. The hydraulic potential distribution and flow lines data showed that the change in flow direction, which was controlled by rainfall amount and antecedent wetness of the soil profile, agreed well with the hydrochemical change across the slope segment during the storm. Hydrograph separation predicted by end‐member mixing analysis (EMMA) using Ca²⁺ and SiO₂ showed that near surface riparian, hillslope soil water and deep riparian groundwater were important in stream flow generation. The evidence of decrease in solutes concentration at a depth of 1 m in the hillslope and 0·6 m in the near surface riparian during peak storm suggested a flushing of high solutes concentration. Most of the solutes accumulated in the deep riparian groundwater zone, which was due to prominent downward flow and agreed well with the residence time. The distinct flow pathways and chemistry between the near surface riparian and deep riparian groundwater zones and the linkage hillslope aquifer and near surface riparian reservoir, which controls rapid flow and solutes flushing during the storm event, are in conflict with the typical assumption that the whole riparian zone resets flow pathways and chemical signature of hillslope soil water, as has been reported in a previous study. Copyright © 2005 John Wiley & Sons, Ltd.     

Investigating mechanisms of stormflow generation by natural tracers and hydrometric data: a small catchment study in the Black Forest,Germany

The importance and interaction of various hydrological pathways and identification of runoff source areas involved in solute transport are still under considerable debate in catchment hydrology. To reveal stormflow generating areas and flow paths, hydrometric behaviour of throughfall, soil water from various depths, runoff, and respective concentrations of the environmental tracers ¹⁸O, Si, K, SO₄ and dissolved organic carbon were monitored for a 14‐week period in a steep headwater catchment in the Black Forest Mountains, Germany. Two stormflow hydrographs were selected and, based on ¹⁸O and Si, chemically separated into three flow components. Their sources were defined using mixing diagrams. Additional information about stormflow generating mechanisms was derived from recession analyses of the basin's complete 5‐year hydrograph record. By providing insight into storage properties and residence times of outflowing reservoirs of the basin, recession analysis proved to be a valuable tool in runoff model conceptualization. Its results agreed well with hydrometric and hydrochemical data. Supported by evaluation of 30 hillslope soil profiles a coherent concept of stormflow generation could be derived: whereas in many steeply sloped basins in the temperate region soil water from hillslopes appears to have an immediate effect on the shape of the stormflow hydrograph, its role at this basin is basically restricted to the recharge of the groundwater reservoir in the near‐channel area. Storm hydrograph peaks appear to be derived from a small direct runoff component supplemented by a fast delivery of baseflow from the groundwater reservoir in the valley bottom. Copyright © 2001 John Wiley & Sons, Ltd.     

Electrical resistivity imaging of hydrologic flow through surface coal mine valley fills with comparison to other landforms

Surface coal mining has altered land cover, near‐surface geologic structure, and hydrologic processes of large areas in central Appalachia, USA. These alterations are associated with changes in water quality such as elevated total‐dissolved solids, which is usually measured via its surrogate, specific conductance (SC). The SC of valley fill effluent streams is a function of fill construction methods, materials, and age; yet hydrologic studies that relate these variables to water quality are sparse due to the difficulty of conducting traditional hydrologic studies in mined landscapes. We used electrical resistivity imaging (ERI) to visualize the subsurface geologic structure and hydrologic flow paths within a valley fill. ERI is a noninvasive geophysical technique that maps spatiotemporal changes in resistivity of the subsurface. We paired ERI with artificial rainfall experiments to track infiltrated water as it moved through the valley fill. Results indicate that ERI can be used to identify subsurface geologic structure and track advancing wetting fronts or preferential flow paths. Our results suggest that the upper portion of the fill contains significant fines, whereas the deeper profile is primarily large rocks and void spaces. Water tended to pond on the surface of compacted areas until it reached preferential flow paths, where it appeared to infiltrate quickly down to >15 m depth in 75 min. ERI applications can improve understanding of how fill construction techniques influence subsurface water movement, and in turn may aid in the development of valley fill construction methods to reduce water quality effects.     

Testing a spatially distributed tracer‐aided runoff model in a snow‐influenced catchment: Effects of multicriteria calibration on streamwater ages

Integrating stable isotope tracers into rainfall‐runoff models allows investigation of water partitioning and direct estimation of travel times and water ages. Tracer data have valuable information content that can be used to constrain models and, in integration with hydrometric observations, test the conceptualization of catchment processes in model structure and parameterization. There is great potential in using tracer‐aided modelling in snow‐influenced catchments to improve understanding of these catchments' dynamics and sensitivity to environmental change. We used the spatially distributed tracer‐aided rainfall‐runoff (STARR) model to simulate the interactions between water storage, flux, and isotope dynamics in a snow‐influenced, long‐term monitored catchment in Ontario, Canada. Multiple realizations of the model were achieved using a combination of single and multiple objectives as calibration targets. Although good simulations of hydrometric targets such as discharge and snow water equivalent could be achieved by local calibration alone, adequate capture of the stream isotope dynamics was predicated on the inclusion of isotope data in the calibration. Parameter sensitivity was highest, and most local, for single calibration targets. With multiple calibration targets, key sensitive parameters were still identifiable in snow and runoff generation routines. Water ages derived from flux tracking subroutines in the model indicated a catchment where runoff is dominated by younger waters, particularly during spring snowmelt. However, resulting water ages were most sensitive to the partitioning of runoff sources from soil and groundwater sources, which was most realistically achieved when isotopes were included in the calibration. Given the paucity of studies where hydrological models explicitly incorporate tracers in snow‐influenced regions, this study using STARR is an important contribution to satisfactorily simulating snowpack dynamics and runoff generation processes, while simultaneously capturing stable isotope variability in snow‐influenced catchments.     

Improved simulation of river water and groundwater exchange in an alluvial plain using the SWAT model

Hydrological interaction between surface and subsurface water systems has a significant impact on water quality, ecosystems and biogeochemistry cycling of both systems. Distributed models have been developed to simulate this function, but they require detailed spatial inputs and extensive computation time. The soil and water assessment tool (SWAT) model is a semi‐distributed model that has been successfully applied around the world. However, it has not been able to simulate the two‐way exchanges between surface water and groundwater. In this study, the SWAT‐landscape unit (LU) model – based on a catena method that routes flow across three LUs (the divide, the hillslope and the valley) – was modified and applied in the floodplain of the Garonne River. The modified model was called SWAT‐LUD. Darcy's equation was applied to simulate groundwater flow. The algorithm for surface water‐level simulation during flooding periods was modified, and the influence of flooding on groundwater levels was added to the model. Chloride was chosen as a conservative tracer to test simulated water exchanges. The simulated water exchange quantity from SWAT‐LUD was compared with the output of a two‐dimensional distributed model, surface–subsurface water exchange model. The results showed that simulated groundwater levels in the LU adjoining the river matched the observed data very well. Additionally, SWAT‐LUD model was able to reflect the actual water exchange between the river and the aquifer. It showed that river water discharge has a significant influence on the surface–groundwater exchanges. The main water flow direction in the river/groundwater interface was from groundwater to river; water that flowed in this direction accounted for 65% of the total exchanged water volume. The water mixing occurs mainly during high hydraulic periods. Flooded water was important for the surface–subsurface water exchange process; it accounted for 69% of total water that flowed from the river to the aquifer. The new module also provides the option of simulating pollution transfer occurring at the river/groundwater interface at the catchment scale. Copyright © 2015 John Wiley & Sons, Ltd.     

Conceptual modelling to assess how the interplay of hydrological connectivity,catchment storage and tracer dynamics controls nonstationary water age estimates

Although catchment storage is an intrinsic control on the rainfall–runoff response of streams, direct measurement remains a major challenge. Coupled models that integrate long‐term hydrometric and isotope tracer data are useful tools that can provide insights into the dynamics of catchment storage and the volumes of water involved. In this study, we use a tracer‐aided hydrological model to characterize catchment storage as a dynamic control on system function related to streamflow generation, which also allows direct estimation of the nonstationarity of water ages. We show that in a wet Scottish upland catchment dominated by runoff generation from riparian peats (histosols) with high water storage, nonstationarity in water age distributions is only clearly detectable during more extreme wet and dry periods. This is explained by the frequency and longevity of hydrological connectivity and the associated relative importance of flow paths contributing younger or older waters to the stream. Generally, these saturated riparian soils represent large mixing zones that buffer the time variance of water age and integrate catchment‐scale partial mixing processes. Although storage simulations depend on model performance, which is influenced by input variability and the degree of isotopic damping in the stream, a longer‐term storage analysis of this model indicates a system that is only sensitive to more extreme hydroclimatic variability. Copyright © 2014 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.     

Contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments

We examined the contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments by conducting detailed hydrochemical observations in five catchments that ranged from zero to second order. End‐member mixing analysis (EMMA) was performed to identify the geographical sources of stream water. Throughfall, hillslope groundwater, shallow bedrock groundwater, and deep bedrock groundwater were identified as end members. The contribution of each end member to storm runoff differed among the catchments because of the differing quantities of riparian groundwater, which was recharged by the bedrock groundwater prior to rainfall events. Among the five catchments, the contribution of throughfall was highest during both baseflow and storm flow in a zero‐order catchment with little contribution from the bedrock groundwater to the riparian reservoir. In zero‐order catchments with some contribution from bedrock groundwater, stream water was dominated by shallow bedrock groundwater during baseflow, but it was significantly influenced by hillslope groundwater during storms. In the first‐order catchment, stream water was dominated by shallow bedrock groundwater during storms as well as baseflow periods. In the second‐order catchment, deeper bedrock groundwater than that found in the zero‐order and first‐order catchments contributed to stream water in all periods, except during large storm events. These results suggest that bedrock groundwater influences the upscaling of storm‐runoff generation processes by affecting the linkages of geomorphic units such as hillslopes, riparian zones, and stream channels. Our results highlight the need for a three‐dimensional approach that considers bedrock groundwater flow when studying the upscaling of storm‐runoff generation processes. Copyright © 2014 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.