Spectral correlation in MODIS water-leaving reflectance retrieval uncertainty.

Spectral remote sensing reflectance, Rrs(λ) (sr-1), is the fundamental quantity used to derive a host of bio-optical and biogeochemical properties of the water column from satellite ocean color measurements. Estimation of uncertainty in those derived geophysical products is therefore dependent on knowledge of the uncertainty in satellite-retrieved Rrs. Furthermore, since the associated algorithms require Rrs at multiple spectral bands, the spectral (i.e., band-to-band) error covariance in Rrs is needed to accurately estimate the uncertainty in those derived properties. This study establishes a derivative-based approach for propagating instrument random noise, instrument systematic uncertainty, and forward model uncertainty into Rrs, as retrieved using NASA's multiple-scattering epsilon (MSEPS) atmospheric correction algorithm, to generate pixel-level error covariance in Rrs. The approach is applied to measurements from Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite and verified using Monte Carlo (MC) analysis. We also make use of this full spectral error covariance in Rrs to calculate uncertainty in phytoplankton pigment chlorophyll-a concentration (chla, mg/m3) and diffuse attenuation coefficient of downwelling irradiance at 490 nm (Kd(490), m-1). Accounting for the error covariance in Rrs generally reduces the estimated relative uncertainty in chla by ∼1-2% (absolute value) in waters with chla < 0.25 mg/m3 where the color index (CI) algorithm is used. The reduction is ∼5-10% in waters with chla > 0.35 mg/m3 where the blue-green ratio (OCX) algorithm is used. Such reduction can be higher than 30% in some regions. For Kd(490), the reduction by error covariance is generally ∼2%, but can be higher than 20% in some regions. The error covariance in Rrs is further verified through forward-calculating chla from MODIS-retrieved and in situ Rrs and comparing estimated uncertainty with observed differences. An 8-day global composite of propagated uncertainty shows that the goal of 35% uncertainty in chla can be achieved over deep ocean waters (chla ≤ 0.1 mg/m3). While the derivative-based approach generates reasonable error covariance in Rrs, some assumptions should be updated as our knowledge improves. These include the inter-band error correlation in top-of-atmosphere reflectance, and uncertainties in the calibration of MODIS 869 nm band, in ancillary data, and in the in situ data used for system vicarious calibration.

Bayesian approach to a generalized inherent optical property model.

Relationships between the absorption and backscattering coefficients of marine optical constituents and ocean color, or remote sensing reflectances Rrs(λ), can be used to predict the concentrations of these constituents in the upper water column. Standard inverse modeling techniques that minimize error between the modeled and observed Rrs(λ) break down when the number of products retrieved becomes similar to, or greater than, the number of different ocean color wavelengths measured. Furthermore, most conventional ocean reflectance inversion approaches, such as the default configuration of NASA's Generalized Inherent Optical Properties algorithm framework (GIOP-DC), require a priori definitions of absorption and backscattering spectral shapes. A Bayesian approach to GIOP is implemented here to address these limitations, where the retrieval algorithm minimizes both the error in retrieved ocean color and the deviation from prior knowledge, calculated using output from a mixture of empirically-derived and best-fit values. The Bayesian approach offers potential to produce an expanded range of parameters related to the spectral shape of absorption and backscattering spectra.

Optimizing retrieval spaces of bio-optical models for remote sensing of ocean color.

We investigated the optimal number of independent parameters required to accurately represent spectral remote sensing reflectances (R rs) by performing principal component analysis on quality controlled in situ and synthetic R rs data. We found that retrieval algorithms should be able to retrieve no more than four free parameters from R rs spectra for most ocean waters. In addition, we evaluated the performance of five different bio-optical models with different numbers of free parameters for the direct inversion of in-water inherent optical properties (IOPs) from in situ and synthetic R rs data. The multi-parameter models showed similar performances regardless of the number of parameters. Considering the computational cost associated with larger parameter spaces, we recommend bio-optical models with three free parameters for the use of IOP or joint retrieval algorithms.

