Identifying four phytoplankton functional types from space: An ecological approach

Deriving maps of phytoplankton taxa based on remote sensing data using bio-optical properties of phytoplankton alone is challenging. A more holistic approach was developed using artificial neural networks, incorporating ecological and geographical knowledge together with ocean color, bio-optical characteristics, and remotely sensed physical parameters. Results show that the combined remote sensing approach could discriminate four major phytoplankton functional types (diatoms, dinoflagellates, coccolithophores, and silicoflagellates) with an accuracy of more than 70%. Models indicate that the most important information for phytoplankton functional type discrimination is spatio-temporal information and sea surface temperature. This approach can supply data for large-scale maps of predicted phytoplankton functional types, and an example is shown.

An intercomparison of bio-optical techniques for detecting dominant phytoplankton size class from satellite remote sensing

Satellite remote sensing of ocean colour is the only method currently available for synoptically measuring wide-area properties of ocean ecosystems, such as phytoplankton chlorophyll biomass. Recently, a variety of bio-optical and ecological methods have been established that use satellite data to identify and differentiate between either phytoplankton functional types (PFTs) or phytoplankton size classes (PSCs). In this study, several of these techniques were evaluated against in situ observations to determine their ability to detect dominant phytoplankton size classes (micro-, nano- and picoplankton). The techniques are applied to a 10-year ocean-colour data series from the SeaWiFS satellite sensor and compared with in situ data (6504 samples) from a variety of locations in the global ocean. Results show that spectral-response, ecological and abundance-based approaches can all perform with similar accuracy. Detection of microplankton and picoplankton were generally better than detection of nanoplankton. Abundance-based approaches were shown to provide better spatial retrieval of PSCs. Individual model performance varied according to PSC, input satellite data sources and in situ validation data types. Uncertainty in the comparison procedure and data sources was considered. Improved availability of in situ observations would aid ongoing research in this field.

iMarNet: an ocean biogeochemistry model intercomparison project within a common physical ocean modelling framework

Ocean biogeochemistry (OBGC) models span a wide variety of complexities, including highly simplified nutrient-restoring schemes, nutrient–phytoplankton–zooplankton–detritus (NPZD) models that crudely represent the marine biota, models that represent a broader trophic structure by grouping organisms as plankton functional types (PFTs) based on their biogeochemical role (dynamic green ocean models) and ecosystem models that group organisms by ecological function and trait. OBGC models are now integral components of Earth system models (ESMs), but they compete for computing resources with higher resolution dynamical setups and with other components such as atmospheric chemistry and terrestrial vegetation schemes. As such, the choice of OBGC in ESMs needs to balance model complexity and realism alongside relative computing cost. Here we present an intercomparison of six OBGC models that were candidates for implementation within the next UK Earth system model (UKESM1). The models cover a large range of biological complexity (from 7 to 57 tracers) but all include representations of at least the nitrogen, carbon, alkalinity and oxygen cycles. Each OBGC model was coupled to the ocean general circulation model Nucleus for European Modelling of the Ocean (NEMO) and results from physically identical hindcast simulations were compared. Model skill was evaluated for biogeochemical metrics of global-scale bulk properties using conventional statistical techniques. The computing cost of each model was also measured in standardised tests run at two resource levels. No model is shown to consistently outperform all other models across all metrics. Nonetheless, the simpler models are broadly closer to observations across a number of fields and thus offer a high-efficiency option for ESMs that prioritise high-resolution climate dynamics. However, simpler models provide limited insight into more complex marine biogeochemical processes and ecosystem pathways, and a parallel approach of low-resolution climate dynamics and high-complexity biogeochemistry is desirable in order to provide additional insights into biogeochemistry–climate interactions.

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.

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed. populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (similar to 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

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.

Variability in phytoplankton community structure in response to the North Atlantic Oscillation and implications for organic carbon flux

The North Atlantic Oscillation (NAO) is a major mode of variability in the North Atlantic, dominating atmospheric and oceanic conditions. Here, we examine the phytoplankton community-structure response to the NAO using the Continuous Plankton Recorder data set. In the Northeast Atlantic, in the transition region between the gyres, variability in the relative influence of subpolar or subtropical-like conditions is reflected in the physical environment. During positive NAO periods, the region experiences subpolar-like conditions, with strong wind stress and deep mixed layers. In contrast, during negative NAO periods, the region shifts toward more subtropical-like conditions. Diatoms dominate the phytoplankton community in positive NAO periods, whereas in negative NAO periods, dinoflagellates outcompete diatoms. The implications for interannual variability in deep ocean carbon flux are examined using data from the Porcupine Abyssal Plain time-series station. Contrary to expectations, carbon flux to 3000 m is enhanced when diatoms are outcompeted by other phytoplankton functional types. Additionally, highest carbon fluxes were not associated with an increase in biomineral content, which implies that ballasting is not playing a dominant role in controlling the flux of material to the deep ocean in this region. In transition zones between gyre systems, phytoplankton populations can change in response to forcing induced by opposing NAO phases.