Retrieving Pigment Concentrations Based on Hyperspectral Measurements of the Phytoplankton Absorption Coefficient in Global Oceans

Phytoplankton communities, which can be easily observed by optical sensors deployed on various types of platforms over diverse temporal and spatial scales, are crucial to marine ecosystems and biogeochemical cycles, and accurate pigment concentrations make it possible to effectively derive information from them. To date, there is no practical approach, however, to retrieving concentrations of detailed pigments from phytoplankton absorption coefficients (aph) with acceptable accuracy and robustness in global oceans. In this study, a novel method, which is a stepwise regression method improved by early stopping (the ES-SR method) based on the derivative of hyperspectral aph, was proposed to retrieve pigment concentrations. This method was developed from an extensive global dataset collected from layers at different depths and contains phytoplankton pigment concentrations and aph. In the case of the logarithm, strong correlations were found between phytoplankton pigment concentrations and the absolute values of the second derivative (aph″)/the fourth derivative (aph4) of aph. According to these correlations, the ES-SR method is effective in obtaining the characteristic wavelengths of phytoplankton pigments for pigment concentration inversion. Compared with the Gaussian decomposition method and principal component regression method, which are based on the derivatives, the ES-SR method implemented on aph″ is the optimum approach with the greatest accuracy for each phytoplankton pigment. More than half of the determination coefficient values (R2log) for all pigments, which were retrieved by performing the ES-SR method on aph″, exceeded 0.7. The values retrieved for all pigments fit well to the one-to-one line with acceptable root mean square error (RMSElog: 0.146–0.508) and median absolute percentage error (MPElog: 8.2–28.5%) values. Furthermore, the poor correlations between the deviations from the values retrieved by the ES-SR method and impact factors related to pigment composition and cell size class show that this method has advantageous robustness. Therefore, the ES-SR method has the potential to effectively monitor phytoplankton community information from hyperspectral optical data in global oceans.

Evaluating semi-analytical algorithms for estimating inherent optical properties in the South China Sea

Using large amounts of bio-optical data collected in the South China Sea (SCS) from 2003 to 2016, this study checks the consistency between well-known semi-analytical algorithms (SAAs)-the quasi-analytical algorithm (QAA) and the default generalized inherent optical property (GIOP-DC)-in retrieving the non-water absorption coefficient (a_nw(λ)), phytoplankton absorption coefficient (a_ph(λ)) and particulate backscattering coefficient (b_bp(λ)) from remote-sensing reflectance (R_rs(λ)) data at 412, 443, 490, 531, and 555 nm. The samples from the SCS are further separated into oligotrophic and mesotrophic water types for the comparison of the SAAs. Several findings are made: First, the values of an_w(λ) derived from the two SAAs deliver similar performance, with R² values ranging from 0.74 to 0.85 and 0.74 to 0.87, implying absolute percent error differences (APDs) from 37.93% to 74.88% and from 32.32% to 71.75% for the QAA and GIOP-DC, respectively. The QAA shows a value of R² between 0.64 and 0.91 and APDs between 43.57% to 83.53%, while the GIOP-DC yields R² between 0.76 to 0.89 and APDs between 44.65% to 79.46% when estimating a_ph(λ). The values of b_bp(λ) derived from the QAA are closer to the in-situ b_bp(λ) values, as indicated by the low values of the normalized centered root-mean-square deviation and normalized standard deviation, which are close to one. Second, a regionally tuned estimation of a_ph(λ) is proposed and recommended for the SCS. This consistency check of inherent optical properties products from SAAs can serve as reference for algorithm selection for further applications, including primary production, carbon, and biogeochemical models of the SCS, and can provide guidance for improving a_ph(λ) estimation.

Evaluation of Satellite-Based Algorithms to Retrieve Chlorophyll-a Concentration in the Canadian Atlantic and Pacific Oceans

