Macrofaunal Functional Diversity Provides Resilience to Nutrient Enrichment in Coastal Sediments

Salt marsh macrofauna: An overview of functions and services

Salt marshes are globally important blue carbon ecosystems, providing essential services such as coastal protection, carbon sequestration, nutrient cycling, and biodiversity support. Among their key inhabitants, macrofauna play critical roles in sustaining ecosystem health and resilience through processes like bioturbation, nutrient cycling, organic matter turnover, and trophic interactions, which in turn support ecosystem services such as fisheries and coastal community livelihoods. Despite their contributions, no comprehensive review has yet focused exclusively on the diverse roles and services of salt marsh macrofauna. This review aims to address this gap by synthesizing current research, supported by a bibliometric analysis revealing significant growth in studies since the year 2000, especially those addressing ecosystem services and climate resilience. We provide an in-depth assessment of macrofaunal functions in bioturbation, nutrient cycling, organic matter dynamics, greenhouse gas regulation, primary and secondary production, and food web interactions. Additionally, we examine the ecosystem services provided, such as provisioning, regulating, and cultural services, and explore the impact of environmental stressors on macrofaunal communities. Finally, this review identifies significant knowledge gaps, offering strategic insights for future research and serving as a vital reference for advancing coastal management and salt marsh conservation strategies.

Effects of in situ experimental warming on metabolic expression in a soft sediment bivalve

Ocean surface temperatures and the frequency and intensity of marine heatwaves are increasing worldwide. Understanding how marine organisms respond and adapt to heat pulses and the rapidly changing climate is crucial for predicting responses of valued species and ecosystems to global warming. Here, we carried out an in situ experiment to investigate sublethal responses to heat spikes of a functionally important intertidal bivalve, the venerid clam Austrovenus stutchburyi. We describe changes in metabolic responses under two warming scenarios (five days and seven days) at two sites (muddy and sandy). Tidal flat warming during every low tide for five days affected the abundance of multiple functional metabolites within this species. The metabolic response was related to pathways such as metabolic energetics, amino acid and lipid metabolism, and accumulation of stress-related metabolites. There was some recovery after cooler weather during the final two days of the experiment. The degree of change was greater in muddy versus sandy sediments. Our findings provide new evidence of the metabolomic response of these important bivalve to heat stress, which could be used for resource managers when implementing strategies to mitigate the impacts of climate change on valuable marine resources.

Climatic, oceanic, freshwater, and local environmental drivers of New Zealand estuarine macroinvertebrates

Understanding the responses of organisms to different environmental drivers is critical for improving ecosystem management and conservation. Estuarine ecosystems are under pressure from multiple anthropogenic stressors (e.g. increasing sediment and nutrient loads, pollution, climate change) that are affecting the functions and services these ecosystems provide. Here, we used long-term estuarine benthic invertebrate monitoring data (∼30 year time-series) to evaluate the responses of macrobenthic invertebrate communities and indicator species to climatic, oceanic, freshwater, and local environmental drivers in New Zealand estuaries. We aimed to improve our ability to predict ecosystem change and understand the effects of multiple environment drivers on benthic communities. Our analyses showed that the abundance and richness of macrobenthic fauna and four indicator taxa (bivalves known to have differing tolerances to sediment mud content: Austrovenus stutchburyi, Macomona liliana, Theora lubrica, and Arthritica bifurca) responded to unique combinations of multiple environmental drivers across sites and times. Macrobenthic responses were highly mixed (i.e., positive and negative) and site-dependent. We also show that responses of macrobenthic fauna were lagged and most strongly related to climatic and oceanic drivers. The way the macrobenthos responded has implications for predicting and understanding the ecological consequences of a rapidly changing environment and how we conserve and manage coastal ecosystems.

Nitrogen enrichment changed the biogeochemical role of sesarmid crabs by shifting their diets in tropical mangrove ecosystems

Sesarmid crabs modulate nutrient dynamics of tropical mangroves through their leaf-eating habit. How N enrichment may alter this regulatory role, and the implications for mangrove nutrient dynamics, remain unclear. Using a mesocosm experiment, we tested how N enrichment could change the microphytobenthos (MPB) communities, thus modifying the crabs' diet and their role in nutrient dynamics. The factorial experiment combined with field investigation revealed a significant increase in the relative abundance of cyanobacteria. Stable isotope analysis suggested that the main carbon source of crabs shifted from leaf litter to cyanobacteria in mesocosms under both high (20×) and low (2×) N enrichment treatments. The significantly lower total cellulase activity of crabs in the mesocosms might explain the decreased carbon assimilation from leaf litter. The changes in the MPB and the microbiome with N enrichment in the presence of crabs may drive significantly higher carbon processing rate in tropical mangroves.

