Microbial carbon use efficiency promotes global soil carbon storage.

Soils store more carbon than other terrestrial ecosystems1,2. How soil organic carbon (SOC) forms and persists remains uncertain1,3, which makes it challenging to understand how it will respond to climatic change3,4. It has been suggested that soil microorganisms play an important role in SOC formation, preservation and loss5-7. Although microorganisms affect the accumulation and loss of soil organic matter through many pathways4,6,8-11, microbial carbon use efficiency (CUE) is an integrative metric that can capture the balance of these processes12,13. Although CUE has the potential to act as a predictor of variation in SOC storage, the role of CUE in SOC persistence remains unresolved7,14,15. Here we examine the relationship between CUE and the preservation of SOC, and interactions with climate, vegetation and edaphic properties, using a combination of global-scale datasets, a microbial-process explicit model, data assimilation, deep learning and meta-analysis. We find that CUE is at least four times as important as other evaluated factors, such as carbon input, decomposition or vertical transport, in determining SOC storage and its spatial variation across the globe. In addition, CUE shows a positive correlation with SOC content. Our findings point to microbial CUE as a major determinant of global SOC storage. Understanding the microbial processes underlying CUE and their environmental dependence may help the prediction of SOC feedback to a changing climate.

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.

Deep learning and process understanding for data-driven Earth system science.

Machine learning approaches are increasingly used to extract patterns and insights from the ever-increasing stream of geospatial data, but current approaches may not be optimal when system behaviour is dominated by spatial or temporal context. Here, rather than amending classical machine learning, we argue that these contextual cues should be used as part of deep learning (an approach that is able to extract spatio-temporal features automatically) to gain further process understanding of Earth system science problems, improving the predictive ability of seasonal forecasting and modelling of long-range spatial connections across multiple timescales, for example. The next step will be a hybrid modelling approach, coupling physical process models with the versatility of data-driven machine learning.

Global patterns of tree wood density.

Wood density is a fundamental property related to tree biomechanics and hydraulic function while playing a crucial role in assessing vegetation carbon stocks by linking volumetric retrieval and a mass estimate. This study provides a high-resolution map of the global distribution of tree wood density at the 0.01° (~1 km) spatial resolution, derived from four decision trees machine learning models using a global database of 28,822 tree-level wood density measurements. An ensemble of four top-performing models combined with eight cross-validation strategies shows great consistency, providing wood density patterns with pronounced spatial heterogeneity. The global pattern shows lower wood density values in northern and northwestern Europe, Canadian forest regions and slightly higher values in Siberia forests, western United States, and southern China. In contrast, tropical regions, especially wet tropical areas, exhibit high wood density. Climatic predictors explain 49%-63% of spatial variations, followed by vegetation characteristics (25%-31%) and edaphic properties (11%-16%). Notably, leaf type (evergreen vs. deciduous) and leaf habit type (broadleaved vs. needleleaved) are the most dominant individual features among all selected predictive covariates. Wood density tends to be higher for angiosperm broadleaf trees compared to gymnosperm needleleaf trees, particularly for evergreen species. The distributions of wood density categorized by leaf types and leaf habit types have good agreement with the features observed in wood density measurements. This global map quantifying wood density distribution can help improve accurate predictions of forest carbon stocks, providing deeper insights into ecosystem functioning and carbon cycling such as forest vulnerability to hydraulic and thermal stresses in the context of future climate change.

Unexpectedly large impact of forest management and grazing on global vegetation biomass.

Carbon stocks in vegetation have a key role in the climate system. However, the magnitude, patterns and uncertainties of carbon stocks and the effect of land use on the stocks remain poorly quantified. Here we show, using state-of-the-art datasets, that vegetation currently stores around 450 petagrams of carbon. In the hypothetical absence of land use, potential vegetation would store around 916 petagrams of carbon, under current climate conditions. This difference highlights the massive effect of land use on biomass stocks. Deforestation and other land-cover changes are responsible for 53-58% of the difference between current and potential biomass stocks. Land management effects (the biomass stock changes induced by land use within the same land cover) contribute 42-47%, but have been underestimated in the literature. Therefore, avoiding deforestation is necessary but not sufficient for mitigation of climate change. Our results imply that trade-offs exist between conserving carbon stocks on managed land and raising the contribution of biomass to raw material and energy supply for the mitigation of climate change. Efforts to raise biomass stocks are currently verifiable only in temperate forests, where their potential is limited. By contrast, large uncertainties hinder verification in the tropical forest, where the largest potential is located, pointing to challenges for the upcoming stocktaking exercises under the Paris agreement.

