This study presents an analysis of the observed inter-annual variability and inter-decadal trends in river discharge across northern Canada for 1964–2013. The 42 rivers chosen for this study span a combined gauged area of 5.26 × 10⁶ km² and are selected based on data availability and quality, gauged area and record length. Inter-annual variability in river discharge is greatest for the eastern Arctic Ocean (coefficient of variation, CV = 16 %) due to the Caniapiscau River diversion into the La Grande Rivière system for enhanced hydropower production. Variability is lowest for the study area as a whole (CV = 7 %). Based on the Mann–Kendall test (MKT), no significant (p > 0.05) trend in annual discharge from 1964 to 2013 is observed in the Bering Sea, western Arctic Ocean, western Hudson and James Bay, and Labrador Sea; for northern Canada as a whole, however, a statistically significant (p < 0.05) decline of 102.8 km³ 25 yr⁻¹ in discharge occurs over the first half of the study period followed by a statistically significant (p < 0.05) increase of 208.8 km³ 25 yr⁻¹ in the latter half. Increasing (decreasing) trends in river discharge to the eastern Hudson and James Bay (eastern Arctic Ocean) are largely explained by the Caniapiscau diversion to the La Grande Rivière system. Strong regional variations in seasonal trends of river discharge are observed, with overall winter (summer) flows increasing (decreasing, with the exception of the most recent decade) partly due to flow regulation and storage for enhanced hydropower production along the Hudson and James Bay, the eastern Arctic Ocean and Labrador Sea. Flow regulation also suppresses the natural variability of river discharge, particularly during cold seasons.

The pan-Arctic domain is undergoing some of Earth’s most rapid and significant changes resulting from anthropogenic and climate-induced alteration of freshwater distribution. Changes in terrestrial freshwater discharge entering the Arctic Basin from pan-Arctic watersheds significantly impact oceanic circulation and sea ice dynamics. Historical streamflow records in high-latitude basins are often discontinuous (seasonal or with large temporal gaps) or sparse (poor spatial coverage), however, making trends from observed records difficult to quantify. Our objectives were to generate a more continuous 90-year record (1981–2070) of spatially distributed freshwater flux for the Arctic Basin (all Arctic draining rivers, including the Yukon), suitable for forcing ocean models, and to analyze the changing simulated trends in freshwater discharge across the domain. We established these data as valid during the historical period (1971–2015) and then used projected futures (preserving uncertainty by running a coupled climate-hydrologic ensemble) to analyze long-term (2021–2070) trends for major Arctic draining rivers. When compared to historic trends reported in the literature, we find that trends are projected to nearly double by 2070, with river discharge to the Arctic Basin increasing by 22% (on average) by 2070. We also find a significant trend toward earlier onset of spring freshet and a general flattening of the average annual hydrograph, with a trend toward decreasing seasonality of Arctic freshwater discharge with climate change and regulation combined. The coupled climate-hydrologic ensemble was then used to force an ocean circulation model to simulate freshwater content and thermohaline circulation. This research provides the marine research community with a daily time series of historic and projected freshwater discharge suitable for forcing sea ice and ocean models. Although important, this work is only a first step in mapping the impacts of climate change on the pan-Arctic region.

Tracer-aided hydrological models are becoming increasingly popular tools as they assist with process understanding and source separation, which facilitates model calibration and diagnosis of model uncertainty (Tetzlaff et al., 2015; Klaus and McDonnell, 2013). Data availability in high-latitude regions, however, proves to be a major challenge associated with this type of application (Tetzlaff et al., 2015). Models require a time series of isotopes in precipitation (δ¹⁸O_ppt) to drive simulations, and throughout much of the world – particularly in sparsely populated high-latitude regions – these data are not widely available. Here we investigate the impact that choice of precipitation isotope product (δ¹⁸O_ppt) has on simulations of streamflow, δ¹⁸O in streamflow (δ¹⁸O_SF), resulting hydrograph separations, and model parameters. In a high-latitude, data-sparse, seasonal basin (Fort Simpson, NWT, Canada), we assess three precipitation isotope products of different spatial and temporal resolutions (i.e. semi-annual static, seasonal KPN43, and daily bias-corrected REMOiso), and apply them to force the isoWATFLOOD tracer-aided hydrologic model. Total simulated streamflow is not significantly impacted by choice of δ¹⁸O_ppt product; however, simulated isotopes in streamflow (δ¹⁸O_SF) and the internal apportionment of water (driven by model parameterization) are impacted. The highest-resolution product (REMOiso) was distinct from the two lower-resolution products (KPN43 and static), but could not be verified as correct due to a lack of daily δ¹⁸O_ppt observations. The resolution of δ¹⁸O_ppt impacts model parameterization and seasonal hydrograph separations, producing notable differences among simulations following large snowmelt and rainfall events when event compositions differ significantly from δ¹⁸O_SF. Capturing and preserving the spatial variability in δ¹⁸O_ppt using distributed tracer-aided models is important because this variability impacts model parameterization. We achieve an understanding of tracer-aided modelling and its application in high-latitude regions with limited δ¹⁸O_ppt observations, and the value such models have in defining modelling uncertainty. In this study, application of a tracer-aided model is able to identify simulations with improved internal process representation, reinforcing the fact that tracer-aided modelling approaches assist with resolving hydrograph component contributions and work towards diagnosing equifinality.