Diversity-stability Relationships Research Articles

Larger nations benefit more than smaller nations from the stabilizing effects of crop diversity.

Crop diversification is increasingly recognized as a strategy to stabilize national food production, yet the benefits of this approach may vary across nations due to the scale dependence of crop diversity and stability. Here we use crop production data from 131 nations from 1961 to 2020 to explore the spatial scale dependence of the crop diversity-stability relationship. Drawing on ecological theory and complementary analytical approaches, we find that as the total national harvested area increases, yield stability increases. Crop diversity stabilizes national yield stability, as does an increase in the number of farms, but these stabilizing effects are weaker in smaller countries. Our findings suggest that enhancing crop diversity at the national level may not provide a de facto universal strategy for increasing yield stability across all countries-with implications for national strategies promoting crop diversification to protect against food system shocks.

A marine heatwave changes the stabilizing effects of biodiversity in kelp forests.

Biodiversity can stabilize ecological communities through biological insurance, but climate and other environmental changes may disrupt this process via simultaneous ecosystem destabilization and biodiversity loss. While changes to diversity-stability relationships (DSRs) and the underlying mechanisms have been extensively explored in terrestrial plant communities, this topic remains largely unexplored in benthic marine ecosystems that comprise diverse assemblages of producers and consumers. By analyzing two decades of kelp forest biodiversity survey data, we discovered changes in diversity, stability, and their relationships at multiple scales (biological organizational levels, spatial scales, and functional groups) that were linked with the most severe marine heatwave ever documented in the North Pacific Ocean. Moreover, changes in the strength of DSRs during/after the heatwave were more apparent among functional groups than both biological organizational levels (population vs. ecosystem levels) and spatial scales (local vs. broad scales). Specifically, the strength of DSRs decreased for fishes, increased for mobile invertebrates and understory algae, and were unchanged for sessile invertebrates during/after the heatwave. Our findings suggest that biodiversity plays a key role in stabilizing marine ecosystems, but the resilience of DSRs to adverse climate impacts primarily depends on the functional identities of ecological communities.

Appropriate livestock grazing alleviates the loss of plant diversity and maintains community resistance in alpine meadows

Alpine meadows constitute one of the major ecosystems on the Qinghai-Tibetan Plateau, with livestock grazing exerting a considerable impact on their biodiversity. However, the degree to which plant diversity influences community stability under different grazing intensities remains unclear in this region. This study conducted controlled grazing experiments across four levels of grazing intensity (no-, low-, medium-, and high-grazing) based on herbage utilization rate to assess the influence of grazing intensities on plant community structure and diversity-stability relationships. We discovered that high-grazing reduced plant diversity and attenuated the temporal stability and resistance of above-ground biomass. No- and low-grazing could alleviate plant biomass loss, with community resistance being optimal under low-grazing. The direct effects of livestock grazing on temporal stability were found to be negligible. Plant characteristics and diversity accounted for a substantial proportion of livestock grazing effects on community resistance (R2 = 0.46), as revealed by piecewise structural equation model analysis. The presence of plant diversity enhances the resistance of alpine meadows against disturbance and accelerates the recovery after grazing. Our results suggest that low-grazing intensity may represent a judicious option for preserving species diversity and community stability on the Qinghai-Tibetan Plateau.

Growing-season drought and nitrogen addition interactively impair grassland ecosystem stability by reducing species diversity, asynchrony, and stability

Aboveground net primary productivity controls the amount of energy available to sustain all living organisms, and its sustainable provision relies on the stability of grassland ecosystems. Human activities leading to global changes, such as increased nitrogen (N) deposition and the more frequent occurrence of extreme precipitation events, with N addition increasing the sensitivity of ecosystem production stability to changes in the precipitation regime. However, whether N addition, in combination with seasonal precipitation increases or severe drought, affects ecosystem stability remains unclear. In this study, we conducted a six-year environmental change monitoring experiment in a semiarid grassland in northern China to test the effects of N addition, seasonal drought, and precipitation increases on the temporal stability of ecosystem productivity. Our study revealed that an interaction between drought and N addition reduced species diversity, species asynchrony, species stability, and thus ecosystem stability. These environmental change drivers (except for precipitation increase) induced a positive relationship between species asynchrony and diversity, whereas N addition interactively with drought and precipitation increase led to a negative relationship between diversity and species stability. Only N addition interactively with drought induced a positive species diversity–ecosystem stability relationship because lower species stability was overcome by increased species asynchrony. Our study is great importance to illustrate that production temporal stability tends to be inhibited with drought, though interactively with nutrient N addition. These findings highlight the primary role of asynchronous dynamics among species in modulating the effects of environmental change on diversity-stability relationships.