The detection and attribution of extreme reductions in vegetation growth across the global land surface.

Negative extreme anomalies in vegetation growth (NEGs) usually indicate severely impaired ecosystem services. These NEGs can result from diverse natural and anthropogenic causes, especially climate extremes (CEs). However, the relationship between NEGs and many types of CEs remains largely unknown at regional and global scales. Here, with satellite-derived vegetation index data and supporting tree-ring chronologies, we identify periods of NEGs from 1981 to 2015 across the global land surface. We find 70% of these NEGs are attributable to five types of CEs and their combinations, with compound CEs generally more detrimental than individual ones. More importantly, we find that dominant CEs for NEGs vary by biome and region. Specifically, cold and/or wet extremes dominate NEGs in temperate mountains and high latitudes, whereas soil drought and related compound extremes are primarily responsible for NEGs in wet tropical, arid and semi-arid regions. Key characteristics (e.g., the frequency, intensity and duration of CEs, and the vulnerability of vegetation) that determine the dominance of CEs are also region- and biome-dependent. For example, in the wet tropics, dominant individual CEs have both higher intensity and longer duration than non-dominant ones. However, in the dry tropics and some temperate regions, a longer CE duration is more important than higher intensity. Our work provides the first global accounting of the attribution of NEGs to diverse climatic extremes. Our analysis has important implications for developing climate-specific disaster prevention and mitigation plans among different regions of the globe in a changing climate.

Human and climate global-scale imprint on sediment transfer during the Holocene.

Accelerated soil erosion has become a pervasive feature on landscapes around the world and is recognized to have substantial implications for land productivity, downstream water quality, and biogeochemical cycles. However, the scarcity of global syntheses that consider long-term processes has limited our understanding of the timing, the amplitude, and the extent of soil erosion over millennial time scales. As such, we lack the ability to make predictions about the responses of soil erosion to long-term climate and land cover changes. Here, we reconstruct sedimentation rates for 632 lakes based on chronologies constrained by 3,980 calibrated 14C ages to assess the relative changes in lake-watershed erosion rates over the last 12,000 y. Estimated soil erosion dynamics were then complemented with land cover reconstructions inferred from 43,669 pollen samples and with climate time series from the Max Planck Institute Earth System Model. Our results show that a significant portion of the Earth surface shifted to human-driven soil erosion rate already 4,000 y ago. In particular, inferred soil erosion rates increased in 35% of the watersheds, and most of these sites showed a decrease in the proportion of arboreal pollen, which would be expected with land clearance. Further analysis revealed that land cover change was the main driver of inferred soil erosion in 70% of all studied watersheds. This study suggests that soil erosion has been altering terrestrial and aquatic ecosystems for millennia, leading to carbon (C) losses that could have ultimately induced feedbacks on the climate system.

Global sensitivities of forest carbon changes to environmental conditions.

The responses of forest carbon dynamics to fluctuations in environmental conditions at a global scale remain elusive. Despite the understanding that favourable environmental conditions promote forest growth, these responses have been challenging to observe across different ecosystems and climate gradients. Based on a global annual time series of aboveground biomass (AGB) estimated from radar satellites between 1992 and 2018, we present forest carbon changes and provide insights on their sensitivities to environmental conditions across scales. Our findings indicate differences in forest carbon changes across AGB classes, with regions with carbon stocks of 50-125 MgC ha-1 depict the highest forest carbon gains and losses, while regions with 125-150 MgC ha-1  have the lowest forest carbon gains and losses in absolute terms. Net forest carbon change estimates show that the arc-of-deforestation and the Congo Basin were the main hotspots of forest carbon loss, while a substantial part of European forest gained carbon during the last three decades. Furthermore, we observe that changes in forest carbon stocks were systematically positively correlated with changes in forest cover fraction. At the same time, it was not necessarily the case with other environmental variables, such as air temperature and water availability at the bivariate level. We also used a model attribution method to demonstrate that atmospheric conditions were the dominant control of forest carbon changes (56% of the total study area) followed by water-related (29% of the total study area) and vegetation (15% of the total study area) conditions. Regionally, we find evidence that carbon gains from long-term forest growth covary with long-term carbon sinks inferred from atmospheric inversions. Our results describe the contributions from the atmosphere, water-related and vegetation conditions to forest carbon changes and provide new insights into the underlying mechanisms of the coupling between forest growth and the global carbon cycle.

