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.

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.

Atmospheric correction over the ocean for hyperspectral radiometers using multi-angle polarimetric retrievals.

We developed a fast and accurate polynomial based atmospheric correction (POLYAC) algorithm for hyperspectral radiometric measurements, which parameterizes the atmospheric path radiances using aerosol properties retrieved from co-located multi-wavelength multi-angle polarimeter (MAP) measurements. This algorithm has been applied to co-located spectrometer for planetary exploration (SPEX) airborne and research scanning polarimeter (RSP) measurements, where SPEX airborne was used as a proxy of hyperspectral radiometers, and RSP as the MAP. The hyperspectral remote sensing reflectance obtained from POLYAC is accurate when compared to Aerosol Robotic Network (AERONET), and Visible Infrared Imaging Radiometer Suite (VIIRS) ocean color products. POLYAC provides a robust alternative atmospheric correction algorithm for hyperspectral or multi-spectral radiometric measurements for scenes involving coastal oceans and/or absorbing aerosols, where traditional atmospheric correction algorithms are less reliable.

SPEX airborne spectropolarimeter calibration and performance.

To improve our understanding of the complex role of aerosols in the climate system and on air quality, measurements are needed of optical and microphysical aerosol. From many studies, it has become evident that a satellite-based multiangle, multiwavelength polarimeter will be essential to provide such measurements. Here, high accuracy (∼0.003) on the degree of linear polarization (DoLP) measurements is important to retrieve aerosol properties with an accuracy needed to advance our understanding of the aerosol effect on climate. SPEX airborne, a multiangle hyperspectral polarimeter, has been developed for observing and characterizing aerosols from NASA's high-altitude research aircraft ER-2. It delivers measurements of radiance and DoLP at visual wavelengths with a spectral resolution of 3 and 7-30 nm, respectively, for radiance and polarization, at nine fixed equidistant viewing angles from -56° to +56° oriented along the ground track, and a swath of 7° oriented across-track. SPEX airborne uses spectral polarization modulation to determine the state of linear polarization of scattered sunlight. This technique has been developed in the Netherlands and has been demonstrated with ground-based instruments. SPEX airborne serves as a demonstrator for a family of space-based SPEX instruments that have the ability to measure and characterize atmospheric aerosol by multiangle hyperspectral polarimetric imaging remotely from a satellite platform. SPEX airborne was calibrated radiometrically and polarimetrically using Jet Propulsion Laboratory (JPL) facilities including the Polarization Stage Generator-2 (PSG-2), which is designed for polarimetric calibration and validation of the Airborne Multiangle SpectroPolarimetric Imager (AirMSPI). Using the PSG-2, the accuracy of the SPEX airborne DoLP measurements in the laboratory setup is found to be 0.002-0.004. Radiometric calibration is realized with an estimated accuracy of 4%. In 2017, SPEX airborne took part in the "Aerosol Characterization from Polarimeters and Lidar" campaign on the ER-2 that included four polarimeters and two lidars. Polarization measurements of SPEX airborne and the coflying Research Scanning Polarimeter (RSP), recorded during the campaign, were compared and display root-mean-square (RMS) differences ranging from 0.004 (at 555 nm) up to 0.02 (at 410 nm). For radiance measurements, excellent agreement between SPEX airborne and RSP is obtained with an RMS difference of ∼4%. The lab- and flight-performance values for polarization are similar to those recently published for AirMSPI, where also an intercomparison with RSP was made using data from field campaigns in 2013. The intercomparison of radiometric and polarimetric data both display negligible bias. The in-flight comparison results provide verification of SPEX airborne's capability to deliver high-quality data.

Intercomparison of airborne multi-angle polarimeter observations from the Polarimeter Definition Experiment.

In early 2013, three airborne polarimeters were flown on the high altitude NASA ER-2 aircraft in California for the Polarimeter Definition Experiment (PODEX). PODEX supported the pre-formulation NASA Aerosol-Cloud-Ecosystem (ACE) mission, which calls for an imaging polarimeter in polar orbit (among other instruments) for the remote sensing of aerosols, oceans, and clouds. Several polarimeter concepts exist as airborne prototypes, some of which were deployed during PODEX as a capabilities test. Two of those instruments to date have successfully produced Level 1 (georegistered, calibrated radiance and polarization) data from that campaign: the Airborne Multiangle Spectropolarimetric Imager (AirMSPI) and the Research Scanning Polarimeter (RSP). We compared georegistered observations of a variety of scene types by these instruments to test whether Level 1 products agreed within stated uncertainties. Initial comparisons found radiometric agreement, but polarimetric biases beyond measurement uncertainties. After subsequent updates to calibration, georegistration, and the measurement uncertainty models, observations from the instruments now largely agree within stated uncertainties. However, the 470 nm reflectance channels have a roughly +6% bias of AirMSPI relative to RSP, beyond expected measurement uncertainties. We also find that observations of dark (ocean) scenes, where polarimetric uncertainty is expected to be largest, do not agree within stated polarimetric uncertainties. Otherwise, AirMSPI and RSP observations are consistent within measurement uncertainty expectations, providing credibility for the subsequent creation of Level 2 (geophysical product) data from these instruments, and comparison thereof. The techniques used in this work can also form a methodological basis for other intercomparisons, for example, of the data gathered during the recent Aerosol Characterization from Polarimeter and Lidar (ACEPOL) field campaign, carried out in October and November of 2017 with four polarimeters (including AirMSPI and RSP).