Remote-sensing reflectance data collected by ocean colour satellites are processed using bio-optical algorithms to retrieve biogeochemical properties of the ocean. One such important property is the concentration of chlorophyll-a, an indicator of phytoplankton biomass that serves a multitude of purposes in various ocean science studies. Here, the performance of two generic chlorophyll-a algorithms (i.e., a band ratio one, Ocean Colour X (OCx), and a semi-analytical one, Garver–Siegel Maritorena (GSM)) was assessed against two large in situ datasets of chlorophyll-a concentration collected between 1999 and 2016 in the Northeast Pacific (NEP) and Northwest Atlantic (NWA) for three ocean colour sensors: Sea-viewing Wide Field-of-view Sensor (SeaWiFS), Moderate Resolution Imaging Spectroradiometer (MODIS), and Visible Infrared Imaging Radiometer Suite (VIIRS). In addition, new regionally-tuned versions of these two algorithms are presented, which reduced the mean error (mg m−3) of chlorophyll-a concentration modelled by OCx in the NWA from −0.40, −0.58 and −0.45 to 0.037, −0.087 and −0.018 for MODIS, SeaWiFS, and VIIRS respectively, and −0.34 and −0.36 to −0.0055 and −0.17 for SeaWiFS and VIIRS in the NEP. An analysis of the uncertainties in chlorophyll-a concentration retrieval showed a strong seasonal pattern in the NWA, which could be attributed to changes in phytoplankton community composition, but no long-term trends were found for all sensors and regions. It was also found that removing the 443 nm waveband for the OCx algorithms significantly improved the results in the NWA. Overall, GSM performed better than the OCx algorithms in both regions for all three sensors but generated fewer chlorophyll-a retrievals than the OCx algorithms.

The Influence of Temperature and Community Structure on Light Absorption by Phytoplankton in the North Atlantic

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.

Factors Regulating the Relationship Between Total and Size-Fractionated Chlorophyll-a in Coastal Waters of the Red Sea

Phytoplankton biomass and size structure are recognized as key ecological indicators. With the aim to quantify the relationship between these two ecological indicators in tropical waters and understand controlling factors, we analyzed the total chlorophyll-a concentration, a measure of phytoplankton biomass, and its partitioning into three size classes of phytoplankton, using a series of observations collected at coastal sites in the central Red Sea. Over a period of 4 years, measurements of flow cytometry, size-fractionated chlorophyll-a concentration, and physical-chemical variables were collected near Thuwal in Saudi Arabia. We fitted a three-component model to the size-fractionated chlorophyll-a data to quantify the relationship between total chlorophyll and that in three size classes of phytoplankton [pico- (<2 μm), nano- (2–20 μm) and micro-phytoplankton (>20 μm)]. The model has an advantage over other more empirical methods in that its parameters are interpretable, expressed as the maximum chlorophyll-a concentration of small phytoplankton (pico- and combined pico-nanophytoplankton, Cpm and Cp,nm, respectively) and the fractional contribution of these two size classes to total chlorophyll-a as it tends to zero (D_p and D_p,n). Residuals between the model and the data (model minus data) were compared with a range of other environmental variables available in the dataset. Residuals in pico- and combined pico-nanophytoplankton fractions of total chlorophyll-a were significantly correlated with water temperature (positively) and picoeukaryote cell number (negatively). We conducted a running fit of the model with increasing temperature and found a negative relationship between temperature and parameters Cpm and Cp,nm and a positive relationship between temperature and parameters D_p and D_p,n. By harnessing the relative red fluorescence of the flow cytometric data, we show that picoeukaryotes, which are higher in cell number in winter (cold) than summer (warm), contain higher chlorophyll per cell than other picophytoplankton and are slightly larger in size, possibly explaining the temperature shift in model parameters, though further evidence is needed to substantiate this finding. Our results emphasize the importance of knowing the water temperature and taxonomic composition of phytoplankton within each size class when understanding their relative contribution to total chlorophyll. Furthermore, our results have implications for the development of algorithms for inferring size-fractionated chlorophyll from satellite data, and for how the partitioning of total chlorophyll into the three size classes may change in a future ocean.

Remote sensing of spatial and temporal patterns of phytoplankton assemblages in the Bohai Sea, Yellow Sea, and east China sea

