Search

International audience; Large proportions of rain water and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time varying processes. CFC‐based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km2 forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from 9 different wells with depths of 2 m to 18 m close to the stream network. Immediately below the water table, CFC‐based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age‐depth relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity.

Flash floods are the highest sediment transporting agent, but are inaccessible for in-situ sampling and have rarely been analyzed by remote sensing technology. Laboratory and field experiments were done to develop linear spectral unmixing (LSU) remote sensing model and evaluate its performance in simulating the suspended sediment concentration (SSC) in flash floods. The models were developed from continuous monitoring in the laboratory and the onsite spectral signature of river bed sediment deposits and flash floods in the Tekeze River and in its tributary, the Tsirare River. The Pearson correlation coefficient was used to determine the variability of correlations between reflectance and SSCs. The coefficient of determination (R2) and root mean square of error (RMSE) were used to evaluate the performance of the generated models. The results found that the Pearson correlation coefficient between SSCs and reflectance varied based on the level of the SSCs, geological colors, and grain sizes. The performance of the LSU model and empirical remote sensing approaches were computed to be R2 = 0.92, and RMSE = ±0.76 g/l in the Tsirare River and R2 = 0.91, and RMSE = ±0.73 g/l in the Tekeze River and R2 = 0.81, RMSE = ±2.65 g/l in the Tsirare river and R2 = 0.76, RMSE = ±10.87 g/l in the Tekeze River, respectively. Hence, the LSU approach of remote sensing was found to be relatively accurate in monitoring and modeling the variability of SSCs that could be applied to the upper Tekeze River basin.

International audience; The knowledge of the sources of base cations in stream water is a prerequisite to assess potential effects of changing environmental conditions such as changing rainfall, weathering or groundwater flows on cation export with stream water. This study use stable Mg isotopes to identify potential sources in the well-studied catchment of Krycklan located on gneissic bedrock covered by quaternary sediments in Sweden. Samples were collected from open filed rain, throughfall, stream, soil, rock and litterfall. The δ26Mg values of these samples was determined and the contributions of different sources to Mg fluxes in the stream were determined from the variation of the Mg isotope and Sr / Mg ratios.The results show an overall variation of 1.10‰ between all samples. In addition, Magnesium isotope ratios varied little in the streamwater and in soil solution, except during snowmelt periods during which a large portion of the annual runoff occurs. Magnesium in the streamwater is explained as a mixture of three pools (open field rain, soil solution and groundwater) with the latter two influenced by catchment processes. Outside the snowmelt period, Mg in streamwater mainly derived from the groundwater, assumed to be mineral weathering signature in this catchment, with a contribution ranging from 12 to 63% to Mg fluxes. Open field rain dominates Mg fluxes in streamwater during spring flood (0 to 78%) and may contribute significantly during larger summer and autumn rainfall events. Soil solution input to streamwater range from 16 to 59% of Mg fluxes in streamwater. Our results demonstrate that δ26Mg values together with Mg concentrations and Sr/Mg ratios can be used to constrain the Mg sources of stream water and quantify weathering release rates.

[Departement_IRSTEA]Eaux [ADD1_IRSTEA]Systèmes aquatiques soumis à des pressions multiples; International audience; Human water use, climate change and land conversion have created a water crisis for billions of individuals and many ecosystems worldwide. Global water stocks and fluxes are estimated empirically and with computer models, but this information is conveyed to policymakers and researchers through water cycle diagrams. Here we compiled a synthesis of the global water cycle, which we compared with 464 water cycle diagrams from around the world. Although human freshwater appropriation now equals half of global river discharge, only 15% of the water cycle diagrams depicted human interaction with water. Only 2% of the diagrams showed climate change or water pollution—two of the central causes of the global water crisis—which effectively conveys a false sense of water security. A single catchment was depicted in 95% of the diagrams, which precludes the representation of teleconnections such as ocean–land interactions and continental moisture recycling. These inaccuracies correspond with specific dimensions of water mismanagement, which suggest that flaws in water diagrams reflect and reinforce the misunderstanding of global hydrology by policymakers, researchers and the public. Correct depictions of the water cycle will not solve the global water crisis, but reconceiving this symbol is an important step towards equitable water governance, sustainable development and planetary thinking in the Anthropocene.

Analysis of hydrological impacts of land use change raises questions about whether, and how much, such impacts are misrepresented because of errors in river flow observations. In this paper, land u ...

For most catchments, there is insufficient data to determine the location of the groundwater surface. For humid climates, it is, therefore, often assumed that the groundwater-surface follows the surface topography. This assumption allows using digital elevation models (DEMs) to estimate the flow directions and catchment boundaries. However, high-resolution elevation data also include many small-scale features that are unlikely to affect the direction of groundwater flow, or only affect it during specific conditions. Furthermore, flow directions may change during events or depending on the water level.The optimal resolution of the DEM for determining groundwater flow directions is not known yet. Therefore, we studied how much DEM derived flow directions and catchment boundaries are affected by the resolution or smoothing of the elevation data for the Krycklan catchment in northern Sweden. We also measured the groundwater levels in two small sub-catchments to determine what DEM resolution best describes the actual groundwater-surface and flow directions.For the topographic analyses, the LiDAR-based elevation data were first smoothed with various filters (e.g., Gaussian filters) and resampled to obtain lower resolution elevation data. We then determined the flow directions for these different DEMs. The aim was to determine where in the catchment the calculated flow directions are most sensitive to the resolution of the topographic data. The results of the topographic analyses show that for some areas, particularly flat areas, ridges, streambanks and locations where the local slope differs from the general slope, the calculated flow directions depend strongly on the resolution and smoothing of the elevation data.To test how well the DEM based groundwater flow directions represent actual flow directions, we installed a dense (5-20 m spacing) network of shallow (1 to 6 m deep) groundwater wells (75 wells in total) in a 1 ha and a 2 ha gauged sub-catchment. The triangular nested design of the groundwater well network allowed us to determine the smaller (5 m) and larger scale (20 m) groundwater gradients. The recorded water levels were augmented and validated by manual measurements during the summers of 2018 and 2019. The high spatial and temporal resolution data allowed us to study the response of the groundwater level and the flow directions to different meteorological situations (e.g., large precipitation events after dry and wet conditions and during a very dry period in summer 2018). These observations indicate that the degree to which the groundwater-surface is a subdued copy of the surface topography varies throughout the year, and provides information on which DEM resolution most accurately represents the groundwater-surface and flow directions.

