Individual species provide multifaceted contributions to the stability of ecosystems

Existing levels of biodiversity and river location may determine changes from small hydropower developments

Global ecosystems are facing anthropogenic threats that affect their ecological functions and biodiversity. However, we still lack an understanding of how biodiversity can mediate the responses of ecosystems or communities to human disturbance across spatial gradients. Here, we examined how existing, spatial patterns of biodiversity influence the ecological effects of small hydropower plants (SHPs) on macroinvertebrates in river ecosystems. This study found that levels of biodiversity (e.g., number of species) can influence the degrees of its alterations by SHPs occurring along elevational gradients. The results of the study reveal that the construction of SHPs has various effects on biodiversity. For example, low-altitude areas with low biodiversity (species richness less than 12) showed a small increase in biodiversity compared to high-altitude areas (species richness more than 12) under SHP disturbances. The increases in the effective habitat area of the river segment could be a driver of the enhanced biodiversity in response to SHP effects. Changes in the numerically dominant species contributed to the overall level of community variation from disturbances. Location-specific strategies may mitigate the effects of SHPs and perhaps other disturbances.

Opposing effects of warming on the stability of above- and belowground productivity in facing an extreme drought event

Climate warming, often accompanied by extreme drought events, could have profound effects on both plant community structure and ecosystem functioning. However, how warming interacts with extreme drought to affect community- and ecosystem-level stability remains a largely open question. Using data from a manipulative experiment with three warming treatments in an alpine meadow that experienced one extreme drought event, we investigated how warming modulates resistance and recovery of community structural and ecosystem functional stability in facing with extreme drought. We found warming decreased resistance and recovery of aboveground net primary productivity (ANPP) and structural resistance but increased resistance and recovery of belowground net primary productivity (BNPP), overall net primary productivity (NPP), and structural recovery. The findings highlight the importance of jointly considering above- and belowground processes when evaluating ecosystem stability under global warming and extreme climate events. The stability of dominant species, rather than species richness and species asynchrony, was identified as a key predictor of ecosystem functional resistance and recovery, except for BNPP recovery. In addition, structural resistance of common species contributed strongly to the resistance changes in BNPP and NPP. Importantly, community structural resistance and recovery dominated the resistance and recovery of BNPP and NPP, but not for ANPP, suggesting the different mechanisms underlie the maintenance of stability of above- versus belowground productivity. This study is among the first to explain that warming modulates ecosystem stability in the face of extreme drought and lay stress on the need to investigate ecological stability at the community level for a more mechanistic understanding of ecosystem stability in response to climate extremes.

Breakdown in seasonal dynamics of subtropical ant communities with land-cover change

Concerns about widespread human-induced declines in insect populations are mounting, yet little is known about how land-use change modifies both the trends and variability of insect communities, particularly in understudied regions. Here, we examine how the seasonal activity patterns of ants-key drivers of terrestrial ecosystem functioning-vary with anthropogenic land-cover change on a subtropical island landscape, and whether differences in temperature or species composition can explain observed patterns. Using trap captures sampled biweekly over 2 years from a biodiversity monitoring network covering Okinawa Island, Japan, we processed 1.2 million individuals and reconstructed activity patterns within and across habitat types. Forest communities exhibited greater temporal variability of activity than those in more developed areas. Using time-series decomposition to deconstruct this pattern, we found that sites with greater human development exhibited ant communities with diminished seasonality, reduced synchrony and higher stochasticity compared with sites with greater forest cover. Our results cannot be explained by variation in regional or site temperature patterns, or by differences in species richness or composition among sites. Our study raises the possibility that disruptions to natural seasonal patterns of functionally key insect communities may comprise an important and underappreciated consequence of global environmental change that must be better understood across Earth's biomes.

Influence of metabolic guilds on a temporal scale in an experimental fermented food derived microbial community

The influence of community diversity, which can be measured at the level of metabolic guilds, on community function is a central question in ecology. Particularly, the long-term temporal dynamic between a community's function and its diversity remains unclear. We investigated the influence of metabolic guild diversity on associated community function by propagating natural microbial communities from a traditionally fermented milk beverage diluted to various levels. Specifically, we assessed the influence of less abundant microbial types, such as yeast, on community functionality and bacterial community compositions over repeated propagation cycles amounting to ∼100 generations. The starting richness of metabolic guilds had a repeatable effect on bacterial community compositions, metabolic profiles, and acidity. The influence of a single metabolic guild, yeast in our study, played a dramatic role on function, but interestingly not on long-term species sorting trajectories of the remaining bacterial community. Our results together suggest an unexpected niche division between yeast and bacterial communities and evidence ecological selection on the microbial communities in our system.

