Integrated global assessment of the natural forest carbon potential

Forests are a substantial terrestrial carbon sink, but anthropogenic changes in land use and climate have considerably reduced the scale of this system1. Remote-sensing estimates to quantify carbon losses from global forests2-5 are characterized by considerable uncertainty and we lack a comprehensive ground-sourced evaluation to benchmark these estimates. Here we combine several ground-sourced6 and satellite-derived approaches2,7,8 to evaluate the scale of the global forest carbon potential outside agricultural and urban lands. Despite regional variation, the predictions demonstrated remarkable consistency at a global scale, with only a 12% difference between the ground-sourced and satellite-derived estimates. At present, global forest carbon storage is markedly under the natural potential, with a total deficit of 226 Gt (model range = 151-363 Gt) in areas with low human footprint. Most (61%, 139 Gt C) of this potential is in areas with existing forests, in which ecosystem protection can allow forests to recover to maturity. The remaining 39% (87 Gt C) of potential lies in regions in which forests have been removed or fragmented. Although forests cannot be a substitute for emissions reductions, our results support the idea2,3,9 that the conservation, restoration and sustainable management of diverse forests offer valuable contributions to meeting global climate and biodiversity targets.

The global biogeography of tree leaf form and habit

Understanding what controls global leaf type variation in trees is crucial for comprehending their role in terrestrial ecosystems, including carbon, water and nutrient dynamics. Yet our understanding of the factors influencing forest leaf types remains incomplete, leaving us uncertain about the global proportions of needle-leaved, broadleaved, evergreen and deciduous trees. To address these gaps, we conducted a global, ground-sourced assessment of forest leaf-type variation by integrating forest inventory data with comprehensive leaf form (broadleaf vs needle-leaf) and habit (evergreen vs deciduous) records. We found that global variation in leaf habit is primarily driven by isothermality and soil characteristics, while leaf form is predominantly driven by temperature. Given these relationships, we estimate that 38% of global tree individuals are needle-leaved evergreen, 29% are broadleaved evergreen, 27% are broadleaved deciduous and 5% are needle-leaved deciduous. The aboveground biomass distribution among these tree types is approximately 21% (126.4 Gt), 54% (335.7 Gt), 22% (136.2 Gt) and 3% (18.7 Gt), respectively. We further project that, depending on future emissions pathways, 17-34% of forested areas will experience climate conditions by the end of the century that currently support a different forest type, highlighting the intensification of climatic stress on existing forests. By quantifying the distribution of tree leaf types and their corresponding biomass, and identifying regions where climate change will exert greatest pressure on current leaf types, our results can help improve predictions of future terrestrial ecosystem functioning and carbon cycling.

Native diversity buffers against severity of non-native tree invasions

Determining the drivers of non-native plant invasions is critical for managing native ecosystems and limiting the spread of invasive species1,2. Tree invasions in particular have been relatively overlooked, even though they have the potential to transform ecosystems and economies3,4. Here, leveraging global tree databases5-7, we explore how the phylogenetic and functional diversity of native tree communities, human pressure and the environment influence the establishment of non-native tree species and the subsequent invasion severity. We find that anthropogenic factors are key to predicting whether a location is invaded, but that invasion severity is underpinned by native diversity, with higher diversity predicting lower invasion severity. Temperature and precipitation emerge as strong predictors of invasion strategy, with non-native species invading successfully when they are similar to the native community in cold or dry extremes. Yet, despite the influence of these ecological forces in determining invasion strategy, we find evidence that these patterns can be obscured by human activity, with lower ecological signal in areas with higher proximity to shipping ports. Our global perspective of non-native tree invasion highlights that human drivers influence non-native tree presence, and that native phylogenetic and functional diversity have a critical role in the establishment and spread of subsequent invasions.

Co-limitation towards lower latitudes shapes global forest diversity gradients

