Reclassifying historical disasters: From single to multi-hazards.

Multi-hazard events, characterized by the simultaneous, cascading, or cumulative occurrence of multiple natural hazards, pose a significant threat to human lives and assets. This is primarily due to the cumulative and cascading effects arising from the interplay of various natural hazards across space and time. However, their identification is challenging, which is attributable to the complex nature of natural hazard interactions and the limited availability of multi-hazard observations. This study presents an approach for identifying multi-hazard events during the past 123 years (1900-2023) using the EM-DAT global disaster database. Leveraging the 'associated hazard' information in EM-DAT, multi-hazard events are detected and assessed in relation to their frequency, impact on human lives and assets, and reporting trends. The interactions between various combinations of natural hazard pairs are explored, reclassifying them into four categories: preconditioned/triggering, multivariate, temporally compounding, and spatially compounding multi-hazard events. The results show, globally, approximately 19 % of the 16,535 disasters recorded in EM-DAT can be classified as multi-hazard events. However, the multi-hazard events recorded in EM-DAT are disproportionately responsible for nearly 59 % of the estimated global economic losses. Conversely, single hazard events resulted in higher fatalities compared to multi-hazard events. The largest proportion of multi-hazard events are associated with floods, storms, and earthquakes. Landslides emerge as the predominant secondary hazards within multi-hazard pairs, primarily triggered by floods, storms, and earthquakes, with the majority of multi-hazard events exhibiting preconditioned/triggering and multivariate characteristics. There is a higher prevalence of multi-hazard events in Asia and North America, whilst temporal overlaps of multiple hazards predominate in Europe. These results can be used to increase the integration of multi-hazard thinking in risk assessments, emergency management response plans and mitigation policies at both national and international levels.

Lexcube: Interactive Visualization of Large Earth System Data Cubes.

Many subsystems of Earth are constantly monitored in space and time and undergo continuous anthropogenic interventions. However, research into this transformation remains largely inaccessible to the public due to the complexity of the Big Data generated by models and Earth observation. To overcome this barrier, we present the Leipzig Explorer of Earth Data Cubes (lexcube.org), an interactive Earth data visualization that allows users to explore terabyte-scale datasets with minimal latency through space, time, variables, and model variants. With over 2800 users and 163,000 API requests since its public release in May 2022, lexcube.org is a novel interactive data cube visualization that embraces the concept of "data cubes," enabling an equal treatment of space and time. We expect this development to be particularly relevant for the emerging exascale Digital Twins of Earth, as interactive visualizations in real-time could remove access barriers and help democratize Earth system sciences.

Leaf-level coordination principles propagate to the ecosystem scale.

Fundamental axes of variation in plant traits result from trade-offs between costs and benefits of resource-use strategies at the leaf scale. However, it is unclear whether similar trade-offs propagate to the ecosystem level. Here, we test whether trait correlation patterns predicted by three well-known leaf- and plant-level coordination theories - the leaf economics spectrum, the global spectrum of plant form and function, and the least-cost hypothesis - are also observed between community mean traits and ecosystem processes. We combined ecosystem functional properties from FLUXNET sites, vegetation properties, and community mean plant traits into three corresponding principal component analyses. We find that the leaf economics spectrum (90 sites), the global spectrum of plant form and function (89 sites), and the least-cost hypothesis (82 sites) all propagate at the ecosystem level. However, we also find evidence of additional scale-emergent properties. Evaluating the coordination of ecosystem functional properties may aid the development of more realistic global dynamic vegetation models with critical empirical data, reducing the uncertainty of climate change projections.

Global variations in critical drought thresholds that impact vegetation.

