Groom et al. Satellite Ocean Colour to end-users. Example applications of the ocean-colour data are presented, focusing on the climate data record and operational applications including water quality and assimilation into numerical models. Current capacity building and training activities pertinent to ocean colour are described and finally a summary of future perspectives is provided.

This paper proposes the use of assimilation of phytoplankton functional types (PFTs) surface chlorophyll for operational forecasting of biogeochemistry on the North‐West European (NWE) Shelf. We explicitly compare the 5‐day forecasting skill of three runs of a physical‐biogeochemical model: (a) a free reference run, (b) a run with daily data assimilation (DA) of total surface chlorophyll (ChlTot), and (c) a run with daily PFTs DA. We show that small total chlorophyll model bias hides comparatively large biases in PFTs chlorophyll, which ChlTot DA fails to correct. This is because the ChlTot DA splits the assimilated total chlorophyll into PFTs by preserving their simulated ratios, rather than taking account of the observed PFT concentrations. Unlike ChlTot DA, PFTs DA substantially improves model representation of PFTs chlorophyll. During forecasting the DA reanalysis skill in representing PFTs chlorophyll degrades toward the free run skill; however, PFTs DA outperforms free run within the whole 5‐day forecasting period. We validated our results with in situ data, and we demonstrated that (in both DA cases) the DA substantially improves the model representation of CO2 fugacity (PFTs DA more than ChlTot DA). ChlTot DA has a positive impact on the representation of silicate, while the PFTs DA seems to have a negative impact. The impact of DA on nitrate and phosphate is not significant. The implications of using a univariate assimilation method, which preserves the phytoplankton stochiometry, and the impact of model biases on the nonassimilated variables are discussed.

We assimilated phytoplankton functional types (PFTs) derived from ocean color into a marine ecosystem model, to improve the simulation of biogeochemical indicators and emerging properties in a shelf sea. Error‐characterized chlorophyll concentrations of four PFTs (diatoms, dinoflagellates, nanoplankton, and picoplankton), as well as total chlorophyll for comparison, were assimilated into a physical‐biogeochemical model of the North East Atlantic, applying a localized Ensemble Kalman filter. The reanalysis simulations spanned the years 1998–2003. The skill of the reference and reanalysis simulations in estimating ocean color and in situ biogeochemical data were compared by using robust statistics. The reanalysis outperformed both the reference and the assimilation of total chlorophyll in estimating the ocean‐color PFTs (except nanoplankton), as well as the not‐assimilated total chlorophyll, leading the model to simulate better the plankton community structure. Crucially, the reanalysis improved the estimates of not‐assimilated in situ data of PFTs, as well as of phosphate and pCO2, impacting the simulation of the air‐sea carbon flux. However, the reanalysis increased further the model overestimation of nitrate, in spite of increases in plankton nitrate uptake. The method proposed here is easily adaptable for use with other ecosystem models that simulate PFTs, for, e.g., reanalysis of carbon fluxes in the global ocean and for operational forecasts of biogeochemical indicators in shelf‐sea ecosystems.

Understanding ecosystem state on the North‐West European (NWE) Shelf is of major importance for both economy and climate research. The purpose of this work is to advance our modeling of in‐water optics on the NWE Shelf, with important implications for how we model primary productivity, as well as for assimilation of water‐leaving radiances. We implement a stand‐alone bio‐optical module into the existing coupled physical‐biogeochemical model configuration. The advantage of the bio‐optical module, when compared to the preexisting light scheme is that it resolves the underwater light spectrally and distinguishes between direct and diffuse downwelling streams. The changed underwater light compares better with both satellite and in situ observations. The module lowered the underwater photosynthetically active radiation, decreasing the simulated primary productivity, but overall, the improved underwater light had relatively limited impact on the phytoplankton seasonal dynamics. We showed that the model skill in representing phytoplankton seasonal cycle (e.g., phytoplankton bloom) can be substantially improved either by assimilation of satellite phytoplankton functional type (PFT) chlorophyll, or by assimilating a novel PFT absorption product. Assimilation of the two PFT products yields similar results, with an important difference in the PFT community structure. Both assimilative runs lead to lower plankton biomass and increase the nutrient concentrations. We discuss some future directions on how to improve our model skill in biogeochemistry without using assimilation, for example, by improving nutrient forcing, retuning the model parameters, and using the bio‐optical module to provide a two‐way physical‐biogeochemical coupling, improving the consistency between model physical and biogeochemical components.

