Anthropogenic disturbance exacerbates resilience loss in the Amazon rainforests

Uncovering the mechanisms that lead to Amazon forest resilience variations is crucial to predict the impact of future climatic and anthropogenic disturbances. Here, we apply a previously used empirical resilience metrics, lag-1 month temporal autocorrelation (TAC), to vegetation optical depth data in C-band (a good proxy of the whole canopy water content) in order to explore how forest resilience variations are impacted by human disturbances and environmental drivers in the Brazilian Amazon. We found that human disturbances significantly increase the risk of critical transitions, and that the median TAC value is ~2.4 times higher in human-disturbed forests than that in intact forests, suggesting a much lower resilience in disturbed forests. Additionally, human-disturbed forests are less resilient to land surface heat stress and atmospheric water stress than intact forests. Among human-disturbed forests, forests with a more closed and thicker canopy structure, which is linked to a higher forest cover and a lower disturbance fraction, are comparably more resilient. These results further emphasize the urgent need to limit deforestation and degradation through policy intervention to maintain the resilience of the Amazon rainforests.

Sub-continental-scale carbon stocks of individual trees in African drylands

The distribution of dryland trees and their density, cover, size, mass and carbon content are not well known at sub-continental to continental scales1-14. This information is important for ecological protection, carbon accounting, climate mitigation and restoration efforts of dryland ecosystems15-18. We assessed more than 9.9 billion trees derived from more than 300,000 satellite images, covering semi-arid sub-Saharan Africa north of the Equator. We attributed wood, foliage and root carbon to every tree in the 0-1,000 mm year-1 rainfall zone by coupling field data19, machine learning20-22, satellite data and high-performance computing. Average carbon stocks of individual trees ranged from 0.54 Mg C ha-1 and 63 kg C tree-1 in the arid zone to 3.7 Mg C ha-1 and 98 kg tree-1 in the sub-humid zone. Overall, we estimated the total carbon for our study area to be 0.84 (±19.8%) Pg C. Comparisons with 14 previous TRENDY numerical simulation studies23 for our area found that the density and carbon stocks of scattered trees have been underestimated by three models and overestimated by 11 models, respectively. This benchmarking can help understand the carbon cycle and address concerns about land degradation24-29. We make available a linked database of wood mass, foliage mass, root mass and carbon stock of each tree for scientists, policymakers, dryland-restoration practitioners and farmers, who can use it to estimate farmland tree carbon stocks from tablets or laptops.

Soil and vegetation water content identify the main terrestrial ecosystem changes

Environmental change is a consequence of many interrelated factors. How vegetation responds to natural and human activity still needs to be well established, quantified, and understood. Recent satellite missions providing hydrologic and ecological indicators enable better monitoring of Earth system changes, yet there is no automatic way to address this issue directly from observations. Here, we develop an observation-based methodology to capture evidence of changes in global terrestrial ecosystems and attribute these changes to natural or anthropogenic activity. We use the longest time record of global microwave L-band soil moisture (SM) and vegetation optical depth (VOD) as satellite data and build spatially-explicit maps of change in soil and vegetation water content and biomass reflecting large ecosystem changes during the last decade, 2010-2020. Regions of prominent trends (from -8% to 9% per year) are observed, especially in humid and semi-arid climates. We further combine such trends with land cover change maps, vegetation greenness, and precipitation variability to assess their relationship with major documented ecosystem changes. Several regions emerge from our results. They cluster changes according to human activity drivers, including deforestation (Amazon, Central Africa) and wildfires (East Australia), artificial reforestation (South-East China), abandonment of farm fields (Central Russia), and climate shifts related to changes in precipitation variability (East Africa, North America, and Central Argentina). Using the high sensitivity of soil and vegetation water content to ecosystem changes, microwave satellite observations enable us to quantify and attribute global vegetation responses to climate or anthropogenic activities as a direct measure of environmental changes and the mechanisms driving them.

