Spatiotemporal variability in hydrochemistry of shallow groundwater in a small pre‐alpine catchment: The importance of landscape elements

Topography and landscape characteristics affect the storage and release of water and, thus, groundwater dynamics and chemistry. Quantification of catchment scale variability in groundwater chemistry and groundwater dynamics may therefore help to delineate different groundwater types and improve our understanding of which parts of the catchment contribute to streamflow. We sampled shallow groundwater from 34 to 47 wells and streamflow at seven locations in a 20‐ha steep mountainous catchment in the Swiss pre‐Alps, during nine baseflow snapshot campaigns. The spatial variability in electrical conductivity, stable water isotopic composition, and major and trace ion concentrations was large and for almost all parameters larger than the temporal variability. Concentrations of copper, zinc, and lead were highest at sites that were relatively dry, whereas concentrations of manganese and iron were highest at sites that had persistent shallow groundwater levels. The major cation and anion concentrations were only weakly correlated to individual topographic or hydrodynamic characteristics. However, we could distinguish four shallow groundwater types based on differences from the catchment average concentrations: riparian zone‐like groundwater, hillslopes and areas with small upslope contributing areas, deeper groundwater, and sites characterized by high magnesium and sulfate concentrations that likely reflect different bedrock material. Baseflow was not an equal mixture of the different groundwater types. For the majority of the campaigns, baseflow chemistry most strongly resembled riparian‐like groundwater for all but one subcatchment. However, the similarity to the hillslope‐type groundwater was larger shortly after snowmelt, reflecting differences in hydrologic connectivity. We expect that similar groundwater types can be found in other catchments with steep hillslopes and wet areas with shallow groundwater levels and recommend sampling of groundwater from all landscape elements to understand groundwater chemistry and groundwater contributions to streamflow.     

The assumption of uniform specific discharge: unsafe at any time?

Nearby catchments in the same landscape are often assumed to have similar specific discharge (runoff per unit catchment area). Five years of streamflow from 14 nested catchments in a 68 km² landscape was used to test this assumption, with the hypothesis that the spatial variability in specific discharge is smaller than the uncertainties in the measurement. The median spatial variability of specific discharge, defined as subcatchment deviation from the catchment outlet, was 33% at the daily scale. This declined to 24% at a monthly scale and 19% at an annual scale. These specific discharge differences are on the same order of magnitude as predicted for major land‐use conversions or a century of climate change. Spatial variability remained when considering uncertainties in specific discharge, and systematic seasonal patterns in specific discharge variation further provide confidence that these differences are more than just errors in the analysis of catchment area, rainfall variability or gauging. Assuming similar specific discharge in nearby catchments can thus lead to spurious conclusions about the effects of disturbance on hydrological and biogeochemical processes. Copyright © 2016 John Wiley & Sons, Ltd.     

Patterns in the linkage of water quantity and quality during low‐flows

This study investigates 72 catchments across the federal state of Baden‐Wuerttemberg, Germany, for changes in water quality during low‐flow events. Data from the state's water quality monitoring network provided seven water quality parameters (water temperature, electrical conductivity, concentrations of chloride, sodium, sulfate, nitrate, and phosphate), which were statistically related to streamflow variability. Water temperature changes during low‐flow showed seasonal dependence. Nitrate concentrations revealed high spatial heterogeneity with about one third of the stations showing decreasing values during low discharge. For most other parameters, concentrations increased during low‐flow. Despite consistent trend directions, the magnitudes of changes with streamflow differed markedly across the state. Both multiple linear regression and a multiple analysis of variances were applied to explain these differences with the help of catchment characteristics. Results indicated that for sulfate and conductivity, geology of the catchments was the most important control, whereas for chloride, sodium, and nitrate, sewage treatment plants had the largest influence. For phosphate, no clear control could be identified. Independent from the applied method, land use was a less important control on river water quality during low‐flow than geology or inflow from sewage treatment plants. These results show that the effects of diffuse and point sources, as well as those of natural and anthropogenic sources differ for different water quality parameters. Overall, a high diversity of potential water quality deterioration signals needs to be considered when the ecological status of rivers is to be protected during low‐flow events.     

Flow dynamics at the continental scale: Streamflow correlation and hydrological similarity

