Flow regimes are critical for determining physical and biological processes in rivers, and their classification and regionalization traditionally seeks to link patterns of flow to physiographic, climate and other information. There are many approaches to, and rationales for, catchment classification, with those focused on streamflow often seeking to relate a particular response characteristic to a physical property or climatic driver. Rationales include such topics as Prediction in Ungauged Basins (PUB), and providing guidance for model selection in poorly understood hydrological systems. The Annual Daily Hydrograph (ADH) is a first-order easily visualized integrated expression of catchment function, and over many years the average ADH is a distinct hydrological signature that differentiate catchments from each other. In this study, we use t-SNE, a state-of-the-art technique of dimensionality reduction, to classify 17110 ADHs for 304 reference catchments in mountainous Western North America. t-SNE is chosen over other conventional methods of dimensionality reduction (e.g. PCA) as it presents greater separability of ADHs, which are projected on a 2D map where the similarities are evaluated according to their map distance. We then utilize a Deep Learning encoder to upgrade the non-parametric t-SNE to a parametric approach, enhancing its capability to address ’unseen’ samples. Results showed that t-SNE successfully clustered ADHs of similar flow regimes on the 2D map and allowed more accurate classification with KNN. In addition, many compact clusters on the 2D map in the coastal Pacific Northwest suggest information redundancy in the local streamflow network. The t-SNE map provides an intuitive way to visualize the similarity of high-dimensional data of ADHs, groups catchments with like characteristics, and avoids the reliance on subjective hydrometric indicators. This article is protected by copyright. All rights reserved.

Reclamation of wetlands, including peatlands, is legally required in the Athabasca oil sands region following bitumen extraction via surface mining which leaves large open pits that are backfilled with saline tailings waste. Six years of hydrochemical data (2013 – 2018) from the Sandhill Fen Watershed (SFW), a 52-ha upland-peatland catchment that was built upon highly saline soft tailings, were used to evaluate salinity and ion patterns and provide insight on its trajectory. In general, electrical conductivity (EC) increased throughout SFW from 1) reduced inflow and outflow, 2) changes in water table positions and 3) increased mixing of site-wide waters. Salinity has increased site-wide over time as EC increased by an average of 1585 and 2313 µS/cm in the wetland and margins, respectively from 2013 to 2018. The uplands were the only region where EC declined by 1747 µS/cm over the six years. There is also evidence of mixing with underlying tailings waste (Na-Cl dominant) as the chemical composition of SFW waters shifted from largely Ca-dominant in 2013 (> 90%) to Na-dominant by 2018 (> 70%). Based on its current conditions, SFW cannot support freshwater peat-forming bryophytes and is most chemically similar to naturally occurring saline fens. A shift in design strategies is recommended (replicating saline instead of freshwater peatlands) to increase the success of these systems. Changes in site-wide average annual electrical conductivity (µS/cm) from 2013 to 2018. • Electrical conductivity (EC) and Na+ increased site-wide from 2013 to 2018. • Decreased inflow and outflow increased EC and ion accumulation. • Higher water tables diluted wetland EC and increased margin EC (ion mobilization). • Increased mixing with tailings waste shifted waters from Ca +2 to Na + dominant. • Waters have become similar to brackish marshes and saline fens.

Water science data are a valuable asset that both underpins the original research project and bolsters new research questions, particularly in view of the increasingly complex water issues facing Canada and the world. Whilst there is general support for making data more broadly accessible, and a number of water science journals and funding agencies have adopted policies that require researchers to share data in accordance with the FAIR (Findable, Accessible, Interoperable, Reusable) principles, there are still questions about effective management of data to protect their usefulness over time. Incorporating data management practices and standards at the outset of a water science research project will enable researchers to efficiently locate, analyze and use data throughout the project lifecycle, and will ensure the data maintain their value after the project has ended. Here, some common misconceptions about data management are highlighted, along with insights and practical advice to assist established and early career water science researchers as they integrate data management best practices and tools into their research. Freely available tools and training opportunities made available in Canada through Global Water Futures, the Portage Network, Gordon Foundation's DataStream, Compute Canada, and university libraries, among others are compiled. These include webinars, training videos, and individual support for the water science community that together enable researchers to protect their data assets and meet the expectations of journals and funders. The perspectives shared here have been developed as part of the Global Water Futures programme's efforts to improve data management and promote the use of common data practices and standards in the context of water science in Canada. Ten best practices are proposed that may be broadly applicable to other disciplines in the natural sciences and can be adopted and adapted globally. This article is protected by copyright. All rights reserved.

