Validation of the NEMO-ERSEM operational ecosystem model for the North West European Continental Shelf

This paper details updates to the Met Office's operational coupled hydrodynamic-ecosystem model from the 7 km Medium-Resolution Continental Shelf – POLCOMS-ERSEM (MRCS-PE) system (Siddorn et al., 2007) to the 7 km Atlantic Margin Model NEMO-ERSEM (AMM7-NE) system. We also provide a validation of the ecosystem component of the new operational system. Comparisons have been made between the model variables and available in situ, satellite and climatological data. The AMM7-NE system has also been benchmarked against the MRCS-PE system. The transition to the new AMM7-NE system was successful and it has been running operationally since March 2012 and has been providing products through MyOcean (http://www.myocean.eu.org) since that time. The results presented herein show the AMM7-NE system performs better than the MRCS-PE system with the most improvement in the model nutrient fields. The problem of nutrient accumulation in the MRCS-PE system appears to be solved in the new AMM7-NE system with nutrient fields improved throughout the domain as discussed in Sect. 4. Improvements in model chlorophyll are also seen but are more modest.

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.

Southern Ocean biogeography and taxonomic resolution: what's in the name?

We examined the taxonomic resolution of zooplankton data required to identify ocean basin scale biogeographic zonation in the Southern Ocean. A 2,154 km transect was completed south of Australia. Sea surface temperature (SST) measured at 1 min intervals showed that seven physical zones were sampled. Zooplankton were collected at a spatial resolution of similar to 9.2 km with a continuous plankton recorder, identified to the highest possible taxonomic resolution and enumerated. Zooplankton assemblage similarity between samples was calculated using the Bray-Curtis index for the taxonomic levels of species, genus, family, order and class after first log(10)(x + 1) (LA) and then presence/absence (PA) transformation of abundance data. Although within and between zone sample similarity increased with decreasing taxonomic resolution, for both data transformations, cluster analysis demonstrated that the biogeographic separation of zones remained at all taxonomic levels when using LA data. ANOSIM confirmed this, detecting significant differences in zooplankton assemblage structure between all seven a priori determined physical zones for all taxonomic levels when using the LA data. In the case of the PA data for the complete data set, and both LA and PA data for a crustacean only data set, no significant differences were detected between zooplankton assemblages in the Polar frontal zone (PFZ) and inter-PFZ at any taxonomic level. Loss of information at resolutions below the species level, particularly in the PA data, prevented the separation of some zones. However, the majority of physical zones were biogeographically distinct from species level to class using both LA and PA transformations. Significant relationships between SST and zooplankton community structure, summarised as NMDS scores, at all taxonomic levels, for both LA and PA transformations, and complete and crustacean only data sets, highlighted the biogeographic relevance of low resolution taxonomic data. The retention of biogeographic information in low taxonomic resolution data shows that data sets collected with different taxonomic resolutions may be meaningfully merged for the post hoc generation of Southern Ocean time series.

Quantifying the impact of environmental variables upon catch-per-unit-effort of the Blue Shark Prionace glauca in the Western English Channel

The effect of environmental variables on blue shark Prionace glauca catch per unit effort (CPUE) in a recreational fishery in the western English Channel, between June and September 1998–2011, was quantified using generalized additive models (GAMs). Sea surface temperature (SST) explained 1·4% of GAM deviance, and highest CPUE occurred at 16·7° C, reflecting the optimal thermal preferences of this species. Surface chlorophyll a concentration (CHL) significantly affected CPUE and caused 27·5% of GAM deviance. Additionally, increasing CHL led to rising CPUE, probably due to higher productivity supporting greater prey biomass. The density of shelf-sea tidal mixing fronts explained 5% of GAM deviance, but was non-significant, with increasing front density negatively affecting CPUE. Time-lagged frontal density significantly affected CPUE, however, causing 12·6% of the deviance in a second GAM and displayed a positive correlation. This outcome suggested a delay between the evolution of frontal features and the subsequent accumulation of productivity and attraction of higher trophic level predators, such as P. glauca.

Indications of a climate effect on Mediterranean fisheries

Using the Food and Agriculture Organization’s (FAO) Mediterranean capture fisheries production dataset in conjunction with global and Mediterranean sea surface temperatures, we investigated trends in fisheries landings and landings per unit of effort of commercially important marine organisms, in relation to temperature oscillations. In addition to the overall warming trend, a temperature shift was detected in the Mediterranean Sea in the late 1990s. Fisheries landings fluctuations were examined for the most abundant commercial species (59 species) and showed significant year-to-year correlations with temperature for nearly 60 % of the cases. From these, the majority (~70 %) were negatively related and showed a reduction of 44 % on average. Increasing trends were found, mainly in the landings of species with short life spans, which seem to have benefited from the increase in water temperature. Τhe effect of oceanic warming is apparent in most species or groups of species sharing ecological (e.g. small and medium pelagic, demersal fish) or taxonomic (e.g. cephalopods, crustaceans) traits. A landings-per-unit-of-effort (LPUE) proxy, using data from the seven Mediterranean European Union member states, also showed significant correlation with temperature fluctuations for six out of the eight species examined, indicating the persistence of temperature influence on landings when the fishing effect is accounted for. The speed of response of marine landings to the warming of the Mediterranean Sea possibly shows both the sensitivity and the vulnerable state of the fish stocks and indicates that climate should be examined together with fisheries as a factor shaping stock fluctuations.

