Changes in Phytoplankton and Zooplankton Communities

D1 - Biological Diversity

D4 - Marine Food Webs

D1.4 - Habitat Distribution

D1.6 - Habitat condition

D4.3 - Abundance/distribution of key trophic groups/species

Plankton form the base of the marine food web and respond rapidly to environmental change, making them important indicators of ecosystem state. Between 2004-2014 plankton communities experienced significant changes in relative abundance, indicating alterations to key aspects of ecosystem functioning. Those changes are widely accepted to be linked to prevailing conditions and may be driven by climate change, nutrient enrichment or other factors.

Area assessed

Printable summary

Background

Plankton (microscopic algae and animals) form the base of the marine food web, making them important indicators for ecosystem state. Changes in plankton communities can affect higher food web levels, such as shellfish, fish and seabirds, since these organisms are supported either directly or indirectly by plankton. Indicators based on plankton lifeforms (i.e. organisms with the same functional traits; Figure 1) can be used to reveal plankton community responses to factors such as nutrient loading from human activities and climate-driven change. When examined in pairs with an ecologically-relevant relationship (Figure 1 for examples) changes in the relative abundance of two lifeforms together (called a lifeform pair) can indicate change in key aspects of ecosystem function, including links between pelagic and benthic communities, energy flows and pathways, and food web interactions.

Figure 1: Plankton lifeform pairs and ecological rationale for their selection

Courtesy of the Integration and Application Network, University of Maryland Center for Environmental Science (ian.umces.edu/symbols/).

For example, changes in the phytoplankton (algae plankton) lifeform can cause changes in the zooplankton (animal plankton) lifeform that feeds on them. At the North-East Atlantic regional scale, plankton community change is strongly linked to prevailing climatic conditions. Pelagic habitats, which are defined based on key water column features, are important to plankton community structure and dynamics (Figure 2).

Because this is a new Indicator Assessment in the first phase of development, no assessment value exists. Instead, the years 2004 to 2008 are used as compared to the last 6-year period (2009 to 2014) to examine changes in lifeform pairs.

Figure 2: Ecohydrodynamic zones (EHDs) in the Greater North Sea, Celtic Seas, and Bay of Biscay and Iberian Coast

EHDs are constructed based on key water column features, which are important to plankton community structure and dynamics. Based on water column structure, there are six predominant EHD are types: permanently mixed throughout the year; permanently stratified throughout the year; regions of freshwater influence (ROFIs); seasonally stratified (for about half the year, including summer); intermittently stratified, and; indeterminate regions (inconsistently alternate between the above levels of stratification). Work is on-going to define echohydrodynamic zones in the Bay of Biscay and Iberian Coast.

Indicators based on plankton lifeforms can be used in some hydrographic conditions to assess community response to sewage pollution (Charvet et al., 1998; Tett et al., 2008), anoxia (Rakocinski, 2012), fishing (Bremner et al., 2004), eutrophication (HELCOM, 2012) and climate change (Beaugrand, 2005). Changes in ecologically-relevant pairs of plankton lifeforms, examined together, can provide information on ecosystem structure and energy flow. The precise combinations of lifeforms comprising each lifeform pair will depend on the habitat and the objective of the indicator, for example as a measure for pelagic habitats, food webs, sea floor integrity or eutrophication. An important advantage of these plankton indicators is that the proposed concepts are relatively easily transferable to other European regional seas (Gowen et al., 2011; Rombouts et al., 2013).

In practice, the use of functional groups, such as lifeforms, is often favoured over indicator species since indices of species abundance are frequently subject to large interannual variation, often due to natural physical dynamics rather than anthropogenic stressors (de Jonge, 2007). Cryptic speciation (species that look the same under a microscope) within the plankton community, alongside the limitations of identifying plankton using routine light microscopy techniques, make the true diversity of plankton communities difficult to assess at a species or genus level. Functional group abundance is often less variable because variability in the abundances of the group’s constituent species averages out. Moreover, indicators based on functional groups have been proved relevant for the description of the community’s structure and biodiversity, and are more easily inter-compared than other indicators based on taxonomy  (Estrada et al., 2004; Mouillot et al., 2006; Gallego et al., 2012; Garmendia et al., 2012).