• Boreal forests evaporate considerably more than higher elevation shrub ecosystems. • Forests exist in a growing season water deficit relying on snowmelt recharge. • ET variability declines with increased vegetation cover. • Majority of growing season water for streamflow is generated at higher elevations. • We propose treeline advance will result in drying of northern catchments. As a result of altitude and latitude amplified climate change, widespread changes in vegetation composition, density and distribution have been observed across northern regions. Despite wide documentation of shrub proliferation and treeline advance, few field-based studies have evaluated the hydrological implications of these changes. Quantification of total evapotranspiration (ET) across a range of vegetation gradients is essential for predicting water yield, yet challenging in cold alpine catchments due to heterogeneous land cover, including both boreal forest and shrub taiga ecosystems. Here, we present six years of surface energy balance components and ET dynamics at three sites along an elevational gradient in a subarctic, alpine catchment near Whitehorse, Yukon Territory, Canada. These sites span a gradient of thermal and vegetation regimes, providing a space-for-time comparison for future ecosystem shifts: 1) a low-elevation boreal white spruce forest (~12–20 m), 2) a mid-elevation subalpine taiga comprised of tall, dense willow ( Salix ) and birch ( Betula ) shrubs (~1–3 m) and 3) a high-elevation subalpine taiga with short, sparse shrub cover (<0.75 m) and moss, lichen, and bare rock. Eddy covariance instrumentation ran year-round at the forest and during the growing season at the two shrub sites. Total ET decreased and interannual variability increased with elevation, with mean May to September ET totals of 349 (±3) mm at the forest, 249 (±10) mm at the tall, dense shrub site, and 240 (±26) mm at the short, sparse shrub site. Comparatively, over the same period, ET:R ratios were the highest and most variable at the forest (2.19 ± 0.37) and similar at the tall, dense shrub (1.22 ± 0.09) and short, sparse shrub (1.14 ± 0.05) sites. Our results suggest that advances in treeline will increase overall ET and lower interannual variability; however, the large growing season water deficit at the forest indicates strong reliance on soil moisture from late fall and snowmelt recharge. In contrast, ET was considerably less at the cooler higher elevation shrub sites , which exhibited similar ET losses over 6 years despite differences in shrub height and abundance. ET rates between the two shrub sites were similar throughout the year, except during the peak growing season. Greater interannual variability in ET at the short, sparse shrub site indicates the reduced influence of vegetation controls on ET. Results suggest that predicted changes in vegetation type and structure in northern regions will have a considerable impact on water partitioning and will vary in a complex way in response to changing precipitation timing, phase and magnitude, growing season length, and vegetation snow and rain interactions.

Abstract As a consequence of increasing temperatures, a rapid increase in shrub vegetation has occurred throughout the circumpolar North and is expected to continue. Rates of shrub expansion are highly variable, both at the regional scale and within local study areas. This study uses repeat airborne LiDAR and field surveys to measure changes in shrub vegetation cover along with landscape-scale variations in a well-studied subarctic headwater catchment in Yukon Territory, Canada. Airborne LiDAR surveys were conducted in August 2007 and 2018, whereas vegetation surveys were conducted in summer 2019. Machine learning classification algorithms were used to predict shrub presence/absence in 2018 based on rasterized LiDAR metrics, with the best-performing model applied to the 2007 LiDAR to create binary shrub cover layers to compare between survey years. Results show a 63.3% total increase in detectable shrub cover >= 0.45 m in height between 2007 and 2018, with an average yearly expansion of 5.8%. These changes were compared across terrain derivatives to quantify the influence of topography on shrub expansion. Terrain comparisons show that shrubs are located in and are preferentially expanding into lower and flatter areas near stream networks, at lower slope positions and with a higher potential for topographic wetness. Overall, the findings from this research reinforce the documented increase in pan-Arctic shrub vegetation in recent years, quantify the variation in shrub expansion over terrain derivatives at the landscape scale, and demonstrate the feasibility of using LiDAR to compare changes in shrub properties over time.

