**Page updated:**
April 13, 2020 **Author:** Curtis Mobley

View PDF

# Atmospheric Transmittances

[Howard Gordon also contributed to this page.]

Equations such as (3) of the Problem Formulation page and (13) of the Normalized Reﬂectances page involve various direct ($T$) and diﬀuse ($t$) atmospheric transmission terms, which require discussion.

### 1 Direct Transmittance

Direct atmospheric transmission is used if only one particular path, or a narrow bundle of nearly colinear paths, connects the source and the observer. This is the case for specular reﬂection, as illustrated in the left panel of Fig. 1. When the sea surface is viewed by the sensor, only a small patch of sea surface is seen as the Sun’s specular reﬂection or direct glint. Other points of the sea surface would be seen as localized patches of Sun glint for other viewing directions, but not by the sensor viewing direction as shown. In each case, the reﬂected radiance is traveling in a very narrow set of directions determined by the Sun’s location and the law of reﬂection (${\mathit{\theta}}_{r}={\mathit{\theta}}_{s}$ for a level surface). The narrow beams of Sun glint for other viewing directions are represented by the dotted and dash lines in the left panel of the ﬁgure. These beams can inﬂuence the direction of interest only via two scatterings: the ﬁrst out of the reﬂected beam and the second into the direction of interest. Such a two-scatter path is shown by the light dashed line. Two scatterings makes the contributions of the unseen beams of specular reﬂection to the direction of interest very small.

Let $\tau $ be the atmospheric optical depth along a vertical path (the nadir viewing direction for the sensor). This $\tau $ includes all eﬀects of atmospheric absorption and attenuation by all atmospheric constituents. For an oﬀ-nadir viewing direction ${\mathit{\theta}}_{v}$, the direct transmittance is then simply

$$T\left({\mathit{\theta}}_{v}\right)=exp\left(-\tau \u2215cos{\mathit{\theta}}_{v}\right)\phantom{\rule{0.3em}{0ex}}.$$ | (1) |

This geometry is analogous to the Lambert-Beer law for radiance propagation of a beam through a homogeneous medium: $L\left(r\right)=L\left(0\right)exp\left(-cr\right)$, where $c$ is the beam attenuation coeﬃcient and $r$ is the distance traveled. In the present case, $\tau =c{r}_{a}$, where ${r}_{a}$ is the distance though the atmosphere on a vertical (nadir-viewing) path, and $r={r}_{a}\u2215cos{\mathit{\theta}}_{v}$ is the atmospheric path length along the viewing direction.

### 2 Diﬀuse Transmittance

For water-leaving radiance ${L}_{w}$, every point of the sea surface is emitting an upward distribution of radiance ${L}_{w}\left(\mathit{\theta},\varphi \right)$, as illustrated in the right panel of Fig. 1. Radiance from all locations and various directions can be scattered into the direction of interest via only one scattering, as illustrated by the dotted line in the right panel of Fig. 1. Radiance scattered out of the beam along the viewing direction can thus be replaced via a single scattering from the radiance emitted by a neighboring point propagating in a diﬀerent direction. The diﬀuse transmittance therefore depends not just on the atmospheric properties and viewing direction, but also on the angular distribution of ${L}_{w}$, which in general is unknown. This situation is analogous to the diﬀuse attenuation ${K}_{d}$ of downwelling irradiance within the water column of a plane parallel ocean that is illuminated at all points of the surface. Because radiance absorbed or scattered out of the path between the source and sensor can be replaced by radiance from other sea surface points and directions, the diﬀuse transmittance $t\left({\mathit{\theta}}_{v}\right)$ will be greater than the direct transmittance $T\left({\mathit{\theta}}_{v}\right)$ (just as the the diﬀuse transmittance $exp\left(-{K}_{d}z\right)$ is always greater than than the beam transmittance $exp\left(-cz\right)$). Accounting for the angular distribution of ${L}_{w}$ and for the scattering processes makes the computation of diﬀuse transmittance more complicated than for direct transmittance.

The diﬀuse transmittance of the water-leaving radiance along a particular viewing direction $\left({\mathit{\theta}}_{v},{\varphi}_{v}\right)$ is by deﬁnition