When examined in ecologically-relevant pairs, lifeforms can provide an indication of changes in various aspects of ecosystem function: the transfer of energy from primary to secondary producers (changes in phytoplankton and zooplankton); the pathway of energy flow and top predators (changes in gelatinous zooplankton and fish larvae); benthic / pelagic coupling, i.e. changes in holoplankton (fully planktonic) and meroplankton (only part of the lifecycle is planktonic, the remainder is benthic) (Figure 1) (Gowen et al., 2011). Both abundance and biomass data can be used to inform lifeform pairs, depending on the lifeform in question and data availability from monitoring programmes. As the knowledge base grows, new pairs can be developed as indicators for pressures other than those currently measured.

The Plankton Index (PI) of lifeform pairs has been developed to track changes in the state of the plankton in marine waters over time. The main features of the method are: (i) the grouping of species of planktonic organisms into functional types or lifeforms; (ii) the display of changes in the abundance of each of these lifeforms using a state-space approach; (iii) calculating a PI to quantify possible changes in the state of the plankton relative to baseline or starting conditions; (iv) relating trends in the PI to trends in human pressure and climate change indices.   See Gowen et al. (2011) and Tett et al. (2008) for further technical information on the method.

The OSPAR Quality Status Report (QSR) 2010 highlighted the potential impact of climate change and other human pressures on plankton communities. Phytoplankton chlorophyll and phytoplankton indicator species are also assessed under the Common Procedure for the Identification of the Eutrophication Status of the OSPAR Maritime Area as parameters in problem and potential problem areas with regard to eutrophication. The QSR 2010 did not include a comparable quantitative assessment of changes in plankton communities to that provided by the plankton lifeforms indicator.

State-Space theory:

Tett et al. (2008) proposed to track changes in the state of the phytoplankton community by means of plots in a state-space and calculating a Plankton Community Index, referred to here as a Plankton Index (PI). The conceptual framework is that ecosystems can be viewed as systems with an instantaneous state defined by values of a set of system state variables which are attributes of the system that change with time in response to each other and external conditions. Building on this approach and plotting plankton lifeform abundance in a multi-dimensional state-space provides a means of monitoring changes in the structure of plankton communities. A state can be defined as a single point in state-space, with co-ordinates provided by the values of the set of state variables, in this case two lifeform abundances, which together are used as a pair.

In the example illustrated in Figure a for diatoms and dinoflagellates, the axes of the two dimensional space are the abundance values of the two lifeforms. Each point represents the state of the ecosystem in terms of the two lifeforms at the time the water sample was collected. Subsequent samples yield additional pairs of abundance values that can be mapped onto the lifeform state-space.

 

Figure a: Mapping the abundance of two lifeforms in state-space

Point A is the ecosystem state at the instant a water sample was taken and is characterised by the abundance values of two lifeforms. Another sample, taken in the same location, yields abundance values that map to a different point in the diatom-dinoflagellate state-space (point B).

The path between the two states is called a trajectory, and the condition of the phytoplankton is defined by the trajectory drawn in the state-space by a set of points. Such trajectories reflect: (i) cyclic and medium-term variability (the higher order consistencies in the plankton that result from seasonal cycles, species succession and interannual variability); and (ii) long-term variability that might result from environmental pressure. The seasonal nature of plankton production and the succession of species in seasonally stratifying seas, result in this trajectory tending in a certain direction (as plankton growth increases in the spring and declines during autumn), such that the trajectory tends towards its starting point. Given roughly constant external pressures, the data collected from a particular location over a period of years form a cloud of points in state-space that can be referred to as a regime. Long-term variability may show a persistent trend of movement away from a starting point in state-space.

