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.

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 and vegetation optical depth 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-20. Regions of prominent trends (from [Formula: see text] 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.

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.

Declining Amazon biomass due to deforestation and subsequent degradation losses exceeding gains.

In the Amazon, deforestation and climate change lead to increased vulnerability to forest degradation, threatening its existing carbon stocks and its capacity as a carbon sink. We use satellite L-Band Vegetation Optical Depth (L-VOD) data that provide an integrated (top-down) estimate of biomass carbon to track changes over 2011-2019. Because the spatial resolution of L-VOD is coarse (0.25°), it allows limited attribution of the observed changes. We therefore combined high-resolution annual maps of forest cover and disturbances with biomass maps to model carbon losses (bottom-up) from deforestation and degradation, and gains from regrowing secondary forests. We show an increase of deforestation and associated degradation losses since 2012 which greatly outweigh secondary forest gains. Degradation accounted for 40% of gross losses. After an increase in 2011, old-growth forests show a net loss of above-ground carbon between 2012 and 2019. The sum of component carbon fluxes in our model is consistent with the total biomass change from L-VOD of 1.3 Pg C over 2012-2019. Across nine Amazon countries, we found that while Brazil contains the majority of biomass stocks (64%), its losses from disturbances were disproportionately high (79% of gross losses). Our multi-source analysis provides a pessimistic assessment of the Amazon carbon balance and highlights the urgent need to stop the recent rise of deforestation and degradation, particularly in the Brazilian Amazon.

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 (EFmax), θcrit, and slope (S) between EF and SM. We find systematic differences of these three variables across biomes. Variation in θcrit, S, and EFmax 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.

Increasing and widespread vulnerability of intact tropical rainforests to repeated droughts.

Intact tropical rainforests have been exposed to severe droughts in recent decades, which may threaten their integrity, their ability to sequester carbon, and their capacity to provide shelter for biodiversity. However, their response to droughts remains uncertain due to limited high-quality, long-term observations covering extensive areas. Here, we examined how the upper canopy of intact tropical rainforests has responded to drought events globally and during the past 3 decades. By developing a long pantropical time series (1992 to 2018) of monthly radar satellite observations, we show that repeated droughts caused a sustained decline in radar signal in 93%, 84%, and 88% of intact tropical rainforests in the Americas, Africa, and Asia, respectively. Sudden decreases in radar signal were detected around the 1997-1998, 2005, 2010, and 2015 droughts in tropical Americas; 1999-2000, 2004-2005, 2010-2011, and 2015 droughts in tropical Africa; and 1997-1998, 2006, and 2015 droughts in tropical Asia. Rainforests showed similar low resistance (the ability to maintain predrought condition when drought occurs) to severe droughts across continents, but American rainforests consistently showed the lowest resilience (the ability to return to predrought condition after the drought event). Moreover, while the resistance of intact tropical rainforests to drought is decreasing, albeit weakly in tropical Africa and Asia, forest resilience has not increased significantly. Our results therefore suggest the capacity of intact rainforests to withstand future droughts is limited. This has negative implications for climate change mitigation through forest-based climate solutions and the associated pledges made by countries under the Paris Agreement.

Climatic and biotic factors influencing regional declines and recovery of tropical forest biomass from the 2015/16 El Niño.

The 2015/16 El Niño brought severe drought and record-breaking temperatures in the tropics. Here, using satellite-based L-band microwave vegetation optical depth, we mapped changes of above-ground biomass (AGB) during the drought and in subsequent years up to 2019. Over more than 60% of drought-affected intact forests, AGB reduced during the drought, except in the wettest part of the central Amazon, where it declined 1 y later. By the end of 2019, only 40% of AGB reduced intact forests had fully recovered to the predrought level. Using random-forest models, we found that the magnitude of AGB losses during the drought was mainly associated with regionally distinct patterns of soil water deficits and soil clay content. For the AGB recovery, we found strong influences of AGB losses during the drought and of [Formula: see text]. [Formula: see text] is a parameter related to canopy structure and is defined as the ratio of two relative height (RH) metrics of Geoscience Laser Altimeter System (GLAS) waveform data-RH25 (25% energy return height) and RH100 (100% energy return height; i.e., top canopy height). A high [Formula: see text] may reflect forests with a tall understory, thick and closed canopy, and/or without degradation. Such forests with a high [Formula: see text] ([Formula: see text] ≥ 0.3) appear to have a stronger capacity to recover than low-[Formula: see text] ones. Our results highlight the importance of forest structure when predicting the consequences of future drought stress in the tropics.

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 CO2 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 by 1.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.

Recent divergence in the contributions of tropical and boreal forests to the terrestrial carbon sink.

Anthropogenic land use and land cover changes (LULCC) have a large impact on the global terrestrial carbon sink, but this effect is not well characterized according to biogeographical region. Here, using state-of-the-art Earth observation data and a dynamic global vegetation model, we estimate the impact of LULCC on the contribution of biomes to the terrestrial carbon sink between 1992 and 2015. Tropical and boreal forests contributed equally, and with the largest share of the mean global terrestrial carbon sink. CO2 fertilization was found to be the main driver increasing the terrestrial carbon sink from 1992 to 2015, but the net effect of all drivers (CO2 fertilization and nitrogen deposition, LULCC and meteorological forcing) caused a reduction and an increase, respectively, in the terrestrial carbon sink for tropical and boreal forests. These diverging trends were not observed when applying a conventional LULCC dataset, but were also evident in satellite passive microwave estimates of aboveground biomass. These datasets thereby converge on the conclusion that LULCC have had a greater impact on tropical forests than previously estimated, causing an increase and decrease of the contributions of boreal and tropical forests, respectively, to the growing terrestrial carbon sink.

Forest management in southern China generates short term extensive carbon sequestration.

Land use policies have turned southern China into one of the most intensively managed forest regions in the world, with actions maximizing forest cover on soils with marginal agricultural potential while concurrently increasing livelihoods and mitigating climate change. Based on satellite observations, here we show that diverse land use changes in southern China have increased standing aboveground carbon stocks by 0.11 ± 0.05 Pg C y-1 during 2002-2017. Most of this regional carbon sink was contributed by newly established forests (32%), while forests already existing contributed 24%. Forest growth in harvested forest areas contributed 16% and non-forest areas contributed 28% to the carbon sink, while timber harvest was tripled. Soil moisture declined significantly in 8% of the area. We demonstrate that land management in southern China has been removing an amount of carbon equivalent to 33% of regional fossil CO2 emissions during the last 6 years, but forest growth saturation, land competition for food production and soil-water depletion challenge the longevity of this carbon sink service.

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.

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 stocks1. 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 CO2 growth-rate anomalies (Pearson's r = 0.86), highlighting the pivotal role of tropical AGC in the global carbon budget.

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 CO2 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-1 (inversions) and by 1.0 PgC yr-1 (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.