The deep porewater of northern peatlands stores large amounts of carbon dioxide (CO2). This store is viewed as a stable feature in the peatland CO2 cycle. Here, we report large and rapid fluctuations in deep porewater CO2 concentration recurring every autumn over four consecutive years in a boreal peatland. Estimates of the vertical diffusion of heat indicate that CO2 diffusion occurs at the turbulent rather than molecular rate. The weakening of porewater thermal stratification in autumn likely increases turbulent diffusion, thus fostering a rapid diffusion of deeper porewater CO2 towards the surface where net losses occur. This phenomenon periodically decreases the peat porewater CO2 store by between 29 and 90 g C m−2 throughout autumn, which is comparable to the peatland’s annual C-sink. Our results establish the need to consider the role of turbulent diffusion in regularly destabilizing the CO2 store in peat porewater., The CO2 stored in deep peat porewater is viewed as a fixed component of peatland C cycling. Here, the authors reveal a hitherto unknown hydro-physical process that results in sudden losses from this CO2 store every autumn.

It has been 12 years since Bishop et al. (2008) wrote the Invited Commentary “Aqua Incognita: the unknown headwaters.” They highlighted that “In most regions, the overwhelming majority of stream length lies beyond the frontiers of any systematic documentation and would have to be represented as a blank space on the assessment map. This means that for the majority of streams that support aquatic life, a systematic understanding is lacking on water quality, habitat, biota, specific discharge, or even how many kilometres of such streams are there. This blank space is so vast that it deserves a name to help us at least to remember that it is there. We propose calling it ‘Aqua Incognita’” (Bishop et al., 2008, p. 1239). We continue to agree with this statement and the need to understand headwater streams better. In this commentary, we want to draw attention to a particular type of headwater stream that is even less frequently examined: headwater streams that flow intermittently, that is, the Aqua Temporaria Incognita. Question 3 of the 23 unsolved problems in hydrology (Blöschl et al., 2019) focuses on ephemeral dryland streams. We argue that this focus needs broadening to headwater temporary streams because they are ubiquitous in all climates. Headwater temporary streams feed larger perennial streams and are particularly sensitive to climate change and other human influences (Jaeger, Olden, & Pelland, 2014; Pumo, Caracciolo, Viola, & Noto, 2016; Reynolds, Shafroth, & LeRoy Poff, 2015). Their effective management and protection, therefore, requires an understanding of both natural and artificial causes of intermittence., We thank the SMIRES (Science and Management of Intermittent Rivers & Ephemeral Streams) COST Action CA15113 for facilitating our many discussions on temporary streams, Willy Bertin and Gaela Le Bechec from SR3A and Judy England from the Environment Agency for information on how and why they use the CrowdWater app, and Helena Ramos Ribeiro and Auria Buchs for collecting the data shown in Figure 2.

The distribution of groundwater ages with depth provides information about subsurface structures and flow dynamics. Upslope measured groundwater age stratifications are commonly used to estimate groundwater recharge rates, whereas downslope measured age stratifications are influenced by recharge conditions, the aquifer structure and interactions between groundwater and surface water. In our study we address the question whether downslope measured groundwater ages from different locations can provide spatial and temporal information about catchment-scale groundwater dynamics and the relationship between groundwater recharge and discharge.We derived an overall groundwater age stratification, representative for the Svartberget subcatchment (0.47 km2) located within the Krycklan study site, by measuring CFCs from 9 different sampling locations with depths of 2 m to 18 m. All sampling locations were downslope and located in basal till which is overlain by ablation till. The CFC-based groundwater age stratification reveals an unexpected pattern, with groundwater ages of already 30 years immediately below the water table. Groundwater ages increase then with depth. We could reproduce the observed groundwater age stratification by using a groundwater flow model and show that the lag of rejuvenation, noticeable in groundwater ages of 30 years at the water table, derives from return flow of groundwater at a subsurface discharge zone that evolves at the interface between the two soil types (basal and ablation till). Furthermore, we demonstrate by varying the infiltration rate how the extent of the discharge zone and the partitioning of the infiltration amount to the two layers change, i.e. young runoff in the upper layer (ablation till) and old groundwater circulation through the deeper layer (basal till).By providing a simple analytical approximations of the observed groundwater age stratification, we show that the extent of the subsurface discharge zone is a powerful indicator of the relation between the recharge and discharge zone, while the vertical gradient of the age-depth relationship provides information about the overall aquifer recharge.

Peat stored in large wetlands plays an important role in the carbon cycle and strongly influences water quality of terrestrial surface water bodies. At the same time, peat is also stored in the direct vicinity of many boreal forest streams. From this strategic position, peat can receive and chemically reset hillslope water before it reaches the stream network. Yet, in contrast to large wetlands, only little spatial information is available on the lateral extent of near-stream peat and even less about its vertical variation. Here, we present field data on peat depth and lateral extent collected from approximately 200 transects (with 12 soil profiles taken per transect) distributed across the entire stream network of the Krycklan boreal catchment in Northern Sweden. This soil profile data revealed a considerable heterogeneity of peat and organic horizon thicknesses. By combining the field data with morphological and geological maps, we show how parent material, stream order and local topography influence near stream peat structures. Furthermore, we discuss potential consequences on surface water quality by linking the detailed peat data set to estimates of lateral, shallow groundwater inflows derived from hydrometric measurements and digital terrain analysis.

