Search

Microplastics are everywhere, including marine sediment. In this study, we evaluated the degradation of polyester, rayon, and cotton sewing threads over nine months when buried in marine sediment in Waitematā Harbour, Auckland, New Zealand. Polyester tensile strength was tested pre- and post-burial to track changes over time. Fourier-transform infrared spectroscopy (FTIR) analysis enabled the examination of the change to the chemical structural integrity of the polyester molecules over time. After one month, rayon and cotton degraded and were invisible to the eye, while visible signs of polyester degradation were apparent after 6 months of burial. This was confirmed by both tensile strength testing and FTIR chemical analysis. While microplastic pollution remains a serious problem, these findings show that at least one type of common plastic does degrade when buried in marine sediments. This likely has implications for seafloor ecosystem functionality and provides hope for plastic circular economy infrastructure., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Marine soft sediments play crucial roles in global biogeochemical cycles and biodiversity. Yet, with organisms often hidden in the sediment, they pose challenges for effective monitoring and management. This study introduces a novel approach utilizing sediment microtopography as a proxy for ecosystem functioning and biodiversity. Combining field sampling, benthic chamber incubations, and advanced Structure-from-Motion photogrammetry techniques, we investigated the relationships between microtopographic features and various ecological parameters across diverse subtidal habitats. Our findings reveal strong associations between sediment microtopography and environmental variables, benthic fluxes, biodiversity metrics, and community functional traits, with microtopography consistently explaining more than 50% of the variance in the data. This research demonstrates the potential of sediment microtopography as a cost-effective and scalable tool for assessing soft-sediment ecosystem dynamics and informing conservation strategies. By providing insights into the links between habitat structure and ecological processes, this study advances our understanding of marine benthic ecology and opens new possibilities for habitat assessment applications worldwide., (© 2024 The Author(s). Ecological Applications published by Wiley Periodicals LLC on behalf of The Ecological Society of America.)

While extensive research has explored the effects of plastic pollution, ecosystem responses remain poorly quantified, especially in field experiments. In this study, we investigated the impact of polyester pollution, a prevalent plastic type, on coastal sediment ecosystem function. Strips of polyester netting were buried into intertidal sediments, and effects on sediment oxygen consumption and polyester additive concentrations were monitored over 72-days. Our results revealed a rapid reduction in the magnitude and variability of sediment oxygen consumption, a crucial ecosystem process, potentially attributed to the loss of the additive di(2-ethylhexyl) phthalate (DEHP) from the polyester material. DEHP concentrations declined by 89% within the first seven days of deployment. However, effects on SOC dissipated after 22 days, indicating a short-term impact and a quick recovery by the ecosystem. Our study provides critical insights into the immediate consequences of plastic pollution on ecosystem metabolism in coastal sediments, contributing to a nuanced understanding of the temporal variation of plastic pollution's multifaceted impacts. Additionally, our research sheds light on the urgent need for comprehensive mitigation strategies to preserve marine ecosystem functionality from plastic pollution impacts., (© 2024. The Author(s).)

As plastic pollution continues to accumulate at the seafloor, concerns around benthic ecosystem functionality heightens. This research demonstrates the systematic effects of polyester microfibers on seafloor organic matter consumption rates, an important benthic ecosystem function connected to multiple reactions and processes. We used a field-based assay to measure the loss of organic matter, both with and without polyester microfiber contamination. We identified sediment organic matter content, mud content, and mean grain size as the main drivers of organic matter consumption, however, polyester microfiber contamination decoupled ecosystem relationships and altered observed organic matter cycling dynamics. Organic matter consumption rates varied across horizontal and vertical spaces, highlighting that consumption and associated plastic effects are dependent on environmental heterogeneity at both small (within sites) and larger (between sites) scales. Our results emphasize the important role habitat heterogeneity plays in seafloor organic matter consumption and the associated effects of plastic pollution on ecosystem function., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Samantha Ladewig reports financial support was provided by The International Yachting Fellowship of Rotarians. Samantha Ladewig reports financial support was provided by Anonymous philanthropic donors, provide through the Microphytes and Microplastics project. If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Environmental management in coastal ecosystems has been challenged by the complex cumulative effects that occur when many small issues result in large ecological shifts. Current environmental management of these spaces focuses on identifying and limiting problematic stressors via a series of assessment techniques. Whilst there is a strong desire among managers to consider complexity in ecological responses to cumulative effects, current approaches for assessing risk focus on breaking down the issues into multiple cause and effect relationships. However, uncertainty arises when data and information for a place are limited, as is commonly the case, and this creates decision paralysis while more information is generated. Here, we discuss how ecological understanding of network interactions in coastal marine ecosystems can be used as a lens to bring together multiple lines of evidence and create actions. We list and describe four characteristics of marine ecosystem interaction networks including the possibility for; 1) indirect effects, 2) effects that emerge as stressor magnitude increases the number of network components implicated, 3) network interactions that amplify these indirect effects, and 4) feedbacks that reinforce or stabilise against indirect effects. We then link these four characteristics to three case studies of common coastal environmental issues to demonstrate how a general understanding of ecological interaction networks can enhance priorities for stressor management that can be applied even when specific data is limited., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier Ltd.)

Bioturbating species play an essential role in regulating nutrient cycling in marine sediments, but their interaction with microplastics (MP) remains poorly understood. Here we investigated the linkage between MP and ecosystem functioning using experimental observations of luminophore distribution in the sediment to parametrize bioturbation coefficients (D b ). this information as fed into a simplified transport-reaction model, allowing us to upscale our experimental results. We found that the composition of bioturbators modulated shifts in the ecosystem functioning under microplastic stress. Maldanid worms (Macroclymenella stewartensis), functionally deep burrowing and upward-conveyor belt feeders, became less active. The D b of M. stewartensis reduced by 25% with the addition of 0.002 g MP cm -2 at surface sediment, causing accumulation of organic matter in the oxic sediment zone and stimulating aerobic respiration by 18%. In contract, the tellinid bivalve Macomona liliana, functionally a surface -deposit feeder that excretes at depth, maintained particle mixing behaviour in MP-contaminated systems. This study provides a mechanistic insight into the impacts of MP and indicates that the functional role of bioturbating species should be involved in assessing the global impact of MP. The model allowed us to understand the broad-scale impact of MP on seafloor habitat., (© 2023. Springer Nature Limited.)

Response diversity increases the potential 'options' for ecological communities to respond to stress (i.e. response capacity). An indicator of community response diversity is the diversity of different traits associated with their capacity to be resistant to stress, to recover and to regulate ecosystem functions. We conducted a network analysis of traits using benthic macroinvertebrate community data from a large-scale field experiment to explore the loss of response diversity along environmental gradients. We elevated sediment nutrient concentrations (a process that occurs with eutrophication) at 24 sites (in 15 estuaries) with varying environmental conditions (water column turbidity and sediment properties). Macroinvertebrate community response capacity to nutrient stress was dependent on the baseline trait network complexity in the ambient community (i.e. non-enriched sediments). The greater the complexity of the baseline network, the less variable the network response to nutrient stress was; in contrast, more variable responses to nutrient stress occurred with simpler networks. Thus, stressors or environmental variables that shift baseline network complexity also shift the capacity for these ecosystems to respond to additional stressors. Empirical studies that explore the mechanisms responsible for loss of resilience are essential to inform our ability to predict changes in ecological states.

