Paired measurements of phytoplankton absorption and backscatter, the inherent optical properties central to the interpretation of ocean colour remote sensing data, are notoriously rare. We present a dataset of Chlorophyll a (Chl a) -specific phytoplankton absorption, scatter and backscatter for 17 different phytoplankton groups, derived from first principles using measured in vivo pigment absorption and a well-validated semi-analytical coated sphere model which simulates the full suite of biophysically consistent phytoplankton optical properties. The optical properties of each simulated phytoplankton cell are integrated over an entire size distribution and are provided at high spectral resolution. The model code is additionally included to enable user access to the complete set of wavelength-dependent, angularly resolved volume scattering functions. This optically coherent dataset of hyperspectral optical properties for a set of globally significant phytoplankton groups has potential for use in algorithm development towards the optimal exploitation of the new age of hyperspectral satellite radiometry.

There is a pressing need for improved bio-optical models of high biomass waters as eutrophication of coastal and inland waters becomes an increasing problem. Seasonal boom conditions in the Southern Benguela and persistent harmful algal production in various inland waters in Southern Africa present valuable opportunities for the development of such modelling capabilities. The phytoplankton-dominated signal of these waters additionally addresses an increased interest in Phytoplankton Functional Type (PFT) analysis. To these ends, an initial validation of a new model of Equivalent Algal Populations (EAP) is presented here. This paper makes a first order comparison of two prominent phytoplankton Inherent Optical Property (IOP) models with the EAP model, which places emphasis on explicit bio-physical modelling of the phytoplankton population as a holistic determinant of inherent optical properties. This emphasis is shown to have an impact on the ability to retrieve the detailed phytoplankton spectral scattering information necessary for PFT applications and to successfully simulate reflectance across wide ranges of physical environments, biomass, and assemblage characteristics.

Synoptic scale knowledge of the size structure of phytoplankton communities can offer insight in to primary ecosystem diversity and biogeochemical variability from operational to the decadal scales. Accordingly, obtaining estimates of size and other phytoplankton functional type descriptors within known confidence limits from remotely sensed data has become a major objective to extend the use of ocean colour data beyond chlorophyll a retrievals. Here, a new forward and inverse modelling structure is proposed to determine information about the cell size of phytoplankton communities using Standard size distributions of two layered spheres to derive a full suite of algal inherent optical properties for a coupled radiative transfer model. This new capability allows explicit quantification of the remote sensing reflectance signal attributable to changes in phytoplankton cell size. Inversion of this model reveals regions within the parameter space where ambiguity may limit potential of inversion algorithms. Validation of the algorithm within the Benguela upwelling system using independent data shows promise for ecosystem applications and further investigation of the interaction between phytoplankton functional types and optical signals. The results here suggest that the utility of assemblage related signals in spectral reflectance is highly sensitive to algal biomass, the presence of other absorbing and scattering constituents and the resultant constituent-specific inherent optical property budget. As such, optimal methods for determining phytoplankton size from (in situ or satellite) ocean colour data will likely rely on appropriately spectrally dense and optimised sensors, well characterised measurement errors including those from atmospheric correction, and an ability to appropriately limit ambiguity within the context of regional inherent optical properties.

There is increasing interdisciplinary interest in phytoplankton community dynamics as the growing environmental problems of water quality (particularly eutrophication) and climate change demand attention. This has led to a pressing need for improved biophysical and causal understanding of Phytoplankton Functional Type (PFT) optical signals, in order for satellite radiometry to be used to detect ecologically relevant phytoplankton assemblage changes. Biophysically and biogeochemically consistent phytoplankton Inherent Optical Property (IOP) models play an important role in achieving this understanding, as the optical effects of phytoplankton assemblage changes can be examined systematically in relation to the bulk optical water-leaving signal. The Equivalent Algal Populations (EAP) model is used here to investigate the source and magnitude of size- and pigment- driven PFT signals in the water-leaving reflectance, as well as the potential to detect these using satellite radiometry. This model places emphasis on the determination of biophysically consistent phytoplankton IOPs, with both absorption and scattering determined by mathematically cogent relationships to the particle complex refractive indices. All IOPs are integrated over an entire size distribution. A distinctive attribute is the model’s comprehensive handling of the spectral and angular character of phytoplankton scattering. Selected case studies and sensitivity analyses reveal that phytoplankton spectral scattering is most useful and the least ambiguous driver of the PFT signal. Key findings are that there is the most sensitivity in phytoplankton backscatter ( $b_{b\varphi }$ ) in the 1–6 $\mathsf{\mu }$ m size range; the backscattering-driven signal in the 520 to 570 nm region is the critical PFT identifier at marginal biomass, and that, while PFT information does appear at blue wavelengths, absorption-driven signals are compromised by ambiguity due to biomass and non-algal absorption. Low signal in the red, due primarily to absorption by water, inhibits PFT detection here. The study highlights the need to quantitatively understand the constraints imposed by phytoplankton biomass and the IOP budget on the assemblage-related signal. A proportional phytoplankton contribution of approximately 40% to the total $b_b$ appears to a reasonable minimum threshold in terms of yielding a detectable optical change in $R_{rs}$ . We hope these findings will provide considerable insight into the next generation of PFT algorithms.

The accurate description of a water body’s volume scattering function (VSF), and hence its phase functions, is critical to the determination of the constituent inherent optical properties (IOPs), the associated spectral water-leaving reflectance, and consequently the retrieval of phytoplankton functional type (PFT) information. The equivalent algal populations (EAP) model has previously been evaluated for phytoplankton-dominated waters, and offers the ability to provide phytoplankton population-specific phase functions, unveiling a new opportunity to further understanding of the causality of the PFT signal. This study presents and evaluates the wavelength dependent, spectrally variable EAP particle phase functions and the subsequent effects on water-leaving reflectance. Comparisons are made with frequently used phase function approximations e.g. the Fournier Forand formulation, as well as with phase functions inferred from measured VSFs in coastal waters. Relative differences in shape and magnitude are quantified. Reflectance modelled with the EAP phase functions is then compared against measured reflectance data from phytoplankton-dominated waters. Further examples of modelled phytoplankton-dominated waters are discussed with reference to choice of phase function for two PFTs (eukaryote and prokaryote) across a range of biomass. Finally a demonstration of the sensitivity of reflectance due to the choice of phase function is presented. The EAP model phase functions account for both spectral and angular variability in phytoplankton backscattering i.e. they display variability which is both spectral and shape-related. It is concluded that phase functions modelled in this way are necessary for investigating the effects of assemblage variability on the ocean colour signal, and should be considered for model closure even in relatively low scattering conditions where phytoplankton dominate the IOPs.