Environmental heterogeneity modulates the effect of plant diversity on the spatial variability of grassland biomass

Plant productivity varies due to environmental heterogeneity, and theory suggests that plant diversity can reduce this variation. While there is strong evidence of diversity effects on temporal variability of productivity, whether this mechanism extends to variability across space remains elusive. Here we determine the relationship between plant diversity and spatial variability of productivity in 83 grasslands, and quantify the effect of experimentally increased spatial heterogeneity in environmental conditions on this relationship. We found that communities with higher plant species richness (alpha and gamma diversity) have lower spatial variability of productivity as reduced abundance of some species can be compensated for by increased abundance of other species. In contrast, high species dissimilarity among local communities (beta diversity) is positively associated with spatial variability of productivity, suggesting that changes in species composition can scale up to affect productivity. Experimentally increased spatial environmental heterogeneity weakens the effect of plant alpha and gamma diversity, and reveals that beta diversity can simultaneously decrease and increase spatial variability of productivity. Our findings unveil the generality of the diversity-stability theory across space, and suggest that reduced local diversity and biotic homogenization can affect the spatial reliability of key ecosystem functions.

Pathways to global-change effects on biodiversity: new opportunities for dynamically forecasting demography and species interactions

In structured populations, persistence under environmental change may be particularly threatened when abiotic factors simultaneously negatively affect survival and reproduction of several life cycle stages, as opposed to a single stage. Such effects can then be exacerbated when species interactions generate reciprocal feedbacks between the demographic rates of the different species. Despite the importance of such demographic feedbacks, forecasts that account for them are limited as individual-based data on interacting species are perceived to be essential for such mechanistic forecasting-but are rarely available. Here, we first review the current shortcomings in assessing demographic feedbacks in population and community dynamics. We then present an overview of advances in statistical tools that provide an opportunity to leverage population-level data on abundances of multiple species to infer stage-specific demography. Lastly, we showcase a state-of-the-art Bayesian method to infer and project stage-specific survival and reproduction for several interacting species in a Mediterranean shrub community. This case study shows that climate change threatens populations most strongly by changing the interaction effects of conspecific and heterospecific neighbours on both juvenile and adult survival. Thus, the repurposing of multi-species abundance data for mechanistic forecasting can substantially improve our understanding of emerging threats on biodiversity.

Covariation of taxonomic and functional facets of β-diversity in Chilean freshwater fish assemblages: Implications for current and future processes of biotic homogenization

The biodiversity of assemblages that experience the introduction and extinction of species may lead to responses in two important facets: The taxonomic and functional diversity. The way in which these facets are associated may reveal important implications and consequences for the conservation of those assemblages. Considering the critical situation of freshwater fishes in continental Chile (30° - 56° S), we analyzed how the taxonomic (TDβ) and functional (FDβ) facets of β-diversity, and their components of turnover and nestedness, are associated. We evaluated changes in β-diversity (ΔTDβ and ΔFDβ), turnover (ΔTDtur and ΔFDtur), and nestedness (ΔTDnes and ΔFDnes) in 20 fish assemblages from their historical (pre-European) to current composition. We also simulated future trends of these changes, assuming that native species with conservation issues would become extinct. Our results show that the fish assemblages studied are in a process of loss of β-diversity, both in taxonomic and functional facets (ΔTDβ = -3.9%; ΔFDβ = -30.4%); also, that these facets are positively correlated in the assemblages studied (r = 0.617; P < 0.05). Both components showed by loss in nestedness (ΔTDnes = -36.9%; ΔFDnes = -60.9%) but gain in turnover (ΔTDtur = 9.2%; ΔFDtur = 12.3%). The functional β-diversity decreased more than the taxonomic (ΔFDβ > ΔTDβ), which was caused chiefly by six exotic species of Salmonidae, whose geographical spread was wider and that at the same time shared several morpho-functional traits. Our forecasts, assuming an intensification in the extinction of Endangered and Vulnerable native species, indicate that the process of homogenization will continue, though at a lower rate. Our study shows that the freshwater ichthyofauna of continental Chile is undergoing biotic homogenization, and that this process involves the facets of taxonomic and functional β-diversity, which are show high correlation between historical and current compositions. Both facets show that process is influenced by nestedness, and while turnover contributes to differentiation (both taxonomic and functional), its importance is overshadowed by nestedness.