To define a regime, an envelope can be drawn about this group of points using a convex hull method. Because of theoretical arguments that the envelope should be doughnut-shaped with a central hole, bounding curves can be fitted outside and inside the cloud of points (Figure b).

 

Figure b. Creating the envelope in three steps, from left (i) to right (iii): An example of a regime defined by the envelope drawn by the convex hull method

The plot is displayed on a logarithmic scale because this is a common method for showing the full extent of seasonal variability. The data are from the ecohydrodynamic zone that best represents the indeterminate area within the Greater North Sea (EHD II Indeterminate, see Figure 2). The colour of each point corresponds to the season within which it was sampled; blue = winter, aqua = spring, yellow = summer, purple = autumn.

The size and shape of the envelope are sensitive to sampling frequency and the total number of samples. Envelopes are made larger by including extreme outer or inner points, and the larger the envelope, the less sensitive it will be to change in the distribution of points in state-space and therefore to detect a change in condition. Conversely, if too many points are excluded the envelope will be small and even minor changes will result in a statistically significant difference. It is therefore desirable to exclude a proportion (p) of points, to eliminate these extremes, and so the 90-percentile was used. Envelopes are therefore drawn around the cloud of points to include a proportion (p=0.9) of the points: with the 5% of points that were most distant from the cloud's centre and the 5% of points that were closest to the centre, excluded.

For a Plankton Index to be calculated, it is necessary to establish starting conditions as the basis for making comparisons. In the example shown in Figure b, data collected from the indeterminate area within the Greater North Sea (Figure 2) between 2004 and 2008, were used to create an envelope. The envelope, thus drawn (Figure c left) defines a domain in state-space that contains a set of trajectories of the diatom-dinoflagellate component of the marine pelagic ecosystem and thus represents the prevailing regime during the starting period.

The next step is to map a new set of data into the starting condition state-space and compare these new data, at least a dozen points (Figure c right). The value of the PI is the proportion of new points that fall inside the envelope (or to be precise, between the inner and outer envelopes). In the example shown in Figure c, 24% or 17 of the 71 new points lie outside, and the PI is 0.76. A value of 1.0 would indicate no change, and a value of zero would show complete change, with all new points plotting outside the starting condition envelope. The envelope was made by excluding 10% of points, so some new points are expected to fall outside: 7.1 (10% of 71), in the case of the example. The exact probability of getting 17 by chance alone when only 7.1 are expected, can be calculated using a binomial series expansion, or by a chi-square calculation (with 1 df and a 1-tail test). The conclusion is that the value of 0.76, is significantly less than the expected value of 0.9, and so the condition of the phytoplankton in this region, as determined by diatoms and dinoflagellates, was statistically significantly (p<0.01) different between the two periods.

 

 

Figure c: An example for the comparison between two periods within the indeterminate area of the Greater North Sea (Figure 2), left) is the reference period of 2004–2008, right) is the comparison period of 2009–2014

The PI value of 0.76 shows a significant change between the two communities with a binomial p-value of 0.0005 (where significance = p-value <0.01).

The time series of the index is produced by comparing the starting condition envelope with data collected from each subsequent year.

Two lifeforms will not be sufficient to describe all of the important characteristics of the plankton. In principle, there is no constraint to adding more lifeforms to the state-space plots (see Tett et al., 2013). The rule is that each additional lifeform must be different to those already used and the axis for each new lifeform must be drawn at right angles to all existing axes. The state-space map must therefore be drawn in as many dimensions as there are state variables (lifeforms) but this becomes complicated when considering the number of lifeforms that may need to be used to fully represent the community structure of the plankton.

  • Two dimensional state-space plots of specific pairs of plankton lifeforms can be used to investigate pelagic habitats and can be combined to provide a holistic plankton indicator to track changes in the condition of the planktonic component of the pelagic ecosystem.
  • Time series of the index will be used to track persistent changes in the condition of the plankton over time and to relate any such trend to trends in human and climate pressures.
 

