Contribution of water-leaving radiances to multiangle, multispectral polarimetric observations over the open ocean: bio-optical model results for case 1 waters

Particle backscattering efficiency as an indicator of suspended particle composition and size in the Pearl River Estuary

Particle backscattering efficiency (Qbbe(λ)) relates closely to particle type and links light backscattering with biogeochemical properties. However, knowledge gaps about Qbbe(λ) persist in the Pearl River Estuary (PRE). We measured Qbbe(λ) within the three transects in the PRE and investigated its relationships with particle composition and size. The results show that particle size influences Qbbe(λ) more than composition. Vertical heterogeneity of Qbbe(λ) was observed; high Qbbe(λ) values spanned vertical columns due to dynamic river-sea interactions, showing distinct layers with higher concentrations at the surface and bottom, but lower mid-layer values unaffected by rivers. Transects influenced by river discharge showed elevated Qbbe(λ) due to high inorganic particle inputs, while coastal currents influenced bottom values via sediment resuspension. This enhanced understanding of Qbbe(λ)'s spatial-temporal variations in the South China Sea is vital for better predicting particle characteristics from optical data.

Scattering direction sampling methods for polarized Monte Carlo simulation of oceanic lidar

Monte Carlo techniques have been widely applied in polarized light simulation. Based on different preconditions, there are two main types of sampling strategies for scattering direction: one is the scalar sampling method; the others are polarized sampling approaches, including the one- and two-point rejection methods. The polarized simulation of oceanic lidar involves a variety of mediums, and an efficient scattering sampling method is the basis for the coupling simulation of the atmosphere and ocean. To determine the optimal scattering sampling method for oceanic lidar simulation, we developed a polarized Monte Carlo model and simulated Mie scattering, Rayleigh scattering, and Petzold average-particle scattering experiments. This simulation model has been validated by comparison with Ramella-Roman's program [Opt. Express13, 4420 (2005)OPEXFF1094-408710.1364/OPEX.13.004420], with differences in reflectance and transmittance Stokes less than 1% in Mie scattering. The simulation results show these scattering sampling methods differ in runtime, scattering angle distributions, and reflectance and transmittance Stokes. Considering the current simulation accuracy of oceanic lidar, the differences in reflectance and transmittance Stokes are acceptable; thus, the runtime becomes the main evaluation factor. The one-point rejection method and scalar sampling method are preferable for the oceanic lidar polarized simulation. Under complex atmosphere-ocean coupling systems, scalar sampling methods may be a better choice since the calculation process of the sampling is independent of the incident Stokes vector.

Atmospheric diffuse transmittance of the linear polarization component of water-leaving radiation

The polarization characteristics of water-leaving radiation contain rich information on oceanic constituents. Determining the atmospheric diffuse transmittance is crucial for extracting the polarization information of water-leaving radiation from the radiation acquired by polarimetry satellites at the top of the atmosphere. However, there is still a lack of understanding of the atmospheric diffuse transmittance of the linear polarization component of water-leaving radiation. Here, we first evaluated the difference between the atmospheric diffuse transmittance of the linear polarization component (T_Q, T_U) and the intensity component (T_I) of the water-leaving radiation based on the Ocean Successive Orders with Atmosphere Advanced radiative transfer model. As a consequence, there were apparent differences between T_Q, T_U and T_I. In the case of a large solar zenith angle and a large viewing zenith angle, the difference between T_Q, T_U and T_I will exceed 1. Meanwhile, compared with T_I, the oceanic constituents had a prominent interference with T_Q and T_U, and the sediment concentration had little interference with T_Q and T_U in low- and medium-turbidity water with respect to the aerosol model, optical thickness, observation geometry, and phytoplankton. Moreover, T_Q and T_U lookup tables were generated for medium- and low-turbidity water, which laid the foundation for extracting the water-leaving radiation polarization information from the satellite observation polarization signal.

Effects of water salinity on the multi-angular polarimetric properties of light reflected from smooth water surfaces