Emerging research shows that microplastic pollution could be impacting seafloor ecosystem function, but this has been primarily demonstrated without environmental and ecological context. This causes uncertainty in the real-world effects of microplastic pollution and leaves out essential information guiding policy and mitigation. In this study, we take a well-supported sampling design and statistical approach commonly employed in benthic ecology to evaluate real-world effects of microplastic pollution on coastal, benthic ecosystem function. We utilised environmental gradients in the Waitemata Harbour of Auckland, New Zealand to evaluate the importance of commonly assessed biological, chemical, and geological sediment variables and the characteristics of microplastic contaminants in driving essential ecosystem functions. Our results showed that models including microplastic terms were more accurate and explained more variability than those without microplastic terms, highlighting that microplastics impact real-world seafloor ecosystem function. Specifically, microplastic fibers significantly influenced oxygen flux (p < 0.03) and the diverse forms of microplastics (i.e., richness) significantly influenced ammonium flux (p < 0.02). Additionally, interactions between microplastic fiber concentrations and mollusc abundances significantly contributed to oxygen flux (p < 0.02). These results provide the first evaluation of in situ relationships between microplastics and ecosystem function. Even more importantly, this study suggests the value of environmental and ecological context for addressing microplastic impacts on benthic ecosystems and argues for further field examination., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier B.V. All rights reserved.)

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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Increasing responsiveness to anthropogenic climate change and the loss of global shellfish ecosystems has heightened interest in the carbon storage and sequestration potential of bivalve-dominated systems. While coastal ecosystems are dynamic zones of carbon transformation and change, current uncertainties and notable heterogeneity in the benthic environment make it difficult to ascertain the climate change mitigation capacity of ongoing coastal restoration projects aimed at revitalizing benthic bivalve populations. In this study we sought to distinguish between direct and indirect effects of subtidal green-lipped mussels (Perna canaliculus) on carbon cycling, and combined published literature with in-situ experiments from restored beds to create a carbon budget for New Zealand's shellfish restoration efforts. A direct summation of biogenic calcification, community respiration, and sediment processes suggests a moderate carbon efflux (+100.1 to 179.6 g C m -2  year -1 ) occurs as a result of recent restoration efforts, largely reflective of the heterotrophic nature of bivalves. However, an examination of indirect effects of restoration on benthic community metabolism and sediment dynamics suggests that beds achieve greater carbon fixation rates and support enhanced carbon burial compared to nearby sediments devoid of mussels. We discuss limitations to our first-order approximation and postulate how the significance of mussel restoration to carbon-related outcomes likely increases over longer timescales. Coastal restoration is often conducted to support the provisioning of many ecosystem services, and we propose here that shellfish restoration not be used as a single measure to offset carbon dioxide emissions, but rather used in tandem with other initiatives to recover a bundle of valued ecosystem services., (© 2022 The Authors. Global Change Biology published by John Wiley & Sons Ltd.)

Coastal vegetated habitats such as seagrasses are known to play a critical role in carbon cycling and the potential to mitigate climate change, as blue carbon habitats have been repeatedly highlighted. However, little information is known about the role of associated macrofauna communities on the dynamics of critical processes of seagrass carbon metabolism (e.g., respiration, turnover, and production). We conducted a field study across a spatial gradient of seagrass meadows involving variable environmental conditions and macrobenthic diversity to investigate (1) the relationship between macrofauna biodiversity and secondary production (i.e., consumer incorporation of organic matter per time unit), and (2) the role of macrofauna communities in seagrass organic carbon metabolism (i.e., respiration and primary production). We show that, although several environmental factors influence secondary production, macrofauna biodiversity controls the range of local seagrass secondary production. We demonstrate that macrofauna respiration rates were responsible for almost 40% of the overall seafloor community respiration. Macrofauna represented on average >25% of the total benthic organic C stocks, high secondary production that is likely to become available to upper trophic levels of the coastal food web. Our findings support the role of macrofauna biodiversity in maintaining productive ecosystems, implying that biodiversity loss due to ongoing environmental change yields less productive seagrass ecosystems. Therefore, the assessment of carbon dynamics in coastal habitats should include associated macrofauna biodiversity elements if we aim to obtain robust estimates of global carbon budgets required to implement management actions for the sustainable functioning of the world's coasts., (© 2022 The Authors. Ecology published by Wiley Periodicals LLC on behalf of The Ecological Society of America.)

Restoration projects are underway internationally in response to global declines in shellfish beds. As diverse biological assemblages underpin a variety of ecosystem services, understanding broader changes in biodiversity associated with mussel restoration becomes increasingly valuable to scientists and restoration practitioners. Studies generally show bivalve beds increase species richness and abundance, but results are scale-dependent and conditional on the mobility of specific communities observed. We examined biodiversity at multiple scales to determine how communities with varying levels of mobility are influenced by subtidal mussel restoration. Significant changes in assemblage structure were observed in both mobile fish and epifaunal communities, with enhanced species richness and total abundance of associated individuals. In contrast, we observed site-dependent effects of bivalve restoration on macrofaunal community structure and composition, with sheltered, harbour mussel bed communities numerically dominated by detritivores accustomed to organically enriched, muddy sediments. Sediment organic matter significantly increased within mussel beds, and distance-based linear models showed that sediment organic matter was an important predictor of macrofaunal assemblage structure on mussel beds, highlighting the significance of benthic-pelagic coupling and biodeposition to soft-sediment communities. This study contributes novel methods and ecological insights on the role of species mobility and site selection in structuring restoration outcomes, better informing future mussel restoration efforts aimed at emphasising functionally-driven ecosystem services., (© 2022. The Author(s).)

The approach of applying stressor load limits or thresholds to aid estuarine management is being explored in many global case studies. However, there is growing concern regarding the influence of multiple stressors and their cumulative effects on the functioning of estuarine ecosystems due to the considerable uncertainty around stressor interactions. Recognising that empirical data limitations hinder parameterisation of detailed models of estuarine ecosystem responses to multiple stressors (suspended sediment, sediment mud and metal content, and nitrogen inputs), an expert driven Bayesian network (BN) was developed and validated. Overall, trends in estuarine condition predicted by the BN model were well supported by field observations, including results that were markedly higher than random (71-84% concordance), providing confidence in the overall model dynamics. The general BN framework was then applied to a case study estuary to demonstrate the model's utility for informing management decisions. Results indicated that reductions in suspended sediment loading were likely to result in improvements in estuarine condition, which was further improved by reductions in sediment mud and metal content, with an increased likelihood of high abundance of ecological communities relative to baseline conditions. Notably, reductions in suspended sediment were also associated with an increased probability of high nuisance macroalgae and phytoplankton if nutrient loading was not also reduced (associated with increased water column light penetration). Our results highlight that if stressor limit setting is to be implemented, limits must incorporate ecosystem responses to cumulative stressors, consider the present and desired future condition of the estuary of interest, and account for the likelihood of unexpected ecological outcomes regardless of whether the experts (or empirical data) suggest a threshold has or has not been triggered., (Copyright © 2021. Published by Elsevier Ltd.)