Streamflow variability in space and time critically affects anthropic water uses and ecosystem services. Unfortunately, spatiotemporal patterns of flow regimes are often unknown, as discharge measurements are usually recorded at a limited number of hydrometric stations unevenly distributed along river networks. Advances in understanding the physical processes that control the spatial patterns of river flows are therefore necessary to predict water availability at ungauged locations or to extrapolate pointwise streamflow observations. This work explores the use of the spatial correlation of river flows as a metric to quantify the similarity between hydrological responses of two catchments. Following a stochastic framework, 340,000 cross‐correlations between pairs of daily streamflows time series are predicted at a seasonal timescale across the contiguous United States using 413 catchments of the MOPEX dataset. Model predictions of streamflow correlation obtained in absence of run‐off information are successfully used to identify catchment outlets sharing similar discharge dynamics and flow regimes across a broad range of geomorphoclimatic conditions, without relying on calibration. The selection of reference streamgauges based on predicted streamflow correlation generally outperforms the selection based on spatial proximity, especially as the density of available gauged sections decreases. Interestingly, correlated outlets share a broad spectrum of hydrological signatures (mean discharge, flow variability, and recession properties), suggesting that catchments forced by analogous frequency and intensity of effective rainfall events might exhibit common geomorphoecological traits leading to similar hydrological responses. The proposed framework provides a physical basis to assist the regionalization of flow dynamics and to interpret the spatial variability of flow regimes along stream networks.     

Quantifying aggregation and change in runoff source in accordance with catchment area increase in a forested headwater catchment

There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5‐km² forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End‐member mixing analysis was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km². The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km², becoming almost constant above 0.1 km². The specific discharge of groundwater showed large variability below 1 km² and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km², it was necessary to consider variations in the discharge rates of both subsurface water and groundwater. Copyright © 2016 John Wiley & Sons, Ltd.     

Contributing sources to baseflow in pre‐alpine headwaters using spatial snapshot sampling

Mountainous headwaters consist of different landscape units including forests, meadows and wetlands. In these headwaters it is unclear which landscape units contribute what percentage to baseflow. In this study, we analysed spatiotemporal differences in baseflow isotope and hydrochemistry to identify catchment‐scale runoff contribution. Three baseflow snapshot sampling campaigns were performed in the Swiss pre‐alpine headwater catchment of the Zwäckentobel (4.25 km²) and six of its adjacent subcatchments. The spatial and temporal variability of δ²H, Ca, DOC, AT, pH, SO₄, Mg and H₄SiO₄ of streamflow, groundwater and spring water samples was analysed and related to catchment area and wetland percentage using bivariate and multivariate methods. Our study found that in the six subcatchments, with variable arrangements of landscape units, the inter‐ and intra catchment variability of isotopic and hydrochemical compositions was small and generally not significant. Stream samples were distinctly different from shallow groundwater. An upper spring zone located near the water divide above 1,400 m and a larger wetland were identified by their distinct spatial isotopic and hydrochemical composition. The upstream wetland percentage was not correlated to the hydrochemical streamflow composition, suggesting that wetlands were less connected and act as passive features with a negligible contribution to baseflow runoff. The isotopic and hydrochemical composition of baseflow changed slightly from the upper spring zone towards the subcatchment outlets and corresponded to the signature of deep groundwater. Our results confirm the need and benefits of spatially distributed snapshot sampling to derive process understanding of heterogeneous headwaters during baseflow. Copyright © 2015 John Wiley & Sons, Ltd.     

Effect of bedrock flow on catchment rainfall‐runoff characteristics and the water balance in forested catchments in Tanzawa Mountains,Japan

In this paper, we examined the role of bedrock groundwater discharge and recharge on the water balance and runoff characteristics in forested headwater catchments. Using rigorous observations of catchment precipitation, discharge and streamwater chemistry, we quantified net bedrock flow rates and contributions to streamwater runoff and the water balance in three forested catchments (second‐order to third‐order catchments) underlain by uniform bedrock in Japan. We found that annual rainfall in 2010 was 3130 mm. In the same period, annual discharge in the three catchments varied from 1800 to 3900 mm/year. Annual net bedrock flow rates estimated by the chloride mass balance method at each catchment ranged from ?1600 to 700 mm/year. The net bedrock flow rates were substantially different in the second‐order and third‐order catchments. During baseflow, discharge from the three catchments was significantly different; conversely, peak flows during large storm events and direct runoff ratios were not significantly different. These results suggest that differences in baseflow discharge rates, which are affected by bedrock flow and intercatchment groundwater transfer, result in the differences in water balance among the catchments. This study also suggests that in these second‐order to third‐order catchments, the drainage area during baseflow varies because of differences between the bedrock drainage area and surface drainage area, but that the effective drainage area during storm flow approaches the surface drainage area. Copyright © 2012 John Wiley & Sons, Ltd.     