Identifying the thresholds of drought that, if crossed, suppress vegetation functioning is vital for accurate quantification of how land ecosystems respond to climate variability and change. We present a globally applicable framework to identify drought thresholds for vegetation responses to different levels of known soil-moisture deficits using four remotely sensed vegetation proxies spanning 2001-2018. The thresholds identified represent critical inflection points for changing vegetation responses from highly resistant to highly vulnerable in response to drought stress, and as a warning signal for substantial vegetation impacts. Drought thresholds varied geographically, with much lower percentiles of soil-moisture anomalies in vegetated areas covered by more forests, corresponding to a comparably stronger capacity to mitigate soil water deficit stress in forested ecosystems. Generally, those lower thresholds are detected in more humid climates. State-of-the-art land models, however, overestimated thresholds of soil moisture (i.e. overestimating drought impacts), especially in more humid areas with higher forest covers and arid areas with few forest covers. Based on climate model projections, we predict that the risk of vegetation damage will increase by the end of the twenty-first century in some hotspots like East Asia, Europe, Amazon, southern Australia and eastern and southern Africa. Our data-based results will inform projections on future drought impacts on terrestrial ecosystems and provide an effective tool for drought management.

A standardized catalogue of spectral indices to advance the use of remote sensing in Earth system research.

Spectral Indices derived from multispectral remote sensing products are extensively used to monitor Earth system dynamics (e.g. vegetation dynamics, water bodies, fire regimes). The rapid increase of proposed spectral indices led to a high demand for catalogues of spectral indices and tools for their computation. However, most of these resources are either closed-source, outdated, unconnected to a catalogue or lacking a common Application Programming Interface (API). Here we present "Awesome Spectral Indices" (ASI), a standardized catalogue of spectral indices for Earth system research. ASI provides a comprehensive machine readable catalogue of spectral indices, which is linked to a Python library. ASI delivers a broad set of attributes for each spectral index, including names, formulas, and source references. The catalogue can be extended by the user community, ensuring that ASI remains current and enabling a wider range of scientific applications. Furthermore, the Python library enables the application of the catalogue to real-world data and thereby facilitates the efficient use of remote sensing resources in multiple Earth system domains.

The global spectrum of plant form and function: enhanced species-level trait dataset.

Here we provide the 'Global Spectrum of Plant Form and Function Dataset', containing species mean values for six vascular plant traits. Together, these traits -plant height, stem specific density, leaf area, leaf mass per area, leaf nitrogen content per dry mass, and diaspore (seed or spore) mass - define the primary axes of variation in plant form and function. The dataset is based on ca. 1 million trait records received via the TRY database (representing ca. 2,500 original publications) and additional unpublished data. It provides 92,159 species mean values for the six traits, covering 46,047 species. The data are complemented by higher-level taxonomic classification and six categorical traits (woodiness, growth form, succulence, adaptation to terrestrial or aquatic habitats, nutrition type and leaf type). Data quality management is based on a probabilistic approach combined with comprehensive validation against expert knowledge and external information. Intense data acquisition and thorough quality control produced the largest and, to our knowledge, most accurate compilation of empirically observed vascular plant species mean traits to date.

Citizen science plant observations encode global trait patterns.

Global maps of plant functional traits are essential for studying the dynamics of the terrestrial biosphere, yet the spatial distribution of trait measurements remains sparse. With the increasing popularity of species identification apps, citizen scientists contribute to growing vegetation data collections. The question emerges whether such opportunistic citizen science data can help map plant functional traits globally. Here we show that we can map global trait patterns by complementing vascular plant observations from the global citizen science project iNaturalist with measurements from the plant trait database TRY. We evaluate these maps using sPlotOpen, a global collection of vegetation plot data. Our results show high correlations between the iNaturalist- and sPlotOpen-based maps of up to 0.69 (r) and higher correlations than to previously published trait maps. As citizen science data collections continue to grow, we can expect them to play a significant role in further improving maps of plant functional traits.

Climatic and soil factors explain the two-dimensional spectrum of global plant trait variation.

Plant functional traits can predict community assembly and ecosystem functioning and are thus widely used in global models of vegetation dynamics and land-climate feedbacks. Still, we lack a global understanding of how land and climate affect plant traits. A previous global analysis of six traits observed two main axes of variation: (1) size variation at the organ and plant level and (2) leaf economics balancing leaf persistence against plant growth potential. The orthogonality of these two axes suggests they are differently influenced by environmental drivers. We find that these axes persist in a global dataset of 17 traits across more than 20,000 species. We find a dominant joint effect of climate and soil on trait variation. Additional independent climate effects are also observed across most traits, whereas independent soil effects are almost exclusively observed for economics traits. Variation in size traits correlates well with a latitudinal gradient related to water or energy limitation. In contrast, variation in economics traits is better explained by interactions of climate with soil fertility. These findings have the potential to improve our understanding of biodiversity patterns and our predictions of climate change impacts on biogeochemical cycles.

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.