A better understanding of the dominant source areas and transport pathways of pesticide losses to surface water is needed for targeting mitigation efforts in a more cost-effective way. To this end, we monitored pesticides in surface water in an agricultural catchment typical of one of the main crop production regions in Sweden. Three small sub-catchments (88-242ha) were selected for water sampling based on a high-resolution digital soil map developed from proximal sensing methods and soil sampling; one sub-catchment had a high proportion of clay soils, another was dominated by coarse sandy soils while the third comprised a mix of soil types. Samples were collected from the stream, from field drains discharging into the stream and from within-field surface runoff during spring and early summer in three consecutive years. These samples were analyzed by LC-MS/MS for 99 compounds, including most of the polar and semi-polar pesticides frequently used in Swedish agriculture. Information on pesticide applications (products, doses and timing) was obtained from annual interviews with the farmers. There were clear and consistent differences in pesticide occurrence in the stream between the three sub-catchments, with both the numbers of detected compounds and concentrations being the largest in the area with a high proportion of clay soils and with very few detections in the sandy sub-catchment. Macropore flow to drains was most likely the dominant loss pathway in the studied area. Many of the compounds that were detected in drainage and stream water samples had not been applied for several years. This suggests that despite the predominant role of fast flow pathways in determining losses to the stream, long-term storage along the transport pathways also occurs, presumably in subsoil horizons where degradation is slow.

Basin infiltration managed aquifer recharge (MAR) is a commonly used method for storing and treating surface water to be used as drinking water. This study examined how the removal of bacteriophage MS2 was affected by the relative age of the sand used for basin infiltration MAR at 4 °C using batch experiments (static and agitated) and column experiments. The sand and the water used in all experiments were characteristic of that used at a basin infiltration MAR scheme in Uppsala, Sweden. Experimental data was fit with numerical models describing the fate and transport of virus in soil—water systems. The “used” sand that had been subjected to intermittent infiltration over a period of eight years had fifteen times the amount of organic carbon than the “new” sand which had yet to be used for infiltration. Results showed that attachment of MS2 to the new sand in batch experiments was relatively irreversible. Attachment to the used sand was reversible. Inactivation of MS2 was slowed when it was attached to the used sand. Results for the column experiments showed that the removal rate of MS2 was significantly lower in columns of used sand than in columns of new sand. Simulations indicated that MS2 would be entirely removed in the infiltration basins with new sand. Less than three log-removal was estimated for basins with used sand. Reduced removal of MS2 by the used sand was deemed to be most likely due to organic coatings on the used sand. Results of this study give deeper insight into the mechanisms responsible for removing virus in infiltration basins and how those mechanisms will change as the sand in the basin is exposed to infiltration water.

Excess nitrogen fertiliser in agricultural soils might be leached to streams and converted to the greenhouse gas nitrous oxide (N2O). To assess the importance of N2O emissions from agricultural streams, concentration dynamics and emissions N2O emissions in streams were investigated in a 32 km2 lowland agricultural catchment located in Sweden. Dissolved N2O concentration was measured at nine occasions between December 2014 and August 2015 at nine stream stations. The stream stations represented sub-catchments with different land use characteristics with agricultural land use ranging from 0 to 63% of the area. Stream N2O percentage saturation ranged 40–2701% and showed large spatial and temporal variations. Statistical analysis using mixed models revealed that N2O concentration was significantly linked to nitrate concentration in the stream water, to the percentage arable land in the sub-catchments as well as to the stream water discharge. Using two empirical equations to estimate the N2O emissions showed that streams were generally a source of N2O to the atmosphere (mean 108 and 175 μg N m−2 h−1 with first and second equation). The catchment scale estimate of N2O stream emissions was compared to the estimate obtained using IPCC guidelines linking N fertilisation inputs and leaching to N2O emissions. The comparison suggested that N2O stream emission calculated using the IPCC methodology might be underestimated. A coarse estimate suggests that N2O stream emissions represent about 4% of the total N2O emissions from N-fertiliser at the catchment scale. Hence while streams covered only 0.1% of the catchment area they were of disproportionate importance as a source of N2O to the atmosphere.

Summary Defining the catchment transit time distribution remains a challenge. Here, we used a new semi-analytical physically-based integrated subsurface flow and advective–dispersive particle movement model to assess the subsurface controls on subsurface water flow paths and transit time distributions. First, we tested the efficacy of the new model for simulation of the observed groundwater dynamics at the well-studied S-transect hillslope (Vastrabacken sub-catchment, Sweden). This system, like many others, is characterized by exponential decline in saturated hydraulic conductivity and porosity with soil depth. The model performed well relative to a tracer-based estimate of transit time distribution as well as observed groundwater depth–discharge relationship within 30 m of the stream. Second, we used the model to assess the effect of changes in the subsurface permeability architecture on flow pathlines and transit time distribution in a set of virtual experiments. Vertical patterns of saturated hydraulic conductivity and porosity with soil depth significantly influenced hillslope transit time distribution. Increasing infiltration rates significantly decreased mean groundwater age, but not the distribution of transit times relative to mean groundwater age. The location of hillslope hydrologic boundaries, including the groundwater divide and no-flow boundary underlying the hillslope, changed the transit time distribution less markedly. These results can guide future decisions on the degree of complexity that is warranted in a physically-based rainfall–runoff model to efficiently and explicitly estimate time invariant subsurface pathlines and transit time distribution.

Understanding how water moves through catchments - from the time it enters as precipitation to when it exits via streamflow - is of fundamental importance to understanding hydrological and biogeoch ...