Ecosystem transpiration and evaporation: Insights from three water flux partitioning methods across FLUXNET sites.

We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants' water use for global water and carbon cycling in a changing climate.

Global covariation of carbon turnover times with climate in terrestrial ecosystems.

The response of the terrestrial carbon cycle to climate change is among the largest uncertainties affecting future climate change projections. The feedback between the terrestrial carbon cycle and climate is partly determined by changes in the turnover time of carbon in land ecosystems, which in turn is an ecosystem property that emerges from the interplay between climate, soil and vegetation type. Here we present a global, spatially explicit and observation-based assessment of whole-ecosystem carbon turnover times that combines new estimates of vegetation and soil organic carbon stocks and fluxes. We find that the overall mean global carbon turnover time is 23(+7)(-4) years (95 per cent confidence interval). On average, carbon resides in the vegetation and soil near the Equator for a shorter time than at latitudes north of 75° north (mean turnover times of 15 and 255 years, respectively). We identify a clear dependence of the turnover time on temperature, as expected from our present understanding of temperature controls on ecosystem dynamics. Surprisingly, our analysis also reveals a similarly strong association between turnover time and precipitation. Moreover, we find that the ecosystem carbon turnover times simulated by state-of-the-art coupled climate/carbon-cycle models vary widely and that numerical simulations, on average, tend to underestimate the global carbon turnover time by 36 per cent. The models show stronger spatial relationships with temperature than do observation-based estimates, but generally do not reproduce the strong relationships with precipitation and predict faster carbon turnover in many semi-arid regions. Our findings suggest that future climate/carbon-cycle feedbacks may depend more strongly on changes in the hydrological cycle than is expected at present and is considered in Earth system models.

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

Turnover concepts in state-of-the-art global vegetation models (GVMs) account for various processes, but are often highly simplified and may not include an adequate representation of the dominant processes that shape vegetation carbon turnover rates in real forest ecosystems at a large spatial scale. Here, we evaluate vegetation carbon turnover processes in GVMs participating in the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP, including HYBRID4, JeDi, JULES, LPJml, ORCHIDEE, SDGVM, and VISIT) using estimates of vegetation carbon turnover rate (k) derived from a combination of remote sensing based products of biomass and net primary production (NPP). We find that current model limitations lead to considerable biases in the simulated biomass and in k (severe underestimations by all models except JeDi and VISIT compared to observation-based average k), likely contributing to underestimation of positive feedbacks of the northern forest carbon balance to climate change caused by changes in forest mortality. A need for improved turnover concepts related to frost damage, drought, and insect outbreaks to better reproduce observation-based spatial patterns in k is identified. As direct frost damage effects on mortality are usually not accounted for in these GVMs, simulated relationships between k and winter length in boreal forests are not consistent between different regions and strongly biased compared to the observation-based relationships. Some models show a response of k to drought in temperate forests as a result of impacts of water availability on NPP, growth efficiency or carbon balance dependent mortality as well as soil or litter moisture effects on leaf turnover or fire. However, further direct drought effects such as carbon starvation (only in HYBRID4) or hydraulic failure are usually not taken into account by the investigated GVMs. While they are considered dominant large-scale mortality agents, mortality mechanisms related to insects and pathogens are not explicitly treated in these models.

Codominant water control on global interannual variability and trends in land surface phenology and greenness.

