Reading tea leaves worldwide: Decoupled drivers of initial litter decomposition mass-loss rate and stabilization

The breakdown of plant material fuels soil functioning and biodiversity. Currently, process understanding of global decomposition patterns and the drivers of such patterns are hampered by the lack of coherent large-scale datasets. We buried 36,000 individual litterbags (tea bags) worldwide and found an overall negative correlation between initial mass-loss rates and stabilization factors of plant-derived carbon, using the Tea Bag Index (TBI). The stabilization factor quantifies the degree to which easy-to-degrade components accumulate during early-stage decomposition (e.g. by environmental limitations). However, agriculture and an interaction between moisture and temperature led to a decoupling between initial mass-loss rates and stabilization, notably in colder locations. Using TBI improved mass-loss estimates of natural litter compared to models that ignored stabilization. Ignoring the transformation of dead plant material to more recalcitrant substances during early-stage decomposition, and the environmental control of this transformation, could overestimate carbon losses during early decomposition in carbon cycle models.

Integrated global assessment of the natural forest carbon potential

Forests are a substantial terrestrial carbon sink, but anthropogenic changes in land use and climate have considerably reduced the scale of this system1. Remote-sensing estimates to quantify carbon losses from global forests2-5 are characterized by considerable uncertainty and we lack a comprehensive ground-sourced evaluation to benchmark these estimates. Here we combine several ground-sourced6 and satellite-derived approaches2,7,8 to evaluate the scale of the global forest carbon potential outside agricultural and urban lands. Despite regional variation, the predictions demonstrated remarkable consistency at a global scale, with only a 12% difference between the ground-sourced and satellite-derived estimates. At present, global forest carbon storage is markedly under the natural potential, with a total deficit of 226 Gt (model range = 151-363 Gt) in areas with low human footprint. Most (61%, 139 Gt C) of this potential is in areas with existing forests, in which ecosystem protection can allow forests to recover to maturity. The remaining 39% (87 Gt C) of potential lies in regions in which forests have been removed or fragmented. Although forests cannot be a substitute for emissions reductions, our results support the idea2,3,9 that the conservation, restoration and sustainable management of diverse forests offer valuable contributions to meeting global climate and biodiversity targets.

The global biogeography of tree leaf form and habit

Understanding what controls global leaf type variation in trees is crucial for comprehending their role in terrestrial ecosystems, including carbon, water and nutrient dynamics. Yet our understanding of the factors influencing forest leaf types remains incomplete, leaving us uncertain about the global proportions of needle-leaved, broadleaved, evergreen and deciduous trees. To address these gaps, we conducted a global, ground-sourced assessment of forest leaf-type variation by integrating forest inventory data with comprehensive leaf form (broadleaf vs needle-leaf) and habit (evergreen vs deciduous) records. We found that global variation in leaf habit is primarily driven by isothermality and soil characteristics, while leaf form is predominantly driven by temperature. Given these relationships, we estimate that 38% of global tree individuals are needle-leaved evergreen, 29% are broadleaved evergreen, 27% are broadleaved deciduous and 5% are needle-leaved deciduous. The aboveground biomass distribution among these tree types is approximately 21% (126.4 Gt), 54% (335.7 Gt), 22% (136.2 Gt) and 3% (18.7 Gt), respectively. We further project that, depending on future emissions pathways, 17-34% of forested areas will experience climate conditions by the end of the century that currently support a different forest type, highlighting the intensification of climatic stress on existing forests. By quantifying the distribution of tree leaf types and their corresponding biomass, and identifying regions where climate change will exert greatest pressure on current leaf types, our results can help improve predictions of future terrestrial ecosystem functioning and carbon cycling.

Native diversity buffers against severity of non-native tree invasions

Determining the drivers of non-native plant invasions is critical for managing native ecosystems and limiting the spread of invasive species1,2. Tree invasions in particular have been relatively overlooked, even though they have the potential to transform ecosystems and economies3,4. Here, leveraging global tree databases5-7, we explore how the phylogenetic and functional diversity of native tree communities, human pressure and the environment influence the establishment of non-native tree species and the subsequent invasion severity. We find that anthropogenic factors are key to predicting whether a location is invaded, but that invasion severity is underpinned by native diversity, with higher diversity predicting lower invasion severity. Temperature and precipitation emerge as strong predictors of invasion strategy, with non-native species invading successfully when they are similar to the native community in cold or dry extremes. Yet, despite the influence of these ecological forces in determining invasion strategy, we find evidence that these patterns can be obscured by human activity, with lower ecological signal in areas with higher proximity to shipping ports. Our global perspective of non-native tree invasion highlights that human drivers influence non-native tree presence, and that native phylogenetic and functional diversity have a critical role in the establishment and spread of subsequent invasions.