Soil compaction is a common consequence of forestry traffic traversing unprotected, moist soils; it decreases porosity and affects hydraulic conductivity even in coarse-textured soils. The aim here was to study root-zone hydrology and vegetation in three microsites (in, between, and beside wheel tracks) 4 to 5 yr after forwarder traffic, on stony and sandy till soils in two clearcuts in northern Sweden. Measurements of soil volumetric water content (VWC), vegetation indicators and one-dimensional hydrological modeling (Hydrus-1D) of wheel tracks and undisturbed soil were conducted. Soil VWC was monitored hourly during 2017 and 2018 in three or four plots along a slope on each site. Soil VWC was also measured once with a portable sensor in 117 plots along two slopes at each site, where the vegetation was recorded and analyzed using Ellenberg indicator indexes. Soil VWC was highest in wheel tracks and lowest between tracks; this was corroborated by the species composition in the wheel tracks (Ellenberg indicator for soil moisture). Bare soil was more frequent in wheel tracks and between tracks than in undisturbed soil. The model simulations indicated that the changed soil hydraulic properties influenced the VWC results in the wheel tracks. However, the differences in average pressure heads in the root zone were small between the microsites and only apparent during dry periods. In the wheel tracks, air-filled porosity was

Stream CO2 emissions contribute significantly to atmospheric climate forcing. While there are strong indications that groundwater inputs sustain these emissions, the specific biogeochemical pathways and timescales involved in this lateral CO2 export are still obscure. Here, via an extensive radiocarbon (14C) characterisation of CO2 and DOC in stream water and its groundwater sources in an old-growth boreal forest, we demonstrate that the 14C-CO2 is consistently in tune with the current atmospheric 14C-CO2 level and shows little association with the 14C-DOC in the same waters. Our findings thus indicate that stream CO2 emissions act as a shortcut that returns CO2 recently fixed by the forest vegetation to the atmosphere. Our results expose a positive feedback mechanism within the C budget of forested catchments, where stream CO2 emissions will be highly sensitive to changes in forest C allocation patterns associated with climate and land-use changes., There is a growing consensus that groundwater inflow supplies most of the C load to streams, but the sources and timescales generating this flux are still unknown. Here, the authors demonstrate that soil respiration, derived from current forest carbon fixation, fuels stream CO2 emissions.

The study evaluates the biogeochemical model ForSAFE-2D, designed to simulate water and chemical transport from the forest to the stream, by simulating the hydrology and the transport of the chemical tracer chloride (Cl−) along a forest hillslope in Northern Sweden. The simulated Cl− exports were in balance with the simulated inputs but measurements suggested a net release of Cl− from the catchment. Underestimated deposition inputs (deposition peaks and possibly dry deposition) were probably and partially responsible for this mismatch. However, we could not exclude that other soil biogeochemical processes omitted in ForSAFE-2D could also contribute to Cl− exports from the catchment. The study showed that ForSAFE-2D is a promising tool to better understand the factors that regulate the chemical export from the forest to the stream. The results also confirmed that there are limitations in using Cl− as a tracer in forest ecosystems.

Studies on hydrology, biogeochemistry, or mineral weathering often rely on assumptions about flow paths, water storage dynamics, and transit times. Testing these assumptions requires detailed hydro ...

Improving the understanding of how stream flow dynamics are influenced by landscape characteristics, such as soils, vegetation and terrain, is a central endeavor of catchment hydrology. Here we inv ...

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 68km(2) landsc ...

Dissolved organic carbon (DOC) is a significant constituent in aquatic ecosystems with concentrations in streams influenced by both temperature and water flow pathway dynamics associated with changes in discharge (streamflow). We investigated the sensitivity of DOC concentrations in 12 high-latitude headwater streams to changes in temperature and discharge using a mathematical model. The implications of differences in sensitivities were explored by using downscaled projections of air temperature and discharge to simulate possible trajectories of DOC concentrations in a changing climate. We found two distinct responses: (i) catchments where stream DOC sensitivity was high to temperature but low to discharge and (ii) catchments where stream DOC sensitivity was low to temperature but high to discharge. Streams with strong seasonal DOC dynamics were more sensitive to temperature changes than nonseasonal systems. In addition, stream DOC sensitivity to discharge was strongly correlated with vertical soil water DOC differences in the near-stream zone. Simulations of possible future changes in DOC concentrations indicated median increases of about 4-24% compared to current levels when using projections of air temperature and discharge but even larger increases were observed for base flow concentrations (13-42%). Streams with high-temperature sensitivity showed the largest increases in DOC concentrations. Our results suggest that future climatic changes could bring significant increases in surface water DOC concentrations in boreal and hemiboreal areas but that the response ultimately is dependent on vertical soil solution DOC differences and soil organic carbon distribution.

Using a 15 year stream record from a northern boreal catchment, we demonstrate that the inter-annual variation in dissolved organic carbon (DOC) concentrations during snowmelt was related to discharge, winter climate and previous DOC export. A short and intense snowmelt gave higher stream water DOC concentrations, as did long winters, while a high previous DOC export during the antecedent summer and autumn resulted in lower concentrations during the following spring. By removing the effect of discharge we could detect that the length of winter affected the modeled soil water DOC concentrations during the following snowmelt period, which in turn affected the concentrations in the stream. Winter climate explained more of the stream water DOC variations than previous DOC export during the antecedent summer and autumn.