We compared stable isotopes of water in plant stem (xylem) water and soil collected over a complete growing season from five well-known long-term study sites in northern/cold regions. These spanned a decreasing temperature gradient from Bruntland Burn (Scotland), Dorset (Canadian Shield), Dry Creek (USA), Krycklan (Sweden), to Wolf Creek (northern Canada). Xylem water was isotopically depleted compared to soil waters, most notably for deuterium. The degree to which potential soil water sources could explain the isotopic composition of xylem water was assessed quantitatively using overlapping polygons to enclose respective data sets when plotted in dual isotope space. At most sites isotopes in xylem water from angiosperms showed a strong overlap with soil water; this was not the case for gymnosperms. In most cases, xylem water composition on a given sampling day could be better explained if soil water composition was considered over longer antecedent periods spanning many months. Xylem water at most sites was usually most dissimilar to soil water in drier summer months, although sites differed in the sequence of change. Open questions remain on why a significant proportion of isotopically depleted water in plant xylem cannot be explained by soil water sources, particularly for gymnosperms. It is recommended that future research focuses on the potential for fractionation to affect water uptake at the soil-root interface, both through effects of exchange between the vapour and liquid phases of soil water and the effects of mycorrhizal interactions. Additionally, in cold regions, evaporation and diffusion of xylem water in winter may be an important process.

Abstract. The interior of western Canada, like many similar cold mid- to high-latitude regions worldwide, is undergoing extensive and rapid climate and environmental change, which may accelerate in the coming decades. Understanding and predicting changes in coupled climate–land–hydrological systems are crucial to society yet limited by lack of understanding of changes in cold-region process responses and interactions, along with their representation in most current-generation land-surface and hydrological models. It is essential to consider the underlying processes and base predictive models on the proper physics, especially under conditions of non-stationarity where the past is no longer a reliable guide to the future and system trajectories can be unexpected. These challenges were forefront in the recently completed Changing Cold Regions Network (CCRN), which assembled and focused a wide range of multi-disciplinary expertise to improve the understanding, diagnosis, and prediction of change over the cold interior of western Canada. CCRN advanced knowledge of fundamental cold-region ecological and hydrological processes through observation and experimentation across a network of highly instrumented research basins and other sites. Significant efforts were made to improve the functionality and process representation, based on this improved understanding, within the fine-scale Cold Regions Hydrological Modelling (CRHM) platform and the large-scale Modélisation Environmentale Communautaire (MEC) – Surface and Hydrology (MESH) model. These models were, and continue to be, applied under past and projected future climates and under current and expected future land and vegetation cover configurations to diagnose historical change and predict possible future hydrological responses. This second of two articles synthesizes the nature and understanding of cold-region processes and Earth system responses to future climate, as advanced by CCRN. These include changing precipitation and moisture feedbacks to the atmosphere; altered snow regimes, changing balance of snowfall and rainfall, and glacier loss; vegetation responses to climate and the loss of ecosystem resilience to wildfire and disturbance; thawing permafrost and its influence on landscapes and hydrology; groundwater storage and cycling and its connections to surface water; and stream and river discharge as influenced by the various drivers of hydrological change. Collective insights, expert elicitation, and model application are used to provide a synthesis of this change over the CCRN region for the late 21st century.

Currently, post-mining landscape plans in the Athabasca Oil Sand Region include large watersheds terminating in pit lakes. In 2012, Base Mine Lake (BML), was constructed with the aim of demonstrating technologies associated with lake reclamation in the region. This paper examines the first 6.5 years of lake-atmosphere energy and carbon exchange. Energetically, BML behaved similar to other northern lakes, storing large quantities of heat in the spring and releasing it in the fall as sensible and latent heat fluxes. At various times a hydrocarbon sheen formed on the lake, which may have suppressed evaporation. However, simple linear relationships failed to statistically quantify the impacts and more comprehensive modelling of the variability may be required. At daily scales, variability in evaporation was well explained by the product of vapour pressure deficit and wind speed as well as the available energy (R2 = 0.74), while sensible heat was explained by the product of wind speed and the difference in air and surface temperature as well as available energy (R2 = 0.85). Spring CH4 fluxes were high, particularly around ice melt, with a maximum flux of 3.3 g m-2 day-1. Otherwise fluxes were low, except during irregular periods. The peak flux of these periods occurred following ~58 h of continuously falling pressure, relating cyclone activity to these large periods of methane emissions. Annually, CO2 and CH4 fluxes were initially high, with median fluxes of 231 mg CO2 m-2 h-1 and 23 mg CH4 m-2 h-1 in 2014. However, the median fluxes reduced quickly and over the least three years of the study (2017 through 2019) the median fluxes declined to 36 mg CO2 m-2 h-1 and 10 mg CH4 m-2 h-1. Overall, BML behaves similar to other boreal lake ecosystems with above average carbon fluxes compared to other constructed reservoirs.

