Altered growth conditions more than reforestation counteracted forest biomass carbon emissions 1990–2020
Trends in global biomass C stocksThe CRAFT model reliably reproduces the observed trends in primary and managed forest biomass C stocks (including both above-ground and belowground biomass) in 1990–2020 with a relative root mean square error (RMSE) of 0.57% between simulated and observed biomass C stocks by the FRA2 at the global level. These low divergences between stock estimates result, however, in global C emissions c. 2 times lower according to the CRAFT simulations than the estimates derived from the FRA (Supplementary Table 2). Still, the CRAFT simulations corroborated the FRA observations while adding information on annual estimations of forest C stocks, rather than the 5-years interval data provided by the FRA (Fig. 1a), and dynamic annual net C emissions (Fig. 1c, d) from managed and primary forests. The five sensitivity analyses carried out on the most uncertain model inputs and assumptions (see ‘Methods’ descriptions and Supplementary Figs. 5–10) confirmed the results presented in Fig.1: The largest deviation derived from the sensitivity analysis considering forest gross instead of net area changes results in a relative RMSE of 1.94% with global C emissions c.2.5 times higher than the FRA estimates (Supplementary Table 2). The simulations from the reference model assumptions yield the best RMSE and closest agreement with the C budgets derived from the FRA, indicating that they are the most optimal.Fig. 1: Global trends in total, primary, and managed forests.a Forest biomass C stocks (GtC); b cumulated change in forest area (Mha; negative values indicate area loss); c cumulated C net emissions (GtC; positive values indicate a C source while negative values indicate a C sink); and d cumulated net change in C-stock densities (tC/ha). See Supplementary Fig. 1 for annual fluxes.Full size imageIn line with the FRA data, we find here that the main trend is a loss of total biomass C stocks following three phases: increase in annual emissions, stagnation and slight recovery of C stocks, resulting in net C emissions from forest biomass (Fig. 1c) by 0.74 GtC or 0.03 GtC/yr between 1990 and 2020, contrasted by an opposite trend of increasing biomass density from 70 to 73 tC/ha in total forest (Fig. 1b, d). These figures are within the range of the estimated sink in forest soil and biomass of 0.1 ± 7.3 GtC/yr in 2001–2019 found by Harris et al.17. Our estimation is also consistent with that of Tubiello et al.3 of 0.11 GtC/yr net C emissions from forest ecosystems. A comparison3 of FRA-derived global forest C emissions with other independent estimates reported in 1990–2015 by National Greenhouse Gas Inventories (NGHGIs)—including the Russian Federation, the USA, China, Indonesia, and India—and by the United Nations Framework Convention on Climate Change for other countries (UNFCCC, 202018) yields a slight difference of c. 18%, although the UNFCCC and NGHGI’s account, by definition, only for emissions from managed land3. Further independent comparisons at the national and macro-regional levels are compiled in Supplementary Table 1 and reveal that C emissions estimated in the present study are in good agreement with other research.Here we find that the net C emissions mostly arise from primary forests, which undergo area loss, but also biomass thickening (Fig. 1b, d). By contrast, in spite of area loss, managed forests act as C-sinks following biomass thickening (Fig. 1b, d). Increasing biomass density is therefore key to counteract net C emissions from forest biomass in 1990–2020. While both harvest rate and burnt area increase globally over the period of observation, the increased forest growth rate that we calculate with CRAFT for both primary and managed forests over 1990–2020 emerges here as the only factor explaining increased biomass density at the global level. This is in line with other research pointing to the relevance of biomass thickening for forest C sequestration19. In addition, our finding that the forest growth rate increased annually by 0.19%, 0.21%, and 0.21% from 1990 to 2020, respectively, for primary, managed and total forests of the world is consistent with Kolby Smith et al.20 who find that also net primary production (NPP) increased annually between 0.10 and 0.25% in the period 1982–2011, as well as with other modeling and remote-sensing studies documenting a global greening trend, i.e., vegetation thickening following increased vegetation growth rate21,22. Note that estimates of annual growth rate increase in 1990–2020 by the sensitivity analyses provide narrow ranges of 0.17–0.19, 0.21–0.23, and 0.20–0.22%, respectively, for primary, managed, and total forests of the world (Supplementary Table 2).Proximate drivers of net C emissionsWe develop six counterfactual scenarios23,24,25 in order to investigate how forest biomass density and forest biomass C stocks would evolve in the hypothetical absence of (i) changes in harvest (CF1); (ii) changes in forest growth rates (CF2); (iii) change in burnt area (CF3); (iv) change in forest area (CF4); (v) harvest (CF5); (vi) burnt area (CF6) (see “Methods” section). The comparison of observed and simulated counterfactual trends allows us to isolate and quantify the influence of these four main drivers on global forest C-stock changes at national resolution (CF1 to 4) as well as to quantify the overall effects of total wood extraction and burnt area (CF5 and 6).At the global level, we find that loss of forest area (CF4) is the main driver of the net C emissions from forest biomass (Fig. 2a). In the absence of changes in area, global forest biomass would act as a cumulative net C sink of c. 26.9 GtC in the study period, creating a difference of 27.6 GtC between the actual and the CF4 C budget. This effect in the absence of area change, however, is a composite of an additional C sink of 30.7 in deforesting countries and an additional C source of 3.8 GtC in reforesting countries. Changes in harvest and burnt area from 1990 to 2020 also drove net C emissions from global forest biomass as emissions drop by c. 5.7 and 1.4 GtC in the respective counterfactual scenarios, thus generating net C-sinks of c. 4.9 and 0.63 GtC (Fig. 2a). These figures are in stark contrast with the estimated total sink of c. 49.1 and 5.4 GtC that would emerge in the hypothetical absence of harvest (CF5) and burnt area (CF6; Fig. 2a), respectively. Only changes in forest growth rates counteract the net C emissions from global forest biomass (CF2; Fig. 2a). In the absence of changes in forest growth rates, global forests would act as net C source of c. 7.4 GtC in 1990–2020, i.e., c. 10 times the actually observed source. This net effect in the absence of growth rate change results from an additional C source of 30.4 in countries experiencing growth rate increase and an additional C sink of 23.0 GtC in countries experiencing growth rate decline.Fig. 2: Counterfactual scenarios (1990–2020) assessing the cumulative impact of: changes in harvest (CF1); changes in forest growth rate (CF2); changes in burnt area (CF3); changes in forest area (CF4); total harvest (CF5); and total fire (CF6) on C-dynamics.Panels (a) and (b) show the global country-level gross and net CF C budgets (GtC) and changes in biomass density (tC/ha), respectively, with negative (red) and positive values (blue) indicating net emissions and sinks, respectively, error bars indicate the range of C budgets estimated across the five sensitivity analyses performed to test the model robustness (see Supplementary Fig. 5 for additional figures showing the net difference between CF and actual C budgets and changes in biomass density, Supplementary Table 3 and Supplementary Fig. 5 for results from sensitivity analyses). Maps show the effects of c CF1; d CF2; e CF3; f CF4; g CF5; h CF6, and are represented as the % of actual biomass C stocks that would be reached in each CF in 2020. Values above 100% (red) indicate that actual change result in net C emissions while values below 100% (blue) indicate that actual change result in a net C sink.Full size imageA sensitivity analysis on the potential underestimation of C-dynamics resulting from the use of net area change data at country level (see “Methods” section and Supplementary Fig. 5) reveals that accounting for gross area changes26 instead of net area change would result in higher global C emissions estimates (4.19 GtC in the sensitivity test versus 0.74 GtC in the reference simulation) but would reveal the same patterns of forest C-dynamic drivers (Supplementary Fig. 5). However, the magnitude of the main drivers would be slightly changed with a lower effect of changes in area (C sink in the hypothetical absence of area changes reaching 20.8 GtC in the sensitivity tests versus 26.9 GtC in the reference assessment) and a higher effect of growth rate changes (C source in the hypothetical absence of growth rate changes reaching 13.1 GtC in the sensitivity tests versus 7.4 GtC in the reference assessment). Generally, the range of results derived from the five sensitivity analyses does not change the relative importance of the individual drivers in any of the scenarios (Fig. 2a, Supplementary Table 3, and Supplementary Fig. 5). However, the sensitivity analyses highlight that the uncertainty is large enough to reverse the cumulated C signal in the absence of changes in harvest (CF1), changes in burnt area (CF3), and the complete absence of burnt areas (CF6). By contrast, the signals of CF2 (no growth rate change), CF4 (no area change), and CF5 (no harvest) are larger than the uncertainty across sensitivity analyses, signaling that our findings on these drivers are most robust.The global trends displayed in Fig. 2a, b are the combined results of diverging national forest dynamics (Fig. 2c, h). In particular, shifts in forest area (CF4) contribute to global net C emissions only in the Global South, excluding Vietnam, India, and Chile (Fig. 2f). The impacts of changes in burnt area and harvest are similarly heterogenous, with considerable effects only in some regions (e.g., Vietnam, Mozambique, Fig. 2c, e). In contrast, changes in forest growth rates are more ubiquitous, mainly positive (leading to C-sinks) for most countries, with a few notable exceptions, mainly in arid or boreal regions (e.g., India, Spain, Argentina, Canada; Fig. 2d). Possible reasons explaining the negative effect of change in forest growth rate are forest degradation, increasing drought, cloudiness, or insect outbreaks15,16,17,18,19. Over the period 1990–2020, the strongest harvest impacts are observed in countries with large area of managed forest and high harvest pressure, mostly located in temperate and subtropical areas (CF5; Fig. 2g), while fire impacts are strong in only a few countries (CF6; Fig. 2h).The fact that we use here country-level data comes both with limitations and advantages. The main limitation associated with national data is that it conceals gross C fluxes in forest biomass dynamics and blurs heterogeneity in growth conditions and anthropogenic management within countries. The country-level resolution aggregates the effects of manifold, partly counteracting processes at the local level—including photosynthesis, maintenance respiration, growth respiration, as well as forest area loss and expansion—on the annual dynamic of primary and managed forest biomass. As a consequence, our optimization of the growth function actually reflects apparent national growth rates resulting from the aggregate of these processes. However, this simplification of forest ecosystem functioning is also an advantage. Our approach reproduces forest biomass dynamics very accurately, which is complementary to most process-based models aimed at depicting biological processes and their abiotic controls27 but providing a wide range of C flux estimations1 and hardly reproducing observation from inventory data1,28,29. By contrast, the strength of the modeling approach implemented here is that it can be run with parsimonious data availability and allows to disentangle the major drivers behind forest C-stock and flux trajectories.Typology of forest biomass changeIn order to identify spatial and temporal patterns of drivers in forest biomass trends, we establish a typology of the main drivers over the period 1990–2020 (Fig. 3b). The typology we established is based on the positive versus negative shift in biomass C stocks, and highlights the most important driver of this shift as assessed through the counterfactual assessment, irrespective of the relative importance of the other drivers shown in Fig. 2. However, as the early separation between increasing and decreasing biomass C stocks in the decision tree (Fig. 3b) may conceal the effect of a major driver counteracting the observed C dynamic, the typology also accounts for possible antagonistic effects by identifying cases in which the main driver of observed C-dynamics is not, in absolute terms, the most important driver (e.g., C stocks increase but the driver with the strongest absolute effect counteracts this positive budget, see also Supplementary Fig. 3). By pinpointing the major drivers of forest change at national levels, such an approach enables to identify major levers for forest conservation.Fig. 3: Main drivers of the net C emissions from forest biomass.a Applied at the national level to the 1990–2020 period; b established according to a Boolean typology using the results from the counterfactual scenario assessment as criteria; c enabling to calculate the sum of net C-sinks and net C sources in each type of forest C-dynamics trajectory identified through the typology, error bars indicating the range of C-sinks and sources by main driver estimated across the five sensitivity analyses, with black and gray bars standing, respectively, for solid and hatched countries (see Supplementary Figs. 6–7 for results from sensitivity analyses). The hatches on the countries (a), typology (b), and bar chart (c) stand for cases in which the driver with the strongest effect actually counteracts the observed carbon budget. The color of the hatches corresponds to the main factor identified by the decision tree algorithm. Abbreviation on the typology: E: C sink driven by forest area Expansion; LH: C sink driven by Lower Harvest; FR: C sink driven by Fire Reduction; EG: C sink driven by Enhanced Growth rate; DG: C source driven by Declining Growth rate; FI: C source driven by Fire Increase; HH: C source driven by Higher Harvest; D: C source driven by Deforestation; NS: non-significant change.Full size imageDeforestation was the dominant driver of net C emissions from forest biomass in most countries of South America and Sub-Saharan Africa, corroborating findings from the literature11,30,31 (Fig. 3a, c). The net C emissions by countries where deforestation is the most significant driver reach c. 21.3 GtC, with only 0.3 GtC of these emissions being counteracted by another major driver (either increased growth rate or lower harvest pressure). These emissions represent c. 92.7% of the 21.9 GtC net emissions arising from all countries acting as net C sources (Fig. 3c). Changes in forest growth rates act as the primary drivers in most countries experiencing a net C sink over the period (Fig. 3a, c). The net C-sinks by countries where changes in forest growth rates are the main driver reach c. 16.4 GtC, with 0.9 GtC of these sinks being counteracted by another major driver (increased harvest pressure in all cases except for Sudan where area loss was the major driver counteracting the C sink). These C-sinks mainly driven by increased growth rate represented c. 77.5% of the 21.1 GtC net sink created by all countries acting as net C-sinks (Fig. 3c).Forest area expansion from 1990 to 2020 is the main driver of forest biomass net C sink in only a few Northern countries but also some Southern countries, namely Vietnam, India, and Chile, in line with findings reported for these countries32,33,34, all together accounting for a net C sink of 3.9 GtC. However, more than half of the C-sinks mainly driven by reforestation are counteracted by another major driver (either declining forest biomass growth rate or increased harvest pressure). Similarly, changes in harvest as well as changes in burnt areas are the main drivers of net C sink or source for a handful of countries in 1990–2020 (Fig. 3a). Finally, declining forest biomass growth rate is the primary driver of net C emissions only in Mongolia and Canada, which is consistent with other studies highlighting slower growth, higher mortality, and insect outbreak events in Canadian forests35,36,37.These highlights derived from the typology remained the same in all sensitivity analyses (Supplementary Figs. 6–7), despite some possible changes in country type identification (Fig. 3a and Supplementary Fig. S6) and amplitude shifts in the attribution of main drivers globally (Fig. 3c and Supplementary Fig. 7). The ranges of values in the attribution of main drivers result from the previously reported differences between the counterfactual and actual C budget estimates across sensitivity analyses (see also Supplementary Tables 2–3) combined with some changes in the type of forest C-dynamics trajectory identified through the typology in countries with large forest biomass stocks: China, India, and Australia (Supplementary Note 1 and Supplementary Fig. 6). However, these shifts do not affect the main conclusions derived from Fig. 3c: in all sensitivity analyses, growth rate changes remain the main driver of global forest biomass C sink with total net C-sinks in countries where increasing growth rate is the main driver (including both solid and hatched countries) ranging from 12.1 to 21.1 GtC, while afforestation always holds the second place of global C sink driver (total net C-sinks in countries where afforestation is the main driver ranging from 2.4 to 7.7 GtC). Similarly, total net C sources by countries where deforestation is the main driver range from −21.9 to −14.0 GtC, thus highlighting that deforestation would by far remain the main driver of forest biomass C emissions across all sensitivity analyses.Implications for forest-based solutionsOur results allow to identify major mechanisms behind observed forest biomass C changes that are immediately relevant for forest-based climate-change-mitigation strategies. We show that deforestation, increasing harvest, and burnt area have driven the net C emissions from forest biomass over the last three decades. Deforestation is the dominant driver, corroborating that protection from deforestation is indispensable1,11,38. On the other hand, forest growth rate is identified as the major driver counteracting net C emissions (Fig. 2a, d). In fact, most of the temperate and boreal countries, with the noteworthy exception of Canada, fall under a type in which enhanced forest growth rate is the major driver of a net C sink (Fig. 3b). Besides, even countries dominated by deforestation in the tropics show significant increases in growth rate (Figs. 2d and 4). These results highlight that enhanced growth rate, rather than reforestation, is the main driver counteracting biomass C emissions in 1990–2020.Fig. 4: Change in forest growth rate and its effects on global carbon stocks.The diagrams show national forest growth rate changes (y-axis) scaled along the cumulated size of the carbon stock in 1990 (x-axis). The area between the graph and the x-axis indicates the C-stock change due to growth rate for total (a), primary (b), and managed forests (c) (see Supplementary Figs. 8–10 for results from sensitivity analyses).Full size imageThese increases in forest growth rate may arise from diverse processes, including climatic and land-use drivers. On the one hand, several studies highlight the effects of environmental drivers—such as warming, atmospheric carbon dioxide (CO2), and nitrogen (N) fertilization1,6,8,11,21,39—on the terrestrial C sink. On the other hand, changes in forest growth rate can also be driven by shifts in forest management practices, such as tree species selection, forest recovery from past degradation and lesser litter grazing12,40,41. Advancing the understanding of the underlying processes of forest growth rate change is key for forging climate-change-mitigation strategies, but it is not straightforward to isolate climatic (e.g., altered CO2 concentration or temperature) from land-use drivers (e.g., non-timber forest uses such as grazing)42. Still, a comparison of trajectories in primary and managed forest growth rate change based on our results allows to derive insights into the interplay of these different drivers (Fig. 4 and Supplementary Fig. 3). From the fact that only 11% of primary forest carbon stocks show declining growth rate trends (Fig. 4c) while a relatively larger carbon stock in managed forest (22%) is affected by declining growth rate trends (Fig. 4b), we can infer that in overall terms—and assuming primary and managed forests of a given country to be similarly affected by climatic drivers —land use is likely to exert a degrading effect on growth rate dynamics. Nevertheless, some countries reveal declining growth rate in primary forest but increasing growth rate in managed forest, thus suggesting that forest management may have an improving effect on forest growth rate in those countries (e.g., USA, Fig. 4b, c, see also Supplementary Fig. 4). In overall terms, this result suggests that globally a reduction of forest use may have the potential to enhance growth rate, thus corroborating previous findings by Quesada et al.14. However, these interpretations warrant a caveat that primary versus managed forest growth rate changes are derived from the FRA data and a state-of-the-art of the literature on changes in primary forest density (see “Methods” section and Supplementary Note 2), the latter being associated with higher uncertainties although the corresponding sensitivity analysis testing suggests these uncertainties to have little impact on the figures displayed here (see Supplementary Tables 2-3 and Supplementary Figs. 5–10).Independent of their origin (management or climate driven), the future trajectories of this driver, forest growth rate, is subject to large uncertainties43,44,45. Research suggests that increasing forest growth rate is a transient phenomenon and might be discontinued in the future46. For instance, several recent studies have pointed toward the saturating effect of CO2 fertilization, which is suspected to be a key process underlying vegetation greening and ensuing thickening21, the risk of increasing mortality and slower growth rate following increasing drought6,47,48, temperature49, and natural disturbances such as insect outbreaks50,51. Even more recently, Duffy et al.52 showed that, in the near-future, temperature increases from business-as-usual trajectories of climate change shall result in a severe reduction, and possibly a reversal, of the terrestrial C sink, despite the remaining unknowns.Therefore, we conclude that, while increasing forest growth rate is the dominant driver counteracting the global net C emissions from forest biomass in the past three decades, it is against a precautionary principle to forge climate strategies that rely on a continuous net C sink effect from the same processes in the future. By contrast, our results suggest that reducing wood harvest (Fig. 2g) and halting deforestation (Fig. 2c) are key strategies to address the challenge of climate-change mitigation. In this context, increasing forest harvest volumes—a strategy often promoted in the course of climate-change-mitigation efforts embraced as the “bioeconomy”—appears to have critical unintended side-effects, despite the potential of wood for substituting some emissions-intensive products and processes53,54,55: by not only reducing the carbon sink function in forests, but also accelerating the overall C turnover rates through rejuvenation of forests and transfer to harvested wood products of lifetimes shorter than those of old-growth forests56,57,58, such strategies result in a critical loss of C sink capacity. Overall, our results plead for a double strategy to enable future forest-based solutions for climate-change mitigation: in the Global South, ending deforestation is the main priority to reverse the net C source toward a net C sink, while in the Global North, lowering wood harvest has the strongest potential to immediately enhance the C sink in forest biomass. More