Marine phytoplankton accounts for roughly half the planetary primary production, and plays significant roles in marine ecosystem functioning, physical and biogeochemical processes, and climate changes. Documenting phytoplankton assemblages' dynamics, particularly their community structure properties, is thus a crucial and also challenging task. A large number of in situ and space-borne observation datasets are collected that cover the marginal seas in the west Pacific, including Bohai Sea, Yellow Sea, and East China Sea. Here, a customized region-specific semi-analytical model is developed in order to detect phytoplankton community structure properties (using phytoplankton size classes, PSCs, as its first-order delegate), and repeatedly tested to assure its reliable performance. Independent in situ validation datasets generate relatively low and acceptable predictive errors (e.g., mean absolute percentage errors, MAPE, are 38.4%, 22.7%, and 34.4% for micro-, nano-, and picophytoplankton estimations, respectively). Satellite synchronization verification also produces comparative predictive errors. By applying this model to long time-series of satellite data, we document the past two-decadal (namely from 1997 to 2017) variation on the PSCs. Satellite-derived records reveal a general spatial distribution rule, namely microphytoplankton accounts for most variation in nearshore regions, when nanophytoplankton dominates offshore water areas, together with a certain high contribution from picophytoplankton. Long time-series of data records indicate a roughly stable tendency during the period of the past twenty years, while there exist periodical changes in a short-term one-year scale. High covariation between marine environment factors and PSCs are further found, with results that underwater light field and sea surface temperature are the two dominant climate variables which exhibit a good ability to multivariate statistically model the PSCs changes in these marginal seas. Specifically, three types of influence induced by underwater light field and sea surface temperature can be generalized to cover different water conditions and regions, and meanwhile a swift response time (approximately < 1 month) of phytoplankton to the changing external environment conditions is found by the wavelet analysis. This study concludes that phytoplankton community structures in the marginal seas remain stable and are year-independent over the past two decades, together with a short-term in-year cycle; this change rule need to be considered in future oceanographic studies.

Retrieving Phytoplankton Size Class from the Absorption Coefficient and Chlorophyll A Concentration Based on Support Vector Machine

The phytoplankton size class (PSC) plays an important role in biogeochemical processes in the ocean. In this study, a regional model of PSCs is proposed to retrieve vertical PSCs from the total minus water absorption coefficient (at-w(λ)) and Chlorophyll a concentration (Chla). The PSC model is developed by first reconstructing phytoplankton absorption and Chla from at-w(λ), and then extracting PSC from them using the support vector machine (SVM). In situ bio-optical data collected in the South China Sea from 2006 to 2013 were used to train the SVM. The proposed PSC model was subsequently validated using an independent PSC dataset from the Northeast South China Sea Cruise in 2015. The results indicate that the PSC model performed better than the three components model, with a value of r2 between 0.35 and 0.66, and the absolute percentage difference between 56% and 181%. On the whole, our PSC model shows a remarkable utility in terms of inferring vertical PSCs from the South China Sea.

The Fundamental Contribution of Phytoplankton Spectral Scattering to Ocean Colour: Implications for Satellite Detection of Phytoplankton Community Structure

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 ϕ ) in the 1–6 μ 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 r s . We hope these findings will provide considerable insight into the next generation of PFT algorithms.

Concentrations of Multiple Phytoplankton Pigments in the Global Oceans Obtained from Satellite Ocean Color Measurements with MERIS

The remote sensing of chlorophyll a concentration from ocean color satellites has been an essential variable quantifying phytoplankton in the past decades, yet estimation of accessory pigments from ocean color remote sensing data has remained largely elusive. In this study, we validated the concentrations of multiple pigments (Cpigs) retrieved from in situ and MEdium Resolution Imaging Spectrometer (MERIS) measured remote sensing reflectance (Rrs(λ)) in the global oceans. A multi-pigment inversion model (MuPI) was used to semi-analytically retrieve Cpigs from Rrs(λ). With a set of globally optimized parameters, the accuracy of the retrievals obtained with MuPI is quite promising. Compared with High-Performance Liquid Chromatography (HPLC) measurements near Bermuda, the concentrations of chlorophyll a, b, c ([Chl-a], [Chl-b], [Chl-c]), photoprotective carotenoids ([PPC]), and photosynthetic carotenoids ([PSC]) can be retrieved from MERIS data with a mean unbiased absolute percentage difference of 38%, 78%, 65%, 36%, and 47%, respectively. The advantage of the MuPI approach is the simultaneous retrievals of [Chl-a] and the accessory pigments [Chl-b], [Chl-c], [PPC], [PSC] from MERIS Rrs(λ) based on a closure between the input and output Rrs(λ) spectra. These results can greatly expand scientific studies of ocean biology and biogeochemistry of the global oceans that are not possible when the only available information is [Chl-a].

An overview of approaches and challenges for retrieving marine inherent optical properties from ocean color remote sensing

