We present long‐term changes in Norwegian lake water quality across regional gradients in atmospheric pollution, air temperature, hydrology, and vegetation using (a) a national representative lake survey carried out in 1995 and 2019 (ThousandLakes), and (b) an annual lake survey from acid‐sensitive catchments (78 lakes, TrendLakes) from 1990 to 2020. Our analysis encompasses all major chemical constituents, for example, anions and cations, dissolved organic matter (DOM), nutrients, iron (Fe), and silicate (SiO₂). During these decades, environmental changes included declines in sulfur (S) and nitrogen (N) deposition, climate warming, and increase in forest biomass. Strong chemical recovery from acidification is found, attributed to large reductions in atmospheric deposition, moderated by catchment processing from land use and climate change. Browning counteracted chemical recovery in some regions, while Ca increased unexpectedly. We suggest that increased weathering, from enhanced terrestrial productivity, is an important driver of increased Ca—substantiated by widespread, substantial increases in SiO₂. Light‐ and nutrient‐limitation has become more prevalent, indicated by higher DOM, lower nitrate (NO₃), and lower NO₃to total phosphorous ratios. Declines in lake NO₃ occurred independently of N deposition, suggesting increased catchment N retention, possibly from increased terrestrial productivity. We conclude that decreased air pollution continues to be a dominant driver of long‐term trends in lake chemistry, but climate‐induced increase in terrestrial weathering processes, governed by increased biomass, is likely to have an increasing impact on future lake acidity, nutrient, and light status, that may cascade along the aquatic continuum from rivers to the coast.

In lakes, the rates of gross primary production (GPP), ecosystem respiration (R), and net ecosystem production (NEP) are often controlled by resource availability. Herein, we explore how catchment vs. within lake predictors of metabolism compare using data from 16 lakes spanning 39°N to 64°N, a range of inflowing streams, and trophic status. For each lake, we combined stream loads of dissolved organic carbon (DOC), total nitrogen (TN), and total phosphorus (TP) with lake DOC, TN, and TP concentrations and high frequency in situ monitoring of dissolved oxygen. We found that stream load stoichiometry indicated lake stoichiometry for C : N and C : P (r² = 0.74 and r² = 0.84, respectively), but not for N : P (r² = 0.04). As we found a strong positive correlation between TN and TP, we only used TP in our statistical models. For the catchment model, GPP and R were best predicted by DOC load, TP load, and load N : P (R² = 0.85 and R² = 0.82, respectively). For the lake model, GPP and R were best predicted by TP concentrations (R² = 0.86 and R² = 0.67, respectively). The inclusion of N : P in the catchment model, but not the lake model, suggests that both N and P regulate metabolism and that organisms may be responding more strongly to catchment inputs than lake resources. Our models predicted NEP poorly, though it is unclear why. Overall, our work stresses the importance of characterizing lake catchment loads to predict metabolic rates, a result that may be particularly important in catchments experiencing changing hydrologic regimes related to global environmental change.

Rising organic charge in northern freshwaters is attributed to increasing levels of dissolved natural organic matter (DNOM) and changes in water chemistry. Organic charge concentration may be determined through charge balance calculations (Org.⁻) or modelled (OAN⁻) using the Oliver and Hruška conceptual models, which are based on the density of weak acid functional sites (SD) present in DNOM. The charge density (CD) is governed by SD as well as protonation and complexation reactions on the functional groups. These models use SD as a key parameter to empirically fit the model to Org.⁻. Utilizing extensive water chemistry datasets, this study shows that spatial and temporal differences in SD and CD are influenced by variations in the humic-to-fulvic ratio of DNOM, organic aluminum (Al) complexation, and the mole fraction of CD to SD, which is governed by acidity. The median SD values obtained for 44 long-term monitored acid-sensitive lakes were 11.1 and 13.9 µEq/mg C for the Oliver and Hruška models, respectively. Over 34 years of monitoring, the CD increased by 70%, likely due to rising pH and declining Al complexation with DNOM. Present-day median SD values for the Oliver and Hruška models in 16 low-order streams are 13.8 and 15.8 µEq/mg C, respectively, and 10.8 and 12.5 µEq/mg C, respectively, in 10 high-order rivers.

Rising organic charge in northern freshwaters is attributed to increasing levels of Dissolved Natural Organic Matter (DNOM) and shifts in water chemistry. Organic charge concentration may be determined by charge balance calculations (Org.-) or modelled (OAN-) using the conceptual Oliver and Hruška models based on the density of weak acid functional sites (SD) present in DNOM. The charge density (CD) is governed by the SD and protonation and complexation reactions on the functional groups. The models use this SD as a key parameter in empirically fitting the model to Org.-. Utilizing extensive datasets of water chemistry, this study shows that spatial and temporal differences in SD and CD are influenced by variations in the humic-to-fulvic ratio of DNOM and the organic aluminum (Al) complexation, as well as the mole fraction of CD to SD governed by the acidity. Site median SD obtained for 44 long-term monitored acid-sensitive lakes was 11.1 and 13.9 µEq/mg C for the Oliver and the Hruška model, respectively. During the 34 years of monitoring the CD increased by 70%, likely due to rising pH and declining Al complexation with DNOM. Present-day median SD values for the Oliver and Hruška models in 16 low-order streams are 13.8 and 15.8 µEq/mg C, respectively, and 10.8 and 12.5 µEq/mg C, respectively, in 10 high order rivers.