Wildfire smoke impacts lake ecosystems.

Wildfire activity is increasing globally. The resulting smoke plumes can travel hundreds to thousands of kilometers, reflecting or scattering sunlight and depositing particles within ecosystems. Several key physical, chemical, and biological processes in lakes are controlled by factors affected by smoke. The spatial and temporal scales of lake exposure to smoke are extensive and under-recognized. We introduce the concept of the lake smoke-day, or the number of days any given lake is exposed to smoke in any given fire season, and quantify the total lake smoke-day exposure in North America from 2019 to 2021. Because smoke can be transported at continental to intercontinental scales, even regions that may not typically experience direct burning of landscapes by wildfire are at risk of smoke exposure. We found that 99.3% of North America was covered by smoke, affecting a total of 1,333,687 lakes ≥10 ha. An incredible 98.9% of lakes experienced at least 10 smoke-days a year, with 89.6% of lakes receiving over 30 lake smoke-days, and lakes in some regions experiencing up to 4 months of cumulative smoke-days. Herein we review the mechanisms through which smoke and ash can affect lakes by altering the amount and spectral composition of incoming solar radiation and depositing carbon, nutrients, or toxic compounds that could alter chemical conditions and impact biota. We develop a conceptual framework that synthesizes known and theoretical impacts of smoke on lakes to guide future research. Finally, we identify emerging research priorities that can help us better understand how lakes will be affected by smoke as wildfire activity increases due to climate change and other anthropogenic activities.

Towards critical white ice conditions in lakes under global warming.

The quality of lake ice is of uppermost importance for ice safety and under-ice ecology, but its temporal and spatial variability is largely unknown. Here we conducted a coordinated lake ice quality sampling campaign across the Northern Hemisphere during one of the warmest winters since 1880 and show that lake ice during 2020/2021 commonly consisted of unstable white ice, at times contributing up to 100% to the total ice thickness. We observed that white ice increased over the winter season, becoming thickest and constituting the largest proportion of the ice layer towards the end of the ice cover season when fatal winter drownings occur most often and light limits the growth and reproduction of primary producers. We attribute the dominance of white ice before ice-off to air temperatures varying around the freezing point, a condition which occurs more frequently during warmer winters. Thus, under continued global warming, the prevalence of white ice is likely to substantially increase during the critical period before ice-off, for which we adjusted commonly used equations for human ice safety and light transmittance through ice.

Current water quality guidelines across North America and Europe do not protect lakes from salinization.

Human-induced salinization caused by the use of road deicing salts, agricultural practices, mining operations, and climate change is a major threat to the biodiversity and functioning of freshwater ecosystems. Yet, it is unclear if freshwater ecosystems are protected from salinization by current water quality guidelines. Leveraging an experimental network of land-based and in-lake mesocosms across North America and Europe, we tested how salinization-indicated as elevated chloride (Cl-) concentration-will affect lake food webs and if two of the lowest Cl- thresholds found globally are sufficient to protect these food webs. Our results indicated that salinization will cause substantial zooplankton mortality at the lowest Cl- thresholds established in Canada (120 mg Cl-/L) and the United States (230 mg Cl-/L) and throughout Europe where Cl- thresholds are generally higher. For instance, at 73% of our study sites, Cl- concentrations that caused a ≥50% reduction in cladoceran abundance were at or below Cl- thresholds in Canada, in the United States, and throughout Europe. Similar trends occurred for copepod and rotifer zooplankton. The loss of zooplankton triggered a cascading effect causing an increase in phytoplankton biomass at 47% of study sites. Such changes in lake food webs could alter nutrient cycling and water clarity and trigger declines in fish production. Current Cl- thresholds across North America and Europe clearly do not adequately protect lake food webs. Water quality guidelines should be developed where they do not exist, and there is an urgent need to reassess existing guidelines to protect lake ecosystems from human-induced salinization.

Forest defoliator outbreaks alter nutrient cycling in northern waters.

Insect defoliators alter biogeochemical cycles from land into receiving waters by consuming terrestrial biomass and releasing biolabile frass. Here, we related insect outbreaks to water chemistry across 12 boreal lake catchments over 32-years. We report, on average, 27% lower dissolved organic carbon (DOC) and 112% higher dissolved inorganic nitrogen (DIN) concentrations in lake waters when defoliators covered entire catchments and reduced leaf area. DOC reductions reached 32% when deciduous stands dominated. Within-year changes in DOC from insect outbreaks exceeded 86% of between-year trends across a larger dataset of 266 boreal and north temperate lakes from 1990 to 2016. Similarly, within-year increases in DIN from insect outbreaks exceeded local, between-year changes in DIN by 12-times, on average. As insect defoliator outbreaks occur at least every 5 years across a wider 439,661 km2 boreal ecozone of Ontario, we suggest they are an underappreciated driver of biogeochemical cycles in forest catchments of this region.

