Interactive effects of multiple stressors vary with consumer interactions, stressor dynamics and magnitude

Fungal identity mediates the impacts of multiple stressors on freshwater ecosystems

Predicting how multiple anthropogenic stressors affect natural ecosystems is a major challenge in ecology. Freshwater ecosystems are threatened worldwide by multiple co-occurring stressors, which can affect aquatic biodiversity, ecosystem functioning and human wellbeing. In stream ecosystems, aquatic fungi play a crucial role in global biogeochemical cycles and food web dynamics, therefore, assessing the functional consequences of fungal biodiversity loss under multiple stressors is crucial. Here, a microcosm approach was used to investigate the effects of multiple stressors (increased temperature and nutrients, drying, and biodiversity loss) on three ecosystem processes: organic matter decomposition, fungal reproduction, and fungal biomass accrual. Net effects of stressors were antagonistic for organic matter decomposition, but additive for fungal reproduction and biomass accrual. Net effects of biodiversity were mainly positive for all processes, even under stress, demonstrating that diversity assures the maintenance of ecosystem processes. Fungal species displayed distinct contributions to each ecosystem process. Furthermore, species with negligible contributions under control conditions changed their role under stress, either enhancing or impairing the communities' performance, emphasizing the importance of fungal species identity. Our study highlights that distinct fungal species have different sensitivities to environmental variability and have different influence on the overall performance of the community. Therefore, preserving high fungal diversity is crucial to maintain fungal species with key ecosystem functions within aquatic communities in face of environmental change.

High-throughput characterization of bacterial responses to complex mixtures of chemical pollutants

Our understanding of how microbes respond to micropollutants, such as pesticides, is almost wholly based on single-species responses to individual chemicals. However, in natural environments, microbes experience multiple pollutants simultaneously. Here we perform a matrix of multi-stressor experiments by assaying the growth of model and non-model strains of bacteria in all 255 combinations of 8 chemical stressors (antibiotics, herbicides, fungicides and pesticides). We found that bacterial strains responded in different ways to stressor mixtures, which could not be predicted simply from their phylogenetic relatedness. Increasingly complex chemical mixtures were both more likely to negatively impact bacterial growth in monoculture and more likely to reveal net interactive effects. A mixed co-culture of strains proved more resilient to increasingly complex mixtures and revealed fewer interactions in the growth response. These results show predictability in microbial population responses to chemical stressors and could increase the utility of next-generation eco-toxicological assays.

Integrating climate model projections into environmental risk assessment: A probabilistic modeling approach

The Society of Environmental Toxicology and Chemistry (SETAC) convened a Pellston workshop in 2022 to examine how information on climate change could be better incorporated into the ecological risk assessment (ERA) process for chemicals as well as other environmental stressors. A major impetus for this workshop is that climate change can affect components of ecological risks in multiple direct and indirect ways, including the use patterns and environmental exposure pathways of chemical stressors such as pesticides, the toxicity of chemicals in receiving environments, and the vulnerability of species of concern related to habitat quality and use. This article explores a modeling approach for integrating climate model projections into the assessment of near- and long-term ecological risks, developed in collaboration with climate scientists. State-of-the-art global climate modeling and downscaling techniques may enable climate projections at scales appropriate for the study area. It is, however, also important to realize the limitations of individual global climate models and make use of climate model ensembles represented by statistical properties. Here, we present a probabilistic modeling approach aiming to combine projected climatic variables as well as the associated uncertainties from climate model ensembles in conjunction with ERA pathways. We draw upon three examples of ERA that utilized Bayesian networks for this purpose and that also represent methodological advancements for better prediction of future risks to ecosystems. We envision that the modeling approach developed from this international collaboration will contribute to better assessment and management of risks from chemical stressors in a changing climate. Integr Environ Assess Manag 2024;20:367-383. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

A synthesis of anthropogenic stress effects on emergence-mediated aquatic-terrestrial linkages and riparian food webs

