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

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.

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 their potential to mitigate climate change as blue carbon habitats have been repeatedly highlighted. However, little 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 while 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 likely becomes 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. Hence, 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 worlds' coasts. This article is protected by copyright. All rights reserved.

Informing the management of multiple stressors on estuarine ecosystems using an expert-based Bayesian Network model

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.

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.

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.

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.

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.

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

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.

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 m² 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 (NO_x) was significantly elevated and related to M. liliana and MPB biochemical properties. Ammoniu (NH₄⁺) 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.

Investigating changes in estuarine ecosystem functioning under future scenarios

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.

Linking changes in species-trait relationships and ecosystem function using a network analysis of traits

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² ], medium [150 g N/m² ], and high [600 g N/m² ]). 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² ) N treatment, functional diversity remained consistent, whereas increasing disturbance to high (600 g N/m² ) N treatment affected the species-trait network structure and caused merging of functional clusters implying a loss of functional trait diversity.

Monitoring for tipping points in the marine environment

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.

Benthic responses to an Antarctic regime shift: food particle size and recruitment biology

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.

The multiple roles of β-diversity help untangle community assembly processes affecting recovery of temperate rocky shores

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.

Long-term environmental monitoring for assessment of change: measurement inconsistencies over time and potential solutions

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.

Changes in the location of biodiversity-ecosystem function hot spots across the seafloor landscape with increasing sediment nutrient loading

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

In situ soft sediment nutrient enrichment: A unified approach to eutrophication field experiments

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.

The role of time and species identities in spatial patterns of species richness and conservation

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.

Optimization of hard clams, polychaetes, physical disturbance and denitrifying bacteria of removing nutrients in marine sediment

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.

Multiple stressors, nonlinear effects and the implications of climate change impacts on marine coastal ecosystems

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.

Same pattern, different mechanism: Locking onto the role of key species in seafloor ecosystem process

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.

Detecting Subtle Shifts in Ecosystem Functioning in a Dynamic Estuarine Environment

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.

Experimenting with ecosystem interaction networks in search of threshold potentials in real-world marine ecosystems

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.

Sensitivity of heterogeneous marine benthic habitats to subtle stressors

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 CO₂ and nutrients in coastal waters.

Counting on β-diversity to safeguard the resilience of estuaries

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.

When small changes matter: the role of cross-scale interactions between habitat and ecological connectivity in recovery

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.

Detecting shifts in ecosystem functioning: the decoupling of fundamental relationships with increased pollutant stress on sandflats

In this study, we investigated the influence of low level contamination by copper, lead and zinc on the functioning of estuarine sandflat ecosystems by comparing the strength and variability of relationships between benthic macrofauna and fluxes (oxygen and nutrients) at three clean and three mildly contaminated sites. Specifically, as indicators of ecosystem functioning, we examined relationships between bivalve biomass, total benthic respiration and ammonium release, and ammonium uptake and benthic primary production. Furthermore, a small amount of organic matter was added to experimental plots at all sites (35 g/0.2 m²) to evaluate stress-on-stress responses relative to controls. Relationships were strongest at the clean sites (steepest slopes, highest r² values, lowest p-values) and weakest at the mildly contaminated sites and in organically enriched plots. Our results suggest that changes in ecosystem functioning may be occurring at mild (<ERL) levels of pollution that are unlikely to cause widespread mortality of macrofauna.

Ecosystem services transcend boundaries: estuaries provide resource subsidies and influence functional diversity in coastal benthic communities

BACKGROUND

Estuaries are highly productive ecosystems that can export organic matter to coastal seas (the 'outwelling hypothesis'). However the role of this food resource subsidy on coastal ecosystem functioning has not been examined.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated the influence of estuarine primary production as a resource subsidy and the influence of estuaries on biodiversity and ecosystem functioning in coastal mollusk-dominated sediment communities. Stable isotope values (δ(13)C, δ(15)N) demonstrated that estuarine primary production was exported to the adjacent coast and contributed to secondary production up to 4 km from the estuary mouth. Further, isotope signatures of suspension feeding bivalves on the adjacent coast (Dosinia subrosea) closely mirrored the isotope values of the dominant bivalves inside the estuaries (Austrovenus stutchburyi), indicating utilization of similar organic matter sources. However, the food subsidies varied between estuaries; with estuarine suspended particulate organic matter (SPOM) dominant at Tairua estuary, while seagrass and fringing vegetation detritus was proportionately more important at Whangapoua estuary, with lesser contributions of estuarine SPOM. Distance from the estuary mouth and the size and density of large bivalves (Dosinia spp.) had a significant influence on the composition of biological traits in the coastal macrobenthic communities, signaling the potential influence of these spatial subsidies on ecosystem functioning.

CONCLUSIONS/SIGNIFICANCE