Response trait diversity and species asynchrony underlie the diversity-stability relationship in Romanian bird communities.

Biodiversity-stability relationships have frequently been studied in ecology, with the recent integration of traits to explain community stability over time. Classical theory underlying the biodiversity-stability relationship posits that different species' responses to the environment should stabilise community-level properties (e.g. biomass or abundance) through compensatory dynamics. However, functional response traits, which aim to predict how species respond to environmental change, are still rarely integrated into studies of ecological stability. Such traits should mechanistically drive community stability, both in terms of community abundance (functional variability) and composition (compositional variability). In turn, whether and how functional or compositional stability scales to affect temporal variation in functional effect traits (a proxy for ecosystem functioning) remains largely unknown, but is key to consistent ecosystem functioning under environmental change. Here, we explore the diversity-stability relationship in bird communities using annual survey data across 98 sites in central Romania, in combination with global trait databases and structural equation models. We show that higher response trait diversity promotes compositional variability directly, and functional variability indirectly via species asynchrony. In turn, functional variability impacts the temporal stability of effect trait diversity. Multiple facets of diversity and community stability differ between natural forests and agricultural or human-dominated survey sites, and the relationship between response diversity and functional variability is mediated by land cover. Further integration of response-and-effect trait frameworks into studies of community stability will enhance understanding of the drivers of biodiversity change, allowing targeted conservation decision-making with a focus on stable ecosystem functioning in the face of global environmental change.

Scale-dependent changes in ecosystem temporal stability over six decades of succession.

A widely assumed, but largely untested, tenet in ecology is that ecosystem stability tends to increase over succession. We rigorously test this idea using 60-year continuous data of old field succession across 480 plots nested within 10 fields. We found that ecosystem temporal stability increased over succession at the larger field scale (γ stability) but not at the local plot scale (α stability). Increased spatial asynchrony among plots within fields increased γ stability, while temporal increases in species stability and decreases in species asynchrony offset each other, resulting in no increase in α stability at the local scale. Furthermore, we found a notable positive diversity-stability relationship at the larger but not local scale, with the increased γ stability at the larger scale associated with increasing functional diversity later in succession. Our results emphasize the importance of spatial scale in assessing ecosystem stability over time and how it relates to biodiversity.

Diversity-stability relationships across organism groups and ecosystem types become decoupled across spatial scales.

The relationship between biodiversity and stability, or its inverse, temporal variability, is multidimensional and complex. Temporal variability in aggregate properties, like total biomass or abundance, is typically lower in communities with higher species diversity (i.e., the diversity-stability relationship [DSR]). At broader spatial extents, regional-scale aggregate variability is also lower with higher regional diversity (in plant systems) and with lower spatial synchrony. However, focusing exclusively on aggregate properties of communities may overlook potentially destabilizing compositional shifts. It is not yet clear how diversity is related to different components of variability across spatial scales, nor whether regional DSRs emerge across a broad range of organisms and ecosystem types. To test these questions, we compiled a large collection of long-term metacommunity data spanning a wide range of taxonomic groups (e.g., birds, fish, plants, invertebrates) and ecosystem types (e.g., deserts, forests, oceans). We applied a newly developed quantitative framework for jointly analyzing aggregate and compositional variability across scales. We quantified DSRs for composition and aggregate variability in local communities and metacommunities. At the local scale, more diverse communities were less variable, but this effect was stronger for aggregate than compositional properties. We found no stabilizing effect of γ-diversity on metacommunity variability, but β-diversity played a strong role in reducing compositional spatial synchrony, which reduced regional variability. Spatial synchrony differed among taxa, suggesting differences in stabilization by spatial processes. However, metacommunity variability was more strongly driven by local variability than by spatial synchrony. Across a broader range of taxa, our results suggest that high γ-diversity does not consistently stabilize aggregate properties at regional scales without sufficient spatial β-diversity to reduce spatial synchrony.

Environmental conditions are the dominant factor influencing stability of terrestrial ecosystems on the Tibetan plateau

