Search

To investigate regional drivers of spatial patterns in macro- and meio-faunal community structure (abundance, biomass and taxonomic diversity) and ecosystem function (sediment community oxygen consumption [SCOC]), we sampled two regions in close proximity on New Zealand's continental margin—the Chatham Rise and the Challenger Plateau. Sites (n= 15) were selected in water depths ranging from 266–1212 m to generate a gradient in sedimentary properties and, in particular, surface pelagic productivity. Both macro- and meio-fauna abundance and biomass was 2–3.5 times higher on the Chatham Rise than on the Challenger Plateau, reflecting regional differences in pelagic primary production. We also found significant inter- and intra-regional differences in macro-fauna taxonomic diversity with two distinctive site groupings in each region. Univariate and multivariate measures of macro-fauna community attributes were most strongly correlated with sediment photosynthetic pigment (explaining 24%–59% of the variation). Sediment pigment content was as equally important in explaining meio-fauna community structure (36%–7%). Unlike community structure, SCOC was most strongly correlated with depth (44%), most likely reflecting temperature effects on benthic metabolism. Our results highlight the importance of a benthic labile food supply in structuring infaunal communities on continental margins and emphasise a tight coupling between pelagic and benthic habitats. [ABSTRACT FROM PUBLISHER]

Mangrove forests can provide important cross-boundary subsidies of organic matter to adjacent habitats through the production, export, decomposition and assimilation of litter. We quantified two of these components in a temperate mangrove forest in Whangamata Harbour, New Zealand: 1) litter production; and 2) decomposition rates as a function of tidal elevation, sediment type and burial depth. Litter traps sampled monthly for a year measured an annual detrital input of 3.24–5.38 t DW ha−1, of which 77% occurred in summer. Decomposition rates depended on litter type, with leaves decomposing faster (63 d to decay by 50%) than pneumatophore and wood material (316 and 460 d, respectively). Buried leaf and wood litter decomposed 1.3–1.4 times slower than litter on the sediment surface; however, tidal elevation and sediment type (mud vs. sand) had no effect. The slow decay of litter (an order of magnitude slower than tropical mangrove litter) suggests that its incorporation into temperate marine food webs may be relatively slow. [ABSTRACT FROM AUTHOR]

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

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

New Zealand has a complex coastal environment spanning a large latitudinal gradient and three water masses. Here we assess whether multivariate analyses of benthic macrofaunal community composition can be a sensitive approach to assessing relative estuarine health across the country, negating the need for regional indices and reducing reliance on reference sites. Community data were used in separate canonical analyses of principal coordinates to create multivariate models of community responses to gradients in mud content and heavy metal contamination. Both models performed well (R 2  = 0.81, 0.71), and were unaffected by regional and estuarine typology differences. The models demonstrate a sensitive and standardized approach to assessing estuarine health that allowed separation of the two stressors. This approach could be applied to other stressors, countries or regions., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

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.

Anthropogenically-modulated reductions in pH, termed ocean acidification, could pose a major threat to the physiological performance, stocks, and biodiversity of calcifiers and may devalue their ecosystem services. Recent debate has focussed on the need to develop approaches to arrest the potential negative impacts of ocean acidification on ecosystems dominated by calcareous organisms. In this study, we demonstrate the role of a discrete (i.e. diffusion) boundary layer (DBL), formed at the surface of some calcifying species under slow flows, in buffering them from the corrosive effects of low pH seawater. The coralline macroalga Arthrocardia corymbosa was grown in a multifactorial experiment with two mean pH levels (8.05 'ambient' and 7.65 a worst case 'ocean acidification' scenario projected for 2100), each with two levels of seawater flow (fast and slow, i.e. DBL thin or thick). Coralline algae grown under slow flows with thick DBLs (i.e., unstirred with regular replenishment of seawater to their surface) maintained net growth and calcification at pH 7.65 whereas those in higher flows with thin DBLs had net dissolution. Growth under ambient seawater pH (8.05) was not significantly different in thin and thick DBL treatments. No other measured diagnostic (recruit sizes and numbers, photosynthetic metrics, %C, %N, %MgCO3) responded to the effects of reduced seawater pH. Thus, flow conditions that promote the formation of thick DBLs, may enhance the subsistence of calcifiers by creating localised hydrodynamic conditions where metabolic activity ameliorates the negative impacts of ocean acidification.

Coastal ecosystems that are characterized by kelp forests encounter daily pH fluctuations, driven by photosynthesis and respiration, which are larger than pH changes owing to ocean acidification (OA) projected for surface ocean waters by 2100. We investigated whether mimicry of biologically mediated diurnal shifts in pH-based for the first time on pH time-series measurements within a kelp forest-would offset or amplify the negative effects of OA on calcifiers. In a 40-day laboratory experiment, the calcifying coralline macroalga, Arthrocardia corymbosa, was exposed to two mean pH treatments (8.05 or 7.65). For each mean, two experimental pH manipulations were applied. In one treatment, pH was held constant. In the second treatment, pH was manipulated around the mean (as a step-function), 0.4 pH units higher during daylight and 0.4 units lower during darkness to approximate diurnal fluctuations in a kelp forest. In all cases, growth rates were lower at a reduced mean pH, and fluctuations in pH acted additively to further reduce growth. Photosynthesis, recruitment and elemental composition did not change with pH, but δ(13)C increased at lower mean pH. Including environmental heterogeneity in experimental design will assist with a more accurate assessment of the responses of calcifiers to OA.

The relationship between species diversity and environmental parameters is poorly understood for the mobile macrofauna of deep-sea habitats due to under-sampling and subsequent lack of accurate taxonomic information. To redress this, cytochrome oxidase c subunit I (COI) DNA sequences were used to estimate species diversity and to compare phoxocephalid amphipod assemblages among 20 stations encompassing a range of environmental conditions. Two regions, east (Chatham Rise) and west (Challenger Plateau) of New Zealand were sampled to depths of 200-1200 m with an epibenthic sled. Using a comparison among identified morphospecies, we found a clear gap in sequence divergences between 6% and 13% and used a 6% threshold to designate molecular operational taxonomic units (MOTUs), as a surrogate to putative species. DNA sequences (n = 297) revealed high total diversity (n = 49 MOTUs), as well as high beta diversity (28 MOTUs found at single location only). Novel phoxocephalid MOTUs were found at most stations, especially on Challenger Plateau and the flanks of Chatham Rise. Analyses of interstation assemblages revealed a major split between regions, indicating minimal overlap in taxon distributions. A cluster of highly similar stations was identified, broadly distributed over the crest of Chatham Rise, in association with elevated food availability, probably resulting from higher surface productivity and relatively shallow depth. Accordingly, multivariate analysis revealed a strong correlation between phoxocephalid assemblages and food supply. This study highlights the value of molecular approaches, in particular COI sequences, for quantifying and comparing diversity in under-sampled and/or under-studied taxa., (© 2012 Blackwell Publishing Ltd.)

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(2)) 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.