Interactions between multiple large macrofauna species and nematode communities — Mechanisms for indirect impacts of trawling disturbance

Highlights •We exposed meiofauna to 7 different large macrofauna species at high and low densities. •Macrofauna presence altered nematode community structure and reduced their abundance. •Macrofauna species had similar effects by reducing the few dominant nematode species. •Meio–macrofauna resource competition and spatial segregation are the main drivers. •Trawling effects on macrofauna affect nematode communities indirectly. Diverse assemblages of infauna in sediments provide important physical and biogeochemical services, but are under increasing pressure by anthropogenic activities, such as benthic trawling. It is known that trawling disturbance has a substantial effect on the larger benthic fauna, with reductions in density and diversity, and changes in community structure, benthic biomass, production, and bioturbation and biogeochemical processes. Largely unknown, however, are the mechanisms by which the trawling impacts on the large benthic macro- and megafauna may influence the smaller meiofauna. To investigate this, a mesocosm experiment was conducted whereby benthic nematode communities from a non-trawled area were exposed to three different densities (absent, low, normal) of 7 large (> 10 mm) naturally co-occurring, bioturbating species which are potentially vulnerable to trawling disturbance. The results showed that total abundances of nematodes were lower if these large macrofauna species were present, but no clear nematode abundance effects could be assigned to the macrofauna density differences. Nematode community structure changed in response to macrofauna presence and density, mainly as a result of the reduced abundance of a few dominant nematode species. Any detectable effects seemed similar for nearly all macrofauna species treatments, supporting the idea that there may be a general indirect, macrofauna-mediated trawling impact on nematode communities. Explanations for these results may be, firstly, competition for food resources, resulting in spatial segregation of the meio- and macrobenthic components. Secondly, different densities of large macrofauna organisms may affect the nematode community structure through different intensities of bioturbatory disturbance or resource competition. These results suggest that removal or reduced densities of larger macrofauna species as a result of trawling disturbance may lead to increased nematode abundance and hints at the validity of interference competition between large macrofauna organisms and the smaller meiofauna, and the energy equivalence hypothesis, where a trade-off is observed between groups of organisms that are dependent on a common source of energy.

Variations in GDGT distributions through the water column in the South East Atlantic Ocean

The TetraEther indeX of 86 carbon atoms (TEX86) temperature proxy is widely used in reconstructions of past sea surface temperature. Most current calibrations are based on surface sediment distributions of the glycerol dialkyl glycerol tetraether lipids (GDGTs) that comprise TEX86 and assume that these GDGTs are exported from the upper mixed layer. However, GDGT export from deeper waters could impact sedimentary GDGT distributions and therefore TEX86 paleothermometry. Here we examine GDGT distributions in suspended particulate matter (SPM) and underlying sediments collected from the Southeast Atlantic Ocean. Our results reveal different GDGT distributions - specifically the ratio between GDGTs bearing 2 vs. 3 cyclopentyl moieties, [2/3] ratios - between surface, subsurface (>50-200 m) and deep water (>200 m) SPM, which suggests the occurrence of in situ (deep) production that is not apparent when considering TEX86. The GDGT distributions in sediments match those of subsurface waters rather than surface waters, suggesting that they have not been preferentially derived from the upper mixed layer; this is consistent with GDGT abundances being highest in shallow subsurface SPM (˜100 to 200 m). It remains unclear what governs the different [2/3] ratios throughout the water column, but it is likely related to a combination of temperature and thaumarchaeotal community structure.

Effect of phytoplankton size classes on bio-optical properties of phytoplankton in the Western Iberian coast: Application of models