Polar seafloor ecosystems are changing rapidly and dramatically, challenging previously held paradigms of extreme dynamical stability. Warming-related declines in polar sea ice are expected to alter fluxes of phytoplankton and under-ice algae to the seafloor. Yet, how changes in food flux cascade through to seafloor communities and functions remains unclear. We leveraged natural spatial and temporal gradients in summertime sea ice extent to better understand the trajectories and implications of climate-related change in McMurdo Sound, Antarctica. McMurdo Sound was expected to be one of the last coastal marine environments on Earth to be affected by planetary warming, but the situation may be changing. Comparing satellite observations of selected coastal sites in McMurdo Sound between 2010-2017 and 2002-2009 revealed more ice-free days per year, and shorter distances to open water during the warmest months each year, in the more recent period. Interdecadal Pacific Oscillation (IPO), Oceanic Niño Index (ONI) and Antarctic Oscillation (AAO) climate indices peaked concurrently between 2014 and 2017 when sea ice breakouts in McMurdo Sound were most spatially and temporally extensive. Increases in sediment chlorophyll a and phaeophytin content (indicating increased deposition of detrital algal food material) were recorded during 2014-2017 at three coastal study sites in McMurdo Sound following the major sea ice breakouts. Soft-sediment seafloor ecosystem metabolism (measured in benthic incubation chambers as dissolved oxygen and inorganic nutrient fluxes) was correlated with sediment algal pigment concentration. Epifaunal invertebrate density, particularly opportunistic sessile suspension feeders, and infaunal community composition also shifted with increased food supply. The ecological characteristics and functions measured at the food-poor sites shifted towards those observed at richer sites at a surprisingly fast pace. These results indicate the sensitivity of the benthos and shed light on Antarctic marine trophic cascades and trajectories of response of iconic high-latitude seafloor habitats to a warming climate., (© 2021 John Wiley & Sons Ltd.)

Coastal sediments, where microplastics (MPs) accumulate, support benthic microalgae (BMA) that contribute to ecosystem functions such as primary production, nutrient recycling and sediment biostabilization. The potential interactions between MPs, BMA and associated properties and functions remain poorly understood. To examine these interactions, a survey of 22 intertidal sites was conducted. MP abundance, size and a suite of MP diversity indices (based on color and shape) were determined from surface sediments alongside biochemical and physical properties. MPs were detected at all sites and dominated by polypropylene (34%), polyester (18%) and polyethylene (11%). Fragment and fiber dominance (16-92% and 6-81% respectively) and color-shape category diversity varied significantly by site. Distance-based linear models demonstrated that estuary-wide, mean grain size and mud were the best predictors of MP abundance-diversity matrices, but variance explained was low (9%). Relationships were improved when the data was split into sandy and muddy habitats. In sandy habitats (<8% mud), physical properties of the bed (mean grain size, mud content and distance from the estuary mouth) were still selected as predictors of MP abundance-diversity (14% variance explained); but a number of bivariate relationships were detected with biochemical properties such as BMA associated pigments and organic matter. In muddy habitats (>8% mud), porewater ammonium was lower when fiber abundance and overall MP diversity were higher. The inclusion of porewater ammonium, organic matter content and pheophytins alongside physical properties explained a greater percentage of the variance in MP abundance-diversity for muddy habitats (21%). The results highlight the importance of examining plastic shapes and MP categories in addition to abundance and emphasize that functionally different habitats should be examined separately to increase our understanding of MP-biota-function relationships., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Ecosystems simultaneously deliver multiple functions that relate to both the activities of resident species and environmental conditions. One of the biggest challenges in multifunctionality assessment is balancing analytical simplicity with ecosystem complexity. As an alternative to index-based approaches, we introduce a multivariate network analysis that uses network theory to assess multifunctionality in terms of the relationships between species' functional traits, environmental characteristics, and functions. We tested our approach in a complex and heterogeneous ecosystem, marine intertidal sandflats. We considered eight ecosystem function, five macrofaunal functional trait groups derived from 36 species, and four environmental characteristics. The indicators of ecosystem functions included the standing stock of primary producers, oxygen production, benthic oxygen consumption, DIN (ammonium and NOx efflux) and phosphate release from the sediments, denitrification, and organic matter degradation at the sediment surface. Trait clusters included functional groups of species that shared combinations of biological traits that affect ecosystem function: small mobile top 2 cm dwellers, suspension feeders, deep-dwelling worms, hard-bodied surface dwellers, and tube-forming worms. Environmental characteristics included sediment organic matter, %mud, %shell hash, and %sediment water content. Our results visualize and quantify how multiple ecosystem elements are connected and contribute to the provision of functions. Small mobile top 2 cm dwellers (among trait clusters) and %mud (among environmental characteristics) were the best predictor for multiple functions. Detailed knowledge of multifunctionality relationships can significantly increase our understanding of the real-world complexity of natural ecosystems. Multivariate network analysis, as a standalone method or applied alongside already existing single index multifunctionality methods, provides means to advance our understanding of how environmental change and biodiversity loss can influence ecosystem performance across multiple dimensions of functionality. Embedding such a detailed yet holistic multifunctionality assessment in environmental decision-making will support the assessment of multiple ecosystem services and social-ecological values., Competing Interests: There are no conflicts of interest., (© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Despite a long history of disturbance-recovery research, we still lack a generalizable understanding of the attributes that drive community recovery potential in seafloor ecosystems. Marine soft-sediment ecosystems encompass a range of heterogeneity from simple low-diversity habitats with limited biogenic structure, to species-rich systems with complex biogenic habitat structure. These differences in biological heterogeneity are a product of natural conditions and disturbance regimes. To search for unifying attributes, we explore whether a set of simple traits can characterize community disturbance-recovery potential using seafloor patch-disturbance experiments conducted in two different soft-sediment landscapes. The two landscapes represent two ends of a spectrum of landscape biotic heterogeneity in order to consider multi-scale disturbance-recovery processes. We consider traits at different levels of biological organization, from the biological traits of individual species, to the traits of species at the landscape scale associated with their occurrence across the landscape and their ability to be dominant. We show that in a biotically heterogeneous landscape (Kawau Bay, New Zealand), seafloor community recovery is stochastic, there is high species turnover, and the landscape-scale traits are good predictors of recovery. In contrast, in a biotically homogeneous landscape (Baltic Sea), the options for recovery are constrained, the recovery pathway is thus more deterministic and the scale of recovery traits important for determining recovery switches to the individual species biological traits within the disturbed patch. Our results imply that these simple, yet sophisticated, traits can be effectively used to characterize community recovery potential and highlight the role of landscapes in providing resilience to patch-scale disturbances., Competing Interests: The authors declare no conflicts of interest., (© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Ecologists have long acknowledged the importance of context dependency related to position along spatial gradients. It is also acknowledged that broad-scale climate patterns can directly and indirectly alter population dynamics. What is not often addressed is whether climate patterns such as the Southern Oscillation interact with population-level temporal patterns and affect the ability of time-series data, such as long-term state of the environment monitoring programmes, to detect change. Monitoring design criteria generally focus on number of data points, sampling frequency and duration, often derived from previous information on species seasonal and multi-year temporal patterns. Our study questioned whether the timing of any changes relative to Southern Oscillation, interacting with species populations dynamics, would also be important. We imposed a series of simulated reductions on macrofaunal abundance data collected regularly over 29 years from two sites, using species selected for observed differences in temporal dynamics. We hypothesized that (1) high within-year sampling frequency would increase detection ability for species with repeatable seasonality cycles and (2) timing of the reduction in abundance relative to the Southern Oscillation was only likely to affect detection ability for long-lived species with multi-year cyclic patterns in abundance. However, regardless of species population dynamics, we found both within-year sampling frequency and the timing of the imposed reduction relative to the Southern Oscillation Index affected detection ability. The latter result, while apparently demonstrating a confounding influence on monitoring, offers the opportunity to improve our ability to detect and interpret analyses of monitoring data, and thus our ability to make recommendations to managers., (© 2021 John Wiley & Sons Ltd.)