Alpine headwaters in subarctic regions are particularly sensitive to climate change, yet there is little information on stream thermal regimes in these areas and how they might respond to global warming. In this paper, we characterize and compare the hydrological and thermal regimes of two subarctic headwater alpine streams within an empirical framework. The streams investigated are located within two adjacent catchments with similar geology, size, elevation and landscape, Granger Creek (GC) and Buckbrush Creek (BB), which are part of the Wolf Creek Research Basin in the Yukon Territory, Canada. Hydrometeorological and high‐resolution stream temperature data were collected throughout summer 2016. Both sites exhibited a flow regime typical of cold alpine headwater catchments influenced by frozen ground and permafrost. Comparatively, GC was characterized by a flashier response with more extreme flows, than BB. In both sites, stream temperature was highly variable and very responsive to short‐term changes in climatic conditions. On average, stream temperature in BB was slightly higher than in GC (respectively 5.8 and 5.7°C), but less variable (average difference between 75th and 25th quantiles of 1.6 and 2.0°C). Regression analysis between mean daily air and stream temperature suggested that a greater relative (to stream flow) groundwater contribution in BB could more effectively buffer atmospheric fluctuations. Heat fluxes were derived and utilized to assess their relative contribution to the energy balance. Overall, non‐advective fluxes followed a daily pattern highly correlated to short‐wave radiation. G1enerally, solar radiation and latent heat were respectively the most important heat source and sink, while air–water interface processes were major factors driving nighttime stream temperature fluctuations.

Abstract Water storage dynamics modulate fluxes within catchments, control the rainfall-runoff response and regulate the velocity of water particles through mixing associated processes. Tracer-aided models are useful tools for tracking the interactions between catchment storage and fluxes, as they can capture both the celerity of the runoff response and the velocity of water particles revealed by tracer dynamics. The phase-space reconstruction of modelled systems can help in this regard; it traces the evolution of a dynamic system from a known initial state as phase trajectories in response to inputs. In this study, we compared the modelled storage-flux dynamics obtained from the application of a spatially distributed tracer-aided hydrological model (STARR) in five contrasting long-term research catchments with varying degrees of snow influence. The models were calibrated using a consistent multivariate methodology based on discharge, isotope composition and snowpack water equivalent. Analysis of extracted modelled storage dynamics gave insights into the system functioning. Large volumes of total stored water needed to be invoked at most sites to reconcile celerity and travel times to match observe discharge and isotope responses. This is because changes in dynamic storage from water balance considerations are small when compared to volume of storage necessary for observed tracer dampening. In the phase-space diagrams, the rates of storage change gave insights into the relative storage volume and seasonal catchment functioning. The storage increase was dominated by hydroclimatic inputs; thus, it presented a stochastic response. Furthermore, depending on the dominance of snow or rainfall inputs, catchments had different seasonal responses in storage dynamics. Decreases in storage were more predictable and reflected the efficiency of catchment drainage, yet at lower storages the influence of ET was also evident. Activation of flow paths due to overland and near-surface flows resulted in non-linearity of catchment functioning largely at high storage states. The storage-discharge relationships generally showed a non-linear distribution, with more scattered states during wettest condition. In turn, all the catchments exhibited an inverse storage effect, with modelled water ages decreasing with increasing storage as lateral flow paths were activated. Insights from this inter-comparison of storage-flux-age dynamics show the benefits of tracer-aided hydrological models in exploring their interactions at well-instrumented sites to better understand hydrological functioning of contrasting catchments.

Mine reclamation in the Athabasca oil sands region Canada, is required by law where companies must reconstruct disturbed landscapes into functioning ecosystems such as forests, wetlands and lakes that existed in the Boreal landscape prior to mining. Winter is a major hydrological factor in this region as snow covers the landscape for 5 to 6 months and is ~25% of the annual precipitation, yet few studies have explored the influence of winter processes on the hydrology of constructed watersheds. One year (2017-2018) of intensive snow hydrology measurements are supplemented with six years (2013-2018) of meteorological measurements from the constructed Sandhill Fen Watershed to: 1) understand snow accumulation and redistribution, snowmelt timing, rate and partitioning, 2) apply a physically-based model for simulating winter processes on hillslopes and 3) evaluate the impact of soil prescriptions and climate change projections on winter processes in reclaimed systems. The 2017-2018 snow season was between November and April and SWE ranged between 40-140 mm. Snow distribution was primarily influenced by topography with little influence of snow trapping from developing vegetation. Snow accumulation was most variable on hillslopes and redistribution was driven by slope position, with SWE greatest at the base of slopes and decreased towards crests. Snowmelt on hillslopes was controlled by slope aspect, as snow declined rapidly on west and south-facing slopes, compared to east and north-facing slopes. Unlike results previously reported on constructed uplands, snowmelt runoff from uplands was much less (~30%), highlighting the influence of different construction materials. Model simulations indicate that antecedent soil moisture and soil temperature have a large influence on partitioning snowmelt over a range of observed conditions. Under a warmer and wetter climate, average annual peak SWE and snow season duration could decline up to 52 % and up to 61 days, respectively while snowmelt runoff ceases completely under the warmest scenarios. Results suggest considerable future variability in snowmelt runoff from hillslopes, yet soil properties can be used to enhance vertical or lateral flows.

