Expanding understanding of optical variability in Lake Superior with a 4-year dataset

Lake Superior is one of the largest freshwater lakes on our planet, but few optical observations have been made to allow for the development and validation of visible spectral satellite remote sensing products. The dataset described here focuses on coincidently observing inherent and apparent optical properties along with biogeochemical parameters. Specifically, we observe remote sensing reflectance, absorption, scattering, backscattering, attenuation, chlorophyll concentration, and suspended particulate matter over the ice-free months of 2013–2016. The dataset substantially increases the optical knowledge of the lake. In addition to visible spectral satellite algorithm development, the dataset is valuable for characterizing the variable light field, particle, phytoplankton, and colored dissolved organic matter distributions, and helpful in food web and carbon cycle investigations. The compiled data can be freely accessed at https://seabass.gsfc.nasa.gov/archive/URI/Mouw/LakeSuperior/

Deriving inherent optical properties from decomposition of hyperspectral non-water absorption

Semi-analytical algorithms (SAAs) developed for multispectral ocean color sensors have benefited from a variety of approaches for retrieving the magnitude and spectral shape of inherent optical properties (IOPs). SAAs generally follow two approaches: 1) simultaneous retrieval of all IOPs, resulting in pre-defined bio-optical models and spectral dependence between IOPs and 2) retrieval of bulk IOPs (absorption and backscattering) first followed by decomposition into separate components, allowing for independent retrievals of some components. Current algorithms used to decompose hyperspectral remotely-sensed reflectance into IOPs follow the first strategy. Here, a spectral deconvolution algorithm for incorporation into the second strategy is presented that decomposes at-w(λ) from in situ measurements and estimates absorption due to phytoplankton (aph(λ)) and colored detrital material (adg(λ)) free of explicit assumptions. The algorithm described here, Derivative Analysis and Iterative Spectral Evaluation of Absorption (DAISEA), provides estimates of aph(λ) and adg(λ) over a spectral range from 350 to 700 nm. Estimated aph(λ) and adg(λ) showed an average normalized root mean square difference of \u3c30% and \u3c20%, respectively, from 350 to 650 nm for the majority of optically distinct environments considered. Estimated Sdg median difference was \u3c20% for all environments considered, while distribution of Sdg uncertainty suggests that biogeochemical variability represented by Sdg can be estimated free of bias. DAISEA results suggest that hyperspectral satellite ocean color data will improve our ability to track biogeochemical processes affiliated with variability in adg(λ) and Sdg free of explicit assumptions

A Satellite Assessment of Environmental Controls of Phytoplankton Community Size Structure

Phytoplankton play a key role as the base of the marine food web and a crucial component in the Earth\u27s carbon cycle. There have been a few regional studies that have utilized satellite‐estimated phytoplankton functional type products in conjunction with other environmental metrics. Here we expand to a global perspective and ask, what are the physical drivers of phytoplankton composition variability? Using a variety of satellite‐observed ocean color products and physical properties spanning 1997–2015, we characterize spatial and temporal variability in phytoplankton community size structure in relation to satellite‐based physical drivers. We consider the relationships globally and by major thermal regimes (cold and warm), dominant size distribution, and chlorophyll concentration variability. Globally, euphotic depth is the most important parameter driving phytoplankton size variability and also over the majority of the high‐latitude ocean and the central gyres. In all other regions, size variability is driven by a balance of light and mode of nutrient delivery. We investigated the relationship between size composition and chlorophyll concentration and the physical drivers through correlation analysis. Changes in size composition over time are regionally varying and explained by temporal shifts in the varying physical conditions. These changes in phytoplankton size composition and the varying underlying physical drivers will ultimately impact carbon export and food web processes in our changing ocean

Global evaluation of particulate organic carbon flux parameterizations and implications for atmospheric pCO\u3csub\u3e2\u3c/sub\u3e

The shunt of photosynthetically derived particulate organic carbon (POC) from the euphotic zone and deep remineralization comprises the basic mechanism of the “biological carbon pump.” POC raining through the “twilight zone” (euphotic depth to 1 km) and “midnight zone” (1 km to 4 km) is remineralized back to inorganic form through respiration. Accurately modeling POC flux is critical for understanding the “biological pump” and its impacts on air‐sea CO2 exchange and, ultimately, long‐term ocean carbon sequestration. Yet commonly used parameterizations have not been tested quantitatively against global data sets using identical modeling frameworks. Here we use a single one‐dimensional physical‐biogeochemical modeling framework to assess three common POC flux parameterizations in capturing POC flux observations from moored sediment traps and thorium‐234 depletion. The exponential decay, Martin curve, and ballast model are compared to data from 11 biogeochemical provinces distributed across the globe. In each province, the model captures satellite‐based estimates of surface primary production within uncertainties. Goodness of fit is measured by how well the simulation captures the observations, quantified by bias and the root‐mean‐square error and displayed using “target diagrams.” Comparisons are presented separately for the twilight zone and midnight zone. We find that the ballast hypothesis shows no improvement over a globally or regionally parameterized Martin curve. For all provinces taken together, Martin\u27s b that best fits the data is [0.70, 0.98]; this finding reduces by at least a factor of 3 previous estimates of potential impacts on atmospheric pCO2 of uncertainty in POC export to a more modest range [−16 ppm, +12 ppm]

