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

Informing ocean color inversion products by seeding with ancillary observations

Ocean reflectance inversion algorithms provide many products used in ecological and biogeochemical models. While a number of different inversion approaches exist, they all use only spectral remote-sensing reflectances (Rrs(λ)) as input to derive inherent optical properties (IOPs) in optically deep oceanic waters. However, information content in Rrs(λ) is limited, so spectral inversion algorithms may benefit from additional inputs. Here, we test the simplest possible case of ingesting optical data ('seeding') within an inversion scheme (the Generalized Inherent Optical Property algorithm framework default configuration (GIOP-DC)) with both simulated and satellite datasets of an independently known or estimated IOP, the particulate backscattering coefficient at 532 nm (bbp(532)). We find that the seeded-inversion absorption products are substantially different and more accurate than those generated by the standard implementation. On global scales, seasonal patterns in seeded-inversion absorption products vary by more than 50% compared to absorption from the GIOP-DC. This study proposes one framework in which to consider the next generation of ocean color inversion schemes by highlighting the possibility of adding information collected with an independent sensor.

Development and Validation of an Empirical Ocean Color Algorithm with Uncertainties: A Case Study with the Particulate Backscattering Coefficient

We explored how algorithm (model) and in situ measurement (observation) uncertainties can effectively be incorporated into empirical ocean color model development and assessment. In this study we focused on methods for deriving the particulate backscattering coefficient at 555 nm, b bp (555) (m-1). We developed a simple empirical algorithm for deriving b bp (555) as a function of a remote sensing reflectance line height (LH) metric. Model training was performed using a high-quality bio-optical dataset that contains coincident in situ measurements of the spectral remote sensing reflectances, R rs (λ) (sr-1), and the spectral particulate backscattering coefficients, b bp (λ). The LH metric used is defined as the magnitude of R rs (555) relative to a linear baseline drawn between R rs (490) and R rs (670). Using an independent validation dataset, we compared the skill of the LH-based model with two other models. We used contemporary validation metrics, including bias and mean absolute error (MAE), that were corrected for model and observation uncertainties. The results demonstrated that measurement uncertainties do indeed impact contemporary validation metrics such as mean bias and MAE. Zeta-scores and z-tests for overlapping confidence intervals were also explored as potential methods for assessing model skill.

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.

Evaluating satellite estimates of particulate backscatter in the global open ocean using autonomous profiling floats

Satellite retrievals of particulate backscattering (bbp) are widely used in studies of ocean ecology and biogeochemistry, but have been historically difficult to validate due to the paucity of available ship-based comparative field measurements. Here we present a comparison of satellite and in situ bbp using observations from autonomous floats (n = 2,486 total matchups across three satellites), which provide bbp at 700 nm. With these data, we quantify how well the three inversion products currently distributed by NASA ocean color retrieve bbp. We find that the median ratio of satellite derived bbp to float bbp ranges from 0.77 to 1.60 and Spearman's rank correlations vary from r = 0.06 to r = 0.79, depending on which algorithm and sensor is used. Model skill degrades with increased spatial variability in remote sensing reflectance, which suggests that more rigorous matchup criteria and factors contributing to sensor noisiness may be useful to address in future work, and/or that we have built in biases in the current widely distributed inversion algorithms.

Approach for Propagating Radiometric Data Uncertainties Through NASA Ocean Color Algorithms

Spectroradiometric satellite observations of the ocean are commonly referred to as "ocean color" remote sensing. NASA has continuously collected, processed, and distributed ocean color datasets since the launch of the Sea-viewing Wide-field-of-view Sensor (SeaWiFS) in 1997. While numerous ocean color algorithms have been developed in the past two decades that derive geophysical data products from sensor-observed radiometry, few papers have clearly demonstrated how to estimate measurement uncertainty in derived data products. As the uptake of ocean color data products continues to grow with the launch of new and advanced sensors, it is critical that pixel-by-pixel data product uncertainties are estimated during routine data processing. Knowledge of uncertainties can be used when studying long-term climate records, or to assist in the development and performance appraisal of bio-optical algorithms. In this methods paper we provide a comprehensive overview of how to formulate first-order first-moment (FOFM) calculus for propagating radiometric uncertainties through a selection of bio-optical models. We demonstrate FOFM uncertainty formulations for the following NASA ocean color data products: chlorophyll-a pigment concentration (Chl), the diffuse attenuation coefficient at 490 nm (K _d,490), particulate organic carbon (POC), normalized fluorescent line height (nflh), and inherent optical properties (IOPs). Using a quality-controlled in situ hyperspectral remote sensing reflectance (R _rs,i ) dataset, we show how computationally inexpensive, yet algebraically complex, FOFM calculations may be evaluated for correctness using the more computationally expensive Monte Carlo approach. We compare bio-optical product uncertainties derived using our test R _rs dataset assuming spectrally-flat, uncorrelated relative uncertainties of 1, 5, and 10%. We also consider spectrally dependent, uncorrelated relative uncertainties in R _rs . The importance of considering spectral covariances in R _rs , where practicable, in the FOFM methodology is highlighted with an example SeaWiFS image. We also present a brief case study of two POC algorithms to illustrate how FOFM formulations may be used to construct measurement uncertainty budgets for ecologically-relevant data products. Such knowledge, even if rudimentary, may provide useful information to end-users when selecting data products or when developing their own algorithms.

Sensitivity of inherent optical properties from ocean reflectance inversion models to satellite instrument wavelength suites

The Earth science community seeks to develop climate data records (CDRs) from satellite measurements of ocean color, a continuous data record which now exceeds twenty years. Space agencies will launch additional instruments in the coming decade that will continue this data record, including the NASA PACE spectrometer. Inherent optical properties (IOPs) quantitatively describe the absorbing and scattering constituents of seawater and can be estimated from satellite-observed spectroradiometric data using semi-analytical algorithms (SAAs). SAAs exploit the contrasting optical signatures of constituent matter at spectral bands observed by satellite sensors. SAA performance, therefore, depends on the spectral resolution of the satellite spectroradiometer. A CDR spanning SeaWiFS, MODIS, OLCI, and PACE, for example, would include IOPs derived using varied wavelength suites if all available wavelengths were considered. Here, we explored differences in derived IOPs that stem simply from the use of (eight) different wavelength suites of input radiometric measurements. Using synthesized data and SeaWiFS Level-3 mission-long composites, we demonstrated equivalent SAA performance for all wavelength suites, but that IOP retrievals vary by several percent across wavelength suites and as a function of water type. The differences equate to roughly ≤ 6, 12, and 7% for a_dg (443), a_ph (443), and b_bp (443), respectively, for waters with C_a ≤ 1 mg m^-3. These values shrink for sensors with similar wavelength suites (e.g., SeaWiFS, MODIS, and MERIS) and rise to substantially larger values for higher C_a waters. Our results also indicate that including 400 nm (in the case of OLCI) influences the derived IOPs, using longer wavelengths (>600 nm) influences the derived IOPs when there is a red signal, and, including additional spectral information shows potential for improved IOP estimation, but not without revisiting SAA parameterizations and execution. While modest in scope, we believe this study contributes to the knowledge base for CDR development. The implication of ignoring such an analysis as CDRs continue to be developed is a prolonged inability to distinguish between algorithmic and environmental contributions to trends and anomalies in the IOP time-series.

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.