Seasonality of epCO2 at different scales along an integrated river continuum within the Dee basin,NE Scotland

Spatial and temporal variations of free CO₂ concentrations in surface waters are mainly controlled by dynamic processes encompassing terrestrial inputs and in‐stream biotic cycling. Free CO₂ can be expressed as ‘excess partial pressure of CO₂’ (epCO₂), indicating supersaturation or under‐saturation with respect to atmospheric CO₂. Seasonal patterns of epCO₂ at sites draining nested upland catchments between 3·40 and 1837 km² were assessed within the River Dee basin in NE Scotland. EpCO₂ values ranged from 0·14 at the lowermost site on the mainstem in autumn to 12·7 on a major tributary during spring. A seasonality index was derived to describe contrasting winter/spring maxima and summer/autumn minima as annual mean epCO₂ values could not clearly distinguish between different sites. Seasonal differences tended to increase downstream as progressive changes in physicochemical conditions enhanced the influence of autotrophic communities. However, perturbations to this continuum occurred as CO₂ inputs from high DOC, heterotrophic tributaries and land‐use changes between open moorland and forest affected downstream continuity. Major tributaries showed reduced differences between seasons compared to the mainstem. Smaller headwaters indicated a lack of seasonality as high connectivity of responsive, peaty soils enabled continual inputs of terrestrially derived free CO₂ to streams concomitant with limited autotrophic CO₂ removal, maintaining epCO₂ > 1 throughout. Seasonality of epCO₂ was mainly driven by the ability of in‐stream biota to consume CO₂ during optimal conditions in summer/autumn. This was confirmed by multiple linear regression analysis, which indicated that, apart from catchment area, baseflow and biotic activity indicators were the best predictors of epCO₂ seasonality characteristics at any particular stage of the river system. Copyright © 2009 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.     

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.     

Using lumped conceptual rainfall-runoff models to simulate daily isotope variability with fractionation in a nested mesoscale catchment

This paper presents 19 months of stable isotope (δ²H and δ¹⁸O) data to enhance understanding of water and solute transport at two spatial scales (2.3 km² and 122 km²) in the agricultural Lunan catchment, Scotland. Daily precipitation and stream isotope data, weekly lake and spring isotope data and monthly groundwater isotope data revealed important insights into flow pathways and mixing of water at both scales. In particular, a deeper groundwater flow path significantly contributes to total streamflow (25-50%). Upstream lake isotope dynamics, susceptible to evaporative fractionation, also appeared to have an important influence on the downstream isotope composition. This unique tracer data set facilitated the conceptualization of a lumped catchment-scale flow-tracer model. The incorporation of hydrological, mixing and fractionation processes based on these data improved simulations of the stream δ²H isotope response at the catchment outlet from 0.37 to 0.56 for the Nash-Sutcliffe statistic. The stable isotope data successfully aided model conceptualization and calibration in the quest for a simple water and solute transport model with improved representation of process dynamics.     

Analysis of hydrological seasonality across northern catchments using monthly precipitation–runoff polygon metrics

Abstract

Seasonality is an important hydrological signature for catchment comparison. Here, the relevance of monthly precipitation–runoff polygons (defined as scatter points of 12 monthly average precipitation–runoff value pairs connected in the chronological monthly sequence) for characterizing seasonality patterns was investigated to describe the hydrological behaviour of 10 catchments spanning a climatic gradient across the northern temperate region. Specifically, the research objectives were to: (a) discuss the extent to which monthly precipitation–runoff polygons can be used to infer active hydrological processes in contrasting catchments; (b) test the ability of quantitative metrics describing the shape, orientation and surface area of monthly precipitation–runoff polygons to discriminate between different seasonality patterns; and (c) examine the value of precipitation–runoff polygons as a basis for catchment grouping and comparison. This study showed that some polygon metrics were as effective as monthly average runoff coefficients for illustrating differences between the 10 catchments. The use of precipitation–runoff polygons was especially helpful to look at the dynamics prevailing in specific months and better assess the coupling between precipitation and runoff and their relative degree of seasonality. This polygon methodology, linked with a range of quantitative metrics, could therefore provide a new simple tool for understanding and comparing seasonality among catchments.

Editor Z.W. Kundzewicz; Associate editor K. Heal

Citation Ali, G., Tetzlaff, D., Kruitbos, L., Soulsby, C., Carey, S., McDonnell, J., Buttle, J., Laudon, H., Seibert, J., McGuire, K., and Shanley, J., 2013. Analysis of hydrological seasonality across northern catchments using monthly precipitation–runoff polygon metrics. Hydrological Sciences Journal, 59 (1), 56–72.     

Storage,mixing, and fluxes of water in the critical zone across northern environments inferred by stable isotopes of soil water