Critical soil moisture thresholds of plant water stress in terrestrial ecosystems

Plant water stress occurs at the point when soil moisture (SM) limits transpiration, defining a critical SM threshold (θ_crit). Knowledge of the spatial distribution of θ_crit is crucial for future projections of climate and water resources. Here, we use global eddy covariance observations to quantify θ_crit and evaporative fraction (EF) regimes. Three canonical variables describe how EF is controlled by SM: the maximum EF (EF_max), θ_crit, and slope (S) between EF and SM. We find systematic differences of these three variables across biomes. Variation in θ_crit, S, and EF_max is mostly explained by soil texture, vapor pressure deficit, and precipitation, respectively, as well as vegetation structure. Dryland ecosystems tend to operate at low θ_crit and show adaptation to water deficits. The negative relationship between θ_crit and S indicates that dryland ecosystems minimize θ_crit through mechanisms of sustained SM extraction and transport by xylem. Our results further suggest an optimal adaptation of local EF-SM response that maximizes growing-season evapotranspiration and photosynthesis.

Exploring complex water stress-gross primary production relationships: Impact of climatic drivers, main effects, and interactive effects

The dominance of vapor pressure deficit (VPD) and soil water content (SWC) for plant water stress is still under debate. These two variables are strongly coupled and influenced by climatic drivers. The impacts of climatic drivers on the relationships between gross primary production (GPP) and water stress from VPD/SWC and the interaction between VPD and SWC are not fully understood. Here, applying statistical methods and extreme gradient boosting models-Shapley additive explanations framework to eddy-covariance observations from the global FLUXNET2015 data set, we found that the VPD-GPP relationship was strongly influenced by climatic interactions and that VPD was more important for plant water stress than SWC across most plant functional types when we removed the effect of main climatic drivers, e.g. air temperature, incoming shortwave radiation and wind speed. However, we found no evidence for a significant influence of elevated CO₂ on stress alleviation, possibly because of the short duration of the records (approximately one decade). Additionally, the interactive effect between VPD and SWC differed from their individual effect. When SWC was high, the SHAP interaction value of SWC and VPD on GPP was decreased with increasing VPD, but when SWC was low, the trend was the opposite. Additionally, we revealed a threshold effect for VPD stress on GPP loss; above the threshold value, the stress on GPP was flattened off. Our results have important implications for independently identifying VPD and SWC limitations on plant productivity, which is meaningful for capturing the magnitude of ecosystem responses to water stress in dynamic global vegetation models.

Uncovering the critical soil moisture thresholds of plant water stress for European ecosystems

Understanding the critical soil moisture (SM) threshold (θ_crit ) of plant water stress and land surface energy partitioning is a basis to evaluate drought impacts and improve models for predicting future ecosystem condition and climate. Quantifying the θ_crit across biomes and climates is challenging because observations of surface energy fluxes and SM remain sparse. Here, we used the latest database of eddy covariance measurements to estimate θ_crit across Europe by evaluating evaporative fraction (EF)-SM relationships and investigating the covariance between vapor pressure deficit (VPD) and gross primary production (GPP) during SM dry-down periods. We found that the θ_crit and soil matric potential threshold in Europe are 16.5% and -0.7 MPa, respectively. Surface energy partitioning characteristics varied among different vegetation types; EF in savannas had the highest sensitivities to SM in water-limited stage, and the lowest in forests. The sign of the covariance between daily VPD and GPP consistently changed from positive to negative during dry-down across all sites when EF shifted from relatively high to low values. This sign of the covariance changed after longer period of SM decline in forests than in grasslands and savannas. Estimated θ_crit from the VPD-GPP covariance method match well with the EF-SM method, showing this covariance method can be used to detect the θ_crit . We further found that soil texture dominates the spatial variability of θ_crit while shortwave radiation and VPD are the major drivers in determining the spatial pattern of EF sensitivities. Our results highlight for the first time that the sign change of the covariance between daily VPD and GPP can be used as an indicator of how ecosystems transition from energy to SM limitation. We also characterized the corresponding θ_crit and its drivers across diverse ecosystems in Europe, an essential variable to improve the representation of water stress in land surface models.

