Using continuous plankton recorder data

The continuous plankton recorder (CPR) survey is the largest multi-decadal plankton monitoring programme in the world. It was initiated in 1931 and by the end of 2004 had counted 207,619 samples and identified 437 phyto- and zooplankton taxa throughout the North Atlantic. CPR data are used extensively by the research community and in recent years have been used increasingly to underpin marine management. Here, we take a critical look at how best to use CPR data. We first describe the CPR itself, CPR sampling, and plankton counting procedures. We discuss the spatial and temporal biases in the Survey, summarise environmental data that have not previously been available, and describe the new data access policy. We supply information essential to using CPR data, including descriptions of each CPR taxonomic entity, the idiosyncrasies associated with counting many of the taxa, the logic behind taxonomic changes in the Survey, the semi-quantitative nature of CPR sampling, and recommendations on choosing the spatial and temporal scale of study. This forms the basis for a broader discussion on how to use CPR data for deriving ecologically meaningful indices based on size, functional groups and biomass that can be used to support research and management. This contribution should be useful for plankton ecologists, modellers and policy makers that actively use CPR data.

The CPR Survey (1948-1997): a gridded database browser of plankton abundance in the North Sea.

The Continuous Plankton Recorder (CPR) survey provides a unique multi- decadal dataset on the abundance of plankton in the North Sea and North Atlantic and is one of only a few monitoring programmes operating at a large spatio- temporal scale. The results of all samples analysed from the survey since 1946 are stored on an Access Database at the Sir Alister Hardy Foundation for Ocean Science (SAHFOS) in Plymouth. The database is large, containing more than two million records (~80 million data points, if zero results are added) for more than 450 taxonomic entities. An open data policy is operated by SAHFOS. However, the data are not on-line and so access by scientists and others wishing to use the results is not interactive. Requests for data are dealt with by the Database Manager. To facilitate access to the data from the North Sea, which is an area of high research interest, a selected set of data for key phytoplankton and zooplankton species has been processed in a form that makes them readily available on CD for research and other applications. A set of MATLAB tools has been developed to provide an interpolated spatio-temporal description of plankton sampled by the CPR in the North Sea, as well as easy and fast access to users in the form of a browser. Using geostatistical techniques, plankton abundance values have been interpolated on a regular grid covering the North Sea. The grid is established on centres of 1 degree longitude x 0.5 degree latitude (~32 x 30 nautical miles). Based on a monthly temporal resolution over a fifty-year period (1948-1997), 600 distribution maps have been produced for 54 zooplankton species, and 480 distribution maps for 57 phytoplankton species over the shorter period 1958-1997. The gridded database has been developed in a user-friendly form and incorporates, as a package on a CD, a set of options for visualisation and interpretation, including the facility to plot maps for selected species by month, year, groups of months or years, long-term means or as time series and contour plots. This study constitutes the first application of an easily accessed and interactive gridded database of plankton abundance in the North Sea. As a further development the MATLAB browser is being converted to a user- friendly Windows-compatible format (WinCPR) for release on CD and via the Web in 2003.

Automatic classification of field-collected dinoflagellates by artificial neural network

Automatic taxonomic categorisation of 23 species of dinoflagellates was demonstrated using field-collected specimens. These dinoflagellates have been responsible for the majority of toxic and noxious phytoplankton blooms which have occurred in the coastal waters of the European Union in recent years and make severe impact on the aquaculture industry. The performance by human 'expert' ecologists/taxonomists in identifying these species was compared to that achieved by 2 artificial neural network classifiers (multilayer perceptron and radial basis function networks) and 2 other statistical techniques, k-Nearest Neighbour and Quadratic Discriminant Analysis. The neural network classifiers outperform the classical statistical techniques. Over extended trials, the human experts averaged 85% while the radial basis network achieved a best performance of 83%, the multilayer perceptron 66%, k-Nearest Neighbour 60%, and the Quadratic Discriminant Analysis 56%.

Defining seasonal marine microbial community dynamics