The latitudinal diversity gradient (LDG) is one of the most recognized global patterns of species richness exhibited across a wide range of taxa. Numerous hypotheses have been proposed in the past two centuries to explain LDG, but rigorous tests of the drivers of LDGs have been limited by a lack of high-quality global species richness data. Here we produce a high-resolution (0.025° × 0.025°) map of local tree species richness using a global forest inventory database with individual tree information and local biophysical characteristics from ~1.3 million sample plots. We then quantify drivers of local tree species richness patterns across latitudes. Generally, annual mean temperature was a dominant predictor of tree species richness, which is most consistent with the metabolic theory of biodiversity (MTB). However, MTB underestimated LDG in the tropics, where high species richness was also moderated by topographic, soil and anthropogenic factors operating at local scales. Given that local landscape variables operate synergistically with bioclimatic factors in shaping the global LDG pattern, we suggest that MTB be extended to account for co-limitation by subordinate drivers.

The 2018 European heatwave led to stem dehydration but not to consistent growth reductions in forests

Heatwaves exert disproportionately strong and sometimes irreversible impacts on forest ecosystems. These impacts remain poorly understood at the tree and species level and across large spatial scales. Here, we investigate the effects of the record-breaking 2018 European heatwave on tree growth and tree water status using a collection of high-temporal resolution dendrometer data from 21 species across 53 sites. Relative to the two preceding years, annual stem growth was not consistently reduced by the 2018 heatwave but stems experienced twice the temporary shrinkage due to depletion of water reserves. Conifer species were less capable of rehydrating overnight than broadleaves across gradients of soil and atmospheric drought, suggesting less resilience toward transient stress. In particular, Norway spruce and Scots pine experienced extensive stem dehydration. Our high-resolution dendrometer network was suitable to disentangle the effects of a severe heatwave on tree growth and desiccation at large-spatial scales in situ, and provided insights on which species may be more vulnerable to climate extremes.

Two large-scale forest scenario modelling approaches for reporting CO2 removal: a comparison for the Romanian forests

BACKGROUND

Forest carbon models are recognized as suitable tools for the reporting and verification of forest carbon stock and stock change, as well as for evaluating the forest management options to enhance the carbon sink provided by sustainable forestry. However, given their increased complexity and data availability, different models may simulate different estimates. Here, we compare carbon estimates for Romanian forests as simulated by two models (CBM and EFISCEN) that are often used for evaluating the mitigation options given the forest-management choices.

RESULTS

The models, calibrated and parameterized with identical or harmonized data, derived from two successive national forest inventories, produced similar estimates of carbon accumulation in tree biomass. According to CBM simulations of carbon stocks in Romanian forests, by 2060, the merchantable standing stock volume will reach an average of 377 m³ ha^-1, while the carbon stock in tree biomass will reach 76.5 tC ha^-1. The EFISCEN simulations produced estimates that are about 5% and 10%, respectively, lower. In addition, 10% stronger biomass sink was simulated by CBM, whereby the difference reduced over time, amounting to only 3% toward 2060.

CONCLUSIONS

This model comparison provided valuable insights on both the conceptual and modelling algorithms, as well as how the quality of the input data may affect calibration and projections of the stock and stock change in the living biomass pool. In our judgement, both models performed well, providing internally consistent results. Therefore, we underline the importance of the input data quality and the need for further data sampling and model improvements, while the preference for one model or the other should be based on the availability and suitability of the required data, on preferred output variables and ease of use.

Changes in global terrestrial live biomass over the 21st century

Live woody vegetation is the largest reservoir of biomass carbon, with its restoration considered one of the most effective natural climate solutions. However, terrestrial carbon fluxes remain the largest uncertainty in the global carbon cycle. Here, we develop spatially explicit estimates of carbon stock changes of live woody biomass from 2000 to 2019 using measurements from ground, air, and space. We show that live biomass has removed 4.9 to 5.5 PgC year^-1 from the atmosphere, offsetting 4.6 ± 0.1 PgC year^-1 of gross emissions from disturbances and adding substantially (0.23 to 0.88 PgC year^-1) to the global carbon stocks. Gross emissions and removals in the tropics were four times larger than temperate and boreal ecosystems combined. Although live biomass is responsible for more than 80% of gross terrestrial fluxes, soil, dead organic matter, and lateral transport may play important roles in terrestrial carbon sink.

Actual European forest management by region, tree species and owner based on 714,000 re-measured trees in national forest inventories

BACKGROUND

European forests have a long record of management. However, the diversity of the current forest management across nations, tree species and owners, is hardly understood. Often when trying to simulate future forest resources under alternative futures, simply the yield table style of harvesting is applied. It is now crucially important to come to grips with actual forest management, now that demand for wood is increasing and the EU Land Use, Land Use Change and Forestry Regulation has been adopted requiring 'continuation of current management practices' as a baseline to set the Forest Reference Level carbon sink.