Tropical tall forests are more sensitive and vulnerable to drought than short forests

Our limited understanding of the impacts of drought on tropical forests significantly impedes our ability in accurately predicting the impacts of climate change on this biome. Here, we investigated the impact of drought on the dynamics of forest canopies with different heights using time-series records of remotely sensed Ku-band vegetation optical depth (Ku-VOD), a proxy of top-canopy foliar mass and water content, and separated the signal of Ku-VOD changes into drought-induced reductions and subsequent non-drought gains. Both drought-induced reductions and non-drought increases in Ku-VOD varied significantly with canopy height. Taller tropical forests experienced greater relative Ku-VOD reductions during drought and larger non-drought increases than shorter forests, but the net effect of drought was more negative in the taller forests. Meta-analysis of in situ hydraulic traits supports the hypothesis that taller tropical forests are more vulnerable to drought stress due to smaller xylem-transport safety margins. Additionally, Ku-VOD of taller forests showed larger reductions due to increased atmospheric dryness, as assessed by vapor pressure deficit, and showed larger gains in response to enhanced water supply than shorter forests. Including the height-dependent variation of hydraulic transport in ecosystem models will improve the simulated response of tropical forests to drought.

Asymmetric responses of ecosystem productivity to rainfall anomalies vary inversely with mean annual rainfall over the conterminous United States

The CONterminous United States (CONUS) presents a large range of climate conditions and biomes where terrestrial primary productivity and its inter-annual variability are controlled regionally by rainfall and/or temperature. Here, the response of ecosystem productivity to those climate variables was investigated across different biomes from 2010 to 2018 using three climate datasets of precipitation, air temperature or drought severity, combined with several proxies of ecosystem productivity: a remote sensing product of aboveground biomass, an net primary productivity (NPP) remote sensing product, an NPP model-based product and four gross primary productivity products. We used an asymmetry index (AI) where positive AI indicates a greater increase of ecosystem productivity in wet years compared to the decline in dry years, and negative AI indicates a greater decline of ecosystem productivity in dry years compared to the increase in wet years. We found consistent spatial patterns of AI across the CONUS for the different products, with negative asymmetries over the Great Plains and positive asymmetries over the southwestern CONUS. Shrubs and, to a lesser extent, evergreen forests show a persistent positive asymmetry, whilst (natural) grasslands appear to have transitioned from positive to negative anomalies during the last decade. The general tendency of dominant negative asymmetry response for ecosystem productivity across the CONUS appears to be influenced by the negative asymmetry of precipitation anomalies. AI was found to be a function of mean rainfall: more positive AIs were found in dry areas where plants are adapted to drought and take advantage of rainfall pulses, and more negative AIs were found in wet areas, with a threshold delineating the two regimes corresponding to a mean annual rainfall of 200-400 mm/year.

Tropical forests did not recover from the strong 2015-2016 El Niño event

Severe drought and extreme heat associated with the 2015-2016 El Niño event have led to large carbon emissions from the tropical vegetation to the atmosphere. With the return to normal climatic conditions in 2017, tropical forest aboveground carbon (AGC) stocks are expected to partly recover due to increased productivity, but the intensity and spatial distribution of this recovery are unknown. We used low-frequency microwave satellite data (L-VOD) to feature precise monitoring of AGC changes and show that the AGC recovery of tropical ecosystems was slow and that by the end of 2017, AGC had not reached predrought levels of 2014. From 2014 to 2017, tropical AGC stocks decreased by1.3 1.2 1.5 Pg C due to persistent AGC losses in Africa (- 0.9 - 1.1 - 0.8 Pg C) and America (- 0.5 - 0.6 - 0.4 Pg C). Pantropically, drylands recovered their carbon stocks to pre-El Niño levels, but African and American humid forests did not, suggesting carryover effects from enhanced forest mortality.