Starting Conditions and Analysis:

The period 2004–2008 was selected to represent the starting conditions because each dataset used here (see section ‘The plankton data’) has sufficient data for that period to construct a reliable starting condition envelope. Additional European datasets will be added in the future as resources become available. This envelope represents starting conditions and not reference conditions. The starting condition envelope was compared with data from the subsequent six-year period (2009–2014). The 2009–2014 period was chosen because it is the most recent period with data and is the same temporal length as the European Union Marine Strategy Framework Directive (MSFD) assessment period.

  • Starting conditions do not represent reference conditions.
 

Spatial Scale:

Because plankton community composition, distribution, and dynamics are closely linked to their environment, the analysis was performed at the scale of the ‘ecohydrodynamic’ zones (EHD) (van Leeuwen et al., 2015). EHDs were determined through analysis of a 50-year hindcast using the General Estuarine Transport Model (GETM) physical model of North Sea hydrodynamics (the model is at lower resolution in the Celtic Seas and is not developed for the Bay of Biscay and Iberian Coast). EHDs are constructed using density stratification, the most important large-scale physical feature in shallow shelf seas (van Leeuwen et al., 2015). Density stratification occurs when the buoyancy of surface waters (influenced by freshwater input or solar heating) is stronger than turbulence and vertical mixing, which limits vertical exchange across the pycnocline (van Leeuwen et al., 2015). The predominant EHD zone types, based on water column structure, are: permanently mixed throughout the year; permanently stratified throughout the year; regions of freshwater influence (ROFIs); seasonally stratified (for about half the year, including summer); intermittently stratified; and indeterminate regions (inconsistently alternate between the above).

At the time of calculating this assessment (2016), the EHD model was more reliable and detailed for the Greater North Sea than for the Celtic Seas. The model had not been developed for the Bay of Biscay and Iberian Coast so that sub-region was considered as a single region of indeterminate status.

  • EHDs are a geographical representation of the conditions which best suit plankton distribution, dynamics and community composition
  • EHDs are constructed based on water column density stratification
 

The Plankton Data:

The assessment has been carried out using phytoplankton and zooplankton data from three sources. The data submitted by the Swedish Meteorological and Hydrological Institute (SMHI) and Plymouth Marine Laboratory (PML, United Kingdom) are fixed-point time series. Four nearshore stations from SMHI were used in the analysis, which fall into the EHD zone that describes the indeterminate regions of the Greater North Sea. Both zooplankton and phytoplankton are sampled once or twice per month all year around. The PML station (L4) is 13 km offshore from Plymouth and is sampled for zooplankton and phytoplankton and a suite of other variables on a weekly basis. Unlike the fixed-point stations, data from the Continuous Plankton Recorder (CPR) survey is collected at a much broader spatial scale through ships-of-opportunity. CPR data are collected offshore and in the open ocean and are best analysed at a monthly time scale. It should be noted that all other data for this assessment come from CPR data collection on major shipping lanes. The CPR survey is coordinated by the Sir Alister Hardy Foundation for Ocean Science (SAHFOS) in the United Kingdom. Data from the different providers were not combined for analysis due to differences in sampling, plankton analysis and enumeration methods. Instead, the datasets were analysed separately. Each dataset has internal QA/QC procedures to ensure consistency and accuracy of the data.

As the SMHI stations are located in the indeterminate EHD of the Greater North Sea their data were averaged for analysis. All CPR samples were averaged for each EHD type in each OSPAR sub-region. In the future, more datasets from other OSPAR Contracting Parties will be incorporated into the analysis.

Figure d: Sample locations

The data used for this assessment were obtained from the United Kingdom (CPR and PML L4) and Sweden (SMHI). Future work will expand the plankton database to include data from additional countries.

Lifeform construction:

Each lifeform pair comprises two ecologically-relevant lifeforms, which contain common functional traits (Table a). The rationale for selecting the lifeform pairs and additional criteria containing supplementary information on lifeforms is listed in Table a.

The database master species list was built by assigning functional traits to each species for a dataset and then adding additional new datasets to expand the master species list (Figure e). The species list for each new dataset was assigned a unique Aphia ID via the World Register of Marine Species (WoRMS). The new dataset’s species list was then compared with the plankton database’s master species list via Aphia IDs and any new species were identified. This process ensures that each species is only entered in the database once. The new species were then manually assigned functional traits by searching the literature; fields were left blank where functional traits for species were unknown. Once traits were assigned, the new species were added to the master species list. Queries were constructed to build lifeforms from the functional trait information (Table b).

A simple method of confidence assessment was used to determine the confidence in each lifeform (Table c). Using expert opinion, each lifeform was evaluated on two characteristics: the ability to identify and speciate organisms in that lifeform using light microscopy and the understanding of the accuracy of determining traits assigned to the lifeform. For example, low confidence is assigned to the lifeform pair ‘diatoms vs autotrophic and mixotrophic dinoflagellates’ because the mixotrophic and autotrophic mode of nutrition of many dinoflagellates species is currently uncertain. Thus the accuracy of assigning the lifeform category is low. Likewise, the lifeform pair ‘carnivorous zooplankton vs non-carnivorous zooplankton’ has a low confidence designation since the feeding habits of many abundant and common zooplankton species remain unknown. Further investigation must also be conducted to decide whether both nuisance and toxin-producing lifeform pairs are required. Only pairs with two high-confidence lifeforms were used in the OSPAR reporting.

  • Species are assigned functional traits and combinations of traits are used to construct lifeforms
  • Lifeform pairs consist of ecologically-relevant lifeforms, with each pair having an ecological rationale for selection
  • A simple confidence assessment method was used to evaluate the confidence in each lifeform. Only pairs containing two high-confidence lifeforms have been used for reporting
 

Figure e: Schematic illustrating the process undertaken to assign functional traits to species, then species to lifeforms

Each species must first be assigned a unique Aphia ID to determine whether it is already present in the master species list.

Table a: Lifeform pairs comprise two ecologically-relevant lifeforms
The ‘Additional criteria’ column contains supplementary information regarding particular lifeforms. A simple confidence rating for each lifeform pair, as described in the text, is also included.

Table b: Each lifeform comprises organisms with particular traits. A query is then used to assign lifeforms to individual species
Lifeform Traits Criteria
Diatoms 'Diatom' only PhytoplanktonType=Diatom
Dinoflagellates 'Dinoflagellate' only PhytoplanktonType=Dinoflagellate
Gelatinous zooplankton 'Gelatinous' only PlanktonType=Zooplankton AND Gelatinous=Y
Fish larvae/eggs 'Fish' only ZooType=Fish
Carnivorous zooplankton 'Carnivore' only PlanktonType=Zooplankton AND Diet=Carnivore
Non-carnivorous zooplankton 'Zooplankton' AND either 'Herbivore', 'Omnivore', OR 'Ambiguous' PlanktonType=Zooplankton AND (Diet=Herbivore OR Omnivore OR Ambiguous)
Crustaceans 'Crustacean' only Crustacean=Y
Large phytoplankton 'Phytoplankton' AND 'Lg' PlanktonType=Phytoplankton AND PhytoplanktonSize=Lg
Small phytoplankton 'Phytoplankton' AND 'Sm' PlanktonType=Phytoplankton AND PhytoplanktonSize=Sm
Phytoplankton 'Phytoplankton' only PlanktonType=Phytoplankton
Autotrophic and mixotrophic dinoflagellates 'Dinoflagellate' AND either 'Auto' OR 'Auto/Mixo' PhytoplanktonType=Dinoflagellate AND (FeedingMech=Auto OR Auto/Mixo)
Pelagic diatoms 'Diatom' AND 'Pelagic' PhytoplanktonType=Diatom AND DiatomDepth=Pelagic
Tychopelagic diatoms 'Diatom' AND 'Tychopelagic' PhytoplanktonType=Diatom AND DiatomDepth=Tychopelagic
Nuisance and toxin-producing diatoms 'Diatom' AND either 'Toxic' OR 'Nuisance' PhytoplanktonType=Diatom AND (HAB = Toxic)
Nuisance and toxin-producing dinoflagellates 'Dinoflagellate' AND either 'Toxic' OR 'Nuisance' PhytoplanktonType=Dinoflagellate AND (HAB = Toxic)
Holoplankton 'Holoplankton' only Habitat=Holoplankton
Meroplankton 'Meroplankton' only Habitat=Meroplankton
Large copepods 'Copepod' AND 'Lg' Copepod=Y AND ZooSize=Lg
Small copepods 'Copepod' AND 'Sm' Copepod=Y AND ZooSize=Sm
Ciliates 'Ciliate' only PhytoplanktonType=Ciliate
Microflagellates 'Dinoflagellate' AND 'Sm' PhytoplanktonType=Dinoflagellate AND PhytoplanktonSize=Sm
Table c: Matrix used to determine the confidence in each lifeform. Only pairs with high confidence were used for this assessment
Easy to ID/speciate Difficult to ID/speciate
Known traits High Low
Unknown traits Low Low