Plastic pollution continues to seep into natural and pristine habitats. Emerging laboratory-based research has evoked concern regarding plastic's impact on ecosystem structure and function, the essence of the ecosystem services that supports our life, wellbeing, and economy. These impacts have yet to be observed in nature where complex ecosystem interaction networks are enveloped in environmental physical and chemical dynamics. Specifically, there is concern that environmental impacts of plastics reach beyond toxicity and into ecosystem processes such as primary production, respiration, carbon and nutrient cycling, filtration, bioturbation, and bioirrigation. Plastics are popularly regarded as recalcitrant carbon molecules, although they have not been fully assessed as such. We hypothesize that plastics can take on similar roles as natural recalcitrant carbon (i.e., lignin and humic substances) in carbon cycling and associated biogeochemistry. In this paper, we review the current knowledge of the impacts of plastic pollution on marine, benthic ecosystem function. We argue for research advancement through (1) employing field experiments, (2) evaluating ecological network disturbances by plastic, and (3) assessing the role of plastics (i.e., a carbon-based molecule) in carbon cycling at local and global scales., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

The distribution of benthic ecosystems, dominated by filter-feeding communities, is highly influenced by the seabed geomorphology. However, the spatial variation in settlement of these species is also affected by near-bottom currents and any changes in light, nutrient concentration and food quality often associated with increases of suspended sediment concentrations within the water column. Detailed predictions of the geographic distribution of filter-feeder species and a deeper understanding of the physical processes influencing their distribution patterns is key for effective management and conservation. To date, predictive distribution modelling has been derived essentially from geomorphological parameters, mainly using spatially limited observations. In this study, seabed mapping, oceanographic modelling, hydrographic records and biological observations are integrated to provide high-resolution prediction of filter-feeder habitat distribution within Queen Charlotte Sound/Tōtaranui and Tory Channel/Kura Te Au, South Island of New Zealand. The aim is to evaluate potential suitable habitat areas for filter-feeders to inform where habitat restoration management should focus efforts to recover communities such as the horse mussel (Atrina zelandica) or the green-lipped mussel (Perna canaliculus), both of which have high economic impact in New Zealand. To accomplish this, Maximum Entropy (MaxEnt) predictive modelling was used to produce Habitat Suitability (HS) maps, using geomorphological parameters and seafloor classification information. Final HS maps also incorporated oceanographic and sediment dynamic information, showing that filter-feeder habitat distribution is highly influenced by the hydrodynamics and sedimentary processes apart from the seafloor geomorphology. Filter-feeder communities inhabit quiescent areas, limited by depth, slope and sediment type; and coincide with regions presenting low near-bottom currents and low turbidity levels. Additionally, the obtained results reveal the effects of the coastal settlements and major marine traffic routes, limiting the suitable habitats to areas with less human impact. This study demonstrates that a multidisciplinary approach is crucial to better predict the spatial distribution of benthic communities, which is key to improve benthic habitat restoration and recovery assessments., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

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., (© 2020 The Authors. Ecological Applications published by Wiley Periodicals LLC on behalf of Ecological Society of America.)

The effects of microplastics (MPs) on the ecological functioning in marine sediments is largely unknown. However, coastal marine sediments and their resident communities play critical roles in the transformation of organic matter and the cycling of nutrients that influence both local and global processes. To investigate how microplastics influence ecosystem functions associated with sediment biogeochemistry, large bivalves and microphytobenthos, we conducted a 31-day laboratory experiment. The experiment tested the role of micro-polyethylene terephthalate (mPETs) at five concentrations (0%, 1%, 3%, 6%, and 8% based on wet weight of top 1 cm sediment). Canonical principle of coordinate analysis (CAP) was applied to assess change on the ecosystem functionality with increasing concentrations of mPETs. Our results highlight that stress effects on ecosystem function are the product of the interaction between Macomona liliana, microphytobenthos and mPETs., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

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.

Background: Species distribution models are commonly used tools to describe diversity patterns and support conservation measures. There is a wide range of approaches to developing SDMs, each highlighting different characteristics of both the data and the ecology of the species or assemblages represented by the data. Yet, signals of species co-occurrences in community data are usually ignored, due to the assumption that such structuring roles of species co-occurrences are limited to small spatial scales and require experimental studies to be detected. Here, our aim is to explore associations among marine sandy-bottom sediment inhabitants and test for the structuring effect of seagrass on co-occurrences among these species across a New Zealand intertidal sandflat, using a joint species distribution model (JSDM)., Results: We ran a JSDM on a total of 27 macrobenthic species co-occurring in 300,000 m 2 of sandflat. These species represented all major taxonomic groups, i.e. polychaetes, bivalves and crustaceans, collected in 400 sampling locations. A number of significant co-occurrences due to shared habitat preferences were present in vegetated areas, where negative and positive correlations were approximately equally common. A few species, among them the gastropods Cominella glandiformis and Notoacmea scapha, co-occurred randomly with other seagrass benthic inhabitants. Residual correlations were less apparent and mostly positive. In bare sand flats shared habitat preferences resulted in many significant co-occurrences of benthic species. Moreover, many negative and positive residual patterns between benthic species remained after accounting for habitat preferences. Some species occurring in both habitats showed similarities in their correlations, such as the polychaete Aglaophamus macroura, which shared habitat preferences with many other benthic species in both habitats, yet no residual correlations remained in either habitat., Conclusions: Firstly, analyses based on a latent variable approach to joint distributions stressed the structuring role of species co-occurrences beyond experimental scales. Secondly, results showed context dependent interactions, highlighted by species having more interconnected networks in New Zealand bare sediment sandflats than in seagrass meadows. These findings stress the critical importance of natural history to modelling, as well as incorporating ecological reality in SDMs.

