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.

The impact of cumulative stressor effects on uncertainty and ecological risk.

To enable environmental management actions to be more effectively prioritized, cumulative effects between multiple stressors need to be accounted for in risk-assessment frameworks. Ecological risk and uncertainty are generally high when multiple stressors occur. In the face of high uncertainty, transparent communication is essential to inform decision-making. The impact of stressor interactions on risk and uncertainty was assessed using generalized linear models for additive and multiplicative effect of six anthropogenic stressors on the abundance of estuarine macrofauna across New Zealand. Models that accounted for multiplicative stressor interactions demonstrated that non-additive effects dominated, had increased explanatory power (6 to 73 % relative increase between models), and thereby reduced the risk of unexpected ecological responses to stress. Secondly, 3D-plots provide important insights in the direction, magnitude and gradients of change, and aid transparency and communication of complex stressor effects. Notably, small changes in a stressor can cause a disproportionally steep gradient of change for a synergistic effect where the tolerance to stressors are lost, and would invoke precautionary management. 3D-plots were able to clearly identify directional shifts where the nature of the interaction changed from antagonistic to synergistic along increasing stressor gradients. For example, increased nitrogen load and exposure caused a shift from positive to negative effect on the abundance of a deposit-feeding polychaete (Magelona). Assessments relying on model coefficient estimates, which provide one effect term, could not capture the complexities observed in 3D-plots and are at risk of mis-identifying interaction types. Finally, visualising model uncertainty demonstrated that although error terms were higher for multiplicative models, they better captured the uncertainty caused by data availability. Together, the steep gradients of change identified in 3D-plots and the higher uncertainty in model predictions in multiplicative models urges more conservative limits to be set for management that account for risk and uncertainty in multiple stressor effects.

Cumulative stressors reduce the self-regulating capacity of coastal ecosystems.

Marine ecosystems are prone to tipping points, particularly in coastal zones where dramatic changes are associated with interactions between cumulative stressors (e.g., shellfish harvesting, eutrophication and sediment inputs) and ecosystem functions. A common feature of many degraded estuaries is elevated turbidity that reduces incident light to the seafloor, resulting from multiple factors including changes in sediment loading, sea-level rise and increased water column algal biomass. To determine whether cumulative effects of elevated turbidity may result in marked changes in the interactions between ecosystem components driving nutrient processing, we conducted a large-scale experiment manipulating sediment nitrogen concentrations in 15 estuaries across a national-scale gradient in incident light at the seafloor. We identified a threshold in incident light that was related to distinct changes in the ecosystem interaction networks (EIN) that drive nutrient processing. Above this threshold, network connectivity was high with clear mechanistic links to denitrification and the role of large shellfish in nitrogen processing. The EIN analyses revealed interacting stressors resulting in a decoupling of ecosystem processes in turbid estuaries with a lower capacity to denitrify and process nitrogen. This suggests that, as turbidity increases with sediment load, coastal areas can be more vulnerable to eutrophication. The identified interactions between light, nutrient processing and the abundance of large shellfish emphasizes the importance of actions that seek to manage multiple stressors and conserve or enhance shellfish abundance, rather than actions focusing on limiting a single stressor.