Quantifying soil water storage, mixing, and release via recharge, transpiration, and evaporation is essential for a better understanding of critical zone processes. Here, we integrate stable isotope (²H and ¹⁸O of soil water, precipitation, and groundwater) and hydrometric (soil moisture) data from 5 long‐term experimental catchments along a hydroclimatic gradient across northern latitudes: Dry Creek (USA), Bruntland Burn (Scotland), Dorset (Canada), Krycklan (Sweden), and Wolf Creek (Canada). Within each catchment, 6 to 11 isotope sampling campaigns occurred at 2 to 4 sampling locations over at least 1 year. Analysis for ²H and ¹⁸O in the bulk pore water was done for >2,500 soil samples either by cryogenic extraction (Dry Creek) or by direct equilibration (other sites). The results showed a similar general pattern that soil water isotope variability reflected the seasonality of the precipitation input signal. However, pronounced differences among sampling locations occurred regarding the isotopic fractionation due to evaporation. We found that antecedent precipitation volumes mainly governed the fractionation signal, temperature and evaporation rates were of secondary importance, and soil moisture played only a minor role in the variability of soil water evaporation fractionation across the hydroclimatic gradient. We further observed that soil waters beneath conifer trees were more fractionated than beneath heather shrubs or red oak trees, indicating higher soil evaporation rates in coniferous forests. Sampling locations closer to streams were more damped and depleted in their stable isotopic composition than hillslope sites, revealing increased subsurface mixing towards the saturated zone and a preferential recharge of winter precipitation. Bulk soil waters generally comprised a high share of waters older than 14 days, which indicates that the water in soil pores are usually not fully replaced by recent infiltration events. The presented stable isotope data of soil water were, thus, a useful tool to track the spatial variability of water fluxes within and from the critical zone. Such data provide invaluable information to improve the representation of critical zone processes in spatially distributed hydrological models.     

Topographic,pedologic and climatic interactions influencing streamflow generation at multiple catchment scales

Dominant flow pathways (DFPs) in mesoscale watersheds are poorly characterized and understood. Here, we make use of a conservative tracer (Gran alkalinity) and detailed information about climatic conditions and physical properties to examine how temporally and spatially variable factors interact to determine DFPs in 12 catchments draining areas from 3.4 to 1829.5 km² (Cairngorms, Scotland). After end‐member mixing was applied to discriminate between near surface and deep groundwater flow pathways, variation partitioning, canonical redundancy analyses and regression models were used to resolve: (i) What is the temporal variability of DFPs in each catchment?; (ii) How do DFPs change across spatial scales and what factors control the differences in hydrological responses?; and (iii) Can a conceptual model be developed to explain the spatiotemporal variability of DFPs as a function of climatic, topographic and soil characteristics? Overall, catchment characteristics were only useful to explain the temporal variability of DFPs but not their spatial variation across scale. The temporal variability of DFPs was influenced most by prevailing hydroclimatic conditions and secondarily soil drainability. The predictability of active DFPs was better in catchments with soils supporting fast runoff generation on the basis of factors such as the cumulative precipitation from the seven previous days, mean daily air temperature and the fractional area covered by Rankers. The best regression model R² was 0.54, thus suggesting that the catchments’ internal complexity was not fully captured by the factors included in the analysis. Nevertheless, this study highlights the utility of combining tracer studies with digital landscape analysis and multivariate statistical techniques to gain insights into the temporal (climatic) and spatial (topographic and pedologic) controls on DFPs. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling phosphorus inputs from agricultural sources and urban areas in river basins

An area-differentiated model approach (MEPhos) for the quantification of mean annual P-inputs from point and diffuse sources is presented. The following pathways are considered: artificial drainage, wash-off, groundwater outflow, soil erosion, rainwater sewers, combined sewer overflows, municipal waste water treatment plants and industrial effluents. Two retention functions for rivers and reservoirs are included in order to model P-sinks within a river basin. This allows a complete record of P-loads in heterogeneous meso- and macroscale river basins and enables validation of modeling results with water quality data on a load basis. The model is applied to the River Ruhr basin (4,485 km²) in Germany, which includes contrasting natural conditions, land use patterns as well as population and industry densities. Based on validated modelling results sub-areas of high P-loads are localized and management options for the reduction of P-inputs to surface waters are proposed taking into account the site conditions of the sub-areas relevant for high P-inputs into surface waters.     

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.     

Calculated wind climatology of the South-Saxonian/North-Czech mountain topography including improved resolution of mountains

A mesoscale model has been applied to calculate climatological means of the surface wind. A reliable average requires more than 40 model runs, which are differentiated by the direction and speed of the geostrophic wind under the assumption of neutral stratification. The frequency distributions of the geostrophic wind have been taken from observations of the 850-hPa winds at the radiosonde station in Prague for a 10-year period. The simulation results have been averaged over all sectors and speed classes of the geostrophic wind according to their frequencies. A comparison of the calculated mean wind speeds with observed ones shows deviations of about 0.4 ms⁻¹ outside the mountains. The representation of steep topography and isolated mountains on the basis of a 3-km horizontal resolution of the simulations needs special treatment in order to reduce the gap of up to 4 ms⁻¹ between observed and simulated mean wind speeds over mountains. Therefore, an empiric speed-up formula has been applied to the isolated mountains that otherwise would fall through the 3-km meshes. The corresponding deviations have been reduced to 1.5 ms⁻¹.