Earlier winter/spring runoff and snowmelt during warmer winters lead to lower summer chlorophyll-a in north temperate lakes.

Winter conditions, such as ice cover and snow accumulation, are changing rapidly at northern latitudes and can have important implications for lake processes. For example, snowmelt in the watershed-a defining feature of lake hydrology because it delivers a large portion of annual nutrient inputs-is becoming earlier. Consequently, earlier and a shorter duration of snowmelt are expected to affect annual phytoplankton biomass. To test this hypothesis, we developed an index of runoff timing based on the date when 50% of cumulative runoff between January 1 and May 31 had occurred. The runoff index was computed using stream discharge for inflows, outflows, or for flows from nearby streams for 41 lakes in Europe and North America. The runoff index was then compared with summer chlorophyll-a (Chl-a) concentration (a proxy for phytoplankton biomass) across 5-53 years for each lake. Earlier runoff generally corresponded to lower summer Chl-a. Furthermore, years with earlier runoff also had lower winter/spring runoff magnitude, more protracted runoff, and earlier ice-out. We examined several lake characteristics that may regulate the strength of the relationship between runoff timing and summer Chl-a concentrations; however, our tested covariates had little effect on the relationship. Date of ice-out was not clearly related to summer Chl-a concentrations. Our results indicate that ongoing changes in winter conditions may have important consequences for summer phytoplankton biomass and production.

Widespread deoxygenation of temperate lakes.

The concentration of dissolved oxygen in aquatic systems helps to regulate biodiversity1,2, nutrient biogeochemistry3, greenhouse gas emissions4, and the quality of drinking water5. The long-term declines in dissolved oxygen concentrations in coastal and ocean waters have been linked to climate warming and human activity6,7, but little is known about the changes in dissolved oxygen concentrations in lakes. Although the solubility of dissolved oxygen decreases with increasing water temperatures, long-term lake trajectories are difficult to predict. Oxygen losses in warming lakes may be amplified by enhanced decomposition and stronger thermal stratification8,9 or oxygen may increase as a result of enhanced primary production10. Here we analyse a combined total of 45,148 dissolved oxygen and temperature profiles and calculate trends for 393 temperate lakes that span 1941 to 2017. We find that a decline in dissolved oxygen is widespread in surface and deep-water habitats. The decline in surface waters is primarily associated with reduced solubility under warmer water temperatures, although dissolved oxygen in surface waters increased in a subset of highly productive warming lakes, probably owing to increasing production of phytoplankton. By contrast, the decline in deep waters is associated with stronger thermal stratification and loss of water clarity, but not with changes in gas solubility. Our results suggest that climate change and declining water clarity have altered the physical and chemical environment of lakes. Declines in dissolved oxygen in freshwater are 2.75 to 9.3 times greater than observed in the world's oceans6,7 and could threaten essential lake ecosystem services2,3,5,11.

Deeper waters are changing less consistently than surface waters in a global analysis of 102 lakes.

Globally, lake surface water temperatures have warmed rapidly relative to air temperatures, but changes in deepwater temperatures and vertical thermal structure are still largely unknown. We have compiled the most comprehensive data set to date of long-term (1970-2009) summertime vertical temperature profiles in lakes across the world to examine trends and drivers of whole-lake vertical thermal structure. We found significant increases in surface water temperatures across lakes at an average rate of + 0.37 °C decade-1, comparable to changes reported previously for other lakes, and similarly consistent trends of increasing water column stability (+ 0.08 kg m-3 decade-1). In contrast, however, deepwater temperature trends showed little change on average (+ 0.06 °C decade-1), but had high variability across lakes, with trends in individual lakes ranging from - 0.68 °C decade-1 to + 0.65 °C decade-1. The variability in deepwater temperature trends was not explained by trends in either surface water temperatures or thermal stability within lakes, and only 8.4% was explained by lake thermal region or local lake characteristics in a random forest analysis. These findings suggest that external drivers beyond our tested lake characteristics are important in explaining long-term trends in thermal structure, such as local to regional climate patterns or additional external anthropogenic influences.

Road Salt Impacts Freshwater Zooplankton at Concentrations below Current Water Quality Guidelines.

