'rtry': An R package to support plant trait data preprocessing.

Plant trait data are used to quantify how plants respond to environmental factors and can act as indicators of ecosystem function. Measured trait values are influenced by genetics, trade-offs, competition, environmental conditions, and phenology. These interacting effects on traits are poorly characterized across taxa, and for many traits, measurement protocols are not standardized. As a result, ancillary information about growth and measurement conditions can be highly variable, requiring a flexible data structure. In 2007, the TRY initiative was founded as an integrated database of plant trait data, including ancillary attributes relevant to understanding and interpreting the trait values. The TRY database now integrates around 700 original and collective datasets and has become a central resource of plant trait data. These data are provided in a generic long-table format, where a unique identifier links different trait records and ancillary data measured on the same entity. Due to the high number of trait records, plant taxa, and types of traits and ancillary data released from the TRY database, data preprocessing is necessary but not straightforward. Here, we present the 'rtry' R package, specifically designed to support plant trait data exploration and filtering. By integrating a subset of existing R functions essential for preprocessing, 'rtry' avoids the need for users to navigate the extensive R ecosystem and provides the functions under a consistent syntax. 'rtry' is therefore easy to use even for beginners in R. Notably, 'rtry' does not support data retrieval or analysis; rather, it focuses on the preprocessing tasks to optimize data quality. While 'rtry' primarily targets TRY data, its utility extends to data from other sources, such as the National Ecological Observatory Network (NEON). The 'rtry' package is available on the Comprehensive R Archive Network (CRAN; https://cran.r-project.org/package=rtry) and the GitHub Wiki (https://github.com/MPI-BGC-Functional-Biogeography/rtry/wiki) along with comprehensive documentation and vignettes describing detailed data preprocessing workflows.

Global patterns of tree wood density.

Wood density is a fundamental property related to tree biomechanics and hydraulic function while playing a crucial role in assessing vegetation carbon stocks by linking volumetric retrieval and a mass estimate. This study provides a high-resolution map of the global distribution of tree wood density at the 0.01° (~1 km) spatial resolution, derived from four decision trees machine learning models using a global database of 28,822 tree-level wood density measurements. An ensemble of four top-performing models combined with eight cross-validation strategies shows great consistency, providing wood density patterns with pronounced spatial heterogeneity. The global pattern shows lower wood density values in northern and northwestern Europe, Canadian forest regions and slightly higher values in Siberia forests, western United States, and southern China. In contrast, tropical regions, especially wet tropical areas, exhibit high wood density. Climatic predictors explain 49%-63% of spatial variations, followed by vegetation characteristics (25%-31%) and edaphic properties (11%-16%). Notably, leaf type (evergreen vs. deciduous) and leaf habit type (broadleaved vs. needleleaved) are the most dominant individual features among all selected predictive covariates. Wood density tends to be higher for angiosperm broadleaf trees compared to gymnosperm needleleaf trees, particularly for evergreen species. The distributions of wood density categorized by leaf types and leaf habit types have good agreement with the features observed in wood density measurements. This global map quantifying wood density distribution can help improve accurate predictions of forest carbon stocks, providing deeper insights into ecosystem functioning and carbon cycling such as forest vulnerability to hydraulic and thermal stresses in the context of future climate change.

Rooting depth and xylem vulnerability are independent woody plant traits jointly selected by aridity, seasonality, and water table depth.

Evolutionary radiations of woody taxa within arid environments were made possible by multiple trait innovations including deep roots and embolism-resistant xylem, but little is known about how these traits have coevolved across the phylogeny of woody plants or how they jointly influence the distribution of species. We synthesized global trait and vegetation plot datasets to examine how rooting depth and xylem vulnerability across 188 woody plant species interact with aridity, precipitation seasonality, and water table depth to influence species occurrence probabilities across all biomes. Xylem resistance to embolism and rooting depth are independent woody plant traits that do not exhibit an interspecific trade-off. Resistant xylem and deep roots increase occurrence probabilities in arid, seasonal climates over deep water tables. Resistant xylem and shallow roots increase occurrence probabilities in arid, nonseasonal climates over deep water tables. Vulnerable xylem and deep roots increase occurrence probabilities in arid, nonseasonal climates over shallow water tables. Lastly, vulnerable xylem and shallow roots increase occurrence probabilities in humid climates. Each combination of trait values optimizes occurrence probabilities in unique environmental conditions. Responses of deeply rooted vegetation may be buffered if evaporative demand changes faster than water table depth under climate change.

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.

The big four of plant taxonomy - a comparison of global checklists of vascular plant names.

Taxonomic checklists used to verify published plant names and identify synonyms are a cornerstone of biological research. Four global authoritative checklists for vascular plants exist: Leipzig Catalogue of Vascular Plants, World Checklist of Vascular Plants, World Flora Online (successor of The Plant List, TPL), and WorldPlants. We compared these four checklists in terms of size and differences across taxa. We matched taxon names of these checklists and TPL against each other, identified differences across checklists, and evaluated the consistency of accepted names linked to individual taxon names. We assessed geographic and phylogenetic patterns of variance. All checklists differed strongly compared with TPL and provided identical information on c. 60% of plant names. Geographically, differences in checklists increased from low to high latitudes. Phylogenetically, we detected strong variability across families. A comparison of name-matching performance on taxon names submitted to the functional trait database TRY, and a check of completeness of accepted names evaluated against an independent, expert-curated checklist of the family Meliaceae, showed a similar performance across checklists. This study raises awareness on the differences in data and approach across these checklists potentially impacting analyses. We propose ideas on the way forward exploring synergies and harmonizing the four global checklists.

Forms of nitrogen inputs regulate the intensity of soil acidification.

Soil acidification induced by reactive nitrogen (N) inputs can alter the structure and function of terrestrial ecosystems. Because different N-transformation processes contribute to the production and consumption of H+ , the magnitude of acidification likely depends on the relative amounts of organic N (ON) and inorganic N (IN) inputs. However, few studies have explicitly measured the effects of N composition on soil acidification. In this study, we first conducted a meta-analysis to test the effects of ON or IN inputs on soil acidification across 53 studies in grasslands. We then compared soil acidification across five different ON:IN ratios and two input rates based on long-term field N addition experiments. The meta-analysis showed that ON had weaker effects on soil acidification than IN when the N addition rate was above 20 g N m-2 year-1 . The field experiment confirmed the findings from meta-analysis: N addition with proportions of ON ≥ 20% caused less soil acidification, especially at a high input rate (30 g N m-2 year-1 ). Structural equation model analysis showed that this result was largely due to a relatively low rate of H+ production from ON as NH3 volatilization and uptake of ON and NH4 + by the dominant grass species Leymus chinensis (which are both lower net contributors to H+ production) result in less NH4 + available for nitrification (which is a higher net contributor to H+ production). These results indicate that the evaluation of soil acidification induced by N inputs should consider N forms and manipulations of relative composition of N inputs may provide an effective approach to alleviate the N-induced soil acidification.

Global beta-diversity of angiosperm trees is shaped by Quaternary climate change.

As Earth's climate has varied strongly through geological time, studying the impacts of past climate change on biodiversity helps to understand the risks from future climate change. However, it remains unclear how paleoclimate shapes spatial variation in biodiversity. Here, we assessed the influence of Quaternary climate change on spatial dissimilarity in taxonomic, phylogenetic, and functional composition among neighboring 200-kilometer cells (beta-diversity) for angiosperm trees worldwide. We found that larger glacial-interglacial temperature change was strongly associated with lower spatial turnover (species replacements) and higher nestedness (richness changes) components of beta-diversity across all three biodiversity facets. Moreover, phylogenetic and functional turnover was lower and nestedness higher than random expectations based on taxonomic beta-diversity in regions that experienced large temperature change, reflecting phylogenetically and functionally selective processes in species replacement, extinction, and colonization during glacial-interglacial oscillations. Our results suggest that future human-driven climate change could cause local homogenization and reduction in taxonomic, phylogenetic, and functional diversity of angiosperm trees worldwide.

Climate-trait relationships exhibit strong habitat specificity in plant communities across Europe.

Ecological theory predicts close relationships between macroclimate and functional traits. Yet, global climatic gradients correlate only weakly with the trait composition of local plant communities, suggesting that important factors have been ignored. Here, we investigate the consistency of climate-trait relationships for plant communities in European habitats. Assuming that local factors are better accounted for in more narrowly defined habitats, we assigned > 300,000 vegetation plots to hierarchically classified habitats and modelled the effects of climate on the community-weighted means of four key functional traits using generalized additive models. We found that the predictive power of climate increased from broadly to narrowly defined habitats for specific leaf area and root length, but not for plant height and seed mass. Although macroclimate generally predicted the distribution of all traits, its effects varied, with habitat-specificity increasing toward more narrowly defined habitats. We conclude that macroclimate is an important determinant of terrestrial plant communities, but future predictions of climatic effects must consider how habitats are defined.

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.

Functional rarity of plants in German hay meadows - Patterns on the species level and mismatches with community species richness.

Functional rarity (FR) - a feature combining a species' rarity with the distinctiveness of its traits - is a promising tool to better understand the ecological importance of rare species and consequently to protect functional diversity more efficiently. However, we lack a systematic understanding of FR on both the species level (which species are functionally rare and why) and the community level (how is FR associated with biodiversity and environmental conditions). Here, we quantify FR for 218 plant species from German hay meadows on a local, regional, and national scale by combining data from 6500 vegetation relevés and 15 ecologically relevant traits. We investigate the association between rarity and trait distinctiveness on different spatial scales via correlation measures and show which traits lead to low or high trait distinctiveness via distance-based redundancy analysis. We test how species richness and FR are correlated, and use boosted regression trees to determine environmental conditions that are driving species richness and FR. On the local scale, only rare species showed high trait distinctiveness while on larger spatial scales rare and common species showed high trait distinctiveness. As infrequent trait attributes (e.g., legumes, low clonality) led to higher trait distinctiveness, we argue that functionally rare species are either specialists or transients. While specialists occupy a particular niche in hay meadows leading to lower rarity on larger spatial scales, transients display distinct but maladaptive traits resulting in high rarity across all spatial scales. More functionally rare species than expected by chance occurred in species-poor communities indicating that they prefer environmental conditions differing from characteristic conditions of species-rich hay meadows. Finally, we argue that functionally rare species are not necessarily relevant for nature conservation because many were transients from surrounding habitats. However, FR can facilitate our understanding of why species are rare in a habitat and under which conditions these species occur.

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.

Global relationships in tree functional traits.

Due to massive energetic investments in woody support structures, trees are subject to unique physiological, mechanical, and ecological pressures not experienced by herbaceous plants. Despite a wealth of studies exploring trait relationships across the entire plant kingdom, the dominant traits underpinning these unique aspects of tree form and function remain unclear. Here, by considering 18 functional traits, encompassing leaf, seed, bark, wood, crown, and root characteristics, we quantify the multidimensional relationships in tree trait expression. We find that nearly half of trait variation is captured by two axes: one reflecting leaf economics, the other reflecting tree size and competition for light. Yet these orthogonal axes reveal strong environmental convergence, exhibiting correlated responses to temperature, moisture, and elevation. By subsequently exploring multidimensional trait relationships, we show that the full dimensionality of trait space is captured by eight distinct clusters, each reflecting a unique aspect of tree form and function. Collectively, this work identifies a core set of traits needed to quantify global patterns in functional biodiversity, and it contributes to our fundamental understanding of the functioning of forests worldwide.

High exposure of global tree diversity to human pressure.

Safeguarding Earth's tree diversity is a conservation priority due to the importance of trees for biodiversity and ecosystem functions and services such as carbon sequestration. Here, we improve the foundation for effective conservation of global tree diversity by analyzing a recently developed database of tree species covering 46,752 species. We quantify range protection and anthropogenic pressures for each species and develop conservation priorities across taxonomic, phylogenetic, and functional diversity dimensions. We also assess the effectiveness of several influential proposed conservation prioritization frameworks to protect the top 17% and top 50% of tree priority areas. We find that an average of 50.2% of a tree species' range occurs in 110-km grid cells without any protected areas (PAs), with 6,377 small-range tree species fully unprotected, and that 83% of tree species experience nonnegligible human pressure across their range on average. Protecting high-priority areas for the top 17% and 50% priority thresholds would increase the average protected proportion of each tree species' range to 65.5% and 82.6%, respectively, leaving many fewer species (2,151 and 2,010) completely unprotected. The priority areas identified for trees match well to the Global 200 Ecoregions framework, revealing that priority areas for trees would in large part also optimize protection for terrestrial biodiversity overall. Based on range estimates for >46,000 tree species, our findings show that a large proportion of tree species receive limited protection by current PAs and are under substantial human pressure. Improved protection of biodiversity overall would also strongly benefit global tree diversity.

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.

Long-term leaf C:N ratio change under elevated CO2 and nitrogen deposition in China: Evidence from observations and process-based modeling.

Climate change, elevating atmosphere CO2 (eCO2) and increased nitrogen deposition (iNDEP) are altering the biogeochemical interactions between plants, microbes and soils, which further modify plant leaf carbon­nitrogen (C:N) stoichiometry and their carbon assimilation capability. Many field experiments have observed large sensitivity of leaf C:N ratio to eCO2 and iNDEP. However, the large-scale pattern of this sensitivity is still unclear, because eCO2 and iNDEP drive leaf C:N ratio toward opposite directions, which are further compounded by the complex processes of nitrogen acquisition and plant-and-microbial nitrogen competition. Here, we attempt to map the leaf C:N ratio spatial variation in the past 5 decades in China with a combination of data-driven model and process-based modeling. These two approaches showed consistent results. Over different regions, we found that leaf C:N ratio had significant but uneven changes between 2 time periods (1960-1989 and 1990-2015): a 5% ± 8% increase for temperate grasslands in northern China, a 3% ± 6% increase for boreal grasslands in western China, and by contrast, a 7% ± 6% decrease for temperate forests in southern China, and a 3% ± 5% decrease for boreal forests in northeastern China. Additionally, the structural equation models indicated that the leaf C:N change was sensitive to ΔNDEP, ΔCO2 and ΔMAT rather than ΔMAP and ecosystem types. Process-based modeling suggested that iNDEP was the main source of soil mineral nitrogen change, dominating leaf C:N ratio change in most areas in China, while eCO2 led to leaf C:N ratio increase in low iNDEP area. This study also indicates that the long-term leaf C:N ratio acclimation was dominated by climate constraint, especially temperature, but was constrained by soil N availability over decade scale.

Comprehensive leaf size traits dataset for seven plant species from digitised herbarium specimen images covering more than two centuries.

BACKGROUND: Morphological leaf traits are frequently used to quantify, understand and predict plant and vegetation functional diversity and ecology, including environmental and climate change responses. Although morphological leaf traits are easy to measure, their coverage for characterising variation within species and across temporal scales is limited. At the same time, there are about 3100 herbaria worldwide, containing approximately 390 million plant specimens dating from the 16th to 21st century, which can potentially be used to extract morphological leaf traits. Globally, plant specimens are rapidly being digitised and images are made openly available via various biodiversity data platforms, such as iDigBio and GBIF. Based on a pilot study to identify the availability and appropriateness of herbarium specimen images for comprehensive trait data extraction, we developed a spatio-temporal dataset on intraspecific trait variability containing 128,036 morphological leaf trait measurements for seven selected species. NEW INFORMATION: After scrutinising the metadata of digitised herbarium specimen images available from iDigBio and GBIF (21.9 million and 31.6 million images for Tracheophyta; accessed date December 2020), we identified approximately 10 million images potentially appropriate for our study. From the 10 million images, we selected seven species (Salix bebbiana Sarg., Alnus incana (L.) Moench, Viola canina L., Salix glauca L., Chenopodium album L., Impatiens capensis Meerb. and Solanum dulcamara L.) , which have a simple leaf shape, are well represented in space and time and have high availability of specimens per species. We downloaded 17,383 images. Out of these, we discarded 5779 images due to quality issues. We used the remaining 11,604 images to measure the area, length, width and perimeter on 32,009 individual leaf blades using the semi-automated tool TraitEx. The resulting dataset contains 128,036 trait records.We demonstrate its comparability to trait data measured in natural environments following standard protocols by comparing trait values from the TRY database. We conclude that the herbarium specimens provide valuable information on leaf sizes. The dataset created in our study, by extracting leaf traits from the digitised herbarium specimen images of seven selected species, is a promising opportunity to improve ecological knowledge about the adaptation of size-related leaf traits to environmental changes in space and time.

An integrated framework of plant form and function: the belowground perspective.

Plant trait variation drives plant function, community composition and ecosystem processes. However, our current understanding of trait variation disproportionately relies on aboveground observations. Here we integrate root traits into the global framework of plant form and function. We developed and tested an overarching conceptual framework that integrates two recently identified root trait gradients with a well-established aboveground plant trait framework. We confronted our novel framework with published relationships between above- and belowground trait analogues and with multivariate analyses of above- and belowground traits of 2510 species. Our traits represent the leaf and root conservation gradients (specific leaf area, leaf and root nitrogen concentration, and root tissue density), the root collaboration gradient (root diameter and specific root length) and the plant size gradient (plant height and rooting depth). We found that an integrated, whole-plant trait space required as much as four axes. The two main axes represented the fast-slow 'conservation' gradient on which leaf and fine-root traits were well aligned, and the 'collaboration' gradient in roots. The two additional axes were separate, orthogonal plant size axes for height and rooting depth. This perspective on the multidimensional nature of plant trait variation better encompasses plant function and influence on the surrounding environment.