where ${L}_{w}$ is the water-leaving radiance at the sea surface and ${L}_{w}^{TOA}$ is the water-leaving radiance that reaches the TOA. Henceforth, all quantities are assumed to be functions of wavelength, and the argument $\lambda $ will be omitted. One way to compute $t\left({\mathit{\theta}}_{v},{\varphi}_{v}\right)$ is to run a coupled ocean-atmosphere radiative transfer model to compute the needed ${L}_{w}\left({\mathit{\theta}}_{v},{\varphi}_{v}\right)$ and ${L}_{w}^{TOA}\left({\mathit{\theta}}_{v},{\varphi}_{v}\right)$ for a wide range of atmospheric and oceanic conditions, Sun and viewing geometries, and wavelengths. The values of $t\left({\mathit{\theta}}_{v},{\varphi}_{v}\right)$ would then be obtained via Eq. (2) and tabulated for later use. Such tables would need to be constructed for all possible aerosol types (i.e, aerosol phase function shapes), aerosol optical thicknesses, water IOPs (or at least for chlorophyll concentrations in Case 1 waters), and viewing geometries. The tables would necessarily be large because of the large number of parameters that can aﬀect $t\left({\mathit{\theta}}_{v},{\varphi}_{v}\right)$.

Yang and Gordon (1997) examined the computation of diﬀuse transmittance and errors therein on the retrieved water-leaving radiance. Using a combination of radiative transfer numerical modeling of the ocean and atmosphere and reciprocity principles, they compared diﬀuse transmittances computed using realistic modeled and measured ${L}_{u}\left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)$ distributions versus diﬀuse transmittances computed on the assumption that ${L}_{u}$ is isotropic. Here ${L}_{u}\left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)$ denotes the upwelling underwater radiance distribution just beneath the sea surface. For a level sea surface, ${L}_{w}\left({\mathit{\theta}}_{v},{\varphi}_{v}\right)={T}_{F}\left({\mathit{\theta}}^{\prime}\right){L}_{u}\left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)\u2215{n}_{w}^{2}$, where ${n}_{w}$ is the index of refraction of the water, the in-water ${\mathit{\theta}}^{\prime}$ is related to the in-air ${\mathit{\theta}}_{v}$ by Snel’s law $sin{\mathit{\theta}}_{v}={n}_{w}sin{\mathit{\theta}}^{\prime}$, and ${T}_{F}\left({\mathit{\theta}}^{\prime}\right)$ is the Fresnel transmittance of the surface from water to air. When ${L}_{u}$ is isotropic, the diﬀuse transmittance is independent of the azimuthal angle and is denoted by ${t}^{\ast}\left({\mathit{\theta}}_{v}\right)$.

Yang and Gordon (1997) (Eq. 3) show via a clever use of reciprocity that the diﬀuse transmittance of radiance ${t}^{\ast}$ along an atmospheric path in the direction of the Sun at solar zenith angle ${\mathit{\theta}}_{o}$ is numerically equal to the diﬀuse transmittance of irradiance from the TOA to a depth just beneath the sea surface, on the assumption that there is no upwelling radiance within the water. That is,

where ${E}_{R}\left({\mathit{\theta}}_{o}\right)$ is the downwelling plane irradiance just beneath the sea surface for an extraterrestrial solar irradiance ${F}_{o}$ incident onto the TOA at angle ${\mathit{\theta}}_{o}$, and ${T}_{F}\left({\mathit{\theta}}_{o}\right)$ is the Fresnel downward transmission of the sea surface for radiance incident at angle ${\mathit{\theta}}_{o}$ from the normal. Since the upwelling radiance ${L}_{u}\left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)$ used to obtain this result is assumed to be isotropic, the azimuthal dependence of ${t}^{\ast}$ is irrevelant, and the desired diﬀuse attenuation for radiance at viewing direction ${\mathit{\theta}}_{v}={\mathit{\theta}}_{o}$ is equal to the value of the irradiance transmission at the same polar angle. The great virtue of Eq. (3) is that it allows the eﬃcient numerical computation of ${t}^{\ast}$ using backward (reverse) Monte Carlo simulation of downwelling irradiance for a given aerosol type and optical thickness.

The retrieved water-leaving radiance is, by Eq. (2),

Let ${L}_{w}^{\ast}$ denote the retrieval when ${t}^{\ast}$ rather than $t$ is used in Eq. (4). The error in the retrieved water-leaving radiance due to using ${t}^{\ast}$ rather than the exact $t$ is

$$\frac{\Delta {L}_{w}}{{L}_{w}}=\frac{{L}_{w}^{\ast}-{L}_{w}}{{L}_{w}}=\frac{t-{t}^{\ast}}{{t}^{\ast}}\phantom{\rule{0.3em}{0ex}}.$$ | (5) |