Results

This assessment reveals change in North-East Atlantic plankton communities between the periods 2004–2008 and 2009–2014. Although lifeform pairs exhibited significant change in some areas of the North-East Atlantic (Figure 3), this does not necessarily imply deterioration of environmental conditions.

 

Figure 3: Change in Plankton Index between the periods 2004–2008 and 2009–2014 for each lifeform pair

Darker blue indicates a more pronounced change. Grey shading represents where there were not enough / well-represented data to determine a Plankton Index. Cells with dots indicate significant change (p<0.01) since the period 2004–2008. Changes in the Plankton Index do not necessarily indicate a deterioration of environmental conditions. They do, however, indicate change from starting conditions.The x-axis represents the ecohydrodynamic zones described in Figure 2.

The ‘holoplankton and meroplankton’ lifeform pair experienced significant change in most areas, suggesting changes in linkage between the benthic and pelagic components of the ecosystem. Changes have also occurred in the ‘small copepod and large copepod’ lifeform pair in many areas, which could indicate possible alterations to food web structure and energy flows. The ‘pelagic and tychopelagic’ (benthic) diatom lifeform pair only underwent significant change in a few areas, which could indicate no important changes in resuspension events in much of the North-East Atlantic.

It is not currently possible to link the changes in lifeform pairs to any particular human pressures, or to link these changes to other biodiversity Indicator Assessments.

The methods and data for this indicator are considered to be of moderate confidence.

Figure 3 shows significant change in many lifeform pairs in many ecohydrodynamic (EHDs) zones between the starting conditions (2004–2008) and the 2009–2014 period. Even though the changes observed are statistically significant, this does not suggest that the changes are symptomatic of a deteriorating environment. It is widely acknowledged that prevailing conditions are the key driver of plankton communities in the North-East Atlantic, particularly at the relatively large EHD scale (see van Leeuwen et al., 2015, 2016).

Progress towards further understanding and interpretation of the results was made during the European Union co-financed EcApRHA project coordinated by OSPAR to address gaps in biodiversity indicator development.

All extended results can be downloaded here.

Confidence Assessment

The plankton lifeform indicator that has been developed and used in this assessment has been published in the peer-reviewed literature and each dataset used has internal quality assurance procedures, however this is its first application to zooplankton data.  This assessment has been undertaken using data collected over a period of 57 years, with >200,000 data collection points across the Greater North Sea, Celtic Seas, Bay of Biscay and Iberian Coast. In-situ data from additional Contracting Parties need to be included. The data and methods used here are therefore considered robust and of moderate confidence.

