Polyester microfiber impacts on coastal sediment organic matter consumption.

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.

Toward a network perspective in coastal ecosystem management.

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.

Ecological network analysis of traits reveals variable response capacity to stress.

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.

Real-world impacts of microplastic pollution on seafloor ecosystem function.

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.

The influence of mussel restoration on coastal carbon cycling.

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.

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

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

Enhancing multiple scales of seafloor biodiversity with mussel restoration.

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.

Positive contribution of macrofaunal biodiversity to secondary production and seagrass carbon metabolism.

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.

Beyond the single index: Investigating ecological mechanisms underpinning ecosystem multifunctionality with network analysis.

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.

Climate cascades affect coastal Antarctic seafloor ecosystem functioning.

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.

The distribution and ecological effects of microplastics in an estuarine ecosystem.

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.

Identifying "vital attributes" for assessing disturbance-recovery potential of seafloor communities.

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.

Sampling frequency, duration and the Southern Oscillation influence the ability of long-term studies to detect sudden change.

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.

Predicting habitat suitability of filter-feeder communities in a shallow marine environment, New Zealand.

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/Totaranui 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.

A call to evaluate Plastic's impacts on marine benthic ecosystem interaction networks.

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.

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

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

The impacts of polyethylene terephthalate microplastics (mPETs) on ecosystem functionality in marine sediment.

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.

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.

Co-occurrence patterns and the large-scale spatial structure of benthic communities in seagrass meadows and bare sand.

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 m2 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.

Effects of Polyester Microfibers on Microphytobenthos and Sediment-Dwelling Infauna.

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 (NOx) was significantly elevated and related to M. liliana and MPB biochemical properties. Ammoniu (NH4+) 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.