Yang and Gordon (1997) found that for viewing angles ${\varphi}_{v}$ perpendicular to the principle plane (the plane of the Sun), the errors in the retrieved ${L}_{w}$ are no more than 4% for viewing angles ${\mathit{\theta}}_{v}\le 45deg$, Sun zenith angles ${\mathit{\theta}}_{o}\le 60deg$, and aerosol optical thicknesses typical of clear atmospheres. The errors in band-ratio algorithms were less; e.g., the error in ${L}_{w}\left(443\right)\u2215{L}_{w}\left(555\right)$ (used to retrieve chlorophyll concentration) is less than 2% except for very clear water and some viewing directions, for which the error in the ratio is about 3%. These parameter ranges covered most of the needs for SeaWiFS and errors of this magnitude were deemed acceptable compared to other errors in the retrieval process (i.e., removal of atmospheric path radiance). In that case, the diﬀuse attenuation can be obtained by pre-computed functions of the form

$${t}^{\ast}\left({\mathit{\theta}}_{v}\right)=A\left({\mathit{\theta}}_{v}\right)exp\left[-B\left({\mathit{\theta}}_{v}\right){\tau}_{a}\right]\phantom{\rule{0.3em}{0ex}}.$$ | (6) |

where $A\left({\mathit{\theta}}_{v}\right)$ and $B\left({\mathit{\theta}}_{v}\right)$ are tabulated for each aerosol type. Thus, for $N$ angles ${\mathit{\theta}}_{v}$, only $2N$ numbers $A\left({\mathit{\theta}}_{v}\right)$ and $B\left({\mathit{\theta}}_{v}\right)$ need to be tabulated for each aerosol type. Moreover, the aerosol type and optical thickness are determined as part of the process to remove the aerosol contribution to the total path radiance. This enables selection of the appropriate $A,B$, and ${\tau}_{a}$, and evaluation of Eq. (4) is operationally feasible. If aerosols are ignored, Eq. (6) reduces to ${t}^{\ast}\left({\mathit{\theta}}_{v}\right)=exp\left[-\frac{1}{2}{\tau}_{R}\u2215cos{\mathit{\theta}}_{v}\right)]$, where ${\tau}_{R}$ is the Rayleigh optical thickness. This is the formula used in the early days for CZCS atmospheric correction. The value of Eq. (6) is that it allows the aerosol optical thickness ${\tau}_{a}$ to be incorporated into the diﬀuse transmittance calculations via a simple exponential and pre-computed $A$ and $B$ values.

However, for larger oﬀ-nadir viewing angles ${\mathit{\theta}}_{v}$, for azimuthal viewing directions near ${\varphi}_{v}=0$ or 180 deg, and for very clear water, the errors $\Delta {L}_{w}\u2215{L}_{w}$ can be as much as 6%. This could be signiﬁcant for the MODIS Aqua sensor, which views a wide range of ${\varphi}_{v}$ directions. Gordon and Franz (2008) therefore re-examined the model of Yang and Gordon (1997) and developed a correction term to ${t}^{\ast}$ so that Eq. (4) becomes

The $\delta \left({\mathit{\theta}}_{v},{\varphi}_{v}\right)$ factor corrects for the bi-directional eﬀects resulting from the use of an isotropic ${L}_{u}$ in the computation of ${t}^{\ast}$, rather than the exact, non-isotropic ${L}_{u}\left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)$. The ${L}_{u}\left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)$ needed for computation of $\delta $ was obtained from ${L}_{u}=\left(f\u2215Q\right)\left({b}_{b}\u2215a\right)$ and chlorophyll-based models for $f\u2215Q$ and ${b}_{b}\u2215a$. Again, the $\delta $ correction term can be pre-computed and tabulated for various aerosol types and water IOPs. Evaluation of the impact of the $\delta $ correction applied to both SeaWiFS and MODIS Aqua data showed that retrieved water-leaving radiances will be in error by no more than $\sim \phantom{\rule{0.3em}{0ex}}1\%$ if ${\mathit{\theta}}_{v}<60deg$ and the $\delta $ correction is omitted. For ${\mathit{\theta}}_{v}>60deg$, i.e. near the edges of scan lines, the use of the $\delta $ correction is warranted.

In current operational practice, the tabulated $A$ and $B$ functions are applied for both downwelling (solar irradiance) and upwelling ${L}_{w}$ paths. The $\delta \left({\mathit{\theta}}_{v}\right)$ correction of Gordon and Franz (2008) is not applied because it is an added complication with no signiﬁcant impact in most instances.