The enduring world forest carbon sink.

The uptake of carbon dioxide (CO2) by terrestrial ecosystems is critical for moderating climate change1. To provide a ground-based long-term assessment of the contribution of forests to terrestrial CO2 uptake, we synthesized in situ forest data from boreal, temperate and tropical biomes spanning three decades. We found that the carbon sink in global forests was steady, at 3.6 ± 0.4 Pg C yr-1 in the 1990s and 2000s, and 3.5 ± 0.4 Pg C yr-1 in the 2010s. Despite this global stability, our analysis revealed some major biome-level changes. Carbon sinks have increased in temperate (+30 ± 5%) and tropical regrowth (+29 ± 8%) forests owing to increases in forest area, but they decreased in boreal (-36 ± 6%) and tropical intact (-31 ± 7%) forests, as a result of intensified disturbances and losses in intact forest area, respectively. Mass-balance studies indicate that the global land carbon sink has increased2, implying an increase in the non-forest-land carbon sink. The global forest sink is equivalent to almost half of fossil-fuel emissions (7.8 ± 0.4 Pg C yr-1 in 1990-2019). However, two-thirds of the benefit from the sink has been negated by tropical deforestation (2.2 ± 0.5 Pg C yr-1 in 1990-2019). Although the global forest sink has endured undiminished for three decades, despite regional variations, it could be weakened by ageing forests, continuing deforestation and further intensification of disturbance regimes1. To protect the carbon sink, land management policies are needed to limit deforestation, promote forest restoration and improve timber-harvesting practices1,3.

Even modest climate change may lead to major transitions in boreal forests.

The sensitivity of forests to near-term warming and associated precipitation shifts remains uncertain1-9. Herein, using a 5-year open-air experiment in southern boreal forest, we show divergent responses to modest climate alteration among juveniles of nine co-occurring North American tree species. Warming alone (+1.6 °C or +3.1 °C above ambient temperature) or combined with reduced rainfall increased the juvenile mortality of all species, especially boreal conifers. Species differed in growth responses to warming, ranging from enhanced growth in Acer rubrum and Acer saccharum to severe growth reductions in Abies balsamea, Picea glauca and Pinus strobus. Moreover, treatment-induced changes in both photosynthesis and growth help explain treatment-driven changes in survival. Treatments in which species experienced conditions warmer or drier than at their range margins resulted in the most adverse impacts on growth and survival. Species abundant in southern boreal forests had the largest reductions in growth and survival due to climate manipulations. By contrast, temperate species that experienced little mortality and substantial growth enhancement in response to warming are rare throughout southern boreal forest and unlikely to rapidly expand their density and distribution. Therefore, projected climate change will probably cause regeneration failure of currently dominant southern boreal species and, coupled with their slow replacement by temperate species, lead to tree regeneration shortfalls with potential adverse impacts on the health, diversity and ecosystem services of regional forests.

Sufficient conditions for rapid range expansion of a boreal conifer.

Unprecedented modern rates of warming are expected to advance boreal forest into Arctic tundra1, thereby reducing albedo2-4, altering carbon cycling4 and further changing climate1-4, yet the patterns and processes of this biome shift remain unclear5. Climate warming, required for previous boreal advances6-17, is not sufficient by itself for modern range expansion of conifers forming forest-tundra ecotones5,12-15,17-20. No high-latitude population of conifers, the dominant North American Arctic treeline taxon, has previously been documented5 advancing at rates following the last glacial maximum (LGM)6-8. Here we describe a population of white spruce (Picea glauca) advancing at post-LGM rates7 across an Arctic basin distant from established treelines and provide evidence of mechanisms sustaining the advance. The population doubles each decade, with exponential radial growth in the main stems of individual trees correlating positively with July air temperature. Lateral branches in adults and terminal leaders in large juveniles grow almost twice as fast as those at established treelines. We conclude that surpassing temperature thresholds1,6-17, together with winter winds facilitating long-distance dispersal, deeper snowpack and increased soil nutrient availability promoting recruitment and growth, provides sufficient conditions for boreal forest advance. These observations enable forecast modelling with important insights into the environmental conditions converting tundra into forest.

Emerging signals of declining forest resilience under climate change.

