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.     

Influence of lowland aquifers and anthropogenic impacts on the isotope hydrology of contrasting mesoscale catchments

We examined the isotope hydrology of eight, contrasting mesoscale (104–488 km²) catchments characterized by a systematic change in the relative importance of upland and lowland areas that reflects the relative distribution of metamorphic and sedimentary rocks. Precipitation and stream water were monitored over a 12‐month period, and stable isotopes were used to examine spatial variations in the hydrometric and tracer dynamics of the catchments. Isotopic tracers were used to examine the temporal dynamics of different runoff sources, and geochemical tracers (alkalinity) were used to identify the geographic sources of runoff. Input–output relationships of isotopic tracers were explored using a gamma function to fit a transit time distribution, which was used to test the hypothesis that the length of mean transit times increased systematically with the cover of sandstone aquifers in the catchments. However, in three catchments, the increased influence of anthropogenic factors, notably reservoir storage, urban runoff and agricultural abstraction for irrigation, prevented reliable transit time estimation. For sites where tentative mean transit time estimates were possible, these varied from around 1.6 years in upland catchments dominated by metamorphic rocks (>75%) and responsive soils to around 4 years in catchments with 34% sandstone cover and freely draining soils. These preliminary results were consistent with inferences of geochemical tracers on the increased role of sedimentary aquifers as runoff sources in lowland areas, but observation from a larger number of sites is needed to confirm this. Copyright © 2012 John Wiley & Sons, Ltd.     

The relative role of soil type and tree cover on water storage and transmission in northern headwater catchments

Soil water storage and stable isotopes dynamics were investigated in dominant soil–vegetation assemblages of a wet northern headwater catchment (3.2 km²) with limited seasonality in precipitation. We determined the relative influence of soil and vegetation cover on storage and transmission processes. Forested and non‐forested sites were compared, on poorly drained histosols in riparian zones and freely draining podzols on steeper hillslopes. Results showed that soil properties exert a much stronger influence than vegetation on water storage dynamics and fluxes, both at the plot and catchment scale. This is mainly linked to the overall energy‐limited climate, restricting evaporation, in conjunction with high soil water storage capacities. Threshold behaviour in runoff responses at the catchment scale was associated with differences in soil water storage and transmission dynamics of different hydropedological units. Linear input–output relationships occurred when runoff was generated predominantly from the permanently wet riparian histosols, which show only small dynamic storage changes. In contrast, nonlinear runoff generation was related to transient periods of high soil wetness on the hillslopes. During drier conditions, more marked differences in soil water dynamics related to vegetation properties emerged, in terms of evaporation and impacts on temporarily increasing dynamic storage potential. Overall, our results suggest that soil type and their influence on runoff generation are dominant over vegetation effects in wet, northern headwater catchments with low seasonality in precipitation. Potential increase of subsurface storage by tree cover (e.g. for flood management) will therefore be spatially distributed throughout the landscape and limited to rare and extreme dry conditions. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of streamflow isotope sampling strategies on the calibration of a tracer-aided rainfall-runoff model

Isotopes are increasingly used in rainfall-runoff models to constrain conceptualisations of internal catchment functioning and reduce model uncertainty. However, there is little guidance on how much tracer data is required to adequately do this, and different studies use data from different sampling strategies. Here, we used a 7-year time series of daily stable water isotope samples of precipitation and streamflow to derive a range of typical stream sampling regimes and investigate how this impacts calibration of a semi-distributed tracer-aided model in terms of flow, deuterium and flux age simulations. Over the 7 years weekly sampling facilitated an almost identical model performance as daily, and there were only slight deteriorations in performance for fortnightly sampling. Monthly sampling resulted in poorer deuterium simulations and greater uncertainty in the derived parameter sets ability to accurately represent catchment functioning, evidenced by unrealistic reductions in the volumes of water available for mixing in the saturation area causing simulated water age decreases. Reducing sampling effort and restricting data collection to 3 years caused reductions in the accuracy of deuterium simulation, though the deterioration did not occur if sampling continued for 5 years. Analysis was also undertaken to consider the effects of reduced sampling effort over the driest and wettest hydrological years to evaluate effects of more extreme conditions. This showed that the model was particularly sensitive to changes in sampling during dry conditions, when the catchment hydrological response is most non-linear. Across all dataset durations, sampling in relation to flow conditions, rather than time, revealed that samples collected at flows >Q50 could provide calibration results comparable to daily sampling. Targeting only extreme high flows resulted in poor deuterium and low flow simulations. This study suggests sufficient characterization of catchment functioning can be obtained through reduced sampling effort over longer timescales and the targeting of flows >Q50.     