A comparison of streamflow,salt and water balances in adjacent farmland and forest catchments in south‐western Victoria,Australia

A study of the hydrologic effects of catchment change from pasture to plantation was carried out in Gatum, south‐western Victoria, Australia. This study describes the hydrologic characteristics of two adjacent catchments: one with 97% grassland and the other one with 62% Eucalyptus globulus plantations. Streamflow from both catchments was intermittent during the 20‐month study period. Monthly streamflow was always greater in the pasture‐dominated catchment compared with the plantation catchment because of lower evapotranspiration in the pasture‐based catchment. This difference in streamflow was also observed even during summer 2010/2011 when precipitation was 74% above average (1954–2012) summer rainfall. Streamflow peaks in the plantation‐based catchment were smaller than in the pasture‐dominated system. Flow duration curves show differences between the pasture and plantation‐dominated catchments and affect both high‐flow and low‐flow periods. Groundwater levels fell (up to 4.4 m) in the plantation catchment during the study period but rose (up to 3.2 m) in the pasture catchment. Higher evapotranspiration in the plantation catchment resulted in falling groundwater levels and greater disconnection of the groundwater system from the stream, resulting in lower baseflow contribution to streamflow. Salt export from each catchment increases with increasing flow and is higher at the pasture catchment, mainly because of the higher flow. Reduced salt loading to streams due to tree planting is generally considered environmentally beneficial in saline areas of south‐eastern Australia, but this benefit is offset by reduced total streamflow. Copyright © 2014 John Wiley & Sons, Ltd.     

Can climate variability information constrain a hydrological model for an ungauged Costa Rican catchment?

Long‐term hydrological data are key to understanding catchment behaviour and for decision making within water management and planning. Given the lack of observed data in many regions worldwide, such as Central America, hydrological models are an alternative for reproducing historical streamflow series. Additional types of information—to locally observed discharge—can be used to constrain model parameter uncertainty for ungauged catchments. Given the strong influence that climatic large‐scale processes exert on streamflow variability in the Central American region, we explored the use of climate variability knowledge as process constraints to constrain the simulated discharge uncertainty for a Costa Rican catchment, assumed to be ungauged. To reduce model uncertainty, we first rejected parameter relationships that disagreed with our understanding of the system. Then, based on this reduced parameter space, we applied the climate‐based process constraints at long‐term, inter‐annual, and intra‐annual timescales. In the first step, we reduced the initial number of parameters by 52%, and then, we further reduced the number of parameters by 3% with the climate constraints. Finally, we compared the climate‐based constraints with a constraint based on global maps of low‐flow statistics. This latter constraint proved to be more restrictive than those based on climate variability (further reducing the number of parameters by 66% compared with 3%). Even so, the climate‐based constraints rejected inconsistent model simulations that were not rejected by the low‐flow statistics constraint. When taken all together, the constraints produced constrained simulation uncertainty bands, and the median simulated discharge followed the observed time series to a similar level as an optimized model. All the constraints were found useful in constraining model uncertainty for an—assumed to be—ungauged basin. This shows that our method is promising for modelling long‐term flow data for ungauged catchments on the Pacific side of Central America and that similar methods can be developed for ungauged basins in other regions where climate variability exerts a strong control on streamflow variability.     

Understanding hydrological processes in glacierized catchments: Evidence and implications of highly variable isotopic and electrical conductivity data