Atmospheric Correction of Satellite Ocean-Color Imagery During the PACE Era.

The Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission will carry into space the Ocean Color Instrument (OCI), a spectrometer measuring at 5nm spectral resolution in the ultraviolet (UV) to near infrared (NIR) with additional spectral bands in the shortwave infrared (SWIR), and two multi-angle polarimeters that will overlap the OCI spectral range and spatial coverage, i. e., the Spectrometer for Planetary Exploration (SPEXone) and the Hyper-Angular Rainbow Polarimeter (HARP2). These instruments, especially when used in synergy, have great potential for improving estimates of water reflectance in the post Earth Observing System (EOS) era. Extending the top-of-atmosphere (TOA) observations to the UV, where aerosol absorption is effective, adding spectral bands in the SWIR, where even the most turbid waters are black and sensitivity to the aerosol coarse mode is higher than at shorter wavelengths, and measuring in the oxygen A-band to estimate aerosol altitude will enable greater accuracy in atmospheric correction for ocean color science. The multi-angular and polarized measurements, sensitive to aerosol properties (e.g., size distribution, index of refraction), can further help to identify or constrain the aerosol model, or to retrieve directly water reflectance. Algorithms that exploit the new capabilities are presented, and their ability to improve accuracy is discussed. They embrace a modern, adapted heritage two-step algorithm and alternative schemes (deterministic, statistical) that aim at inverting the TOA signal in a single step. These schemes, by the nature of their construction, their robustness, their generalization properties, and their ability to associate uncertainties, are expected to become the new standard in the future. A strategy for atmospheric correction is presented that ensures continuity and consistency with past and present ocean-color missions while enabling full exploitation of the new dimensions and possibilities. Despite the major improvements anticipated with the PACE instruments, gaps/issues remain to be filled/tackled. They include dealing properly with whitecaps, taking into account Earth-curvature effects, correcting for adjacency effects, accounting for the coupling between scattering and absorption, modeling accurately water reflectance, and acquiring a sufficiently representative dataset of water reflectance in the UV to SWIR. Dedicated efforts, experimental and theoretical, are in order to gather the necessary information and rectify inadequacies. Ideas and solutions are put forward to address the unresolved issues. Thanks to its design and characteristics, the PACE mission will mark the beginning of a new era of unprecedented accuracy in ocean-color radiometry from space.

Calibration and validation of Airborne Multiangle SpectroPolarimetric Imager (AirMSPI) polarization measurements.

The Airborne Multiangle SpectroPolarimetric Imager (AirMSPI), a precursor to the future Multi-Angle Imager for Aerosols satellite instrument, is a remote-sensing instrument for the characterization of atmospheric aerosols and clouds. To help discriminate between different aerosol particle types, which is crucial to improve our understanding of their impact on climate and air quality, AirMSPI acquires imagery over multiple view angles in the ultraviolet, visible, and near-infrared, and it employs dual photoelastic modulators (PEMs) to target an uncertainty requirement of ±0.005 in the degree of linear polarization (DoLP) at selected wavelengths. Laboratory polarimetric calibrations using a second-generation Polarization State Generator-2 (PSG-2) and validation measurements at 0

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.

Remote sensing of aerosols with small satellites in formation flight.

Determination of aerosol optical properties with orbital passive remote sensing is a difficult task, as observations often have limited information. Multi-angle instruments, such as the Multi-angle Imaging SpectroRadiometer (MISR) and the POlarization and Directionality of the Earth's Reflectances (POLDER), seek to address this by making information rich multi-angle observations, which can be used to better retrieve aerosol optical properties. The paradigm for such instruments is that each angle view is made from one platform, with, for example, a gimbaled sensor or multiple fixed view angle sensors. This restricts the observing geometry to a plane within the scene Bidirectional Reflectance Distribution Function (BRDF ) observed at the top of the atmosphere (TOA). New technological developments, however, support sensors on small satellites flying in formation, which could be a beneficial alternative. Such sensors may have only one viewing direction each, but the agility of small satellites allows one to control this direction and change it over time. When such agile satellites are flown in formation and their sensors pointed to the same location at approximately the same time, they could sample a distributed set of geometries within the scene BRDF . In other words, observations from multiple satellites can take a variety of view zenith and azimuth angles, and are not restricted to one azimuth plane as is the case with a single multi-angle instrument. It is not known, however, if this is as potentially capable as a multi-angle platform for the purposes of aerosol remote sensing. Using a systems engineering tool coupled with an information content analysis technique, we investigate the feasibility of such an approach for the remote sensing of aerosols. These tools test the mean results of all geometries encountered in an orbit. We find that small satellites in formation are equally capable as multi-angle platforms for aerosol remote sensing, as long as their calibration accuracies and measurement uncertainties are equivalent. As long as the viewing geometries are dispersed throughout the BRDF , it appears the quantity of view angles determines the information content of the observations, not the specific observation geometry. Given the smoothly varying nature of BRDF 's observed at the TOA, this is reasonable, and supports the viability of aerosol remote sensing with small satellites flying in formation. The incremental improvement in information content that we found with number of view angles also supports the concept of a resilient mission comprised of multiple satellites that are continuously replaced as they age or fail.

