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.

Towards a dynamic list of Amazonian tree species

To provide an empirical foundation for estimates of the Amazonian tree diversity, we recently published a checklist of 11,675 tree species recorded to date in the region (ter Steege H, et al. (2016) The discovery of the Amazonian tree flora with an updated checklist of all known tree taxa. Scientific Reports 6:29549). From this total of plant records compiled from public databases and literature, widely used in studies on the Amazonian plant diversity, only 6,727 tree species belong to the first taxonomically-vetted checklist published for the region (Cardoso D, et al. (2017) Amazon plant diversity revealed by a taxonomically verified species list. PNAS 114:10695-10700). The striking difference in these two numbers spurred us to evaluate both lists, in order to release an improved Amazonian tree list; to discuss species inclusion criteria; and to highlight the ecological importance of verifying the occurrence of “non-Amazonian” trees in the region through the localization and identification of specimens. A number of species in the 2016 checklist that are not trees, non-native, synonyms, or misspellings were removed and corresponded to about 23% of the names. Species not included in the taxonomically-vetted checklist but verified by taxonomists to occur in Amazonia as trees were retained. Further, the inclusion of recently recorded/new species (after 2016), and recent taxonomic changes added up to an updated checklist including 10,071 species recorded for the Amazon region and shows the dynamic nature of establishing an authoritative checklist of Amazonian tree species. Completing and improving this list is a long-term, high-value commitment that will require a collaborative approach involving ecologists, taxonomists, and practitioners.

Variation in stem mortality rates determines patterns of above-ground biomass in Amazonian forests: implications for dynamic global vegetation models

Understanding the processes that determine above-ground biomass (AGB) in Amazonian forests is important for predicting the sensitivity of these ecosystems to environmental change and for designing and evaluating dynamic global vegetation models (DGVMs). AGB is determined by inputs from woody productivity [woody net primary productivity (NPP)] and the rate at which carbon is lost through tree mortality. Here, we test whether two direct metrics of tree mortality (the absolute rate of woody biomass loss and the rate of stem mortality) and/or woody NPP, control variation in AGB among 167 plots in intact forest across Amazonia. We then compare these relationships and the observed variation in AGB and woody NPP with the predictions of four DGVMs. The observations show that stem mortality rates, rather than absolute rates of woody biomass loss, are the most important predictor of AGB, which is consistent with the importance of stand size structure for determining spatial variation in AGB. The relationship between stem mortality rates and AGB varies among different regions of Amazonia, indicating that variation in wood density and height/diameter relationships also influences AGB. In contrast to previous findings, we find that woody NPP is not correlated with stem mortality rates and is weakly positively correlated with AGB. Across the four models, basin-wide average AGB is similar to the mean of the observations. However, the models consistently overestimate woody NPP and poorly represent the spatial patterns of both AGB and woody NPP estimated using plot data. In marked contrast to the observations, DGVMs typically show strong positive relationships between woody NPP and AGB. Resolving these differences will require incorporating forest size structure, mechanistic models of stem mortality and variation in functional composition in DGVMs.

Variation in stem mortality rates determines patterns of above-ground biomass in Amazonian forests: implications for dynamic global vegetation models

Understanding the processes that determine above-ground biomass (AGB) in Amazonian forests is important for predicting the sensitivity of these ecosystems to environmental change and for designing and evaluating dynamic global vegetation models (DGVMs). AGB is determined by inputs from woody productivity [woody net primary productivity (NPP)] and the rate at which carbon is lost through tree mortality. Here, we test whether two direct metrics of tree mortality (the absolute rate of woody biomass loss and the rate of stem mortality) and/or woody NPP, control variation in AGB among 167 plots in intact forest across Amazonia. We then compare these relationships and the observed variation in AGB and woody NPP with the predictions of four DGVMs. The observations show that stem mortality rates, rather than absolute rates of woody biomass loss, are the most important predictor of AGB, which is consistent with the importance of stand size structure for determining spatial variation in AGB. The relationship between stem mortality rates and AGB varies among different regions of Amazonia, indicating that variation in wood density and height/diameter relationships also influences AGB. In contrast to previous findings, we find that woody NPP is not correlated with stem mortality rates and is weakly positively correlated with AGB. Across the four models, basin-wide average AGB is similar to the mean of the observations. However, the models consistently overestimate woody NPP and poorly represent the spatial patterns of both AGB and woody NPP estimated using plot data. In marked contrast to the observations, DGVMs typically show strong positive relationships between woody NPP and AGB. Resolving these differences will require incorporating forest size structure, mechanistic models of stem mortality and variation in functional composition in DGVMs.