The response of evapotranspiration (ET) to warming is of critical importance to the water and carbon cycle of the boreal biome, a mosaic of land cover types dominated by forests and peatlands. The effect of warming-induced vapour pressure deficit (VPD) increases on boreal ET remains poorly understood because peatlands are not specifically represented as plant functional types in Earth system models. Here we show that peatland ET increases more than forest ET with increasing VPD using observations from 95 eddy covariance tower sites. At high VPD of more than 2 kPa, peatland ET exceeds forest ET by up to 30%. Future (2091–2100) mid-growing season peatland ET is estimated to exceed forest ET by over 20% in about one-third of the boreal biome for RCP4.5 and about two-thirds for RCP8.5. Peatland-specific ET responses to VPD should therefore be included in Earth system models to avoid biases in water and carbon cycle projections.

Peatlands and forests cover large areas of the boreal biome and are critical for global climate regulation. They also regulate regional climate through heat and water vapour exchange with the atmosphere. Understanding how land-atmosphere interactions in peatlands differ from forests may therefore be crucial for modelling boreal climate system dynamics and for assessing climate benefits of peatland conservation and restoration. To assess the biophysical impacts of peatlands and forests on peak growing season air temperature and humidity, we analysed surface energy fluxes and albedo from 35 peatlands and 37 evergreen needleleaf forests - the dominant boreal forest type - and simulated air temperature and vapour pressure deficit (VPD) over hypothetical homogeneous peatland and forest landscapes. We ran an evapotranspiration model using land surface parameters derived from energy flux observations and coupled an analytical solution for the surface energy balance to an atmospheric boundary layer (ABL) model. We found that peatlands, compared to forests, are characterized by higher growing season albedo, lower aerodynamic conductance, and higher surface conductance for an equivalent VPD. This combination of peatland surface properties results in a ∼20% decrease in afternoon ABL height, a cooling (from 1.7 to 2.5 °C) in afternoon air temperatures, and a decrease in afternoon VPD (from 0.4 to 0.7 kPa) for peatland landscapes compared to forest landscapes. These biophysical climate impacts of peatlands are most pronounced at lower latitudes (∼45°N) and decrease toward the northern limit of the boreal biome (∼70°N). Thus, boreal peatlands have the potential to mitigate the effect of regional climate warming during the growing season. The biophysical climate mitigation potential of peatlands needs to be accounted for when projecting the future climate of the boreal biome, when assessing the climate benefits of conserving pristine boreal peatlands, and when restoring peatlands that have experienced peatland drainage and mining. © 2020 The Author(s). Published by IOP Publishing Ltd. (Less)

Abstract. Permafrost strongly controls hydrological processes in cold regions, and our understanding of how changes in seasonal and perennial frozen ground disposition and linked storage dynamics affects runoff generation processes remains limited. Storage dynamics and water redistribution are influenced by the seasonal variability and spatial heterogeneity of frozen ground, snow accumulation and melt. Stable isotopes provide a potentially useful technique to quantify the dynamics of water sources, flow paths and ages; yet few studies have employed isotope data in permafrost-influenced catchments. Here, we applied the conceptual model STARR (Spatially distributed Tracer-Aided Rainfall-Runoff model), which facilitates fully distributed simulations of hydrological storage dynamics and runoff processes, isotopic composition and water ages. We adapted this model to a subarctic catchment in Yukon Territory, Canada, with a time-variable implementation of field capacity to include the influence of thaw dynamics. A multi-criteria calibration based on stream flow, snow water equivalent and isotopes was applied to three years of data. The integration of isotope data in the spatially distributed model provided the basis to quantify spatio-temporal dynamics of water storage and ages, emphasizing the importance of thaw layer dynamics in mixing and damping the melt signal. By using the model conceptualisation of spatially and temporally variant storage, this study demonstrates the ability of tracer-aided modelling to capture thaw layer dynamics that cause mixing and damping of the isotopic melt signal.