Forest ecosystems depend on their capacity to withstand and recover from natural and anthropogenic perturbations (that is, their resilience)1. Experimental evidence of sudden increases in tree mortality is raising concerns about variation in forest resilience2, yet little is known about how it is evolving in response to climate change. Here we integrate satellite-based vegetation indices with machine learning to show how forest resilience, quantified in terms of critical slowing down indicators3-5, has changed during the period 2000-2020. We show that tropical, arid and temperate forests are experiencing a significant decline in resilience, probably related to increased water limitations and climate variability. By contrast, boreal forests show divergent local patterns with an average increasing trend in resilience, probably benefiting from warming and CO2 fertilization, which may outweigh the adverse effects of climate change. These patterns emerge consistently in both managed and intact forests, corroborating the existence of common large-scale climate drivers. Reductions in resilience are statistically linked to abrupt declines in forest primary productivity, occurring in response to slow drifting towards a critical resilience threshold. Approximately 23% of intact undisturbed forests, corresponding to 3.32 Pg C of gross primary productivity, have already reached a critical threshold and are experiencing a further degradation in resilience. Together, these signals reveal a widespread decline in the capacity of forests to withstand perturbation that should be accounted for in the design of land-based mitigation and adaptation plans.

Overwintering fires in boreal forests.

Forest fires are usually viewed within the context of a single fire season, in which weather conditions and fuel supply can combine to create conditions favourable for fire ignition-usually by lightning or human activity-and spread1-3. But some fires exhibit 'overwintering' behaviour, in which they smoulder through the non-fire season and flare up in the subsequent spring4,5. In boreal (northern) forests, deep organic soils favourable for smouldering6, along with accelerated climate warming7, may present unusually favourable conditions for overwintering. However, the extent of overwintering in boreal forests and the underlying factors influencing this behaviour remain unclear. Here we show that overwintering fires in boreal forests are associated with hot summers generating large fire years and deep burning into organic soils, conditions that have become more frequent in our study areas in recent decades. Our results are based on an algorithm with which we detect overwintering fires in Alaska, USA, and the Northwest Territories, Canada, using field and remote sensing datasets. Between 2002 and 2018, overwintering fires were responsible for 0.8 per cent of the total burned area; however, in one year this amounted to 38 per cent. The spatiotemporal predictability of overwintering fires could be used by fire management agencies to facilitate early detection, which may result in reduced carbon emissions and firefighting costs.

Exchangeable manganese regulates carbon storage in the humus layer of the boreal forest.

The huge carbon stock in humus layers of the boreal forest plays a critical role in the global carbon cycle. However, there remains uncertainty about the factors that regulate below-ground carbon sequestration in this region. Notably, based on evidence from two independent but complementary methods, we identified that exchangeable manganese is a critical factor regulating carbon accumulation in boreal forests across both regional scales and the entire boreal latitudinal range. Moreover, in a novel fertilization experiment, manganese addition reduced soil carbon stocks, but only after 4 y of additions. Our results highlight an underappreciated mechanism influencing the humus carbon pool of boreal forests.

General reversal of N-decomposition relationship during long-term decomposition in boreal and temperate forests.

Decomposition of dead organic matter is fundamental to carbon (C) and nutrient cycling in terrestrial ecosystems, influencing C fluxes from the biosphere to the atmosphere. Theory predicts and evidence strongly supports that the availability of nitrogen (N) limits litter decomposition. Positive relationships between substrate N concentrations and decomposition have been embedded into ecosystem models. This decomposition paradigm, however, relies on data mostly from short-term studies analyzing controls on early-stage decomposition. We present evidence from three independent long-term decomposition investigations demonstrating that the positive N-decomposition relationship is reversed and becomes negative during later stages of decomposition. First, in a 10-y decomposition experiment across 62 woody species in a temperate forest, leaf litter with higher N concentrations exhibited faster initial decomposition rates but ended up a larger recalcitrant fraction decomposing at a near-zero rate. Second, in a 5-y N-enrichment experiment of two tree species, leaves with experimentally enriched N concentrations had faster decomposition initial rates but ultimately accumulated large slowly decomposing fractions. Measures of amino sugars on harvested litter in two experiments indicated that greater accumulation of microbial residues in N-rich substrates likely contributed to larger slowly decomposing fractions. Finally, a database of 437 measurements from 120 species in 45 boreal and temperate forest sites confirmed that higher N concentrations were associated with a larger slowly decomposing fraction. These results challenge the current treatment of interactions between N and decomposition in many ecosystems and Earth system models and suggest that even the best-supported short-term controls of biogeochemical processes might not predict long-term controls.

Warming-induced tree growth may help offset increasing disturbance across the Canadian boreal forest.