Salinity is an important environmental factor regulating the aquatic system structure of lakes and other water bodies. Changes in salinity, which can be caused by human activities, can adversely impact the life of water organisms. The refractive index, which can be directly related to water salinity, also controls the polarimetric properties of light reflected from the water surface. In this study, polarimetric measurements of smooth water surfaces with different salinity content were performed at different viewing zenith angles in the wavelength range of 450-1000 nm in the specular reflection directions. The results show that the light reflected from the water surface (defined as reflectance factor) in one measurement direction can be replaced by the reflectance factor derived from polarimetric measurements, and if the polarizer absorptance is considered, the average relative difference is less than 3%. The degree of linear polarization (DOLP) was used to retrieve the refractive indices of water with different salinities based on the Fresnel reflection coefficient. The inverted refractive indices not only have high accuracy (uncertainty from 0.9% to 1.8%) but also have a very strong relationship with the water salinity content. Our study shows the possibility of estimating the variation in water salinity using multi-angular polarimetric measurements.

Use of A Neural Network-Based Ocean Body Radiative Transfer Model for Aerosol Retrievals from Multi-Angle Polarimetric Measurements

For aerosol retrieval from multi-angle polarimetric (MAP) measurements over the ocean it is important to accurately account for the contribution of the ocean-body to the top-of-atmosphere signal, especially for wavelengths <500 nm. Performing online radiative transfer calculations in the coupled atmosphere ocean system is too time consuming for operational retrieval algorithms. Therefore, mostly lookup-tables of the ocean body reflection matrix are used to represent the lower boundary in an atmospheric radiative transfer model. For hyperspectral measurements such as those from Spectro-Polarimeter for Planetary Exploration (SPEXone) on the NASA Plankton, Aerosol, Cloud and ocean Ecosystem (PACE) mission, also the use of look-up tables is unfeasible because they will become too big. In this paper, we propose a new method for aerosol retrieval over ocean from MAP measurements using a neural network (NN) to model the ocean body reflection matrix. We apply the NN approach to synthetic SPEXone measurements and also to real data collected by SPEX airborne during the Aerosol Characterization from Polarimeter and Lidar (ACEPOL) campaign. We conclude that the NN approach is well capable for aerosol retrievals over ocean, introducing no significant error on the retrieved aerosol properties

Sensitivity study on the contribution of scattering by randomly oriented nonspherical hydrosols to linear polarization in clear to semi-turbid shallow waters

The influence of hydrosol nonsphericity on the polarization characteristics of light under water is investigated by combining accurate single-scattering models for randomly oriented spheroidal scatterers with a radiative transfer model that employs Stokes formalism and considers refraction of direct unpolarized solar radiation and 100% linearly polarized radiation at the air-water interface followed by single scattering. Variations in what we call the "linear polarization phase function" (the degree of linear polarization as a function of scattering angle and the angle of linear polarization as a function of scattering angle) are examined for a wide range of spheroid aspect ratios and complex refractive indices of hydrosols. Implications for polarization-sensitive marine organisms and for remote sensing of the marine environment are discussed.

Characterization of radiance from the ocean surface by hyperspectral imaging

A novel snapshot hyperspectral imager is introduced for ocean color (OC) applications and its capabilities are demonstrated. The instrument provides hyperspectral radiance images with a wide field-of-view (FOV) and short exposure time, which is valuable for the direct characterization of the wind-roughened surface in various illumination conditions and wind speeds. Uncertainties in the total(L_t), sky (L_s) and derived water-leaving (L_w)radiances at viewing angles of 20-60° are determined as a function of wind speed together with associated correlation coefficients and variances of the sea surface reflectance coefficient ρ. Estimated L_w uncertainties can partially explain the inaccuracy of satellite retrievals in the blue bands in the coastal waters. It is shown that in above-water measurements in no-glint conditions with viewing and azimuth angles of 40° and 90°, respectively, for both L_t(λ) and L_s(λ) the impact of FOV is minimal at least up to measured W = 5.7 m/s for full-angle FOV of 4° and larger. Implications of uncertainties for the derivation of water leaving radiance in above-water ship-borne and AERONET-OC measurements are discussed.

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.

Measuring and Modeling the Polarized Upwelling Radiance Distribution in Clear and Coastal Waters

