An optimization approach has been developed for simultaneous retrieval of aerosol properties and normalized water-leaving radiance (nLw) from multispectral, multiangular, and polarimetric observations over ocean. The main features of the method are (1) use of a simplified bio-optical model to estimate nLw, followed by an empirical refinement within a specified range to improve its accuracy; (2) improved algorithm convergence and stability by applying constraints on the spatial smoothness of aerosol loading and Chlorophyll a (Chl a) concentration across neighboring image patches and spectral constraints on aerosol optical properties and nLw across relevant bands; and (3) enhanced Jacobian calculation by modeling and storing the radiative transfer (RT) in aerosol/Rayleigh mixed layer, pure Rayleigh-scattering layers, and ocean medium separately, then coupling them to calculate the field at the sensor. This approach avoids unnecessary and time-consuming recalculations of RT in unperturbed layers in Jacobian evaluations. The Markov chain method is used to model RT in the aerosol/Rayleigh mixed layer and the doubling method is used for the uniform layers of the atmosphere–ocean system. Our optimization approach has been tested using radiance and polarization measurements acquired by the Airborne Multiangle SpectroPolarimetric Imager (AirMSPI) over the AERONET USC_SeaPRISM ocean site (6 February 2013) and near the AERONET La Jolla site (14 January 2013), which, respectively, reported relatively high and low aerosol loadings. Validation of the results is achieved through comparisons to AERONET aerosol and ocean color products. For comparison, the USC_SeaPRISM retrieval is also performed by use of the Generalized Retrieval of Aerosol and Surface Properties algorithm (Dubovik et al., 2011). Uncertainties of aerosol and nLw retrievals due to random and systematic instrument errors are analyzed by truth-in/truth-out tests with three Chl a concentrations, five aerosol loadings, three different types of aerosols, and nine combinations of solar incidence and viewing geometries.

Tomography aims to recover a three-dimensional (3D) density map of a medium or an object. In medical imaging, it is extensively used for diagnostics via X-ray computed tomography (CT). We define and derive a tomography of cloud droplet distributions via passive remote sensing. We use multi-view polarimetric images to fit a 3D polarized radiative transfer (RT) forward model. Our motivation is 3D volumetric probing of vertically-developed convectively-driven clouds that are ill-served by current methods in operational passive remote sensing. Current techniques are based on strictly 1D RT modeling and applied to a single cloudy pixel, where cloud geometry defaults to that of a plane-parallel slab. Incident unpolarized sunlight, once scattered by cloud-droplets, changes its polarization state according to droplet size. Therefore, polarimetric measurements in the rainbow and glory angular regions can be used to infer the droplet size distribution. This work defines and derives a framework for a full 3D tomography of cloud droplets for both their mass concentration in space and their distribution across a range of sizes. This 3D retrieval of key microphysical properties is made tractable by our novel approach that involves a restructuring and differentiation of an open-source polarized 3D RT code to accommodate a special two-step optimization technique. Physically-realistic synthetic clouds are used to demonstrate the methodology with rigorous uncertainty quantification.

This paper presents the physical basis of the Earth Polychromatic Imaging Camera (EPIC) cloud product algorithms and an initial evaluation of their performance. Since June 2015, EPIC has been providing observations of the sunlit side of the Earth with its 10 spectral channels ranging from the UV to the near-infrared. A suite of algorithms has been developed to generate the standard EPIC Level 2 cloud products that include cloud mask, cloud effective pressure/height, and cloud optical thickness. The EPIC cloud mask adopts the threshold method and utilizes multichannel observations and ratios as tests. Cloud effective pressure/height is derived with observations from the O₂ A-band (780 and 764 nm) and B-band (680 and 688 nm) pairs. The EPIC cloud optical thickness retrieval adopts a single-channel approach in which the 780 and 680 nm channels are used for retrievals over ocean and over land, respectively. Comparison with co-located cloud retrievals from geosynchronous earth orbit (GEO) and low earth orbit (LEO) satellites shows that the EPIC cloud product algorithms are performing well and are consistent with theoretical expectations. These products are publicly available at the Atmospheric Science Data Center at the NASA Langley Research Center for climate studies and for generating other geophysical products that require cloud properties as input.