Satellite discrimination of Karenia mikimotoi and Phaeocystis harmful algal blooms in European coastal waters: Merged classification of ocean colour data

The detection of dense harmful algal blooms (HABs) by satellite remote sensing is usually based on analysis of chlorophyll-a as a proxy. However, this approach does not provide information about the potential harm of bloom, nor can it identify the dominant species. The developed HAB risk classification method employs a fully automatic data-driven approach to identify key characteristics of water leaving radiances and derived quantities, and to classify pixels into “harmful”, “non-harmful” and “no bloom” categories using Linear Discriminant Analysis (LDA). Discrimination accuracy is increased through the use of spectral ratios of water leaving radiances, absorption and backscattering. To reduce the false alarm rate the data that cannot be reliably classified are automatically labelled as “unknown”. This method can be trained on different HAB species or extended to new sensors and then applied to generate independent HAB risk maps; these can be fused with other sensors to fill gaps or improve spatial or temporal resolution. The HAB discrimination technique has obtained accurate results on MODIS and MERIS data, correctly identifying 89% of Phaeocystis globosa HABs in the southern North Sea and 88% of Karenia mikimotoi blooms in the Western English Channel. A linear transformation of the ocean colour discriminants is used to estimate harmful cell counts, demonstrating greater accuracy than if based on chlorophyll-a; this will facilitate its integration into a HAB early warning system operating in the southern North Sea.

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.

Monitoring of a quasi-stationary eddy in the Bay of Biscay by means of satellite, in situ and model results

The presence of a quasi-stationary anticyclonic eddy within the southeastern Bay of Biscay (centred around 44°30′N-4°W) has been reported on various occasions in the bibliography. The analysis made in this study for the period 2003–2010, by using in situ and remote sensing measurements and model results shows that this mesoscale coherent structure is present almost every year from the end of winter-beginning of spring, to the beginning of fall. During this period it remains in an area limited to the east by the Landes Plateau, to the west by Le Danois Bank and Torrelavega canyon and to the northwest by the Jovellanos seamount. All the observations and analysis made in this contribution, suggest that this structure is generated between Capbreton and Torrelavega canyons. Detailed monitoring from in situ and remote sensing data of an anticyclonic quasi-stationary eddy, in 2008, shows the origin of this structure from a warm water current located around 43°42′N-3°30′W in mid-January. This coherent structure is monitored until August around the same area, where it has a marked influence on the Sea Level Anomaly, Sea Surface Temperature and surface Chlorophyll-a concentration. An eddy tracking method, applied to the outputs of a numerical model, shows that the model is able to reproduce this type of eddy, with similar 2D characteristics and lifetimes to that suggested by the observations and previous works. This is the case, for instance, of the simulated MAY04 eddy, which was generated in May 2004 around Torrelavega canyon and remained quasi-stationary in the area for 4 months. The diameter of this eddy ranged from 40 to 60 km, its azimuthal velocity was less than 20 cm s−1, its vertical extension reached 3000–3500 m depth during April and May and it was observed to interact with other coherent structures.

On the temporal consistency of chlorophyll products derived from three ocean-colour sensors

Satellite ocean-colour sensors have life spans lasting typically five-to-ten years. Detection of long-term trends in chlorophyll-a concentration (Chl-a) using satellite ocean colour thus requires the combination of different ocean-colour missions with sufficient overlap to allow for cross-calibration. A further requirement is that the different sensors perform at a sufficient standard to capture seasonal and inter-annual fluctuations in ocean colour. For over eight years, the SeaWiFS, MODIS-Aqua and MERIS ocean-colour sensors operated in parallel. In this paper, we evaluate the temporal consistency in the monthly Chl-a time-series and in monthly inter-annual variations in Chl-a among these three sensors over the 2002–2010 time period. By subsampling the monthly Chl-a data from the three sensors consistently, we found that the Chl-a time-series and Chl-a anomalies among sensors were significantly correlated for >90% of the global ocean. These correlations were also relatively insensitive to the choice of three Chl-a algorithms and two atmospheric-correction algorithms. Furthermore, on the subsampled time-series, correlations between Chl-a and time, and correlations between Chl-a and physical variables (sea-surface temperature and sea-surface height) were not significantly different for >92% of the global ocean. The correlations in Chl-a and physical variables observed for all three sensors also reflect previous theories on coupling between physical processes and phytoplankton biomass. The results support the combining of Chl-a data from SeaWiFS, MODIS-Aqua and MERIS sensors, for use in long-term Chl-a trend analysis, and highlight the importance of accounting for differences in spatial sampling among sensors when combining ocean-colour observations.