Our study demonstrated that the locations where ecosystem services like productivity are generated are not necessarily where the services are utilized. Further, we identified indirect positive effects of the nutrient subsidies on biodiversity (the estuarine subsidies influenced the bivalves, which in turn affected the diversity and functional trait composition of the coastal sediment macrofaunal communities). These findings highlight the importance of integrative ecosystem-based management that maintains the connectivity of estuarine and coastal ecosystems.

Interaction networks in coastal soft-sediments highlight the potential for change in ecological resilience

Recent studies emphasize the role of indirect relationships and feedback loops in maintaining ecosystem resilience. Environmental changes that impact on the organisms involved in these processes have the potential to initiate threshold responses and fundamentally shift the interactions within an ecosystem. However, empirical studies are hindered by the difficulty of designing appropriate manipulative experiments to capture this complexity. Here we employ structural equation modeling to define and test the architecture of ecosystem interaction networks. Using survey data from 19 estuaries we investigate the interactions between biological (abundance of large bioturbating macrofauna, microphytobenthos, and detrital matter) and physical (sediment grain size) processes. We assess the potential for abrupt changes in the architecture of the network and the strength of interactions to occur across environmental gradients. Our analysis identified a potential threshold in the relationship between sediment mud content and benthic chlorophyll a, at -12 microg/g, using quantile regression. Below this threshold, the interaction network involved different variables and fewer feedbacks than above. This approach has potential to improve our empirical understanding of thresholds in ecological systems and our ability to design manipulative experiments that test how and when a threshold will be passed. It can also be used to indicate to resource managers that a particular system has the potential to exhibit threshold responses to environmental change, emphasizing precautionary management and facilitating a better understanding of how persistent multiple stressors threaten the resilience and long-term use of natural ecosystems.

Disturbance of sandflats by thin terrigenous sediment deposits: consequences for primary production and nutrient cycling

The contributions of biodiversity to ecosystem functioning are increasingly recognized by ecologists, with biodiversity loss considered a significant threat to the maintenance of life-supporting ecosystem goods and services. Although ecologists have increased the amount of realism in biodiversity-ecosystem functioning (BEF) experiments, effects on functioning are rarely investigated in the field in conjunction with disturbances that affect biodiversity. Here, effects on functioning were investigated in situ following experimental disturbance of a realistic type and magnitude. Experimental deposits of terrigenous sediment (5 mm thickness) were applied at three intertidal sites in Whangapoua Harbour (Coromandel Peninsula, New Zealand), where soil erosion associated with logging is a serious problem. Experimental disturbance by terrigenous sediment impacted macrofaunal abundance and community composition (suspension and deposit feeders), and there were coincident shifts in ecosystem functioning at all three sites. Relative to controls, sediment-treated plots had higher rates of ammonium efflux (despite 50% fewer macrofaunal individuals) and lower rates of gross primary production and photosynthetic efficiency (despite similar sediment chlorophyll a concentrations). Judging from nutrient ratios in sediment pore water, microphyte primary producers living on the sediment surface in control plots were likely nitrogen limited (the normal situation for marine waters), whereas microphytes in sediment-treated plots were likely phosphate limited. Gross primary production and photosynthetic efficiency were significantly and positively correlated with ammonium uptake in control plots (r2 = 0.463 and 0.856, respectively) but not in treated plots (P > 0.05, r2 < 0.3). We suggest that the higher content of charged particles (clay) in the experimental deposits bound up phosphate, limiting its bio-availability, and shifted functional relationships between sandflat macrofauna and sediment primary productivity.

beta-diversity and species accumulation in antarctic coastal benthos: influence of habitat, distance and productivity on ecological connectivity

High Antarctic coastal marine environments are comparatively pristine with strong environmental gradients, which make them important places to investigate biodiversity relationships. Defining how different environmental features contribute to shifts in beta-diversity is especially important as these shifts reflect both spatio-temporal variations in species richness and the degree of ecological separation between local and regional species pools. We used complementary techniques (species accumulation models, multivariate variance partitioning and generalized linear models) to assess how the roles of productivity, bio-physical habitat heterogeneity and connectivity change with spatial scales from metres to 100's of km. Our results demonstrated that the relative importance of specific processes influencing species accumulation and beta-diversity changed with increasing spatial scale, and that patterns were never driven by only one factor. Bio-physical habitat heterogeneity had a strong influence on beta-diversity at scales <290 km, while the effects of productivity were low and significant only at scales >40 km. Our analysis supports the emphasis on the analysis of diversity relationships across multiple spatial scales and highlights the unequal connectivity of individual sites to the regional species pool. This has important implications for resilience to habitat loss and community homogenisation, especially for Antarctic benthic communities where rates of recovery from disturbance are slow, there is a high ratio of poor-dispersing and brooding species, and high biogenic habitat heterogeneity and spatio-temporal variability in primary production make the system vulnerable to disturbance. Consequently, large areas need to be included within marine protected areas for effective management and conservation of these special ecosystems in the face of increasing anthropogenic disturbance.