Here we describe, the longest microbial time-series analyzed to date using high-resolution 16S rRNA tag pyrosequencing of samples taken monthly over 6 years at a temperate marine coastal site off Plymouth, UK. Data treatment effected the estimation of community richness over a 6-year period, whereby 8794 operational taxonomic units (OTUs) were identified using single-linkage preclustering and 21 130 OTUs were identified by denoising the data. The Alphaproteobacteria were the most abundant Class, and the most frequently recorded OTUs were members of the Rickettsiales (SAR 11) and Rhodobacteriales. This near-surface ocean bacterial community showed strong repeatable seasonal patterns, which were defined by winter peaks in diversity across all years. Environmental variables explained far more variation in seasonally predictable bacteria than did data on protists or metazoan biomass. Change in day length alone explains >65% of the variance in community diversity. The results suggested that seasonal changes in environmental variables are more important than trophic interactions. Interestingly, microbial association network analysis showed that correlations in abundance were stronger within bacterial taxa rather than between bacteria and eukaryotes, or between bacteria and environmental variables.

Pattern of benthic biomass size spectra from shallow waters in the East China Seas

Benthic biomass size spectra (BSS) and normalized biomass size spectra were constructed, and benthic secondary production was estimated by a size spectrum equation in the shallow waters in the East China Sea, ranging latitudinally from 40A degrees N to 29A degrees N. The BSS patterns were bimodal, two biomass peaks corresponding to meiofauna and macrofauna, respectively, separated by a trough of low biomass at 8-256 mu g individual dry weight which varied in position with median sediment particle size. The BSS also displayed bimodality within meiofauna size ranges, which in most stations was due to the relative proportions of nematodes and other meiofauna taxa. Re-analysis of data from sites in the UK, South Africa, and Antarctic showed a similar bimodality in the adult species body size distribution within the meiofauna size range. Macrofaunal production estimated by the size spectrum equation was very similar to the results of Brey90 empirical equation. However, these production values were much lower than those calculated by Brey01. Different individual dry-to-wet conversion ratios, temperature deviation, and macrofauna taxonomic composition might be responsible for the between-model differences. The macrofaunal P/B ratios calculated by this equation ranged from 0.3 to 3.4 which were in accordance with values from Northern Hemisphere mid-latitudes. Meiofaunal production estimates will need further empirical support.

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.

Minor impact of ocean acidification to the composition of the active microbial community in an Arctic sediment

Effects of ocean acidification on the composition of the active bacterial and archaeal community within Arctic surface sediment was analysed in detail using 16S rRNA 454 pyrosequencing. Intact sediment cores were collected and exposed to one of five different pCO(2) concentrations [380 (present day), 540, 750, 1120 and 3000 atm] and RNA extracted after a period of 14 days exposure. Measurements of diversity and multivariate similarity indicated very little difference between pCO(2) treatments. Only when the highest and lowest pCO(2) treatments were compared were significant differences evident, namely increases in the abundance of operational taxonomic units most closely related to the Halobacteria and differences to the presence/absence structure of the Planctomycetes. The relative abundance of members of the classes Planctomycetacia and Nitrospira increased with increasing pCO(2) concentration, indicating that these groups may be able to take advantage of changing pH or pCO(2) conditions. The modest response of the active microbial communities associated with these sediments may be due to the low and fluctuating pore-water pH already experienced by sediment microbes, a result of the pH buffering capacity of marine sediments, or due to currently unknown factors. Further research is required to fully understand the impact of elevated CO2 on sediment physicochemical parameters, biogeochemistry and microbial community dynamics.

On the roles of cell size and trophic strategy in North Atlantic diatom and dinoflagellate communities

We have examined the inter- and intra-group seasonal succession of 113 diatom and dinoflagellate taxa, as surveyed by the Continuous Plankton Recorder (CPR) in the North Atlantic, by grouping taxa according to two key functional traits: cell size (mg C cell21) and trophic strategy (photoautotrophy, mixotrophy, or heterotrophy). Mixotrophic dinoflagellates follow photoautotrophic diatoms but precede their obligate heterotrophic counterparts in the succession because of the relative advantages afforded by photosynthesizing when light and nutrients are available in spring. The mean cell size of the sampled diatoms is smallest in the summer, likely because of the higher specific nutrient affinity of smaller relative to larger cells. Contrastingly, we hypothesize that mixotrophy diminishes the size selection based on nutrient limitation and accounts for the lack of a seasonal size shift among surveyed dinoflagellates. Relatively small, heterotrophic dinoflagellates (mg C cell21 , 1023) peak after other, larger dinoflagellates, in part because of the increased abundance of their small prey during nutrientdeplete summer months. The largest surveyed diatoms (mg C cell21 . 1022) bloom later than others, and we hypothesize that this may be because of their relatively slow maximum potential growth rates and high internal nutrient storage, as well as to the slower predation of these larger cells. The new trait database and analysis presented here helps translate the taxonomic information of the CPR survey into metrics that can be directly compared with trait-based models.