The upwelling spectral radiance distribution is polarized, and this polarization varies with the optical properties of the water body. Knowledge of the polarized, upwelling, bidirectional radiance distribution function (BRDF) is important for generating consistent, long-term data records for ocean color because the satellite sensors from which the data are derived are sensitive to polarization. In addition, various studies have indicated that measurement of the polarization of the radiance leaving the ocean can used to determine particle characteristics (Tonizzo et al., 2007; Ibrahim et al., 2016; Chami et al., 2001). Models for the unpolarized BRDF (Morel et al., 2002; Lee et al., 2011) have been validated (Voss et al., 2007; Gleason et al., 2012), but variations in the polarization of the upwelling radiance due to the sun angle, viewing geometry, dissolved material, and suspended particles have not been systematically documented. In this work, we simulated the upwelling radiance distribution using a Monte Carlo-based radiative transfer code and measured it using a set of fish-eye cameras with linear polarizing filters. The results of model-data comparisons from three field experiments in clear and turbid coastal conditions showed that the degree of linear polarization (DOLP) of the upwelling light field could be determined by the model with an absolute error of ±0.05 (or 5% when the DOLP was expressed in %). This agreement was achieved even with a fixed scattering Mueller matrix, but did require in situ measurements of the other inherent optical properties, e.g., scattering coefficient, absorption coefficient, etc. This underscores the difficulty that is likely to be encountered using the particle scattering Mueller matrix (as indicated through the remote measurement of the polarized radiance) to provide a signature relating to the properties of marine particles beyond the attenuation/absorption coefficient.

Variability of the reflectance coefficient of skylight from the ocean surface and its implications to ocean color

The value and spectral dependence of the reflectance coefficient (ρ) of skylight from wind-roughened ocean surfaces is critical for determining accurate water leaving radiance and remote sensing reflectances from shipborne, AERONET-Ocean Color and satellite observations. Using a vector radiative transfer code, spectra of the reflectance coefficient and corresponding radiances near the ocean surface and at the top of the atmosphere (TOA) are simulated for a broad range of parameters including flat and windy ocean surfaces with wind speeds up to 15 m/s, aerosol optical thicknesses of 0-1 at 440nm, wavelengths of 400-900 nm, and variable Sun and viewing zenith angles. Results revealed a profound impact of the aerosol load and type on the spectral values of ρ. Such impacts, not included yet in standard processing, may produce significant inaccuracies in the reflectance spectra retrieved from above-water radiometry and satellite observations. Implications for satellite cal/val activities as well as potential changes in measurement and data processing schemes are discussed.

Simultaneous polarimeter retrievals of microphysical aerosol and ocean color parameters from the "MAPP" algorithm with comparison to high-spectral-resolution lidar aerosol and ocean products

We present an optimal-estimation-based retrieval framework, the microphysical aerosol properties from polarimetry (MAPP) algorithm, designed for simultaneous retrieval of aerosol microphysical properties and ocean color bio-optical parameters using multi-angular total and polarized radiances. Polarimetric measurements from the airborne NASA Research Scanning Polarimeter (RSP) were inverted by MAPP to produce atmosphere and ocean products. The RSP MAPP results are compared with co-incident lidar measurements made by the NASA High-Spectral-Resolution Lidar HSRL-1 and HSRL-2 instruments. Comparisons are made of the aerosol optical depth (AOD) at 355 and 532 nm, lidar column-averaged measurements of the aerosol lidar ratio and Ångstrøm exponent, and lidar ocean measurements of the particulate hemispherical backscatter coefficient and the diffuse attenuation coefficient. The measurements were collected during the 2012 Two-Column Aerosol Project (TCAP) campaign and the 2014 Ship-Aircraft Bio-Optical Research (SABOR) campaign. For the SABOR campaign, 73% RSP MAPP retrievals fall within ±0.04 AOD at 532 nm as measured by HSRL-1, with an R value of 0.933 and root-mean-square deviation of 0.0372. For the TCAP campaign, 53% of RSP MAPP retrievals are within 0.04 AOD as measured by HSRL-2, with an R value of 0.927 and root-mean-square deviation of 0.0673. Comparisons with HSRL-2 AOD at 355 nm during TCAP result in an R value of 0.959 and a root-mean-square deviation of 0.0694. The RSP retrievals using the MAPP optimal estimation framework represent a key milestone on the path to a combined lidar+polarimeter retrieval using both HSRL and RSP measurements.

Airborne and shipborne polarimetric measurements over open ocean and coastal waters: intercomparisons and implications for spaceborne observations