The spatial and temporal characterization of geochemical tracers over Alpine glacierized catchments is particularly difficult, but fundamental to quantify groundwater, glacier melt, and rain water contribution to stream runoff. In this study, we analysed the spatial and temporal variability of δ²H and electrical conductivity (EC) in various water sources during three ablation seasons in an 8.4‐km² glacierized catchment in the Italian Alps, in relation to snow cover and hydro‐meteorological conditions. Variations in the daily streamflow range due to melt‐induced runoff events were controlled by maximum daily air temperature and snow covered area in the catchment. Maximum daily streamflow decreased with increasing snow cover, and a threshold relation was found between maximum daily temperature and daily streamflow range. During melt‐induced runoff events, stream water EC decreased due to the contribution of glacier melt water to stream runoff. In this catchment, EC could be used to distinguish the contribution of subglacial flow (identified as an end member, enriched in EC) from glacier melt water to stream runoff, whereas spring water in the study area could not be considered as an end member. The isotopic composition of snow, glacier ice, and melt water was not significantly correlated with the sampling point elevation, and the spatial variability was more likely affected by postdepositional processes. The high spatial and temporal variability in the tracer signature of the end members (subglacial flow, rain water, glacier melt water, and residual winter snow), together with small daily variability in stream water δ²H dynamics, are problematic for the quantification of the contribution of the identified end members to stream runoff, and call for further research, possibly integrated with other natural or artificial tracers.     

Hydrologic response to valley‐scale structure in alpine headwaters

Few systematic studies of valley‐scale geomorphic drivers of streamflow regimes in complex alpine headwaters have compared response between catchments. As a result, little guidance is available for regional‐scale hydrological research and monitoring efforts that include assessments of ecosystem function. Physical parameters such as slope, elevation range, drainage area and bedrock geology are often used to stratify differences in streamflow response between sampling sites within an ecoregion. However, these metrics do not take into account geomorphic controls on streamflow specific to glaciated mountain headwaters. The coarse‐grained nature of depositional features in alpine catchments suggests that these landforms have little water storage capacity because hillslope runoff moves rapidly just beneath the rock mantle before emerging in fluvial networks. However, recent studies show that a range of depositional features, including talus slopes, protalus ramparts and ‘rock‐ice’ features may have more storage capacity than previously thought. To better evaluate potential differences in streamflow response among basins with extensive coarse depositional features and those without, we examined the relationships between streamflow discharge, stable isotopes, water temperature and the amplitude of the diurnal signal at five basin outlets. We also quantified the percentages of colluvial channel length measured along the stepped longitudinal profile. Colluvial channels, characterized by the presence of surficial, coarse‐grained depositional features, presented sediment‐rich, transport‐limited morphologies that appeared to have a cumulative effect on the timing and volume of flow downstream. Measurements taken from colluvial channels flowing through depositional landforms showed median recession constants (K_r) of 0.9–0.95, δ¹⁸O values of ≥?14.5 and summer diurnal amplitudes ≤0.8 as compared with more typical surface water recession constant values of 0.7, δ¹⁸O ≤ ?13.5 and diurnal amplitudes >2.0. Our results demonstrated strong associations between the percentage of colluvial channel length within a catchment and moderated streamflow regimes, water temperatures, diurnal signals and depleted δ¹⁸O related to groundwater influx. Copyright © 2014 John Wiley & Sons, Ltd.     

Projecting climate change impacts on stream flow regimes with tracer‐aided runoff models ‐ preliminary assessment of heterogeneity at the mesoscale

The northern mid‐high latitudes form a region that is sensitive to climate change, and many areas already have seen – or are projected to see – marked changes in hydroclimatic drivers on catchment hydrological function. In this paper, we use tracer‐aided conceptual runoff models to investigate such impacts in a mesoscale (749 km²) catchment in northern Scotland. The catchment encompasses both sub‐arctic montane sub‐catchments with high precipitation and significant snow influence and drier, warmer lowland sub‐catchments. We used downscaled HadCM3 General Circulation Model outputs through the UKCP09 stochastic weather generator to project the future climate. This was based on synthetic precipitation and temperature time series generated from three climate change scenarios under low, medium and high greenhouse gas emissions. Within an uncertainty framework, we examined the impact of climate change at the monthly, seasonal and annual scales and projected impacts on flow regimes in upland and lowland sub‐catchments using hydrological models with appropriate process conceptualization for each landscape unit. The results reveal landscape‐specific sensitivity to climate change. In the uplands, higher temperatures result in diminishing snow influence which increases winter flows, with a concomitant decline in spring flows as melt reduces. In the lowlands, increases in air temperatures and re‐distribution of precipitation towards autumn and winter lead to strongly reduced summer flows despite increasing annual precipitation. The integration at the catchment outlet moderates these seasonal extremes expected in the headwaters. This highlights the intimate connection between hydrological dynamics and catchment characteristics which reflect landscape evolution. It also indicates that spatial variability of changes in climatic forcing combined with differential landscape sensitivity in large heterogeneous catchments can lead to higher resilience of the integrated runoff response. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of forestry on summertime low flows and physical fish habitat in snowmelt‐dominant headwater catchments of the Pacific Northwest