Evaluating the SEVIRI Fire Thermal Anomaly Detection Algorithm across the Central African Republic Using the MODIS Active Fire Product

Satellite-based remote sensing of active fires is the only practical way to consistently and continuously monitor diurnal fluctuations in biomass burning from regional, to continental, to global scales. Failure to understand, quantify, and communicate the performance of an active fire detection algorithm, however, can lead to improper interpretations of the spatiotemporal distribution of biomass burning, and flawed estimates of fuel consumption and trace gas and aerosol emissions. This work evaluates the performance of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) Fire Thermal Anomaly (FTA) detection algorithm using seven months of active fire pixels detected by the Moderate Resolution Imaging Spectroradiometer (MODIS) across the Central African Republic (CAR). Results indicate that the omission rate of the SEVIRI FTA detection algorithm relative to MODIS varies spatially across the CAR, ranging from 25% in the south to 74% in the east. In the absence of confounding artifacts such as sunglint, uncertainties in the background thermal characterization, and cloud cover, the regional variation in SEVIRI's omission rate can be attributed to a coupling between SEVIRI's low spatial resolution detection bias (i.e., the inability to detect fires below a certain size and intensity) and a strong geographic gradient in active fire characteristics across the CAR. SEVIRI's commission rate relative to MODIS increases from 9% when evaluated near MODIS nadir to 53% near the MODIS scene edges, indicating that SEVIRI errors of commission at the MODIS scene edges may not be false alarms but rather true fires that MODIS failed to detect as a result of larger pixel sizes at extreme MODIS scan angles. Results from this work are expected to facilitate (i) future improvements to the SEVIRI FTA detection algorithm; (ii) the assimilation of the SEVIRI and MODIS active fire products; and (iii) the potential inclusion of SEVIRI into a network of geostationary sensors designed to achieve global diurnal active fire monitoring.

On the Front Line: frontal zones as priority at-sea conservation areas for mobile marine vertebrates

1.Identifying priority areas for marine vertebrate conservation is complex because species of conservation concern are highly mobile, inhabit dynamic habitats and are difficult to monitor. 2.Many marine vertebrates are known to associate with oceanographic fronts – physical interfaces at the transition between water masses – for foraging and migration, making them important candidate sites for conservation. Here, we review associations between marine vertebrates and fronts and how they vary with scale, regional oceanography and foraging ecology. 3.Accessibility, spatiotemporal predictability and relative productivity of front-associated foraging habitats are key aspects of their ecological importance. Predictable mesoscale (10s–100s km) regions of persistent frontal activity (‘frontal zones’) are particularly significant. 4.Frontal zones are hotspots of overlap between critical habitat and spatially explicit anthropogenic threats, such as the concentration of fisheries activity. As such, they represent tractable conservation units, in which to target measures for threat mitigation. 5.Front mapping via Earth observation (EO) remote sensing facilitates identification and monitoring of these hotspots of vulnerability. Seasonal or climatological products can locate biophysical hotspots, while near-real-time front mapping augments the suite of tools supporting spatially dynamic ocean management. 6.Synthesis and applications. Frontal zones are ecologically important for mobile marine vertebrates. We surmise that relative accessibility, predictability and productivity are key biophysical characteristics of ecologically significant frontal zones in contrasting oceanographic regions. Persistent frontal zones are potential priority conservation areas for multiple marine vertebrate taxa and are easily identifiable through front mapping via EO remote sensing. These insights are useful for marine spatial planning and marine biodiversity conservation, both within Exclusive Economic Zones and in the open oceans.

