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

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

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

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

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

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

Who is contributing where? Predicting ecosystem service multifunctionality for shellfish species through ecological principles

A key challenge in environmental management is determining how to manage multiple ecosystem services (ES) simultaneously, to ensure efficient and sustainable use of the environment and its resources. In marine environments, the spatial assessment of ES is lagging as a result of data-scarcity and modelling complexity. Applying mechanistic models to link ecological processes with ecosystem functions and services to assess areas of high ES potential can bridge this gap and accommodate assessments of functional differences between service providers. Here, we applied an ecosystem principles approach to assess ES potential for food provision, water quality regulation, nitrogen removal, and sediment stabilisation, provided by two estuarine bivalves (Austrovenus stutchburyi and Paphies australis) that differ in habitat association (broad and narrow distributions), to gain insight into the utility of these models for local-scale management. Maps of individual ES displayed differing patterns related to habitat associations of the species providing them, with variation in the quantities of services being delivered and locations of importance. Areas of importance for the provision of multiple services (number of services provided and their combined intensity per species) were assessed using hotspot analyses, which suggested that areas of high shellfish density at the harbour entrances were important for ES multifunctionality. A targeted management approach that includes environmental context, rather than a focus solely on the protection of high-density shellfish areas, is required to sustain the provision of individual ES.

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.

Sea Ice Dynamics Drive Benthic Microbial Communities in McMurdo Sound, Antarctica

Climate change is driving dramatic variability in sea ice dynamics, a key driver in polar marine ecosystems. Projected changes in Antarctica suggest that regional warming will force dramatic shifts in sea ice thickness and persistence, altering sea ice-associated primary production and deposition to the seafloor. To improve our understanding of the impacts of sea ice change on benthic ecosystems, we directly compared the benthic microbial communities underlying first-year sea ice (FYI) and multi-year sea ice (MYI). Using two tractable coastal habitats in McMurdo Sound, Antarctica, where FYI (Cape Evans) and MYI (New Harbour) prevail, we show that the structure and composition of the benthic microbial communities reflect the legacy of sea ice dynamics. At Cape Evans, an enrichment of known heterotrophic algal polysaccharide degrading taxa (e.g., Flavobacteriaceae, unclassified Gammaproteobacteria, and Rubritaleaceae) and sulfate-reducing bacteria (e.g., Desulfocapsaceae) correlated with comparatively higher chlorophyll a (14.2±0.8μgg^-1) and total organic carbon content (0.33%±0.04), reflecting increased productivity and seafloor deposition beneath FYI. Conversely, at New Harbour, an enrichment of known archaeal (e.g., Nitrosopumilaceae) and bacterial (e.g., Woeseiaceae and Nitrospiraceae) chemoautotrophs was common in sediments with considerably lower chlorophyll a (1.0±0.24μgg^-1) and total organic carbon content (0.17%±0.01), reflecting restricted productivity beneath MYI. We also report evidence of a submarine discharge of sub-permafrost brine from Taylor Valley into New Harbour. By comparing our two study sites, we show that under current climate-warming scenarios, changes to sea ice productivity and seafloor deposition are likely to initiate major shifts in benthic microbial communities, with heterotrophic organic matter degradation processes becoming increasingly important. This study provides the first assessment of how legacy sea ice conditions influence benthic microbial communities in Antarctica, contributing insight into sea ice-benthic coupling and ecosystem functioning in a polar environment.

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.

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

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

Invasive ecosystem engineers threaten benthic nitrogen cycling by altering native infaunal and biofouling communities