Marine Zooplankton of Southern Britain - Part 1: Radiolaria, Heliozoa, Foraminifera, Ciliophora, Cnidaria, Ctenophora, Platyhelminthes, Nemertea, Rotifera and Mollusca. Occasional Publication of the Marine Biological Association 25

This series of three guides (of which this is Part 1) collates taxonomic identification information for the zooplankton groups recorded off south-west Britain , primarily for local identification and training purposes. However, because prevailing currents also bring oceanic zooplankton into the English Channel , the range of species sampled off Plymouth covers the majority found over the shallower parts of northern European continental shelf (excluding the Mediterranean Sea ), so the guides should be more widely useful and hopefully make tackling zooplankton identification easier for a wider audience. The commonest truly planktonic species and the most widely studied groups are covered in most detail, but some information is also included on benthic, epibenthic and parasitic species that are sampled occasionally. For all groups there is at least information on their morphology, guidance on their identification and bibliographies giving identification resources.

Marine Zooplankton of Southern Britain - Part 2: Arachnida, Pycnogonida, Cladocera, Facetotecta, Cirripedia and Copepoda. Occasional Publication of the Marine Biological Association 26

This series of three guides (of which this is Part 2) collates taxonomic identification information for the zooplankton groups recorded off south-west Britain , primarily for local identification and training purposes. However, because prevailing currents also bring oceanic zooplankton into the English Channel , the range of species sampled off Plymouth covers the majority found over the shallower parts of northern European continental shelf (excluding the Mediterranean Sea ), so the guides should be more widely useful and hopefully make tackling zooplankton identification easier for a wider audience. The commonest truly planktonic species and the most widely studied groups are covered in most detail, but some information is also included on benthic, epibenthic and parasitic species that are sampled occasionally. For all groups there is at least information on their morphology, guidance on their identification and bibliographies giving identification resources.

A new procedure to optimize the selection of groups in a classification tree

Agglomerative cluster analyses encompass many techniques, which have been widely used in various fields of science. In biology, and specifically ecology, datasets are generally highly variable and may contain outliers, which increase the difficulty to identify the number of clusters. Here we present a new criterion to determine statistically the optimal level of partition in a classification tree. The criterion robustness is tested against perturbated data (outliers) using an observation or variable with values randomly generated. The technique, called Random Simulation Test (RST), is tested on (1) the well-known Iris dataset [Fisher, R.A., 1936. The use of multiple measurements in taxonomic problems. Ann. Eugenic. 7, 179–188], (2) simulated data with predetermined numbers of clusters following Milligan and Cooper [Milligan, G.W., Cooper, M.C., 1985. An examination of procedures for determining the number of clusters in a data set. Psychometrika 50, 159–179] and finally (3) is applied on real copepod communities data previously analyzed in Beaugrand et al. [Beaugrand, G., Ibanez, F., Lindley, J.A., Reid, P.C., 2002. Diversity of calanoid copepods in the North Atlantic and adjacent seas: species associations and biogeography. Mar. Ecol. Prog. Ser. 232, 179–195]. The technique is compared to several standard techniques. RST performed generally better than existing algorithms on simulated data and proved to be especially efficient with highly variable datasets.

Monitoring changes in phytoplankton abundance and composition in the Northwest Atlantic: a comparison of results obtained by Continuous Plankton Recorder sampling and colour satellite imagery

Phytoplankton abundance in the NW Atlantic was measured by continuous plankton recorder (CPR) sampling along tracks between Iceland and the western Scotian Shelf from 1998 to 2006, when sea-surface chlorophyll (SSChl) measurements were also being made by ocean colour satellite imagery using the SeaWiFS sensor. Seasonal and inter-annual changes in phytoplankton abundance were examined using data collected by both techniques, averaged over each of four shelf regions and four deep ocean regions. CPR sampling had gaps (missing months) in all regions and in the four deep ocean regions satellite observations were too sparse between November and February to be of use. Average seasonal cycles of SSChl were similar to those of total diatom abundance in seven regions, to those of the phytoplankton colour index in six regions, but were not similar to those of total dinoflagellate abundance anywhere. Large inter-annual changes in spring bloom dynamics were captured by both samplers in shelf regions. Changes in annual (or 8 months) averages of SSChl did not generally follow those of the CPR indices within regions and multi-year averages of SSChl, and the three CPR indices were generally higher in shelf than in deep ocean regions. Remote sensing and CPR sampling provide complementary ways of monitoring phytoplankton in the ocean: the former has superior temporal and spatial coverage and temporal resolution, and the latter provides better taxonomic information.

