Evaluation of climate-related carbon turnover processes in global vegetation models for boreal and temperate forests

The evolution, complexity and diversity of models of long-term forest dynamics

To assess the impacts of climate change on vegetation from stand to global scales, models of forest dynamics that include tree demography are needed. Such models are now available for 50 years, but the currently existing diversity of model formulations and its evolution over time are poorly documented. This hampers systematic assessments of structural uncertainties in model-based studies.We conducted a meta-analysis of 28 models, focusing on models that were used in the past five years for climate change studies. We defined 52 model attributes in five groups (basic assumptions, growth, regeneration, mortality and soil moisture) and characterized each model according to these attributes. Analyses of model complexity and diversity included hierarchical cluster analysis and redundancy analysis.Model complexity evolved considerably over the past 50 years. Increases in complexity were largest for growth processes, while complexity of modelled establishment processes increased only moderately. Model diversity was lowest at the global scale, and highest at the landscape scale. We identified five distinct clusters of models, ranging from very simple models to models where specific attribute groups are rendered in a complex manner and models that feature high complexity across all attributes.Most models in use today are not balanced in the level of complexity with which they represent different processes. This is the result of different model purposes, but also reflects legacies in model code, modelers' preferences, and the 'prevailing spirit of the epoch'. The lack of firm theories, laws and 'first principles' in ecology provides high degrees of freedom in model development, but also results in high responsibilities for model developers and the need for rigorous model evaluation. Synthesis. The currently available model diversity is beneficial: convergence in simulations of structurally different models indicates robust projections, while convergence of similar models may convey a false sense of certainty. The existing model diversity-with the exception of global models-can be exploited for improved projections based on multiple models. We strongly recommend balanced further developments of forest models that should particularly focus on establishment and mortality processes, in order to provide robust information for decisions in ecosystem management and policymaking.

A data-driven estimate of litterfall and forest carbon turnover and the drivers of their inter-annual variabilities in forest ecosystems across China

Strong influences of climate and land-cover changes on terrestrial ecosystems urgently need to re-estimate forest carbon turnover time (τ_forest), i.e., the residence time of carbon (C) in the living forest carbon reservoir in China, to reduce uncertainties in ecosystem carbon sinks under ongoing climate change. However, in absence of accurate carbon loss (e.g., forest litterfall), τ_forest estimate based on the non-steady-state assumption (NSSA) in forest ecosystems across China is still unclear. In this study, thus, we first compiled a litterfall dataset with 1025 field observations, and applied a Random Forest (RF) algorithm with the linkage of gridded environmental variables to predict litterfall from 2000 to 2019 with a fine spatial resolution of 1 km and a temporal resolution of one year. Finally, τ_forest was also estimated with the data-driven litterfall product. Results showed that RF algorithm could well predict the spatial and temporal patterns of forest litterfall with a model efficiency of 0.58 and root mean square error of 78.7 g C m^-2 year^-1. Mean litterfall was 205.4 ± 1.1 Tg C year^-1 (mean ± standard error) with a significant increasing trend of 0.65 ± 0.14 Tg C year^-2 from 2000 to 2019 (p < 0.01), indicating an increasing carbon loss from litterfall. Mean τ_forest was 26.2 ± 0.1 years with a significant decreasing trend of -0.11 ± 0.02 years (p < 0.01) from 2000 to 2019. Climate change dominated the inter-annual variability of τ_forest in high latitude areas, and land-cover change dominated the regions with intensive human activities. These findings suggested that carbon loss from vegetation to the atmosphere becomes more quickly in recent decades, with significant implication for vegetation carbon cycling-climate feedbacks. Meanwhile, the developed litterfall and τ_forest datasets can serve as a benchmark for biogeochemical models to accurately estimate global carbon cycling.

Evidence for conifer sucrose transporters' functioning in the light-dependent adjustment of sugar allocation

Sucrose is the central unit of carbon and energy in plants. Active intercellular transport of sucrose is mediated by sucrose transporters (SUTs), genes for which have been found in the genomes of all land plants. However, they have only been assigned functions in angiosperm species. Here, we cloned two types of SUTs from two gymnosperms, the conifers Cedrus deodara (Roxb. G. Don) and Pinus massoniana Lambert, and analyzed their sucrose transport activities. Uptake of the fluorescent sucrose-analog esculin into tobacco epidermis cells expressing the conifer SUT confirmed their transport ability. To determine their function in planta, we investigated their mRNA abundance in relation to photosynthesis and sugar levels in leaves, inner bark, wood and roots. Combined with measurements of protein abundance and immunolocalization of C. deodara SUTs, our results suggest a role for CdSUT1G and CdSUT2 in supporting phloem transport under varying environmental conditions. The implications of these findings regarding conifer physiology and SUT evolution are discussed.

