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.

The global spectrum of plant form and function.

Earth is home to a remarkable diversity of plant forms and life histories, yet comparatively few essential trait combinations have proved evolutionarily viable in today's terrestrial biosphere. By analysing worldwide variation in six major traits critical to growth, survival and reproduction within the largest sample of vascular plant species ever compiled, we found that occupancy of six-dimensional trait space is strongly concentrated, indicating coordination and trade-offs. Three-quarters of trait variation is captured in a two-dimensional global spectrum of plant form and function. One major dimension within this plane reflects the size of whole plants and their parts; the other represents the leaf economics spectrum, which balances leaf construction costs against growth potential. The global plant trait spectrum provides a backdrop for elucidating constraints on evolution, for functionally qualifying species and ecosystems, and for improving models that predict future vegetation based on continuous variation in plant form and function.

Climate extremes and the carbon cycle.

The terrestrial biosphere is a key component of the global carbon cycle and its carbon balance is strongly influenced by climate. Continuing environmental changes are thought to increase global terrestrial carbon uptake. But evidence is mounting that climate extremes such as droughts or storms can lead to a decrease in regional ecosystem carbon stocks and therefore have the potential to negate an expected increase in terrestrial carbon uptake. Here we explore the mechanisms and impacts of climate extremes on the terrestrial carbon cycle, and propose a pathway to improve our understanding of present and future impacts of climate extremes on the terrestrial carbon budget.

Plants, microorganisms, and soil temperatures contribute to a decrease in methane fluxes on a drained Arctic floodplain.

As surface temperatures are expected to rise in the future, ice-rich permafrost may thaw, altering soil topography and hydrology and creating a mosaic of wet and dry soil surfaces in the Arctic. Arctic wetlands are large sources of CH4 , and investigating effects of soil hydrology on CH4 fluxes is of great importance for predicting ecosystem feedback in response to climate change. In this study, we investigate how a decade-long drying manipulation on an Arctic floodplain influences CH4 -associated microorganisms, soil thermal regimes, and plant communities. Moreover, we examine how these drainage-induced changes may then modify CH4 fluxes in the growing and nongrowing seasons. This study shows that drainage substantially lowered the abundance of methanogens along with methanotrophic bacteria, which may have reduced CH4 cycling. Soil temperatures of the drained areas were lower in deep, anoxic soil layers (below 30 cm), but higher in oxic topsoil layers (0-15 cm) compared to the control wet areas. This pattern of soil temperatures may have reduced the rates of methanogenesis while elevating those of CH4 oxidation, thereby decreasing net CH4 fluxes. The abundance of Eriophorum angustifolium, an aerenchymatous plant species, diminished significantly in the drained areas. Due to this decrease, a higher fraction of CH4 was alternatively emitted to the atmosphere by diffusion, possibly increasing the potential for CH4 oxidation and leading to a decrease in net CH4 fluxes compared to a control site. Drainage lowered CH4 fluxes by a factor of 20 during the growing season, with postdrainage changes in microbial communities, soil temperatures, and plant communities also contributing to this reduction. In contrast, we observed CH4 emissions increased by 10% in the drained areas during the nongrowing season, although this difference was insignificant given the small magnitudes of fluxes. This study showed that long-term drainage considerably reduced CH4 fluxes through modified ecosystem properties.

Linking plant and ecosystem functional biogeography.

Classical biogeographical observations suggest that ecosystems are strongly shaped by climatic constraints in terms of their structure and function. On the other hand, vegetation function feeds back on the climate system via biosphere-atmosphere exchange of matter and energy. Ecosystem-level observations of this exchange reveal very large functional biogeographical variation of climate-relevant ecosystem functional properties related to carbon and water cycles. This variation is explained insufficiently by climate control and a classical plant functional type classification approach. For example, correlations between seasonal carbon-use efficiency and climate or environmental variables remain below 0.6, leaving almost 70% of variance unexplained. We suggest that a substantial part of this unexplained variation of ecosystem functional properties is related to variations in plant and microbial traits. Therefore, to progress with global functional biogeography, we should seek to understand the link between organismic traits and flux-derived ecosystem properties at ecosystem observation sites and the spatial variation of vegetation traits given geoecological covariates. This understanding can be fostered by synergistic use of both data-driven and theory-driven ecological as well as biophysical approaches.

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.

Effects of climate extremes on the terrestrial carbon cycle: concepts, processes and potential future impacts.

Extreme droughts, heat waves, frosts, precipitation, wind storms and other climate extremes may impact the structure, composition and functioning of terrestrial ecosystems, and thus carbon cycling and its feedbacks to the climate system. Yet, the interconnected avenues through which climate extremes drive ecological and physiological processes and alter the carbon balance are poorly understood. Here, we review the literature on carbon cycle relevant responses of ecosystems to extreme climatic events. Given that impacts of climate extremes are considered disturbances, we assume the respective general disturbance-induced mechanisms and processes to also operate in an extreme context. The paucity of well-defined studies currently renders a quantitative meta-analysis impossible, but permits us to develop a deductive framework for identifying the main mechanisms (and coupling thereof) through which climate extremes may act on the carbon cycle. We find that ecosystem responses can exceed the duration of the climate impacts via lagged effects on the carbon cycle. The expected regional impacts of future climate extremes will depend on changes in the probability and severity of their occurrence, on the compound effects and timing of different climate extremes, and on the vulnerability of each land-cover type modulated by management. Although processes and sensitivities differ among biomes, based on expert opinion, we expect forests to exhibit the largest net effect of extremes due to their large carbon pools and fluxes, potentially large indirect and lagged impacts, and long recovery time to regain previous stocks. At the global scale, we presume that droughts have the strongest and most widespread effects on terrestrial carbon cycling. Comparing impacts of climate extremes identified via remote sensing vs. ground-based observational case studies reveals that many regions in the (sub-)tropics are understudied. Hence, regional investigations are needed to allow a global upscaling of the impacts of climate extremes on global carbon-climate feedbacks.

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.

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.

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.

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.

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.

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.

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.

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.

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.