Assimilation of remotely-sensed optical properties to improve marine biogeochemistry modelling

In this paper we evaluate whether the assimilation of remotely-sensed optical data into a marine ecosystem model improves the simulation of biogeochemistry in a shelf sea. A localized Ensemble Kalman filter was used to assimilate weekly diffuse light attenuation coefficient data, Kd(443) from SeaWiFs, into an ecosystem model of the western English Channel. The spatial distributions of (unassimilated) surface chlorophyll from satellite, and a multivariate time series of eighteen biogeochemical and optical variables measured in situ at one long-term monitoring site were used to evaluate the system performance for the year 2006. Assimilation reduced the root mean square error and improved the correlation with the assimilated Kd(443) observations, for both the analysis and, to a lesser extent, the forecast estimates, when compared to the reference model simulation. Improvements in the simulation of (unassimilated) ocean colour chlorophyll were less evident, and in some parts of the Channel the simulation of this data deteriorated. The estimation errors for the (unassimilated) in situ data were reduced for most variables with some exceptions, e.g. dissolved nitrogen. Importantly, the assimilation adjusted the balance of ecosystem processes by shifting the simulated food web towards the microbial loop, thus improving the estimation of some properties, e.g. total particulate carbon. Assimilation of Kd(443) outperformed a comparative chlorophyll assimilation experiment, in both the estimation of ocean colour data and in the simulation of independent in situ data. These results are related to relatively low error in Kd(443) data, and because it is a bulk optical property of marine ecosystems. Assimilation of remotely-sensed optical properties is a promising approach to improve the simulation of biogeochemical and optical variables that are relevant for ecosystem functioning and climate change studies.

How many Coccolithovirus genotypes does it take to terminate an Emiliania huxleyi bloom?

Giant viruses are known to be significant mortality agents of phytoplankton, often being implicated in the terminations of large Emiliania huxleyi blooms. We have previously shown the high temporal variability of E. huxleyi-infecting coccolithoviruses (EhVs) within a Norwegian fjord mesocosm. In the current study we investigated EhV dynamics within a naturally-occurring E. huxleyi bloom in the Western English Channel. Using denaturing gradient gel electrophoresis and marker gene sequencing, we uncovered a spatially highly dynamic Coccolithovirus population that was associated with a genetically stable E. huxleyi population as revealed by the major capsid protein gene (mcp) and coccolith morphology motif (CMM), respectively. Coccolithoviruses within the bloom were found to be variable with depth and unique virus populations were detected at different stations sampled indicating a complex network of EhV-host infections. This ultimately will have significant implications to the internal structure and longevity of ecologically important E. huxleyi blooms.

Comparison of new and primary production models using SeaWiFS data in contrasting hydrographic zones of the northern North Atlantic.

The accuracy of two satellite models of marine primary (PP) and new production (NP) were assessed against 14C and 15N uptake measurements taken during six research cruises in the northern North Atlantic. The wavelength resolving model (WRM) was more accurate than the Vertical General Production Model (VGPM) for computation of both PP and NP. Mean monthly satellite maps of PP and NP for both models were generated from 1997 to 2010 using SeaWiFS data for the Irminger basin and North Atlantic. Intra- and inter-annual variability of the two models was compared in six hydrographic zones. Both models exhibited similar spatio-temporal patterns: PP and NP increased from April to June and decreased by August. Higher values were associated with the East Greenland Current (EGC), Iceland Basin (ICB) and the Reykjanes Ridge (RKR) and lower values occurred in the Central Irminger Current (CIC), North Irminger Current (NIC) and Southern Irminger Current (SIC). The annual PP and NP over the SeaWiFS record was 258 and 82 gC m-2 yr-1 respectively for the VGPM and 190 and 41 gC m-2 yr-1 for the WRM. Average annual cumulative sum in the anomalies of NP for the VGPM were positively correlated with the North Atlantic Oscillation (NAO) in the EGC, CIC and SIC and negatively correlated with the multivariate ENSO index (MEI) in the ICB. By contrast, cumulative sum of the anomalies of NP for the WRM were significantly correlated with NAO only in the EGC and CIC. NP from both VGPM and WRM exhibited significant negative correlations with Arctic Oscillation (AO) in all hydrographic zones. The differences in estimates of PP and NP in these hydrographic zones arise principally from the parameterisation of the euphotic depth and the SST dependence of photo-physiological term in the VGPM, which has a greater sensitivity to variations in temperature than the WRM. In waters of 0 to 5C PP using the VGPM was 43% higher than WRM, from 5 to 10C the VGPM was 29% higher and from 10 to 15C the VGPM was 27% higher.