METHODS

Based on a large dataset of 714,000 re-measured trees in National Forest inventories from 13 regions, we are now able to analyse actual forest harvesting.

CONCLUSIONS

From this large set of repeated tree measurements we can conclude that there is no such thing as yield table harvesting in Europe. We found general trends of increasing harvest probability with higher productivity of the region and the species, but with important deviations related to local conditions like site accessibility, state of the forest resource (like age), specific subsidies, importance of other forest services, and ownership of the forest. As a result, we find a huge diversity in harvest regimes. Over the time period covered in our inventories, the average harvest probability over all regions was 2.4% yr-1 (in number of trees) and the mortality probability was 0.4% yr-1. Our study provides underlying and most actual data that can serve as a basis for quantifying 'continuation of current forest management'. It can be used as a cornerstone for the base period as required for the Forest Reference Level for EU Member States.

Understanding the implications of the EU-LULUCF regulation for the wood supply from EU forests to the EU

BACKGROUND

In June 2018, the European Parliament and Council of the European Union adopted a legislative regulation for incorporating greenhouse gas emissions and removals from Land Use, Land Use Change and Forestry (EU-LULUCF) under its 2030 Climate and Energy Framework. The LULUCF regulation aim to incentivise EU Member States to decrease greenhouse gas emissions and increase removals in the LULUCF sector. The regulation, however, does not set a target for increasing the LULUCF carbon sink, but rather includes a 'no net debit' target for LULUCF (Forests and Agricultural soils). For Managed Forest Land (MFL) an accounting framework with capped credits for additional mitigation against a set forest reference level (FRL) was agreed for 2021-2030. The FRL gives the projected future carbon sink in the two compliance periods 2021-2025 and 2026-2030 under "continuation of forest management practices as they were in the reference period 2000-2009". This FRL was disputed by some Member States as it was perceived to put a limit on their future wood harvesting from MFL. Here we simulated with the EFISCEN European forest model the "continuation of forest management practices" and determined the corresponding wood harvest for 26 EU countries under progressing age classes.

RESULTS

The simulations showed that under "continuation of forest management practices" the harvest (wood removals) in the 26 EU countries as a whole can increase from 420 million m³/year in 2000-2009 to 560 million m³/year in 2050 due to progressing age classes. This implies there is a possibility to increase absolute wood harvests without creating debits compared to the forest reference level. However, the manner in which 'continuation of forest management' developed with a progressing age class development over time, meant that in some countries the future harvesting exceeded 90% of the increment. Since this generally is considered to be unsustainable we additionally set a harvesting cut-off as max 90% of increment to be harvested for each individual country as a possible interpretation of sustainability criteria that are included in the regulation. Using this additional limit the projected harvest will only increase to 493 million m³/year.

CONCLUSIONS

The worry from Member States (MS) that the FRL will prevent any additional harvesting seems unwarranted. Due to differences between Member States concerning the state of their forest resources, the FRL as a baseline for harvesting works out very differently for the different Member States. The FRL may have other unforeseen consequences which we discuss. Under all scenarios the living forest biomass sink shows a decline. This can be counteracted through incentivising measures under Climate Smart Forestry.

Species distribution models predict temporal but not spatial variation in forest growth

Bioclimate envelope models have been widely used to illustrate the discrepancy between current species distributions and their potential habitat under climate change. However, the realism and correct interpretation of such projections has been the subject of considerable discussion. Here, we investigate whether climate suitability predictions correlate to tree growth, measured in permanent inventory plots and inferred from tree‐ring records. We use the ensemble classifier RandomForest and species occurrence data from ~200,000 inventory plots to build species distribution models for four important European forestry species: Norway spruce, Scots pine, European beech, and pedunculate oak. We then correlate climate‐based habitat suitability with volume measurements from ~50‐year‐old stands, available from ~11,000 inventory plots. Secondly, habitat projections based on annual historical climate are compared with ring width from ~300 tree‐ring chronologies. Our working hypothesis is that habitat suitability projections from species distribution models should to some degree be associated with temporal or spatial variation in these growth records. We find that the habitat projections are uncorrelated with spatial growth records (inventory plot data), but they do predict interannual variation in tree‐ring width, with an average correlation of .22. Correlation coefficients for individual chronologies range from values as high as .82 or as low as −.31. We conclude that tree responses to projected climate change are highly site‐specific and that local suitability of a species for reforestation is difficult to predict. That said, projected increase or decrease in climatic suitability may be interpreted as an average expectation of increased or reduced growth over larger geographic scales.