Large projected increases in forest disturbance pose a major threat to future wood fiber supply and carbon sequestration in the cold-limited, Canadian boreal forest ecosystem. Given the large sensitivity of tree growth to temperature, warming-induced increases in forest productivity have the potential to reduce these threats, but research efforts to date have yielded contradictory results attributed to limited data availability, methodological biases, and regional variability in forest dynamics. Here, we apply a machine learning algorithm to an unprecedented network of over 1 million tree growth records (1958 to 2018) from 20,089 permanent sample plots distributed across both Canada and the United States, spanning a 16.5 °C climatic gradient. Fitted models were then used to project the near-term (2050 s time period) growth of the six most abundant tree species in the Canadian boreal forest. Our results reveal a large, positive effect of increasing thermal energy on tree growth for most of the target species, leading to 20.5 to 22.7% projected gains in growth with climate change under RCP 4.5 and 8.5. The magnitude of these gains, which peak in the colder and wetter regions of the boreal forest, suggests that warming-induced growth increases should no longer be considered marginal but may in fact significantly offset some of the negative impacts of projected increases in drought and wildfire on wood supply and carbon sequestration and have major implications on ecological forecasts and the global economy.

Increasing wildfires threaten historic carbon sink of boreal forest soils.

Boreal forest fires emit large amounts of carbon into the atmosphere primarily through the combustion of soil organic matter1-3. During each fire, a portion of this soil beneath the burned layer can escape combustion, leading to a net accumulation of carbon in forests over multiple fire events4. Climate warming and drying has led to more severe and frequent forest fires5-7, which threaten to shift the carbon balance of the boreal ecosystem from net accumulation to net loss1, resulting in a positive climate feedback8. This feedback will occur if organic-soil carbon that escaped burning in previous fires, termed 'legacy carbon', combusts. Here we use soil radiocarbon dating to quantitatively assess legacy carbon loss in the 2014 wildfires in the Northwest Territories of Canada2. We found no evidence for the combustion of legacy carbon in forests that were older than the historic fire-return interval of northwestern boreal forests9. In forests that were in dry landscapes and less than 60 years old at the time of the fire, legacy carbon that had escaped burning in the previous fire cycle was combusted. We estimate that 0.34 million hectares of young forests (<60 years) that burned in the 2014 fires could have experienced legacy carbon combustion. This implies a shift to a domain of carbon cycling in which these forests become a net source-instead of a sink-of carbon to the atmosphere over consecutive fires. As boreal wildfires continue to increase in size, frequency and intensity7, the area of young forests that experience legacy carbon combustion will probably increase and have a key role in shifting the boreal carbon balance.

Metatranscriptomics captures dynamic shifts in mycorrhizal coordination in boreal forests.

Carbon storage and cycling in boreal forests-the largest terrestrial carbon store-is moderated by complex interactions between trees and soil microorganisms. However, existing methods limit our ability to predict how changes in environmental conditions will alter these associations and the essential ecosystem services they provide. To address this, we developed a metatranscriptomic approach to analyze the impact of nutrient enrichment on Norway spruce fine roots and the community structure, function, and tree-microbe coordination of over 350 root-associated fungal species. In response to altered nutrient status, host trees redefined their relationship with the fungal community by reducing sugar efflux carriers and enhancing defense processes. This resulted in a profound restructuring of the fungal community and a collapse in functional coordination between the tree and the dominant Basidiomycete species, and an increase in functional coordination with versatile Ascomycete species. As such, there was a functional shift in community dominance from Basidiomycetes species, with important roles in enzymatically cycling recalcitrant carbon, to Ascomycete species that have melanized cell walls that are highly resistant to degradation. These changes were accompanied by prominent shifts in transcriptional coordination between over 60 predicted fungal effectors, with more than 5,000 Norway spruce transcripts, providing mechanistic insight into the complex molecular dialogue coordinating host trees and their fungal partners. The host-microbe dynamics captured by this study functionally inform how these complex and sensitive biological relationships may mediate the carbon storage potential of boreal soils under changing nutrient conditions.

Shedding light on the increased carbon uptake by a boreal forest under diffuse solar radiation across multiple scales.