Comprehensive polarimetric closure is demonstrated using observations from two in-situ polarimeters and Vector Radiative Transfer (VRT) modeling. During the Ship-Aircraft Bio-Optical Research (SABOR) campaign, the novel CCNY HyperSAS-POL polarimeter was mounted on the bow of the R/V Endeavor and acquired hyperspectral measurements from just above the surface of the ocean, while the NASA GISS Research Scanning Polarimeter was deployed onboard the NASA LaRC's King Air UC-12B aircraft. State-of-the-art, ancillary measurements were used to characterize the atmospheric and marine contributions in the VRT model, including those of the High Spectral Resolution Lidar (HSRL), the AErosol RObotic NETwork for Ocean Color (AERONET-OC), a profiling WETLabs ac-9 spectrometer and the Multi-spectral Volume Scattering Meter (MVSM). An open-ocean and a coastal scene are analyzed, both affected by complex aerosol conditions. In each of the two cases, it is found that the model is able to accurately reproduce the Stokes components measured simultaneously by each polarimeter at different geometries and viewing altitudes. These results are mostly encouraging, considering the different deployment strategies of RSP and HyperSAS-POL, which imply very different sensitivities to the atmospheric and ocean contributions, and open new opportunities in above-water polarimetric measurements. Furthermore, the signal originating from each scene was propagated to the top of the atmosphere to explore the sensitivity of polarimetric spaceborne observations to changes in the water type. As expected, adding polarization as a measurement capability benefits the detection of such changes, reinforcing the merits of the full-Stokes treatment in modeling the impact of atmospheric and oceanic constituents on remote sensing observations.

Estimation of chlorophyll concentration in waters near Hokkaido using the linear combination method

An inversion algorithm is implemented to retrieve the surface chlorophyll a (Chl) concentration using satellite observation data from the MODIS instrument. The algorithm employs a simple and flexible index (LCI) to combine with the Chl without explicit correction for aerosol scattering. To investigate the sensitivity of LCI to Chl and other influence factors, an oceanic radiative transfer model coupled with a comprehensive bio-optical module is developed. It is studied that the LCI is significantly linear to Chl and not sensitive to other influence factors, except in very low oceanic salinity or scattering angle conditions, where over a 12% relative difference of derived Chl exists. Inversion results show the retrieved Chl are highly consistent with the MODIS operational data products in waters near Hokkaido, with the correlation coefficient, root mean square deviation, and average percentage difference of 0.9702, 0.3756 mg m^-3, and 13.89%, respectively. Investigation of the validity of this algorithm with a variety of atmospheric conditions indicates that the residual influence of atmosphere on the LCI index, after Rayleigh scattering correction, is generally within ± 0.001, allowing the retrieval error of Chl at less than 25% in most cases. A good comparison between retrieval and in situ measurements is also identified and implies that the retrieval accuracy via the LCI method depends on the linear combination coefficients used and the bio-optical module selected, while effects of polarization can be ignored.

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.

Equivalence of internal and external mixture schemes of single scattering properties in vector radiative transfer

Polarized radiation fields in a turbid medium are influenced by single-scattering properties of scatterers. It is common that media contain two or more types of scatterers, which makes it essential to properly mix single-scattering properties of different types of scatterers in the vector radiative transfer theory. The vector radiative transfer solvers can be divided into two basic categories: the stochastic and deterministic methods. The stochastic method is basically the Monte Carlo method, which can handle scatterers with different scattering properties explicitly. This mixture scheme is called the external mixture scheme in this paper. The deterministic methods, however, can only deal with a single set of scattering properties in the smallest discretized spatial volume. The single-scattering properties of different types of scatterers have to be averaged before they are input to deterministic solvers. This second scheme is called the internal mixture scheme. The equivalence of these two different mixture schemes of scattering properties has not been demonstrated so far. In this paper, polarized radiation fields for several scattering media are solved using the Monte Carlo and successive order of scattering (SOS) methods and scattering media contain two types of scatterers: Rayleigh scatterers (molecules) and Mie scatterers (aerosols). The Monte Carlo and SOS methods employ external and internal mixture schemes of scatterers, respectively. It is found that the percentage differences between radiances solved by these two methods with different mixture schemes are of the order of 0.1%. The differences of Q/I, U/I, and V/I are of the order of 10^-5∼10^-4, where I, Q, U, and V are the Stokes parameters. Therefore, the equivalence between these two mixture schemes is confirmed to the accuracy level of the radiative transfer numerical benchmarks. This result provides important guidelines for many radiative transfer applications that involve the mixture of different scattering and absorptive particles.

Vector radiative transfer model for coupled atmosphere and ocean systems including inelastic sources in ocean waters