Quantifying the impact of environmental variables upon catch-per-unit-effort of the Blue Shark Prionace glauca in the Western English Channel

The effect of environmental variables on blue shark Prionace glauca catch per unit effort (CPUE) in a recreational fishery in the western English Channel, between June and September 1998–2011, was quantified using generalized additive models (GAMs). Sea surface temperature (SST) explained 1·4% of GAM deviance, and highest CPUE occurred at 16·7° C, reflecting the optimal thermal preferences of this species. Surface chlorophyll a concentration (CHL) significantly affected CPUE and caused 27·5% of GAM deviance. Additionally, increasing CHL led to rising CPUE, probably due to higher productivity supporting greater prey biomass. The density of shelf-sea tidal mixing fronts explained 5% of GAM deviance, but was non-significant, with increasing front density negatively affecting CPUE. Time-lagged frontal density significantly affected CPUE, however, causing 12·6% of the deviance in a second GAM and displayed a positive correlation. This outcome suggested a delay between the evolution of frontal features and the subsequent accumulation of productivity and attraction of higher trophic level predators, such as P. glauca.

Operational monitoring and forecasting of bathing water quality through exploiting satellite Earth observation and models: The AlgaRisk demonstration service

Coastal zones and shelf-seas are important for tourism, commercial fishing and aquaculture. As a result the importance of good water quality within these regions to support life is recognised worldwide and a number of international directives for monitoring them now exist. This paper describes the AlgaRisk water quality monitoring demonstration service that was developed and operated for the UK Environment Agency in response to the microbiological monitoring needs within the revised European Union Bathing Waters Directive. The AlgaRisk approach used satellite Earth observation to provide a near-real time monitoring of microbiological water quality and a series of nested operational models (atmospheric and hydrodynamic-ecosystem) provided a forecast capability. For the period of the demonstration service (2008–2013) all monitoring and forecast datasets were processed in near-real time on a daily basis and disseminated through a dedicated web portal, with extracted data automatically emailed to agency staff. Near-real time data processing was achieved using a series of supercomputers and an Open Grid approach. The novel web portal and java-based viewer enabled users to visualise and interrogate current and historical data. The system description, the algorithms employed and example results focussing on a case study of an incidence of the harmful algal bloom Karenia mikimotoi are presented. Recommendations and the potential exploitation of web services for future water quality monitoring services are discussed.

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.

Are Prorocentrum hoffmannianum and Prorocentrum belizeanum (Dinophyceae, Prorocentrales), the same species? An integration of morphological and molecular data.

The taxonomic assignment of Prorocentrum species is based on morphological characteristics; however, morphological variability has been found for several taxa isolated from different geographical regions. In this study, we evaluated species boundaries of Prorocentrum hoffmannianum and Prorocentrum belizeanum based on morphological and molecular data. A detailed morphological analysis was done, concentrating on the periflagellar architecture. Molecular analyses were performed on partial Small Sub-Unit (SSU) rDNA, partial Large Sub-Unit (LSU) rDNA, complete Internal Transcribed Spacer Regions (ITS1-5.8S-ITS2), and partial cytochrome b (cob) sequences. We concatenated the SSU-ITS-LSU fragments and constructed a phylogenetic tree using Bayesian Inference (BI) and maximum likelihood (ML) methods. Morphological analyses indicated that the main characters, such as cell size and number of depressions per valve, normally used to distinguish P. hoffmannianum from P. belizeanum, overlapped. No clear differences were found in the periflagellar area architecture. Prorocentrum hoffmannianum and P. belizeanum were a highly supported monophyletic clade separated into three subclades, which broadly corresponded to the sample collection regions. Subtle morphological overlaps found in cell shape, size, and ornamentation lead us to conclude that P. hoffmanianum and P. belizeanum might be considered conspecific. The molecular data analyses did not separate P. hoffmannianum and P. belizeanum into two morphospecies, and thus, we considered them to be the P. hoffmannianum species complex because their clades are separated by their geographic origin. These geographic and genetically distinct clades could be referred to as ribotypes: (A) Belize, (B) Florida-Cuba, (C1) India, and (C2) Australia.