Chlorophyll-a satellite products are routinely used in oceanography, providing a synoptic and global view of phytoplankton abundance. However, these products lack information on the community structure of the phytoplankton, which is crucial for ecological modelling and ecosystem studies. To assess the usefulness of existing methods to differentiate phytoplankton functional types (PFT) or phytoplankton size classes from satellite data, in-situ phytoplankton samples collected in the Western Iberian coast, on the North-East Atlantic, were analysed for pigments and absorption spectra. Water samples were collected in five different locations, four of which were located near the shore and another in an open-ocean, seamount region. Three different modelling approaches for deriving phytoplankton size classes were applied to the in situ data. Approaches tested provide phytoplankton size class information based on the input of pigments data (Brewin et al., 2010), absorption spectra data (Ciotti et al., 2002) or both (Uitz et al., 2008). Following Uitz et al. (2008), results revealed high variability in microphytoplankton chlorophyll-specific absorption coefficients, ranging from 0.01 to 0.09 m2 (mg chl)− 1 between 400 and 500 nm. This spectral analysis suggested, in one of the regions, the existence of small cells (< 20 μm) in the fraction of phytoplankton presumed to be microphytoplankton (based on diagnostic pigments). Ciotti et al. (2002) approach yielded the highest differences between modelled and measured absorption spectra for the locations where samples had high variability in community structure and cell size. The Brewin et al. (2010) pigment-based model was adjusted and a set of model coefficients are presented and recommended for future studies in offshore water of the Western Iberian coast.

Volcanic ash supply to the surface ocean – remote sensing of biological responses and their wider biogeochemical significance

Transient micronutrient enrichment of the surface ocean can enhance phytoplankton growth rates and alter microbial community structure with an ensuing spectrum of biogeochemical feedbacks. Strong phytoplankton responses to micronutrients supplied by volcanic ash have been reported recently. Here we: (i) synthesize findings from these recent studies; (ii) report the results of a new remote sensing study of ash fertilization; and (iii) calculate theoretical bounds of ash-fertilized carbon export. Our synthesis highlights that phytoplankton responses to ash do not always simply mimic that of iron amendment; the exact mechanisms for this are likely biogeochemically important but are not yet well understood. Inherent optical properties of ash-loaded seawater suggest rhyolitic ash biases routine satellite chlorophyll-a estimation upwards by more than an order of magnitude for waters with <0.1 mg chlorophyll-a m-3, and less than a factor of 2 for systems with >0.5 mg chlorophyll-a m-3. For this reason post-ash-deposition chlorophyll-a changes in oligotrophic waters detected via standard Case 1 (open ocean) algorithms should be interpreted with caution. Remote sensing analysis of historic events with a bias less than a factor of 2 provided limited stand-alone evidence for ash-fertilization. Confounding factors were poor coverage, incoherent ash dispersal, and ambiguity ascribing biomass changes to ash supply over other potential drivers. Using current estimates of iron release and carbon export efficiencies, uncertainty bounds of ash-fertilized carbon export for 3 events are presented. Patagonian iron supply to the Southern Ocean from volcanic eruptions is less than that of windblown dust on thousand year timescales but can dominate supply at shorter timescales. Reducing uncertainties in remote sensing of phytoplankton response and nutrient release from ash are avenues for enabling assessment of the oceanic response to large-scale transient nutrient enrichment.

Spatio-temporal variability in ammonia oxidation and ammonia oxidising bacteria and archaea in coastal sediments of the Western English Channel

The abundance of ammonia-oxidising bacterial (AOB) and ammonia-oxidising archaeal (AOA) (amoA) genes and ammonia oxidation rates were compared bimonthly from July 2008 to May 2011 in 4 contrasting coastal sediments in the western English Channel. Despite a higher abundance of AOA amoA genes within all sediments and at all time-points, rates of ammonia oxidation correlated with AOB and not AOA amoA gene abundance. Sediment type was a major factor in determining both AOB amoA gene abundance and AOB community structure, possibly due to deeper oxygen penetration into the sandier sediments, increasing the area available for ammonia oxidation. Decreases in AOB amoA gene abundance were evident during summer and autumn, with maximum abundance and ammonia oxidation rates occurring in winter and early spring. PCR-DGGE of AOB amoA genes indicated that no seasonal changes to community composition occurred; however, a gradual movement in community composition occurred at 3 of the sites studied. The lack of correlation between AOA amoA gene abundance and ammonium oxidation rates, or any other environmental variable measured, may be related to the higher spatial variation amongst measurements, obscuring temporal trends, or the bimonthly sampling, which may have been too infrequent to capture temporal variability in the deposition of fresh organic matter. Alternatively, AOA may respond to changing substrate concentrations by an increase or decrease in transcript rather than gene abundance.

