Microbial spatial variability: An example from the Celtic Sea

In July 2004, dominant populations of microbial ultraplankton (<5 μm), in the surface of the Celtic Sea (between UK and Eire), were repeatedly mapped using flow cytometry, at 1.5 km resolution over a region of diameter 100 km. The numerically dominant representatives of all basic functional types were enumerated including one group of phototrophic bacteria (Syn), two groups of phytoplankton (PP, NP), three groups of heterotrophic bacterioplankton (HB) and the regionally dominant group of heterotrophic protists (HP). The distributions of all organisms showed strong spatial variability with little relation to variability in physical fields such as salinity and temperature. Furthermore, there was little agreement between distributions of different organisms. The only linear correlation consistently explaining more than 50% of the variance between any pairing of the organism groups enumerated is between two different groups of HB. Specifically, no linear, or non-linear, relationship is found between any pairings of SYB, PP or HB groups with their protist predators HP. Looking for multiple dependencies, factor analysis reveals three groupings: Syn, PP and low nucleic acid content HB (LNA); high nucleic acid content HB (HNA); HP and NP. Even the manner in which the spatial variability of Syn, PP and HB abundance varies as a function of lengthscale (represented by a semivariogram) differs significantly from that for HP. In summary, although all microbial planktonic groups enumerated are present and numerically dominant throughout the region studied, at face value the relationships between them seem weak. Nevertheless, the behaviour of a simple, illustrative ecological model, with strongly interacting phototrophs and heterotrophs, with stochastic forcing, is shown to be consistent with the observed poor correlations and differences in how spatial variability varies with lengthscale. Thus, our study suggests that a comparison of microbial abundances alone may not discern strong underlying trophic interactions. Specific knowledge of these processes, in particular grazing, will be required to explain the causes of the observed microbial spatial variability and its resulting consequences for the functioning of the ecosystem.

Response of an Arctic Sediment Nitrogen Cycling Community to Increased CO2

Ocean acidification influences sediment/water nitrogen fluxes, possibly by impacting on the microbial process of ammonia oxidation. To investigate this further, undisturbed sediment cores collected from Ny Alesund harbour (Svalbard) were incubated with seawater adjusted to CO2 concentrations of 380, 540, 760, 1,120 and 3,000 μatm. DNA and RNA were extracted from the sediment surface after 14 days' exposure and the abundance of bacterial and archaeal ammonia oxidising (amoA) genes and transcripts quantified using quantitative polymerase chain reaction. While there was no change to the abundance of bacterial amoA genes, an increase to 760 μatm pCO2 reduced the abundance of bacterial amoA transcripts by 65 %, and this was accompanied by a shift in the composition of the active community. In contrast, archaeal amoA gene and transcript abundance both doubled at 3,000 μatm, with an increase in species richness also apparent. This suggests that ammonia oxidising bacteria and archaea in marine sediments have different pH optima, and the impact of elevated CO2 on N cycling may be dependent on the relative abundances of these two major microbial groups. Further evidence of a shift in the balance of key N cycling groups was also evident: the abundance of nirS-type denitrifier transcripts decreased alongside bacterial amoA transcripts, indicating that NO3 − produced by bacterial nitrification fuelled denitrification. An increase in the abundance of Planctomycete-specific 16S rRNA, the vastmajority of which grouped with known anammox bacteria, was also apparent at 3,000 μatm pCO2. This could indicate a possible shift from coupled nitrification–denitrification to anammox activity at elevated CO2.

Satellite estimates of net community production indicate predominance of net autotrophy in the Atlantic Ocean

There is ongoing debate as to whether the oligotrophic ocean is predominantly net autotrophic and acts as a CO2 sink, or net heterotrophic and therefore acts as a CO2 source to the atmosphere. This quantification is challenging, both spatially and temporally, due to the sparseness of measurements. There has been a concerted effort to derive accurate estimates of phytoplankton photosynthesis and primary production from satellite data to fill these gaps; however there have been few satellite estimates of net community production (NCP). In this paper, we compare a number of empirical approaches to estimate NCP from satellite data with in vitro measurements of changes in dissolved O2 concentration at 295 stations in the N and S Atlantic Ocean (including the Antarctic), Greenland and Mediterranean Seas. Algorithms based on power laws between NCP and particulate organic carbon production (POC) derived from 14C uptake tend to overestimate NCP at negative values and underestimate at positive values. An algorithm that includes sea surface temperature (SST) in the power function of NCP and 14C POC has the lowest bias and root-mean square error compared with in vitro measured NCP and is the most accurate algorithm for the Atlantic Ocean. Nearly a 13 year time series of NCP was generated using this algorithm with SeaWiFS data to assess changes over time in different regions and in relation to climate variability. The North Atlantic subtropical and tropical Gyres (NATL) were predominantly net autotrophic from 1998 to 2010 except for boreal autumn/winter, suggesting that the northern hemisphere has remained a net sink for CO2 during this period. The South Atlantic subtropical Gyre (SATL) fluctuated from being net autotrophic in austral spring-summer, to net heterotrophic in austral autumn–winter. Recent decadal trends suggest that the SATL is becoming more of a CO2 source. Over the Atlantic basin, the percentage of satellite pixels with negative NCP was 27%, with the largest contributions from the NATL and SATL during boreal and austral autumn–winter, respectively. Variations in NCP in the northern and southern hemispheres were correlated with climate indices. Negative correlations between NCP and the multivariate ENSO index (MEI) occurred in the SATL, which explained up to 60% of the variability in NCP. Similarly there was a negative correlation between NCP and the North Atlantic Oscillation (NAO) in the Southern Sub-Tropical Convergence Zone (SSTC),which explained 90% of the variability. There were also positive correlations with NAO in the Canary Current Coastal Upwelling (CNRY) and Western Tropical Atlantic (WTRA)which explained 80% and 60% of the variability in each province, respectively. MEI and NAO seem to play a role in modifying phases of net autotrophy and heterotrophy in the Atlantic Ocean.