Periods of summertime low flows are often critical for fish. This study quantified the impacts of forest clear‐cutting on summertime low flows and fish habitat and how they evolved through time in two snowmelt‐dominant headwater catchments in the southern interior of British Columbia, Canada. A paired‐catchment analysis was applied to July–September water yield, the number of days each year with flow less than 10% of mean annual discharge, and daily streamflow for each calendar day. The postharvest time series were divided into treatment periods of approximately 6–10 years, which were analysed independently to evaluate how the effects of forestry changed through time. An instream flow assessment using a physical habitat simulation‐style approach was used to relate streamflow to the availability of physical habitat for resident rainbow trout. About two decades after the onset of logging and as the extent of logging increased to approximately 50% of the catchments, reductions in daily summertime low flows became more significant for the July–September yield (43%) and for the analysis by calendar day (11–68%). Reductions in summertime low flows were most pronounced in the catchment with the longest postharvest time series. On the basis of the temporal patterns of response, we hypothesize that the delayed reductions in late‐summer flow represent the combined effects of a persistent advance in snowmelt timing in combination with at least a partial recovery of transpiration and interception loss from the regenerating forests. These results indicate that asymptotic hydrological recovery as time progresses following logging is not suitable for understanding the impacts of forest harvesting on summertime low flows. Additionally, these reductions in streamflow corresponded to persistent decreases in modelled fish habitat availability that typically ranged from 20% to 50% during the summer low‐flow period in one of the catchments, suggesting that forest harvest may have substantial delayed effects on rearing salmonids in headwater streams.     

Hydrological resilience to forest fire in the subarctic Canadian shield

Understanding the role of forest fires on water budgets of subarctic Precambrian Shield catchments is important because of growing evidence that fire activity is increasing. Most research has focused on assessing impacts on individual landscape units, so it is unclear how changes manifest at the catchment scale enough to alter water budgets. The objective of this study was to determine the water budget impact of a forest fire that partially burned a ~450 km² subarctic Precambrian Shield basin. Water budget components were measured in a pair of catchments: one burnt and another unburnt. Burnt and unburnt areas had comparable net radiation, but thaw was deeper in burned areas. There were deeper snow packs in burns. Differences in streamflow between the catchments were within measurement uncertainty. Enhanced winter streamflow from the burned watershed was evident by icing growth at the streamflow gauge location, which was not observed in the unburned catchment. Wintertime water chemistry was also clearly elevated in dissolved organics, and organic-associated nutrients. Application of a framework to assess hydrological resilience of watersheds to wildfire reveal that watersheds with both high bedrock and open water fractions are more resilient to hydrological change after fire in the subarctic shield, and resilience decreases with increasingly climatically wet conditions. This suggests significant changes in runoff magnitude, timing and water chemistry of many Shield catchments following wildfire depend on pre-fire land cover distribution, the extent of the wildfire and climatic conditions that follow the fire.     

Basin-scale analysis of rainfall and runoff in Peru (1969–2004): Pacific,Titicaca and Amazonas drainages

Abstract

We present a procedure for estimating Q95 low flows in both gauged and ungauged catchments where Q95 is the flow that is exceeded 95% of the time. For each step of the estimation procedure, a number of alternative methods was tested on the Austrian data set by leave-one-out cross-validation, and the method that performed best was used in the final procedure. To maximise the accuracy of the estimates, we combined relevant sources of information including long streamflow records, short streamflow records, and catchment characteristics, according to data availability. Rather than deriving a single low flow estimate for each catchment, we estimated lower and upper confidence limits to allow local information to be incorporated in a practical application of the procedure. The components of the procedure consist of temporal (climate) adjustments for short records; grouping catchments into eight seasonality-based regions; regional regressions of low flows with catchment characteristics; spatial adjustments for exploiting local streamflow data; and uncertainty assessment. The results are maps of lower and upper confidence limits of low flow discharges for 21 000 sub-catchments in Austria.     