Predicting the effects of invasive ecosystem engineering species in new bioregions has proved elusive. In part this is because separating biological effects from purely physical mechanisms has been little studied and yet could help predict potentially damaging bioinvasions. Here we tested the effects of a large bio-engineering fanworm Sabella spallanzanii (Sabella) versus worm-like structures (mimics) on gas and nutrient fluxes in a marine soft bottom sediment. Experimental plots of sediment in Hauraki Gulf (New Zealand) were used to test the hypothesis that ecosystem engineers negatively influence benthic ecosystem function through autogenic mechanisms, facilitating activity by biofouling organisms and competitive exclusion of native infauna. Enhanced physical structure associated with Sabella and mimics increased nitrogen fluxes, community metabolism and reduced denitrification from 23 μmol m⁻² h⁻¹ to zero at densities greater than 25 m². Sabella plots on average had greater respiration (29%), NH₄ release (33%), and greater NO₃ release (52%) compared to mimics, suggesting allogenic (biological) mechanisms occur, but play a secondary role to autogenic (physical) mechanisms. The dominance of autogenic mechanisms indicates that bio-engineers are likely to cause significant impacts when established, regardless of fundamental differences in recipient regions or identity of the introduced bio-engineer. In the case of Sabella spallanzanii, compromised denitrification has the potential to tip the balance of net solute and gas exchanges and cause further ecological degradation in an already eutrophic system.

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

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

In situ response of Antarctic under-ice primary producers to experimentally altered pH

Elevated atmospheric CO₂ concentrations are contributing to ocean acidification (reduced seawater pH and carbonate concentrations), with potentially major ramifications for marine ecosystems and their functioning. Using a novel in situ experiment we examined impacts of reduced seawater pH on Antarctic sea ice-associated microalgal communities, key primary producers and contributors to food webs. pH levels projected for the following decades-to-end of century (7.86, 7.75, 7.61), and ambient levels (7.99), were maintained for 15 d in under-ice incubation chambers. Light, temperature and dissolved oxygen within the chambers were logged to track diurnal variation, with pH, O₂, salinity and nutrients assessed daily. Uptake of CO₂ occurred in all treatments, with pH levels significantly elevated in the two extreme treatments. At the lowest pH, despite the utilisation of CO₂ by the productive microalgae, pH did not return to ambient levels and carbonate saturation states remained low; a potential concern for organisms utilising this under-ice habitat. However, microalgal community biomass and composition were not significantly affected and only modest productivity increases were noted, suggesting subtle or slightly positive effects on under-ice algae. This in situ information enables assessment of the influence of future ocean acidification on under-ice community characteristics in a key coastal Antarctic habitat.

Sediment carbon and nutrient fluxes from cleared and intact temperate mangrove ecosystems and adjacent sandflats

The loss of mangrove ecosystems is associated with numerous impacts on coastal and estuarine function, including sediment carbon and nutrient cycling. In this study we compared in situ fluxes of carbon dioxide (CO₂) from the sediment to the atmosphere, and fluxes of dissolved inorganic nutrients and oxygen across the sediment-water interface, in intact and cleared mangrove and sandflat ecosystems in a temperate estuary. Measurements were made 20 and 25months after mangrove clearance, in summer and winter, respectively. Sediment CO₂ efflux was over two-fold higher from cleared than intact mangrove ecosystems at 20 and 25months after mangrove clearance. The higher CO₂ efflux from the cleared site was explained by an increase in respiration of dead root material along with sediment disturbance following mangrove clearance. In contrast, sediment CO₂ efflux from the sandflat site was negligible (≤9.13±1.18mmolm^-2d^-1), associated with lower sediment organic matter content. The fluxes of inorganic nutrients (NH₄⁺, NO_x and PO₄^3-) from intact and cleared mangrove sediments were low (≤20.37±18.66μmolm^-2h^-1). The highest NH₄⁺ fluxes were measured at the sandflat site (69.21±13.49μmolm^-2h^-1). Lower inorganic nutrient fluxes within the cleared and intact mangrove sites compared to the sandflat site were associated with lower abundance of larger burrowing macrofauna. Further, a higher fraction of organic matter, silt and clay content in mangrove sediments may have limited nutrient exchange.

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.

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.

Effects of Detrital Subsidies on Soft-Sediment Ecosystem Function Are Transient and Source-Dependent