Individual species and site dynamics are the main drivers of spatial scaling of stability in aquatic communities

Introduction

Any measure of ecological stability scales with the spatial and temporal extent of the data on which it is based. The magnitude of stabilization effects at increasing spatial scale is determined by the degree of synchrony between local and regional species populations.

Methods

We applied two recently developed approaches to quantify these stabilizing effects to time series records from three aquatic monitoring data sets differing in environmental context and organism type.

Results and Discussion

We found that the amount and general patterns of stabilization with increasing spatial scale only varied slightly across the investigated species groups and systems. In all three data sets, the relative contribution of stabilizing effects via asynchronous dynamics across space was higher than compensatory dynamics due to differences in biomass fluctuations across species and populations. When relating the stabilizing effects of individual species and sites to species and site-specific characteristics as well as community composition and aspects of spatial biomass distribution patterns, however, we found that the effects of single species and sites showed large differences and were highly context dependent, i.e., dominant species can but did not necessarily have highly stabilizing or destabilizing effects on overall community biomass. The sign and magnitude of individual contributions depended on community structure and the spatial distribution of biomass and species in space. Our study therefore provides new insights into the mechanistic understanding of ecological stability patterns across scales in natural species communities.

Stability of multi-layer ecosystems

Community structure is reported to play a critical role in ecosystem stability, which indicates the ability of a system to return to equilibrium after perturbations. However, current studies rely on the assumption that the community consists of only a single-layer structure. In practice, multi-layer structures are common in ecosystems, e.g. the distributions of both microorganisms in strata and fish in the ocean usually stratify into multi-layer structures. Here we use multi-layer networks to model species interactions within each layer and between different layers, and study the stability of multi-layer ecosystems under different interaction types. We show that competitive interactions within each layer have a more substantial stabilizing effect in multi-layer ecosystems relative to their impact in single-layer ecosystems. Surprisingly, between different layers, we find that competition between species destabilizes the ecosystem. We further provide a theoretical analysis of the stability of multi-layer ecosystems and confirm the robustness of our findings for different connectances between layers, numbers of species in each layer, and numbers of layers. Our work provides a general framework for investigating the stability of multi-layer ecosystems and uncovers the double-sided role of competitive interactions in the stability of multi-layer ecosystems.

Impacts of extreme climatic events on trophic network complexity and multidimensional stability

Untangling the relationship between network complexity and ecological stability under climate change is an arduous challenge for theoretical and empirical ecology. Even more so, when considering extreme climatic events. Here, we studied the effects of extreme climatic events (heatwaves) on the complexity of realistic freshwater ecosystems using topological and quantitative trophic network metrics. Next, we linked changes in network complexity with the investigation of four stability components (temporal stability, resistance, resilience, and recovery) of community's functional, compositional, and energy flux stability. We found reduction in topological network complexity to be correlated with reduction of functional and compositional resistance. However, temperature-driven increase in link-weighted network complexity increased functional and energy flux recovery and resilience, but at the cost of increased compositional instability. Overall, we propose an overarching approach to elucidate the effects of climate change on multidimensional stability through the lens of network complexity, providing helpful insights for preserving ecosystems stability under climate change.

Temporal Dynamics of Rare and Abundant Soil Bacterial Taxa from Different Fertilization Regimes Under Various Environmental Disturbances