Abstract. High latitude environments store approximately half of the global organic carbon pool in peatlands, organic soils and permafrost while large arctic rivers convey an estimated 18–50 Tg C a−1 to the Arctic Ocean. Warming trends associated with climate change affect dissolved organic carbon (DOC) export from terrestrial to riverine environments. However, there is limited consensus as to whether exports will increase or decrease due to complex interactions between climate, soils, vegetation, and associated production, mobilization and transport processes. A large body of research has focused on large river system DOC and DOM lability and observed trends conserved across years, whereas investigation at smaller watershed scales show that thermokarst and fire have a transient impact on hydrologically-mediated solute transport. This study, located in the Wolf Creek Research Basin situated ~ 20 km south of Whitehorse, YT, Canada, utilises a nested design to assess seasonal and annual patterns of DOC and DOM composition across diverse landscape types (headwater, wetland, lake) and watershed scales. Peak DOC concentration and export occurred during freshet per most northern watersheds, however, peaks were lower than a decade ago at the headwater site Granger Creek. DOM composition was most variable during freshet with high A254, SUVA254 and low FI and BIX. DOM composition was relatively insensitive to flow variation during summer and fall. The influence of increasing watershed scale and downstream mixing of landscape contributions was an overall dampening of DOC concentrations and optical indices with increasing groundwater contribution. Forecasted vegetation shifts, permafrost thaw and other changes due to climate change may alter DOM sources from predominantly organic soils to decomposing vegetation, and facilitate transport through deeper flow pathways with an enhanced groundwater role.

Abstract. A set of hydrometeorological data is presented in this paper, which can be used to characterize the hydrometeorology and climate of a subarctic mountain basin and has proven particularly useful for forcing hydrological models and assessing their performance in capturing hydrological processes in subarctic alpine environments. The forcing dataset includes daily precipitation, hourly air temperature, humidity, wind, solar and net radiation, soil temperature, and geographical information system data. The model performance assessment data include snow depth and snow water equivalent, streamflow, soil moisture, and water level in a groundwater well. This dataset was recorded at different elevation bands in Wolf Creek Research Basin, near Whitehorse, Yukon Territory, Canada, representing forest, shrub tundra, and alpine tundra biomes from 1993 through 2014. Measurements continue through 2018 and are planned for the future at this basin and will be updated to the data website. The database presented and described in this article is available for download at https://doi.org/10.20383/101.0113.

Abstract. High-latitude environments store approximately half of the global organic carbon pool in peatlands, organic soils and permafrost, while large Arctic rivers convey an estimated 18–50 Tg C a−1 to the Arctic Ocean. Warming trends associated with climate change affect dissolved organic carbon (DOC) export from terrestrial to riverine environments. However, there is limited consensus as to whether exports will increase or decrease due to complex interactions between climate, soils, vegetation, and associated production, mobilization and transport processes. A large body of research has focused on large river system DOC and dissolved organic matter (DOM) lability and observed trends conserved across years, whereas investigation at smaller watershed scales show that thermokarst and fire have a transient impact on hydrologically mediated solute transport. This study, located in the Wolf Creek Research Basin situated ∼20 km south of Whitehorse, YT, Canada, utilizes a nested design to assess seasonal and annual patterns of DOC and DOM composition across diverse landscape types (headwater, wetland and lake) and watershed scales. Peak DOC concentration and export occurred during freshet, as is the case in most northern watersheds; however, peaks were lower than a decade ago at the headwater site Granger Creek. DOM composition was most variable during freshet with high A254 and SUVA254 and low FI and BIX. DOM composition was relatively insensitive to flow variation during summer and fall. The influence of increasing watershed scale and downstream mixing of landscape contributions was an overall dampening of DOC concentrations and optical indices with increasing groundwater contribution. Forecasted vegetation shifts, enhanced permafrost and seasonal thaw, earlier snowmelt, increased rainfall and other projected climate-driven changes will alter DOM sources and transport pathways. The results from this study support a projected shift from predominantly organic soils (high aromaticity and less fresh) to decomposing vegetation (more fresh and lower aromaticity). These changes may also facilitate flow and transport via deeper flow pathways and enhance groundwater contributions to runoff.