Widespread use of NaCl for road deicing has caused increased chloride concentrations in lakes near urban centers and areas of high road density. Chloride can be toxic, and water quality guidelines have been created to regulate it and protect aquatic life. However, these guidelines may not adequately protect organisms in low-nutrient, soft water lakes such as those underlain by the Precambrian Shield. We tested this hypothesis by conducting laboratory experiments on six Daphnia species using a soft water culture medium. We also examined temporal changes in cladoceran assemblages in the sediments of two small lakes on the Canadian Shield: one near a highway and the other >3 km from roads where salt is applied in the winter. Our results showed that Daphnia were sensitive to low chloride concentrations with decreased reproduction and increased mortality occurring between 5 and 40 mg Cl-/L. Analysis of cladoceran remains in lake sediments revealed changes in assemblage composition that coincided with the initial application of road salt in this region. In contrast, there were no changes detected in the remote lake. We found that 22.7% of recreational lakes in Ontario have chloride concentrations between 5 and 40 mg/L suggesting that cladoceran zooplankton in these lakes may already be experiencing negative effects of chloride.

Storm impacts on phytoplankton community dynamics in lakes.

In many regions across the globe, extreme weather events such as storms have increased in frequency, intensity, and duration due to climate change. Ecological theory predicts that such extreme events should have large impacts on ecosystem structure and function. High winds and precipitation associated with storms can affect lakes via short-term runoff events from watersheds and physical mixing of the water column. In addition, lakes connected to rivers and streams will also experience flushing due to high flow rates. Although we have a well-developed understanding of how wind and precipitation events can alter lake physical processes and some aspects of biogeochemical cycling, our mechanistic understanding of the emergent responses of phytoplankton communities is poor. Here we provide a comprehensive synthesis that identifies how storms interact with lake and watershed attributes and their antecedent conditions to generate changes in lake physical and chemical environments. Such changes can restructure phytoplankton communities and their dynamics, as well as result in altered ecological function (e.g., carbon, nutrient and energy cycling) in the short- and long-term. We summarize the current understanding of storm-induced phytoplankton dynamics, identify knowledge gaps with a systematic review of the literature, and suggest future research directions across a gradient of lake types and environmental conditions.

Stochastic disturbance regimes alter patterns of ecosystem variability and recovery.

Altered ecosystem variability is an important ecological response to disturbance yet understanding of how various attributes of disturbance regimes affect ecosystem variability is limited. To improve the framework for understanding the disturbance regime attributes that affect ecosystem variability, we examine how the introduction of stochasticity to disturbance parameters (frequency, severity and extent) alters simulated recovery when compared to deterministic outcomes from a spatially explicit simulation model. We also examine the agreement between results from empirical studies and deterministic and stochastic configurations of the model. We find that stochasticity in disturbance frequency and spatial extent leads to the greatest increase in the variance of simulated dynamics, although stochastic severity also contributes to departures from the deterministic case. The incorporation of stochasticity in disturbance attributes improves agreement between empirical and simulated responses, with 71% of empirical responses correctly classified by stochastic configurations of the model as compared to 47% using the purely deterministic model. By comparison, only 2% of empirical responses were correctly classified by the deterministic model and misclassified by stochastic configurations of the model. These results indicate that stochasticity in the attributes of a disturbance regime alters the patterns and classification of ecosystem variability, suggesting altered recovery dynamics. Incorporating stochastic disturbance processes into models may thus be critical for anticipating the ecological resilience of ecosystems.

Terrestrial loads of dissolved organic matter drive inter-annual carbon flux in subtropical lakes during times of drought.

Lentic ecosystems are important agents of local and global carbon cycling, but their contribution varies along gradients of dissolved organic matter (DOM) and productivity. We investigated how contrasting summer and autumn precipitation can shape annual and inter-annual variation in ecosystem carbon (C) flux (gross primary production (GPP), ecosystem respiration (ER), and CO2 efflux) in two subtropical lakes differing substantially in trophic state and water color. Instrumented buoys recorded time series of free-water DO, terrestrial DOM (tDOM), chlorophyll a, water temperature profiles, and meteorological measurements over five years (2009-2011 and 2014-2015). Reduced precipitation caused immediate and prolonged effects on C flux in both lakes. During the drought year (2014) GPP and ER declined by 60 to 80% and both lakes were either CO2 sinks or neutral. In the subsequent wet year (2015), GPP and ER increased by 40 to 110%, and both lakes shifted to strong net CO2 emitters. Higher ecosystem R resulted from larger GPP while higher tDOM contributed to a dramatic increase in dissolved inorganic carbon, which intensified CO2 emission in both lakes. C flux was more responsive in the clear mesotrophic lake, declining by approximately 40% in the cumulative GPP and ER, and increasing by >400% in CO2 efflux whereas changes in the oligotrophic colored lake were more modest (approximately 30% and 300% for metabolic declines and efflux increases, respectively). Temporal variation and magnitude of C flux were governed by tDOM-mediated changes in epilimnetic nutrient levels and hypolimnetic light availability. This study demonstrated terrestrial loads of DOM strongly influence the inter-annual response and sensitivity of ecosystem C flux to variation in inter-annual precipitation. Our findings have important implications for predicting the trend, magnitude, duration, and sensitivity of the response of C flux in subtropical lakes/reservoirs to future changes in precipitation patterns under altered climatic conditions.