Both respiration and photosynthesis determine the scaling of plankton metabolism in the oligotrophic ocean

Despite its importance to ocean–climate interactions, the metabolic state of the oligotrophic ocean has remained controversial for 415 years. Positions in the debate are that it is either hetero- or autotrophic, which suggests either substantial unaccounted for organic matter inputs, or that all available photosynthesis (P) estimations (including 14C) are biased. Here we show the existence of systematic differences in the metabolic state of the North (heterotrophic) and South (autotrophic) Atlantic oligotrophic gyres, resulting from differences in both P and respiration (R). The oligotrophic ocean is neither auto- nor heterotrophic, but functionally diverse. Our results show that the scaling of plankton metabolism by generalized P:R relationships that has sustained the debate is biased, and indicate that the variability of R, and not only of P, needs to be considered in regional estimations of the ocean’s metabolic state.

Nanoplankton and picoplankton in the Western English Channel: abundance and seasonality from 2007–2013

The nano- and picoplankton community at Station L4 in the Western English Channel was studied between 2007 and 2013 by flow cytometry to quantify abundance and investigate seasonal cycles within these communities. Nanoplankton included both photosynthetic and heterotrophic eukaryotic single-celled organisms while the picoplankton included picoeukaryote phytoplankton, Synechococcus sp. cyanobacteria and heterotrophic bacteria. A Box–Jenkins Transfer Function climatology analysis of surface data revealed that Synechococcus sp., cryptophytes, and heterotrophic flagellates had bimodal annual cycles. Nanoeukaryotes and both high and low nucleic acid-containing bacteria (HNA and LNA, respectively) groups exhibited unimodal annual cycles. Phaeocystis sp., whilst having clearly defined abundance maxima in spring was not detectable the rest of the year. Coccolithophores exhibited a weak seasonal cycle, with abundance peaks in spring and autumn. Picoeukaryotes did not exhibit a discernable seasonal cycle at the surface. Timings of maximum group abundance varied through the year. Phaeocystis sp. and heterotrophic flagellates peaked in April/May. Nanoeukaryotes and HNA bacteria peaked in June/July and had relatively high abundance throughout the summer. Synechococcus sp., cryptophytes and LNA bacteria all peaked from mid to late September. The transfer function model techniques used represent a useful means of identifying repeating annual cycles in time series data with the added ability to detect trends and harmonic terms at different time scales from months to decades.

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.

Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles

Global ocean biogeochemistry models currently employed in climate change projections use highly simplified representations of pelagic food webs. These food webs do not necessarily include critical pathways by which ecosystems interact with ocean biogeochemistry and climate. Here we present a global biogeochemical model which incorporates ecosystem dynamics based on the representation of ten plankton functional types (PFTs); six types of phytoplankton, three types of zooplankton, and heterotrophic bacteria. We improved the representation of zooplankton dynamics in our model through (a) the explicit inclusion of large, slow-growing zooplankton, and (b) the introduction of trophic cascades among the three zooplankton types. We use the model to quantitatively assess the relative roles of iron vs. grazing in determining phytoplankton biomass in the Southern Ocean High Nutrient Low Chlorophyll (HNLC) region during summer. When model simulations do not represent crustacean macrozooplankton grazing, they systematically overestimate Southern Ocean chlorophyll biomass during the summer, even when there was no iron deposition from dust. When model simulations included the developments of the zooplankton component, the simulation of phytoplankton biomass improved and the high chlorophyll summer bias in the Southern Ocean HNLC region largely disappeared. Our model results suggest that the observed low phytoplankton biomass in the Southern Ocean during summer is primarily explained by the dynamics of the Southern Ocean zooplankton community rather than iron limitation. This result has implications for the representation of global biogeochemical cycles in models as zooplankton faecal pellets sink rapidly and partly control the carbon export to the intermediate and deep ocean.