Run to the hills: Forest growth responsiveness to drought increased at higher elevation during the late 20th century

Climate warming is expected to enhance forest growth in cold-limited biomes while triggering reductions in drought-limited biomes. However, as temperature raises, it is unclear how temperature- and drought-growth couplings shift across elevation gradients in different biomes. We still lack comprehensive analyses on how altitude modulates the influence of temperature and drought on tree growth during the second half of the 20th century when climate warming accelerated. We compared the worldwide responses of tree growth (RWI, ring-width indices) to two of its major climatic constraints, growing-season minimum temperatures and drought (SPEI index), across biomes and elevation gradients during two periods with different warming rates (1960-1980 vs. 1980-2000). We found a decrease in the correlations of minimum temperatures with growth, but a strengthening of drought-growth relationships. However, these patterns varied across biomes because correlations between growth and temperature decreased in temperate forests and woodland shrubland, while correlations between growth and SPEI increased in boreal forests and decreased in temperate forests. Differences in growth responsiveness to climate between the two periods were more marked for mid-latitude forests situated between 1200 and 1600 m. The slopes of the relationships between growth-temperature correlations and elevation decreased in late spring and midsummer. The slopes of the relationships between growth-drought correlations and elevation increased in temperate forests and woodland shrubland suggesting that drought impacts are "climbing" in these biomes. Temperature controls on forest growth are relaxing as the climate warms, while drought is becoming a more significant constraint for tree growth, particularly for mid-elevation forests and in drought-prone woodland and shrubland. The strengthening of drought-growth coupling should be considered in vegetation models to reduce the uncertainty on forest climate mitigation.

Data-driven estimates of global litter production imply slower vegetation carbon turnover

Accurate quantification of vegetation carbon turnover time (τ_veg ) is critical for reducing uncertainties in terrestrial vegetation response to future climate change. However, in the absence of global information of litter production, τ_veg could only be estimated based on net primary productivity under the steady-state assumption. Here, we applied a machine-learning approach to derive a global dataset of litter production by linking 2401 field observations and global environmental drivers. Results suggested that the observation-based estimate of global natural ecosystem litter production was 44.3 ± 0.4 Pg C year^-1 . By contrast, land-surface models (LSMs) overestimated the global litter production by about 27%. With this new global litter production dataset, we estimated global τ_veg (mean value 10.3 ± 1.4 years) and its spatial distribution. Compared to our observation-based τ_veg , modelled τ_veg tended to underestimate τ_veg at high latitudes. Our empirically derived gridded datasets of litter production and τ_veg will help constrain global vegetation models and improve the prediction of global carbon cycle.

Excessive positive response of model-simulated land net primary production to climate changes over circumboreal forests

Land carbon cycle components in an Earth system model (ESM) play a crucial role in the projections of forest ecosystem responses to climate/environmental changes. Evaluating models from the viewpoint of observations is essential for an improved understanding of model performance and for identifying uncertainties in their outputs. Herein, we evaluated the land net primary production (NPP) for circumboreal forests simulated with 10 ESMs in Phase 5 of the Coupled Model Intercomparison Project by comparisons with observation-based indexes for forest productivity, namely, the composite version 3G of the normalized difference vegetation index (NDVI3g) and tree-ring width index (RWI). These indexes show similar patterns in response to past climate change over the forests, i.e., a one-year time lag response and smaller positive responses to past climate changes in comparison with the land NPP simulated by the ESMs. The latter showed overly positive responses to past temperature and/or precipitation changes in comparison with the NDVI3g and RWI. These results indicate that ESMs may overestimate the future forest NPP of circumboreal forests (particularly for inland dry regions, such as inner Alaska and Canada, and eastern Siberia, and for hotter, southern regions, such as central Europe) under the expected increases in both average global temperature and precipitation, which are common to all current ESMs.

Rainfall manipulation experiments as simulated by terrestrial biosphere models: Where do we stand?