Solar radiation is scattered by cloud cover, aerosols and other particles in the atmosphere, all of which are affected by global changes. Furthermore, the diffuse fraction of solar radiation is increased by more frequent forest fires and likewise would be if climate interventions such as stratospheric aerosol injection were adopted. Forest ecosystem studies predict that an increase in diffuse radiation would result in higher productivity, but ecophysiological data are required to identify the processes responsible within the forest canopy. In our study, the response of a boreal forest to direct, diffuse and heterogeneous solar radiation conditions was examined during the daytime in the growing season to determine how carbon uptake is affected by radiation conditions at different scales. A 10-year data set of ecosystem, shoot and forest floor vegetation carbon and water-flux data was examined. Ecosystem-level carbon assimilation was higher under diffuse radiation conditions in comparison with direct radiation conditions at equivalent total photosynthetically active radiation (PAR). This was driven by both an increase in shoot and forest floor vegetation photosynthetic rate. Most notably, ecosystem-scale productivity was strongly related to the absolute amount of diffuse PAR, since it integrates both changes in total PAR and diffuse fraction. This finding provides a gateway to explore the processes by which absolute diffuse PAR enhances productivity, and the long-term persistence of this effect under scenarios of higher global diffuse radiation.

Rapid migration of Mongolian oak into the southern Asian boreal forest.

The migration of trees induced by climatic warming has been observed at many alpine treelines and boreal-tundra ecotones, but the migration of temperate trees into southern boreal forest remains less well documented. We conducted a field investigation across an ecotone of temperate and boreal forests in northern Greater Khingan Mountains of northeast China. Our analysis demonstrates that Mongolian oak (Quercus mongolica), an important temperate tree species, has migrated rapidly into southern boreal forest in synchrony with significant climatic warming over the past century. The average rate of migration is estimated to be 12.0 ± 1.0 km decade-1 , being slightly slower than the movement of isotherms (14.7 ± 6.4 km decade-1 ). The migration rate of Mongolian oak is the highest observed among migratory temperate trees (average rate 4.0 ± 1.0 km decade-1 ) and significantly higher than the rates of tree migration at boreal-tundra ecotones (0.9 ± 0.4 km decade-1 ) and alpine treelines (0.004 ± 0.003 km decade-1 ). Compared with the coexisting dominant boreal tree species, Dahurian larch (Larix gmelinii), temperate Mongolian oak is observed to have significantly lower capacity for light acquisition, comparable water-use efficiency but stronger capacity to utilize nutrients especially the most limiting nutrient, nitrogen. In the context of climatic warming, and in addition to a high seed dispersal capacity and potential thermal niche differences, the advantage of nutrient utilization, reflected by foliar elementomes and stable nitrogen isotope ratios, is also likely a key mechanism for Mongolian oak to coexist with Dahurian larch and facilitate its migration toward boreal forest. These findings highlight a rapid deborealization of southern Asian boreal forest in response to climatic warming.

Pixel walking along the boreal forest-Arctic tundra ecotone: Large scale ground-truthing of satellite-derived greenness (NDVI).

In this Technical Advance, we describe a novel method to improve ecological interpretation of remotely sensed vegetation greenness measurements that involved sampling 24,395 Landsat pixels (30 m) across 639 km of Alaska's central Brooks Range. The method goes well beyond the spatial scale of traditional plot-based sampling and thereby more thoroughly relates ground-based observations to satellite measurements. Our example dataset illustrates that, along the boreal-Arctic boundary, vegetation with the greatest Landsat Normalized Difference Vegetation Index (NDVI) is taller than 1 m, woody, and deciduous; whereas vegetation with lower NDVI tends to be shorter, evergreen, or non-woody. The field methods and associated analyses advance efforts to inform satellite data with ground-based vegetation observations using field samples collected at spatial scales that closely match the resolution of remotely sensed imagery.

Hybridization mediated range expansion and climate change resilience in two keystone tree species of boreal forests.

Current global climate change is expected to affect biodiversity negatively at all scales leading to mass biodiversity loss. Many studies have shown that the distribution of allele frequencies across a species' range is often influenced by specific genetic loci associated with local environmental variables. This association reflects local adaptation and allele changes at those loci could thereby contribute to the evolutionary response to climate change. However, predicting how species will adapt to climate change from this type of data alone remains challenging. In the present study, we combined exome capture sequences and environmental niche reconstruction, to test multiple methods for assessing local adaptation and climate resilience in two widely distributed conifers, Norway spruce and Siberian spruce. Both species are keystone species of the boreal forest and share a vast hybrid zone. We show that local adaptation in conifers can be detected through allele frequency variation, population-level ecological preferences, and historical niche movement. Moreover, we integrated genetic and ecological information into genetic offset predictive models to show that hybridization plays a central role in expanding the niche breadth of the two conifer species and may help both species to cope better with future changing climates. This joint genetic and ecological analysis also identified spruce populations that are at risk under current climate change.