Conclusion

This assessment indicates that there is variability in the plankton community for all lifeforms, which is in accordance with the published scientific literature on plankton dynamics.

While the indicator assessment shows there is change, further work is needed to draw conclusions on the magnitude, direction and the key pressures or environmental factors driving change in lifeform pairs. Interpretation of the results and further refinement of the methodology have still to take place. An extensive peer-reviewed research base, however, suggests that prevailing oceanographic and climatic conditions are the overall driver of plankton change in the North-East Atlantic.

This assessment indicates spatial and temporal variability in the plankton community for all lifeforms. This is in accordance with expert knowledge on plankton dynamics from the scientific literature (McQuatters-Gollop et al., 2015). The plankton lifeform indicator has been published in the peer-reviewed literature (Tett et al., 2008) and has been developed at the sub-regional scale using data from the pan-European Continuous Plankton Recorder survey which has documented quality assurance procedures and has supported over 1000 scientific publications (McQuatters-Gollop et al., 2015). The data and methods used here are therefore considered robust and of high confidence.

Conclusions on the results cannot yet be made as additional work is needed to investigate drivers of change. Interpretation of the results and further refinement to the methodology will occur at a later time but is resource dependent.

Knowledge Gaps

Further scientific research is needed to examine the magnitude and direction of change in the Plankton Index with respect to each lifeform pair, as well as the ecological consequences of such change, for each lifeform pair in each ecohydrodynamic zone. It is also necessary to investigate links between change in lifeform pairs and human and climatic pressures. If changes due to prevailing conditions (such as natural variability and climate change) can be separated from those caused by human pressures in each region, this will help to inform management decision-making by allowing the application of regionally-targeted management measures only where needed.

The following work needs to be taken forward to address the knowledge gaps and improve the assessment provided by this indicator:

  • Determination of the appropriate time period to represent starting conditions;
  • Identification and testing of possible assessment value options;
  • Refinement of spatial scale of ecohydrodynamic zones in the Bay of Biscay and Iberian Coast;
  • Investigation into appropriate Harmful Algal Blooms (HABs) lifeform pair;
  • Incorporation of additional plankton datasets as data and resources become available (including other OSPAR regions, especially to address EU Member States’ MSFD reporting requirements);
  • Incorporate historical data, and based on these potentially revisit the methodology around setting envelopes;
  • Examination of drivers of change in plankton lifeforms in each area, which will allow apportionment of change to climate and anthropogenic pressures;
  • Determination of the traits for many plankton species – for some species even information about basic biological characteristics, such as diet, is not yet known;
  • Understanding the role of pico and nanoplankton (small plankton composed by cells with sizes between respectively 0.2 μm and 2 μm and 2 μm and 20 µm) in the ecosystem. These small size categories are difficult to measure routinely and are thus mainly ignored. However, they make up a significant proportion of the plankton biomass relevant to marine food webs. Their ecosystem role needs further investigation so that they can be included in an appropriate new or existing lifeform; and
  • Establishment of monitoring for pico and nanoplakton to provide data on a relatively data-sparse component of the plankton community.

Beaugrand, G., (2005). Monitoring pelagic ecosystems using plankton indicators. ICES Journal of Marine Science, 62: 333-338.

Bremner, J., Frid, C.L.J. and Rogers, S.I., (2004). Biological traits of the North Sea benthos – Does fishing affect benthic ecosystem function? . In: P. Barnes, Thomas, J. (Editor), Benthic habitats and the effects of fishing. American Fisheries Society, Bethesda, MD, pp. 477-489.

Charvet, S., Roger, M.C., Faessel, B. and Lafont, M., (1998). Biomonitoring of freshwater ecosystems by the use of biological traits. Annals of Limnology – International Journal of Limnology, 34: 455-464.

de Jonge, V.N., (2007). Toward the application of ecological concepts in EU coastal water management. Marine Pollution Bulletin, 55: 407-414.