Disturbance to conserved bacterial communities in the cold water gorgonian coral Eunicella verrucosa

The bacterial communities associated with healthy and diseased colonies of the cold-water gorgonian coral Eunicella verrucosa at three sites off the south-west coast of England were compared using denaturing gradient gel electrophoresis (DGGE) and clone libraries. Significant differences in community structure between healthy and diseased samples were discovered, as were differences in the level of disturbance to these communities at each site; this correlated with depth and sediment load. The majority of cloned sequences from healthy coral tissue affiliated with the Gammaproteobacteria. The stability of the bacterial community and dominance of specific genera found across visibly healthy colonies suggest the presence of a specific microbial community. Affiliations included a high proportion of Endozoicomonas sequences, which were most similar to sequences found in tropical corals. This genus has been found in a number of invertebrates and is suggested to have a role in coral health and in the metabolisation of dimethylsulfoniopropionate (DMSP) produced by zooxanthellae. However, screening of colonies for the presence of zooxanthellae produced a negative result. Diseased colonies showed a decrease in affiliated clones and an increase in clones related to potentially harmful/transient microorganisms but no increase in a particular pathogen. This study demonstrates that a better understanding of these bacterial communities, the factors that affect them and their role in coral health and disease will be of critical importance in predicting future threats to temperate gorgonian communities.

Trait-based representation of diatom functional diversity in a plankton functional type model of the eutrophied southern North Sea

We introduce a trait-based description of diatom functional diversity to an existing plankton functional type (PFT) model, implemented for the eutrophied coastal ecosystem in the Southern Bight of the North Sea. The trait-based description represents a continuum of diatom species, each characterized by a distinct cell volume, and includes size dependence of four diatom traits: the maximum growth rate, the half-saturation constants for nutrient uptake, the photosynthetic efficiency, and the relative affinity of copepods for diatoms. Through competition under seasonally varying forcing, the fitness of each diatom varies throughout time, and the outcome of competition results in a changing community structure. The predicted seasonal change in mean cell volume of the community is supported by field observations: smaller diatoms, which are more competitive in terms of resource acquisition, prevail during the first spring bloom, whereas the summer bloom is dominated by larger species which better resist grazing. The size-based model is used to determine the ecological niche of diatoms in the area and identifies a range of viable sizes that matches observations. The general trade-off between small, competitive diatoms and large, grazing-resistant species is a convenient framework to study patterns in diatom functional diversity. PFT models and trait-based approaches constitute promising complementary tools to study community structure in marine ecosystems.

There are no whole truths in meta-analyses: all their truths are half truths

In a recent letter, Thomsen & Wernberg (2015) rean-alyzed data compiled for our recent paper (Lyonset al., 2014). In that paper, we examined the effectsof macroalgal blooms and macroalgal mats on sevenimportant measures of community structure and eco-system functioning and explored several ecologicaland methodological factors that might explain someof the variation in the observed effects. Thomsen &Wernberg (2015) re-analyzed two small subsets of the data, focusing on experimental studies examining effects of blooms/mats on invertebrate abundance.Their analyses revealed two interesting patterns.First, they showed that macroalgal blooms reducedthe abundance of communities that Thomsen andWernberg categorized as ‘mainly infauna’, whileincreasing the abundance of communities categorized as ‘mainly epifauna’. Second, they showed that theimpacts of macroalgal blooms on ‘mainly infauna’communities increased with algal density in experiments that included multiple levels of algal density.These findings, as well as the conclusions that Thomsen & Wernberg (2015) draw from them, are largely consistent with our own expectations and interpretations. However, we also feel that some caution is required when interpreting the results of their analyses.

The assessment of a global marine ecosystem model on the basis of emergent properties and ecosystem function: a case study with ERSEM