ForestClim-Bioclimatic variables for microclimate temperatures of European forests

Microclimate research gained renewed interest over the last decade and its importance for many ecological processes is increasingly being recognized. Consequently, the call for high-resolution microclimatic temperature grids across broad spatial extents is becoming more pressing to improve ecological models. Here, we provide a new set of open-access bioclimatic variables for microclimate temperatures of European forests at 25 × 25 m² resolution.

Macroevolutionary decline in mycorrhizal colonization and chemical defense responsiveness to mycorrhization

Arbuscular mycorrhizal fungi (AMF) have evolved associations with roots of 60% plant species, but the net benefit for plants vary broadly from mutualism to parasitism. Yet, we lack a general understanding of the evolutionary and ecological forces driving such variation. To this end, we conducted a comparative phylogenetic experiment with 24 species of Plantago, encompassing worldwide distribution, to address the effect of evolutionary history and environment on plant growth and chemical defenses in response to AMF colonization. We demonstrate that different species within one plant genus vary greatly in their ability to associate with AMF, and that AMF arbuscule colonization intensity decreases monotonically with increasing phylogenetic branch length, but not with concomitant changes in pedological and climatic conditions across species. Moreover, we demonstrate that species with the highest colonization levels are also those that change their defensive chemistry the least. We propose that the costs imposed by high AMF colonization in terms of reduced changes in secondary chemistry might drive the observed macroevolutionary decline in mycorrhization.

Global maps of soil temperature

Research in global change ecology relies heavily on global climatic grids derived from estimates of air temperature in open areas at around 2 m above the ground. These climatic grids do not reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions occur and most terrestrial species reside. Here, we provide global maps of soil temperature and bioclimatic variables at a 1-km2 resolution for 0-5 and 5-15 cm soil depth. These maps were created by calculating the difference (i.e. offset) between in situ soil temperature measurements, based on time series from over 1200 1-km2 pixels (summarized from 8519 unique temperature sensors) across all the world's major terrestrial biomes, and coarse-grained air temperature estimates from ERA5-Land (an atmospheric reanalysis by the European Centre for Medium-Range Weather Forecasts). We show that mean annual soil temperature differs markedly from the corresponding gridded air temperature, by up to 10°C (mean = 3.0 ± 2.1°C), with substantial variation across biomes and seasons. Over the year, soils in cold and/or dry biomes are substantially warmer (+3.6 ± 2.3°C) than gridded air temperature, whereas soils in warm and humid environments are on average slightly cooler (-0.7 ± 2.3°C). The observed substantial and biome-specific offsets emphasize that the projected impacts of climate and climate change on near-surface biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments. The global soil-related bioclimatic variables provided here are an important step forward for any application in ecology and related disciplines. Nevertheless, we highlight the need to fill remaining geographic gaps by collecting more in situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications.

Priming effects in soils across Europe

Land use is a key factor driving changes in soil carbon (C) cycle and contents worldwide. The priming effect (PE)-CO₂ emissions from changed soil organic matter decomposition in response to fresh C inputs-is one of the most unpredictable phenomena associated with C cycling and related nutrient mobilization. Yet, we know very little about the influence of land use on soil PE across contrasting environments. Here, we conducted a continental-scale study to (i) determine the PE induced by ¹³ C-glucose additions to 126 cropland and seminatural (forests and grasslands) soils from 22 European countries; (ii) compare PE magnitude in soils under various crop types (i.e., cereals, nonpermanent industrial crops, and orchards); and (iii) model the environmental factors influencing PE. On average, PEs were negative in seminatural (with values ranging between -60 and 26 µg C g^-1 soil after 35 days of incubation; median = -11) and cropland (from -55 to 27 µC g^-1 soil; median = -4.3) soils, meaning that microbial communities preferentially switched from soil organic C decomposition to glucose mineralization. PE was significantly less negative in croplands compared with seminatural ecosystems and not influenced by the crop type. PE was driven by soil basal respiration (reflecting microbial activity), microbial biomass C, and soil organic C, which were all higher in seminatural ecosystems compared with croplands. This cross European experimental and modeling study elucidated that PE intensity is dependent on land use and allowed to clarify the factors regulating this important C cycling process.

ForestTemp - Sub-canopy microclimate temperatures of European forests