Effect of phytoplankton size classes on bio-optical properties of phytoplankton in the Western Iberian coast: Application of models

Chlorophyll-a satellite products are routinely used in oceanography, providing a synoptic and global view of phytoplankton abundance. However, these products lack information on the community structure of the phytoplankton, which is crucial for ecological modelling and ecosystem studies. To assess the usefulness of existing methods to differentiate phytoplankton functional types (PFT) or phytoplankton size classes from satellite data, in-situ phytoplankton samples collected in the Western Iberian coast, on the North-East Atlantic, were analysed for pigments and absorption spectra. Water samples were collected in five different locations, four of which were located near the shore and another in an open-ocean, seamount region. Three different modelling approaches for deriving phytoplankton size classes were applied to the in situ data. Approaches tested provide phytoplankton size class information based on the input of pigments data (Brewin et al., 2010), absorption spectra data (Ciotti et al., 2002) or both (Uitz et al., 2008). Following Uitz et al. (2008), results revealed high variability in microphytoplankton chlorophyll-specific absorption coefficients, ranging from 0.01 to 0.09 m2 (mg chl)− 1 between 400 and 500 nm. This spectral analysis suggested, in one of the regions, the existence of small cells (< 20 μm) in the fraction of phytoplankton presumed to be microphytoplankton (based on diagnostic pigments). Ciotti et al. (2002) approach yielded the highest differences between modelled and measured absorption spectra for the locations where samples had high variability in community structure and cell size. The Brewin et al. (2010) pigment-based model was adjusted and a set of model coefficients are presented and recommended for future studies in offshore water of the Western Iberian coast.

Challenges and difficulties in assessing the Environmental Status under the requirements of the Ecosystem Approach in North-African countries, illustrated by eutrophication assessment

Marine ecosystems provide many ecosystem goods and services. However, these ecosystems and the benefits they create for humans are subject to competing uses and increasing pressures. As a consequence of the increasing threats to the marine environment, several regulations require applying an ecosystem-based approach for managing the marine environment. Within the Mediterranean Sea, in 2008, the Contracting Parties of the Mediterranean Action Plan decided to progressively apply the Ecosystem Approach (EcAp) with the objective of achieving Good Environmental Status (GES) for 2018. To assess the Environmental Status, the EcAp proposes 11 Ecological Objectives, each of which requires a set of relevant indicators to be integrated. Progress towards the EcAp entails a gradual and important challenge for North-African countries, and efforts have to be initiated to propose and discuss methods. Accordingly, to enhance the capacity of North-African countries to implement EcAp and particularly to propose and discuss indicators and methods to assess GES, the aim of this manuscript is to identify the practical problems and gaps found at each stage of the Environmental Status assessment process. For this purpose, a stepwise method has been proposed to assess the Environmental Status using Ecologic Objective 5-Eutrophication as example.

Basking sharks and oceanographic fronts: quantifying associations in the north‐east Atlantic