Identifying the relative importance of climatic and other environmental controls on the interannual variability and trends in global land surface phenology and greenness is challenging. Firstly, quantifications of land surface phenology and greenness dynamics are impaired by differences between satellite data sets and phenology detection methods. Secondly, dynamic global vegetation models (DGVMs) that can be used to diagnose controls still reveal structural limitations and contrasting sensitivities to environmental drivers. Thus, we assessed the performance of a new developed phenology module within the LPJmL (Lund-Potsdam-Jena managed Lands) DGVM with a comprehensive ensemble of three satellite data sets of vegetation greenness and ten phenology detection methods, thereby thoroughly accounting for observational uncertainties. The improved and tested model allows us quantifying the relative importance of environmental controls on interannual variability and trends of land surface phenology and greenness at regional and global scales. We found that start of growing season interannual variability and trends are in addition to cold temperature mainly controlled by incoming radiation and water availability in temperate and boreal forests. Warming-induced prolongations of the growing season in high latitudes are dampened by a limited availability of light. For peak greenness, interannual variability and trends are dominantly controlled by water availability and land-use and land-cover change (LULCC) in all regions. Stronger greening trends in boreal forests of Siberia than in North America are associated with a stronger increase in water availability from melting permafrost soils. Our findings emphasize that in addition to cold temperatures, water availability is a codominant control for start of growing season and peak greenness trends at the global scale.

Influence of physiological phenology on the seasonal pattern of ecosystem respiration in deciduous forests.

Understanding the environmental and biotic drivers of respiration at the ecosystem level is a prerequisite to further improve scenarios of the global carbon cycle. In this study we investigated the relevance of physiological phenology, defined as seasonal changes in plant physiological properties, for explaining the temporal dynamics of ecosystem respiration (RECO) in deciduous forests. Previous studies showed that empirical RECO models can be substantially improved by considering the biotic dependency of RECO on the short-term productivity (e.g., daily gross primary production, GPP) in addition to the well-known environmental controls of temperature and water availability. Here, we use a model-data integration approach to investigate the added value of physiological phenology, represented by the first temporal derivative of GPP, or alternatively of the fraction of absorbed photosynthetically active radiation, for modeling RECO at 19 deciduous broadleaved forests in the FLUXNET La Thuile database. The new data-oriented semiempirical model leads to an 8% decrease in root mean square error (RMSE) and a 6% increase in the modeling efficiency (EF) of modeled RECO when compared to a version of the model that does not consider the physiological phenology. The reduction of the model-observation bias occurred mainly at the monthly time scale, and in spring and summer, while a smaller reduction was observed at the annual time scale. The proposed approach did not improve the model performance at several sites, and we identified as potential causes the plant canopy heterogeneity and the use of air temperature as a driver of ecosystem respiration instead of soil temperature. However, in the majority of sites the model-error remained unchanged regardless of the driving temperature. Overall, our results point toward the potential for improving current approaches for modeling RECO in deciduous forests by including the phenological cycle of the canopy.

Characterizing the Response of Vegetation Cover to Water Limitation in Africa Using Geostationary Satellites.

Hydrological interactions between vegetation, soil, and topography are complex, and heterogeneous in semi-arid landscapes. This along with data scarcity poses challenges for large-scale modeling of vegetation-water interactions. Here, we exploit metrics derived from daily Meteosat data over Africa at ca. 5 km spatial resolution for ecohydrological analysis. Their spatial patterns are based on Fractional Vegetation Cover (FVC) time series and emphasize limiting conditions of the seasonal wet to dry transition: the minimum and maximum FVC of temporal record, the FVC decay rate and the FVC integral over the decay period. We investigate the relevance of these metrics for large scale ecohydrological studies by assessing their co-variation with soil moisture, and with topographic, soil, and vegetation factors. Consistent with our initial hypothesis, FVC minimum and maximum increase with soil moisture, while the FVC integral and decay rate peak at intermediate soil moisture. We find evidence for the relevance of topographic moisture variations in arid regions, which, counter-intuitively, is detectable in the maximum but not in the minimum FVC. We find no clear evidence for wide-spread occurrence of the "inverse texture effect" on FVC. The FVC integral over the decay period correlates with independent data sets of plant water storage capacity or rooting depth while correlations increase with aridity. In arid regions, the FVC decay rate decreases with canopy height and tree cover fraction as expected for ecosystems with a more conservative water-use strategy. Thus, our observation-based products have large potential for better understanding complex vegetation-water interactions from regional to continental scales.