Ecosystem models are often assessed using quantitative metrics of absolute ecosystem state, but these model-data comparisons are disproportionately vulnerable to discrepancies in the location of important circulation features. An alternative method is to demonstrate the models capacity to represent ecosystem function; the emergence of a coherent natural relationship in a simulation indicates that the model may have an appropriate representation of the ecosystem functions that lead to the emergent relationship. Furthermore, as emergent properties are large-scale properties of the system, model validation with emergent properties is possible even when there is very little or no appropriate data for the region under study, or when the hydrodynamic component of the model differs significantly from that observed in nature at the same location and time. A selection of published meta-analyses are used to establish the validity of a complex marine ecosystem model and to demonstrate the power of validation with emergent properties. These relationships include the phytoplankton community structure, the ratio of carbon to chlorophyll in phytoplankton and particulate organic matter, the ratio of particulate organic carbon to particulate organic nitrogen and the stoichiometric balance of the ecosystem. These metrics can also inform aspects of the marine ecosystem model not available from traditional quantitative and qualitative methods. For instance, these emergent properties can be used to validate the design decisions of the model, such as the range of phytoplankton functional types and their behaviour, the stoichiometric flexibility with regards to each nutrient, and the choice of fixed or variable carbon to nitrogen ratios.

Feeding selectivity of bivalve larvae on natural plankton assemblages in the Western English Channel

Meroplankton, including bivalve larvae, are an important and yet understudied component of coastal marine food webs. Understanding the baseline of meroplankton ecology is imperative to establish and predict their sensitivity to local and global marine stressors. Over an annual cycle (October 2009–September 2010), bivalve larvae were collected from the Western Channel Observatory time series station L4 (50°15.00′N, 4°13.02′W). The morphologically similar larvae were identified by analysis of the 18S nuclear small subunit ribosomal RNA gene, and a series of incubation experiments were conducted to determine larval ingestion rates on natural plankton assemblages. Complementary gut content analysis was performed using a PCR-based method for detecting prey DNA both from field-collected larvae and those from the feeding experiments. Molecular identification of bivalve larvae showed the community composition to change over the course of the sampling period with domination by Phaxas in winter and higher diversity in autumn. The larvae selected for nanoeukaryotes (2–20 µm) including coccolithophores (<20 µm) which together comprised >75 % of the bivalve larvae diet. Additionally, a small percentage of carbon ingested originated from heterotrophic ciliates (<30 µm). The molecular analysis of bivalve larvae gut content provided increased resolution of identification of prey consumed and demonstrated that the composition of prey consumed established through bottle incubations conferred with that established from in situ larvae. Despite changes in bivalve larvae community structure, clearance rates of each prey type did not change significantly over the course of the experiment, suggesting different bivalve larvae species may consume similar prey.

Next Generation Sequencing Reveals the Hidden Diversity of Zooplankton Assemblages

Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel next generation sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify richness and diversity of a mixed zooplankton assemblage from a productive time series site in the Western English Channel. Methodology/Principle Findings Plankton net hauls (200 µm) were taken at the Western Channel Observatory station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,041 sequences were obtained for all samples. The sequences clustered into 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 135 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 58 taxonomic groups. Conclusions Metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and hard-to-identify meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for elucidating the true diversity and species richness of zooplankton communities. While this approach allows for broad diversity assessments of plankton it may become increasingly attractive in future if sequence reference libraries of accurately identified individuals are better populated.

An unusually large phytoplankton spring bloom drives rapid changes in benthic diversity and ecosystem function