Expanding understanding of optical variability in Lake Superior with a 4-year dataset

Lake Superior is one of the largest freshwater lakes on our planet, but few optical observations have been made to allow for the development and validation of visible spectral satellite remote sensing products. The dataset described here focuses on coincidently observing inherent and apparent optical properties along with biogeochemical parameters. Specifically, we observe remote sensing reflectance, absorption, scattering, backscattering, attenuation, chlorophyll concentration, and suspended particulate matter over the ice-free months of 2013-2016. The dataset substantially increases the optical knowledge of the lake. In addition to visible spectral satellite algorithm development, the dataset is valuable for characterizing the variable light field, particle, phytoplankton, and colored dissolved organic matter distributions, and helpful in food web and carbon cycle investigations. The compiled data can be freely accessed at https://seabass.gsfc.nasa.gov/archive/URI/Mouw/LakeSuperior/

Hit or miss? Impact of time series resolution on resolving phytoplankton dynamics at hourly, weekly, and satellite remote sensing frequencies

Characterizing marine phytoplankton community variability is crucial to designing sampling strategies and interpreting time series. Satellite remote sensing, microscopy sampling, and flow through imaging systems have widely different resolutions: from weekly or monthly with microscopy sampling to daily when no cloud cover or glint is present with polar-orbiting satellites, and hourly for autonomous imaging instruments. To improve our understanding of data robustness against sampling resolution at different taxonomic levels, we analyze 2 yr of data from an Imaging FlowCytobot with hourly resolution and resample it to daily, satellite-temporal, and weekly microscopy sampling resolution. We show that weekly and satellite-temporal resolutions are sufficient to resolve general community composition but that the randomness of satellite-temporal resolution can result in overrepresenting or underrepresenting certain categories. While the yearly phytoplankton biomass bloom is detected in late winter by all four resolutions, category-specific yearly blooms are generally consistent in timing but often underestimated or missed by the weekly and satellite-temporal resolutions, introducing a bias in year-to-year comparisons. A minimum of biweekly sampling, particularly during known bloom periods, would lower the bias in such categories. Similarly, sampling time should be considered as daily variations are category-specific. Overall, morning and low tide sampling tended to have higher biomass. We provide tables for categories detected by the IFCB in Narragansett Bay with their major bloom characteristics and recorded daily variability to inform future sampling designs. These results provide tools to interpret past and future time series, including possible detection of specific taxonomic groups with targeted satellite algorithms

Optimization and assessment of phytoplankton size class algorithms for ocean color data on the Northeast U.S. continental shelf

The size structure of phytoplankton communities influences important ecological and biogeochemical processes, including the transfer of energy through marine food webs. A variety of algorithms have been developed to estimate phytoplankton size classes (PSCs) from satellite ocean color data. However, many of these algorithms were developed for application to the global ocean, and their performance in more productive, optically complex coastal and continental shelf regions warrants evaluation. In this study, several existing PSC models were applied in the Northeast U.S. continental shelf (NES) region and compared with in situ PSC estimates derived from a local HPLC pigment data set. The effect of regional re-parameterization and incorporation of sea surface temperature (SST) into existing abundance-based model frameworks was investigated and model performance was assessed using an independent data set. Abundance-based model re-parameterization alone did not result in significant improvement in model performance compared with other models. However, the inclusion of SST led to a consistent reduction in model error for all size classes. Of two absorption-based algorithms tested, the best performing approach displayed similar performance metrics to the regional SST-dependent abundance-based model. The SST-dependent model and the absorption-based method were applied to monthly composites of the NES region for April and September 2019 and qualitatively compared. The results highlight the benefit of considering SST in abundance-based models and the applicability of absorption-based PSC methods in optically complex regions

Bio-optical Properties of Cyanobacteria Blooms in Western Lake Erie