Anthropogenic stress alters the linkage between aquatic and terrestrial ecosystems in various ways. Here, we review the contemporary literature on how alterations in aquatic systems through environmental pollution, invasive species and hydromorphological changes carry-over to terrestrial ecosystems and the food webs therein. We consider both the aquatic insect emergence and flooding as pathways through which stressors can propagate from the aquatic to the terrestrial system. We specifically synthesize and contextualize results on the roles of pollutants in the emergence pathway and their top-down consequences. Our review revealed that the emergence and flooding pathway are only considered in isolation and that the overall effects of invasive species or pollutants on food webs at the water-land interface require further attention. While very few recent studies looked at invasive species, a larger number of studies focused on metal transfer compared to pesticides, pharmaceuticals or PCBs, and multiple stress studies up to now left aquatic-terrestrial linkages unconsidered. Recent research on pollutants and emergence used aquatic-terrestrial mesocosms to elucidate the effects of aquatic stressors such as the mosquito control agent Bti, metals or pesticides to understand the effects on riparian spiders. Quality parameters, such as the structural and functional composition of emergent insect communities, the fatty acid profiles, yet also the composition of pollutants transferred to land prove to be important for the effects on riparian spiders. Process-based models including quality of emergence are useful to predict the resulting top-down directed food web effects in the terrestrial recipient ecosystem. In conclusion, we present and recommend a combination of empirical and modelling approaches in order to understand the complexity of aquatic-terrestrial stressor propagation and its spatial and temporal variation.

A class of their own? Water-soluble polymer pollution impacting a freshwater host-pathogen system

While the inclusion of synthetic polymers such as primary microplastics within personal care products have been widely restricted under EU/UK Law, water-soluble polymers (WSPs) have so far slipped the net of global chemical regulation despite evidence that these could be polluting wastewater effluents at concentrations greatly exceeding those of microplastics. Polyvinyl alcohol (PVA) and polyvinylpyrrolidone (PVP) represent WSPs with common industry and household uses, down-the-drain disposal and a direct route to wastewater treatment plants, conveying high risk of environmental leaching into freshwater ecosystems. The current study is the first investigating the impacts of predicted environmental concentrations of these WSPs on life-history traits of two freshwater species also constituting a disease model (fish - Poecilia reticulata and parasite - Gyrodactylus turnbulli). Single effects of WSPs on fish as well as their interactive effects with infection of the ectoparasite were determined over a 45-day exposure. Generally, WSPs reduced fish growth and increased routine metabolic rate of fish implying a depleted energetic budget, however these effects were dose, exposure time and polymer dependent. Parasitic infection alone caused a significant reduction in fish growth and enhanced fish routine metabolic rate. In contrast, a non-additive effect on metabolic rate was evident in fish experiencing simultaneous infection and WSP exposure, suggesting a protective effect of the two WSPs for fish also exposed to a metazoan ectoparasite. Off-host parasite survival was significantly lowered by both WSPs; however, parasite counts of infected fish also exposed to WSP were not significantly different from the control, implying more complex mechanisms may underpin this stressor interaction. Distinct detrimental impacts were inflicted on both organisms implying environmental leaching of WSPs may be causing significant disruption to interspecies interactions within freshwater ecosystems. Additionally, these results could contribute to sustainable development in industry, as we conclude PVA represents a less harmful alternative to PVP.

A process framework for integrating stressor-response functions into cumulative effects models

Stressor-response (SR) functions quantify ecological responses to natural environmental variation or anthropogenic stressors. They are also core drivers of cumulative effects (CE) models, which are increasingly recognized as essential management tools to grapple with the diffuse footprint of human impacts. Here, we provide a process framework for the identification, development, and integration of SR functions into CE models, and highlight their consequential properties, behaviour, criteria for selecting appropriate stressors and responses, and general approaches for deriving them. Management objectives (and causal effect pathways) will determine the ultimate stressor and target response variables of interest (i.e., individual growth/survival, population size, community structure, ecosystem processes), but data availability will constrain whether proxies need to be used for the target stressor or response variables. Available data and confidence in underlying mechanisms will determine whether empirical or mechanistic (theoretical) SR functions are optimal. Uncertainty in underlying SR functions is often the primary source of error in CE modelling, and monitoring outcomes through adaptive management to iteratively refine parameterization of SR functions is a key element of model application. Dealing with stressor interactions is an additional challenge, and in the absence of known or suspected interaction mechanisms, controlling main effects should remain the primary focus. Indicators of suspected interaction presence (i.e., much larger or smaller responses to stressor reduction than expected during monitoring) should be confirmed through adaptive management cycles or targeted stressor manipulations. Where possible, management decisions should selectively take advantage of interactions to strategically mitigate stressor impacts (i.e., by using antagonisms to suppress stressor impacts, and by using synergisms to efficiently reduce them).