In 2012, the Western English Channel experienced an unusually large and long-lived phytoplankton spring bloom. When compared with data from the past 20 years, average phytoplankton biomass at Station L4 (part of the Western Channel Observatory) was approximately 3× greater and lasted 50% longer than any previous year. Regular (mostly weekly) box core samples were collected from this site before, during and after the bloom to determine its impact on macrofaunal abundance, diversity, biomass, community structure and function. The spring bloom of 2012 was shown to support a large and rapid response in the majority of benthic taxa and functional groups. However, key differences in the precise nature of this response, as well as in its timing, was observed between different macrofauna feeding groups. Deposit feeders responded almost instantly at the start of the bloom, primarily thorough an increase in abundance. Suspension feeders and opportunistic/predatory/carnivorous taxa responded slightly more slowly and primarily with an increase in biomass. At the end of the bloom a rapid decline in macrobenthic abundance, diversity and biomass closely followed the decline in phytoplankton biomass. With suspension feeders showing evidence of this decline a few weeks before deposit feeders, it was concluded that this collapse in benthic communities was driven primarily by food availability and competition. However, it is possible that environmental hypoxia and the presence of toxic benthic cyanobacteria could also have contributed to this decline. This study shows evidence for strong benthic–pelagic coupling at L4; a shallow (50 m), coastal, fine-sand habitat. It also demonstrates that in such habitats, it is not just planktonic organisms that demonstrate clear community phenology. Different functional groups within the benthic assemblage will respond to the spring bloom in specific manner, with implications for key ecosystem functions and processes, such as secondary production and bioturbation. Only by taking integrated benthic and pelagic observations over such fine temporal scales (weekly) was the current study able to identify the intimate structure of the benthic response. Similar studies from other habitats and under different bloom conditions are urgently needed to fully appreciate the strength of benthic–pelagic coupling in shallow coastal environments.

Modelling the future biogeography of North Atlantic zooplankton communities in response to climate change

Advances in habitat and climate modelling allow us to reduce uncertainties of climate change impacts on species distribution. We evaluated the impacts of future climate change on community structure, diversity, distribution and phenology of 14 copepod species in the North Atlantic. We developed and validated habitat models for key zooplankton species using continuous plankton recorder (CPR) survey data collected at mid latitudes of the North Atlantic. Generalized additive models (GAMs) were applied to relate the occurrence of species to environmental variables. Models were projected to future (2080–2099) environmental conditions using coupled hydroclimatix–biogeochemical models under the Intergovernmental Panel on Climate Change (IPCC) A1B climate scenario, and compared to present (2001–2020) conditions. Our projections indicated that the copepod community is expected to respond substantially to climate change: a mean poleward latitudinal shift of 8.7 km per decade for the overall community with an important species range variation (–15 to 18 km per decade); the species seasonal peak is expected to occur 12–13 d earlier for Calanus finmarchicus and C. hyperboreus; and important changes in community structure are also expected (high species turnover of 43–79% south of the Oceanic Polar Front). The impacts of the change expected by the end of the century under IPCC global warming scenarios on copepods highlight poleward shifts, earlier seasonal peak and changes in biodiversity spatial patterns that might lead to alterations of the future North Atlantic pelagic ecosystem. Our model and projections are supported by a temporal validation undertaken using the North Atlantic climate regime shift that occurred in the 1980s: the habitat model built in the cold period (1970–1986) has been validated in the warm period (1987–2004).

Relationships between biodiversity and the stability of marine ecosystems: comparisons at a European scale using meta-analysis.

The relationship between biodiversity and stability of marine benthic assemblages was investigated using existing data sets (n = 28) covering various spatial (m-km) and temporal (1973-2006) scales in different benthic habitats (emergent rock, rock pools and sedimentary habitats) through meta-analyses. Assemblage stability was estimated by measuring temporal variances of species richness, total abundance (density or % cover) and community species composition and abundance structure (using multivariate analyses). Positive relationships between temporal variability in species number and richness were generally observed at both quadrat (<1 m2) and site (100 m2) scales, while no relationships were observed by multivariate analyses. Positive relationships were also observed at the scale of site between temporal variability in species number and variability in community structure with evenness estimates. This implies that the relationship between species richness or evenness and species richness variability is slightly positive and depends on the scale of observation, suggesting that biodiversity per se is important for the stability of ecosystems. Changes within community assemblages in terms of structure are, however, generally independent of biodiversity, suggesting no effect of diversity, but the potential impact of individual species, and/or environmental factors. Except for sedimentary and rock pool habitats, no relationship was observed between temporal variation of the aggregated variable of total abundances and diversity at either scale. Overall our results emphasise that relationships depend on scale of measurements, type of habitats and the marine systems (North Atlantic and Mediterranean) considered.