NZTD - The New Zealand Trait Database for shallow-water marine benthic invertebrates

Macrobenthic traits, for example feeding mode, life history, morphology, are increasingly used for determining responses of macrobenthic fauna to environmental change and influences on ecosystem functioning. Yet, trait information is scarce or non-existent in several parts of the world, such as New Zealand. This deficit makes collecting trait data a difficult and time-consuming task, limiting its potential use in trait-based assessments. Here, we present the New Zealand Trait Database (NZTD) for marine benthic invertebrates, the first comprehensive assessment of macrobenthic traits in New Zealand. The NZTD provides trait information for more than 700 macrobenthic taxa, categorised by 18 traits and 77 trait modalities. The NZTD includes five freely downloadable datasets, (1) the macrobenthic trait dataset, with outcomes from a fuzzy coding procedure, (2) the trait source information, (3) the references by taxa, (4) the full references list, and (5) the full taxa list used in the NZTD. Establishing the NZTD closes the trait knowledge gap in New Zealand and facilitates future research applying trait-based approaches to New Zealand's coastal macrofauna.

Crustaceans functional diversity in mangroves and adjacent mudflats of the Persian Gulf and Gulf of Oman

The effects of mangrove vegetation on functional features of crustacean assemblages were investigated using the Biological Traits Analysis (BTA). The study was carried out at four major sites in the arid mangrove ecosystem of the Persian Gulf and Gulf of Oman. Samples of Crustacea and associated environmental variables were taken seasonally (February 2018 and June 2019) from two habitats: a vegetated area with both mangrove trees and pneumatophores, and an adjacent mudflat. In each site, functional traits for the species were assigned using seven categories based on bioturbation, adult mobility, feeding habits and life-strategy traits. The results showed that the crabs (i.e., Opusia indica, Nasima dotilliformis and Ilyoplax frater) were widely distributed across all sites and habitats. The vegetated habitats supported higher taxonomic diversity than the mudflats, which highlights the importance of mangrove structural complexity for crustacean assemblages. Species inhabiting vegetated habitats were characterized by stronger presence of conveyor building species, detritivore, predator, grazer, lecithotrophic larval development, body size of 50-100 mm, and swimmer trait modalities. While, mudflat habitat enhanced the occurrences of surface deposit feeder, planktotrophic larval development, body size of <5 mm, and life span of 2-5 years. The results of our study showed that taxonomic diversity increased from the mudflats to the mangrove vegetated habitats. However, functional diversity did not differ between habitats. The significant differences in species and functional trait composition were observed between vegetated habitats and adjacent mudflats, stressing that different habitats may harbor different species and trait sets, likely as a result of habitat complexity. The use of taxonomic and functional attributes generate complementary information that can helps us to reach more efficient conclusions in terms of biodiversity conservation and ecosystem functionality in mangrove ecosystems.

Estimating Effects of Sea Level Rise on Benthic Biodiversity and Ecosystem Functioning in a Large Meso-Tidal Coastal Lagoon

Estuaries are among the world's most productive ecosystems, but due to their geographic location, they are at the forefront of anthropogenic pressures. Sea level rise (SLR) is one major consequence of climate change that poses a threat to estuaries with extensive intertidal habitats. The ecological implications of intertidal habitat loss have been largely overlooked despite their likely significance. We aimed to address this knowledge gap by investigating how benthic macroinvertebrate communities and their contributions to ecosystem function are likely to respond to SLR. Based on a spatially extensive dataset (119 sites) from a large coastal lagoon, depth, sediment chlorophyll concentrations, mud content, and average current speed were identified as the main drivers of community compositional turnover. Shifts in benthic community structure and associated functional implications were then evaluated using depth as a proxy for SLR. Three main macrofaunal groups representing intertidal, shallow subtidal, and deep subtidal habitats were identified. Functional trait analysis indicated low functional redundancy for a key intertidal suspension-feeding bivalve (Austrovenus stutchburyi) and the lack of a shallow subtidal functional replacement should intertidal habitats become inundated. These findings strongly suggest SLR and the associated environmental changes will alter estuarine macroinvertebrate communities, with implications for future ecosystem function and resilience.