Ocean color measured from satellites provides daily global, synoptic views of spectral waterleaving reflectances that can be used to generate estimates of marine inherent optical properties (IOPs). These reflectances, namely the ratio of spectral upwelled radiances to spectral downwelled irradiances, describe the light exiting a water mass that defines its color. IOPs are the spectral absorption and scattering characteristics of ocean water and its dissolved and particulate constituents. Because of their dependence on the concentration and composition of marine constituents, IOPs can be used to describe the contents of the upper ocean mixed layer. This information is critical to further our scientific understanding of biogeochemical oceanic processes, such as organic carbon production and export, phytoplankton dynamics, and responses to climatic disturbances. Given their importance, the international ocean color community has invested significant effort in improving the quality of satellite-derived IOP products, both regionally and globally. Recognizing the current influx of data products into the community and the need to improve current algorithms in anticipation of new satellite instruments (e.g., the global, hyperspectral spectroradiometer of the NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission), we present a synopsis of the current state of the art in the retrieval of these core optical properties. Contemporary approaches for obtaining IOPs from satellite ocean color are reviewed and, for clarity, separated based their inversion methodology or the type of IOPs sought. Summaries of known uncertainties associated with each approach are provided, as well as common performance metrics used to evaluate them. We discuss current knowledge gaps and make recommendations for future investment for upcoming missions whose instrument characteristics diverge sufficiently from heritage and existing sensors to warrant reassessing current approaches.

Retrieval of phytoplankton cell size from chlorophyll a specific absorption and scattering spectra of phytoplankton

Phytoplankton cell size is an important property that affects diverse ecological and biogeochemical processes, and analysis of the absorption and scattering spectra of phytoplankton can provide important information about phytoplankton size. In this study, an inversion method for extracting quantitative phytoplankton cell size data from these spectra was developed. This inversion method requires two inputs: chlorophyll a specific absorption and scattering spectra of phytoplankton. The average equivalent-volume spherical diameter (ESD_v) was calculated as the single size approximation for the log-normal particle size distribution (PSD) of the algal suspension. The performance of this method for retrieving cell size was assessed using the datasets from cultures of 12 phytoplankton species. The estimations of a(λ) and b(λ) for the phytoplankton population using ESD_v had mean error values of 5.8%-6.9% and 7.0%-10.6%, respectively, compared to the a(λ) and b(λ) for the phytoplankton populations using the log-normal PSD. The estimated values of C_iESD_v were in good agreement with the measurements, with r²=0.88 and relative root mean square error (NRMSE)=25.3%, and relatively good performances were also found for the retrieval of ESD_v with r²=0.78 and NRMSE=23.9%.

Remote estimation of phytoplankton size fractions using the spectral shape of light absorption

Phytoplankton size structure plays an important role in ocean biogeochemical processes. The light absorption spectra of phytoplankton provide a great potential for retrieving phytoplankton size structure because of the strong dependence on the packaging effect caused by phytoplankton cell size and on different pigment compositions related to phytoplankton taxonomy. In this study, we investigated the variability in light absorption spectra of phytoplankton in relation to the size structure. Based on this, a new approach was proposed for estimating phytoplankton size fractions. Our approach use the spectral shape of the normalized phytoplankton absorption coefficient (a(ph)(λ)) through principal component analysis (PCA). Values of a(ph)(λ) were normalized to remove biomass effects, and PCA was conducted to separate the spectral variance of normalized a(ph)(λ) into uncorrelated principal components (PCs). Spectral variations captured by the first four PC modes were used to build relationships with phytoplankton size fractions. The results showed that PCA had powerful ability to capture spectral variations in normalized a(ph)(λ), which were significantly related to phytoplankton size fractions. For both hyperspectral a(ph)(λ) and multiband a(ph)(λ), our approach is applicable. We evaluated our approach using wide in situ data collected from coastal waters and the global ocean, and the results demonstrated a good and robust performance in estimating phytoplankton size fractions in various regions. The model performance was further evaluated by a(ph)(λ) derived from in situ remote sensing reflectance (R(rs)(λ)) with a quasi-analytical algorithm. Using R(rs)(λ) only at six bands, accurate estimations of phytoplankton size fractions were obtained, with R(2) values of 0.85, 0.61, and 0.76, and root mean-square errors of 0.130, 0.126, and 0.112 for micro-, nano-, and picophytoplankton, respectively. Our approach provides practical basis for remote estimation of phytoplankton size structure using a(ph)(λ) derived from satellite observations or rapid field instrument measurements in the future.

Discrimination of phytoplankton functional groups using an ocean reflectance inversion model

