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

Time-dependent polarized radiative transfer in an atmosphere-ocean system exposed to external illumination

Time-dependent polarized radiative transfer in an atmosphere-ocean system exposed to external illumination is numerically investigated. The specular reflection and transmission effects based on the relative refractive index between the atmosphere and water are considered. A modified Monte Carlo (MMC) algorithm combined with time shift and superposition principle, which significantly improves the computational efficiency of the traditional Monte Carlo (TMC) method, is developed to simulate the time-dependent polarized radiative transfer process. The accuracy and computational superiority of the MMC for polarized radiative transfer in the atmosphere-ocean system are validated, and the time-resolved polarized radiative signals are discussed.

Radiative Transfer Modeling of Phytoplankton Fluorescence Quenching Processes

We report the first radiative transfer model that is able to simulate phytoplankton fluorescence with both photochemical and non-photochemical quenching included. The fluorescence source term in the inelastic radiative transfer equation is proportional to both the quantum yield and scalar irradiance at excitation wavelengths. The photochemical and nonphotochemical quenching processes change the quantum yield based on the photosynthetic active radiation. A sensitivity study was performed to demonstrate the dependence of the fluorescence signal on chlorophyll a concentration, aerosol optical depths and solar zenith angles. This work enables us to better model the phytoplankton fluorescence, which can be used in the design of new space-based sensors that can provide sufficient sensitivity to detect the phytoplankton fluorescence signal. It could also lead to more accurate remote sensing algorithms for the study of phytoplankton physiology.

Evaluation of Weighting Average Functions as a Simplification of the Radiative Transfer Simulation in Vertically Inhomogeneous Eutrophic Waters

Current water color remote sensing algorithms typically do not consider the vertical variations of phytoplankton. Ecolight with a radiative transfer program was used to simulate the underwater light field of vertical inhomogeneous waters based on the optical properties of a eutrophic lake (i.e., Lake Chaohu, China). Results showed that the vertical distribution of chlorophyll-a (Chla(z)) can considerably affect spectrum shape and magnitude of apparent optical properties (AOPs), including subsurface remote sensing reflectance in water (rrs(λ, z)) and the diffuse attenuation coefficient (Kx(λ, z)). The vertical variations of Chla(z) changed the spectrum shapes of rrs(λ, z) at the green and red wavelengths with a maximum value at approximately 590 nm, and changed the Kx(λ, z) from blue to red wavelength range with no obvious spectral variation. The difference between rrs(λ, z) at depth z m and its asymptotic value (Δrrs(λ, z)) could reach to ~78% in highly stratified waters. Diffuse attenuation coefficient of downwelling plane irradiance (Kd(λ, z)) had larger vertical variations, especially near water surface, in highly stratified waters. Three weighting average functions performed well in less stratified waters, and the weighting average function proposed by Zaneveld et al., (2005) performed best in highly stratified waters. The total contribution of the first three layers to rrs(λ, 0−) was approximately 90%, but the contribution of each layer in the water column to rrs(λ, 0−) varied with wavelength, vertical distribution of Chla(z) profiles, concentration of suspended particulate inorganic matter (SPIM), and colored dissolved organic matter (CDOM). A simple stratified remote sensing reflectance model considering the vertical distribution of phytoplankton was built based on the contribution of each layer to rrs(λ, 0−).

Ocean Color Analytical Model Explicitly Dependent on the Volume Scattering Function

An analytical radiative transfer (RT) model for remote sensing reflectance that includes the bidirectional reflectance distribution function (BRDF) is described. The model, called ZTT (Zaneveld-Twardowski-Tonizzo), is based on the restatement of the RT equation by Zaneveld (1995) in terms of light field shape factors. Besides remote sensing geometry considerations (solar zenith angle, viewing angle, and relative azimuth), the inputs are Inherent Optical Properties (IOPs) absorption a and backscattering bb coefficients, the shape of the particulate volume scattering function (VSF) in the backward direction, and the particulate backscattering ratio. Model performance (absolute error) is equivalent to full RT simulations for available high quality validation data sets, indicating almost all residual errors are inherent to the data sets themselves, i.e., from the measurements of IOPs and radiometry used as model input and in match up assessments, respectively. Best performance was observed when a constant backward phase function shape based on the findings of Sullivan and Twardowski (2009) was assumed in the model. Critically, using a constant phase function in the backward direction eliminates a key unknown, providing a path toward inversion to solve for a and bb. Performance degraded when using other phase function shapes. With available data sets, the model shows stronger performance than current state-of-the-art look-up table (LUT) based BRDF models used to normalize reflectance data, formulated on simpler first order RT approximations between rrs and bb/a or bb/(a + bb) (Morel et al., 2002; Lee et al., 2011). Stronger performance of ZTT relative to LUT-based models is attributed to using a more representative phase function shape, as well as the additional degrees of freedom achieved with several physically meaningful terms in the model. Since the model is fully described with analytical expressions, errors for terms can be individually assessed, and refinements can be readily made without carrying out the gamut of full RT computations required for LUT-based models. The ZTT model is invertible to solve for a and bb from remote sensing reflectance, and inversion approaches are being pursued in ongoing work. The focus here is with development and testing of the in-water forward model, but current ocean color remote sensing approaches to cope with an air-sea interface and atmospheric effects would appear to be transferable. In summary, this new analytical model shows good potential for future ocean color inversion with low bias, well-constrained uncertainties (including the VSF), and explicit terms that can be readily tuned. Emphasis is put on application to the future NASA Plankton, Aerosol, Cloud, and ocean Ecosystem (PACE) mission.

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.

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.