Ecological research heavily relies on coarse-gridded climate data based on standardized temperature measurements recorded at 2 m height in open landscapes. However, many organisms experience environmental conditions that differ substantially from those captured by these macroclimatic (i.e. free air) temperature grids. In forests, the tree canopy functions as a thermal insulator and buffers sub-canopy microclimatic conditions, thereby affecting biological and ecological processes. To improve the assessment of climatic conditions and climate-change-related impacts on forest-floor biodiversity and functioning, high-resolution temperature grids reflecting forest microclimates are thus urgently needed. Combining more than 1200 time series of in situ near-surface forest temperature with topographical, biological and macroclimatic variables in a machine learning model, we predicted the mean monthly offset between sub-canopy temperature at 15 cm above the surface and free-air temperature over the period 2000-2020 at a spatial resolution of 25 m across Europe. This offset was used to evaluate the difference between microclimate and macroclimate across space and seasons and finally enabled us to calculate mean annual and monthly temperatures for European forest understories. We found that sub-canopy air temperatures differ substantially from free-air temperatures, being on average 2.1°C (standard deviation ± 1.6°C) lower in summer and 2.0°C higher (±0.7°C) in winter across Europe. Additionally, our high-resolution maps expose considerable microclimatic variation within landscapes, not captured by the gridded macroclimatic products. The provided forest sub-canopy temperature maps will enable future research to model below-canopy biological processes and patterns, as well as species distributions more accurately.

The global distribution and environmental drivers of aboveground versus belowground plant biomass

A poor understanding of the fraction of global plant biomass occurring belowground as roots limits our understanding of present and future ecosystem function and carbon pools. Here we create a database of root-mass fractions (RMFs), an index of plant below- versus aboveground biomass distributions, and generate quantitative, spatially explicit global maps of RMFs in trees, shrubs and grasses. Our analyses reveal large gradients in RMFs both across and within vegetation types that can be attributed to resource availability. High RMFs occur in cold and dry ecosystems, while low RMFs dominate in warm and wet regions. Across all vegetation types, the directional effect of temperature on RMFs depends on water availability, suggesting feedbacks between heat, water and nutrient supply. By integrating our RMF maps with existing aboveground plant biomass information, we estimate that in forests, shrublands and grasslands, respectively, 22%, 47% and 67% of plant biomass exists belowground, with a total global belowground fraction of 24% (20-28%), that is, 113 (90-135) Gt carbon. By documenting the environmental correlates of root biomass allocation, our results can inform model projections of global vegetation dynamics under current and future climate scenarios.

Building a global database of soil microbial biomass and function: a call for collaboration

Global analyses are emerging as valuable complements to local and regional scale studies in ecology and are useful for examining many of the major environmental issues that we face today. Soil ecology has significantly benefited from these developments, with recent syntheses unearthing interesting, unexpected biogeographic patterns in belowground biotic communities. However, some questions still remain unanswered, and the accuracy of these studies is inevitably limited by the extent of the data they draw upon. This is a particular problem in global ecology because most datasets used exhibit geographic bias in sample distribution. Here, we work towards addressing this problem with an open call for collaboration on a planned global analysis of soil phospholipid fatty acid and potential enzyme activity measurements. We summarize the current extent of our dataset, outline the planned analyses, and provide information for prospective collaborators who would like to contribute or learn more.

The Genome of Peronospora belbahrii Reveals High Heterozygosity, a Low Number of Canonical Effectors, and TC-Rich Promoters

Along with Plasmopara destructor, Peronosopora belbahrii has arguably been the economically most important newly emerging downy mildew pathogen of the past two decades. Originating from Africa, it has started devastating basil production throughout the world, most likely due to the distribution of infested seed material. Here, we present the genome of this pathogen and results from comparisons of its genomic features to other oomycetes. The assembly of the nuclear genome was around 35.4 Mbp in length, with an N₅₀ scaffold length of around 248 kbp and an L₅₀ scaffold count of 46. The circular mitochondrial genome consisted of around 40.1 kbp. From the repeat-masked genome, 9,049 protein-coding genes were predicted, out of which 335 were predicted to have extracellular functions, representing the smallest secretome so far found in peronosporalean oomycetes. About 16% of the genome consists of repetitive sequences, and, based on simple sequence repeat regions, we provide a set of microsatellites that could be used for population genetic studies of P. belbahrii. P. belbahrii has undergone a high degree of convergent evolution with other obligate parasitic pathogen groups, reflecting its obligate biotrophic lifestyle. Features of its secretome, signaling networks, and promoters are presented, and some patterns are hypothesized to reflect the high degree of host specificity in Peronospora species. In addition, we suggest the presence of additional virulence factors apart from classical effector classes that are promising candidates for future functional studies.