We assessed the hydrological implications of climate effects on vegetation phenology in northern environments by fusion of data from remote-sensing and local catchment monitoring. Studies using satellite data have shown earlier and later dates for the start (SOS) and end of growing seasons (EOS), respectively, in the Northern Hemisphere over the last 3 decades. However, estimates of the change greatly depend on the satellite data utilized. Validation with experimental data on climate-vegetation-hydrology interactions requires long-term observations of multiple variables which are rare and usually restricted to small catchments. In this study, we used two NDVI (normalized difference vegetation index) products (at ~25 & 0.5 km spatial resolutions) to infer SOS and EOS for six northern catchments, and then investigated the likely climate impacts on phenology change and consequent effects on catchment water yield, using both assimilated data (GLDAS: global land data assimilation system) and direct catchment observations. The major findings are: (1) The assimilated air temperature compared well with catchment observations (regression slopes and R2 close to 1), whereas underestimations of summer rainstorms resulted in overall underestimations of precipitation (regression slopes of 0.3-0.7, R2 ≥ 0.46). (2) The two NDVI products inferred different vegetation phenology characteristics. (3) Increased mean pre-season temperature significantly influenced the advance of SOS and delay of EOS. The precipitation influence was weaker, but delayed SOS corresponding to increased pre-season precipitation at most sites can be related to later snow melting. (4) Decreased catchment streamflow over the last 15 years could be related to the advance in SOS and extension of growing seasons. Greater streamflow reductions in the cold sites than the warm ones imply stronger climate warming impacts on vegetation and hydrology in colder northerly environments. The methods used in this study have potential for better understanding interactions between vegetation, climate and hydrology in observation-scarce regions.

Abstract Northern peatlands contain up to 20% of the ∼3000 Pg of global soil organic carbon. Carbon-rich peatlands cover upwards of 65% of the landscape in northern Canada where resource extraction activities disturb both the carbon pools and the future carbon sequestration capacity of the landscape. Previous estimates of the carbon losses from this disturbance predict a complete loss of the region’s peatland carbon pool. Mining industries operating in these sensitive environments have recently begun constructing closure landscapes which are intended to develop carbon cycle processes similar to undisturbed northern peatlands. This study investigates eddy covariance fluxes of carbon dioxide (CO2) at one of Canada’s first fully constructed boreal plains watersheds, the Sandhill Fen Watershed. During the first three years since inception, only the lowland region had an annual net ecosystem exchange of CO2 (NEE) indicative of increasing CO2 sink potential. The lowland region was characterized by saturated salvaged peat soils, standing water, thriving communities of Typha and Carex spp. and was a net CO2 sink of 77 g C m−2 in the third year. At the same time the upland and the midland regions characterized by moist salvaged peat soils and a mix of herbaceous, shrub and planted Picea glauca and Pinus banksiana remained net sources of CO2. Despite similar rates of gross primary production, ecosystem and plot-level respiration rates in the lowland were significantly lower than in the midland region, likely due to very low reduction potentials within the lowland’s saturated soils. With no other significant outflows of carbon, the lowland of the Sandhill Fen Watershed may be in the early stages of organic matter accumulation. Due to limited oxidation of the salvaged peat substrate in the lowland region, wetland reclamation employing these techniques may reduce the disturbance loss of the carbon pool in the boreal plains.

Abstract. Permafrost strongly controls hydrological processes in cold regions. Our understanding of how changes in seasonal and perennial frozen ground disposition and linked storage dynamics affect runoff generation processes remains limited. Storage dynamics and water redistribution are influenced by the seasonal variability and spatial heterogeneity of frozen ground, snow accumulation and melt. Stable isotopes are potentially useful for quantifying the dynamics of water sources, flow paths and ages, yet few studies have employed isotope data in permafrost-influenced catchments. Here, we applied the conceptual model STARR (the Spatially distributed Tracer-Aided Rainfall–Runoff model), which facilitates fully distributed simulations of hydrological storage dynamics and runoff processes, isotopic composition and water ages. We adapted this model for a subarctic catchment in Yukon Territory, Canada, with a time-variable implementation of field capacity to include the influence of thaw dynamics. A multi-criteria calibration based on stream flow, snow water equivalent and isotopes was applied to 3 years of data. The integration of isotope data in the spatially distributed model provided the basis for quantifying spatio-temporal dynamics of water storage and ages, emphasizing the importance of thaw layer dynamics in mixing and damping the melt signal. By using the model conceptualization of spatially and temporally variable storage, this study demonstrates the ability of tracer-aided modelling to capture thaw layer dynamics that cause mixing and damping of the isotopic melt signal.