Using stable isotopes to estimate travel times in a data‐sparse Arctic catchment: Challenges and possible solutions

Use of isotopes to quantify the temporal dynamics of the transformation of precipitation into run‐off has revealed fundamental new insights into catchment flow paths and mixing processes that influence biogeochemical transport. However, catchments underlain by permafrost have received little attention in isotope‐based studies, despite their global importance in terms of rapid environmental change. These high‐latitude regions offer limited access for data collection during critical periods (e.g., early phases of snowmelt). Additionally, spatio‐temporal variable freeze–thaw cycles, together with the development of an active layer, have a time variant influence on catchment hydrology. All of these characteristics make the application of traditional transit time estimation approaches challenging. We describe an isotope‐based study undertaken to provide a preliminary assessment of travel times at Siksik Creek in the western Canadian Arctic. We adopted a model–data fusion approach to estimate the volumes and isotopic characteristics of snowpack and meltwater. Using samples collected in the spring/summer, we characterize the isotopic composition of summer rainfall, melt from snow, soil water, and stream water. In addition, soil moisture dynamics and the temporal evolution of the active layer profile were monitored. First approximations of transit times were estimated for soil and streamwater compositions using lumped convolution integral models and temporally variable inputs including snowmelt, ice thaw, and summer rainfall. Comparing transit time estimates using a variety of inputs revealed that transit time was best estimated using all available inflows (i.e., snowmelt, soil ice thaw, and rainfall). Early spring transit times were short, dominated by snowmelt and soil ice thaw and limited catchment storage when soils are predominantly frozen. However, significant and increasing mixing with water in the active layer during the summer resulted in more damped steam water variation and longer mean travel times (~1.5 years). The study has also highlighted key data needs to better constrain travel time estimates in permafrost catchments.     

Fine scale variability of hyporheic hydrochemistry in salmon spawning gravels with contrasting groundwater-surface water interactions

There is increasing realisation of the importance of groundwater–surface water (GW–SW) interactions in understanding freshwater ecology. A study that assessed the influence of local GW–SW interactions on shallow (<250 mm) hyporheic water quality at two contrasting salmon spawning locations in Scotland, UK is reported. At a groundwater-dominated site, continuous logging sensors revealed that hyporheic dissolved oxygen (DO) concentrations changed rapidly in response to changing hydrological conditions. Low volume (25 ml) spot samples revealed fine-scale spatial variability (<0.05 m) consistent with a vertically shifting boundary layer between source waters. At a surface-water-dominated location, hyporheic water was typically characterised by high DO and electrical conductivity values, characteristic of surface water. Small reductions in DO at this site are hypothesised to be associated with short residence hyporheic discharge. A comparison between in-situ (logging DO sensor data) and ex-situ (small volume sampling) methods revealed good agreement, potentially allowing deployment of the two methods in stratified sampling programmes. This study demonstrates that hyporheic water quality varies over fine spatial and temporal scales and that future studies need to design sampling strategies that consider the scales appropriate to both the ecology and the hyporheic processes of interest.     

Isotope hydrology and water sources in a heavily urbanized stream

Complex networks of both natural and engineered flow paths control the hydrology of streams in major cities through spatio-temporal variations in connection and disconnection of diverse water sources. We used spatially extensive and temporally intensive sampling of water stable isotopes to disentangle the hydrological sources of the heavily urbanized Panke catchment (~220 km²) in the north of Berlin, Germany. The isotopic data enabled us to partition stream water sources across the catchment using a Bayesian mixing analysis. The upper part of the catchment streamflow is dominated by groundwater (~75%) from gravel aquifers. In dry summer periods, streamflow becomes intermittent in the upper catchment, possibly as a result of local groundwater abstractions. Storm drainage dominates the responses to precipitation events. Although such events can dramatically change the isotopic composition of the upper stream network, storm drainage only accounts for 10%–15% of annual streamflow. Moving downstream, subtle changes in sources and isotope signatures occur as catchment characteristics vary and the stream is affected by different tributaries. However, effluents from a wastewater treatment plant (WWTP), serving 700,000 people, dominate stream flow in the lower catchment (~90% of annual runoff) where urbanization effects are more dramatic. The associated increase in sealed surfaces downstream also reduces the relative contribution of groundwater to streamflow. The volume and isotopic composition of storm runoff is again dominated by urban drainage, though in the lower catchment, still only about 10% of annual runoff comes from storm drains. The study shows the potential of stable water isotopes as inexpensive tracers in urban catchments that can provide a more integrated understanding of the complex hydrology of major cities. This offers an important evidence base for guiding the plans to develop and re-develop urban catchments to protect, restore, and enhance their ecological and amenity value.     