Cumulative effects of multiple stressors impact an endangered seagrass population and fish communities

Coastal ecosystems are becoming increasingly threatened by human activities and there is growing appreciation that management must consider the impacts of multiple stressors. Cumulative effects assessments (CEAs) have become a popular tool for identifying the distribution and intensity of multiple human stressors in coastal ecosystems. Few studies, however, have demonstrated strong correlations between CEAs and change in ecosystem condition, questioning its management use. Here, we apply a CEA to the endangered seagrass Posidonia australis in Pittwater, NSW, Australia, using spatial data on known stressors to seagrass related to foreshore development, water quality, vessel traffic and fishing. We tested how well cumulative effects scores explained changes in P. australis extent measured between 2005 and 2019 using high-resolution aerial imagery. A negative correlation between P. australis and estimated cumulative effects scores was observed (R2 = 22 %), and we identified a threshold of cumulative effects above which losses of P. australis became more likely. Using baited remote underwater video, we surveyed fishes over P. australis and non-vegetated sediments to infer and quantify how impacts of cumulative effects to P. australis extent would flow on to fish assemblages. P. australis contained a distinct assemblage of fish, and on non-vegetated sediments the abundance of sparids, which are of importance to fisheries, increased with closer proximity to P. australis. Our results demonstrate the negative impact of multiple stressors on P. australis and the consequences for fish biodiversity and fisheries production across much of the estuary. Management actions aimed at reducing or limiting cumulative effects to low and moderate levels will help conserve P. australis and its associated fish biodiversity and productivity.

Ecosystem-level responses to multiple stressors using a time-dynamic food-web model: The case of a re-oligotrophicated coastal embayment (Saronikos Gulf, E Mediterranean)

Multiple stressors may combine in unexpected ways to alter the structure of ecological systems, however, our current ability to evaluate their ecological impact is limited due to the lack of information concerning historic trophic interactions and ecosystem dynamics. Saronikos Gulf is a heavily exploited embayment in the E Mediterranean that has undergone significant ecological alterations during the last 20 years including a shift from long-standing eutrophic to oligotrophic conditions in the mid-2000's. Here we used a historical Ecopath food-web model of Saronikos Gulf (1998-2000) and fitted the time-dynamic module Ecosim to biomass and catch time series for the period 2001-2020. We then projected the model forward in time from 2021 to 2050 under 8 scenarios to simulate ecosystem responses to the individual and combined effect of sea surface temperature increase, primary productivity shifts and fishing effort release. Incorporating trophic interactions, climate warming, fishing and primary production improved model fit, depicting that both fishing and the environment have historically influenced ecosystem dynamics. Retrospective simulations of the model captured historical biomass and catch trends of commercially important stocks and reproduced successfully the marked recovery of marine resources 10 years after re-oligotrophication. In future scenarios increasing temperature had a detrimental impact on most functional groups, increasing and decreasing productivity had a positive and negative effect on all respectively, while fishing reductions principally benefited top predators. Combined stressors produced synergistic or antagonistic effects depending on the direction and magnitude of change of each stressor in isolation while their overall impact seemed to be strongly mediated via food-web interactions. Such holistic approaches advance of our mechanistic understanding of ecosystems enabling us to develop more effective management strategies in the face of a rapidly changing marine environment.

Winter is coming: Interactions of multiple stressors in winter and implications for the natural world

Winter is a key driver of ecological processes in freshwater, marine and terrestrial ecosystems, particularly in higher latitudes. Species have evolved various adaptive strategies to cope with food limitations and the cold and dark wintertime. However, human-induced climate change and other anthropogenic stressors are impacting organisms in winter in unpredictable ways. In this paper, we show that global change experiments investigating multiple stressors have predominantly been conducted during summer months. However, effects of anthropogenic stressors sometimes differ between winter and other seasons, necessitating comprehensive investigations. Here, we outline a framework for understanding the different effects of anthropogenic stressors in winter compared to other seasons and discuss the primary mechanisms that will alter ecological responses of organisms (microbes, animals and plants). For instance, while the magnitude of some anthropogenic stressors can be greater in winter than in other seasons (e.g. some pollutants), others may alleviate natural winter stress (e.g. warmer temperatures). These changes can have immediate, delayed or carry-over effects on organisms during winter or later seasons. Interactions between stressors may also vary with season. We call for a renewed research direction focusing on multiple stressor effects on winter ecology and evolution to fully understand, and predict, how ecosystems will fare under changing winters. We also argue the importance of incorporating the interactions of anthropogenic stressors with winter into ecological risk assessments, management and conservation efforts.