Ecosystem stability is essential for the sustainable provision of diverse ecosystem services. However, the factors that maintain ecosystem stability and their relative importance on the Tibetan Plateau, a region sensitive to climate change, remain unclear. Here, we combined data from ground-based biodiversity surveys at 143 sites from 2019 to 2021 with the temporal stability of ecosystems derived from remote sensing data from 2000 to 2020 to disentangle mechanisms of diversity–stability relationships. We further quantified the impact of biodiversity (taxonomic, functional, and phylogenetic diversity) and environmental context (spatial location, climate, and soil conditions) on temporal stability. Our results show that the stability of a typical ecosystem on the Tibetan Plateau is mainly regulated by environmental factors, and the environmental context can directly affect the stability of the ecosystem rather than indirectly through biodiversity. These findings are critical for adaptation measures and prioritizing conservation areas for future climate change scenarios.

Species asynchrony stabilizes productivity over 20 years in Northeast China.

The stability of forest productivity can reflect the functioning of forest ecosystems. It is a crucial topic to understand the relationship between biodiversity and ecosystem functions in ecology. Although previous studies have made great progress in understanding the effects of diversity, species asynchrony, and other factors on community biomass and productivity, few studies have explored how these factors affect the temporal stability of productivity. In this study, we hypothesized that diversity, species asynchrony, and topography would directly or indirectly impact the temporal stability of productivity. To test this hypothesis, we used a multiple regression model and a piecewise structural equation model based on the inventory data over 20 years (5-year intervals) from 1992 to 2012 at Jingouling Forest Farm in Northeast China. Our results showed that species asynchrony was the main driving factor affecting the temporal stability of productivity. Structural diversity significantly decreased community stability, while species diversity had a nonsignificant effect on it. We found the combination of a multiple regression model and a piecewise structural equation model is an effective method for evaluating the factors that influence community stability. The effect of species asynchrony is crucial for understanding the ecological mechanisms underlying the diversity-stability relationship in mixed forests.

Mechanisms underpinning community stability along a latitudinal gradient: Insights from a niche-based approach.

At large scales, the mechanisms underpinning stability in natural communities may vary in importance due to changes in species composition, mean abundance, and species richness. Here we link species characteristics (niche positions) and community characteristics (richness and abundance) to evaluate the importance of stability mechanisms in 156 butterfly communities monitored across three European countries and spanning five bioclimatic regions. We construct niche-based hierarchical structural Bayesian models to explain first differences in abundance, population stability, and species richness between the countries, and then explore how these factors impact community stability both directly and indirectly (via synchrony and population stability). Species richness was partially explained by the position of a site relative to the niches of the species pool, and species near the centre of their niche had higher average population stability. The differences in mean abundance, population stability, and species richness then influenced how much variation in community stability they explained across the countries. We found, using variance partitioning, that community stability in Finnish communities was most influenced by community abundance, whereas this aspect was unimportant in Spain with species synchrony explaining most variation; the UK was somewhat intermediate with both factors explaining variation. Across all countries, the diversity-stability relationship was indirect with species richness reducing synchrony which increased community stability, with no direct effects of species richness. Our results suggest that in natural communities, biogeographical variation observed in key drivers of stability, such as population abundance and species richness, leads to community stability being limited by different factors and that this can partially be explained due to the niche characteristics of the European butterfly assemblage.

How to measure response diversity

Abstract The insurance effect of biodiversity—that diversity stabilises aggregate ecosystem properties—is mechanistically underlain by inter‐ and intraspecific trait variation in organismal responses to the environment. This variation, termed response diversity, is therefore a potentially critical determinant of ecological stability. However, response diversity has yet to be widely quantified, possibly due to difficulties in its measurement. Even when it has been measured, approaches have varied. Here, we review methods for measuring response diversity and from them distil a methodological framework for quantifying response diversity from experimental and/or observational data, which can be practically applied in laboratory and field settings across a range of taxa. Previous empirical studies on response diversity most commonly invoke response traits as proxies aimed at capturing species' ecological responses to the environment. Our approach, which is based on environment‐dependent ecological responses to any biotic or abiotic environmental variable, is conceptually simple and robust to any form of environmental response, including nonlinear responses. Given its derivation from empirical data on species' ecological responses, this approach should more directly reflect response diversity than the trait‐based approach dominant in the literature. By capturing even subtle inter‐ or intraspecific variation in environmental responses, and environment dependencies in response diversity, we hope this framework will motivate tests of the diversity–stability relationship from a new perspective, and provide an approach for mapping, monitoring and conserving this critical dimension of biodiversity.

Assortative mating for reproductive timing affects population recruitment and resilience in a quantitative genetic model.