Ocean reflectance inversion models (ORMs) provide a mechanism for inverting the color of the water observed by a satellite into marine inherent optical properties (IOPs), which can then be used to study phytoplankton community structure. Most ORMs effectively separate the total signal of the collective phytoplankton community from other water column constituents; however, few have been shown to effectively identify individual contributions by multiple phytoplankton groups over a large range of environmental conditions. We evaluated the ability of an ORM to discriminate between Noctiluca miliaris and diatoms under conditions typical of the northern Arabian Sea. We: (1) synthesized profiles of IOPs that represent bio-optical conditions for the Arabian Sea; (2) generated remote-sensing reflectances from these profiles using Hydrolight; and (3) applied the ORM to the synthesized reflectances to estimate the relative concentrations of diatoms and N. miliaris. By comparing the estimates from the inversion model with those from synthesized vertical profiles, we identified those conditions under which the ORM performs both well and poorly. Even under perfectly controlled conditions, the absolute accuracy of ORM retrievals degraded when further deconstructing the derived total phytoplankton signal into subcomponents. Although the absolute magnitudes maintained biases, the ORM successfully detected whether or not Noctiluca miliaris appeared in the simulated water column. This quantitatively calls for caution when interpreting the absolute magnitudes of the retrievals, but qualitatively suggests that the ORM provides a robust mechanism for identifying the presence or absence of species.

Novel method for quantifying the cell size of marine phytoplankton based on optical measurements

Phytoplankton size is important for the pelagic food web and oceanic ecosystems. However, the size of phytoplankton is difficult to quantify because of methodological constraints. To address this limitation, we have exploited the phytoplankton package effect to develop a new method for estimating the mean cell size of individual phytoplankton populations. This method was validated using a data set that contained simultaneous measurements of phytoplankton absorption and cell size distributions from 13 phytoplankton species. Comparing with existing methods, our method is more efficient with good accuracy, and it could potentially be applied in current in situ optical instruments.

Multivariate approach for the retrieval of phytoplankton size structure from measured light absorption spectra in the Mediterranean Sea (BOUSSOLE site)

Models based on the multivariate partial least squares (PLS) regression technique are developed for the retrieval of phytoplankton size structure from measured light absorption spectra (BOUSSOLE site, northwestern Mediterranean Sea). PLS-models trained with data from the Mediterranean Sea showed good accuracy in retrieving, over the nine-year BOUSSOLE time series, the concentrations of total chlorophyll a [Tchl a], of the sum of seven diagnostic pigments and of pigments associated with micro, nano, and picophytoplankton size classes separately. PLS-models trained using either total particle or phytoplankton absorption spectra performed similarly, and both reproduced seasonal variations of biomass and size classes derived by high performance liquid chromatography. Satisfactory retrievals were also obtained using PLS-models trained with a data set including various locations of the world's oceans, with however a lower accuracy. These results open the way to an application of this method to absorption spectra derived from hyperspectral and field satellite radiance measurements.

Generalized ocean color inversion model for retrieving marine inherent optical properties

Ocean color measured from satellites provides daily, global estimates of marine inherent optical properties (IOPs). Semi-analytical algorithms (SAAs) provide one mechanism for inverting the color of the water observed by the satellite into IOPs. While numerous SAAs exist, most are similarly constructed and few are appropriately parameterized for all water masses for all seasons. To initiate community-wide discussion of these limitations, NASA organized two workshops that deconstructed SAAs to identify similarities and uniqueness and to progress toward consensus on a unified SAA. This effort resulted in the development of the generalized IOP (GIOP) model software that allows for the construction of different SAAs at runtime by selection from an assortment of model parameterizations. As such, GIOP permits isolation and evaluation of specific modeling assumptions, construction of SAAs, development of regionally tuned SAAs, and execution of ensemble inversion modeling. Working groups associated with the workshops proposed a preliminary default configuration for GIOP (GIOP-DC), with alternative model parameterizations and features defined for subsequent evaluation. In this paper, we: (1) describe the theoretical basis of GIOP; (2) present GIOP-DC and verify its comparable performance to other popular SAAs using both in situ and synthetic data sets; and, (3) quantify the sensitivities of their output to their parameterization. We use the latter to develop a hierarchical sensitivity of SAAs to various model parameterizations, to identify components of SAAs that merit focus in future research, and to provide material for discussion on algorithm uncertainties and future emsemble applications.

Impact of computational methods and spectral models on the retrieval of optical properties via spectral optimization