Assessing the suitability of lakes and reservoirs for recreation using Landsat 8.

Water clarity has long been used as a visual indicator of the condition of water quality. The clarity of waters is generally valued for esthetic and recreational purposes. Water clarity is often assessed using a Secchi disk attached to a measured line and lowered to a depth where it can be no longer seen. We have applied an approach which uses atmospherically corrected Landsat 8 data to estimate the water clarity in freshwater bodies by using the quasi-analytical algorithm (QAA) and Contrast Theory to predict Secchi depths for more than 270 lakes and reservoirs across the continental US. We found that incorporating Landsat 8 spectral data into methodologies created to retrieve the inherent optical properties (IOP) of coastal waters was effective at predicting in situ measures of the clarity of inland water bodies. The predicted Secchi depths were used to evaluate the recreational suitability for swimming and recreation using an assessment framework developed from public perception of water clarity. Results showed approximately 54% of the water bodies in our dataset were classified as "marginally suitable to suitable" with approximately 31% classed as "eminently suitable" and approximately 15% classed as "totally unsuitable-unsuitable". The implications are that satellites engineered for terrestrial applications can be successfully used with traditional ocean color algorithms and methods to measure the water quality of freshwater environments. Furthermore, operational land-based satellite sensors have the temporal repeat cycles, spectral resolution, wavebands, and signal-to-noise ratios to be repurposed to monitor water quality for public use and trophic status of complex inland waters.

Optimal estimation framework for ocean color atmospheric correction and pixel-level uncertainty quantification.

Ocean color (OC) remote sensing requires compensation for atmospheric scattering and absorption (aerosol, Rayleigh, and trace gases), referred to as atmospheric correction (AC). AC allows inference of parameters such as spectrally resolved remote sensing reflectance (Rrs(λ);sr-1) at the ocean surface from the top-of-atmosphere reflectance. Often the uncertainty of this process is not fully explored. Bayesian inference techniques provide a simultaneous AC and uncertainty assessment via a full posterior distribution of the relevant variables, given the prior distribution of those variables and the radiative transfer (RT) likelihood function. Given uncertainties in the algorithm inputs, the Bayesian framework enables better constraints on the AC process by using the complete spectral information compared to traditional approaches that use only a subset of bands for AC. This paper investigates a Bayesian inference research method (optimal estimation [OE]) for OC AC by simultaneously retrieving atmospheric and ocean properties using all visible and near-infrared spectral bands. The OE algorithm analytically approximates the posterior distribution of parameters based on normality assumptions and provides a potentially viable operational algorithm with a reduced computational expense. We developed a neural network RT forward model look-up table-based emulator to increase algorithm efficiency further and thus speed up the likelihood computations. We then applied the OE algorithm to synthetic data and observations from the moderate resolution imaging spectroradiometer (MODIS) on NASA's Aqua spacecraft. We compared the Rrs(λ) retrieval and its uncertainty estimates from the OE method with in-situ validation data from the SeaWiFS bio-optical archive and storage system (SeaBASS) and aerosol robotic network for ocean color (AERONET-OC) datasets. The OE algorithm improved Rrs(λ) estimates relative to the NASA standard operational algorithm by improving all statistical metrics at 443, 555, and 667 nm. Unphysical negative Rrs(λ), which often appears in complex water conditions, was reduced by a factor of 3. The OE-derived pixel-level Rrs(λ) uncertainty estimates were also assessed relative to in-situ data and were shown to have skill.

Satellites quantify the spatial extent of cyanobacterial blooms across the United States at multiple scales.