Quantitative models that simulate the inheritance and evolution of fitness-linked traits offer a method for predicting how environmental or anthropogenic perturbations can affect the dynamics of wild populations. Random mating between individuals within populations is a key assumption of many such models used in conservation and management to predict the impacts of proposed management or conservation actions. However, recent evidence suggests that non-random mating may be underestimated in wild populations and play an important role in diversity-stability relationships. Here we introduce a novel individual-based quantitative genetic model that incorporates assortative mating for reproductive timing, a defining attribute of many aggregate breeding species. We demonstrate the utility of this framework by simulating a generalized salmonid lifecycle, varying input parameters, and comparing model outputs to theoretical expectations for several eco-evolutionary, population dynamic scenarios. Simulations with assortative mating systems resulted in more resilient and productive populations than those that were randomly mating. In accordance with established ecological and evolutionary theory, we also found that decreasing the magnitude of trait correlations, environmental variability, and strength of selection each had a positive effect on population growth. Our model is constructed in a modular framework so that future components can be easily added to address pressing issues such as the effects of supportive breeding, variable age structure, differential selection by sex or age, and fishery interactions on population growth and resilience. With code published in a public Github repository, model outputs may easily be tailored to specific study systems by parameterizing with empirically generated values from long-term ecological monitoring programs.

Species portfolio effects dominate seasonal zooplankton stabilization within a large temperate lake.

Portfolio effects (PEs) in ecology refer to the suite of phenomenon where the temporal variation of aggregate ecosystem properties (i.e., abundance) is lower than that of their ecosystem components. An example of this is where differential responses of species to environmental variation generate stability at higher levels of ecological organization (e.g., local community, metapopulation, metacommunity). Most of the research examining such PEs has focused on spatial or interannual variation of ecosystems; however, as global change continues to alter seasonality and ecosystem functioning, understanding the underlying food web structures that help maintain stability at multiple spatial and temporal scales is critical to managing ecological systems. Recent advances investigating diversity-stability relationships has led to the development of frameworks that incorporate a metacommunity perspective which allows for the partitioning of PEs across organizational scales (i.e., local community, metapopulation, cross-community, metacommunity) from local population dynamics (total). This partitioning yields insights into the mechanisms that generate observed PEs in nature. Here, we employed one of these recently developed frameworks on a temporally (1986-1999, 2008-2019) and spatially (five sampling stations, local communities) extensive data set of zooplankton abundance (e.g., density) within a large temperate lake to investigate how temporal (seasonal) and spatial (among site) PEs influence stability within the zooplankton metacommunity. We found that seasonal asynchrony of different zooplankton species within local communities and across communities generated the vast majority of stabilization, while spatial (i.e., metapopulation) dynamics were more synchronous and contributed little to overall system stability. Furthermore, significantly positive diversity-asynchrony relationships at the total, local- and cross-community scales were found as asynchrony was positively correlated with local Shannon diversity. Last, a comparison of PEs over the time periods, during which significant local and global changes (i.e., climate warming, invasive species) have occurred suggests that PEs may be eroding, as increasingly synchronous dynamics and declining diversity in recent years have led to a rise in metacommunity variability. We end by arguing for the critical importance of understanding seasonally driven stabilizing mechanisms as local and global changes threaten to fundamentally alter seasonal signals with potentially strong implications for the structures that lend stability to ecosystems.

Biodiversity stabilizes plant communities through statistical-averaging effects rather than compensatory dynamics

Understanding the relationship between biodiversity and ecosystem stability is a central goal of ecologists. Recent studies have concluded that biodiversity increases community temporal stability by increasing the asynchrony between the dynamics of different species. Theoretically, this enhancement can occur through either increased between-species compensatory dynamics, a fundamentally biological mechanism; or through an averaging effect, primarily a statistical mechanism. Yet it remains unclear which mechanism is dominant in explaining the diversity-stability relationship. We address this issue by mathematically decomposing asynchrony into components separately quantifying the compensatory and statistical-averaging effects. We applied the new decomposition approach to plant survey and experimental data from North American grasslands. We show that statistical averaging, rather than compensatory dynamics, was the principal mediator of biodiversity effects on community stability. Our simple decomposition approach helps integrate concepts of stability, asynchrony, statistical averaging, and compensatory dynamics, and suggests that statistical averaging, rather than compensatory dynamics, is the primary means by which biodiversity confers ecological stability.

Aridity determines the effects of warming on community stability in Inner Mongolian grassland