Bitumen extraction via surface mining in the Athabasca Oil Sands Region results in permanent alteration of boreal forests and wetlands. As part of their legal requirements, oil companies must reclaim disturbed landscapes into functioning ecosystems. Despite considerable work establishing upland forests, only two pilot wetland-peatland systems integrated within a watershed have been constructed to date. Peatland reclamation is challenging as it requires complete reconstruction with few guidelines or previous work in this region. Furthermore, the variable sub-humid climate and salinity of tailings materials present additional challenges. In 2012, Syncrude Canada Ltd. constructed a 52-ha pilot upland-wetland system, the Sandhill Fen Watershed, which was designed with a pump and underdrain system to provide freshwater and enhance drainage to limit salinization from underlying soft tailings materials that have elevated electrical conductivity (EC) and Na+. The objective of this research is to evaluate the hydrochemical response of a constructed wetland to variations in hydrology and water management with respect to water sources, flow pathways and major chemical transformations in the three years following commissioning. Results suggest that active water management practices in 2013 kept EC relatively low, with most wetland sites <1000 μS/cm with Na+ concentrations <250 mg/L. With limited management in 2014 and 2015, the EC increased in the wetland to >1000 μS/cm in 2014 and >2000 μS/cm in 2015. The most notable change was the emergence of several Na+ enriched zones in the margins. Here, Na+ concentrations were two to three times higher than other sites. Stable isotopes of water support that the Na+ enriched areas arise from underlying process-affected water in the tailings, providing evidence of its upward transport and seepage under a natural hydrologic regime. In future years, salinity is expected to evolve in its flow pathways and diffusion, yet the timeline and extent of these changes are uncertain.

We thank Audrey Innes for isotope analysis at University of Aberdeen for Bruntland Burn and Krycklan sites, Johannes Tiwari (SLU) for isotope sampling in Krycklan, Pernilla Lofvenius (SLU) for providing PET data for Krycklan (via SITES), and Jeff McDonnell and Kim Janzen (University of Saskatchewan) for soil water isotope analysis for the Dorset and Wolf Creek sites. The Krycklan part was funded by the KAW Branch-Point project. We acknowledge the funding from the European Research Council (ERC, project GA 335910 VeWa). We thank the Editor and three anonymous reviewers for their critical comments during the peer-review process.

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

Time series remote sensing vegetation indices derived from SPOT 5 data are compared with vegetation structure and eddy covariance flux data at 15 dry to wet reclamation and reference sites within the Oil Sands region of Alberta, Canada. This comprehensive analysis examines the linkages between indicators of ecosystem function and change trajectories observed both at the plot level and within pixels. Using SPOT imagery, we find that higher spatial resolution datasets (e.g. 10 m) improves the relationship between vegetation indices and structural measurements compared with interpolated (lower resolution) pixels. The simple ratio (SR) vegetation index performs best when compared with stem density-based indicators (R2 = 0.65; p < 0.00), while the normalised difference vegetation index (NDVI) and soil adjusted vegetation index (SAVI) are most comparable to foliage indicators (leaf area index (LAI) and canopy cover (R2 = 0.52-0.78; p > 0.02). Fluxes (net ecosystem production (NEP) and gross ecosystem production (GEP)) are most related to NDVI and SAVI when these are interpolated to larger 20 m × 20 m pixels (R2 = 0.44-0.50; p < 0.00). As expected, decreased sensitivity of NDVI is problematic for sites with LAI > 3 m2 m-2, making this index more appropriate for newly regenerating reclamation areas. For sites with LAI < 3 m2 m-2, trajectories of vegetation change can be mapped over time and are within 2.7% and 3.3% of annual measured LAI changes observed at most sites. This study demonstrates the utility of remote sensing in combination with field and eddy covariance data for monitoring and scaling of reclaimed and reference site productivity within and beyond the Oil Sands Region of western Canada.

Surface mining in northern Alberta transforms wetlands and forests into open pits, tailings ponds, and overburden. As part of their license to operate, mine operators are required to reclaim this altered landscape to a predisturbance capacity. In 2012, Syncrude Canada Limited constructed one of the first of two reclaimed wetlands, the Sandhill Fen Watershed (SFW), to evaluate wetland reclamation strategies. SFW is a 52-ha system atop soft-tailings that includes an inflow/outflow pump system, underdrains, upland hummocks, and a fen lowland. In this study, water table dynamics of the fen lowland were evaluated in the 2 years following commissioning (2014–2015) to assess whether this newly constructed watershed has hydrological conditions that facilitate hydric soils with water table regimes similar to reference systems. Results indicate that the location and hydrophysical properties of placed materials control water table responses to both water management and precipitation. This differential water table response in the SFW lowland drove lateral fluxes between adjacent landforms, suggesting periods of intermittent water supply from uplands to wetlands along hummock margins. As in natural systems, the lowland fen exhibited several lateral flow reversals over the 2 years depending upon water level. Water tables on-average were greater than those observed in natural analogues. Comparison during these first 2 years following commissioning contribute to the increasing insight as to how construction and management practices support reclamation postmining.