Microfibers often dominate sediment microplastic samples, but little is known about their ecological effects on benthic organisms and functions. Polyethylene terephthalate) (PET) microfibers were added to 36 sediment chambers at six concentrations (0-0.5 g kg -1 sediment) to assess the effects on microphytobenthos (MPB), a key deposit-feeding bivalve, Macomona liliana , and sediment nutrient pools. MPB photosynthesis was promoted in 18 chambers through a 12 h light/dark cycle. Another 18 chambers were maintained under dark conditions to inhibit photosynthesis. After 35 days of MPB growth and stabilization, four M. liliana were added to each chamber for a further 40 days. MPB biomass and composition were examined alongside M. liliana biochemical and behavioral properties and porewater dissolved inorganic nutrient concentrations. Increasing microfibers resulted in lower MPB biomass, fewer diatom-associated fatty acids (FAs), and an increase in cyanobacteria. The changes in MPB coincided with up to 75% lower energy reserves and reduced burrowing activity in M. liliana . In the light, nitrate + nitrate (NO x ) was significantly elevated and related to M. liliana and MPB biochemical properties. Ammoniu (NH 4 + ) concentrations increased but were variable in both the light and the dark. Our results suggest that increasing microfiber concentrations influence the interactions between M. liliana and MPB and affect biogeochemical processing in coastal marine sediments.

Estuaries are subject to disturbance by land-based sediment and nutrient inputs, resulting in changes to the ecosystems and the functions and services that they support. Spatial mapping tools that identify how functional hotspots in the estuary may shift in location and intensity under different disturbance scenarios highlight to managers the trajectory of change and the value of active management and restoration, but to date these tools are only available in the most intensively researched ecosystems. Using empirical data derived from long-term monitoring and multi-habitat field experiments we developed future scenarios representing different impacts of environmental degradation on estuarine ecosystem functions that are important for supporting ecosystem services. We used the spatial prioritization software Zonation in a novel fashion to assess effects of different disturbance scenarios on critical soft-sediment ecosystem processes (nutrient fluxes and sediment erodibility measures) that are influenced by macrofaunal communities and local environment conditions. We compared estimates of current conditions with three scenarios linked to changes in land-use and resulting downstream impacts on estuarine ecosystems to determine how disturbance influences the distribution of high value areas for ecosystem function. Scenarios investigated the implications of habitat degradation associated with sediment deposition and declines in large sediment-dwelling animal abundance whose behavior has important influences on ecosystem function. Our analyses demonstrate decreases in the majority of ecosystem processes under scenarios associated with disturbances. These results suggest that it is important to restore biodiversity and ecosystem function and that the application of Zonation in this context offers a simple, rapid and cost-effective way of identifying priority actions and locations for restoration, and how these shift due to multiple impacts., (© 2020 by the Ecological Society of America.)

A major challenge in ecology and environmental management is linking changes in community composition to ecosystem functions. We developed the network analysis of traits (NAT) to show changes in community network structure based on the changes in the composition and connectivity between clusters of species that share traits that imply shifts in functional diversity. We tested the application of NAT on a 113 species found on an intertidal sandflat that was subject to experimental nitrogen addition (control [0 g N/m 2 ], medium [150 g N/m 2 ], and high [600 g N/m 2 ]). This allowed us to directly link mechanistic changes in community composition and function with the trait-space network patterns revealed by NAT. We demonstrate that under medium (150 g N/m 2 ) N treatment, functional diversity remained consistent, whereas increasing disturbance to high (600 g N/m 2 ) N treatment affected the species-trait network structure and caused merging of functional clusters implying a loss of functional trait diversity., (© 2019 by the Ecological Society of America.)

