Stressors increase the impacts of coastal macrofauna biodiversity loss on ecosystem multifunctionality

There is substantial evidence that biodiversity underpins ecosystem functioning, but it is unclear how these relationships change with multiple stressors in complex real-world settings. Coastal zones are affected by numerous stressors (for example, sediment input and nutrient runoff from land) and the cumulative effects of these stressors may result in pronounced and unexpected changes in the functioning of ecosystems. To investigate the cumulative effects of turbidity and elevated nutrients on coastal biodiversity-ecosystem functioning relationships, we performed a large-scale field experiment manipulating in situ sediment porewater ammonium concentrations and measured multiple ecosystem functions related to carbon fixation and mineralisation in 15 estuaries with varying levels of turbidity. The results indicated that the benthic macrofauna diversity (species richness, abundance, and functional richness) declined with increased porewater ammonium concentrations and there were clear thresholds in light levels at the seafloor in relation to the biodiversity-ecosystem function relationships. Multifunctionality indices (an integrated index of all measured functions) in moderately turbid and turbid estuaries (daily mean seafloor PAR  420 µmol m−2 s−1) however remained relatively constant, suggesting that they were more resilient against the nutrient-induced biodiversity loss. Our results demonstrate that ecosystems already stressed by stressors that alter functional performance (turbidity) may be more prone to loss of overall functioning if biodiversity is reduced by another stressor (nutrient enrichment), highlighting the potential snowballing effects of cumulative change.