Inelastic scattering plays an important role in ocean optics. The main inelastic scattering mechanisms include Raman scattering, fluorescence by colored dissolved organic matter (FDOM), and fluorescence by chlorophyll. This paper reports an implementation of all three inelastic scattering mechanisms in the exact vector radiative transfer model for coupled atmosphere and ocean Systems (CAOS). Simulation shows that FDOM contributes to the water radiation field in the broad visible spectral region, while chlorophyll fluorescence is limited in a narrow band centered at 685 nm. This is consistent with previous findings in the literature. The fluorescence distribution as a function of depth and viewing angle is presented. The impacts of fluorescence to the degree of linear polarization (DoLP) and orientation of the polarization ellipse (OPE) are studied. The DoLP is strongly influenced by inelastic scattering at wavelengths with strong inelastic scattering contribution. The OPE is less affected by inelastic scattering but it has a noticeable impact, in terms of the angular region of positive polarization, in the backward direction. This effect is more apparent for deeper water depth.

Polarized transfer functions of the ocean surface for above-surface determination of the vector submarine light field

A method is developed to determine the underwater polarized light field from above sea surface observations. A hybrid approach combining vector radiative transfer simulations and the Monte Carlo method is used to determine the transfer functions of polarized light for wind-driven ocean surfaces. Transfer functions for surface-reflected skylight and upward transmission of light through the sea surface are presented for many common viewing and solar geometries for clear-sky conditions. Sensitivity of reflection matrices to environmental conditions is examined and can vary up to 50% due to wind speed, 25% due to atmospheric aerosol load, and 10% due to radiometer field-of-view. Scalar transmission is largely independent of water type and varies a few percent with wind speed, while polarized components can change up to 10% in high winds. Considerations for determining the water-leaving radiance (scalar or vector) are discussed.

Radiative transfer simulations of the two-dimensional ocean glint reflectance and determination of the sea surface roughness

An optimized discrete-ordinate radiative transfer model (DISORT3) with a pseudo-two-dimensional bidirectional reflectance distribution function (BRDF) is used to simulate and validate ocean glint reflectances at an infrared wavelength (1036 nm) by matching model results with a complete set of BRDF measurements obtained from the NASA cloud absorption radiometer (CAR) deployed on an aircraft. The surface roughness is then obtained through a retrieval algorithm and is used to extend the simulation into the visible spectral range where diffuse reflectance becomes important. In general, the simulated reflectances and surface roughness information are in good agreement with the measurements, and the diffuse reflectance in the visible, ignored in current glint algorithms, is shown to be important. The successful implementation of this new treatment of ocean glint reflectance and surface roughness in DISORT3 will help improve glint correction algorithms in current and future ocean color remote sensing applications.

OSOAA: a vector radiative transfer model of coupled atmosphere-ocean system for a rough sea surface application to the estimates of the directional variations of the water leaving reflectance to better process multi-angular satellite sensors data over the ocean

In this study, we present a radiative transfer model, so-called OSOAA, that is able to predict the radiance and degree of polarization within the coupled atmosphere-ocean system in the presence of a rough sea surface. The OSOAA model solves the radiative transfer equation using the successive orders of scattering method. Comparisons with another operational radiative transfer model showed a satisfactory agreement within 0.8%. The OSOAA model has been designed with a graphical user interface to make it user friendly for the community. The radiance and degree of polarization are provided at any level, from the top of atmosphere to the ocean bottom. An application of the OSOAA model is carried out to quantify the directional variations of the water leaving reflectance and degree of polarization for phytoplankton and mineral-like dominated waters. The difference between the water leaving reflectance at a given geometry and that obtained for the nadir direction could reach 40%, thus questioning the Lambertian assumption of the sea surface that is used by inverse satellite algorithms dedicated to multi-angular sensors. It is shown as well that the directional features of the water leaving reflectance are weakly dependent on wind speed. The quantification of the directional variations of the water leaving reflectance obtained in this study should help to correctly exploit the satellite data that will be acquired by the current or forthcoming multi-angular satellite sensors.

Polarized reflectance and transmittance properties of windblown sea surfaces