Correction: Kernel methods and their derivatives: Concept and perspectives for the earth system sciences.

[This corrects the article DOI: 10.1371/journal.pone.0235885.].

A unified vegetation index for quantifying the terrestrial biosphere.

Empirical vegetation indices derived from spectral reflectance data are widely used in remote sensing of the biosphere, as they represent robust proxies for canopy structure, leaf pigment content, and, subsequently, plant photosynthetic potential. Here, we generalize the broad family of commonly used vegetation indices by exploiting all higher-order relations between the spectral channels involved. This results in a higher sensitivity to vegetation biophysical and physiological parameters. The presented nonlinear generalization of the celebrated normalized difference vegetation index (NDVI) consistently improves accuracy in monitoring key parameters, such as leaf area index, gross primary productivity, and sun-induced chlorophyll fluorescence. Results suggest that the statistical approach maximally exploits the spectral information and addresses long-standing problems in satellite Earth Observation of the terrestrial biosphere. The nonlinear NDVI will allow more accurate measures of terrestrial carbon source/sink dynamics and potentials for stabilizing atmospheric CO2 and mitigating global climate change.

Kernel methods and their derivatives: Concept and perspectives for the earth system sciences.

Kernel methods are powerful machine learning techniques which use generic non-linear functions to solve complex tasks. They have a solid mathematical foundation and exhibit excellent performance in practice. However, kernel machines are still considered black-box models as the kernel feature mapping cannot be accessed directly thus making the kernels difficult to interpret. The aim of this work is to show that it is indeed possible to interpret the functions learned by various kernel methods as they can be intuitive despite their complexity. Specifically, we show that derivatives of these functions have a simple mathematical formulation, are easy to compute, and can be applied to various problems. The model function derivatives in kernel machines is proportional to the kernel function derivative and we provide the explicit analytic form of the first and second derivatives of the most common kernel functions with regard to the inputs as well as generic formulas to compute higher order derivatives. We use them to analyze the most used supervised and unsupervised kernel learning methods: Gaussian Processes for regression, Support Vector Machines for classification, Kernel Entropy Component Analysis for density estimation, and the Hilbert-Schmidt Independence Criterion for estimating the dependency between random variables. For all cases we expressed the derivative of the learned function as a linear combination of the kernel function derivative. Moreover we provide intuitive explanations through illustrative toy examples and show how these same kernel methods can be applied to applications in the context of spatio-temporal Earth system data cubes. This work reflects on the observation that function derivatives may play a crucial role in kernel methods analysis and understanding.

Inferring causation from time series in Earth system sciences.

The heart of the scientific enterprise is a rational effort to understand the causes behind the phenomena we observe. In large-scale complex dynamical systems such as the Earth system, real experiments are rarely feasible. However, a rapidly increasing amount of observational and simulated data opens up the use of novel data-driven causal methods beyond the commonly adopted correlation techniques. Here, we give an overview of causal inference frameworks and identify promising generic application cases common in Earth system sciences and beyond. We discuss challenges and initiate the benchmark platform causeme.net to close the gap between method users and developers.

Publisher Correction: Stand age and species richness dampen interannual variation of ecosystem-level photosynthetic capacity.

In the version of this Article originally published, the wrong Supplementary Information pdf was uploaded, in which the figures did not correspond with those mentioned in the main text and the R code was not presented properly. This has now been replaced.

Global trait-environment relationships of plant communities.