Previous studies indicate that cyanobacterial harmful algal bloom (cyanoHAB) frequency, extent, and magnitude have increased globally over the past few decades. However, little quantitative capability is available to assess these metrics of cyanoHABs across broad geographic scales and at regular intervals. Here, the spatial extent was quantified from a cyanobacteria algorithm applied to two European Space Agency satellite platforms-the MEdium Resolution Imaging Spectrometer (MERIS) onboard Envisat and the Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3. CyanoHAB spatial extent was defined for each geographic area as the percentage of valid satellite pixels that exhibited cyanobacteria above the detection limit of the satellite sensor. This study quantified cyanoHAB spatial extent for over 2,000 large lakes and reservoirs across the contiguous United States (CONUS) during two time periods: 2008-2011 via MERIS and 2017-2020 via OLCI when cloud-, ice-, and snow-free imagery was available. Approximately 56% of resolvable lakes were glaciated, 13% were headwater, isolated, or terminal lakes, and the rest were primarily drainage lakes. Results were summarized at national-, regional-, state-, and lake-scales, where regions were defined as nine climate regions which represent climatically consistent states. As measured by satellite, changes in national cyanoHAB extent did have a strong increase of 6.9% from 2017 to 2020 (|Kendall's tau (τ)| = 0.56; gamma (γ) = 2.87 years), but had negligible change (|τ| = 0.03) from 2008 to 2011. Two of the nine regions had moderate (0.3 ≤ |τ| < 0.5) increases in spatial extent from 2017 to 2020, and eight of nine regions had negligible (|τ| < 0.2) change from 2008 to 2011. Twelve states had a strong or moderate increase from 2017 to 2020 (|τ| ≥ 0.3), while only one state had a moderate increase and two states had a moderate decrease from 2008 to 2011. A decrease, or no change, in cyanoHAB spatial extent did not indicate a lack of issues related to cyanoHABs. Sensitivity results of randomly omitted daily CONUS scenes confirm that even with reduced data availability during a short four-year temporal assessment, the direction and strength of the changes in spatial extent remained consistent. We present the first set of national maps of lake cyanoHAB spatial extent across CONUS and demonstrate an approach for quantifying past and future changes at multiple spatial scales. Results presented here provide water quality managers information regarding current cyanoHAB spatial extent and quantify rates of change.

Satellites for long-term monitoring of inland U.S. lakes: The MERIS time series and application for chlorophyll-a.

Lakes and other surface fresh waterbodies provide drinking water, recreational and economic opportunities, food, and other critical support for humans, aquatic life, and ecosystem health. Lakes are also productive ecosystems that provide habitats and influence global cycles. Chlorophyll concentration provides a common metric of water quality, and is frequently used as a proxy for lake trophic state. Here, we document the generation and distribution of the complete MEdium Resolution Imaging Spectrometer (MERIS; Appendix A provides a complete list of abbreviations) radiometric time series for over 2300 satellite resolvable inland bodies of water across the contiguous United States (CONUS) and more than 5,000 in Alaska. This contribution greatly increases the ease of use of satellite remote sensing data for inland water quality monitoring, as well as highlights new horizons in inland water remote sensing algorithm development. We evaluate the performance of satellite remote sensing Cyanobacteria Index (CI)-based chlorophyll algorithms, the retrievals for which provide surrogate estimates of phytoplankton concentrations in cyanobacteria dominated lakes. Our analysis quantifies the algorithms' abilities to assess lake trophic state across the CONUS. As a case study, we apply a bootstrapping approach to derive a new CI-to-chlorophyll relationship, ChlBS, which performs relatively well with a multiplicative bias of 1.11 (11%) and mean absolute error of 1.60 (60%). While the primary contribution of this work is the distribution of the MERIS radiometric timeseries, we provide this case study as a roadmap for future stakeholders' algorithm development activities, as well as a tool to assess the strengths and weaknesses of applying a single algorithm across CONUS.

Satellite remote sensing to assess cyanobacterial bloom frequency across the United States at multiple spatial scales.

Cyanobacterial blooms can have negative effects on human health and local ecosystems. Field monitoring of cyanobacterial blooms can be costly, but satellite remote sensing has shown utility for more efficient spatial and temporal monitoring across the United States. Here, satellite imagery was used to assess the annual frequency of surface cyanobacterial blooms, defined for each satellite pixel as the percentage of images for that pixel throughout the year exhibiting detectable cyanobacteria. Cyanobacterial frequency was assessed across 2,196 large lakes in 46 states across the continental United States (CONUS) using imagery from the European Space Agency's Ocean and Land Colour Instrument for the years 2017 through 2019. In 2019, across all satellite pixels considered, annual bloom frequency had a median value of 4% and a maximum value of 100%, the latter indicating that for those satellite pixels, a cyanobacterial bloom was detected by the satellite sensor for every satellite image considered. In addition to annual pixel-scale cyanobacterial frequency, results were summarized at the lake- and state-scales by averaging annual pixel-scale results across each lake and state. For 2019, average annual lake-scale frequencies also had a maximum value of 100%, and Oregon and Ohio had the highest average annual state-scale frequencies at 65% and 52%. Pixel-scale frequency results can assist in identifying portions of a lake that are more prone to cyanobacterial blooms, while lake- and state-scale frequency results can assist in the prioritization of sampling resources and mitigation efforts. Satellite imagery is limited by the presence of snow and ice, as imagery collected in these conditions are quality flagged and discarded. Thus, annual bloom frequencies within nine climate regions were investigated to determine whether missing data biased results in climate regions more prone to snow and ice, given that their annual summaries would be weighted toward the summer months when cyanobacterial blooms tend to occur. Results were unbiased by the time period selected in most climate regions, but a large bias was observed for the Northwest Rockies and Plains climate region. Moderate biases were observed for the Ohio Valley and the Southeast climate regions. Finally, a clustering analysis was used to identify areas of high and low cyanobacterial frequency across CONUS based on average annual lake-scale cyanobacterial frequencies for 2019. Several clusters were identified that transcended state, watershed, and eco-regional boundaries. Combined with additional data, results from the clustering analysis may offer insight regarding large-scale drivers of cyanobacterial blooms.

Bayesian retrieval of optically relevant properties from hyperspectral water-leaving reflectances.

Current methods to retrieve optically relevant properties from ocean color observations do not explicitly make use of prior knowledge about property distributions. Here we implement a simplified Bayesian approach that takes into account prior probability distributions on two sets of five optically relevant parameters, and conduct a retrieval of these parameters using hyperspectral simulated water-leaving reflectances. We focus specifically on the ability of the model to distinguish between two optically similar phytoplankton taxa, diatoms and Noctiluca scintillans. The inversion retrieval gives most-likely concentrations and uncertainty estimates, and we find that the model is able to probabilistically predict the occurrence of Noctiluca scintillans blooms using these metrics. We discuss how this method can be expanded to include a priori covariances between different parameters, and show the effect of varying measurement uncertainty and spectral resolution on Noctiluca scintillans bloom predictions.

Comparing level-2 and level-3 satellite ocean color retrieval validation methodologies.

Many ocean color data applications leverage global spatially composited level-3 (L3) satellite data because of their regular Earth-grid frame of reference. However, ocean color satellite retrieval performance is routinely evaluated on level-2 (L2) data at the native satellite swath resolution and geometries. This study assesses how accurately binned and gridded L3 data represent L2 satellite data products via satellite-to-in situ match-up activities. L2 and L3 satellite data retrievals of the photosynthetic pigment chlorophyll-a are compared with a common in situ dataset, revealing similar L2 and L3 satellite-to-in situ performance for both MODIS-Aqua and VIIRS-SNPP. This agreement implies that L2 validation results are generally applicable to L3 data. However, uncertainties are introduced during the generation of L3 data from L2 data. L3 data comparisons introduce a wider temporal window between the time of in situ measurement and the time of the satellite observation, which can unintentionally reflect on the quality of the satellite retrieval or algorithm performance. The choice of L3 map projection may introduce additional uncertainty by spatially distorting the true location of the satellite retrievals. Each manipulation of satellite data beyond the instrument's native spatiotemporal reference (L2) reduces the applicability of L2 validation results to higher data processing levels.

CHLOROPHYLL ALGORITHMS FOR OCEAN COLOR SENSORS - OC4, OC5 & OC6.

A high degree of consistency and comparability among chlorophyll algorithms is necessary to meet the goals of merging data from concurrent overlapping ocean color missions for increased coverage of the global ocean and to extend existing time series to encompass data from recently launched missions and those planned for the near future, such as PACE, OLCI, HawkEye, EnMAP and SABIA-MAR. To accomplish these goals, we developed 65 empirical ocean color (OC) maximum band ratio (MBR) algorithms for 25 satellite instruments using the largest available and most globally representative database of coincident in situ chlorophyll a and remote sensing reflectances. Excellent internal consistency was achieved across these OC 'Version -7' algorithms, as demonstrated by a median regression slope and coefficient of determination (R2) of 0.985 and 0.859, respectively, between 903 pairwise comparisons of OC-modeled chlorophyll. SeaWiFS and MODIS-Aqua satellite-to-in situ match-up results indicated equivalent, and sometimes superior, performance to current heritage chlorophyll algorithms. During the past forty years of ocean color research the violet band (412 nm) has rarely been used in empirical algorithms to estimate chlorophyll concentrations in oceanic surface water. While the peak in chlorophyll-specific absorption coincides with the 443 nm band present on most ocean color sensors, the magnitude of chlorophyll-specific absorption at 412 nm can reach upwards of ~70% of that at 443 nm. Nearly one third of total chlorophyll-specific absorption between 400 and 700 nm occurs below 443 nm, suggesting that bands below 443 nm, such as the 412 nm band present on most ocean color sensors, may also be useful in detecting chlorophyll under certain conditions and assumptions. The 412 nm band is also the brightest band (that is, with the most dominant magnitude) in remotely sensed reflectances retrieved by heritage passive ocean color instruments when chlorophyll is less than ~0.1 mg m-3, which encompasses ~24% of the global ocean. To attempt to exploit this additional spectral information, we developed two new families of OC algorithms, the OC5 and OC6 algorithms, which include the 412 nm band in the MBR. By using this brightest band in MBR empirical chlorophyll algorithms, the highest possible dynamic range of MBR may be achieved in these oligotrophic areas. The terms oligotrophic, mesotrophic, and eutrophic get frequent use in the scientific literature to designate trophic status; however, quantitative definitions in terms of chlorophyll levels are arbitrarily defined. We developed a new, reproducible, bio-optically based index for trophic status based on the frequency of the brightest, maximum band in the MBR for the OC6_SEAWIFS algorithm, along with remote sensing reflectances from the entire SeaWiFS mission. This index defines oligotrophic water as chlorophyll less than ~0.1 mg m-3, eutrophic water as chlorophyll above 1.67 mg m-3 and mesotrophic water as chlorophyll between 0.1 and 1.67 mg m-3. Applying these criteria to the 40-year mean global ocean chlorophyll data set revealed that oligotrophic, mesotrophic, and eutrophic water occupy ~24%, 67%, and 9%, respectively, of the area of the global ocean on average.

Approach for identifying optically shallow pixels when processing ocean-color imagery.

Ocean-color remote sensing is routinely used to derive marine geophysical parameters from sensor-observed spectral water-leaving radiances. However, in clear geometrically shallow regions, traditional ocean-color algorithms can be confounded by light reflected from the seafloor. Such regions are typically referred to as "optically shallow". When performing spatiotemporal analyses of ocean color datasets, optically shallow features such as submerged sand banks and coral reefs can lead to unexpected regional biases. Most contemporary approaches mask or flag suspected optically shallow pixels based on ancillary bathymetric data. However, the extent to which seafloor reflectance contaminates the water-leaving radiance is dependent on bathymetry, water clarity and seafloor albedo. In this paper, an approach for flagging optically shallow pixels has been developed that considers all three of these variables. In the method, the optical depth of the water column at 547 nm, ζ(547), is predicted from bathymetric data and estimated water-column optical properties. If ζ(547) is less then the pre-defined threshold, a pixel is flagged as potentially optically shallow. Radiative transfer modeling was used to identify a conservative threshold value of ζ(547) = 20 for a bright sand seafloor. In addition, pixels in waters shallower than 5 m are also flagged. We also examined how varying bathymetric datasets may affect the optically shallow flag using MODIS data. It is anticipated that the optically shallow flag will benefit end-users when quality controlling derived ocean color products. Further, the flag may prove useful as a mechanism for switching between optically deep and shallow algorithms during ocean color processing.

Performance metrics for the assessment of satellite data products: an ocean color case study.

Performance assessment of ocean color satellite data has generally relied on statistical metrics chosen for their common usage and the rationale for selecting certain metrics is infrequently explained. Commonly reported statistics based on mean squared errors, such as the coefficient of determination (r2), root mean square error, and regression slopes, are most appropriate for Gaussian distributions without outliers and, therefore, are often not ideal for ocean color algorithm performance assessment, which is often limited by sample availability. In contrast, metrics based on simple deviations, such as bias and mean absolute error, as well as pair-wise comparisons, often provide more robust and straightforward quantities for evaluating ocean color algorithms with non-Gaussian distributions and outliers. This study uses a SeaWiFS chlorophyll-a validation data set to demonstrate a framework for satellite data product assessment and recommends a multi-metric and user-dependent approach that can be applied within science, modeling, and resource management communities.

Retrieval of aerosol properties and water-leaving reflectance from multi-angular polarimetric measurements over coastal waters.

Ocean color remote sensing is an important tool to monitor water quality and biogeochemical conditions of ocean. Atmospheric correction, which obtains water-leaving radiance from the total radiance measured by satellite-borne or airborne sensors, remains a challenging task for coastal waters due to the complex optical properties of aerosols and ocean waters. In this paper, we report a research algorithm on aerosol and ocean color retrieval with emphasis on coastal waters, which uses coupled atmosphere and ocean radiative transfer model to fit polarized radiance measurements at multiple viewing angles and multiple wavelengths. Ocean optical properties are characterized by a generalized bio-optical model with direct accounting for the absorption and scattering of phytoplankton, colored dissolved organic matter (CDOM) and non-algal particles (NAP). Our retrieval algorithm can accurately determine the water-leaving radiance and aerosol properties for coastal waters, and may be used to improve the atmospheric correction when apply to a hyperspectral ocean color instrument.

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.

Mobile device application for monitoring cyanobacteria harmful algal blooms using Sentinel-3 satellite Ocean and Land Colour Instruments.

Cyanobacterial harmful algal blooms (cyanoHAB) cause human and ecological health problems in lakes worldwide. The timely distribution of satellite-derived cyanoHAB data is necessary for adaptive water quality management and for targeted deployment of water quality monitoring resources. Software platforms that permit timely, useful, and cost-effective delivery of information from satellites are required to help managers respond to cyanoHABs. The Cyanobacteria Assessment Network (CyAN) mobile device application (app) uses data from the European Space Agency Copernicus Sentinel-3 satellite Ocean and Land Colour Instrument (OLCI) in near realtime to make initial water quality assessments and quickly alert managers to potential problems and emerging threats related to cyanobacteria. App functionality and satellite data were validated with 25 state health advisories issued in 2017. The CyAN app provides water quality managers with a user-friendly platform that reduces the complexities associated with accessing satellite data to allow fast, efficient, initial assessments across lakes.

Satellite monitoring of cyanobacterial harmful algal bloom frequency in recreational waters and drinking source waters.

Cyanobacterial harmful algal blooms (cyanoHAB) cause extensive problems in lakes worldwide, including human and ecological health risks, anoxia and fish kills, and taste and odor problems. CyanoHABs are a particular concern in both recreational waters and drinking source waters because of their dense biomass and the risk of exposure to toxins. Successful cyanoHAB assessment using satellites may provide an indicator for human and ecological health protection, In this study, methods were developed to assess the utility of satellite technology for detecting cyanoHAB frequency of occurrence at locations of potential management interest. The European Space Agency's MEdium Resolution Imaging Spectrometer (MERIS) was evaluated to prepare for the equivalent series of Sentine1-3 Ocean and Land Colour Imagers (OLCI) launched in 2016 as part of the Copernicus program. Based on the 2012 National Lakes Assessment site evaluation guidelines and National Hydrography Dataset, the continental United States contains 275,897 lakes and reservoirs >1 hectare in area. Results from this study show that 5.6 % of waterbodies were resolvable by satellites with 300 m single-pixel resolution and 0.7 % of waterbodies were resolvable when a three by three pixel (3×3-pixel) array was applied based on minimum Euclidian distance from shore. Satellite data were spatially joined to U.S. public water surface intake (PWSI) locations, where single-pixel resolution resolved 57% of the PWSI locations and a 3×3-pixel array resolved 33% of the PWSI locations. Recreational and drinking water sources in Florida and Ohio were ranked from 2008 through 2011 by cyanoHAB frequency above the World Health Organization's (WHO) high threshold for risk of 100,000 cells mL-1. The ranking identified waterbodies with values above the WHO high threshold, where Lake Apopka, FL (99.1 %) and Grand Lake St. Marys, OH (83 %) had the highest observed bloom frequencies per region. The method presented here may indicate locations with high exposure to cyanoHABs and therefore can be used to assist in prioritizing management resources and actions for recreational and drinking water sources.

Implementation of an analytical Raman scattering correction for satellite ocean-color processing.

Raman scattering of photons by seawater molecules is an inelastic scattering process. This effect can contribute significantly to the water-leaving radiance signal observed by space-borne ocean-color spectroradiometers. If not accounted for during ocean-color processing, Raman scattering can cause biases in derived inherent optical properties (IOPs). Here we describe a Raman scattering correction (RSC) algorithm that has been integrated within NASA's standard ocean-color processing software. We tested the RSC with NASA's Generalized Inherent Optical Properties algorithm (GIOP). A comparison between derived IOPs and in situ data revealed that the magnitude of the derived backscattering coefficient and the phytoplankton absorption coefficient were reduced when the RSC was applied, whilst the absorption coefficient of colored dissolved and detrital matter remained unchanged. Importantly, our results show that the RSC did not degrade the retrieval skill of the GIOP. In addition, a time-series study of oligotrophic waters near Bermuda showed that the RSC did not introduce unwanted temporal trends or artifacts into derived IOPs.

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.

Retrieving marine inherent optical properties from satellites using temperature and salinity-dependent backscattering by seawater.

Time-series of marine inherent optical properties (IOPs) from ocean color satellite instruments provide valuable data records for studying long-term time changes in ocean ecosystems. Semi-analytical algorithms (SAAs) provide a common method for estimating IOPs from radiometric measurements of the marine light field. Most SAAs assign constant spectral values for seawater absorption and backscattering, assume spectral shape functions of the remaining constituent absorption and scattering components (e.g., phytoplankton, non-algal particles, and colored dissolved organic matter), and retrieve the magnitudes of each remaining constituent required to match the spectral distribution of measured radiances. Here, we explore the use of temperature- and salinity-dependent values for seawater backscattering in lieu of the constant spectrum currently employed by most SAAs. Our results suggest that use of temperature- and salinity-dependent seawater spectra elevate the SAA-derived particle backscattering, reduce the non-algal particles plus colored dissolved organic matter absorption, and leave the derived absorption by phytoplankton unchanged.