Quantitative estimates of the land-atmosphere exchange of gaseous elemental mercury (GEM) are biased by the measurement technique employed, because no standard method or scale in space and time are agreed upon. Here we present concurrent GEM exchange measurements over a boreal peatland using a novel relaxed eddy accumulation (REA) system, a rectangular Teflon (R) dynamic flux chamber (DFC) and a DFC designed according to aerodynamic considerations (Aero-DFC). During four consecutive days the DFCs were placed alternately on two measurement plots in every cardinal direction around the REA sampling mast. Spatial heterogeneity in peat surface characteristics (0-34 cm) was identified by measuring total mercury in eight peat cores (57 +/- 8 ng g(-1), average SE), vascular plant coverage (32-52%), water table level (4.5-14.1 cm) and dissolved gaseous elemental mercury concentrations (28-51 pg L-1) in the peat water. The GEM fluxes measured by the DFCs showed a distinct diel pattern, but no spatial difference in the average fluxes was detected (ANOVA, alpha = 0.05). Even though the correlation between the Teflon DFC and Aero-DFC was significant (r = 0.76, p < 0.05) the cumulative flux of the Aero-DFC was a factor of three larger. The average flux of the Aero-DFC (1.9 ng m(-2) h(-1)) and REA (2 ng m(-2) h(-1)) were in good agreement. The results indicate that the novel REA design is in agreement for cumulative flux estimates with the Aero-DFC, which incorporates the effect of atmospheric turbulence. The comparison was performed over a fetch with spatially rather homogenous GEM flux dynamics under fairly consistent weather conditions, minimizing the effect of weather influence on the data from the three measurement systems. However, in complex biomes with heterogeneous surface characteristics where there can be large spatial variability in GEM gas exchange, the small footprint of chambers ( < 0.2 m(2)) makes for large coefficients of variation. Thus many chamber measurement replications are needed to establish a credible biome GEM flux estimate, even for a single point in time. Dynamic flux chambers will, however, be able to resolve systematic differences between small scale features, such as experimentally manipulated plots or small scale spatial heterogeneity.

Water’s journey from rain to stream by Harald Grip and Allan Rodhe (1985, in Swedish: Vattnets vag fran regn till back) was one of the first textbooks to present groundwater contributions as a major feature of runoff generation, with implications for water quality and management. Three decades later, we have the privilege of presenting a special issue of Hydrological Processes, ‘Runoff Generation in a Nordic Light: 30 Years with Water’s Journey from Rain to Stream’ that seeks to introduce the book to a larger audience and continue the journey of ideas that the authors set in motion with their book.

Headwater streams export CO2 as lateral downstream export and vertical evasion from the stream surface. CO2 in boreal headwater streams generally originates from adjacent terrestrial areas, so determining the sources and rate of CO2 transport along the hillslope–riparian–stream continuum could improve estimates of CO2 export via the aquatic pathway, especially by quantifying evasion at higher temporal resolutions. Continuous measurements of dissolved CO2 concentrations and water table were made along the hillslope–riparian–stream continuum in the Västrabäcken sub-catchment of the Krycklan catchment, Sweden. Daily water and CO2 export from the hillslope and riparian zone were estimated over one hydrological year (October 2012–September 2013) using a flow-concentration model and compared with measured lateral downstream CO2 export. Total water export over the hydrological year from the hillslope was 230 mm yr−1 compared with 270 mm yr−1 from the riparian zone. This corresponds well (proportional to the relative upslope contributing area) to the annual catchment runoff of 265 mm yr−1. Total CO2 export from the riparian zone to the stream was 3.0 g CO2-C m−2 yr−1. A hotspot for riparian CO2 export was observed at 30–50 cm depth (accounting for 71 % of total riparian export). Seasonal variability was high with export peaks during the spring flood and autumn storm events. Downstream lateral CO2 export (determined from stream water dissolved CO2 concentrations and discharge) was 1.2 g CO2-C m−2 yr−1. Subtracting downstream lateral export from riparian export (3.0 g CO2-C m−2 yr−1) gives 1.8 g CO2-C m−2 yr−1 which can be attributed to evasion losses (accounting for 60 % of export via the aquatic pathway). The results highlight the importance of terrestrial CO2 export, especially from the riparian zone, for determining catchment aquatic CO2 losses and the importance of the CO2 evasion component to carbon export via the aquatic conduit.

The current generation of catchment travel time distribution (TTD) research, integrating nearly three decades of work since publication of Water's Journey from Rain to Stream, seeks to represent the full distribution in catchment travel times and its temporal variability. Here, we compare conceptualizations of increasing complexity with regards to mixing of water storages and evaluate how these assumptions influence time-variable TTD estimates for two catchments with contrasting climates: the Gardsjon catchment in Sweden and the Marshall Gulch catchment in Arizona, USA. Our results highlight that, as long as catchment TTDs cannot be measured directly but need to be inferred from input-output signals of catchments, the inferred catchment TTDs depend strongly on the underlying assumptions of mixing within a catchment. Furthermore, we found that the conceptualization of the evapotranspiration flux strongly influences the inferred travel times of stream discharge. For the wet and forested Gardsjon catchment in Sweden, we inferred that evapotranspiration most likely resembles a completely mixed sample of the water stored in the catchment; however, for the drier Marshall Gulch catchment in Arizona, evapotranspiration predominantly contained the younger water stored in the catchment. For the Marshall Gulch catchment, this higher probability for young water in evapotranspiration resulted in older water in the stream compared to travel times inferred with assumptions of complete mixing. New observations that focus on the TTD of the evapotranspiration flux and the actual travel time of water through a catchment are necessary to improve identification of mixing and consequently travel times of stream water. Copyright © 2014 John Wiley & Sons, Ltd.

Forestry operations can increase the export of mercury (both total and methyl) to surface waters. However, little is known about the relative contribution of different forestry practices. We address this question using a paired-catchment study that distinguishes the effects of site preparation from the antecedent logging. Runoff water from three catchments, two harvested and one untreated control, was sampled biweekly during one year prior to logging, two years after logging, and three years after site preparation. The logging alone did not significantly increase the concentrations of either total or methyl-mercury in runoff, but export increased by 50-70% in one of the harvested catchments as a consequence of increased runoff volume. The combined effects of logging and site preparation increased total and methyl-mercury concentrations by 30-50% relative to preharvest conditions in both treated catchments. The more pronounced concentration effect after site preparation compared to logging could be related to site preparation being conducted during summer. This caused more soil disturbance than logging, which was done during winter with snow covering the ground. The results suggest that the cumulative impact of forest harvest on catchment mercury outputs depends on when and how forestry operations are implemented.