Spatial and seasonal variation in macrozoobenthic density, biomass and community composition in a major tropical intertidal area, the Bijagós Archipelago, West-Africa

The coastal intertidal ecosystem of the Bijagós Archipelago, Guinea-Bissau, one of the largest and most important in West Africa, sustains a considerable proportion of the migratory shorebird populations of the East Atlantic Flyway and operates as a nursery area for benthic fish in the region. The macrozoobenthos in these mudflats constitute the main food source for both groups so that spatial and temporal variation in their abundance and community composition is likely to influence the abundance and distribution of fish and birds. In this study we described the spatial and temporal dynamics in the density, biomass, and community composition of macrozoobenthos across six intertidal flats in three islands of the Bijagós Archipelago. Overall, the Bijagós Archipelago was characterised by a highly species-rich macrozoobenthic community, with ca. 88 taxa identified across all sites, reaching a mean density of 1871 ± 58.3 ind.m-2 (mean ± SE) and mean biomass of 5.65 ± 0.41 g of AFDM.m-2 (ash-free dry mass per m2), values much lower than what was described for nearby intertidal areas, namely the Band d'Arguin, Mauritania. Density and biomass of major macrozoobenthos classes (Bivalvia, Polychaeta, Malacostraca and Gastropoda) differed across sites and months, displaying an overall increase in density towards the final months of the dry season (March and April). Similarly, community composition also differed significantly between sites and throughout the season. The site with most distinct community composition (Adonga) supported low diversity and high abundance of a few bivalve species, whilst other two sites that hosted the most diverse communities, were also the most similar between them (Anrumai and Abu). These spatial and temporal patterns constitute an important baseline to improve knowledge of this intertidal ecosystem and will contribute towards a better understanding of the spatial and temporal distribution patterns of their consumers.

Altered nitrogen transformation pathways and a legacy of sediment organic matter enrichment

