Global Retrieval of Diatom Abundance Based on Phytoplankton Pigments and Satellite Data

Phytoplankton thermal trait parameterization alters community structure and biogeochemical processes in a modeled ocean

Phytoplankton exhibit diverse physiological responses to temperature which influence their fitness in the environment and consequently alter their community structure. Here, we explored the sensitivity of phytoplankton community structure to thermal response parameterization in a modelled marine phytoplankton community. Using published empirical data, we evaluated the maximum thermal growth rates (μmax ) and temperature coefficients (Q10 ; the rate at which growth scales with temperature) of six key Phytoplankton Functional Types (PFTs): coccolithophores, cyanobacteria, diatoms, diazotrophs, dinoflagellates, and green algae. Following three well-documented methods, PFTs were either assumed to have (1) the same μmax and the same Q10 (as in to Eppley, 1972), (2) a unique μmax but the same Q10 (similar to Kremer et al., 2017), or (3) a unique μmax and a unique Q10 (following Anderson et al., 2021). These trait values were then implemented within the Massachusetts Institute of Technology biogeochemistry and ecosystem model (called Darwin) for each PFT under a control and climate change scenario. Our results suggest that applying a μmax and Q10 universally across PFTs (as in Eppley, 1972) leads to unrealistic phytoplankton communities, which lack diatoms globally. Additionally, we find that accounting for differences in the Q10 between PFTs can significantly impact each PFT's competitive ability, especially at high latitudes, leading to altered modeled phytoplankton community structures in our control and climate change simulations. This then impacts estimates of biogeochemical processes, with, for example, estimates of export production varying by ~10% in the Southern Ocean depending on the parameterization. Our results indicate that the diversity of thermal response traits in phytoplankton not only shape community composition in the historical and future, warmer ocean, but that these traits have significant feedbacks on global biogeochemical cycles.

Plankton Imagery Data Inform Satellite-Based Estimates of Diatom Carbon

Estimating the biomass of phytoplankton communities via remote sensing is a key requirement for understanding global ocean ecosystems. Of particular interest is the carbon associated with diatoms given their unequivocal ecological and biogeochemical roles. Satellite-based algorithms often rely on accessory pigment proxies to define diatom biomass, despite a lack of validation against independent diatom biomass measurements. We used imaging-in-flow cytometry to quantify diatom carbon in the western North Atlantic, and compared results to those obtained from accessory pigment-based approximations. Based on this analysis, we offer a new empirical formula to estimate diatom carbon concentrations from chlorophyll a. Additionally, we developed a neural network model in which we integrated chlorophyll a and environmental information to estimate diatom carbon distributions in the western North Atlantic. The potential for improving satellite-based diatom carbon estimates by integrating environmental information into a model, compared to models that are based solely on chlorophyll a, is discussed.

Decadal changes in global phytoplankton compositions influenced by biogeochemical variables

The global environmental changes owing to natural and anthropogenic influences are challenging the structure and functioning of the ocean ecosystem. The complex processes interacting within the physical, chemical, and biological environment at different spatio-temporal scales and their impact on the ocean ecosystem processes are yet to be investigated. A long term trend on phytoplankton biomass in terms of Chlorophyll-a concentration (Chl-a), phytoplankton compositions and the processes that control the variability is required for understanding the ocean ecosystem. This study investigated decadal trends (2002-2015) of phytoplankton composition and biogeochemical parameters over the Global Ocean (GO), Southern Ocean (SO), and the Arctic Ocean (AO) using ocean color remote sensing and assimilated data from the National Aeronautics and Space Administration (NASA) Ocean biogeochemical model. The results revealed the dominance of larger cell phytoplankton mainly diatoms throughout the SO and AO; however, the coccolithophores dominate in the remaining part of the GO. Analysis of nutrients showed that nitrate is not a limiting factor for the variability of phytoplankton biomass in the SO and AO. The low nitrate concentration influenced in the rest of the GO. The photosynthetically available radiation (PAR) limiting the phytoplankton biomass and composition in the SO and AO. Although the SO is known as the high nutrient low chlorophyll (HNLC) region of the GO, the low iron concentration along with the PAR co-limits the growth of phytoplankton biomass. Trend analysis showed that an increase in Chl-a and diatoms in the SO and AO. In contrast, it declined significantly in the other regions of the GO, in response to the consistent increase in sea surface temperature. The results indicated that, shifting of phytoplankton community from regional to global scale have a greater implication for climate change and marine ecosystem.