Understanding the mechanisms linking oceanographic processes and marine vertebrate habitat use is critical to effective management of populations of conservation concern. The basking shark Cetorhinus maximus has been shown to associate with oceanographic fronts – physical interfaces at the transitions between water masses – to exploit foraging opportunities resulting from aggregation of zooplankton. However, the scale, significance and variability of these observed associations have not yet been established. Here, we quantify the influence of mesoscale (10s – 100s km) frontal activity on habitat use over timescales of weeks to months. We use animal-mounted archival tracking with composite front mapping via Earth Observation (EO) remote sensing to provide an oceanographic context to individual shark movements. We investigate levels of association with fronts occurring over two spatio-temporal scales, (i) broad-scale seasonally persistent frontal zones and (ii) contemporaneous mesoscale thermal and chl-a fronts. Using random walk simulations and logistic regression within an iterative generalised linear mixed modelling (GLMM) framework, we find that seasonal front frequency is a significant predictor of shark presence. Temporally-matched oceanographic metrics also indicate that sharks demonstrate a preference for productive regions, and associate with contemporaneous thermal and chl-a fronts more frequently than could be expected at random. Moreover, we highlight the importance of cross-frontal temperature change and persistence, which appear to interact to affect the degree of prey aggregation along thermal fronts. These insights have clear implications for understanding the preferred habitats of basking sharks in the context of anthropogenic threat management and marine spatial planning in the northeast Atlantic.

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.

Seasonality and interannual variability of the European Slope Current from 20 years of altimeter data compared with in situ measurements

The European Slope Current (SC) is a major section of the warm poleward flow from the Atlantic to the Arctic, which also moderates the exchange of heat, salt, nutrients and carbon between the deep ocean and the European shelf seas. The mean structure of the geostrophic flow, seasonality, interannual variability and long-term trend of SC are appraised with an unprecedented continuous 20-year satellite altimeter dataset. Comparisons with long term in situ data showed a maximum correlation of r2=0.51 between altimeter and Acoustic Doppler Current Profilers (ADCP), with similar results for drogued buoy data. Mean geostrophic currents were appraised more comprehensively than previous attempts, and the paths of 4 branches of the North Atlantic Current (NAC) and positions of 5 eddies in the region were derived quantitatively. A consistent seasonal cycle in the flow of the SC was found at all 8 sections along the European shelf slope, with maximum poleward flow in the winter and minimum in the summer. The seasonal difference in the altimetry current speed amounted to ~8-10 cm s-1 at the northern sections, but only ~5 cm s-1 on the Bay of Biscay slopes. This extended altimeter dataset indicates significant regional and seasonal variations, and has revealed new insights into the interannual variability of the SC. It is shown that there is a peak poleward flow at most positions along a ~2000 km stretch of the continental slope from Portugal to Scotland during 1995-1997, but this did not clearly relate to the extreme negative North Atlantic Oscillation (NAO) in the winter of 1995-1996. The speed of the SC exhibited a long term decreasing trend of ~1% per year. By contrast the NAC showed no significant trend over the 20-year period. Major changes in the NAC occurred three times, and these changes followed decreases in the NAO index.

Basking sharks and oceanographic fronts: quantifying associations in the north-east Atlantic

Understanding the mechanisms linking oceanographic processes and marine vertebrate habitat use is critical to effective management of populations of conservation concern. The basking shark Cetorhinus maximus has been shown to associate with oceanographic fronts – physical interfaces at the transitions between water masses – to exploit foraging opportunities resulting from aggregation of zooplankton. However, the scale, significance and variability of these observed associations have not yet been established. Here, we quantify the influence of mesoscale (10s – 100s km) frontal activity on habitat use over timescales of weeks to months. We use animal-mounted archival tracking with composite front mapping via Earth Observation (EO) remote sensing to provide an oceanographic context to individual shark movements. We investigate levels of association with fronts occurring over two spatio-temporal scales, (i) broad-scale seasonally persistent frontal zones and (ii) contemporaneous mesoscale thermal and chl-a fronts. Using random walk simulations and logistic regression within an iterative generalised linear mixed modelling (GLMM) framework, we find that seasonal front frequency is a significant predictor of shark presence. Temporally-matched oceanographic metrics also indicate that sharks demonstrate a preference for productive regions, and associate with contemporaneous thermal and chl-a fronts more frequently than could be expected at random. Moreover, we highlight the importance of cross-frontal temperature change and persistence, which appear to interact to affect the degree of prey aggregation along thermal fronts. These insights have clear implications for understanding the preferred habitats of basking sharks in the context of anthropogenic threat management and marine spatial planning in the northeast Atlantic.