The regional variability in tundra and boreal carbon dioxide (CO₂) fluxes can be high, complicating efforts to quantify sink‐source patterns across the entire region. Statistical models are increasingly used to predict (i.e., upscale) CO₂ fluxes across large spatial domains, but the reliability of different modeling techniques, each with different specifications and assumptions, has not been assessed in detail. Here, we compile eddy covariance and chamber measurements of annual and growing season CO₂ fluxes of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) during 1990–2015 from 148 terrestrial high‐latitude (i.e., tundra and boreal) sites to analyze the spatial patterns and drivers of CO₂ fluxes and test the accuracy and uncertainty of different statistical models. CO₂ fluxes were upscaled at relatively high spatial resolution (1 km²) across the high‐latitude region using five commonly used statistical models and their ensemble, that is, the median of all five models, using climatic, vegetation, and soil predictors. We found the performance of machine learning and ensemble predictions to outperform traditional regression methods. We also found the predictive performance of NEE‐focused models to be low, relative to models predicting GPP and ER. Our data compilation and ensemble predictions showed that CO₂ sink strength was larger in the boreal biome (observed and predicted average annual NEE −46 and −29 g C m⁻² yr⁻¹, respectively) compared to tundra (average annual NEE +10 and −2 g C m⁻² yr⁻¹). This pattern was associated with large spatial variability, reflecting local heterogeneity in soil organic carbon stocks, climate, and vegetation productivity. The terrestrial ecosystem CO₂ budget, estimated using the annual NEE ensemble prediction, suggests the high‐latitude region was on average an annual CO₂ sink during 1990–2015, although uncertainty remains high.

The complete landscape surface of the active Mackenzie River Delta (13,135 km²) was manually partitioned into discrete lakes (3331 km²), channels (1744 km²), wetlands (1614 km²), and dry floodplain area (6446 km²) via GIS analysis of digital topographic maps recently available for the system. The census total of lakes (49,046) is almost twice as large as prior estimates. Using this new information, total lake volume in the delta during the post river flooding period is estimated as 5.4 km³. Total floodwater storage in the delta lakes and floodplain at peak water levels is estimated at 25.8 km³ and thus is equivalent to about 47% of Mackenzie River flow (55.4 km³ yr⁻¹) during the high‐discharge period of delta breakup. During this period the stored river water can be envisioned in the form of a thin layer of water (2.3 m thick on average) spread out over 11,200 km² of lakes and flooded vegetation and exposed to 24 h d⁻¹ solar irradiance. Consequently, this temporarily stored water has significant potential to affect the composition of river water flowing to the Beaufort Shelf as it recedes to the river channels after the flood peak.

Using a whole-watershed approach and a combination of historical, contemporary, modeled and paleolimnological datasets, we show that the High Arctic’s largest lake by volume (Lake Hazen) has succumbed to climate warming with only a ~1 °C relative increase in summer air temperatures. This warming deepened the soil active layer and triggered large mass losses from the watershed’s glaciers, resulting in a ~10 times increase in delivery of glacial meltwaters, sediment, organic carbon and legacy contaminants to Lake Hazen, a >70% decrease in lake water residence time, and near certainty of summer ice-free conditions. Concomitantly, the community assemblage of diatom primary producers in the lake shifted dramatically with declining ice cover, from shoreline benthic to open-water planktonic species, and the physiological condition of the only fish species in the lake, Arctic Char, declined significantly. Collectively, these changes place Lake Hazen in a biogeochemical, limnological and ecological regime unprecedented within the past ~300 years.

The Canadian Arctic has vast freshwater resources, and fish are important in the diet of many Northerners. Mercury is a contaminant of concern because of its potential toxicity and elevated bioaccumulation in some fish populations. Over the last decade, significant advances have been made in characterizing the cycling and fate of mercury in these freshwater environments. Large amounts of new data on concentrations, speciation and fluxes of Hg are provided and summarized for water and sediment, which were virtually absent for the Canadian Arctic a decade ago. The biogeochemical processes that control the speciation of mercury remain poorly resolved, including the sites and controls of methylmercury production. Food web studies have examined the roles of Hg uptake, trophic transfer, and diet for Hg bioaccumulation in fish, and, in particular, advances have been made in identifying determinants of mercury levels in lake-dwelling and sea-run forms of Arctic char. In a comparison of common freshwater fish species that were sampled across the Canadian Arctic between 2002 and 2009, no geographic patterns or regional hotspots were evident. Over the last two to four decades, Hg concentrations have increased in some monitored populations of fish in the Mackenzie River Basin while other populations from the Yukon and Nunavut showed no change or a slight decline. The different Hg trends indicate that the drivers of temporal change may be regional or habitat-specific. The Canadian Arctic is undergoing profound environmental change, and preliminary evidence suggests that it may be impacting the cycling and bioaccumulation of mercury. Further research is needed to investigate climate change impacts on the Hg cycle as well as biogeochemical controls of methylmercury production and the processes leading to increasing Hg levels in some fish populations in the Canadian Arctic.