Global climate change has emerged as a critical environmental problem. Different types of climate extremes drive soil microbial communities to alternative states, leading to a series of consequences for soil microbial ecosystems and related functions. An effective method is urgently needed for buffering microbial communities to tackle environmental disturbances. Here, we conducted a series of mesocosm experiments in which the organic (NOF) and chemical fertilizer (NCF) long-term-amended soil microbiotas were subjected to environmental disturbances that included drought, flooding, freeze-thaw cycles, and heat. We subsequently tracked the temporal dynamics of rare and abundant bacterial taxa in NOF and NCF treatment soils to assess the efficiencies of organic amendments in recovery of soil microbiome. Our results revealed that freeze-thaw cycles and drought treatments showed weaker effects on bacterial communities than flooding and heat. The turnover between rare and abundant taxa occurred in postdisturbance succession of flooding and heat treatments, indicating that new equilibria were tightly related to the rare taxa in both NCF and NOF treatment soils. The Bayesian fits of modeling for the microbiome recovery process revealed that the stability of abundant taxa in NOF was higher than that in NCF soil. In particular, the NOF treatment soil reduced the divergence from the initial bacterial community after weak perturbations occurred. Together, we demonstrated that long-term organic input is an effective strategy to enhance the thresholds for transition to alternative states via enhancing the stability of abundant bacterial species. These findings provide a basis for the sustainable development of agricultural ecosystems in response to changing climates. IMPORTANCE Different climate extremes are expected to be a major threat to crop production, and the soil microbiome has been known to play a crucial role in agricultural ecosystems. In recent years, we have known that organic amendments are an effective method for optimizing the composition and functioning of the soil microbial community and maintaining the health of the soil ecosystem. However, the effects of organic fertilization on buffering bacterial communities against environmental disturbances and the underlying mechanisms are still unclear. We conducted a series of mesocosm experiments and showed that organic fertilizers had additional capacities in promoting the soil microbiome to withstand climate change effects. Our study provides both mechanistic insights as well as a direct guide for the sustainable development of agricultural ecosystems in response to climate change.

Organic fertilization enhances the resistance and resilience of soil microbial communities under extreme drought

INTRODUCTION

The soil bacterial microbiome plays a crucial role in ecosystem functioning. The composition and functioning of the microbiome are tightly controlled by the physicochemical surrounding. Therefore, the microbiome is responsive to management, such as fertilization, and to climate change, such as extreme drought. It remains a challenge to retain microbiome functioning under drought.

OBJECTIVES

This work aims to reveal if fertilization with organic fertilizer, can enhance resistance and resilience of bacterial communities and their function in extreme drought and subsequent rewetting compared with conventional fertilizers.

METHODS

In soil mesocosms, we induced a long-term drought for 80 days with subsequent rewetting for 170 days to follow bacterial community dynamics in organic (NOF) and chemical (NCF) fertilization regimes.

RESULTS

Our results showed that bacterial diversity was higher with NOF than with NCF during drought. In particular, the ecological resilience and recovery of bacterial communities under NOF were higher than in NCF. We found these bacterial community features to enhance pathogen-inhibiting functions in NOF compared to NCF during late recovery. The other soil ecology functional analyses revealed that bacterial biomass recovered in the early stage after rewetting, while soil respiration increased continuously following prolonged time after rewetting.

CONCLUSION

Together, our study indicates that organic fertilization can enhance the stability of the soil microbiome and ensures that specific bacterial-driven ecosystem functions recover after rewetting. This may provide the basis for more sustainable agricultural practices to counterbalance negative climate change-induced effects on soil functioning.

Anthropogenic pressures enhance the deleterious effects of extreme storms on rocky shore communities

Climate change is not only changing the mean values of environmental parameters that modulate ecosystems, but also the regime of disturbances. Among them, extreme events have a key role in structuring biological communities. Ecosystems are frequently suffering multiple anthropogenic pressures which can cause effects that are not additive. Thus, the effects of extreme events need to be studied in combination with other pressures to adequately evaluate their consequences. We performed a manipulative approach in two rocky shores in the Mediterranean with contrasting levels of anthropogenic pressure (mainly eutrophication) simulating storms with different disturbance regimes in the intertidal and subtidal zones. In the short-term, an extreme storm had a greater impact on the species assemblage than other disturbance regimes, being especially notable in the area suffering from a high anthropogenic pressure. In this area, the species assemblages that suffered from an extreme storm took a longer time to recover than the ones affected by other disturbance regimes and were generally more affected after the disturbance. The intertidal zone, having more variable environmental conditions than the subtidal zone, was more resistant and able to recover from extreme storms. Our results suggest that the effects of extreme events on biological communities could be strengthened when co-occurring with anthropogenic pressures, especially ecosystems adapted to less variable environmental conditions. Thus, limiting other anthropogenic pressures that ecosystems are suffering is crucial to maintain the natural resistance and recovery capacity of ecosystems towards extreme events such as storms.

Nutrients and herbivores impact grassland stability across spatial scales through different pathways