Seasonal bias in global ocean color observations

In this study, we identify a seasonal bias in the ocean color satellite-derived remote sensing reflectances (Rrs(λ);sr-1) at the ocean color validation site, Marine Optical BuoY. The seasonal bias in Rrs(λ) is present to varying degrees in all ocean color satellites examined, including the Visible Infrared Imaging Radiometer Suite, Sea-Viewing Wide Field-of-View Sensor, and Moderate Resolution Imaging Spectrometer. The relative bias in Rrs has spectral dependence. Products derived from Rrs(λ) are affected by the bias to varying degrees, with particulate backscattering varying up to 50% over a year, chlorophyll varying up to 25% over a year, and absorption from phytoplankton or dissolved material varying by up to 15%. The propagation of Rrs(λ) bias into derived products is broadly confirmed on regional and global scales using Argo floats and data from the cloud-aerosol lidar with orthogonal polarization instrument aboard the cloud-aerosol lidar and infrared pathfinder satellite. The artifactual seasonality in ocean color is prominent in areas of low biomass (i.e., subtropical gyres) and is not easily discerned in areas of high biomass. While we have eliminated several candidates that could cause the biases in Rrs(λ), there are still outstanding questions regarding potential contributions from atmospheric corrections. Specifically, we provide evidence that the aquatic bidirectional reflectance distribution function may in part cause the observed seasonal bias, but this does not preclude an additional effect of the aerosol estimation. Our investigation highlights the contributions that atmospheric correction schemes can make in introducing biases in Rrs(λ), and we recommend more simulations to discern these influence Rrs(λ) biases. Community efforts are needed to find the root cause of the seasonal bias because all past, present, and future data are, or will be, affected until a solution is implemented.

Evaluation of diagnostic pigments to estimate phytoplankton size classes

Phytoplankton accessory pigments are commonly used to estimate phytoplankton size classes, particularly during development and validation of biogeochemical models and satellite ocean color-based algorithms. The diagnostic pigment analysis (DPA) is based on bulk measurements of pigment concentrations and relies on assumptions regarding the presence of specific pigments in different phytoplankton taxonomic groups. Three size classes are defined by the DPA: picoplankton, nanoplankton, and microplankton. Until now, the DPA has not been evaluated against an independent approach that provides phytoplankton size calculated on a per-cell basis. Automated quantitative cell imagery of microplankton and some nanoplankton, used in combination with conventional flow cytometry for enumeration of picoplankton and nanoplankton, provide a novel opportunity to perform an independent evaluation of the DPA. Here, we use a data set from the North Atlantic Ocean that encompasses all seasons and a wide range of chlorophyll concentrations (0.18-5.14 mg m^-3). Results show that the DPA overestimates microplankton and picoplankton when compared to cytometry data, and subsequently underestimates the contribution of nanoplankton to total biomass. In contrast to the assumption made by the DPA that the microplankton size class is largely made up of diatoms and dinoflagellates, imaging-in-flow cytometry shows significant presence of diatoms and dinoflagellates in the nanoplankton size class. Additionally, chlorophyll b is commonly attributed solely to picoplankton by the DPA, but Chl b-containing phytoplankton are observed with imaging in both nanoplankton and microplankton size classes. We suggest revisions to the DPA equations and application of uncertainties when calculating size classes from diagnostic pigments.