There is a growing use of remote sensing observations for detecting and quantifying freshwater cyanobacteria populations, yet the inherent optical properties of these communities in natural settings, fundamental to bio-optical algorithms, are not well known. Toward bridging this knowledge gap, we measured a full complement of optical properties in western Lake Erie during cyanobacteria blooms in the summers of 2013 and 2014. Our measurements focus attention on the optical uniqueness of cyanobacteria blooms, which have consequences for remote sensing and bio-optical modeling. We found the cyanobacteria blooms in the western basin during our field work were dominated by Microcystis, while the waters in the adjacent central basin were dominated by Planktothrix. Chlorophyll concentrations ranged from 1 to over 135 μg/L across the study area with the highest concentrations associated with Microcystis in the western basin. We observed large, amorphous colonial Microcystis structures in the bloom area characterized by high phytoplankton absorption and high scattering coefficients with a mean particle backscatter ratio at 443 nm \u3e 0.03, which is higher than other plankton types and more comparable to suspended inorganic sediments. While our samples contained mixtures of both, our analysis suggests high contributions to the measured scatter and backscatter coefficients from cyanobacteria. Our measurements provide new insights into the optical properties of cyanobacteria blooms, and indicate that current semi-analytic models are likely to have problems resolving a closed solution in these types of waters as many of our observations are beyond the range of existing model components. We believe that different algorithm or model approaches are needed for these conditions, specifically for phytoplankton absorption and particle backscatter components. From a remote sensing perspective, this presents a challenge not only in terms of a need for new algorithms, but also for determining when to apply the best algorithm for a given situation. These results are new in the sense that they represent a complete description of the optical properties of freshwater cyanobacteria blooms, and are likely to be representative of bloom conditions for other systems containing Microcystis cells and colonies

Characterizing CDOM Spectral Variability Across Diverse Regions and Spectral Ranges

Satellite remote sensing of colored dissolved organic matter (CDOM) has focused on CDOM absorption (aCDOM) at a reference wavelength, as its magnitude provides insight into the underwater light field and large-scale biogeochemical processes. CDOM spectral slope, SCDOM, has been treated as a constant or semiconstant parameter in satellite retrievals of aCDOM despite significant regional and temporal variabilities. SCDOM and other optical metrics provide insights into CDOM composition, processing, food web dynamics, and carbon cycling. To date, much of this work relies on fluorescence techniques or aCDOM in spectral ranges unavailable to current and planned satellite sensors (e.g., \u3c 300 nm). In preparation for anticipated future hyperspectral satellite missions, we take the first step here of exploring global variability in SCDOM and fit deviations in the aCDOM spectra using the recently proposed Gaussian decomposition method. From this, we investigate if global variability in retrieved SCDOM and Gaussian components is significant and regionally distinct. We iteratively decreased the spectral range considered and analyzed the number, location, and magnitude of fitted Gaussian components to understand if a reduced spectral range impacts information obtained within a common spectral window. We compared the fitted slope from the Gaussian decomposition method to absorption-based indices that indicate CDOM composition to determine the ability of satellite-derived slope to inform the analysis and modeling of large-scale biogeochemical processes. Finally, we present implications of the observed variability for remote sensing of CDOM characteristics via SCDOM

Microplankton fraction (Sfm) estimated form satellite output (OC-CCI), links to netCDF files

We provide satellite output that estimates phytoplankton size class as percent microplankton (Sfm, > 20 μm). The calculation of Sfm requires satellite spectral remote sensing reflectance (Rrs(λ)), chlorophyll concentration ([Chl]), and absorption due to dissolved and detrital matter (adg(λ)). These are taken from 4km, monthly, Ocean Color Climate Change Initiative (OC-CCI, version 3.0, www.esa-oceancolour-cci.org) products which are globally merged Sea-Viewing Wide Field-of-View Sensor (SeaWiFS), Medium Resolution Imaging Spectrometer (MERIS), Moderate Resolution Imaging Spectroradiometer (MODIS-Aqua), and Visible Infrared Imaging Radiometer Suite (VIIRS) imagery for a continuous record from 1997 through 2015. This is an absorption-based approach from Mouw and Yoder (2010) (doi:10.1029/2010JC006337) where the chlorophyll-specific absorption spectra for phytoplankton size class extremes, pico- (0.2-2 μm) and microplankton (> 20 μm), are weighted by Sfm. Sfm is estimated from a look-up table containing simulated [Chl], adg(443), Rrs(λ), and Sfm. For a given pixel, satellite-estimated [Chl] and adg(443), are used to narrow the search space within the look-up table. Of the remaining options, the closest simulated Rrs(λ) to the satellite-observed Rrs(λ) is selected and the associated Sfm is assigned