Wildfire-induced increases in photosynthesis in boreal forest ecosystems of North America.

Observations of the annual cycle of atmospheric CO2 in high northern latitudes provide evidence for an increase in terrestrial metabolism in Arctic tundra and boreal forest ecosystems. However, the mechanisms driving these changes are not yet fully understood. One proposed hypothesis is that ecological change from disturbance, such as wildfire, could increase the magnitude and change the phase of net ecosystem exchange through shifts in plant community composition. Yet, little quantitative work has evaluated this potential mechanism at a regional scale. Here we investigate how fire disturbance influences landscape-level patterns of photosynthesis across western boreal North America. We use Alaska and Canadian large fire databases to identify the perimeters of wildfires, a Landsat-derived land cover time series to characterize plant functional types (PFTs), and solar-induced fluorescence (SIF) from the Orbiting Carbon Observatory-2 (OCO-2) as a proxy for photosynthesis. We analyze these datasets to characterize post-fire changes in plant succession and photosynthetic activity using a space-for-time approach. We find that increases in herbaceous and sparse vegetation, shrub, and deciduous broadleaf forest PFTs during mid-succession yield enhancements in SIF by 8-40% during June and July for 2- to 59-year stands relative to pre-fire controls. From the analysis of post-fire land cover changes within individual ecoregions and modeling, we identify two mechanisms by which fires contribute to long-term trends in SIF. First, increases in annual burning are shifting the stand age distribution, leading to increases in the abundance of shrubs and deciduous broadleaf forests that have considerably higher SIF during early- and mid-summer. Second, fire appears to facilitate a long-term shift from evergreen conifer to broadleaf deciduous forest in the Boreal Plain ecoregion. These findings suggest that increasing fire can contribute substantially to positive trends in seasonal CO2 exchange without a close coupling to long-term increases in carbon storage.

The biological controls of soil carbon accumulation following wildfire and harvest in boreal forests: A review.

Boreal forests are frequently subjected to disturbances, including wildfire and clear-cutting. While these disturbances can cause soil carbon (C) losses, the long-term accumulation dynamics of soil C stocks during subsequent stand development is controlled by biological processes related to the balance of net primary production (NPP) and outputs via heterotrophic respiration and leaching, many of which remain poorly understood. We review the biological processes suggested to influence soil C accumulation in boreal forests. Our review indicates that median C accumulation rates following wildfire and clear-cutting are similar (0.15 and 0.20 Mg ha-1 year-1, respectively), however, variation between studies is extremely high. Further, while many individual studies show linear increases in soil C stocks through time after disturbance, there are indications that C stock recovery is fastest early to mid-succession (e.g. 15-80 years) and then slows as forests mature (e.g. >100 years). We indicate that the rapid build-up of soil C in younger stands appears not only driven by higher plant production, but also by a high rate of mycorrhizal hyphal production, and mycorrhizal suppression of saprotrophs. As stands mature, the balance between reductions in plant and mycorrhizal production, increasing plant litter recalcitrance, and ectomycorrhizal decomposers and saprotrophs have been highlighted as key controls on soil C accumulation rates. While some of these controls appear well understood (e.g. temporal patterns in NPP, changes in aboveground litter quality), many others remain research frontiers. Notably, very little data exists describing and comparing successional patterns of root production, mycorrhizal functional traits, mycorrhizal-saprotroph interactions, or C outputs via heterotrophic respiration and dissolved organic C following different disturbances. We argue that these less frequently described controls require attention, as they will be key not only for understanding ecosystem C balances, but also for representing these dynamics more accurately in soil organic C and Earth system models.

Eastern Canadian boreal forest soil and foliar chemistry show evidence of resilience to long-term nitrogen addition.