Diversity loss from multiple interacting disturbances is regime-dependent

Anthropogenic activities expose many ecosystems to multiple novel disturbances simultaneously. Despite this, how biodiversity responds to simultaneous disturbances remains unclear, with conflicting empirical results on their interactive effects. Here, we experimentally test how one disturbance (an invasive species) affects the diversity of a community over multiple levels of another disturbance regime (pulse mortality). Specifically, we invade stably coexisting bacterial communities under four different pulse frequencies, and compare their final resident diversity to uninvaded communities under the same pulse mortality regimes. Our experiment shows that the disturbances synergistically interact, such that the invader significantly reduces resident diversity at high pulse frequency, but not at low. This work therefore highlights the need to study simultaneous disturbance effects over multiple disturbance regimes as well as to carefully document unmanipulated disturbances, and may help explain the conflicting results seen in previous multiple-disturbance work.

Mean species responses predict effects of environmental change on coexistence

Environmental change research is plagued by the curse of dimensionality: the number of communities at risk and the number of environmental drivers are both large. This raises the pressing question if a general understanding of ecological effects is achievable. Here, we show evidence that this is indeed possible. Using theoretical and simulation-based evidence for bi- and tritrophic communities, we show that environmental change effects on coexistence are proportional to mean species responses and depend on how trophic levels on average interact prior to environmental change. We then benchmark our findings using relevant cases of environmental change, showing that means of temperature optima and of species sensitivities to pollution predict concomitant effects on coexistence. Finally, we demonstrate how to apply our theory to the analysis of field data, finding support for effects of land use change on coexistence in natural invertebrate communities.

Chemical Mixtures and Multiple Stressors: Same but Different?

Ecosystems are strongly influenced by multiple anthropogenic stressors, including a wide range of chemicals and their mixtures. Studies on the effects of multiple stressors have largely focussed on nonchemical stressors, whereas studies on chemical mixtures have largely ignored other stressors. However, both research areas face similar challenges and require similar tools and methods to predict the joint effects of chemicals or nonchemical stressors, and frameworks to integrate multiple chemical and nonchemical stressors are missing. We provide an overview of the research paradigms, tools, and methods commonly used in multiple stressor and chemical mixture research and discuss potential domains of cross-fertilization and joint challenges. First, we compare the general paradigms of ecotoxicology and (applied) ecology to explain the historical divide. Subsequently, we compare methods and approaches for the identification of interactions, stressor characterization, and designing experiments. We suggest that both multiple stressor and chemical mixture research are too focused on interactions and would benefit from integration regarding null model selection. Stressor characterization is typically more costly for chemical mixtures. While for chemical mixtures comprehensive classification systems at suborganismal level have been developed, recent classification systems for multiple stressors account for environmental context. Both research areas suffer from rather simplified experimental designs that focus on only a limited number of stressors, chemicals, and treatments. We discuss concepts that can guide more realistic designs capturing spatiotemporal stressor dynamics. We suggest that process-based and data-driven models are particularly promising to tackle the challenge of prediction of effects of chemical mixtures and nonchemical stressors on (meta-)communities and (meta-)food webs. We propose a framework to integrate the assessment of effects for multiple stressors and chemical mixtures. Environ Toxicol Chem 2023;42:1915-1936. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Adverse effects of the pesticide chlorpyrifos on the physiology of a damselfly only occur at the cold and hot extremes of a temperature gradient