Detrital subsidies from marine macrophytes are prevalent in temperate estuaries, and their role in structuring benthic macrofaunal communities is well documented, but the resulting impact on ecosystem function is not understood. We conducted a field experiment to test the effects of detrital decay on soft-sediment primary production, community metabolism and nutrient regeneration (measures of ecosystem function). Twenty four (2 m²) plots were established on an intertidal sandflat, to which we added 0 or 220 g DW m^-2 of detritus from either mangroves (Avicennia marina), seagrass (Zostera muelleri), or kelp (Ecklonia radiata) (n = 6 plots per treatment). Then, after 4, 17 and 46 d we measured ecosystem function, macrofaunal community structure and sediment properties. We hypothesized that (1) detrital decay would stimulate benthic primary production either by supplying nutrients to the benthic macrophytes, or by altering the macrofaunal community; and (2) ecosystem responses would depend on the stage and rate of macrophyte decay (a function of source). Avicennia detritus decayed the slowest with a half-life (t₅₀) of 46 d, while Zostera and Ecklonia had t₅₀ values of 28 and 2.6 d, respectively. However, ecosystem responses were not related to these differences. Instead, we found transient effects (up to 17 d) of Avicennia and Ecklonia detritus on benthic primary production, where initially (4 d) these detrital sources suppressed primary production, but after 17 d, primary production was stimulated in Avicennia plots relative to controls. Other ecosystem function response variables and the macrofaunal community composition were not altered by the addition of detritus, but did vary with time. By sampling ecosystem function temporally, we were able to capture the in situ transient effects of detrital subsidies on important benthic ecosystem functions.

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.

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.

Sedimentary environment influences the effect of an infaunal suspension feeding bivalve on estuarine ecosystem function

The suspension feeding bivalve Austrovenus stutchburyi is a key species on intertidal sandflats in New Zealand, affecting the appearance and functioning of these systems, but is susceptible to several environmental stressors including sedimentation. Previous studies into the effect of this species on ecosystem function have been restricted in space and time, limiting our ability to infer the effect of habitat change on functioning. We examined the effect of Austrovenus on benthic primary production and nutrient dynamics at two sites, one sandy, the other composed of muddy-sand to determine whether sedimentary environment alters this key species' role. At each site we established large (16 m²) plots of two types, Austrovenus addition and removal. In winter and summer we deployed light and dark benthic chambers to quantify oxygen and nutrient fluxes and measured sediment denitrification enzyme activity to assess denitrification potential. Rates of gross primary production (GPP) and ammonium uptake were significantly increased when Austrovenus was added, relative to removed, at the sandy site (GPP, 1.5 times greater in winter and summer; ammonium uptake, 8 times greater in summer; 3-factor analysis of variance (ANOVA), p<0.05). Denitrification potential was also elevated in Austrovenus addition plots at the sandy site in summer (by 1.6 times, p<0.1). In contrast, there was no effect of Austrovenus treatment on any of these variables at the muddy-sand site, and overall rates tended to be lower at the muddy-sand site, relative to the sandy site (e.g. GPP was 2.1 to 3.4 times lower in winter and summer, respectively, p<0.001). Our results suggest that the positive effects of Austrovenus on system productivity and denitrification potential is limited at a muddy-sand site compared to a sandy site, and reveal the importance of considering sedimentary environment when examining the effect of key species on ecosystem function.

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.

Forecasting the limits of resilience: integrating empirical research with theory

Despite the increasing evidence of drastic and profound changes in many ecosystems, often referred to as regime shifts, we have little ability to understand the processes that provide insurance against such change (resilience). Modelling studies have suggested that increased variance may foreshadow a regime shift, but this requires long-term data and knowledge of the functional links between key processes. Field-based research and ground-truthing is an essential part of the heuristic that marries theoretical and empirical research, but experimental studies of resilience are lagging behind theory, management and policy requirements. Empirically, ecological resilience must be understood in terms of community dynamics and the potential for small shifts in environmental forcing to break the feedbacks that support resilience. Here, we integrate recent theory and empirical data to identify ways we might define and understand potential thresholds in the resilience of nature, and thus the potential for regime shifts, by focusing on the roles of strong and weak interactions, linkages in meta-communities, and positive feedbacks between these and environmental drivers. The challenge to theoretical and field ecologists is to make the shift from hindsight to a more predictive science that is able to assist in the implementation of ecosystem-based management.