Estrada, M., Henriksen, P., Gasol, J.M., Casamayor, E.O. and Pedrós-Alió, C., (2004). Diversity of Planktonic Photoautotrophic Microorganisms Along a Salinity Gradient as Depicted by Microscopy, Flow Cytometry, Pigment Analysis and DNA-based Methods. FEMS Microbiology Ecology, 49: 281-293.

Gallego, I., Davidson, T.A., Jeppesen, E., Perez-Martinez, C., Sanchez-Castillo, P., Juan, M., Fuentes-Rodriguez, F., Leon, D., Penalver, P., Toja, J. and Casas, J.J., (2012). Taxonomic or ecological approaches? Searching for phytoplankton surrogates in the determination of richness and assemblage composition in ponds. Ecological Indicators, 18: 575-585.

Garmendia, M., Borja, Á., Franco, J. and Revilla, M., (2012). Phytoplankton composition indicators for the assessment of eutrophication in marine waters: Present state and challenges within the European directives. Marine Pollution Bulletin, 66: 7-16.

Gowen, R.J., McQuatters-Gollop, A., Tett, P., Best, M., Bresnan, E., Castellani, C., Cook, K., Forster, R., Scherer, C. and Mckinney, A., (2011). The Development of UK Pelagic (Plankton) Indicators and Targets for the MSFD. Advice to Defra, Belfast, UK, 41 pp.

HELCOM, (2012). Development of the HELCOM core-set indicators Part B, GES 8/2012/7b. HELCOM, Brussels,

McQuatters-Gollop, A., Edwards, M., Helaouët, P., Johns, D.G., Owens, N.J.P., Raitsos, D.E., Schroeder, D., Skinner, J. and Stern, R.F., (2015). The Continuous Plankton Recorder survey: how can long-term phytoplankton datasets deliver Good Environmental Status? . Estuarine, Coastal and Shelf Science, 162: 88-97.

Mouillot, D., Spatharis, S., Reizopoulou, S., Laugier, T., Sabetta, L., Basset, A. and Chi, T.D., (2006). Alternatives to taxonomic-based approaches to assess changes in transitional water communities. Aquatic Conservation-Marine and Freshwater Ecosystems, 16: 469-482.

Rakocinski, C.F., (2012). Evaluating macrobenthic process indicators in relation to organic enrichment and hypoxia. Ecological Indicators, 13: 1-12.

Rombouts, I., Beaugrand, G., Fizzala, X., Gaill, F., Greenstreet, S.P.R., Lamare, S., Le Loc’h, F., McQuatters-Gollop, A., Mialet, B., Niquil, N., Percelay, J., Renaud, F., Rossberg, A.G. and Féral, J.P., (2013). Food web indicators under the Marine Strategy Framework Directive: from complexity to simplicity? Ecological Indicators, 29: 246–254.

Tett, P., Carreira, C., Mills, D.K., van Leeuwen, S., Foden, J., Bresnan, E. and Gowen, R.J., (2008). Use of a phytoplankton community index to assess the health of coastal waters. ICES Journal of Marine Science, 65: 1475-1482.

Tett, P., Gowen, R., Painting, S., Elliott, M., Forster, R., Mills, D., Bresnan, E., Capuzzo, E., Fernandes, T. and Foden, J. (2013) Framework for understanding marine ecosystem health. Marine Ecology Progress Series, 494, 1-27. https://doi.org/10.3354/meps10539

van Leeuwen, S., Tett, P., Mills, D. and Molen, J.v.d., (2015). Stratified and nonstratified areas in the North Sea: Long-term variability and biological and policy implications. Journal of Geophysical Research Oceans, 120: 4670-4686.

van Leeuwen, S.M., Le Quesne, W.F. and Parker, E.R., (2016). Potential future fisheries yields in shelf waters: a model study of the effects of climate change and ocean acidification. Biogeosciences, 13: 441-454.