The boreal forest is one of the world's largest terrestrial biome and plays crucial roles in global biogeochemical cycles, such as carbon (C) sequestration in vegetation and soil. However, the impacts of decades of N deposition on N-limited ecosystems, like the eastern Canadian boreal forest, remain unclear. For 13 years, N deposition was simulated by periodically adding ammonium nitrate on soils of two boreal coniferous forests (i.e., balsam fir and black spruce) of eastern Canada, at low (LN) and high (HN) rates, corresponding to 3 and 10 times the ambient N deposition, respectively. We show that more than a decade of N addition had no strong effects on mineral soil C, N, P, and cation concentrations and on foliar total Ca, K, Mg, and Mn concentrations. In organic soil, C stock was not affected by N addition while N stock increased, and exchangeable Ca2+ and Mg2+ decreased at the balsam fir site under HN treatment. At both sites, LN treatment had nearly no impact on foliage and soil chemistry but foliar N and N:P significantly increased under HN treatment, potentially leading to foliar nutrient imbalance. Overall, our work indicates that, in the eastern Canadian boreal forest, soil and foliar nutrient concentrations and stocks are resilient to increasing N deposition potentially because, in the context of N limitation, extra N would be rapidly immobilized by soil micro-organisms and vegetation. These findings could improve modeling future boreal forest soil C stocks and biomass growth and could help in planning forest management strategies in eastern Canada.

Postfire Phosphorus Enrichment Mitigates Nitrogen Loss in Boreal Forests.

Net nitrogen mineralization (Nmin) and nitrification regulate soil N availability and loss after severe wildfires in boreal forests experiencing slow vegetation recovery. Yet, how microorganisms respond to postfire phosphorus (P) enrichment to alter soil N transformations remains unclear in N-limited boreal forests. Here, we investigated postfire N-P interactions using an intensive regional-scale sampling of 17 boreal forests in the Greater Khingan Mountains (Inner Mongolia-China), a laboratory P-addition incubation, and a continental-scale meta-analysis. We found that postfire soils had an increased risk of N loss by accelerated Nmin and nitrification along with low plant N demand, especially during the early vegetation recovery period. The postfire N/P imbalance created by P enrichment acts as a "N retention" strategy by inhibiting Nmin but not nitrification in boreal forests. This strategy is attributed to enhanced microbial N-use efficiency and N immobilization. Importantly, our meta-analysis found that there was a greater risk of N loss in boreal forest soils after fires than in other climatic zones, which was consistent with our results from the 17 soils in the Greater Khingan Mountains. These findings demonstrate that postfire N-P interactions play an essential role in mitigating N limitation and maintaining nutrient balance in boreal forests.

Increasing fire and the decline of fire adapted black spruce in the boreal forest.

Intensifying wildfire activity and climate change can drive rapid forest compositional shifts. In boreal North America, black spruce shapes forest flammability and depends on fire for regeneration. This relationship has helped black spruce maintain its dominance through much of the Holocene. However, with climate change and more frequent and severe fires, shifts away from black spruce dominance to broadleaf or pine species are emerging, with implications for ecosystem functions including carbon sequestration, water and energy fluxes, and wildlife habitat. Here, we predict that such reductions in black spruce after fire may already be widespread given current trends in climate and fire. To test this, we synthesize data from 1,538 field sites across boreal North America to evaluate compositional changes in tree species following 58 recent fires (1989 to 2014). While black spruce was resilient following most fires (62%), loss of resilience was common, and spruce regeneration failed completely in 18% of 1,140 black spruce sites. In contrast, postfire regeneration never failed in forests dominated by jack pine, which also possesses an aerial seed bank, or broad-leaved trees. More complete combustion of the soil organic layer, which often occurs in better-drained landscape positions and in dryer duff, promoted compositional changes throughout boreal North America. Forests in western North America, however, were more vulnerable to change due to greater long-term climate moisture deficits. While we find considerable remaining resilience in black spruce forests, predicted increases in climate moisture deficits and fire activity will erode this resilience, pushing the system toward a tipping point that has not been crossed in several thousand years.

Analysis of taiga and tundra lake browning trends from 2002 to 2021 using MODIS data.

Lakes in taiga and tundra regions may be silently undergoing changes due to global warming. One of those changes is browning in lake color. The browning interacts with the carbon cycle, ecosystem dynamics, and water quality in freshwater systems. However, spatiotemporal variabilities of browning in these regions have not been well documented. Using MODIS remote sensing reflectance at near ultraviolet wavelengths from 2002 to 2021 on the Google Earth Engine platform, we quantified long-term browning trends across 7616 lakes (larger than 10 km2) in taiga and tundra biomes. These lakes showed an overall decreased trend in browning (Theil-Sen Slope = 0.00015), with ∼36% of these lakes showing browning trends, and ∼1% of these lakes showing statistically significant (p-value <0.05) browning trends. The browning trends more likely occurred in small lakes in high latitude, low ground ice content regions, where air temperature increased and precipitation decreased. While temperature is projected to increase in response to climate change, our results provide one means to understand how biogeochemical cycles and ecological dynamics respond to climate change.