The effects of habitat loss, fragmentation, and community homogenization on resilience in estuaries

When changes in the frequency and extent of disturbance outstrip the recovery potential of resident communities, the selective removal of species contributes to habitat loss and fragmentation across landscapes. The degree to which habitat change is likely to influence community resilience will depend on metacommunity structure and connectivity. Thus ecological connectivity is central to understanding the potential for cumulative effects to impact upon diversity. The importance of these issues to coastal marine communities, where the prevailing concept of open communities composed of highly dispersive species is being challenged, indicates that these systems may be more sensitive to cumulative impacts than previously thought. We conducted a disturbance-recovery experiment across gradients of community type and environmental conditions to assess the roles of ecological connectivity and regional variations in community structure on the recovery of species richness, total abundance, and community composition in Mahurangi Harbour, New Zealand. After 394 days, significant differences in recovery between sites were apparent. Statistical models explaining a high proportion of the variability (R2 > 0.92) suggested that community recovery rates were controlled by a combination of physical and ecological features operating across spatial scales, affecting successional processes. The dynamic and complex interplay of ecological and environmental processes we observed driving patch recovery across the estuarine landscape are integral to recovery from disturbances in heterogeneous environments. This link between succession/recovery, disturbance, and heterogeneity confirms the utility of disturbance-recovery experiments as assays for cumulative change due to fragmentation and habitat change in estuaries.

Bioturbators enhance ecosystem function through complex biogeochemical interactions

Predicting the consequences of species loss is critically important, given present threats to biological diversity such as habitat destruction, overharvesting and climate change. Several empirical studies have reported decreased ecosystem performance (for example, primary productivity) coincident with decreased biodiversity, although the relative influence of biotic effects and confounding abiotic factors has been vigorously debated. Whereas several investigations focused on single trophic levels (for example, grassland plants), studies of whole systems have revealed multiple layers of feedbacks, hidden drivers and emergent properties, making the consequences of species loss more difficult to predict. Here we report functionally important organisms and considerable biocomplexity in a sedimentary seafloor habitat, one of Earth's most widespread ecosystems. Experimental field measurements demonstrate how the abundance of spatangoid urchins--infaunal (in seafloor sediment) grazers/deposit feeders--is positively related to primary production, as their activities change nutrient fluxes and improve conditions for production by microphytobenthos (sedimentatry microbes and unicellular algae). Declines of spatangoid urchins after trawling are well documented, and our research linking these bioturbators to important benthic-pelagic fluxes highlights potential ramifications for productivity in coastal oceans.

Dredging-induced nutrient release from sediments to the water column in a southeastern saltmarsh tidal creek

Dredging is a large-scale anthropogenic disturbance agent in coastal and estuarine habitats that can profoundly affect water quality. We examined the impact of a small-scale dredging operation in a salt marsh in South Carolina by comparing nutrient levels (NH(4)(+), NO(x), PO(4)(-)) and total suspended solid concentrations before and during dredging activities. Nutrient enrichment was evaluated within the context of tidal, seasonal, and inter-annual variability by using long-term water chemistry data provided by the North Inlet-Winyah Bay National Estuarine Research Reserve. The conditions of the dredging permit (i.e., its relatively small scale), the season chosen for the work (fall-winter), the nature of the sediments dredged (coarse-grained), and the amount of natural variability in the estuary's water chemistry (even on a daily time-scale) all minimized the impact of the dredging activities. Results of this study will add to the limited body of empirical data that should be considered in evaluating future dredging permit applications related to shallow estuarine waterways.