The NSERC Canadian Lake Pulse Network: A national assessment of lake health providing science for water management in a changing climate.

The distribution and quality of water resources vary dramatically across Canada, and human impacts such as land-use and climate changes are exacerbating uncertainties in water supply and security. At the national level, Canada has no enforceable standards for safe drinking water and no comprehensive water-monitoring program to provide detailed, timely reporting on the state of water resources. To provide Canada's first national assessment of lake health, the NSERC Canadian Lake Pulse Network was launched in 2016 as an academic-government research partnership. LakePulse uses traditional approaches for limnological monitoring as well as state-of-the-art methods in the fields of genomics, emerging contaminants, greenhouse gases, invasive pathogens, paleolimnology, spatial modelling, statistical analysis, and remote sensing. A coordinated sampling program of about 680 lakes together with historical archives and a geomatics analysis of over 80,000 lake watersheds are used to examine the extent to which lakes are being altered now and in the future, and how this impacts aquatic ecosystem services of societal importance. Herein we review the network context, objectives and methods.

Widespread diminishing anthropogenic effects on calcium in freshwaters.

Calcium (Ca) is an essential element for almost all living organisms. Here, we examined global variation and controls of freshwater Ca concentrations, using 440 599 water samples from 43 184 inland water sites in 57 countries. We found that the global median Ca concentration was 4.0 mg L-1 with 20.7% of the water samples showing Ca concentrations ≤ 1.5 mg L-1, a threshold considered critical for the survival of many Ca-demanding organisms. Spatially, freshwater Ca concentrations were strongly and proportionally linked to carbonate alkalinity, with the highest Ca and carbonate alkalinity in waters with a pH around 8.0 and decreasing in concentrations towards lower pH. However, on a temporal scale, by analyzing decadal trends in >200 water bodies since the 1980s, we observed a frequent decoupling between carbonate alkalinity and Ca concentrations, which we attributed mainly to the influence of anthropogenic acid deposition. As acid deposition has been ameliorated, in many freshwaters carbonate alkalinity concentrations have increased or remained constant, while Ca concentrations have rapidly declined towards or even below pre-industrial conditions as a consequence of recovery from anthropogenic acidification. Thus, a paradoxical outcome of the successful remediation of acid deposition is a globally widespread freshwater Ca concentration decline towards critically low levels for many aquatic organisms.

A comparative evaluation of five common pairwise tests of species association.

Researchers have long viewed patterns of species association as key to understanding the processes that structure communities. Community-level tests of species association have received the most attention; however, pairwise species associations may offer greater opportunity for linking patterns to specific mechanisms. Although several tests of pairwise association have been developed, there remain gaps in our understanding of their performance. Consequently, it is unclear whether these methods reliably detect patterns of association, or if any one method is superior. We maximized association patterns for single species pairs in synthetic community matrices and examined how accurately five pairwise association tests found that pair, while not finding others (i.e., type I and II error rates). All tests are more likely to miss patterns of association than to falsely detect them. When we maximized association for a species pair that included one or more rare or common species, tests were frequently unable to identify that pair as significantly associated. Consequently, these tests are best suited for identifying significant associations between pairs of species that occur in an intermediate number of samples; for such pairs, three of the five tests considered here detected 100% of the pairs for which we maximized associations.

Citizen science shows systematic changes in the temperature difference between air and inland waters with global warming.

Citizen science projects have a long history in ecological studies. The research usefulness of such projects is dependent on applying simple and standardized methods. Here, we conducted a citizen science project that involved more than 3500 Swedish high school students to examine the temperature difference between surface water and the overlying air (Tw-Ta) as a proxy for sensible heat flux (QH). If QH is directed upward, corresponding to positive Tw-Ta, it can enhance CO2 and CH4 emissions from inland waters, thereby contributing to increased greenhouse gas concentrations in the atmosphere. The students found mostly negative Tw-Ta across small ponds, lakes, streams/rivers and the sea shore (i.e. downward QH), with Tw-Ta becoming increasingly negative with increasing Ta. Further examination of Tw-Ta using high-frequency temperature data from inland waters across the globe confirmed that Tw-Ta is linearly related to Ta. Using the longest available high-frequency temperature time series from Lake Erken, Sweden, we found a rapid increase in the occasions of negative Tw-Ta with increasing annual mean Ta since 1989. From these results, we can expect that ongoing and projected global warming will result in increasingly negative Tw-Ta, thereby reducing CO2 and CH4 transfer velocities from inland waters into the atmosphere.