Generation of random sea surfaces using wave variance spectra and Fourier transforms is formulated in a way that guarantees conservation of wave energy and fully resolves wave height and slope variances. Monte Carlo polarized ray tracing, which accounts for multiple scattering between light rays and wave facets, is used to compute effective Mueller matrices for reflection and transmission of air- or water-incident polarized radiance. Irradiance reflectances computed using a Rayleigh sky radiance distribution, sea surfaces generated with Cox-Munk statistics, and unpolarized ray tracing differ by 10%-18% compared with values computed using elevation- and slope-resolving surfaces and polarized ray tracing. Radiance reflectance factors, as used to estimate water-leaving radiance from measured upwelling and sky radiances, are shown to depend on sky polarization, and improved values are given.

Polarization impacts on the water-leaving radiance retrieval from above-water radiometric measurements

Above-water measurements of water-leaving radiance are widely used for water-quality monitoring and ocean-color satellite data validation. Reflected skylight in above-water radiometry needs to be accurately estimated prior to derivation of water-leaving radiance. Up-to-date methods to estimate reflection of diffuse skylight on rough sea surfaces are based on radiative transfer simulations and sky radiance measurements. But these methods neglect the polarization state of the incident skylight, which is generally highly polarized. In this paper, the effects of polarization on the sea surface reflectance and the subsequent water-leaving radiance estimation are investigated. We show that knowledge of the polarization field of the diffuse skylight significantly improves above-water radiometry estimates, in particular in the blue part of the spectrum where the reflected skylight is dominant. A newly developed algorithm based on radiative transfer simulations including polarization is described. Its application to the standard Aerosol Robotic Network-Ocean Color and hyperspectral radiometric measurements of the 1.5-year dataset acquired at the Long Island Sound site demonstrates the noticeable importance of considering polarization for water-leaving radiance estimation. In particular it is shown, based on time series of collocated data acquired in coastal waters, that the azimuth range of measurements leading to good-quality data is significantly increased, and that these estimates are improved by more than 12% at 413 nm. Full consideration of polarization effects is expected to significantly improve the quality of the field data utilized for satellite data validation or potential vicarious calibration purposes.

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.

Observation of non-principal plane neutral points in the in-water upwelling polarized light field

Neutral points are specific directions in the light field where the three Stokes parameters Q, U, V, and thus the degree of polarization simultaneously go to zero. We have made the first measurement of non-principal-plane neutral points in the upwelling light field in natural waters. These neutral points are located at approximately 40°- 80° nadir angle and between 120° - 160° azimuth to the sun which is well off of the principal plane. Calculations show that the neutral point positions are very sensitive to the balance in the incident light between the partially polarized skylight and the direct solar beam.

Polarized radiance fields under a dynamic ocean surface: a three-dimensional radiative transfer solution

The hybrid matrix operator, Monte Carlo (HMOMC) method previously reported [Appl. Opt.47, 1063-1071 (2008)] is improved by neglecting higher-order terms in the coupling of the matrix operators and by introducing a dual grid scheme. The computational efficiency for solving the vector radiative transfer equation in a full 3D coupled atmosphere-surface-ocean system is substantially improved, and, thus, large-scale simulations of the radiance distribution become feasible. The improved method is applied to the computation of the polarized radiance field under realistic surface waves simulated by the power spectral density method. To the authors' best knowledge, this is the first time that the polarized radiance field under a dynamic ocean surface and the underwater image of an object above such an ocean surface have been reported.

Polarized light in coastal waters: hyperspectral and multiangular analysis

Measurements of the underwater polarized light field were performed at different stations, atmospheric conditions and water compositions using a newly developed hyperspectral and multiangular polarimeter during a recent cruise in the coastal areas of New York Harbor - Sandy Hook, NJ region (USA). Results are presented for waters with chlorophyll concentrations 1.3-4.8 microg/l and minerals concentrations 2.0- 3.9 mg/l. Angular and spectral variations of the degree of polarization are found to be consistent with theory. Maximum values of the degree of polarization do not exceed 0.4 and the position of the maximum is close to 100 masculine scattering angle. Normalized radiances and degrees of polarization are compared with simulated ones obtained with a Monte Carlo radiative transfer code for the atmosphere-ocean system and show satisfactory agreement.

Near-infrared active polarimetric and multispectral laboratory demonstrator for target detection

We report on the design and exploitation of a real-field laboratory demonstrator combining active polarimetric and multispectral functions. Its building blocks, including a multiwavelength pulsed optical parametric oscillator at the emission side and a hyperspectral imager with polarimetric capability at the reception side, are described. The results obtained with this demonstrator are illustrated on some examples and discussed. In particular it is found that good detection performances rely on joint use of intensity and polarimetric images, with these images exhibiting complementary signatures in most cases.

