Quantifying the resilience of vegetated ecosystems is key to constraining both present-day and future global impacts of anthropogenic climate change. Here we apply both empirical and theoretical resilience metrics to remotely-sensed vegetation data in order to examine the role of water availability and variability in controlling vegetation resilience at the global scale. We find a concise global relationship where vegetation resilience is greater in regions with higher water availability. We also reveal that resilience is lower in regions with more pronounced inter-annual precipitation variability, but find less concise relationships between vegetation resilience and intra-annual precipitation variability. Our results thus imply that the resilience of vegetation responds differently to water deficits at varying time scales. In view of projected increases in precipitation variability, our findings highlight the risk of ecosystem degradation under ongoing climate change.

The character and health of ecosystems worldwide is tightly coupled to changes in Earth’s climate. Theory suggests that ecosystem resilience—the ability of ecosystems to resist and recover from external shocks such as droughts and fires—can be inferred from their natural variability. Here, we quantify vegetation resilience globally with complementary metrics based on two independent long-term satellite records. We first empirically confirm that the recovery rates from large perturbations can be closely approximated from internal vegetation variability across vegetation types and climate zones. On the basis of this empirical relationship, we quantify vegetation resilience continuously and globally from 1992 to 2017. Long-term vegetation resilience trends are spatially heterogeneous, with overall increasing resilience in the tropics and decreasing resilience at higher latitudes. Shorter-term trends, however, reveal a marked shift towards a global decline in vegetation resilience since the early 2000s, particularly in the equatorial rainforest belt.

Snow meltwaters account for most of the yearly water budgets of many catchments in High Mountain Asia (HMA). We examine trends in snow water equivalent (SWE) using passive microwave data (1987 to 2009). We find an overall decrease in SWE in HMA, despite regions of increased SWE in the Pamir, Kunlun Shan, Eastern Himalaya, and Eastern Tien Shan. Although the average decline in annual SWE across HMA (contributing area, 2641 × 10³ km²) is low (average, −0.3%), annual SWE losses conceal distinct seasonal and spatial heterogeneities across the study region. For example, the Tien Shan has seen both strong increases in winter SWE and sharp declines in spring and summer SWE. In the majority of catchments, the most negative SWE trends are found in mid-elevation zones, which often correspond to the regions of highest snow-water storage and are somewhat distinct from glaciated areas. Negative changes in SWE storage in these mid-elevation zones have strong implications for downstream water availability.

Potential climate tipping points pose a growing risk for societies, and policy is calling for improved anticipation of them. Satellite remote sensing can play a unique role in identifying and anticipating tipping phenomena across scales. Where satellite records are too short for temporal early warning of tipping points, complementary spatial indicators can leverage the exceptional spatial-temporal coverage of remotely sensed data to detect changing resilience of vulnerable systems. Combining Earth observation with Earth system models can improve process-based understanding of tipping points, their interactions, and potential tipping cascades. Such fine-resolution sensing can support climate tipping point risk management across scales.

Concerns have been raised that the resilience of vegetated ecosystems may be negatively impacted by ongoing anthropogenic climate and land-use change at the global scale. Several recent studies present global vegetation resilience trends based on satellite data using diverse methodological set-ups. Here, upon a systematic comparison of data sets, spatial and temporal pre-processing, and resilience estimation methods, we propose a methodology that avoids different biases present in previous results. Nevertheless, we find that resilience estimation using optical satellite vegetation data is broadly problematic in dense tropical and high-latitude boreal forests, regardless of the vegetation index chosen. However, for wide parts of the mid-latitudes—especially with low biomass density—resilience can be reliably estimated using several optical vegetation indices. We infer a spatially consistent global pattern of resilience gain and loss across vegetation indices, with more regions facing declining resilience, especially in Africa, Australia and central Asia.

Many widely used observational data sets are comprised of several overlapping instrument records. While data inter-calibration techniques often yield continuous and reliable data for trend analysis, less attention is generally paid to maintaining higher-order statistics such as variance and autocorrelation. A growing body of work uses these metrics to quantify the stability or resilience of a system under study and potentially to anticipate an approaching critical transition in the system. Exploring the degree to which changes in resilience indicators such as the variance or autocorrelation can be attributed to non-stationary characteristics of the measurement process – rather than actual changes in the dynamical properties of the system – is important in this context. In this work we use both synthetic and empirical data to explore how changes in the noise structure of a data set are propagated into the commonly used resilience metrics lag-one autocorrelation and variance. We focus on examples from remotely sensed vegetation indicators such as vegetation optical depth and the normalized difference vegetation index from different satellite sources. We find that time series resulting from mixing signals from sensors with varied uncertainties and covering overlapping time spans can lead to biases in inferred resilience changes. These biases are typically more pronounced when resilience metrics are aggregated (for example, by land-cover type or region), whereas estimates for individual time series remain reliable at reasonable sensor signal-to-noise ratios. Our work provides guidelines for the treatment and aggregation of multi-instrument data in studies of critical transitions and resilience.