Diel Surface Temperature Range Scales with Lake Size.

Ecological and biogeochemical processes in lakes are strongly dependent upon water temperature. Long-term surface warming of many lakes is unequivocal, but little is known about the comparative magnitude of temperature variation at diel timescales, due to a lack of appropriately resolved data. Here we quantify the pattern and magnitude of diel temperature variability of surface waters using high-frequency data from 100 lakes. We show that the near-surface diel temperature range can be substantial in summer relative to long-term change and, for lakes smaller than 3 km2, increases sharply and predictably with decreasing lake area. Most small lakes included in this study experience average summer diel ranges in their near-surface temperatures of between 4 and 7°C. Large diel temperature fluctuations in the majority of lakes undoubtedly influence their structure, function and role in biogeochemical cycles, but the full implications remain largely unexplored.

Determining the probability of cyanobacterial blooms: the application of Bayesian networks in multiple lake systems.

A Bayesian network model was developed to assess the combined influence of nutrient conditions and climate on the occurrence of cyanobacterial blooms within lakes of diverse hydrology and nutrient supply. Physicochemical, biological, and meteorological observations were collated from 20 lakes located at different latitudes and characterized by a range of sizes and trophic states. Using these data, we built a Bayesian network to (1) analyze the sensitivity of cyanobacterial bloom development to different environmental factors and (2) determine the probability that cyanobacterial blooms would occur. Blooms were classified in three categories of hazard (low, moderate, and high) based on cell abundances. The most important factors determining cyanobacterial bloom occurrence were water temperature, nutrient availability, and the ratio of mixing depth to euphotic depth. The probability of cyanobacterial blooms was evaluated under different combinations of total phosphorus and water temperature. The Bayesian network was then applied to quantify the probability of blooms under a future climate warming scenario. The probability of the "high hazardous" category of cyanobacterial blooms increased 5% in response to either an increase in water temperature of 0.8°C (initial water temperature above 24°C) or an increase in total phosphorus from 0.01 mg/L to 0.02 mg/L. Mesotrophic lakes were particularly vulnerable to warming. Reducing nutrient concentrations counteracts the increased cyanobacterial risk associated with higher temperatures.

The jellification of north temperate lakes.

Calcium (Ca) concentrations are decreasing in softwater lakes across eastern North America and western Europe. Using long-term contemporary and palaeo-environmental field data, we show that this is precipitating a dramatic change in Canadian lakes: the replacement of previously dominant pelagic herbivores (Ca-rich Daphnia species) by Holopedium glacialis, a jelly-clad, Ca-poor competitor. In some lakes, this transformation is being facilitated by increases in macro-invertebrate predation, both from native (Chaoborus spp.) and introduced (Bythotrephes longimanus) zooplanktivores, to which Holopedium, with its jelly coat, is relatively invulnerable. Greater representation by Holopedium within cladoceran zooplankton communities will reduce nutrient transfer through food webs, given their lower phosphorus content relative to daphniids, and greater absolute abundances may pose long-term problems to water users. The dominance of jelly-clad zooplankton will likely persist while lakewater Ca levels remain low.

Synchronous dynamics of zooplankton competitors prevail in temperate lake ecosystems.

Although competing species are expected to exhibit compensatory dynamics (negative temporal covariation), empirical work has demonstrated that competitive communities often exhibit synchronous dynamics (positive temporal covariation). This has led to the suggestion that environmental forcing dominates species dynamics; however, synchronous and compensatory dynamics may appear at different length scales and/or at different times, making it challenging to identify their relative importance. We compiled 58 long-term datasets of zooplankton abundance in north-temperate and sub-tropical lakes and used wavelet analysis to quantify general patterns in the times and scales at which synchronous/compensatory dynamics dominated zooplankton communities in different regions and across the entire dataset. Synchronous dynamics were far more prevalent at all scales and times and were ubiquitous at the annual scale. Although we found compensatory dynamics in approximately 14% of all combinations of time period/scale/lake, there were no consistent scales or time periods during which compensatory dynamics were apparent across different regions. Our results suggest that the processes driving compensatory dynamics may be local in their extent, while those generating synchronous dynamics operate at much larger scales. This highlights an important gap in our understanding of the interaction between environmental and biotic forces that structure communities.