Changes in rainfall amounts and patterns have been observed and are expected to continue in the near future with potentially significant ecological and societal consequences. Modelling vegetation responses to changes in rainfall is thus crucial to project water and carbon cycles in the future. In this study, we present the results of a new model-data intercomparison project, where we tested the ability of 10 terrestrial biosphere models to reproduce the observed sensitivity of ecosystem productivity to rainfall changes at 10 sites across the globe, in nine of which, rainfall exclusion and/or irrigation experiments had been performed. The key results are as follows: (a) Inter-model variation is generally large and model agreement varies with timescales. In severely water-limited sites, models only agree on the interannual variability of evapotranspiration and to a smaller extent on gross primary productivity. In more mesic sites, model agreement for both water and carbon fluxes is typically higher on fine (daily-monthly) timescales and reduces on longer (seasonal-annual) scales. (b) Models on average overestimate the relationship between ecosystem productivity and mean rainfall amounts across sites (in space) and have a low capacity in reproducing the temporal (interannual) sensitivity of vegetation productivity to annual rainfall at a given site, even though observation uncertainty is comparable to inter-model variability. (c) Most models reproduced the sign of the observed patterns in productivity changes in rainfall manipulation experiments but had a low capacity in reproducing the observed magnitude of productivity changes. Models better reproduced the observed productivity responses due to rainfall exclusion than addition. (d) All models attribute ecosystem productivity changes to the intensity of vegetation stress and peak leaf area, whereas the impact of the change in growing season length is negligible. The relative contribution of the peak leaf area and vegetation stress intensity was highly variable among models.

Organizing principles for vegetation dynamics

Plants and vegetation play a critical-but largely unpredictable-role in global environmental changes due to the multitude of contributing processes at widely different spatial and temporal scales. In this Perspective, we explore approaches to master this complexity and improve our ability to predict vegetation dynamics by explicitly taking account of principles that constrain plant and ecosystem behaviour: natural selection, self-organization and entropy maximization. These ideas are increasingly being used in vegetation models, but we argue that their full potential has yet to be realized. We demonstrate the power of natural selection-based optimality principles to predict photosynthetic and carbon allocation responses to multiple environmental drivers, as well as how individual plasticity leads to the predictable self-organization of forest canopies. We show how models of natural selection acting on a few key traits can generate realistic plant communities and how entropy maximization can identify the most probable outcomes of community dynamics in space- and time-varying environments. Finally, we present a roadmap indicating how these principles could be combined in a new generation of models with stronger theoretical foundations and an improved capacity to predict complex vegetation responses to environmental change.

Plant respiration: Controlled by photosynthesis or biomass?

Two simplifying hypotheses have been proposed for whole-plant respiration. One links respiration to photosynthesis; the other to biomass. Using a first-principles carbon balance model with a prescribed live woody biomass turnover, applied at a forest research site where multidecadal measurements are available for comparison, we show that if turnover is fast the accumulation of respiring biomass is low and respiration depends primarily on photosynthesis; while if turnover is slow the accumulation of respiring biomass is high and respiration depends primarily on biomass. But the first scenario is inconsistent with evidence for substantial carry-over of fixed carbon between years, while the second implies far too great an increase in respiration during stand development-leading to depleted carbohydrate reserves and an unrealistically high mortality risk. These two mutually incompatible hypotheses are thus both incorrect. Respiration is not linearly related either to photosynthesis or to biomass, but it is more strongly controlled by recent photosynthates (and reserve availability) than by total biomass.

Constraining modelled global vegetation dynamics and carbon turnover using multiple satellite observations

The response of land ecosystems to future climate change is among the largest unknowns in the global climate-carbon cycle feedback. This uncertainty originates from how dynamic global vegetation models (DGVMs) simulate climate impacts on changes in vegetation distribution, productivity, biomass allocation, and carbon turnover. The present-day availability of a multitude of satellite observations can potentially help to constrain DGVM simulations within model-data integration frameworks. Here, we use satellite-derived datasets of the fraction of absorbed photosynthetic active radiation (FAPAR), sun-induced fluorescence (SIF), above-ground biomass of trees (AGB), land cover, and burned area to constrain parameters for phenology, productivity, and vegetation dynamics in the LPJmL4 DGVM. Both the prior and the optimized model accurately reproduce present-day estimates of the land carbon cycle and of temporal dynamics in FAPAR, SIF and gross primary production. However, the optimized model reproduces better the observed spatial patterns of biomass, tree cover, and regional forest carbon turnover. Using a machine learning approach, we found that remaining errors in simulated forest carbon turnover can be explained with bioclimatic variables. This demonstrates the need to improve model formulations for climate effects on vegetation turnover and mortality despite the apparent successful constraint of simulated vegetation dynamics with multiple satellite observations.

Pervasive decreases in living vegetation carbon turnover time across forest climate zones

