Experiments suggest that biodiversity enhances the ability of ecosystems to maintain multiple functions, such as carbon storage, productivity, and the buildup of nutrient pools (multifunctionality). However, the relationship between biodiversity and multifunctionality has never been assessed globally in natural ecosystems. We report here on a global empirical study relating plant species richness and abiotic factors to multifunctionality in drylands, which collectively cover 41% of Earth's land surface and support over 38% of the human population. Multifunctionality was positively and significantly related to species richness. The best-fitting models accounted for over 55% of the variation in multifunctionality and always included species richness as a predictor variable. Our results suggest that the preservation of plant biodiversity is crucial to buffer negative effects of climate change and desertification in drylands.

The biogeochemical cycles of carbon (C), nitrogen (N) and phosphorus (P) are interlinked by primary production, respiration and decomposition in terrestrial ecosystems. It has been suggested that the C, N and P cycles could become uncoupled under rapid climate change because of the different degrees of control exerted on the supply of these elements by biological and geochemical processes. Climatic controls on biogeochemical cycles are particularly relevant in arid, semi-arid and dry sub-humid ecosystems (drylands) because their biological activity is mainly driven by water availability. The increase in aridity predicted for the twenty-first century in many drylands worldwide may therefore threaten the balance between these cycles, differentially affecting the availability of essential nutrients. Here we evaluate how aridity affects the balance between C, N and P in soils collected from 224 dryland sites from all continents except Antarctica. We find a negative effect of aridity on the concentration of soil organic C and total N, but a positive effect on the concentration of inorganic P. Aridity is negatively related to plant cover, which may favour the dominance of physical processes such as rock weathering, a major source of P to ecosystems, over biological processes that provide more C and N, such as litter decomposition. Our findings suggest that any predicted increase in aridity with climate change will probably reduce the concentrations of N and C in global drylands, but increase that of P. These changes would uncouple the C, N and P cycles in drylands and could negatively affect the provision of key services provided by these ecosystems.

Summary 1.Habitat filtering (HF, trait convergence) and niche differentiation (ND, trait divergence) are known to impact upon plant community structure. Both processes integrate individual responses to the abiotic environment and biotic interactions. Thus, it is difficult to clearly identify the underlying abiotic and biotic factors that ultimately impact community structure by looking at community-level patterns of trait divergence or convergence alone. 2. We used a functional trait-based and multiscale approach to assess how biotic interactions and aridity determine the functional structure of semi-arid shrublands sampled along a large aridity gradient in Spain. At the regional scale, we investigated functional differences among species (axes of specialization) to identify important traits for community assembly. At the community scale, we evaluated the relative impact of HF and ND on community structure using a null model approach. Finally, at the plant neighbourhood scale, we evaluated the impact of biotic interactions on community structure by investigating the spatial patterns of trait aggregation. 3. The shrub species surveyed can be separated along four axes of specialization based on their above-ground architecture and leaf morphology. Our community scale analysis suggested that the functional structure of semi-arid communities was clearly non-random, HF and ND acting independently on different traits to determine community structure along the aridity gradient. At the plant neighbourhood scale, the spatial distribution of species was also clearly not random, suggesting that competition and facilitation impacted on the observed changes in the functional diversity of shrubland communities along the aridity gradient. 4. Synthesis: Our results demonstrated that HF and ND acted simultaneously on independent traits to jointly determine community structure. Most importantly, our multiscale approach suggested that competition and facilitation interplayed with aridity to determine this structure. Competition appeared to be constant along the aridity gradient and explained the high functional diversity observed in semi-arid shrublands. Facilitation affected subordinate and rare species and, thus, may act to enhance the biodiversity of these ecosystems. Finally, the framework employed in our study allows moving forward from the examination of patterns to the development of mechanistic traitbased approaches to study plant community assembly.

The stability of ecological communities is critical for the stable provisioning of ecosystem services, such as food and forage production, carbon sequestration, and soil fertility. Greater biodiversity is expected to enhance stability across years by decreasing synchrony among species, but the drivers of stability in nature remain poorly resolved. Our analysis of time series from 79 datasets across the world showed that stability was associated more strongly with the degree of synchrony among dominant species than with species richness. The relatively weak influence of species richness is consistent with theory predicting that the effect of richness on stability weakens when synchrony is higher than expected under random fluctuations, which was the case in most communities. Land management, nutrient addition, and climate change treatments had relatively weak and varying effects on stability, modifying how species richness, synchrony, and stability interact. Our results demonstrate the prevalence of biotic drivers on ecosystem stability, with the potential for environmental drivers to alter the intricate relationship among richness, synchrony, and stability.