Inter‐comparison of hydro‐climatic regimes across northern catchments: synchronicity,resistance and resilience

The higher mid‐latitudes of the Northern Hemisphere are particularly sensitive to climate change as small differences in temperature determine frozen ground status, precipitation phase, and the magnitude and timing of snow accumulation and melt. An international inter‐catchment comparison program, North‐Watch, seeks to improve our understanding of the sensitivity of northern catchments to climate change by examining their hydrological and biogeochemical responses. The catchments are located in Sweden (Krycklan), Scotland (Mharcaidh, Girnock and Strontian), the United States (Sleepers River, Hubbard Brook and HJ Andrews) and Canada (Catamaran, Dorset and Wolf Creek). This briefing presents the initial stage of the North‐Watch program, which focuses on how these catchments collect, store and release water and identify ‘types’ of hydro‐climatic catchment response. At most sites, a 10‐year data of daily precipitation, discharge and temperature were compiled and evaporation and storage were calculated. Inter‐annual and seasonal patterns of hydrological processes were assessed via normalized fluxes and standard flow metrics. At the annual‐scale, relations between temperature, precipitation and discharge were compared, highlighting the role of seasonality, wetness and snow/frozen ground. The seasonal pattern and synchronicity of fluxes at the monthly scale provided insight into system memory and the role of storage. We identified types of catchments that rapidly translate precipitation into runoff and others that more readily store water for delayed release. Synchronicity and variance of rainfall–runoff patterns were characterized by the coefficient of variation (cv) of monthly fluxes and correlation coefficients. Principal component analysis (PCA) revealed clustering among like catchments in terms of functioning, largely controlled by two components that (i) reflect temperature and precipitation gradients and the correlation of monthly precipitation and discharge and (ii) the seasonality of precipitation and storage. By advancing the ecological concepts of resistance and resilience for catchment functioning, results provided a conceptual framework for understanding susceptibility to hydrological change across northern catchments. Copyright © 2010 John Wiley & Sons, Ltd.     

Using high‐resolution isotope data and alternative calibration strategies for a tracer‐aided runoff model in a nested catchment

Testing hydrological models over different spatio‐temporal scales is important for both evaluating diagnostics and aiding process understanding. High‐frequency (6‐hr) stable isotope sampling of rainfall and runoff was undertaken during 3‐week periods in summer and winter within 12 months of daily sampling in a 3.2‐km² catchment in the Scottish Highlands. This was used to calibrate and test a tracer‐aided model to assess the (a) information content of high‐resolution data, (b) effect of different calibration strategies on simulations and inferred processes, and (c) model transferability to <1‐km² subcatchment. The 6‐hourly data were successfully incorporated without loss of model performance, improving the temporal resolution of the modelling, and making it more relevant to the time dynamics of the isotope and hydrometric response. However, this added little new information due to old‐water dominance and riparian mixing in this peatland catchment. Time variant results, from differential split sample testing, highlighted the importance of calibrating to a wide range of hydrological conditions. This also provided insights into the nonstationarity of catchment mixing processes, in relation to storage and water ages, which varied markedly depending on the calibration period. Application to the nested subcatchment produced equivalent parameterization and performance, highlighting similarity in dominant processes. The study highlighted the utility of high‐resolution data in combination with tracer‐aided models, applied at multiple spatial scales, as learning tools to enhance process understanding and evaluation of model behaviour across nonstationary conditions. This helps reveal more fully the catchment response in terms of the different mechanistic controls on both wave celerites and particle velocities.