Under globally accelerating rates of ecosystem degradation, maintaining ecosystem function is a priority to avoid loss of valuable ecosystem services. Two factors are important: changes to the disturbance regime (stresses imposed) and resilience of biodiversity and ecosystem functions (the ecosystem's capacity to respond to change). Various attributes at different scales of ecological organisation confer resilience (from individual species to communities at landscape scales), but it is critical to understand how these attributes interact to inform how ecosystem function changes with disturbances that vary in intensity, spatial extent, and frequency. Individual species attributes influence their resistance, while attributes at the landscape-scale influence recovery of communities and function. Understanding resilience to disturbances requires defining the characteristics of a resilient community at multiple scales., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Increasingly studies are reporting sudden and dramatic changes in the structure and function of communities or ecosystems. The prevalence of these reports demonstrates the importance for management of being able to detect whether these have happened and, preferably, whether they are likely to occur. Ecological theory provides the rationale for why such changes occur and a variety of statistical indicators of approach that have generic properties have been developed. However, whether the theory has successfully translated into monitoring programmes is unknown. We searched the literature for guidelines that would drive design of monitoring programmes able to detect past and approaching tipping points and analysed marine monitoring programmes in New Zealand. We found very few guidelines in the ecological, environmental or monitoring literature, although both simulation and marine empirical studies suggest that within-year sampling increases the likelihood of detecting approaching tipping points. The combination of the need to monitor both small and medium scale temporal dynamics of multiple variables to detect tipping points meant that few marine monitoring programmes in New Zealand were fit for that purpose. Interestingly, we found many marine examples of studies detecting past and approaching TP with fewer data than was common in the theoretical literature. We, therefore, suggest that utilizing ecological knowledge is of paramount importance in designing and analyzing time-series monitoring for tipping points and increasing the certainty for short-term or infrequent datasets of whether a tipping point has occurred. As monitoring plays an important role in management of tipping points by providing supporting information for other locations about when and why a tipping point may occur, we believe that monitoring for tipping points should be promoted., (Copyright © 2018. Published by Elsevier Ltd.)

Polar ecosystems are bellwether indicators of climate change and offer insights into ecological resilience. In this study, we describe contrasting responses to an apparent regime shift of two very different benthic communities in McMurdo Sound, Antarctica. We compared species-specific patterns of benthic invertebrate abundance and size between the west (low productivity) and east (higher productivity) sides of McMurdo Sound across multiple decades (1960s-2010) to depths of 60 m. We present possible factors associated with the observed changes. A massive and unprecedented shift in sponge recruitment and growth on artificial substrata observed between the 1980s and 2010 contrasts with lack of dramatic sponge settlement and growth on natural substrata, emphasizing poorly understood sponge recruitment biology. We present observations of changes in populations of sponges, bryozoans, bivalves, and deposit-feeding invertebrates in the natural communities on both sides of the sound. Scientific data for Antarctic benthic ecosystems are scant, but we gather multiple lines of evidence to examine possible processes in regional-scale oceanography during the eight years in which the sea ice did not clear out of the southern portion of McMurdo Sound. We suggest that large icebergs blocked currents and advected plankton, allowed thicker multi-year ice, and reduced light to the benthos. This, in addition to a possible increase in iron released from rapidly melting glaciers, fundamentally shifted the quantity and quality of primary production in McMurdo Sound. A hypothesized shift from large to small food particles is consistent with increased recruitment and growth of sponges on artificial substrata, filter-feeding polychaetes, and some bryozoans, as well as reduced populations of bivalves and crinoids that favor large particles, and echinoderms Sterechinus neumayeri and Odontaster validus that predominantly feed on benthic diatoms and large phytoplankton mats that drape the seafloor after spring blooms. This response of different guilds of filter feeders to a hypothesized shift from large to small phytoplankton points to the enormous need for and potential value of holistic monitoring programs, particularly in pristine ecosystems, that could yield both fundamental ecological insights and knowledge that can be applied to critical conservation concerns as climate change continues., (© 2018 The Authors. Ecological Applications published by Wiley Periodicals, Inc. on behalf of Ecological Society of America.)

Metacommunity theory highlights the potential of β-diversity as a useful link to empirical research, especially in diverse systems where species exhibit a range of stage-dependent dispersal characteristics. To investigate the importance of different components and scales of β-diversity in community assembly, we conducted a large-scale disturbance experiment and compared relative recovery across multiple sites and among plots within sites on the rocky shore. Six sites were spread along 80 km of coastline and, at each site, five plots were established, matching disturbed and undisturbed quadrats. Recovery was not complete at any of the sites after 1 year for either epibenthos (mostly composed of macroalgae and, locally, mussels) or infauna. Significant differences in recovery among sites were observed for epibenthos but not for infauna, suggesting that different community assembly processes were operating. This was supported by epibenthos in the recovering plots having higher species turnover than in undisturbed sediment, and recovery well predicted by local diversity, while infaunal recovery was strongly influenced by the epibenthic community's habitat complexity. However, infaunal community recovery did not simply track formation of habitat by recovering epibenthos, but appeared to be overlain by within-site and among-site aspects of infaunal β-diversity. These results suggest that documenting changes in the large plants and animals alone will be a poor surrogate for rocky shore community assembly processes. No role for ecological connectivity (negative effect of among-site β-diversity) in driving recovery was observed, suggesting a low risk of effects from multiple disturbances propagating along the coast, but a limited resilience at the site scale to large-scale disturbances such as landslides or oil spills., Competing Interests: The authors have no competing interests.

The importance of long-term environmental monitoring and research for detecting and understanding changes in ecosystems and human impacts on natural systems is widely acknowledged. Over the last decades, a number of critical components for successful long-term monitoring have been identified. One basic component is quality assurance/quality control protocols to ensure consistency and comparability of data. In Norway, the authorities require environmental monitoring of the impacts of the offshore petroleum industry on the Norwegian continental shelf, and in 1996, a large-scale regional environmental monitoring program was established. As a case study, we used a sub-set of data from this monitoring to explore concepts regarding best practices for long-term environmental monitoring. Specifically, we examined data from physical and chemical sediment samples and benthic macroinvertebrate assemblages from 11 stations from six sampling occasions during the period 1996-2011. Despite the established quality assessment and quality control protocols for this monitoring program, we identified several data challenges, such as missing values and outliers, discrepancies in variable and station names, changes in procedures without calibration, and different taxonomic resolution. Furthermore, we show that the use of different laboratories over time makes it difficult to draw conclusions with regard to some of the observed changes. We offer recommendations to facilitate comparison of data over time. We also present a new procedure to handle different taxonomic resolution, so valuable historical data is not discarded. These topics have a broader relevance and application than for our case study.

In coastal ecosystems, climate change affects multiple environmental factors, yet most predictive models are based on simple cause-and-effect relationships. Multiple stressor scenarios are difficult to predict because they can create a ripple effect through networked ecosystem functions. Estuarine ecosystem function relies on an interconnected network of physical and biological processes. Estuarine habitats play critical roles in service provision and represent global hotspots for organic matter processing, nutrient cycling and primary production. Within these systems, we predicted functional changes in the impacts of land-based stressors, mediated by changing light climate and sediment permeability. Our in-situ field experiment manipulated sea level, nutrient supply, and mud content. We used these stressors to determine how interacting environmental stressors influence ecosystem function and compared results with data collected along elevation gradients to substitute space for time. We show non-linear, multi-stressor effects deconstruct networks governing ecosystem function. Sea level rise altered nutrient processing and impacted broader estuarine services ameliorating nutrient and sediment pollution. Our experiment demonstrates how the relationships between nutrient processing and biological/physical controls degrade with environmental stress. Our results emphasise the importance of moving beyond simple physically-forced relationships to assess consequences of climate change in the context of ecosystem interactions and multiple stressors.

Declining biodiversity and loss of ecosystem function threatens the ability of habitats to contribute ecosystem services. However, the form of the relationship between biodiversity and ecosystem function (BEF) and how relationships change with environmental change is poorly understood. This limits our ability to predict the consequences of biodiversity loss on ecosystem function, particularly in real-world marine ecosystems that are species rich, and where multiple ecosystem functions are represented by multiple indicators. We investigated spatial variation in BEF relationships across a 300 000 m 2 intertidal sandflat by nesting experimental manipulations of sediment pore water nitrogen concentration into sites with contrasting macrobenthic community composition. Our results highlight the significance of many different elements of biodiversity associated with environmental characteristics, community structure, functional diversity, ecological traits or particular species (ecosystem engineers) to important functions of coastal marine sediments (benthic oxygen consumption, ammonium pore water concentrations and flux across the sediment-water interface). Using the BEF relationships developed from our experiment, we demonstrate patchiness across a landscape in functional performance and the potential for changes in the location of functional hot and cold spots with increasing nutrient loading that have important implications for mapping and predicating change in functionality and the concomitant delivery of ecosystem services., (© 2017 The Author(s).)

Adding fertiliser to sediments is an established way of studying the effects of eutrophication but a lack of consistent methodology, reporting on enrichment levels, or guidance on application rates precludes rigorous synthesis and meta-analysis. We developed a simple enrichment technique then applied it to 28 sites across an intertidal sandflat. Fertiliser application rates of 150 and 600gNm(-2) resulted in pore water ammonium concentrations respectively 1-110 and 4-580×ambient, with greater elevations observed in deeper (5-7cm) than surface (0-2cm) sediments. These enrichment levels were similar to eutrophic estuaries and were maintained for at least seven weeks. The high between-site variability could be partially explained by the sedimentary environment and macrofaunal community (42%), but only at the high application rate. We suggest future enrichment studies should be conducted in situ across large environmental gradients to incorporate real world complexity and increase generality of conclusions., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Many conservation actions are justified on the basis of managing biodiversity. Biodiversity, in terms of species richness, is largely the product of rare species. This is problematic because the intensity of sampling needed to characterize communities and patterns of rarity or to justify the use of surrogates has biased sampling in favor of space over time. However, environmental fluctuations interacting with community dynamics lead to temporal variations in where and when species occur, potentially affecting conservation planning by generating uncertainty about results of species distribution modeling (including range determinations), selection of surrogates for biodiversity, and the proportion of biodiversity composed of rare species. To have confidence in the evidence base for conservation actions, one must consider whether temporal replication is necessary to produce broad inferences. Using approximately 20 years of macrofaunal data from tidal flats in 2 harbors, we explored variation in the identity of rare, common, restricted range, and widespread species over time and space. Over time, rare taxa were more likely to increase in abundance or occurrence than to remain rare or disappear and to exhibit temporal patterns in their occurrence. Space-time congruency in ranges (i.e., spatially widespread taxa were also temporally widespread) was observed only where samples were collected across an environmental gradient. Fifteen percent of the taxa in both harbors changed over time from having spatially restricted ranges to having widespread ranges. Our findings suggest that rare species can provide stability against environmental change, because the majority of species were not random transients, but that selection of biodiversity surrogates requires temporal validation. Rarity needs to be considered both spatially and temporally, as species that occur randomly over time are likely to play a different role in ecosystem functioning than those exhibiting temporal structure (e.g., seasonality). Moreover, temporal structure offers the opportunity to place management and conservation activities within windows of maximum opportunity., (© 2016 Society for Conservation Biology.)

Marine organisms are known to play important roles in transforming nutrients in sediments, however, guidelines to optimize sediment restoration are not available. We conducted a laboratory mesocosm experiment to investigate the role of hard clams, polychaetes, the degree of physical disturbance and denitrifying bacterial concentrations in removing total nitrogen (TN), total phosphorus (TP), and total organic carbon (TOC) in marine sediments. Response surface methodology was employed to analyze the results of initial experiments and in a subsequent experiment identified optimal combinations of parameters. Balancing the TN, TP, TOC removal efficiency, our model predicted 39% TN removal, 33% TP removal, and 42% TOC removal for a 14-day laboratory bioremediation trial using hard clams biomass of 1.2kgm(-2), physical disturbance depth of 16.4cm, bacterial density of 0.18Lm(-2), and polychaetes biomass of 0.16kgm(-2), respectively. These results emphasize the value of combining different species in field-based bioremediation., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Global climate change will undoubtedly be a pressure on coastal marine ecosystems, affecting not only species distributions and physiology but also ecosystem functioning. In the coastal zone, the environmental variables that may drive ecological responses to climate change include temperature, wave energy, upwelling events and freshwater inputs, and all act and interact at a variety of spatial and temporal scales. To date, we have a poor understanding of how climate-related environmental changes may affect coastal marine ecosystems or which environmental variables are likely to produce priority effects. Here we use time series data (17 years) of coastal benthic macrofauna to investigate responses to a range of climate-influenced variables including sea-surface temperature, southern oscillation indices (SOI, Z4), wind-wave exposure, freshwater inputs and rainfall. We investigate responses from the abundances of individual species to abundances of functional traits and test whether species that are near the edge of their tolerance to another stressor (in this case sedimentation) may exhibit stronger responses. The responses we observed were all nonlinear and some exhibited thresholds. While temperature was most frequently an important predictor, wave exposure and ENSO-related variables were also frequently important and most ecological variables responded to interactions between environmental variables. There were also indications that species sensitive to another stressor responded more strongly to weaker climate-related environmental change at the stressed site than the unstressed site. The observed interactions between climate variables, effects on key species or functional traits, and synergistic effects of additional anthropogenic stressors have important implications for understanding and predicting the ecological consequences of climate change to coastal ecosystems., (© 2015 John Wiley & Sons Ltd.)

Seafloor biodiversity is a key mediator of ecosystem functioning, but its role is often excluded from global budgets or simplified to black boxes in models. New techniques allow quantification of the behavior of animals living below the sediment surface and assessment of the ecosystem consequences of complex interactions, yielding a better understanding of the role of seafloor animals in affecting key processes like primary productivity. Combining predictions based on natural history, behavior of key benthic species and environmental context allow assessment of differences in functioning and process, even when the measured ecosystem property in different systems is similar. Data from three sedimentary systems in New Zealand illustrate this. Analysis of the behaviors of the infaunal ecosystem engineers in each system revealed three very different mechanisms driving ecosystem function: density and excretion, sediment turnover and surface rugosity, and hydraulic activities and porewater bioadvection. Integrative metrics of ecosystem function in some cases differentiate among the systems (gross primary production) and in others do not (photosynthetic efficiency). Analyses based on behaviors and activities revealed important ecosystem functional differences and can dramatically improve our ability to model the impact of stressors on ecosystem and global processes.

Spatial variation in the composition of communities is the product of many biotic and environmental interactions. A neglected factor in the analysis of community distribution patterns is the multi-scale nature of the data, which has implications for understanding ecological processes and the development of conservation and environmental management practice. Drawing on recently established multivariate spatial analyses, we investigate whether including relationships between spatial structure and abiotic variables enable us to better discern patterns of species and communities across scales. Data comprised 1200 macrozoobenthic samples collected over an array of distances (30 cm to 1 km) in three New Zealand harbours, as well as commonly used abiotic variables, such as sediment characteristics and chlorophyll a concentrations, measured at the same scales. Moran's eigenvector mapping was used to extract spatial scales at which communities were structured. Benthic communities, representing primarily bivalves, polychaetes and crustaceans, were spatially structured at four spatial scales, i.e. >100 m, 50-100 m, 50-15 m, and < 15 m. A broad selection of abiotic variables contributed to the large-scale variation, whereas a more limited set explained part of the fine-scale community structure. Across all scales, less than 30% of the variation in spatial structure was captured by our analysis. The large number of species (48) making up the 10 highest species scores based on redundancy analyses illustrate the variability of species-scale associations. Our results emphasise that abiotic variables and biodiversity are related at all scales investigated and stress the importance of assessing the relationship between environmental variables and the abundance and distribution of biological assemblages across a range of different scales.

Identifying the effects of stressors before they impact ecosystem functioning can be challenging in dynamic, heterogeneous 'real-world' ecosystems. In aquatic systems, for example, reductions in water clarity can limit the light available for photosynthesis, with knock-on consequences for secondary consumers, though in naturally turbid wave-swept estuaries, detecting the effects of elevated turbidity can be difficult. The objective of this study was to investigate the effects of shading on ecosystem functions mediated by sandflat primary producers (microphytobenthos) and deep-dwelling surface-feeding macrofauna (Macomona liliana; Bivalvia, Veneroida, Tellinidae). Shade cloths (which reduced incident light intensity by ~80%) were deployed on an exposed, intertidal sandflat to experimentally stress the microphytobenthic community associated with the sediment surface. After 13 weeks, sediment properties, macrofauna and fluxes of oxygen and inorganic nutrients across the sediment-water interface were measured. A multivariate metric of ecosystem function (MF) was generated by combining flux-based response variables, and distance-based linear models were used to determine shifts in the drivers of ecosystem function between non-shaded and shaded plots. No significant differences in MF or in the constituent ecosystem function variables were detected between the shaded and non-shaded plots. However, shading reduced the total explained variation in MF (from 64% in non-shaded plots to 15% in shaded plots) and affected the relative influence of M. liliana and other explanatory variables on MF. This suggests that although shade stress may shift the drivers of ecosystem functioning (consistent with earlier investigations of shading effects on sandflat interaction networks), ecosystem functions appear to have a degree of resilience to those changes.

Earth is in the midst of a biodiversity crisis that is impacting the functioning of ecosystems and the delivery of valued goods and services. However, the implications of large scale species losses are often inferred from small scale ecosystem functioning experiments with little knowledge of how the dominant drivers of functioning shift across scales. Here, by integrating observational and manipulative experimental field data, we reveal scale-dependent influences on primary productivity in shallow marine habitats, thus demonstrating the scalability of complex ecological relationships contributing to coastal marine ecosystem functioning. Positive effects of key consumers (burrowing urchins, Echinocardium cordatum) on seafloor net primary productivity (NPP) elucidated by short-term, single-site experiments persisted across multiple sites and years. Additional experimentation illustrated how these effects amplified over time, resulting in greater primary producer biomass (sediment chlorophyll a content) in the longer term, depending on climatic context and habitat factors affecting the strengths of mutually reinforcing feedbacks [corrected]. The remarkable coherence of results from small and large scales is evidence of real-world ecosystem function scalability and ecological self-organisation. This discovery provides greater insights into the range of responses to broad-scale anthropogenic stressors in naturally heterogeneous environmental settings.

The effective application of biodiversity-ecosystem function (BEF) research to societal needs amid the Anthropocene represents the next grand challenge for ecology. Biodiversity knowledge that is most meaningful to society must reconcile insights derived from theory with detailed experiments and broad-scale trends. This perspective requires science that addresses high species richness, redundancy, and natural variability, which simplified 'model systems' cannot mimic. Here, we illustrate solutions of biodiversity knowledge to management and societal problems that combine BEF with scaling experiments, analysis of BEF along environmental gradients, and mapping technologies. We primarily draw examples from biophysical interactions in seafloor environments, which cover 70% of the Earth and add significantly to global ecosystem functions and services., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Thresholds profoundly affect our understanding and management of ecosystem dynamics, but we have yet to develop practical techniques to assess the risk that thresholds will be crossed. Combining ecological knowledge of critical system interdependencies with a large-scale experiment, we tested for breaks in the ecosystem interaction network to identify threshold potential in real-world ecosystem dynamics. Our experiment with the bivalves Macomona liliana and Austrovenus stutchburyi on marine sandflats in New Zealand demonstrated that reductions in incident sunlight changed the interaction network between sediment biogeochemical fluxes, productivity, and macrofauna. By demonstrating loss of positive feedbacks and changes in the architecture of the network, we provide mechanistic evidence that stressors lead to break points in dynamics, which theory predicts predispose a system to a critical transition.

It is important to understand the consequences of low level disturbances on the functioning of ecological communities because of the pervasiveness and frequency of this type of environmental change. In this study we investigated the response of a heterogeneous, subtidal, soft-sediment habitat to small experimental additions of organic matter and calcium carbonate to examine the sensitivity of benthic ecosystem functioning to changes in sediment characteristics that relate to the environmental threats of coastal eutrophication and ocean acidification. Our results documented significant changes between key biogeochemical and sedimentary variables such as gross primary production, ammonium uptake and dissolved reactive phosphorus flux following treatment additions. Moreover, the application of treatments affected relationships between macrofauna communities, sediment characteristics (e.g., chlorophyll a content) and biogeochemical processes (oxygen and nutrient fluxes). In this experiment organic matter and calcium carbonate showed persistent opposing effects on sedimentary processes, and we demonstrated that highly heterogeneous sediment habitats can be surprisingly sensitive to subtle perturbations. Our results have important biological implications in a world with relentless anthropogenic inputs of atmospheric CO2 and nutrients in coastal waters.

Habitat-modifying organisms that impact other organisms and local functioning are important in determining ecosystem resilience. However, it is often unclear how the outcome of interactions performed by key species varies depending on the spatial and temporal disturbance context which makes the prediction of disturbance-driven regime shifts difficult. We investigated the strength and generality of effects of the filter feeding cockle Cerastoderma edule on its ambient intertidal benthic physical and biological environment. By comparing the magnitude of the effect of experimental cockle removal between a non-cohesive and a sheltered cohesive sediment in two different periods of the year, we show that the outcome of cockle interference effects relates to differences in physical disturbance, and to temporal changes in suspended sediment load and ontogenetic changes in organism traits. Interference effects were only present in the cohesive sediments, though the effects varied seasonally. Cockle presence decreased only the density of surface-dwelling species suggesting that interference effects were particularly mediated by bioturbation of the surface sediments. Furthermore, density reductions in the presence of cockles were most pronounced during the season when larvae and juveniles were present, suggesting that these life history stages are most vulnerable to interference competition. We further illustrate that cockles may enhance benthic microalgal biomass, most likely through the reduction of surface-dwelling grazing species, especially in periods with high sediment load and supposedly also high bioturbation rates. Our results emphasize that the physical disturbance of the sediment may obliterate biotic interactions, and that temporal changes in environmental stressors, such as suspended sediments, may affect the outcome of key species interference effects at the local scale. Consequently, natural processes and anthropogenic activities that change bed shear stress and sediment dynamics in coastal soft-sediment systems will affect cockle-mediated influences on ecosystem properties and therefore the resilience of these systems to environmental change.

Coastal ecosystems are often stressed by non-point source and cumulative effects that can lead to local-scale community homogenisation and a concomitant loss of large-scale ecological connectivity. Here we investigate the use of β-diversity as a measure of both community heterogeneity and ecological connectivity. To understand the consequences of different environmental scenarios on heterogeneity and connectivity, it is necessary to understand the scale at which different environmental factors affect β-diversity. We sampled macrofauna from intertidal sites in nine estuaries from New Zealand's North Island that represented different degrees of stress derived from land-use. We used multiple regression models to identify relationships between β-diversity and local sediment variables, factors related to the estuarine and catchment hydrodynamics and morphology and land-based stressors. At local scales, we found higher β-diversity at sites with a relatively high total richness. At larger scales, β-diversity was positively related to γ-diversity, suggesting that a large regional species pool was linked with large-scale heterogeneity in these systems. Local environmental heterogeneity influenced β-diversity at both local and regional scales, although variables at the estuarine and catchment scales were both needed to explain large scale connectivity. The estuaries expected a priori to be the most stressed exhibited higher variance in community dissimilarity between sites and connectivity to the estuary species pool. This suggests that connectivity and heterogeneity metrics could be used to generate early warning signals of cumulative stress.

Interaction between the diversity of local communities and the degree of connectivity between them has the potential to influence local recovery rates and thus profoundly affect community dynamics in the face of the cumulative impacts that occur across regions. Although such complex interactions have been modeled, field experiments in natural ecosystems to investigate the importance of interactions between local and regional processes are rare, especially so in coastal marine seafloor habitats subjected to many types of disturbance. We conducted a defaunation experiment at eight subtidal sites, incorporating manipulation of habitat structure, to test the relative importance of local habitat features and colonist supply in influencing macrobenthic community recovery rate. Our sites varied in community composition, habitat characteristics, and hydrodynamic conditions, and we conducted the experiment in two phases, exposing defaunated plots to colonists during periods of either high or low larval colonist supply. In both phases of the experiment, five months after disturbance, we were able to develop models that explained a large proportion of variation in community recovery rate between sites. Our results emphasize that the connectivity to the regional species pool influences recovery rate, and although local habitat effects were important, the strength of these effects was affected by broader-scale site characteristics and connectivity. Empirical evidence that cross-scale interactions are important in disturbance-recovery dynamics emphasizes the complex dynamics underlying seafloor community responses to cumulative disturbance.