Satellite-observed pantropical carbon dynamics

Changes in terrestrial tropical carbon stocks have an important role in the global carbon budget. However, current observational tools do not allow accurate and large-scale monitoring of the spatial distribution and dynamics of carbon stocks¹. Here, we used low-frequency L-band passive microwave observations to compute a direct and spatially explicit quantification of annual aboveground carbon (AGC) fluxes and show that the tropical net AGC budget was approximately in balance during 2010 to 2017, the net budget being composed of gross losses of -2.86 PgC yr^-1 offset by gross gains of -2.97 PgC yr^-1 between continents. Large interannual and spatial fluctuations of tropical AGC were quantified during the wet 2011 La Niña year and throughout the extreme dry and warm 2015-2016 El Niño episode. These interannual fluctuations, controlled predominantly by semiarid biomes, were shown to be closely related to independent global atmospheric CO₂ growth-rate anomalies (Pearson's r = 0.86), highlighting the pivotal role of tropical AGC in the global carbon budget.

Uncertainty in Soil Moisture Retrievals: an Ensemble Approach using SMOS L-Band Microwave Data

The uncertainty of surface soil moisture (SM) retrievals from satellite brightness temperature (TB) observations depends primarily on the choice of radiative transfer model (RTM) parameters, prior SM information and TB inputs. This paper studies the sensitivity of several established and experimental SM retrieval products from the Soil Moisture Ocean Salinity (SMOS) mission to these choices at 11 reference sites, located in 7 watersheds across the United States (US). Different RTM parameter sets cause large biases between retrievals. Whereas typical RTM parameter sets are calibrated for SM retrievals, it is shown that a parameter set carefully optimized for TB forward modeling can also be used for retrieving SM. It is also shown that the inclusion of dynamic prior SM estimates in a Bayesian retrieval scheme can strongly improve SM retrievals, regardless of the choice of RTM parameters. The second part of this paper evaluates the ensemble uncertainty metrics for SM retrievals obtained by propagating a wide range of RTM parameters through the RTM, and the relationship with time series metrics obtained by comparing SM retrievals with in situ data. As expected for bounded variables, the total spread in the ensemble SM retrievals is smallest for wet and dry SM values and highest for intermediate SM values. After removal of the strong long-term SM bias associated with the RTM parameter values for individual ensemble members, the remaining anomaly ensemble SM spread shows higher values when SM deviates further from its long-term mean SM. This reveals higher-order biases (e.g. differences in variances) in the retrieval error, which should be considered when characterizing retrieval error. The time-average anomaly ensemble SM spread of 0.037 m3/m3 approximates the actual time series unbiased root-mean-square-difference of 0.042 m3/m3 between ensemble mean retrievals and in situ data across the reference sites.

Changes in rainfall distribution promote woody foliage production in the Sahel

Dryland ecosystems comprise a balance between woody and herbaceous vegetation. Climate change impacts rainfall timing, which may alter the respective contributions of woody and herbaceous plants on the total vegetation production. Here, we apply 30 years of field-measured woody foliage and herbaceous mass from Senegal and document a faster increase in woody foliage mass (+17 kg ha-1 yr-1) as compared to herbaceous mass (+3 kg ha-1 yr-1). Annual rainfall trends were partitioned into core wet-season rains (+0.7 mm yr-1), supporting a weak but periodic (5-year cycles) increase in herbaceous mass, and early/late rains (+2.1 mm yr-1), explaining the strongly increased woody foliage mass. Satellite observations confirm these findings for the majority of the Sahel, with total herbaceous/woody foliage mass increases by 6%/20%. We conclude that the rainfall recovery in the Sahel does not benefit herbaceous vegetation to the same extent as woody vegetation, presumably favoured by increased early/late rains.

Impact of the 2015/2016 El Niño on the terrestrial carbon cycle constrained by bottom-up and top-down approaches

Evaluating the response of the land carbon sink to the anomalies in temperature and drought imposed by El Niño events provides insights into the present-day carbon cycle and its climate-driven variability. It is also a necessary step to build confidence in terrestrial ecosystems models' response to the warming and drying stresses expected in the future over many continents, and particularly in the tropics. Here we present an in-depth analysis of the response of the terrestrial carbon cycle to the 2015/2016 El Niño that imposed extreme warming and dry conditions in the tropics and other sensitive regions. First, we provide a synthesis of the spatio-temporal evolution of anomalies in net land–atmosphere CO₂ fluxes estimated by two in situ measurements based on atmospheric inversions and 16 land-surface models (LSMs) from TRENDYv6. Simulated changes in ecosystem productivity, decomposition rates and fire emissions are also investigated. Inversions and LSMs generally agree on the decrease and subsequent recovery of the land sink in response to the onset, peak and demise of El Niño conditions and point to the decreased strength of the land carbon sink: by 0.4–0.7 PgC yr⁻¹ (inversions) and by 1.0 PgC yr⁻¹ (LSMs) during 2015/2016. LSM simulations indicate that a decrease in productivity, rather than increase in respiration, dominated the net biome productivity anomalies in response to ENSO throughout the tropics, mainly associated with prolonged drought conditions.

This article is part of a discussion meeting issue ‘The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications’.

Coupling of ecosystem-scale plant water storage and leaf phenology observed by satellite

Plant water storage is fundamental to the functioning of terrestrial ecosystems by participating in plant metabolism, nutrient and sugar transport, and maintenance of the integrity of the hydraulic system of the plant. However, a global view of the size and dynamics of the water pools stored in plant tissues is still lacking. Here, we report global patterns of seasonal variations in ecosystem-scale plant water storage and their relationship with leaf phenology, based on space-borne measurements of L-band vegetation optical depth. We find that seasonal variations in plant water storage are highly synchronous with leaf phenology for the boreal and temperate forests, but asynchronous for the tropical woodlands, where the seasonal development of plant water storage lags behind leaf area by up to 180 days. Contrasting patterns of the time lag between plant water storage and terrestrial groundwater storage are also evident in these ecosystems. A comparison of the water cycle components in seasonally dry tropical woodlands highlights the buffering effect of plant water storage on the seasonal dynamics of water supply and demand. Our results offer insights into ecosystem-scale plant water relations globally and provide a basis for an improved parameterization of eco-hydrological and Earth system models.

Satellite passive microwaves reveal recent climate-induced carbon losses in African drylands

The African continent is facing one of the driest periods in the past three decades as well as continued deforestation. These disturbances threaten vegetation carbon (C) stocks and highlight the need for improved capabilities of monitoring large-scale aboveground carbon stock dynamics. Here we use a satellite dataset based on vegetation optical depth derived from low-frequency passive microwaves (L-VOD) to quantify annual aboveground biomass-carbon changes in sub-Saharan Africa between 2010 and 2016. L-VOD is shown not to saturate over densely vegetated areas. The overall net change in drylands (53% of the land area) was -0.05 petagrams of C per year (Pg C yr^-1) associated with drying trends, and a net change of -0.02 Pg C yr^-1 was observed in humid areas. These trends reflect a high inter-annual variability with a very dry year in 2015 (net change, -0.69 Pg C) with about half of the gross losses occurring in drylands. This study demonstrates, first, the applicability of L-VOD to monitor the dynamics of carbon loss and gain due to weather variations, and second, the importance of the highly dynamic and vulnerable carbon pool of dryland savannahs for the global carbon balance, despite the relatively low carbon stock per unit area.