Plant functional traits directly affect ecosystem functions. At the species level, trait combinations depend on trade-offs representing different ecological strategies, but at the community level trait combinations are expected to be decoupled from these trade-offs because different strategies can facilitate co-existence within communities. A key question is to what extent community-level trait composition is globally filtered and how well it is related to global versus local environmental drivers. Here, we perform a global, plot-level analysis of trait-environment relationships, using a database with more than 1.1 million vegetation plots and 26,632 plant species with trait information. Although we found a strong filtering of 17 functional traits, similar climate and soil conditions support communities differing greatly in mean trait values. The two main community trait axes that capture half of the global trait variation (plant stature and resource acquisitiveness) reflect the trade-offs at the species level but are weakly associated with climate and soil conditions at the global scale. Similarly, within-plot trait variation does not vary systematically with macro-environment. Our results indicate that, at fine spatial grain, macro-environmental drivers are much less important for functional trait composition than has been assumed from floristic analyses restricted to co-occurrence in large grid cells. Instead, trait combinations seem to be predominantly filtered by local-scale factors such as disturbance, fine-scale soil conditions, niche partitioning and biotic interactions.

Ecosystem functioning is enveloped by hydrometeorological variability.

Terrestrial ecosystem processes, and the associated vegetation carbon dynamics, respond differently to hydrometeorological variability across timescales, and so does our scientific understanding of the underlying mechanisms. Long-term variability of the terrestrial carbon cycle is not yet well constrained and the resulting climate-biosphere feedbacks are highly uncertain. Here we present a comprehensive overview of hydrometeorological and ecosystem variability from hourly to decadal timescales integrating multiple in situ and remote-sensing datasets characterizing extra-tropical forest sites. We find that ecosystem variability at all sites is confined within a hydrometeorological envelope across sites and timescales. Furthermore, ecosystem variability demonstrates long-term persistence, highlighting ecological memory and slow ecosystem recovery rates after disturbances. However, simulation results with state-of-the-art process-based models do not reflect this long-term persistent behaviour in ecosystem functioning. Accordingly, we develop a cross-time-scale stochastic framework that captures hydrometeorological and ecosystem variability. Our analysis offers a perspective for terrestrial ecosystem modelling and paves the way for new model-data integration opportunities in Earth system sciences.

Stand age and species richness dampen interannual variation of ecosystem-level photosynthetic capacity.

The total uptake of carbon dioxide by ecosystems via photosynthesis (gross primary productivity, GPP) is the largest flux in the global carbon cycle. A key ecosystem functional property determining GPP is the photosynthetic capacity at light saturation (GPPsat), and its interannual variability (IAV) is propagated to the net land-atmosphere exchange of CO2. Given the importance of understanding the IAV in CO2 fluxes for improving the predictability of the global carbon cycle, we have tested a range of alternative hypotheses to identify potential drivers of the magnitude of IAV in GPPsat in forest ecosystems. Our results show that while the IAV in GPPsat within sites is closely related to air temperature and soil water availability fluctuations, the magnitude of IAV in GPPsat is related to stand age and biodiversity (R2 = 0.55, P < 0.0001). We find that the IAV of GPPsat is greatly reduced in older and more diverse forests, and is higher in younger forests with few dominant species. Older and more diverse forests seem to dampen the effect of climate variability on the carbon cycle irrespective of forest type. Preserving old forests and their diversity would therefore be beneficial in reducing the effect of climate variability on Earth's forest ecosystems.

Stand age and species richness dampen interannual variation of ecosystem-level photosynthetic capacity.

The total uptake of carbon dioxide by ecosystems via photosynthesis (gross primary productivity, GPP) is the largest flux in the global carbon cycle. A key ecosystem functional property determining GPP is the photosynthetic capacity at light saturation (GPPsat), and its interannual variability (IAV) is propagated to the net land-atmosphere exchange of CO2. Given the importance of understanding the IAV in CO2 fluxes for improving the predictability of the global carbon cycle, we have tested a range of alternative hypotheses to identify potential drivers of the magnitude of IAV in GPPsat in forest ecosystems. Our results show that while the IAV in GPPsat within sites is closely related to air temperature and soil water availability fluctuations, the magnitude of IAV in GPPsat is related to stand age and biodiversity (R2 = 0.55, P < 0.0001). We find that the IAV of GPPsat is greatly reduced in older and more diverse forests, and is higher in younger forests with few dominant species. Older and more diverse forests seem to dampen the effect of climate variability on the carbon cycle irrespective of forest type. Preserving old forests and their diversity would therefore be beneficial in reducing the effect of climate variability on Earth's forest ecosystems.