Ecotoxicological studies considerably improved realism by assessing the toxicity of pollutants at different temperatures. Nevertheless, they may miss key interaction patterns between pollutants and temperature by typically considering only part of the natural thermal gradient experienced by species and ignoring daily temperature fluctuations (DTF). We therefore tested in a common garden laboratory experiment the effects of the pesticide chlorpyrifos across a range of mean temperatures and DTF on physiological traits (related to oxidative stress and bioenergetics) in low- and high-latitude populations of Ischnura elegans damselfly larvae. As expected, the impact of chlorpyrifos varied along the wide range of mean temperatures (12-34 °C). None of the physiological traits (except the superoxide anion levels) were affected by chlorpyrifos at the intermediate mean temperatures (20-24 °C). Instead, most of them were negatively affected by chlorpyrifos (reduced activity levels of the antioxidant defense enzymes superoxide dismutase [SOD], catalase [CAT] and peroxidase [PER], and a reduced energy budget) at the very high (≥28 °C) or extreme high temperatures (≥32 °C), and to lesser extent at the lower mean temperatures (≤16 °C). Notably, at the lower mean temperatures the negative impact of chlorpyrifos was often only present or stronger under DTF. Although the chlorpyrifos effects on the physiological traits greatly depended on the experimentally imposed thermal gradient, patterns were mainly consistent across the natural latitude-associated thermal gradient, indicating the generality of our results. The thermal patterns in chlorpyrifos-induced physiological responses contributed to the observed toxicity patterns in life history (reduced survival and growth at low and high mean temperatures). Taken together, our results underscore the importance of evaluating pesticide toxicity along a temperature gradient and of taking a mechanistic approach with a focus on physiology, to improve our understanding of the combined effects of pollutants and temperature in natural populations.

The Asymmetric Response Concept explains ecological consequences of multiple stressor exposure and release

Our capacity to predict trajectories of ecosystem degradation and recovery is limited, especially when impairments are caused by multiple stressors. Recovery may be fast or slow and either complete or partial, sometimes result in novel ecosystem states or even fail completely. Here, we introduce the Asymmetric Response Concept (ARC) that provides a basis for exploring and predicting the pace and magnitude of ecological responses to, and release from, multiple stressors. The ARC holds that three key mechanisms govern population, community and ecosystem trajectories. Stress tolerance is the main mechanism determining responses to increasing stressor intensity, whereas dispersal and biotic interactions predominantly govern responses to the release from stressors. The shifting importance of these mechanisms creates asymmetries between the ecological trajectories that follow increasing and decreasing stressor intensities. This recognition helps to understand multiple stressor impacts and to predict which measures will restore communities that are resistant to restoration.

Tipping points and interactive effects of chronic human disturbance and acute heat stress on coral diversity

Multiple anthropogenic stressors co-occur ubiquitously in natural ecosystems. However, multiple stressor studies often produce conflicting results, potentially because the nature and direction of stressor interactions depends upon the strength of the underlying stressors. Here, we first examine how coral α- and β-diversities vary across sites spanning a gradient of chronic local anthropogenic stress before and after a prolonged marine heatwave. Developing a multiple stressor framework that encompasses non-discrete stressors, we then examine interactions between the continuous and discrete stressors. We provide evidence of additive effects, antagonistic interactions (with heatwave-driven turnover in coral community composition diminishing as the continuous stressor increased), and tipping points (at which the response of coral Hill-richness to stressors changed from additive to near synergistic). We show that community-level responses to multiple stressors can vary, and even change qualitatively, with stressor intensity, underscoring the importance of examining complex, but realistic continuous stressors to understand stressor interactions and their ecological impacts.

Predicting effects of multiple interacting global change drivers across trophic levels

Global change encompasses many co-occurring anthropogenic drivers, which can act synergistically or antagonistically on ecological systems. Predicting how different global change drivers simultaneously contribute to observed biodiversity change is a key challenge for ecology and conservation. However, we lack the mechanistic understanding of how multiple global change drivers influence the vital rates of multiple interacting species. We propose that reaction norms, the relationships between a driver and vital rates like growth, mortality, and consumption, provide insights to the underlying mechanisms of community responses to multiple drivers. Understanding how multiple drivers interact to affect demographic rates using a reaction-norm perspective can improve our ability to make predictions of interactions at higher levels of organization-that is, community and food web. Building on the framework of consumer-resource interactions and widely studied thermal performance curves, we illustrate how joint driver impacts can be scaled up from the population to the community level. A simple proof-of-concept model demonstrates how reaction norms of vital rates predict the prevalence of driver interactions at the community level. A literature search suggests that our proposed approach is not yet used in multiple driver research. We outline how realistic response surfaces (i.e., multidimensional reaction norms) can be inferred by parametric and nonparametric approaches. Response surfaces have the potential to strengthen our understanding of how multiple drivers affect communities as well as improve our ability to predict when interactive effects emerge, two of the major challenges of ecology today.