Spatial and temporal variability in Scotian Shelf zooplankton communities

Zooplankton are indicators of the response of marine ecosystems to environmental variability. The relationships between zooplankton communities on the Scotian Shelf and hydrographic and geographic regions of the Scotian Shelf in the 1990s and 2000s were described using complementary data sets, each resolving different space and time scales. The Atlantic Zone Monitoring Program (AZMP) sampled Scotian Shelf zooplankton from the whole water column twice per year at stations along three cross-shelf transects and semi-monthly at a fixed station on the inshore central shelf, while Continuous Plankton Recorder (CPR) samples were collected from near-surface waters approximately monthly on an along-shelf transect. Variability patterns were compared among these three data sets to identify robust spatial and interannual trends. Stations were clustered based on taxonomic composition, and spatial clusters were compared to hydrographic boundaries and bathymetry to determine whether temporal changes in community composition were driven by changes in water mass distributions on the shelf. This project identifies zooplankton community and abundance shifts that may affect fish recruitment in the northwest Atlantic and contributes to development of ecosystem-based fisheries management on the Scotian Shelf.

US NOAA Fisheries and UK SAHFOS CPR surveys: comparison of methods and data

The US National Oceanic and Atmospheric Administration (NOAA) Fisheries Continuous Plankton Recorder (CPR) Survey has sampled four routes: Boston–Nova Scotia (1961–present), New York toward Bermuda (1976–present), Narragansett Bay–Mount Hope Bay–Rhode Island Sound (1998–present) and eastward of Chesapeake Bay (1974–1980). NOAA involvement began in 1974 when it assumed responsibility for the existing Boston–Nova Scotia route from what is now the UK's Sir Alister Hardy Foundation for Ocean Science (SAHFOS). Training, equipment and computer software were provided by SAHFOS to ensure continuity for this and standard protocols for any new routes. Data for the first 14 years of this route were provided to NOAA by SAHFOS. Comparison of collection methods; sample processing; and sample identification, staging and counting techniques revealed near-consistency between NOAA and SAHFOS. One departure involved phytoplankton counting standards. This has since been addressed and the data corrected. Within- and between-survey taxonomic and life-stage names and their consistency through time were, and continue to be, an issue. For this, a cross-reference table has been generated that contains the SAHFOS taxonomic code, NOAA taxonomic code, NOAA life-stage code, National Oceanographic Data Center (NODC) taxonomic code, Integrated Taxonomic Information System (ITIS) serial number and authority and consistent use/route. This table is available for review/use by other CPR surveys. Details of the NOAA and SAHFOS comparison and analytical techniques unique to NOAA are presented.

Genetic analysis of Continuous Plankton Recorder (CPR) samples

Genetic analysis of Continuous Plankton Recorder (CPR) samples is enabling greater taxonomic resolution and the study of plankton population structure. Here, we present some results from the genetic analysis of CPR samples collected in the North Sea and north-eastern Atlantic that reveal the impacts of climate on benthic-pelagic coupling and the food web. We show that pronounced changes in the North Sea meroplankton are related to an increased abundance and spatial distribution of the larvae of the benthic echinoderm, Echinocardium cordatum. Key stages of reproduction in E. cordatum, gametogenesis and spawning, are influenced by winter and spring sea temperature (January-May). A stepwise increase in sea temperature after 1987, which has created warmer conditions earlier in the year, together with increased summer phytoplankton, may benefit the reproduction and survival of this benthic species. Competition between the larvae of E. cordatum and other holozooplanlcton taxa may now be altering the trophodynamics of the summer pelagic ecosystem. In the north-eastern Atlantic the genetic analysis of fish larvae sampled by the CPR has revealed an unprecedented increase in the abundance of juvenile snake pipefish, Entelurus aequoreiis since 2002. We argue that increased sea surface temperatures in winter and spring when the eggs of E. aqueoreus, which are brooded by the male, are developing and the young larvae are growing in the plankton are a likely cause. The increased abundance of this species in Atlantic and adjacent European seas already appears to be influencing the marine food web.