Water-leaving contribution to polarized radiation field over ocean.

The top-of-atmosphere (TOA) radiation field from a coupled atmosphere-ocean system (CAOS) includes contributions from the atmosphere, surface, and water body. Atmospheric correction of ocean color imagery is to retrieve water-leaving radiance from the TOA measurement, from which ocean bio-optical properties can be obtained. Knowledge of the absolute and relative magnitudes of water-leaving signal in the TOA radiation field is important for designing new atmospheric correction algorithms and developing retrieval algorithms for new ocean biogeochemical parameters. In this paper we present a systematic sensitivity study of water-leaving contribution to the TOA radiation field, from 340 nm to 865 nm, with polarization included. Ocean water inherent optical properties are derived from bio-optical models for two kinds of waters, one dominated by phytoplankton (PDW) and the other by non-algae particles (NDW). In addition to elastic scattering, Raman scattering and fluorescence from dissolved organic matter in ocean waters are included. Our sensitivity study shows that the polarized reflectance is minimized for both CAOS and ocean signals in the backscattering half plane, which leads to numerical instability when calculating water leaving relative contribution, the ratio between polarized water leaving and CAOS signals. If the backscattering plane is excluded, the water-leaving polarized signal contributes less than 9% to the TOA polarized reflectance for PDW in the whole spectra. For NDW, the polarized water leaving contribution can be as much as 20% in the wavelength range from 470 to 670 nm. For wavelengths shorter than 452 nm or longer than 865 nm, the water leaving contribution to the TOA polarized reflectance is in general smaller than 5% for NDW. For the TOA total reflectance, the water-leaving contribution has maximum values ranging from 7% to 16% at variable wavelengths from 400 nm to 550 nm from PDW. The water leaving contribution to the TOA total reflectance can be as large as 35% for NDW, which is in general peaked at 550 nm. Both the total and polarized reflectances from water-leaving contributions approach zero in the ultraviolet and near infrared bands. These facts can be used as constraints or guidelines when estimating the water leaving contribution to the TOA reflectance for new atmospheric correction algorithms for ocean color imagery.

Analysis of fine-mode aerosol retrieval capabilities by different passive remote sensing instrument designs.

Remote sensing of aerosol optical properties is difficult, but multi-angle, multi-spectral, polarimetric instruments have the potential to retrieve sufficient information about aerosols that they can be used to improve global climate models. However, the complexity of these instruments means that it is difficult to intuitively understand the relationship between instrument design and retrieval success. We apply a Bayesian statistical technique that relates instrument characteristics to the information contained in an observation. Using realistic simulations of fine size mode dominated spherical aerosols, we investigate three instrument designs. Two of these represent instruments currently in orbit: the Multiangle Imaging SpectroRadiometer (MISR) and the POLarization and Directionality of the Earths Reflectances (POLDER). The third is the Aerosol Polarimetry Sensor (APS), which failed to reach orbit during recent launch, but represents a viable design for future instruments. The results show fundamental differences between the three, and offer suggestions for future instrument design and the optimal retrieval strategy for current instruments. Generally, our results agree with previous validation efforts of POLDER and airborne prototypes of APS, but show that the MISR aerosol optical thickness uncertainty characterization is possibly underestimated.

Analysis of shipboard aerosol optical thickness measurements from multiple sunphotometers aboard the R/V Ronald H. Brown during the Aerosol Characterization Experiment--Asia.

Marine sunphotometer measurements collected aboard the R/V Ronald H. Brown during the Aerosol Characterization Experiment-Asia (ACE-Asia) are used to evaluate the ability of complementary instrumentation to obtain the best possible estimates of aerosol optical thickness and Angstrom exponent from ships at sea. A wide range of aerosol conditions, including clean maritime conditions and highly polluted coastal environments, were encountered during the ACE-Asia cruise. The results of this study suggest that shipboard hand-held sunphotometers and fast-rotating shadow-band radiometers (FRSRs) yield similar measurements and uncertainties if proper measurement protocols are used and if the instruments are properly calibrated. The automated FRSR has significantly better temporal resolution (2 min) than the hand-held sunphotometers when standard measurement protocols are used, so it more faithfully represents the variability of the local aerosol structure in polluted regions. Conversely, results suggest that the hand-held sunphotometers may perform better in clean, maritime air masses for unknown reasons. Results also show that the statistical distribution of the Angstrom exponent measurements is different when the distributions from hand-held sunphotometers are compared with those from the FRSR and that the differences may arise from a combination of factors.