Predicting Aquifer Response Time for Application in Catchment Modeling

It is well established that changes in catchment land use can lead to significant impacts on water resources. Where land‐use changes increase evapotranspiration there is a resultant decrease in groundwater recharge, which in turn decreases groundwater discharge to streams. The response time of changes in groundwater discharge to a change in recharge is a key aspect of predicting impacts of land‐use change on catchment water yield. Predicting these impacts across the large catchments relevant to water resource planning can require the estimation of groundwater response times from hundreds of aquifers. At this scale, detailed site‐specific measured data are often absent, and available spatial data are limited. While numerical models can be applied, there is little advantage if there are no detailed data to parameterize them. Simple analytical methods are useful in this situation, as they allow the variability in groundwater response to be incorporated into catchment hydrological models, with minimal modeling overhead. This paper describes an analytical model which has been developed to capture some of the features of real, sloping aquifer systems. The derived groundwater response timescale can be used to parameterize a groundwater discharge function, allowing groundwater response to be predicted in relation to different broad catchment characteristics at a level of complexity which matches the available data. The results from the analytical model are compared to published field data and numerical model results, and provide an approach with broad application to inform water resource planning in other large, data‐scarce catchments.     

The dynamics of natural pipe hydrological behaviour in blanket peat

Natural soil pipes are found in peatlands, but little is known about their hydrological role. This paper presents the most complete set of pipe discharge data to date from a deep blanket peatland in Northern England. In a 17.4‐ha catchment, we identified 24 perennially flowing and 60 ephemerally flowing pipe outlets. Eight pipe outlets along with the catchment outlet were continuously gauged over an 18‐month period. The pipes in the catchment were estimated to produce around 13.7% of annual streamflow, with individual pipes often producing large peak flows (maximum peak of 3.8 l s^?1). Almost all pipes, whether ephemerally or perennially flowing, shallow or deep (outlets > 1 m below the peat surface), showed increased discharge within a mean of 3 h after rainfall commencement and were dominated by stormflow, indicating good connectivity between the peatland surface and the pipes. However, almost all pipes had a longer period between the hydrograph peak and the return to base flow compared with the stream (mean of 23.9 h for pipes, 19.7 h for stream). As a result, the proportion of streamflow produced by the pipes at any given time increased at low flows and formed the most important component of stream discharge for the lowest 10% of flows. Thus, a small number of perennially flowing pipes became more important to the stream system under low‐flow conditions and probably received water via matrix flow during periods between storms. Given the importance of pipes to streamflow in blanket peatlands, further research is required into their wider role in influencing stream water chemistry, water temperature and fluvial carbon fluxes, as well as their role in altering local hydrochemical cycling within the peat mass itself. Enhanced piping within peatlands caused by environmental change may lead to changes in the streamflow regime with larger low flows and more prolonged drainage of the peat. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessment of the impact of pools on the water isotopic signature of a boreal patterned peatland

Recent studies have shown that boreal peatlands exhibit considerable chemical variability but without clear spatial pattern. This chemical heterogeneity illustrates the complex hydrological behaviour of peatlands, particularly patterned fen. Isotopic, chemical and physical tracers were used to describe the hydrological behaviour of a small boreal headwater catchment (13 ha) during the snow‐free period with a special emphasis on the downstream patterned fen. Results showed that shallow pools were mixed every day during the summer, particularly during nights or discharge periods. Despite large water storage capacities in pools, which should induce large buffer effect, hydrological behaviour of patterned fen is more similar to a piston flow process. This is probably because of the division of the fen into successive small cascading streamflow reservoirs. The consequences were a rapid change of the chemical signature throughout the fen, particularly upstream. A spatial pattern was observed downstream in early summer. The isotopic signature passed from an upstream depleted and homogeneous signature to a progressively enriched downstream signature. However, this pattern was not identified during the wetter period (late summer), probably because the discharge, which dominated the water budget, decreased the surface water residence time and flushed a large proportion of stored surface water. We developed for this patterned fen a conceptual model of the surface flow to explain these particular mixing effects and the implications on the dynamics of the chemical signature. To further our understanding of similar boreal headwater catchments, future work should include the development of a multiple mixed‐reservoir model. Copyright © 2015 John Wiley & Sons, Ltd.