Nutrients and herbivores are well-known drivers of grassland diversity and stability in local communities. However, whether they interact to impact the stability of aboveground biomass and whether these effects depend on spatial scales remain unknown. It is also unclear whether nutrients and herbivores impact stability via different facets of plant diversity including species richness, evenness, and changes in community composition through time and space. We used a replicated experiment adding nutrients and excluding herbivores for 5 years in 34 global grasslands to explore these questions. We found that both nutrient addition and herbivore exclusion alone reduced stability at the larger spatial scale (aggregated local communities; gamma stability), but through different pathways. Nutrient addition reduced gamma stability primarily by increasing changes in local community composition over time, which was mainly driven by species replacement. Herbivore exclusion reduced gamma stability primarily by decreasing asynchronous dynamics among local communities (spatial asynchrony). Their interaction weakly increased gamma stability by increasing spatial asynchrony. Our findings indicate that disentangling the processes operating at different spatial scales may improve conservation and management aiming at maintaining the ability of ecosystems to reliably provide functions and services for humanity.

Predators mitigate the destabilising effects of heatwaves on multitrophic stream communities

Amidst the global extinction crisis, climate change will expose ecosystems to more frequent and intense extreme climatic events, such as heatwaves. Yet, whether predator species loss-a prevailing characteristic of the extinction crisis-will exacerbate the ecological consequences of extreme climatic events remains largely unknown. Here, we show that the loss of predator species can interact with heatwaves to moderate the compositional stability of ecosystems. We exposed multitrophic stream communities, with and without a dominant predator species, to realistic current and future heatwaves and found that heatwaves destabilised algal communities by homogenising them in space. However, this happened only when the predator was absent. Additional heatwave impacts on multiple aspects of stream communities, including changes to the structure of algal and macroinvertebrate communities, as well as total algal biomass and its temporal variability, were not apparent during heatwaves and emerged only after the heatwaves had passed. Taken together, our results suggest that the ecological consequences of heatwaves can amplify over time as their impacts propagate through biological interaction networks, but the presence of predators can help to buffer such impacts. These findings underscore the importance of conserving trophic structure, and highlight the potential for species extinctions to amplify the effects of climate change and extreme events.

Facing up to climate change: Community composition varies with aspect and surface temperature in the rocky intertidal

Marine rocky intertidal organisms are amongst those most affected by climate change with regional distributional changes observed for many species. Although often ascribed to increased sea surface temperatures, precise assessment of the local habitat conditions underpinning observed and predicted changes in community assembly is lacking. Here we examine how aspect (i.e. north-south orientation) affects intertidal community composition and how rock surface temperatures and stress responses of two dominant grazer species (Patella spp.) elucidate emergent differences. We quantified year-round temperature variation and surveyed intertidal community composition on paired natural rock gullies with Equator- (EF) and Pole-facing (PF) surfaces. We also investigated variation in limpet (Patella spp.) reproductive phenology and osmotic stress. Average annual temperatures were 0.8 °C (1.6 °C at low tide) higher, with six-fold more frequent extremes (i.e. > 30 °C) on EF than PF surfaces. Intertidal community composition varied with aspect across trophic levels with greater overall species richness, abundance of primary producers and grazers on PF-surfaces, and greater barnacle abundance on EF-surfaces. Although species richness of organisms from different biogeographical origins ('Boreal' or 'Lusitanian') did not vary, the Lusitanian limpet Patella depressa exhibited earlier reproductive development on EF-surfaces and both limpet species exhibited greater thermal stress on EF-surfaces. We argue that our study system provides a good model for understanding how temperature variation at local scales can affect community composition, as well as ecophysiological and ecological responses to climate change and so better inform and predict regional range shifts over coming decades.

Global patterns of resilience decline in vertebrate populations

Maintaining the resilience of natural populations, their ability to resist and recover from disturbance, is crucial to prevent biodiversity loss. However, the lack of appropriate data and quantitative tools has hampered our understanding of the factors determining resilience on a global scale. Here, we quantified the temporal trends of two key components of resilience-resistance and recovery-in >2000 population time-series of >1000 vertebrate species globally. We show that the number of threats to which a population is exposed is the main driver of resilience decline in vertebrate populations. Such declines are driven by a non-uniform loss of different components of resilience (i.e. resistance and recovery). Increased anthropogenic threats accelerating resilience loss through a decline in the recovery ability-but not resistance-of vertebrate populations. These findings suggest we may be underestimating the impacts of global change, highlighting the need to account for the multiple components of resilience in global biodiversity assessments.