Thermal Performance Curves in a Polluted World: Too Cold and Too Hot Temperatures Synergistically Increase Pesticide Toxicity

Ecotoxicological studies typically cover only a limited part of the natural thermal range of populations and ignore daily temperature fluctuations (DTFs). Therefore, we may miss important stressor interaction patterns and have poor knowledge on how pollutants affect thermal performance curves (TPCs), which is needed to improve insights into the fate of populations to warming in a polluted world. We tested the single and combined effects of pesticide exposure and DTFs on the TPCs of low- and high-latitude populations of Ischnura elegans damselfly larvae. While chlorpyrifos did not have any effect at the intermediate mean temperatures (20-24 °C), it became toxic (reflecting synergisms) at lower (≤16 °C, reduced growth) and especially at higher (≥28 °C, reduced survival and growth) mean temperatures, resulting in more concave-shaped TPCs. Remarkably, these toxicity patterns were largely consistent at both latitudes and hence across a natural thermal gradient. Moreover, DTFs magnified the pesticide-induced survival reductions at 34 °C. The TPC perspective allowed us to identify different toxicity patterns and interaction types (mainly additive vs synergistic) across the thermal gradient. This highlights the importance of using thermal gradients to make more realistic predictions about the impact of pesticides in a warming world and of warming in a polluted world.

Mercury Contamination Challenges the Behavioral Response of a Keystone Species to Arctic Climate Change

Combined effects of multiple, climate change-associated stressors are of mounting concern, especially in Arctic ecosystems. Elevated mercury (Hg) exposure in Arctic animals could affect behavioral responses to changes in foraging landscapes caused by climate change, generating interactive effects on behavior and population resilience. We investigated this hypothesis in little auks (Alle alle), a keystone Arctic seabird. We compiled behavioral data for 44 birds across 5 years using accelerometers while also quantifying blood Hg and environmental conditions. Warm sea surface temperature (SST) and low sea ice coverage reshaped time activity budgets (TABs) and diving patterns, causing decreased resting, increased flight, and longer dives. Mercury contamination was not associated with TABs. However, highly contaminated birds lengthened interdive breaks when making long dives, suggesting Hg-induced physiological limitations. As dive durations increased with warm SST, subtle toxicological effects threaten to increasingly constrain diving and foraging efficiency as climate change progresses, with ecosystem-wide repercussions.

Perceived multiple stressor effects depend on sample size and stressor gradient length

Multiple stressors are continuously deteriorating surface waters worldwide, posing many challenges for their conservation and restoration. Combined effect types of multiple stressors range from single-stressor dominance to complex interactions. Identifying prevalent combined effect types is critical for environmental management, as it helps to prioritise key stressors for mitigation. However, it remains unclear whether observed single and combined stressor effects reflect true ecological processes unbiased by sample size and length of stressor gradients. Therefore, we examined the role of sample size and stressor gradient lengths in 158 paired-stressor response cases with over 120,000 samples from rivers, lakes, transitional and marine ecosystems around the world. For each case, we split the overall stressor gradient into two partial gradients (lower and upper) and investigated associated changes in single and combined stressor effects. Sample size influenced the identified combined effect types, and stressor interactions were less likely for cases with fewer samples. After splitting gradients, 40 % of cases showed a change in combined effect type, 30 % no change, and 31 % showed a loss in stressor effects. These findings suggest that identified combined effect types may often be statistical artefacts rather than representing ecological processes. In 58 % of cases, we observed changes in stressor effect directions after the gradient split, suggesting unimodal stressor effects. In general, such non-linear responses were more pronounced for organisms at higher trophic levels. We conclude that observed multiple stressor effects are not solely determined by ecological processes, but also strongly depend on sampling design. Observed effects are likely to change when sample size and/or gradient length are modified. Our study highlights the need for improved monitoring programmes with sufficient sample size and stressor gradient coverage. Our findings emphasize the importance of adaptive management, as stress reduction measures or further ecosystem degradation may change multiple stressor-effect relationships, which will then require associated changes in management strategies.