A global database of soil nematode abundance and functional group composition

As the most abundant animals on earth, nematodes are a dominant component of the soil community. They play critical roles in regulating biogeochemical cycles and vegetation dynamics within and across landscapes and are an indicator of soil biological activity. Here, we present a comprehensive global dataset of soil nematode abundance and functional group composition. This dataset includes 6,825 georeferenced soil samples from all continents and biomes. For geospatial mapping purposes these samples are aggregated into 1,933 unique 1-km pixels, each of which is linked to 73 global environmental covariate data layers. Altogether, this dataset can help to gain insight into the spatial distribution patterns of soil nematode abundance and community composition, and the environmental drivers shaping these patterns.

Global distribution of earthworm diversity

Soil organisms, including earthworms, are a key component of terrestrial ecosystems. However, little is known about their diversity, their distribution, and the threats affecting them. We compiled a global dataset of sampled earthworm communities from 6928 sites in 57 countries as a basis for predicting patterns in earthworm diversity, abundance, and biomass. We found that local species richness and abundance typically peaked at higher latitudes, displaying patterns opposite to those observed in aboveground organisms. However, high species dissimilarity across tropical locations may cause diversity across the entirety of the tropics to be higher than elsewhere. Climate variables were found to be more important in shaping earthworm communities than soil properties or habitat cover. These findings suggest that climate change may have serious implications for earthworm communities and for the functions they provide.

Understanding climate change from a global analysis of city analogues

Combating climate change requires unified action across all sectors of society. However, this collective action is precluded by the ‘consensus gap’ between scientific knowledge and public opinion. Here, we test the extent to which the iconic cities around the world are likely to shift in response to climate change. By analyzing city pairs for 520 major cities of the world, we test if their climate in 2050 will resemble more closely to their own current climate conditions or to the current conditions of other cities in different bioclimatic regions. Even under an optimistic climate scenario (RCP 4.5), we found that 77% of future cities are very likely to experience a climate that is closer to that of another existing city than to its own current climate. In addition, 22% of cities will experience climate conditions that are not currently experienced by any existing major cities. As a general trend, we found that all the cities tend to shift towards the sub-tropics, with cities from the Northern hemisphere shifting to warmer conditions, on average ~1000 km south (velocity ~20 km.year^-1), and cities from the tropics shifting to drier conditions. We notably predict that Madrid’s climate in 2050 will resemble Marrakech’s climate today, Stockholm will resemble Budapest, London to Barcelona, Moscow to Sofia, Seattle to San Francisco, Tokyo to Changsha. Our approach illustrates how complex climate data can be packaged to provide tangible information. The global assessment of city analogues can facilitate the understanding of climate change at a global level but also help land managers and city planners to visualize the climate futures of their respective cities, which can facilitate effective decision-making in response to on-going climate change.

Comparative genomics of the nonlegume Parasponia reveals insights into evolution of nitrogen-fixing rhizobium symbioses

Nodules harboring nitrogen-fixing rhizobia are a well-known trait of legumes, but nodules also occur in other plant lineages, with rhizobia or the actinomycete Frankia as microsymbiont. It is generally assumed that nodulation evolved independently multiple times. However, molecular-genetic support for this hypothesis is lacking, as the genetic changes underlying nodule evolution remain elusive. We conducted genetic and comparative genomics studies by using Parasponia species (Cannabaceae), the only nonlegumes that can establish nitrogen-fixing nodules with rhizobium. Intergeneric crosses between Parasponia andersonii and its nonnodulating relative Trema tomentosa demonstrated that nodule organogenesis, but not intracellular infection, is a dominant genetic trait. Comparative transcriptomics of P. andersonii and the legume Medicago truncatula revealed utilization of at least 290 orthologous symbiosis genes in nodules. Among these are key genes that, in legumes, are essential for nodulation, including NODULE INCEPTION (NIN) and RHIZOBIUM-DIRECTED POLAR GROWTH (RPG). Comparative analysis of genomes from three Parasponia species and related nonnodulating plant species show evidence of parallel loss in nonnodulating species of putative orthologs of NIN, RPG, and NOD FACTOR PERCEPTION Parallel loss of these symbiosis genes indicates that these nonnodulating lineages lost the potential to nodulate. Taken together, our results challenge the view that nodulation evolved in parallel and raises the possibility that nodulation originated ∼100 Mya in a common ancestor of all nodulating plant species, but was subsequently lost in many descendant lineages. This will have profound implications for translational approaches aimed at engineering nitrogen-fixing nodules in crop plants.