Universal scaling of robustness of ecosystem services to species loss

Ensuring reliable supply of services from nature is key to the sustainable development and well-being of human societies. Varied and frequently complex relationships between biodiversity and ecosystem services have, however, frustrated our capacity to quantify and predict the vulnerability of those services to species extinctions. Here, we use a qualitative Boolean modelling framework to identify universal drivers of the robustness of ecosystem service supply to species loss. These drivers comprise simple features of the networks that link species to the functions they perform that, in turn, underpin a service. Together, they define what we call network fragility. Using data from >250 real ecological networks representing services such as pollination and seed-dispersal, we demonstrate that network fragility predicts remarkably well the robustness of empirical ecosystem services. We then show how to quantify contributions of individual species to ecosystem service robustness, enabling quantification of how vulnerability scales from species to services. Our findings provide general insights into the way species, functional traits, and the links between them together determine the vulnerability of ecosystem service supply to biodiversity loss.

Species diversity and food web structure jointly shape natural biological control in agricultural landscapes

Land-use change and agricultural intensification concurrently impact natural enemy (e.g., parasitoid) communities and their associated ecosystem services (ESs), i.e., biological pest control. However, the extent to which (on-farm) parasitoid diversity and food webs mediate landscape-level influences on biological control remains poorly understood. Here, drawing upon a 3-year study of quantitative parasitoid-hyperparasitoid trophic networks from 25 different agro-landscapes, we assess the cascading effects of landscape composition, species diversity and trophic network structure on ecosystem functionality (i.e., parasitism, hyperparasitism). Path analysis further reveals cascaded effects leading to biological control of a resident crop pest, i.e., Aphis gossypii. Functionality is dictated by (hyper)parasitoid diversity, with its effects modulated by food web generality and vulnerability. Non-crop habitat cover directly benefits biological control, whereas secondary crop cover indirectly lowers hyperparasitism. Our work underscores a need to simultaneously account for on-farm biodiversity and trophic interactions when investigating ESs within dynamic agro-landscapes.

Extreme flood disturbance effects on multiple dimensions of river invertebrate community stability

Multidimensional analysis of community stability has recently emerged as an overarching approach to evaluating ecosystem response to disturbance. However, the approach has previously been applied only in experimental and modelling studies. We applied this concept to an 18-year time series (2000-2017) of macroinvertebrate community dynamics from a southeast Alaskan river to further develop and test the approach in relation to the effects of two extreme flood events occurring in 2005 (event 1) and 2014 (event 2). Five components of stability were calculated for pairs of pre- or post-event years. Individual components were tested for differences between pre- and post-event time periods. Stability components' pairwise correlations were assessed and ellipsoids of stability were developed for each time period and compared to a null model derived from the permuted dataset. Only one stability component demonstrated a significant difference between time periods. In contrast, 80% of moderate and significant correlations between stability components were degraded post-disturbance and significant changes to the form of stability ellipsoids were observed. Ellipsoids of stability for all periods after the initial disturbance (2005) were not different to the null model. Our results illustrate that the dimensionality of stability approach can be applied to natural ecosystem time-series data. The major increase in dimensionality of stability observed following disturbance potentially indicates significant shifts in the processes which drive stability following disturbance. This evidence improves our understanding of community response beyond what is possible through analysis of individual stability components.

Effects of multiple stressors on the dimensionality of ecological stability

Ecological stability is a multidimensional construct. Investigating multiple stability dimensions is key to understand how ecosystems respond to disturbance. Here, we evaluated the single and combined effects of common agricultural stressors (insecticide, herbicide and nutrients) on four dimensions of stability (resistance, resilience, recovery and invariability) and on the overall dimensionality of stability (DS) using the results of a freshwater mesocosm experiment. Functional recovery and resilience to pesticides were enhanced in nutrient-enriched systems, whereas compositional recovery was generally not achieved. Pesticides did not affect compositional DS, whereas functional DS was significantly increased by the insecticide only in non-enriched systems. Stressor interactions acted non-additively on single stability dimensions as well as on functional DS. Moreover we demonstrate that pesticides can modify the correlation between functional and compositional aspects of stability. Our study shows that different disturbance types, and their interactions, require specific management actions to promote ecosystem stability.