Spectral optimization algorithm (SOA) is a well-accepted scheme for the retrieval of water constituents from the measurement of ocean color radiometry. It defines an error function between the input and output remote sensing reflectance spectrum, with the latter modeled with a few variables that represent the optically active properties, while the variables are solved numerically by minimizing the error function. In this paper, with data from numerical simulations and field measurements as input, we evaluate four computational methods for minimization (optimization) for their efficiency and accuracy on solutions, and illustrate impact of bio-optical models on the retrievals. The four optimization routines are the Levenberg-Marquardt (LM), the Generalized Reduced Gradient (GRG), the Downhill Simplex Method (Amoeba), and the Simulated Annealing-Downhill Simplex (i.e. SA + Amoeba, hereafter abbreviated as SAA). The Garver-Siegel-Maritorena SOA model is used as a base to test these computational methods. It is observed that 1) LM is the fastest method, but SAA has the largest number of valid retrievals; 2) the quality of final solutions are strongly influenced by the forms of spectral models (or eigen functions); and 3) dynamically-varying eigen functions are necessary to obtain smaller errors for both reflectance spectrum and retrievals. Results of this study provide helpful guidance for the selection of a computational method and spectral models if an SOA scheme is to be used to process ocean color images.

Algal biogeography: metagenomics shows distribution of a picoplanktonic pelagophyte

How can we determine the distribution of uncultured marine microorganisms? Targeted metagenomics has provided the complete chloroplast genome sequence, and the distribution, for a picoplanktonic pelagophyte alga.

Particle backscattering as a function of chlorophyll and phytoplankton size structure in the open-ocean

Using an extensive database of in situ observations we present a model that estimates the particle backscattering coefficient as a function of the total chlorophyll concentration in the open-ocean (Case-1 waters). The parameters of the model include a constant background component and the chlorophyll-specific backscattering coefficients associated with small (<20 μm) and large (>20 μm) phytoplankton. The new model performed with similar accuracy when compared with a traditional power-law function, with the additional benefit of providing information on the role of phytoplankton size. The observed spectral-dependency (γ) of model parameters was consistent with past observations, such that γ associated with the small phytoplankton population was higher than that of large phytoplankton. Furthermore, γ associated with the constant background component suggests this component is likely attributed to submicron particles. We envisage that the model would be useful for improving Case-1 ocean-colour models, assimilating light into multi-phytoplankton ecosystem models and improving estimates of phytoplankton size structure from remote sensing.

Model of phytoplankton absorption based on three size classes

Using the phytoplankton size-class model of Brewin et al. [Ecol. Model.221, 1472 (2010)], the two-population absorption model of Sathyendranath et al. [Int. J. Remote. Sens.22, 249 (2001)] and Devred et al. [J. Geophys. Res.111, C03011 (2006)] is extended to three populations of phytoplankton, namely, picophytoplankton, nanophytoplankton, and microphytoplankton. The new model infers total and size-dependent phytoplankton absorption as a function of the total chlorophyll-a concentration. A main characteristic of the model is that all the parameters that describe it have biological or optical interpretation. The three-population model performs better than the two-population model at retrieving total phytoplankton absorption. Accounting for the contributions of picophytoplankton and nanophytoplankton, rather than the combination of both as in the two-population model, improved significantly the retrieval of phytoplankton absorption at low chlorophyll-a concentrations. Class-dependent specific absorption of phytoplankton derived using the model compares well with previously published models. However, the model presented in this paper provides the specific absorption of three size classes and is applicable to a continuum of chlorophyll-a concentrations. Absorption obtained from remotely sensed chlorophyll-a using our model compares well with in situ absorption measurements.

Retrieval of phytoplankton size from bio-optical measurements: theory and applications

The absorption coefficient of a substance distributed as discrete particles in suspension is less than that of the same material dissolved uniformly in a medium-a phenomenon commonly referred to as the flattening effect. The decrease in the absorption coefficient owing to flattening effect depends on the concentration of the absorbing pigment inside the particle, the specific absorption coefficient of the pigment within the particle, and on the diameter of the particle, if the particles are assumed to be spherical. For phytoplankton cells in the ocean, with diameters ranging from less than 1 µm to more than 100 µm, the flattening effect is variable, and sometimes pronounced, as has been well documented in the literature. Here, we demonstrate how the in vivo absorption coefficient of phytoplankton cells per unit concentration of its major pigment, chlorophyll a, can be used to determine the average cell size of the phytoplankton population. Sensitivity analyses are carried out to evaluate the errors in the estimated diameter owing to potential errors in the model assumptions. Cell sizes computed for field samples using the model are compared qualitatively with indirect estimates of size classes derived from high performance liquid chromatography data. Also, the results are compared quantitatively against measurements of cell size in laboratory cultures. The method developed is easy-to-apply as an operational tool for in situ observations, and has the potential for application to remote sensing of ocean colour data.