Improving national-scale carbon stock inventories using knowledge on land use history

National-scale inventories of soil organic carbon (SOC) and forest floor carbon (FFC) stocks have a high uncertainty. Inventories are often based on the interpolation of sampled information, often using a number of covariables to help such interpolation. The rationale for the choice of these covariables is not always documented, despite the fact that many local-scale studies have identified the factors explaining spatial variability of SOC and FFC stocks. These studies indicate, among others the importance of long-term land use history. Despite this, information on the effects of land use history has never been used to explain variability of carbon stocks in national-scale inventories. We designed an alternative method to improve national-scale inventories of SOC and FCC for the Dutch sand area that takes stock of the findings of detailed case studies. Determinants for SOC and FFC stocks derived from landscape-scale case studies were used to map national-scale spatial variability and to calculate national totals. The resulting national-scale spatial distribution was compared with the SOC stock map from the current Dutch greenhouse gas inventory. Using land use history to explain SOC variability decreased the error of the SOC stock estimate in 60 % of the area. The error in FFC stocks decreased in half of the forest area after including soil fertility, tree species, and forest age as explanatory factors. Estimates with reduced uncertainty will make land use and land management a more attractive and acceptable mitigation option to reduce emissions of greenhouse gases for the LULUCF sector.

Ecologically implausible carbon response?

Magnani et al. present a very strong correlation between mean lifetime net ecosystem production (NEP, defined as the net rate of carbon (C) accumulation in ecosystems) and wet nitrogen (N) deposition. For their data in the range 4.9-9.8 kg N ha(-1) yr(-1), on which the correlation largely depends, the response is approximately 725 kg C per kg N in wet deposition. According to the authors, the maximum N wet deposition level of 9.8 kg N ha(-1) yr(-1) is equivalent to a total deposition of 15 kg N ha(-1 )yr(-1), implying a net sequestration near 470 kg C per kg N of total deposition. We question the ecological plausibility of the relationship and show, from a multi-factor analysis of European forest measurements, how interactions with site productivity and environment imply a much smaller NEP response to N deposition.

Europe's terrestrial biosphere absorbs 7 to 12% of European anthropogenic CO2 emissions

Most inverse atmospheric models report considerable uptake of carbon dioxide in Europe's terrestrial biosphere. In contrast, carbon stocks in terrestrial ecosystems increase at a much smaller rate, with carbon gains in forests and grassland soils almost being offset by carbon losses from cropland and peat soils. Accounting for non-carbon dioxide carbon transfers that are not detected by the atmospheric models and for carbon dioxide fluxes bypassing the ecosystem carbon stocks considerably reduces the gap between the small carbon-stock changes and the larger carbon dioxide uptake estimated by atmospheric models. The remaining difference could be because of missing components in the stock-change approach, as well as the large uncertainty in both methods. With the use of the corrected atmosphere- and land-based estimates as a dual constraint, we estimate a net carbon sink between 135 and 205 teragrams per year in Europe's terrestrial biosphere, the equivalent of 7 to 12% of the 1995 anthropogenic carbon emissions.

CARBON SINKS IN TEMPERATE FORESTS

▪ Abstract  In addition to being scientifically exciting, commercially important, and environmentally essential, temperate forests have also become a key diplomatic item in international climate negotiations as potential sinks for carbon. This review presents the methods used to estimate carbon sequestration, identifies the constraints and opportunities for carbon sequestration in temperate forests, addresses the issues raised by the monitoring of carbon sequestration, and analyzes uncertainties pertaining to the sequestration of carbon by temperate forests. This review serves a dual purpose: It aims at informing policy makers about carbon sequestration in temperate forests and at making forest ecologists, biogeochemists, and atmospheric scientists aware of the structure of an international agreement to reduce CO2 and other greenhouse gas emissions and some of the real, still answered scientific questions that it poses.