Ongoing climate warming affects biodiversity and ecosystem functions and services. However, the mechanisms underlying the effects of climate warming on biodiversity-stability relationships are unclear. Whether these effects depend on environmental factors such as aridity remains an open question. Here, we conducted a factorial field experiment to investigate the effects of warming on nine grasslands in a transect spanning ∼1500 km with increasing aridity in Inner Mongolia, China. Our results showed that warming-induced decrease in species richness did not alter community variability across nine grasslands. However, species richness was negatively related to community variability under warming conditions, but not under control conditions. This indicated that warming regulated the species diversity-stability relationship. Warming increased species synchrony, species variability, and dominant species variability. These outcomes lead to increases in temporal variability of community productivity. Climate fluctuation increased species variability and species synchrony within communities, thereby reducing community stability. Importantly, the size of warming effects on community variability via biotic mechanisms grew with increasing aridity. Our results indicate that species diversity plays a crucial role in regulating community stability under climate warming. Our findings reveal the generality of biotic mechanisms responsible for grassland stability under future climate warming, and highlight that the geographic variation in aridity determines the effect sizes of warming on grassland stability. We suggest that ongoing and future climate warming might disrupt grasslands stability with high aridity, thereby reducing ecosystem functions and services.

Examining the diversity, stability and functioning of marine fish communities across a latitudinal gradient

AbstractAimAs anthropogenic stressors on the biosphere intensify, understanding how communities respond to disturbances is critical. Biodiversity is often thought to promote the stability of communities over time and enhance ecosystem functioning. However, results have been inconsistent, and the multifaceted linkages among diversity, stability and functioning under acute disturbances remain poorly understood. We experimentally tested the responses of marine fish communities to disturbance (i.e., acute habitat loss) across a diversity gradient spanning 35 degrees of latitude in the western Atlantic Ocean to assess the diversity–stability relationship and the interplay between diversity, stability, and fish biomass recovery (as a proxy for function) in marine fish communities.LocationWestern Atlantic Ocean [Maine, Massachusetts, North Carolina, Florida (USA), Belize and Panama].Time period2016–2017.Major taxa studiedSmall, bottom‐dwelling (‘cryptobenthic’) fishes.MethodsWe experimentally tested the response of marine fish communities to disturbance across a diversity gradient on human‐made dock pilings. We holistically sampled cryptobenthic fish communities, then we imposed a severe disturbance by removing all benthic epifauna. We then compared the community stability, defined as the constancy in community composition, on disturbed and undisturbed pilings after one year.ResultsDiversity showed a negative effect on community stability at both the regional (across docks) and local (within docks) scales. Similarly, local diversity was negatively correlated with ecosystem function. These effects are exacerbated by the habitat loss imposed via our experimental treatment.Main conclusionsOur results suggest that habitat loss may re‐shuffle diverse, tropical communities more intensively than species‐poor, temperate communities, which impacts biomass recovery, our proxy of functioning. Contrary to ecological theory, in small‐bodied, benthos‐associated vertebrate communities, biodiversity may neither promote stability nor functioning, suggesting that human disturbances may be particularly impactful in tropical, high‐diversity ecosystems.

Tree diversity depending on environmental gradients promotes biomass stability via species asynchrony in China’s forest ecosystems

There is mounting evidence that biodiversity promotes ecological stability in changing environments. However, understanding diversity–stability relationships and their underlying mechanisms across large-scale tree diversity and natural environmental gradients are still controversial and largely lacking. We used thirty-nine 0.12 ha long-term permanent forest plots spanning China’s various forest types to test the effects of multiple abiotic (climate, soil, age and topography) and biotic factors (taxonomic and structural diversity, functional diversity and community-mean traits, and species asynchrony) on biomass stability and its components (mean biomass and biomass variability) over time. We used multiple analytical methods to identify the best explanatory variables and complicated causal relationships for community biomass stability. Our results showed that species richness increased biomass stability by promoting species asynchrony. Structural and functional diversity had a weaker effect on biomass stability. Forest age and structural diversity increased mean biomass and biomass variability significantly and simultaneously. Communities dominated by tree species with high wood density had high biomass stability. Soil nitrogen enhanced biomass stability directly and indirectly through its effects on mean biomass. Soil nitrogen to phosphorus ratio increased biomass stability via increasing species asynchrony. Precipitation indirectly increased biomass stability by affecting tree diversity. Moreover, the direct and indirect effects of soil nutrients on biomass stability were greater than that of climate variables. Our results suggest that species asynchrony is the main mechanism proposed to explain the stabilizing effect of diversity on community biomass, supporting two mechanisms, namely, the biodiversity insurance hypothesis and complementary dynamics. Soil and climate factors also play an important role in shaping diversity–stability relationships. Our results provide a new insight into how tree diversity affects ecosystem stability across diverse community types and large-scale environmental gradients.