A vector radiative transfer model for coupled atmosphere and ocean systems based on successive order of scattering method

A vector radiative transfer model has been developed for coupled atmosphere and ocean systems based on the Successive Order of Scattering (SOS) Method. The emphasis of this study is to make the model easy-to-use and computationally efficient. This model provides the full Stokes vector at arbitrary locations which can be conveniently specified by users. The model is capable of tracking and labeling different sources of the photons that are measured, e.g. water leaving radiances and reflected sky lights. This model also has the capability to separate florescence from multi-scattered sunlight. The delta - fit technique has been adopted to reduce computational time associated with the strongly forward-peaked scattering phase matrices. The exponential - linear approximation has been used to reduce the number of discretized vertical layers while maintaining the accuracy. This model is developed to serve the remote sensing community in harvesting physical parameters from multi-platform, multi-sensor measurements that target different components of the atmosphere-oceanic system.

Investigation of the variations in the water leaving polarized reflectance from the POLDER satellite data over two biogeochemical contrasted oceanic areas

The biogeochemical characterization of marine particles suspended in sea water, is of fundamental importance in many areas of ocean science. Previous studies based on theoretical calculations and field measurements have demonstrated the importance of the use of the polarized light field in the retrieval of the suspended marine particles properties. However, because of the weakness of the water leaving polarized signal and of the limited number of appropriate spatial sensors, such measurements have never been exploited from space. Here we show that the marine polarized remote sensing reflectance, as detected from the POLarization and Directionality of the Earth's Reflectances (POLDER) sensor, can be measured from space over bright waters and in absence of aerosols. This feasibility study is carried out over two oceanic areas characterized by different nature of the bulk particulate assemblage: the Barents sea during an intense coccolithophore bloom, and the Rio de la Plata estuary waters dominated by suspended sediments. The retrieved absolute values of the degree of polarization, P, its angular pattern, and its behavior with the scattering level are consistent with theory and field measurements. Radiative transfer simulations confirm the sensitivity of the POLDER-2 P values to the nature of the particulate assemblage. These preliminary results are very promising for the assessment of the bulk particle composition from remote sensing of the polarized signal, at least over highly scattering waters.

Retrieval of chlorophyll fluorescence from reflectance spectra through polarization discrimination: modeling and experiments

The polarization discrimination technique we recently developed, shows that it is possible to separate the elastic scattering and the chlorophyll fluorescence signal from the water-leaving radiance by making use of the fact that the elastically scattered components are partially polarized, while the fluorescence signal is unpolarized. The technique has been shown to be applicable to a wide range of water conditions. We present an extension of experimental and analytical results, which serve to define the scope of this technique and its range of applicability. A new analysis, based on vector radiative transfer computations, and on laboratory and field measurements on eastern Long Island and in the Chesapeake Bay, shows that the technique is generally effective for both open ocean and coastal waters, but that it is limited if the ocean bottom albedo and/or multiple scattering due to very high mineral particle concentrations result in depolarizing the water-leaving radiance. In addition, we show that in contrast with the polarization-based retrieval, the traditional method of extracting fluorescence height using the baseline method can give significant errors, particularly for coastal waters where it strongly overestimates the fluorescence values.

Polarization-discrimination technique to maximize the lidar signal-to-noise ratio for daylight operations

The impact and potential of a polarization-selection technique to reduce the sky background signal for linearly polarized monostatic elastic backscatter lidar measurements are examined. Taking advantage of naturally occurring polarization properties in scattered skylight, we devised a polarization-discrimination technique in which both the lidar transmitter and the receiver track and minimize detected sky background noise while maintaining maximum lidar signal throughput. Lidar elastic backscatter measurements, carried out continuously during daylight hours at 532 nm, show as much as a factor of square root 10 improvement in the signal-to-noise ratio (SNR) over conventional unpolarized schemes. For vertically pointing lidars, the largest improvements are limited to the early morning and late afternoon hours, while for lidars scanning azimuthally and in elevation at angles other than vertical, significant improvements are achievable over more extended time periods with the specific times and improvement factors depending on the specific angle between the lidar and the solar axes. The resulting diurnal variations in SNR improvement sometimes show an asymmetry with the solar angle that analysis indicates can be attributed to changes in observed relative humidity that modifies the underlying aerosol microphysics and observed optical depth.