Forests play a major role in the global carbon cycle. Previous studies on the capacity of forests to sequester atmospheric CO2 have mostly focused on carbon uptake, but the roles of carbon turnover time and its spatiotemporal changes remain poorly understood. Here, we used long-term inventory data (1955 to 2018) from 695 mature forest plots to quantify temporal trends in living vegetation carbon turnover time across tropical, temperate, and cold climate zones, and compared plot data to 8 Earth system models (ESMs). Long-term plots consistently showed decreases in living vegetation carbon turnover time, likely driven by increased tree mortality across all major climate zones. Changes in living vegetation carbon turnover time were negatively correlated with CO2 enrichment in both forest plot data and ESM simulations. However, plot-based correlations between living vegetation carbon turnover time and climate drivers such as precipitation and temperature diverged from those of ESM simulations. Our analyses suggest that forest carbon sinks are likely to be constrained by a decrease in living vegetation carbon turnover time, and accurate projections of forest carbon sink dynamics will require an improved representation of tree mortality processes and their sensitivity to climate in ESMs.

Underestimated ecosystem carbon turnover time and sequestration under the steady state assumption: A perspective from long-term data assimilation

It is critical to accurately estimate carbon (C) turnover time as it dominates the uncertainty in ecosystem C sinks and their response to future climate change. In the absence of direct observations of ecosystem C losses, C turnover times are commonly estimated under the steady state assumption (SSA), which has been applied across a large range of temporal and spatial scales including many at which the validity of the assumption is likely to be violated. However, the errors associated with improperly applying SSA to estimate C turnover time and its covariance with climate as well as ecosystem C sequestrations have yet to be fully quantified. Here, we developed a novel model-data fusion framework and systematically analyzed the SSA-induced biases using time-series data collected from 10 permanent forest plots in the eastern China monsoon region. The results showed that (a) the SSA significantly underestimated mean turnover times (MTTs) by 29%, thereby leading to a 4.83-fold underestimation of the net ecosystem productivity (NEP) in these forest ecosystems, a major C sink globally; (b) the SSA-induced bias in MTT and NEP correlates negatively with forest age, which provides a significant caveat for applying the SSA to young-aged ecosystems; and (c) the sensitivity of MTT to temperature and precipitation was 22% and 42% lower, respectively, under the SSA. Thus, under the expected climate change, spatiotemporal changes in MTT are likely to be underestimated, thereby resulting in large errors in the variability of predicted global NEP. With the development of observation technology and the accumulation of spatiotemporal data, we suggest estimating MTTs at the disequilibrium state via long-term data assimilation, thereby effectively reducing the uncertainty in ecosystem C sequestration estimations and providing a better understanding of regional or global C cycle dynamics and C-climate feedback.

Twentieth century redistribution in climatic drivers of global tree growth

Energy and water limitations of tree growth remain insufficiently understood at large spatiotemporal scales, hindering model representation of interannual or longer-term ecosystem processes. By assessing and statistically scaling the climatic drivers from 2710 tree-ring sites, we identified the boreal and temperate land areas where tree growth during 1930-1960 CE responded positively to temperature (20.8 ± 3.7 Mio km²; 25.9 ± 4.6%), precipitation (77.5 ± 3.3 Mio km²; 96.4 ± 4.1%), and other parameters. The spatial manifestation of this climate response is determined by latitudinal and altitudinal temperature gradients, indicating that warming leads to geographic shifts in growth limitations. We observed a significant (P < 0.001) decrease in temperature response at cold-dry sites between 1930-1960 and 1960-1990 CE, and the total temperature-limited area shrunk by -8.7 ± 0.6 Mio km². Simultaneously, trees became more limited by atmospheric water demand almost worldwide. These changes occurred under mild warming, and we expect that continued climate change will trigger a major redistribution in growth responses to climate.

Estimation of Above-Ground Biomass over Boreal Forests on Siberia Using Updated In Situ, ALOS-2 PALSAR-2, and RADARSAT-2 Data

The estimation of above-ground biomass (AGB) in boreal forests is of special concern as it constitutes the highest carbon pool in the northern hemisphere. In particularly, monitoring of the forests in the Russian Federation is important as some regions have not been inventoried for many years. This study explores the combination of multi-frequency, multi-polarization, and multi-temporal radar data as one key approach to provide an accurate estimate of forest biomass. The data from L-band Advanced Land Observing Satellite 2 (ALOS-2) Phased Array L-Band Synthetic Aperture Radar 2 (PALSAR-2), together with C-band RADARSAT-2 data, were applied for AGB estimation. Backscatter coefficients from L- and C-band radar were used independently and in combination with a non-parametric model to retrieve AGB data for a boreal forest in Siberia (Krasnoyarskiy Kray). AGB estimation was performed using the random forests machine learning algorithm. The results demonstrated that high estimation accuracies can be achieved at a spatial resolution of 0.25 ha. When the L-band data alone were used for the retrieval, a corrected root-mean-square error (RMSEcor) of 29.4 t ha−1 was calculated. A marginal decrease in RMSEcor was observed when only the filtered L-band backscatter data, without ratio and texture, were used (29.1 t ha−1). The inclusion of the C-band data reduced the over and underestimation; the bias was reduced from 5.5 t ha−1 to 4.7 t ha−1; and a RMSEcor of 30.2 t ha−1 was calculated.