Bacterial Recruitment to Carnivorous Pitcher Plant Communities: Identifying Sources Influencing Plant Microbiome Composition and Function.

Processes influencing recruitment of diverse bacteria to plant microbiomes remain poorly understood. In the carnivorous pitcher plant Sarracenia purpurea model system, individual pitchers open to collect rainwater, invertebrates and a diverse microbial community, and this detrital food web is sustained by captured insect prey. This study examined how potential sources of bacteria affect the development of the bacterial community within pitchers, how the host plant tissue affects community development and how established vs. assembling communities differ. In a controlled greenhouse experiment, seven replicate pitchers were allocated to five treatments to exclude specific bacterial sources or host tissue: milliQ water only, milliQ + insect prey, rainwater + prey, established communities + prey, artificial pitchers with milliQ + prey. Community composition and functions were examined over 8–40 weeks using bacterial gene sequencing and functional predictions, measurements of cell abundance, hydrolytic enzyme activity and nutrient transformations. Distinct community composition and functional differences between artificial and real pitchers confirm an important influence of host plant tissue on community development, but also suggest this could be partially related to host nutrient uptake. Significant recruitment of bacteria to pitchers from air was evident from many taxa common to all treatments, overlap in composition between milliQ, milliQ + prey, and rainwater + prey treatments, and few taxa unique to milliQ only pitchers. Community functions measured as hydrolytic enzyme (chitinase, protease) activity suggested a strong influence of insect prey additions and were linked to rapid transformation of insect nutrients into dissolved and inorganic sources. Bacterial taxa found in 6 of 7 replicate pitchers within treatments, the “core microbiome” showed tighter successional trajectories over 8 weeks than all taxa. Established pitcher community composition was more stable over 8 weeks, suggesting a diversity-stability relationship and effect of microinvertebrates on bacteria. This study broadly demonstrates that bacterial composition in host pitcher plants is related to both stochastic and specific bacterial recruitment and host plants influence microbial selection and support microbiomes through capture of insect prey.

High-order correlations in species interactions lead to complex diversity-stability relationships for ecosystems.

How ecosystems maintain stability is an active area of research. Inspired by applications of random matrix theory in nuclear physics, May showed decades ago that in an ecosystem model with many randomly interacting species, increasing species diversity decreases the stability of the ecosystem. There have since been many additions to May's efforts, one being an improved understanding the effect of mutualistic, competitive, or predator-prey-like correlations between pairs of species. Here we extend a random matrix technique developed in the context of spin-glass theory to study the effect of high-order correlations among species interactions. The resulting analytically solvable models include next-to-nearest-neighbor correlations in the species interaction network, such as the enemy of my enemy is my friend, as well as higher-order correlations. We find qualitative differences from May and others' models, including nonmonotonic diversity-stability relationships. Furthermore, inclusion of particular next-to-nearest-neighbor correlations in predator-prey as opposed to mutualist-competitive networks causes the former to transition to being more stable at higher species diversity. We discuss potential applicability of our results to microbiota engineering and to the ecology of interpredator interactions, such as cub predation between lions and hyenas as well as companionship between humans and dogs.

Crop diversity effects on temporal agricultural production stability across European regions

Stabilizing agricultural production is fundamental to food security. At the national level, increasing the effective diversity of cultivated crops has been found to increase temporal production stability, i.e., the year-to-year stability of total caloric production of all crops combined. Here, we specifically investigated these effects at the regional level for the European Union and tested the effect of crop diversity in relation to agricultural inputs, soil properties, climate instability, and time on caloric, protein, and fat stability, as we hypothesized that the effect of diversity is context dependent. We further investigated these relationships for specific countries. We found that greater crop diversity was consistently associated with an increase in production stability, particularly in regions with large areas equipped for irrigation and low soil type diversity. For instance, in Spain and Italy, crop diversity showed the strongest positive effect among all predictors, while on the European level, the stabilizing effect of nitrogen use was substantially higher. In Germany, the crop diversity-stability relationship was weak, suggesting that crops react similarly to climatic, economic, and political factors or are grown in the same periods. With this study, we substantiate previous findings that crop diversity stabilizes agricultural caloric production and extend these with regard to protein and fat. The results elucidate the key drivers that enhance production stability for different European countries and regions, which is of key importance for a comparably productive agricultural region like Europe.