Boreal headwater streams have been identified as hot spots for evasion of greenhouse gases (GHGs). This study was the first to systematically determine the concentrations of CO2 and CH4 in hemiboreal headwater streams. The use of a headspace sampling method focusing on GHGs in combination with a statistically representative selection of more than 200 streams across two regions in Sweden was the basis for defining the base flow concentrations of CO2 and CH4. All streams were supersaturated relative to the atmosphere in CO2 and the majority in CH4 for the 82% of streams in which CH4 was detected. The spatial variability in both CO2 and CH4 was high but positively related to total organic carbon, mean annual temperature, and proportion of peatland in the catchment. There were, however, regional differences in the spatial controls, which are something that predictive models need to consider. The data set allowed for comparison between a headspace and an alkalinity-based method for determining CO2. More than 50% of the streams contained no alkalinity which made the alkalinity-based determination of CO2 impossible. In addition, half of the streams with alkalinity had alkalinities low enough (

Community perception of forest–water relationship was gathered using participatory rural appraisal tools in four watersheds of the Blue Nile River Basin in Ethiopia. These were compared and contrasted with the observational records of forest cover and water flow. Upstream and downstream communities were assessed separately to check for differences in perception based on location within a watershed. The key result of the study was that people in the study watersheds had a range of perceptions about the forest–water relationship which were watershed specific. The perceptions were generally consistent with observational evidence from the same watersheds. This study highlighted the need for locale-specific approaches to land and water management in the Basin, as well as the potential value of using community perceptions to complement the observational records which can have spatial and temporal limitations.

Forest disturbance such as clear-cutting has been identified as an important factor for increasing dissolved organic carbon (DOC) concentrations in boreal streams. We used a long-term data set of s ...

There is still a need for catchment hydrological and transport models that properly integrate the effects of preferential flows while accounting for differences in velocities and celerities. A modeling methodology is presented here which uses particle tracking methods to simulate both flow and transport in multiple pathways in a single consistent solution. Water fluxes and storages are determined by the volume and density of particles and transport is attained by labeling the particles with information that may be tracked throughout the lifetime of that particle in the catchment. The methodology allows representation of preferential flows through the use of particle velocity distributions, and mixing between pathways can be achieved with pathway transition probabilities. A transferable 3-D modeling methodology is presented for the first time and applied to a unique step-shift isotope experiment that was carried out at the 0.63 ha G1 catchment in Gardsjon, Sweden. This application highlights the importance of combining flow and transport in hydrological representations, and the importance of pathway velocity distributions and interactions in obtaining a satisfactory representation of the observations.

Forty-five years (1960–2004) of hydrological data from 12 watersheds in the Abbay Basin, Ethiopia, were tested for possible trends over the entire time series and differences in medians (step-wise changes) between three sub-periods. The classification of the sub-periods was based on the major political changes in 1975 and 1991. Variables investigated were rainfall (P), total flow (Q t), high flow (Q h), low flow (Q l), low flow index (LFI) and run-off coefficient (C). Data were checked for outliers, errors and homogeneity. Trend was tested after serial and cross-correlation tests. The data for each variable were serially uncorrelated from 1 to 10 lag years. There were five globally significant trends out of 50 test cases and 36 significant step-wise changes out of 180 tests. The majority of the significant changes were watershed specific. Run-off coefficient was the single variable showing a consistently increasing trend and stood for ca. 25 % of the total significant trends and step-wise changes. Half of these changes occurred after 1991. We concluded that despite the land use policy changes in 1975 and 1991, as well as the long-term soil degradation, the hydrological regime was quite stable over the 45-year period, with the exception of an increase in the run-off coefficient in the latter part of the run-off record in some watersheds.

Riparian zones (RZ) are a major factor controlling water chemistry in forest streams. Base cations' (BC) concentrations, fluxes, and cycling in the RZ merit attention because a changing climate and increased forest harvesting could have negative consequences, including re-acidification, for boreal surface waters. We present a two-year study of BC and silica (Si) flow-weighted concentrations from 13 RZ and 14 streams in different landscape elements of a boreal catchment in northern Sweden. The spatial variation in BC and Si dynamics in both RZ and streams was explained by differences in landscape element type, with highest concentrations in silty sediments and lowest concentrations in peat-dominated wetland areas. Temporal stability in BC and Si concentrations in riparian soil water, remarkably stable Mg/Ca ratios, and homogeneous mineralogy suggest that patterns found in the RZ are a result of a distinct mineralogical upslope signal in groundwater. Stream water Mg/Ca ratios indicate that the signal is subsequently maintained in the streams. Flow-weighted concentrations of Ca, Mg, and Na in headwater streams were represented by the corresponding concentrations in the RZ, which were estimated using the Riparian Flow-Concentration Integration Model (RIM) approach. Stream and RZ flow-weighted concentrations differed for K and Si, suggesting a stronger biogeochemical influence on these elements, including K recirculation by vegetation and retention of Si within the RZ. Potential increases in groundwater levels linked to forest harvesting or changes in precipitation regimes would tend to reduce BC concentrations from RZ to streams, potentially leading to episodic acidification.

[1] Land cover changes can have significant impacts on hydrological regime. The objective of this study was to detect possible hydrological changes of four watersheds in the Blue Nile Basin using a model-based method for hydrological change detection. The four watersheds, Birr, Upper-Didesa, Gilgel Abbay, and Koga range in size from 260 to 1800 km2. The changes were assessed based on model parameters, model residuals, and in the overall function of the watersheds in transferring rainfall into runoff. The entire time series (1960–2004) was divided into three periods based on political and land management policy changes. A conceptual rainfall-runoff model, the HBV (Hydrologiska Byrans Vattenbalansavdelning) model, was used for the analysis, and suitable parameter sets for each period were found based on a Monte Carlo approach. The values of six out of nine parameters changed significantly between the periods. Model residuals also showed significant changes between the three periods in three of the four watersheds. On the other hand, the overall functioning of the watersheds in processing rainfall to runoff changed little. So even though the individual parameters and model residuals were changing, the integrated functioning of the watersheds showed minimal changes. This study demonstrated the value of using different approaches for detecting hydrological change and highlighted the sensitivity of the outcome to the applied modeling and statistical methods.

[1] Aquatic CO2 concentrations are highly variable and strongly linked to discharge, but until recently, measurements have been largely restricted to low-frequency manual sampling. Using new in situ CO2 sensors, we present concurrent, high-frequency (

We investigated the effects of forest clear-cutting on snow accumulation, melt dynamics and associated stream responses in a paired catchment experiment in northern Sweden. Two boreal first-order s ...

International audience; Predicting hydrological catchment behaviour based on measurable (and preferably widely available) catchment characteristics has been one of the main goals of hydrological modelling. Residence time distributions provide synoptic information about catchment functioning and can be useful metrics to predict their behaviours. Moreover, residence time distributions highlight a wide range of characteristic scales (spatial and temporal) and mixing processes. However, catchment-specific heterogeneity means that the link between residence time distributions and catchment characteristics is complex. Investigating this link for a wide range of catchments could reveal the role of topography, geology, land-use, climate and other factors in controlling catchment hydrology. Meaningful comparison is often challenging given the diversity of data and model structures and formats. To address this need, we are introducing a new virtual platform called Catchment virtual Observatory for Sharing flow and transport models outputs (COnSOrT). The goal of COnSOrT is to promote catchment intercomparison by sharing calibrated model outputs. Compiling commensurable results in COnSOrT will help evaluate model performance, quantify inter-catchment controls on hydrology, and identify research gaps and priorities in catchment science. Researchers interested in sharing or using calibrated model results are invited to participate in the virtual observatory. Participants may test post-processing methods on a wide range of catchment environments to evaluate the generality of their findings.

Groundwater flowing from hillslopes through riparian (near-stream) soils often undergoes chemical transformations that can substantially influence stream water chemistry. We used landscape analysis to predict total organic carbon (TOC) concentration profiles and groundwater levels measured in the riparian zone (RZ) of a 67 km2 catchment in Sweden. TOC exported laterally from 13 riparian soil profiles was then estimated based on the riparian flow–concentration integration model (RIM). Much of the observed spatial variability of riparian TOC concentrations in this system could be predicted from groundwater levels and the topographic wetness index (TWI). Organic riparian peat soils in forested areas emerged as hotspots exporting large amounts of TOC. These TOC fluxes were subject to considerable temporal variations caused by a combination of variable flow conditions and changing soil water TOC concentrations. Mineral riparian gley soils, on the other hand, were related to rather small TOC export rates and were characterized by relatively time-invariant TOC concentration profiles. Organic and mineral soils in RZs constitute a heterogeneous landscape mosaic that potentially controls much of the spatial variability of stream water TOC. We developed an empirical regression model based on the TWI to move beyond the plot scale and to predict spatially variable riparian TOC concentration profiles for RZs underlain by glacial till.

Several major articles from the past decade and beyond conclude the impact of reforestation or afforestation on water yield is negative: additional forest cover will reduce and removing forests will raise downstream water availability. A second group of authors argue the opposite: planting additional forests should raise downstream water availability and intensify the hydrologic cycle. Obtaining supporting evidence for this second group of authors has been more difficult due to the larger scales at which the positive effects of forests on the water cycle may be seen. We argue that forest cover is inextricably linked to precipitation. Forest-driven evapotranspiration removed from a particular catchment contributes to the availability of atmospheric moisture vapor and its cross-continental transport, raising the likelihood of precipitation events and increasing water yield, in particular in continental interiors more distant from oceans. Seasonal relationships heighten the importance of this phenomenon. We review the arguments from different scales and perspectives. This clarifies the generally beneficial relationship between forest cover and the intensity of the hydrologic cycle. While evidence supports both sides of the argument – trees can reduce runoff at the small catchment scale – at larger scales, trees are more clearly linked to increased precipitation and water availability. Progressive deforestation, land conversion from forest to agriculture and urbanization have potentially negative consequences for global precipitation, prompting us to think of forest ecosystems as global public goods. Policy-making attempts to measure product water footprints, estimate the value of ecosystem services, promote afforestation, develop drought mitigation strategies and otherwise manage land use must consider the linkage of forests to the supply of precipitation.

This paper explores the flow paths and turnover times within a catchment characterized by the transmissivity feedback mechanism where there is a strong increase in the saturated hydraulic conductiv ...

Although water storage is an important variable to understand the hydrological functioning of a catchment, it is challenging to estimate the total water storage in a catchment. Catchment water stor ...

We bring together three decades of research from a boreal catchment to facilitate an improved mechanistic understanding of surface water dissolved organic carbon (DOC) regulation across multiple scales. The Krycklan Catchment Study encompasses 15 monitored nested research catchments, ranging from 3 to 6900 ha in size, as well as a set of monitored transects of forested and wetland soils. We show that in small homogenous catchments, hydrological functioning provides a first order control on the temporal variability of stream water DOC. In larger, more heterogeneous catchments, stream water DOC dynamics are regulated by the combined effect of hydrological mechanisms and the proportion of major landscape elements, such as wetland and forested areas. As a consequence, streams with heterogeneous catchments undergo a temporal switch in the DOC source. In a typical boreal catchment covered by 10-20% wetlands, DOC originates predominantly from wetland sources during low flow conditions. During high flow, the major source of DOC is from forested areas of the catchment. We demonstrate that by connecting knowledge about DOC sources in the landscape with detailed hydrological process understanding, an improved representation of stream water DOC regulation can be provided. The purpose of this study is to serve as a framework for appreciating the role of regulating mechanisms, connectivity and scaling for understanding the pattern and dynamics of surface water DOC across complex landscapes. The results from this study suggest that the sensitivity of stream water DOC in the boreal landscape ultimately depends on changes within individual landscape elements, the proportion and connectivity of these affected landscape elements, and how these changes are propagated downstream.

Short-term increases in soil solution nitrate (NO3-) concentration are often observed after forest harvest, even in N-limited systems. We model NO3- leaching below the rooting zone as a function of site productivity. Using national forest inventories and published estimates of N attenuation in rivers and the riparian zone, we estimate effects of stem-only harvesting on NO3- leaching to groundwater, surface waters and the marine environment. Stem-only harvesting is a minor contributor to NO3- pollution of Swedish waters. Effects in surface waters are rapidly diluted downstream, but can be locally important for shallow well-waters as well as for the total amount of N reaching the sea. Harvesting adds approximately 8 Gg NO3-N to soil waters in Sweden, with local concentrations up to 7 mg NO3-N 1(-1). Of that, similar to 3.3 Gg reaches the marine environment. This is similar to 3% of the overall Swedish N load to the Baltic. (C) 2010 Elsevier Ltd. All rights reserved.

Small streams with catchment areas

To determine if forestry affects stream water dissolved organic carbon (DOC) concentrations, we conducted high frequency water sampling at a clear-cut catchment experiment in northern Sweden 1 year after harvesting. The overall finding was that harvesting significantly increased stream water DOC in these boreal forest catchments, at least during the growing season. The results indicate a DOC concentration increase of up to 50% during early summer on the two harvested catchments relative to the two control catchments. The analysis supports the hypothesis that a raised groundwater level following harvesting caused the increased DOC concentration during both hydrological episodes and low flow conditions. Harvesting resulted in a 70% increase in DOC export due to the combined effect of runoff and DOC concentration during the June–October study period. Given the extent of forestry activity in the boreal landscape, these results demonstrate that tree harvesting will affect the water quality of the region.

This paper defines landscape-scale patterns in the character of natural organic matter (NOM) and tests for relationships to catchment soil, vegetation and topography. The drainage network of a boreal catchment, subcatchment size 0.12–78 km2, in Northern Sweden was sampled in August 2002 during a period of stable low water flow. The NOM was characterized with UV/Vis spectroscopy, fluorescence, XAD-8 fractionation (%humic substances), gel permeation chromatography (apparent molecular weight), and elemental composition (C:N). The largest spatial variation was found for C:N, absorbance ratio, and specific visible absorptivity. The lowest variation was in fluorescence index, %humic substances and molecular retention time. But the variation in total organic carbon (TOC), iron and aluminium concentration was more than twice that of C:N. Between headwater and downstream sites no significant changes were distinguished in the NOM character. At stream reaches, junctions and lakes little change (

Summary Two different but complementary modelling approaches for reproducing the observed dynamics of total organic carbon (TOC) in a boreal stream are presented. One is based on a regression analysis, while the other is based on riparian soil conditions using a convolution of flow and concentration. Both approaches are relatively simple to establish and help to identify gaps in the process understanding of the TOC transport from soils to catchments runoff. The largest part of the temporal variation of stream TOC concentrations (4–46 mg L−1) in a forested headwater stream in the boreal zone in northern Sweden may be described using a four-parameter regression equation that has runoff and transformed air temperature as sole input variables. Runoff is assumed to be a proxy for soil wetness conditions and changing flow pathways which in turn caused most of the stream TOC variation. Temperature explained a significant part of the observed inter-annual variability. Long-term riparian hydrochemistry in soil solutions within 4 m of the stream also captures a surprisingly large part of the observed variation of stream TOC and highlights the importance of riparian soils. The riparian zone was used to reproduce stream TOC with the help of a convolution model based on flow and average riparian chemistry as input variables. There is a significant effect of wetting of the riparian soil that translates into a memory effect for subsequent episodes and thus contributes to controlling stream TOC concentrations. Situations with high flow introduce a large amount of variability into stream water TOC that may be related to memory effects, rapid groundwater fluctuations and other processes not identified so far. Two different climate scenarios for the region based on the IPCC scenarios were applied to the regression equation to test what effect the expected increase in precipitation and temperature and resulting changes in runoff would have on stream TOC concentrations assuming that the soil conditions remain unchanged. Both scenarios resulted in a mean increase of stream TOC concentrations of between 1.5 and 2.5 mg L−1 during the snow free season, which amounts to approximately 15% more TOC export compared to present conditions. Wetter and warmer conditions in the late autumn led to a difference of monthly average TOC of up to 5 mg L−1, suggesting that stream TOC may be particularly susceptible to climate variability during this season.

Criteria are needed for distinguishing naturally acid water from that acidified by air pollution, especially in the organic-rich waters of northern Sweden. The Steady-State Water Chemistry Model (SSWC) was augmented to include organic acidity so that it could predict pre-industrial pH in organic-rich waters. The resulting model predictions of pre-industrial ANC and pH were then tested against diatom predictions of pre-industrial pH and alkalinity in 58 lakes from N. Sweden (after alkalinity was converted to ANC using the CBALK method). The SSWC Model's predictions of pre-industrial lake pH in N. Sweden did not correspond well with the diatom predictions, even when accounting for the uncertainty in the diatom model. This was due to the SSWC's sensitivity to short-term fluctuations in contemporary water chemistry. Thus the SSWC Model is not suitable for judging the acidification of individual lakes in areas such as northern Sweden where the degree of chronic acidification is small, or without a good average value of contemporary water chemistry. These results should be considered when assessing the accuracy of critical loads calculated using SSWC.