Drivers and mechanisms of tree mortality in moist tropical forests

Tree mortality rates appear to be increasing in moist tropical forests (MTFs) with significant carbon cycle consequences. Here, we review the state of knowledge regarding MTF tree mortality, create a conceptual framework with testable hypotheses regarding the drivers, mechanisms and interactions that may underlie increasing MTF mortality rates, and identify the next steps for improved understanding and reduced prediction. Increasing mortality rates are associated with rising temperature and vapor pressure deficit, liana abundance, drought, wind events, fire and, possibly, CO₂ fertilization-induced increases in stand thinning or acceleration of trees reaching larger, more vulnerable heights. The majority of these mortality drivers may kill trees in part through carbon starvation and hydraulic failure. The relative importance of each driver is unknown. High species diversity may buffer MTFs against large-scale mortality events, but recent and expected trends in mortality drivers give reason for concern regarding increasing mortality within MTFs. Models of tropical tree mortality are advancing the representation of hydraulics, carbon and demography, but require more empirical knowledge regarding the most common drivers and their subsequent mechanisms. We outline critical datasets and model developments required to test hypotheses regarding the underlying causes of increasing MTF mortality rates, and improve prediction of future mortality under climate change.

Uncertainty Quantification of Extratropical Forest Biomass in CMIP5 Models over the Northern Hemisphere

Simplified representations of processes influencing forest biomass in Earth system models (ESMs) contribute to large uncertainty in projections. We evaluate forest biomass from eight ESMs outputs archived in the Coupled Model Intercomparison Project Phase 5 (CMIP5) using the biomass data synthesized from radar remote sensing and ground-based observations across northern extratropical latitudes. ESMs exhibit large biases in the forest distribution, forest fraction, and mass of carbon pools that contribute to uncertainty in forest total biomass (biases range from −20 Pg C to 135 Pg C). Forest total biomass is primarily positively correlated with precipitation variations, with surface temperature becoming equally important at higher latitudes, in both simulations and observations. Relatively small differences in forest biomass between the pre-industrial period and the contemporary period indicate uncertainties in forest biomass were introduced in the pre-industrial model equilibration (spin-up), suggesting parametric or structural model differences are a larger source of uncertainty than differences in transient responses. Our findings emphasize the importance of improved (1) models of carbon allocation to biomass compartments, (2) distribution of vegetation types in models, and (3) reproduction of pre-industrial vegetation conditions, in order to reduce the uncertainty in forest biomass simulated by ESMs.

Detecting early warning signals of tree mortality in boreal North America using multiscale satellite data

Increasing tree mortality from global change drivers such as drought and biotic infestations is a widespread phenomenon, including in the boreal zone where climate changes and feedbacks to the Earth system are relatively large. Despite the importance for science and management communities, our ability to forecast tree mortality at landscape to continental scales is limited. However, two independent information streams have the potential to inform and improve mortality forecasts: repeat forest inventories and satellite remote sensing. Time series of tree-level growth patterns indicate that productivity declines and related temporal dynamics often precede mortality years to decades before death. Plot-level productivity, in turn, has been related to satellite-based indices such as the Normalized difference vegetation index (NDVI). Here we link these two data sources to show that early warning signals of mortality are evident in several NDVI-based metrics up to 24 years before death. We focus on two repeat forest inventories and three NDVI products across western boreal North America where productivity and mortality dynamics are influenced by periodic drought. These data sources capture a range of forest conditions and spatial resolution to highlight the sensitivity and limitations of our approach. Overall, results indicate potential to use satellite NDVI for early warning signals of mortality. Relationships are broadly consistent across inventories, species, and spatial resolutions, although the utility of coarse-scale imagery in the heterogeneous aspen parkland was limited. Longer-term NDVI data and annually remeasured sites with high mortality levels generate the strongest signals, although we still found robust relationships at sites remeasured at a typical 5 year frequency. The approach and relationships developed here can be used as a basis for improving forest mortality models and monitoring systems.