We present a model that estimates the spectral phytoplankton absorption coefficient (aph(λ)) of four phytoplankton groups (picophytoplankton, nanophytoplankton, dinoflagellates, and diatoms) as a function of the total chlorophyll-a concentration (C) and sea surface temperature (SST). Concurrent data on aph(λ) (at 12 visible wavelengths), C and SST, from the surface layer (<20 m depth) of the North Atlantic Ocean, were partitioned into training and independent validation data, the validation data being matched with satellite ocean-colour observations. Model parameters (the chlorophyll-specific phytoplankton absorption coefficients of the four groups) were tuned using the training data and found to compare favourably (in magnitude and shape) with results of earlier studies. Using the independent validation data, the new model was found to retrieve total aph(λ) with a similar performance to two earlier models, using either in situ or satellite data as input. Although more complex, the new model has the advantage of being able to determine aph(λ) for four phytoplankton groups and of incorporating the influence of SST on the composition of the four groups. We integrate the new four-population absorption model into a simple model of ocean colour, to illustrate the influence of changes in SST on phytoplankton community structure, and consequently, the blue-to-green ratio of remote-sensing reflectance. We also present a method of propagating error through the model and illustrate the technique by mapping errors in group-specific aph(λ) using a satellite image. We envisage the model will be useful for ecosystem model validation and assimilation exercises and for investigating the influence of temperature change on ocean colour.

In this work we produced a long‐term reanalysis of the phytoplankton community structure in the Mediterranean Sea and used it to define ecoregions. These were based on the spatial variability of the phytoplankton type fractions and their influence on selected carbon fluxes. A regional ocean color product of four phytoplankton functional types (PFTs; diatoms, dinoflagellates, nanophytoplankton, and picophytoplankton) was assimilated into a coupled physical‐biogeochemical model of the Mediterranean Sea (Proudman Oceanographic Laboratory Coastal Ocean Modelling System‐European Regional Seas Ecosystem Model, POLCOMS–ERSEM) by using a 100‐member ensemble Kalman filter, in a reanalysis simulation for years 1998–2014. The reanalysis outperformed the reference simulation in representing the assimilated ocean color PFT fractions to total chlorophyll, although the skill for the ocean color PFT concentrations was not improved significantly. The reanalysis did not impact noticeably the reference simulation of not assimilated in situ observations, with the exception of a slight bias reduction for the situ PFT concentrations, and a deterioration of the phosphate simulation. We found that the Mediterranean Sea can be subdivided in three PFT‐based ecoregions, derived from the spatial variability of the PFT fraction dominance or relevance. Picophytoplankton dominates the largest part of open ocean waters; microphytoplankton dominates in a few, highly productive coastal spots near large‐river mouths; nanophytoplankton is relevant in intermediate‐productive coastal and Atlantic‐influenced waters. The trophic and carbon sedimentation efficiencies are highest in the microphytoplankton ecoregion and lowest in the picophytoplankton and nanophytoplankton ecoregions. The reanalysis and regionalization offer new perspectives on the variability of the structure and functioning of the phytoplankton community and related biogeochemical fluxes, with foreseeable applications in Blue Growth of the Mediterranean Sea.

Characterization of chlorophyll and sea surface temperature (SST) structural heterogeneity using their scaling properties can provide a useful tool to estimate the relative importance of key physical and biological drivers. Seasonal, annual, and also instantaneous spatial distributions of chlorophyll and SST, determined from satellite measurements, in seven different coastal and shelf‐sea regions around the UK have been studied. It is shown that multifractals provide a very good approximation to the scaling properties of the data: in fact, the multifractal scaling function is well approximated by universal multifractal theory. The consequence is that all of the statistical information about data structure can be reduced to being described by two parameters. It is further shown that also bathymetry scales in the studied regions as multifractal. The SST and chlorophyll multifractal structures are then explained as an effect of bathymetry and turbulence.

In this study we demonstrate the many strengths of scale analysis: we use it to evaluate the Nucleus for European Modelling of the Ocean model skill in representing sea surface temperature (SST) in the Southern Ocean by comparing three model resolutions: 1/12°, 1/4°, and 1°. We show that while 4–5 times resolution scale is sufficient for each model resolution to reproduce the magnitude of satellite Earth Observation (EO) SST spatial variability to within ±10%, the representation of ∼100‐km SST variability patterns is substantially (e.g., ∼50% at 750 km) improved by increasing model resolution from 1° to 1/12°. We also analyzed the dominant scales of the SST model input drivers (short‐wave radiation, air‐sea heat fluxes, wind stress components, wind stress curl, and bathymetry) variability with the purpose of determining the optimal SST model input driver resolution. The SST magnitude of variability is shown to scale with two power law regimes separated by a scaling break at ∼200‐km scale. The analysis of the spatial and temporal scales of dominant SST driver impact helps to interpret this scaling break as a separation between two different dynamical regimes: the (relatively) fast SST dynamics below ∼200 km governed by eddies, fronts, Ekman upwelling, and air‐sea heat exchange, while above ∼200 km the SST variability is dominated by long‐term (seasonal and supraseasonal) modes and the SST geography.