Observation-based assessment of secondary water effects on seasonal vegetation decay across Africa.

Local studies and modeling experiments suggest that shallow groundwater and lateral redistribution of soil moisture, together with soil properties, can be highly important secondary water sources for vegetation in water-limited ecosystems. However, there is a lack of observation-based studies of these terrain-associated secondary water effects on vegetation over large spatial domains. Here, we quantify the role of terrain properties on the spatial variations of dry season vegetation decay rate across Africa obtained from geostationary satellite acquisitions to assess the large-scale relevance of secondary water effects. We use machine learning based attribution to identify where and under which conditions terrain properties related to topography, water table depth, and soil hydraulic properties influence the rate of vegetation decay. Over the study domain, the machine learning model attributes about one-third of the spatial variations of vegetation decay rates to terrain properties, which is roughly equally split between direct terrain effects and interaction effects with climate and vegetation variables. The importance of secondary water effects increases with increasing topographic variability, shallower groundwater levels, and the propensity to capillary rise given by soil properties. In regions with favorable terrain properties, more than 60% of the variations in the decay rate of vegetation are attributed to terrain properties, highlighting the importance of secondary water effects on vegetation in Africa. Our findings provide an empirical assessment of the importance of local-scale secondary water effects on vegetation over Africa and help to improve hydrological and vegetation models for the challenge of bridging processes across spatial scales.

A climate risk analysis of Earth's forests in the 21st century.

Earth's forests harbor extensive biodiversity and are currently a major carbon sink. Forest conservation and restoration can help mitigate climate change; however, climate change could fundamentally imperil forests in many regions and undermine their ability to provide such mitigation. The extent of climate risks facing forests has not been synthesized globally nor have different approaches to quantifying forest climate risks been systematically compared. We combine outputs from multiple mechanistic and empirical approaches to modeling carbon, biodiversity, and disturbance risks to conduct a synthetic climate risk analysis for Earth's forests in the 21st century. Despite large uncertainty in most regions we find that some forests are consistently at higher risk, including southern boreal forests and those in western North America and parts of the Amazon.

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.

Memory effects of climate and vegetation affecting net ecosystem CO2 fluxes in global forests.

Forests play a crucial role in the global carbon (C) cycle by storing and sequestering a substantial amount of C in the terrestrial biosphere. Due to temporal dynamics in climate and vegetation activity, there are significant regional variations in carbon dioxide (CO2) fluxes between the biosphere and atmosphere in forests that are affecting the global C cycle. Current forest CO2 flux dynamics are controlled by instantaneous climate, soil, and vegetation conditions, which carry legacy effects from disturbances and extreme climate events. Our level of understanding from the legacies of these processes on net CO2 fluxes is still limited due to their complexities and their long-term effects. Here, we combined remote sensing, climate, and eddy-covariance flux data to study net ecosystem CO2 exchange (NEE) at 185 forest sites globally. Instead of commonly used non-dynamic statistical methods, we employed a type of recurrent neural network (RNN), called Long Short-Term Memory network (LSTM) that captures information from the vegetation and climate's temporal dynamics. The resulting data-driven model integrates interannual and seasonal variations of climate and vegetation by using Landsat and climate data at each site. The presented LSTM algorithm was able to effectively describe the overall seasonal variability (Nash-Sutcliffe efficiency, NSE = 0.66) and across-site (NSE = 0.42) variations in NEE, while it had less success in predicting specific seasonal and interannual anomalies (NSE = 0.07). This analysis demonstrated that an LSTM approach with embedded climate and vegetation memory effects outperformed a non-dynamic statistical model (i.e. Random Forest) for estimating NEE. Additionally, it is shown that the vegetation mean seasonal cycle embeds most of the information content to realistically explain the spatial and seasonal variations in NEE. These findings show the relevance of capturing memory effects from both climate and vegetation in quantifying spatio-temporal variations in forest NEE.