Estuaries are ecologically valuable ecosystems that process nutrients through complex biogeochemical processes. Here we identify drivers and inhibitors of nitrogen removal in unvegetated intertidal sandflats at six sites in Manukau Harbour (37° 2.00'S 174° 42.00'E) to quantify the exchange of solutes across the sediment-water interface, with nitrogen removal rates (NRR) measured at two contrasting sites (PI and CB) near and far, respectively, from an historical wastewater treatment plant. Solute fluxes were paired with sediment and macrofauna community data to identify drivers of ecosystem function. Fluxes of oxygen, nitrogen and phosphorous were found to vary among sites, with differences attributed to variation in sediment characteristics (grain size, chlorophyll a, organic content) and macrofauna community structure. Mean NRR was found to vary between sites (PI = 32.2 vs CB = 217.9 μmol N₂/m²/h), with bioturbating macrofauna (bivalves Austrovenus stutchburyi and Macomona liliana), microphytobenthic biomass, and exposure to nutrients likely key contributing drivers.

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) decreased with the loss of macrofauna biodiversity. Functioning in low-turbidity estuaries (daily mean 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.

Biological traits approaches in benthic marine ecology: Dead ends and new paths

Biological traits analysis (BTA) links community structure to both ecological functions and response to environmental drivers through species' attributes. In consequence, it has become a popular approach in marine benthic studies. However, BTA will reach a dead end if the scientific community does not acknowledge its current shortcomings and limitations: (a) uncertainties related to data origins and a lack of standardized reporting of trait information; (b) knowledge gaps on the role of multiple interacting traits on driving the organisms' responses to environmental variability; (c) knowledge gaps regarding the mechanistic links between traits and functions; (d) a weak focus on the spatial and temporal variability that is inherent to the trait expression of species; and, last but not least, (e) the large reliance on expert knowledge due to an enormous knowledge gap on the basic ecology of many benthic species. BTA will only reach its full potential if the scientific community is able to standardize and unify the reporting and storage of traits data and reconsider the importance of baseline observational and experimental studies to fill knowledge gaps on the mechanistic links between biological traits, functions, and environmental variability. This challenge could be assisted by embracing new technological advances in marine monitoring, such as underwater camera technology and artificial intelligence, and making use of advanced statistical approaches that consider the interactive nature and spatio-temporal variability of biological systems. The scientific community has to abandon some dead ends and explore new paths that will improve our understanding of individual species, traits, and the functioning of benthic ecosystems.

Long-term response of marine benthic fauna to thin-layer capping with powdered activated carbon in the Grenland fjords, Norway

The Grenland fjords in Norway have a long history of contamination by large emissions of dioxins and mercury. As a possible sediment remediation method in situ, thin-layer capping with powdered activated carbon (AC) mixed with clay was applied at two test sites at 30 m and 95 m depth in the Grenland fjords. This study presents long-term effects of the AC treatment on the benthic community structure, i.e. nine years after capping. Capping with AC significantly reduced the number of species, their abundance and biomass at the two test sites, compared to uncapped reference sites. At the more shallow site, the dominant brittle star species Amphiura filiformis disappeared shortly after capping and did not re-establish nine years after capping. At the deeper site, the AC treatment also caused long-lasting negative effects on the benthic community, but some recovery was observed after nine years. Ecological indices used to assess environmental status did not capture the impaired benthic communities caused by the capping. The present study is the first documentation of negative effects of powdered AC on marine benthic communities on a decadal scale. Our results show that the benefits of reduced contaminant bioavailability from capping with AC should be carefully weighed against the cost of long-term detrimental effects on the benthic community. More research is needed to develop a thin-layer capping material that is efficient at sequestering contaminants without being harmful to benthic species.

Nitrogen and phosphorus enrichment cause declines in invertebrate populations: a global meta-analysis

Human-driven changes in nitrogen (N) and phosphorus (P) inputs are modifying biogeochemical cycles and the trophic state of many habitats worldwide. These alterations are predicted to continue to increase, with the potential for a wide range of impacts on invertebrates, key players in ecosystem-level processes. Here, we present a meta-analysis of 1679 cases from 207 studies reporting the effects of N, P, and combined N + P enrichment on the abundance, biomass, and richness of aquatic and terrestrial invertebrates. Nitrogen and phosphorus additions decreased invertebrate abundance in terrestrial and aquatic ecosystems, with stronger impacts under combined N + P additions. Likewise, N and N + P additions had stronger negative impacts on the abundance of tropical than temperate invertebrates. Overall, the effects of nutrient enrichment did not differ significantly among major invertebrate taxonomic groups, suggesting that changes in biogeochemical cycles are a pervasive threat to invertebrate populations across ecosystems. The effects of N and P additions differed significantly among invertebrate trophic groups but N + P addition had a consistent negative effect on invertebrates. Nutrient additions had weaker or inconclusive impacts on invertebrate biomass and richness, possibly due to the low number of case studies for these community responses. Our findings suggest that N and P enrichment affect invertebrate community structure mainly by decreasing invertebrate abundance, and these effects are dependent on the habitat and trophic identity of the invertebrates. These results highlight the important effects of human-driven nutrient enrichment on ecological systems and suggest a potential driver for the global invertebrate decline documented in recent years.

Environmental DNA metabarcoding reveals estuarine benthic community response to nutrient enrichment - Evidence from an in-situ experiment

Nutrient loading is a major threat to estuaries and coastal environments worldwide, therefore, it is critical that we have good monitoring tools to detect early signs of degradation in these ecologically important and vulnerable ecosystems. Traditionally, bottom-dwelling macroinvertebrates have been used for ecological health assessment but recent advances in environmental genomics mean we can now characterize less visible forms of biodiversity, offering a more holistic view of the ecosystem and potentially providing early warning signals of disturbance. We carried out a manipulative nutrient enrichment experiment (0, 150 and 600 g N fertilizer m^-2) in two estuaries in New Zealand to assess the effects of nutrient loading on benthic communities. After seven months of enrichment, environmental DNA (eDNA) metabarcoding was used to examine the response of eukaryotic (18S rRNA), diatom only (rbcL) and bacterial (16S rRNA) communities. Multivariate analyses demonstrated changes in eukaryotic, diatom and bacterial communities in response to nutrient enrichment at both sites, despite differing environmental conditions. These patterns aligned with changes in macrofaunal communities identified using traditional morphological techniques, confirming concordance between disturbance indicators detected by eDNA and current monitoring approaches. Clear shifts in eukaryotic and bacterial indicator taxa were seen in response to nutrient loading while changes in diatom only communities were more subtle. Community changes were discernible between 0 and 150 g N m^-2 treatments, suggesting that estuary health assessment tools could be developed to detect early signs of degradation. Increasing variation in community structure associated with nutrient loading could also be used as an indicator of stress or approaching tipping points. This work represents a first step towards the development of molecular-based estuary monitoring tools, which could provide a more holistic and standardized approach to ecosystem health assessment with faster turn-around times and lower costs.

Unraveling ecosystem functioning in intertidal soft sediments: the role of density-driven interactions

Although they only occupy a relatively small portion of the surface of the planet, coastal habitats are some of the most productive and valued ecosystems in the world. Among these habitats, tidal flats are an important component of many harbours and estuaries, but their deterioration due to human activities poses a serious threat to biodiversity and ecosystem function. Benthic communities are usually arranged in patches dominated by key species with overlapping distributions. Understanding the ecological consequences of interactions between these species in transition zones where their habitats overlap is necessary in order to quantify their contribution to overall ecosystem functioning and to scale-up and generalize results. Spatial transition in abundance and the interaction of multiple factors that drive ecosystem function are complex processes that require real-world research. Through a multi-site mensurative experiment, we show that transition areas drive non-linear effects on biogeochemical fluxes that have important implications for quantifying overall functioning. In our study the main drivers of ecosystem function were the abundance of two large but functionally very different species rather than biodiversity per se. Furthermore, we demonstrate that the use of the biogenic features created by specific infaunal species at the sediment–water interface is a better predictor of ecosystem functioning than the density of the species per se, making this approach particularly appealing for large scale, mapping and monitoring studies.

Mapping the estuarine ecosystem service of pollutant removal using empirically validated boosted regression tree models

Humans rely on the natural environment and benefit from the goods and services provided by natural ecosystems. Quantification and mapping of ecosystem services (ES) is required to better protect valued ES benefits under pressure from anthropogenic activities. The removal of excess nitrogen, a recognized catchment-derived pollutant, by biota in estuarine soft sediments is an important ES that potentially ameliorates the development of eutrophication symptoms. Here, we quantified estuarine benthic sediment characteristics and denitrification enzyme activity (DEA), a proxy of inorganic N removal, at 109 sites in four estuaries to develop a general ("global") model for predicting DEA. Our initial global model for linking DEA and environmental characteristics had good explanatory power, with sediment mud content having the strongest influence on DEA (60%), followed by sediment organic matter content (≈35%) and sediment chlorophyll a content (≈5%). Predicted and empirically evaluated DEA values in a fifth estuary (Whitford, n = 90 validation sites) were positively correlated (r = 0.77), and the fit and certainty of the model (based on two types of uncertainty measures) increased further after the validation sites were incorporated into it. The model tended to underpredict DEA at the upper end of its range (at the muddier, more organically enriched sites), and the relative roles of the three environmental predictors differed in Whitford relative to the four previously sampled estuaries (reducing the explained deviance relative to the initial global model). Our detailed quantification of DEA and methodological description for producing empirically validated maps, complete with uncertainty information, represents an important first step in the construction of nutrient pollution removal ES maps for use in coastal marine spatial management. This technique can likely be adapted to map other ecosystem functions and ES proxies worldwide.

Nitrogen enrichment increases greenhouse gas emissions from emerged intertidal sandflats

Unvegetated, intertidal sandflats play a critical role in estuarine carbon and nutrient dynamics. However, these ecosystems are under increasing threat from anthropogenic stressors, especially nitrogen enrichment. While research in this area typically focuses on sediment-water exchanges of carbon and nutrients during tidal inundation, there remain significant gaps in our understanding of GHG (Greenhouse Gas) fluxes during tidal emergence. Here we use in situ benthic chambers to quantify GHG fluxes during tidal emergence and investigate the impact of nitrogen enrichment on these fluxes. Our results demonstrate significant differences in magnitude and direction of GHG fluxes between emerged and submerged flats, demonstrating the importance of considering tidal state when estimating GHG emissions from intertidal flats. These responses were related to differences in microphytobenthic and macrofaunal activity, illustrating the important role of ecology in mediating fluxes from intertidal flats. Our results further demonstrate that nitrogen enrichment of 600 gN m-2 was associated with, on average, a 1.65x increase in CO2 uptake under light (photosynthetically active) conditions and a 1.35x increase in CO2 emission under dark conditions, a 3.8x increase in CH4 emission and a 15x increase in N2O emission overall. This is particularly significant given the large area intertidal flats cover globally, and their increasing exposure to anthropogenic stressors.

Shifts in trait-based and taxonomic macrofauna community structure along a 27-year time-series in the south-eastern North Sea

Current research revealed distinct changes in ecosystem functions, and thus in ecosystem stability and resilience, caused by changes in community structure and diversity loss. Benthic species play an important role in benthic-pelagic coupling, such as through the remineralization of deposited organic material, and changes to benthic community structure and diversity have associated with changes in ecosystem functioning, ecosystem stability and resilience. However, the long-term variability of traits and functions in benthic communities is largely unknown. By using abundance and bioturbation potential of macrofauna samples, taken along a transect from the German Bight towards the Dogger Bank in May 1990 and annually from 1995 to 2017, we analysed the taxonomic and trait-based macrofauna long-term community variability and diversity. Taxonomic and trait-based diversity remained stable over time, while three different regimes were found, characterised by changes in taxonomic and trait-based community structure. Min/max autocorrelation factor analysis revealed the climatic variables sea surface temperature (SST) and North Atlantic Oscillation Index (NAOI), nitrite, and epibenthic abundance as most important environmental drivers for taxonomic and trait-based community changes.

Biodiversity breakpoints along stress gradients in estuaries and associated shifts in ecosystem interactions

Denitrification in coastal sediments can provide resilience to eutrophication in estuarine ecosystems, but this key ecosystem function is impacted directly and indirectly by increasing stressors. The erosion and loading of fine sediments from land, resulting in sedimentation and elevated sediment muddiness, presents a significant threat to coastal ecosystems worldwide. Impacts on biodiversity with increasing sediment mud content are relatively well understood, but corresponding impacts on denitrification are uncharacterised. Soft sediment ecosystems have a network of interrelated biotic and abiotic ecosystem components that contribute to microbial nitrogen cycling, but these components (especially biodiversity measures) and their relationships with ecosystem functions are sensitive to stress. With a large dataset spanning broad environmental gradients this study uses interaction network analysis to present a mechanistic view of the ecological interactions that contribute to microbial nitrogen cycling, showing significant changes above and below a stressor (mud) threshold. Our models demonstrate that positive biodiversity effects become more critical with a higher level of sedimentation stress, and show that effective ecosystem management for resilience requires different action under different scenarios.

The application of bioturbators for aquatic bioremediation: Review and meta-analysis

Human activities introduce significant contamination into aquatic systems that impact biodiversity and ecosystem function. Many contaminants accumulate, and remediation options are now required worldwide. One method for bioremediation involves the application of macrofauna to stimulate microbial ecosystem processes including contaminant removal. However, if we are to confidently apply such a technique, we need clarity on the effect of bioturbators on different contaminants and how these vary under different environmental scenarios. Here we used a systematic review and meta-analysis to analyse current knowledge on the activities of bioturbating macrofauna in contaminated sediments and quantify how bioturbation-bioremediation changes depend on the taxonomic group, the aquatic ecosystem and important environmental variables. Three common contaminant classes were reviewed and analysed: metals, nutrients (i.e. ammonia and phosphorous) and polycyclic aromatic hydrocarbons (PAH). In addition, meta-regressions were calculated to estimate the effect of environmental and experimental design variables on effect sizes. Meta-analytic results revealed that deeper burrowing and more active sediment surface animals (e.g. polychaetes) increased metal release from sediments, nutrients and oxygen uptake by microbial fractions in comparison to bioturbators that inhabit shallower depths in sediments. In addition, there was a different effect of bioturbators on response variables in different aquatic systems. Finally, bioturbator effects on nutrient and metal release appeared modulated by context-specific variables such as temperature, pH, sediment grain size, animal density and experimental duration. Our findings highlight critical knowledge gaps such as field applications, less studied macrobenthic fauna and the incorporation of molecular approaches. Our results provide the first quantitative synthesis of the effects of bioturbators on contaminant fate and the variables that need to be considered for the optimization of this method as a viable approach for sediment remediation and contaminant management in aquatic systems.