It looks like you're using Internet Explorer 6: Features on this site are not supported by that browser version. Please upgrade to the latest version of Internet Explorer.

You are viewing a version of the Ocean Optics Web Book site specially-formatted for printing.

Ocean Optics Web Book

Skip to main content
Ocean Optics Web Book

Atmospheric Correction

Normalized Reflectances

Page updated: May 25, 2016
Principal author: Curtis Mobley
 

[Howard Gordon and David Antoine also contributed to this page.]

Ocean-color remote sensing algorithms usually work with remote-sensing reflectances or normalized water-leaving reflectances. The calculation and interpretation of those quantities are discussed in detail on this page.

The ratio of water-leaving radiance $ L_{\rm {w}}$ to incident sky irradiance $ E_{\rm {d}}$ is an apparent optical property (AOP) that has only weak dependence on external parameters such as solar zenith angle and sky conditions, but which is strongly correlated with water-column inherent optical properties (IOPs). However, the remote-sensing reflectance $ R_{\rm {rs}} \equiv L_{\rm {w}}/E_{\rm {d}}$ still does depend somewhat on the atmospheric and other conditions at the time of measurement and thus, strictly speaking, is tied to the particular time and location of the observation. (Further discussion and examples are given on the Reflectances page.)

It would be desirable to have an AOP that completely removes the effects of solar zenith angle, viewing direction, atmospheric conditions, and sea state, while retaining a strong dependence on the water IOPs. It would then be possible to compare this AOP for measurements made at different times and/or locations, and thereby extract information about the differences in the water columns for the different measurements. Even for measurements made at the same time and location, normalization to a common set of conditions is needed, e.g., when comparing in situ measurements having different viewing directions. Such an AOP is obtained via the concept of the normalized water-leaving reflectance.

Normalized Radiances and Reflectances

Let $ L_{\rm {u}}(z, \theta_{\rm {s}}, \theta_{\rm {v}}, \phi)$ be the in-water, upwelling radiance at depth $ z$ , for a Sun zenith angle of $ \theta_s$ and a viewing direction of $ \theta_{\rm {v}}, \phi$ . Polar viewing direction $ \theta_{\rm {v}} = 0$ indicates a direction looking at the nadir, detecting radiance traveling toward the zenith. The azimuthal angle $ \phi$ is measured relative to the Sun's azimuthal direction. $ L_{\rm {w}}(\theta_{\rm {s}}, \theta_{\rm {v}}, \phi)$ denotes the corresponding water-leaving radiance, which is measured just above the sea surface. These radiances depend strongly on wavelength, which is not shown. In practice, $ L_{\rm {u}}$ is measured by instruments in the water. $ L_{\rm {w}}$ can be obtained by atmospheric correction of a radiance measured at the top of the atmosphere, from an above-surface measurement at sea level after correction for surface reflectance, or from an in-water measurement of $ L_{\rm {u}}$ extrapolated upward through the sea surface.

One goal of normalization is to transform a satellite-based measurement of top-of-the-atmosphere radiance $ L_{\rm {t}}$ into something that can be compared with a standard measurement made in situ, in the ocean, for whatever Sun zenith angle, viewing direction, atmospheric conditions, and wave state occurred at the time of the satellite measurement. Let this standard in situ measurement be the nadir-viewing radiance measured just below the sea surface, $ L_{\rm {u}}(0^-, \theta_{\rm {v}} = 0)$ . Depth $ z = 0^-$ refers to a location in the water just below the sea surface; $ 0^+$ refers to a location in the air just above the sea surface. Dividing $ L_{\rm {u}}(0^-, \theta_{\rm {v}} = 0)$ by the downwelling plane irradiance within the water, $ E_{\rm {d}}(0^-)$ , gives the in-water Remote Sensing Ratio $ RSR$ :

$\displaystyle RSR \equiv \frac{L_{\rm {u}}(0^-, \theta_{\rm {v}} = 0)}{E_{\rm {d}}(0^-)} \,.$ (1)

The division of $ L_{\rm {u}}(0^-, \theta_{\rm {v}} = 0)$ by $ E_{\rm {d}}(0^-)$ removes the "zeroth order" effect of solar zenith angle $ \theta_{\rm {s}}$ and the "first order" atmospheric effects (including aerosol effects) on the magnitude of $ L_{\rm {u}}(0^-, \theta_{\rm {v}} = 0)$ . We now want to transform $ L_{\rm {t}}$ into something comparable.

In the early days of satellite remote sensing, it was sometimes assumed that the upwelling underwater radiance distribution is isotropic. Under that assumption, $ RSR$ is approximately what you would get if the Sun were at the zenith and there were no atmosphere (i.e., the sky were black). This was the origin of statements that the normalized water-leaving radiance is the water-leaving radiance "which would exit the sea surface if the Sun were at the zenith and if the atmosphere were absent" (Gordon et al. (1988), page 10,910).

Fig. [*] compares upwelling and water-leaving radiances for "no atmosphere" or "black sky" vs. realistic sky conditions. The curves of this figure were generated using HydroLight with the Sun placed at the zenith. The sky was either black (a collimated incident sky radiance) or had a diffuse radiance angular distribution typical of a clear sky. The water IOPs were determined using the same "new Case 1" bio-optical model as for the simulations of the previous section. Runs were made for chlorophyll values of 0.05, 0.5 and 5 $ \unit{mg~m^{-3}}$ . The runs were at a wavelength of 550 nm and the sea surface was level. Each radiance is normalized by its value at the nadir-viewing direction to isolate the differences in the shapes of the curves. The upper set of curves shows the shape of the upwelling radiances $ L_{\rm {u}}$ just below the sea surface as functions of the in-water, off-nadir viewing angle $ \theta_{\rm {v}}'$ at right angles ( $ \phi = 90~\unit{deg}$ ) to the solar plane. The lower set of curves (those curving downward in the figure) shows the corresponding water-leaving radiances $ L_{\rm {w}}$ as functions of the in-air, off-nadir viewing angle $ \theta_{\rm {v}}$ . For a level sea surface, these in-water and in-air angles are related by Snel's law $ \sin \theta_{\rm {v}} = n_w \sin \theta_{\rm {v}}'$ , where $ n_w \approx 1.34$ is the water index of refraction. Although the differences in the black-sky and real-sky radiances are only a few percent over the range of angles relevant to most remote sensing ( $ \theta_{\rm {v}} \lesssim 60$ deg), differences of this magnitude are significant given the high accuracy requirements for retrieved water-leaving radiances in ocean remote sensing.

Figure: Comparisons of $ L_{\rm{u}}(z=0^-, \theta_{\rm{s}}=0, \theta_{\rm{v}}', \phi=90)$ and $ L_{\rm{w}}(\theta_{\rm{s}}=0, \theta_{\rm{v}}, \phi=90)$ for a zenith Sun in a black sky (no atmosphere; dashed lines) vs. a zenith Sun in a typical atmosphere (solid lines). The colors identify the chlorophyll concentrations. The black line at an ordinate value of 1 corresponds to an isotropic radiance distribution.
Image a03fe7417dc3b33aac6de78722781e94

The observation that the upwelling radiance distribution depends on the sky radiance distribution indicates that the idea of removing the atmosphere is too extreme.Morel and Gentili (1996) (page 4852) therefore revised the definition of normalized water-leaving radiance to be "...the radiance that could be measured by a nadir-viewing instrument, if the Sun were at the zenith in the absence of any atmospheric loss, and when the Earth is at its mean distance from the Sun." The distinction between "atmosphere were absent" (i.e., a vacuum) and "absence of any atmospheric loss" (i.e., no attenuation by the atmosphere) is fundamental to understanding the developments of this section.

The normalization proceeds as follows. The first step is to account for solar zenith angle and atmospheric attenuation effects on $ L_{\rm {w}}(\theta_{\rm {s}}, \theta_{\rm {v}}, \phi)$ via (e.g., Gordon and Clark (1981), page 10,910; or Gordon and Wang (1994b), Eq. 4):

$\displaystyle [L_{\rm {w}}(\theta_{\rm {v}}, \phi)]_{\rm {N}} \equiv \frac{L_{\rm {w}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)}{\cos \theta_{\rm {s}} \,t(\theta_{\rm {s}})} \,,$    

where $ t(\theta_{\rm {s}})$ is the atmospheric diffuse transmittance for irradiance in the Sun's direction for the given atmospheric conditions. Recent papers include an explicit factor to correct $ L_{\rm {w}}$ for the Earth-Sun distance at the time of measurement:

$\displaystyle [L_{\rm {w}}(\theta_{\rm {v}}, \phi)]_N \equiv \left ( \frac{R}{R_{\rm {o}}} \right )^2 \frac{ L_{\rm {w}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)}{\cos \theta_{\rm {s}} \,t(\theta_{\rm {s}})} \,.$ (2)

Here $ R$ is the Earth-Sun distance at the time of measurement, and $ R_{\rm {o}}$ is the mean Earth-Sun distance. The $ (R/R_{\rm {o}})^2$ factor corrects the $ L_{\rm {w}}$ measurement to what it would be at the mean Earth-Sun distance. (The solar irradiance at the TOA varies by about 8% over the course of a year due to the Earth's elliptical orbit.)

$ [L_{\rm {w}} (\theta_{\rm {v}}, \phi)]_{\rm {N}}$ is called the normalized water-leaving radiance. It is the water-leaving radiance that would occur if the earth were at the mean Earth-Sun distance, the Sun were at the zenith, and the atmosphere were non-attenuating. Note that although the factors of $ (R/R_{\rm {o}})^2$ , $ \cos \theta_{\rm {s}}$ , and $ t(\theta_{\rm {s}})$ largely remove the effects of Earth-Sun distance, solar zenith angle, and atmospheric attenuation, respectively, on the measured $ L_{\rm {w}}$ , the normalized water-leaving radiance still refers to a particular viewing direction and depends on the sky angular radiance distribution at the time of observation.

Multiplying the $ [L_{\rm {w}} (\theta_{\rm {v}}, \phi)]_{\rm {N}}$ of Eq. ([*]) by a factor of $ \pi/F_{\rm {o}}$ , where $ \pi$ has units of steradian and $ F_{\rm {o}}$ is the extraterrestrial solar irradiance at the mean Earth-Sun distance, gives the nondimensional normalized water-leaving reflectance $ [\rho_{\rm {w}}]_{\rm {N}}$ (Gordon and Wang (1994b), page 7756):

$\displaystyle [\rho_{\rm {w}}(\theta_{\rm {v}}, \phi)]_{\rm {N}} \equiv \frac{\pi}{F_{\rm {o}}} [L_{\rm {w}}(\theta_{\rm {v}}, \phi)]_{\rm {N}} = \pi \frac{ \left ( \frac{R}{R_{\rm {o}}} \right )^2 L_{\rm {w}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)}{F_{\rm {o}} \, \cos \theta_{\rm {s}} \,t(\theta_{\rm {s}})} \,.$ (3)

(Radiant energy is a physical quantity that propagates through space and that can leave the water. Radiance is a derived physical quantity that likewise can leave the water, so it makes sense to speak of the "water-leaving radiance." Reflectance, however, is a property of a surface. Reflectance does not propagate through space and it cannot leave a surface, so it does not make sense to speak of the "water-leaving reflectance." However, the term "water-leaving reflectance" is well established shorthand for "reflectance based on the water-leaving radiance and the incident irradiance," and we will use the term even if it is linguistically somewhat improper.)

The remote-sensing reflectance $ R_{\rm {rs}}$ is usually defined as

$\displaystyle R_{\rm {rs}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi) \equiv \frac{L_{\rm {w}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)}{E_{\rm {d}}(0^+, \theta_{\rm {s}})} \,.$ (4)

In this definition, both $ L_{\rm {w}}$ and $ E_{\rm {d}}$ are values for the Earth-Sun distance at the time of measurement. However, this $ R_{\rm {rs}}$ is numerically the same as what would be obtained if both $ L_{\rm {w}}$ and $ E_{\rm {d}}$ were corrected to the mean Earth-Sun distance by $ (R/R_{\rm {o}})^2$ factors applied to each, because the correction factors on $ L_{\rm {w}}$ and $ E_{\rm {d}}$ cancel out. That is to say, the irradiance at the sea surface for Earth-Sun distance $ R$ is

$\displaystyle E_{\rm {d}}(0^+, \theta_{\rm {s}}) = F_{\rm {o}} \left ( \frac{R_{\rm {o}}}{R} \right )^2 \cos \theta_{\rm {s}} \, t(\theta_{\rm {s}}) \,.$ (5)

It thus follows that

$\displaystyle [\rho_{\rm {w}}(\theta_{\rm {v}}, \phi)]_{\rm {N}} = \pi \frac{ \left ( \frac{R}{R_{\rm {o}}} \right )^2 L_{\rm {w}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)}{F_{\rm {o}} \cos \theta_{\rm {s}} \, t(\theta_{\rm {s}})} = \pi \frac{L_{\rm {w}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)}{E_{\rm {d}}(\theta_{\rm {s}})} = \pi R_{\rm {rs}}(\theta_{\rm {v}}, \phi) \,.$ (6)

Another way to view $ [\rho_{\rm {w}}]_{\rm {N}}$ is to think of it as the bidirectional reflectance distribution function (BRDF) of the ocean normalized by the BRDF of a perfectly reflecting Lambertian surface. The BRDF of a surface as measured in the laboratory is the radiance reflected by the surface divided by the incident plane irradiance onto the surface. The BRDF of a Lambertian surface whose irradiance reflectance is $ \mathcal{R}$ is $ \mathcal{R}/\pi$ , with units of inverse steradian. For a perfect Lambertian reflector, $ \mathcal{R} = 1$ , and

$\displaystyle [\rho_{\rm {w}}]_{\rm {N}} = \frac{BRDF_{\rm {ocean}}}{BRDF_{\rm {Lamb}}} = \frac{L_{\rm {w}} / E_{\rm {d}}}{1/\pi} = \pi R_{\rm {rs}}\,.$    

This makes it clear that the $ \pi$ carries units of solid angle, so that $ [\rho_{\rm {w}}]_{\rm {N}}$ is nondimensional.


The BRDF Effect

As noted above, the normalizations contained in $ [\rho_{\rm {w}}(\theta_{\rm {v}}, \phi)]_{\rm {N}}$ or $ R_{\rm {rs}}(\theta_{\rm {v}}, \phi)$ remove the effects of solar zenith angle, atmospheric attenuation, and Earth-Sun distance on a measured radiance $ L_{\rm {w}}$ . However, $ [\rho_{\rm {w}}(\theta_{\rm {v}}, \phi)]_{\rm {N}}$ still refers to a particular viewing direction $ (\theta_{\rm {v}}, \phi)$ . This dependence ties $ [\rho_{\rm {w}}(\theta_{\rm {v}}, \phi)]_{\rm {N}}$ to the angular distribution of the upwelling underwater radiance and to the transmittance through the sea surface from water to air (which depends on the wave state, i.e., on the wind speed). The upwelling underwater radiance in turn depends on the angular distribution of the incident sky radiance, surface transmittance from air to water, and to the absorbing and scattering properties of the water body (the scattering phase function in particular). The dependence of the upwelling radiance distribution on the sky radiance distribution, viewing geometry, and water optical properties is commonly called the BRDF effect. The final step is to remove the BRDF effect to the greatest extent possible.

The BRDF effect was studied by Morel and colleagues in a series of papers, Morel and Gentili (1991), Morel and Gentili (1993), and Morel and Gentili (1996), culminating in Morel et al. (2002). They used numerical radiative transfer models to compute correction factors that would transform a measurement made for a particular Sun zenith angle, viewing direction, wind speed, atmospheric conditions, and water IOPs into a measurement that corresponds to a zenith Sun and nadir viewing direction for a typical marine atmosphere and for Case 1 water with a given chlorophyll value. The correction involves three separate factors, $ \mathfrak{R}$ , $ f$ , and $ Q$ , as follows.

Let $ \mathfrak{R}(\theta_{\rm {v}}',W)$ be a nondimensional factor that accounts for all transmission and reflection effects by the wind-blown sea surface when $ E_{\rm {d}}(0^+)$ is transmitted downward through the surface to give $ E_{\rm {d}}(0^-)$ , and $ L_{\rm {u}}(0^-,\theta_{\rm {v}}',\phi)$ is transmitted upward through the surface to give $ L_{\rm {w}}(0^+,\theta_{\rm {v}},\phi)$ . Polar angle $ \theta_{\rm {v}}'$ (measured from the nadir) is the underwater angle that is refracted by the surface into the above-surface viewing direction $ \theta_{\rm {v}}$ of the water-leaving radiance $ L_{\rm {w}}(\theta_{\rm {v}},\phi_{\rm {v}})$ . $ W$ is the wind speed. $ \mathfrak{R}(\theta_{\rm {v}}',W)$ depends on the wind speed (i.e., the surface wave state) and the water index of refraction via the Snel's law mapping of $ \theta_{\rm {v}}'$ to $ \theta_{\rm {v}}$ . However, Gordon (2005) showed that the dependence of $ \mathfrak{R}(\theta_{\rm {v}}',W)$ on wind speed is very weak, and usually $ \mathfrak{R}$ can be computed with adequate accuracy over a wide range of viewing angles using $ W = 0$ . The detailed derivation of $ \mathfrak{R}(\theta_{\rm {v}}',W)$ in Morel and Gentili (1996) (Eq. 5 and Appendix D) shows that, strictly speaking, $ \mathfrak{R}$ also depends on the solar zenith angle and the angular distribution of the incident atmospheric radiance, which affect how much incident irradiance is transmitted through the sea surface. Likewise, $ \mathfrak{R}$ depends on the water IOPs via the in-water irradiance reflectance $ R(0^+) = E_{\rm {u}}(0^+)/E_{\rm {d}}(0^+)$ . However, these dependencies are weak compared to the directional ( $ \theta_{\rm {v}}'$ ) dependency and so, for compactness and consistency with Morel's notation, are not shown as arguments.

Let $ \mathfrak{R}_{\rm {o}}(W)$ be the reference value of $ \mathfrak{R}(\theta_{\rm {v}}',W)$ corresponding to transmission normal to the mean sea surface ( $ \theta_{\rm {v}} = \theta_{\rm {v}}' = 0$ ). Multiplying the $ [L_{\rm {w}}(\theta_{\rm {v}})]_{\rm {N}}$ of Eq. ([*]) by $ \mathfrak{R}_{\rm {o}}(W)/\mathfrak{R}(\theta_{\rm {v}}',W)$ corrects for surface-transmission effects for the actual viewing direction $ \theta_{\rm {v}}'$ and wind speed $ W$ . For $ \theta_{\rm {v}}' \approx 0$ , $ \mathfrak{R}_{\rm {o}} \approx 0.53$ ; $ \mathfrak{R}_{\rm {o}}$ is greater than 0.52 out to angles of 50 deg. Note that Figure 4 of Morel et al. (2002), which shows ratios of $ \mathfrak{R}_{\rm {o}}(W)/\mathfrak{R}(\theta_{\rm {v}}',W)$ , was incorrectly computed; see the discussion and revised figures in Gordon (2005).

The dimensionless factor $ f$ is defined by $ E_{\rm {u}}(0^-)/E_{\rm {d}}(0^-) \equiv f \times (b_{\rm {b}}/a)$ , where $ a$ and $ b_{\rm {b}}$ are the water absorption and backscatter coefficients respectively. This factor parameterizes how downwelling irradiance in the water is converted to upwelling irradiance by backscatter and reduced by absorption. That is to say, $ f$ relates the irradiance reflectance within the water to the most relevant IOPs. $ f$ values are in the range of 0.3 to 0.5 (Morel and Gentili (1996), Fig. 2).

The factor $ Q \equiv E_{\rm {u}}(0^-)/L_{\rm {u}}(0^-)$ (units of steradian) describes the angular distribution of the upwelling radiance. $ Q = \pi~\unit{sr}$ for an isotropic upwelling radiance distribution; actual in-water radiance distributions typically have $ Q$ values in the range of 3 to 6 sr (Morel and Gentili (1996), Fig. 3).

In practice, $ f$ and $ Q$ are combined to give a term that has less variability than the individual factors. The combined factor $ f/Q \equiv L_{\rm {u}}(0^-,\theta_{\rm {v}}',\phi)/E_{\rm {d}}(0^-)$ describes how the downwelling irradiance just beneath the sea surface is reflected back upward as upwelling radiance in the direction ( $ \theta_{\rm {v}}',\phi$ ). The $ f/Q$ term thus describes both the efficiency of conversion of downwelling irradiance into upwelling radiance, and the angular (non-isotropic) distribution of the upwelling underwater radiance that generates the water-leaving radiance. $ f/Q$ values are typically in the range to 0.07 to 0.15 (Morel and Gentili (1996), Fig. 6). Let $ f_{\rm {o}}/Q_{\rm {o}}$ refer to the ratio for the nadir viewing direction and Sun at the zenith.

Multiplying $ [L_{\rm {w}} (\theta_{\rm {v}}, \phi)]_{\rm {N}}$ by $ (f_{\rm {o}}/Q_{\rm {o}})/(f/Q)$ corrects for the difference of the actual angular distribution of the upwelling radiance and what that distribution would be for the Sun at the zenith, nadir viewing, for the particular atmospheric and oceanic conditions used to compute $ f$ and $ Q$.

Applying these BRDF corrections to the $ [L_{\rm {w}} (\theta_{\rm {v}}, \phi)]_{\rm {N}}$ of Eq. ([*]) gives (Morel et al. (2002), Eq. 13):

$\displaystyle [L_{\rm {w}}]_{\rm {N}}^{\rm {ex}} \equiv [L_{\rm {w}}(\theta_{\rm {v}},\phi)]_{\rm {N}} \, \frac{\mathfrak{R}_{\rm {o}}(W)}{\mathfrak{R}(\theta_{\rm {v}}', W)} \,\frac{f_{\rm {o}}({\rm {ATM}},W,\rm {IOP})}{Q_{\rm {o}}({\rm {ATM}},W,\rm {IOP})} \, \left [ \frac{f(\theta_{\rm {s}},{\rm {ATM}},W,\rm {IOP})}{Q(\theta_{\rm {s}},\theta_{\rm {v}}',\phi,{\rm {ATM}},W,\rm {IOP})} \right ]^{-1} \,.$ (7)

Morel et al. call $ [L_{\rm {w}}(\lambda)]_{\rm {N}}^{\rm {ex}}$ the "exact normalized water-leaving radiance." The arguments "ATM" and "IOP" refer to the specific set of atmospheric conditions and water inherent optical properties used to compute $ f$ and $ Q$ . (As previously noted, these arguments are omitted from the $ \mathfrak{R}$ terms because the ATM and IOP dependencies are small for $ \mathfrak{R}$ .) This $ [L_{\rm {w}}(\lambda)]_{\rm {N}}^{\rm {ex}}$ is equivalent to the normalization seen in Eq. (2) of Franz et al. (2007), except for a sensor-specific correction factor for out-of-band wavelength response. It should be noted that the $ f_{\rm {o}}/Q_{\rm {o}}$ factor has arguments of ( $ {\rm {ATM}},W,{\rm {IOP}}$ ) because those values correspond to the same atmospheric and oceanic conditions as the $ f/Q$ factor; the difference is that $ f_{\rm {o}}/Q_{\rm {o}}$ corresponds to $ \theta_{\rm {s}} = 0$ and $ \theta_{\rm {v}}' = 0$ . Values of the $ (\mathfrak{R}_o/\mathfrak{R})(f_{\rm {o}}/Q_{\rm {o}})/(f/Q)$ product are typically in the 0.6 to 1.2 range, depending on the IOPs, solar zenith angle, atmospheric conditions, wind speed, and wavelength.

It is noted that the wind speed at the time of observation, $ W$ , is shown in both the $ \mathfrak{R}_{\rm {o}}$ and $ (f_{\rm {o}}/Q_{\rm {o}})$ terms in Eq. ([*]). It can be argued that these terms should use a reference value of $ W = 0$ , so that all quantities are referred to a level sea surface. However, as presently implemented by the OBPG, the same wind speed $ W$ is used in the reference terms.

The atmospheric conditions used to compute the $ f/Q$ factors were typical of marine atmospheres and are summarized as follows (Morel and Gentili (1996); Morel et al. (2002)):

  • The atmosphere was modeled as 50 layers, each 1 km thick
  • The aerosol optical thickness was $ \tau_a = 0.2$ at 550 nm.
  • The tropospheric aerosols had a relative humidity of 70% and were distributed in the upper 45 atmospheric layers
  • The marine aerosols had a relative humidity of 90% and were distributed in the lower 5 atmospheric layers
  • The aerosol phase functions were modeled using Mie theory and the data of Shettle and Fenn (1979).
  • The sea surface was modeled using a Gaussian distribution of wave slopes based on the empirical wind-speed, wave-slope data of Cox and Munk (1954) (page 847) for the given wind speed $ W$ . $ W$ was taken to be 0 for the $ Q$ calculations, although there are still some residual capillary waves according to the empirical slope data of Cox and Munk.

The water optical properties were obtained from bio-optical models for Case 1 water, for which the IOPs can all be parameterized by the chlorophyll concentration $ Chl$ . In particular, the scattering phase function was computed as a chlorophyll-weighted sum of phase functions for "small" and "large" particles, which themselves were computed using T-matrix theory for non-spherical particles with different size distributions. The details of the IOP bio-optical models are given in Morel et al. (2002).

The radiative transfer calculations were carried out using a Monte Carlo code for calculation of the $ \mathfrak{R}$ factors, and using HydroLight for the $ f/Q$ calculations. The $ \mathfrak{R}$ are tabulated for an exactly level sea surface and for wind speeds of $ W = 0, 4, 10$ and $ 16~\unit{m~s^{-1}}$ ; the $ W = 0$ values include a small amount of residual capillary waves because the Cox-Munk mean square sea surface slopes are not extactly zero for a zero wind speed. The $ f/Q$ calculations were done both with and without Raman scattering by water. None of the calculations included polarization. These codes were run for

  • 7 wavelengths (412.5, 442.5, 490, 510, 560, 620, and 660 nm)
  • 6 chlorophyll values ($ Chl$ = 0.03, 0.1, 0.3, 1.0, 3.0 and $ 10~\unit{mg~m^{-3}}$ ). The water was homogeneous.
  • 6 solar zenith angles ( $ \theta_{\rm {s}}$ = 0, 15, 30, 45, 60 and 75 deg)
  • 13 azimuthal angles ($ \phi$ = 0 to 180 deg by steps of 15 deg)
  • 17 nadir angles ( $ \theta_{\rm {v}}'$ = 1.078, 3.411, 6.289, 9.278, 12.300, 15.330, 18.370, 21.410, 24.450, 27.500, 30.540, 33.590, 36.640, 39.690, 42.730, 45.780, and 48.830 deg)
These runs give a total of $ 7 \times 6 \times 6 \times 13 \times 17 = 55,629 ~f/Q$ values, which are organized into $ 7 \times 6 \times 6 = 252$ tables, each with 13 columns and 17 rows. A separate table gives values of $ \mathfrak{R}(\theta_{\rm {v}}',W)$ at $ \theta_{\rm {v}}'$ increments of 1 deg and for wind speeds of $ W = 0, 4, 10$ and $ 16~\unit{m~s^{-2}}$ . The tables including Raman effects are available at BRDF tables.

Use of these tables requires the chlorophyll concentration, which is not a priori known. In the initial study of Morel and Gentili (1996), a band-ratio algorithm was used to obtain an initial guess for $ Chl$ . That value was then used in the tables to obtain the BRDF correction factors. Those factors give a new estimate of $ [L_{\rm {w}}(\lambda)]_{\rm {N}}^{\rm {ex}}$ , which can be used to obtain an new value of $ Chl$ , and so on. However, the tables are not tied to how the $ Chl$ value is obtained, so in practice any algorithm can be used to obtain $ Chl$ values from the reflectances in an iterative process.

Finally, the exact normalized water-leaving radiance of Eq. ([*]) is used in Eq. ([*]) to obtain the exact normalized water-leaving reflectance:

$\displaystyle [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}} \equiv$ $\displaystyle ~ \frac{\pi}{F_{\rm {o}}} [L_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ (8)
$\displaystyle =$ $\displaystyle ~\left\{ \frac{\pi}{F_{\rm {o}} \cos \theta_{\rm {s}} \,t(\theta_{\rm {s}})} \frac{\mathfrak{R}_{\rm {o}}(W)}{\mathfrak{R}(\theta_{\rm {v}}',W)} \, \frac{f_{\rm {o}}({\rm {ATM}},W,\rm {IOP})}{Q_{\rm {o}}({\rm {ATM}},W,\rm {IOP})} \, \left [ \frac{f(\theta_{\rm {s}},{\rm {ATM}},W,\rm {IOP})}{Q(\theta_{\rm {s}},\theta_{\rm {v}}',\phi,{\rm {ATM}},W,\rm {IOP})} \right ]^{-1} \right \} \times$    
  $\displaystyle ~~~L_{\rm {w}}(\theta_{\rm {s}}, \theta_{\rm {v}}, \phi)\,.$ (9)

The quantities in brackets can all be obtained from pre-computed look-up-tables given the Sun and viewing geometry, wavelength, atmospheric conditions used to obtain $ L_{\rm {w}}(\theta_{\rm {s}}, \theta_{\rm {v}}, \phi)$ from $ L_{\rm {t}}(\theta_{\rm {s}}, \theta_{\rm {v}}, \phi)$ (used to determine $ t(\theta_{\rm {s}})$ ), and the chlorophyll concentration. The chlorophyll concentration determines the IOPs according to the bio-optical models for Case 1 water used in the Morel et al. calculations.

Although Morel and others call $ [L_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ the "exact" normalized water-leaving radiance, and $ [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ the "exact" normalized water-leaving reflectance, it must be remembered that these quantities are exact only if the atmosphere and ocean have exactly the same absorbing and scattering properties as used in the model simulations upon which the BRDF correction factors are based. That will of course not in general be the case. In Case 1 waters, the differences between the Morel et al. Case 1 IOP model and the actual ocean IOPs are often small enough that the BRDF-corrected quantities are sufficiently accurate for remote sensing. However, the differences can become large in Case 2 waters. Research therefore continues on ways to improve the BRDF correction, both to extend its validity to Case 2 water and to remove the need to estimate the chlorophyll concentration in order to use the look-up tables (Lee et al. (2011), Fan et al. (2016)). Although further improvements can be anticipated, the BRDF correction as described above remains the OPBG operational algorithm at the time of this writing.

The OBPG works with radiance to get to $ [L_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ . However, when doing atmospheric correction on TOA radiances, the various look-up-tables used for Rayleigh correction, etc. are in terms of reflectances $ [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ . Equation ([*]) allows easy conversion from one to the other, depending on which quantity is most convenient for a given step of the atmospheric correction process.

It is to be noted that the "remote-sensing reflectance" reported by the NASA OBPG as a standard Level 2 product of the NASA ocean color satellites such as MODIS and VIIRS is the exact normalized water-leaving reflectance of Eq. ([*]), divided by $ \pi$ :

$\displaystyle R_{\rm {rs}}({\rm {NASA}}) = \frac{[\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}}}{\pi} = \frac{[L_{\rm {w}}]_{\rm {N}}^{\rm {ex}}}{F_{\rm {o}}} \,.$ (10)

However, the "remote-sensing reflectance" computed by HydroLight is $ R_{\rm {rs}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)$ as defined by Eq. ([*]). A given HydroLight run computes $ R_{\rm {rs}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)$ for all viewing directions $ \theta_{\rm {v}}, \phi$ for the given $ \theta_{\rm {s}}$ and other conditions of wind speed, IOPs, and atmospheric radiance distribution. Thus the HydroLight $ R_{\rm {rs}}$ outputs for various viewing directions incorporate the BRDF effects. (Indeed, comparison of HydroLight-computed $ R_{\rm {rs}}(\theta_{\rm {s}},\theta_{\rm {v}}, \phi)$ values with HydroLight values for a zenith Sun and nadir viewing direction is how the $ f/Q$ BRDF factors were determined in the Morel et al. studies.) Thus there is no need for an explicit BRDF correction to the HydroLight-computed $ R_{\rm {rs}}$ . If a HydroLight run has the Sun at the zenith, then the nadir-viewing HydroLight $ R_{\rm {rs}}$ times $ \pi$ corresponds to $ [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$:

$\displaystyle [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}} = \pi R_{\rm {rs}}({\rm {HydroLight}};\theta_{\rm {s}} = 0,\theta_{\rm {v}} = 0) \,.$ (11)

Recalling Eq. ([*]), this gives

$\displaystyle R_{\rm {rs}}({\rm {NASA}}) = R_{\rm {rs}}({\rm {HydroLight}}; \theta_{\rm {s}} = 0,\theta_{\rm {v}} = 0)\,.$ (12)

There should be only a very small difference in HydroLight's $ \pi R_{\rm {rs}}$ and $ [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ , attributable to any differences in the angular distributions of the sky radiances used to compute the Morel BRDF factors and as used in HydroLight. Morel et al. (2002) (page 6295) notes that these differences are negligible.

It should be noted that HydroLight $ R_{\rm {rs}}$ values are valid for whatever IOPs were used in the run; there is no restriction to homogeneous Case 1 water and no need to estimate a chlorophyll concentration as must be done when applying the Morel BRDF factors. Thus HydroLight can give a very general $ [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ via Eq. ([*]), without the assumptions made by Morel et al. when developing the BRDF correction factors seen in Eq. ([*]).

Figure [*] illustrates the magnitudes of the corrections to $ \pi R_{\rm {rs}}(\theta_{\rm {s}}, \theta_{\rm {v}}, \phi)$ values. HydroLight was first run to generate remote-sensing reflectances for a Sun zenith angle of $ \theta_{\rm {s}} = 50$ deg, with the Sun in a clear sky. The water IOPs and chlorophyll values were the same as those used for Fig. [*]. The runs included Raman scatter by water and fluorescence by chlorophyll and CDOM. The dashed lines in Fig. [*] show the values of $ \pi R_{\rm {rs}}(\theta_{\rm {s}} = 50, \theta_{\rm {v}} = 30, \phi = 90)$ . These viewing angles correspond to a sensor viewing the ocean at an off-nadir angle of 30 deg at right angles to the solar plane. The HydroLight runs were then repeated with the Sun at the zenith. The solid lines in the figures show the resulting values of $ [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ as determined by Eq. ([*]). Depending on the water IOPs, Sun zenith angle, viewing direction, and wavelength, the normalization can change a spectrum by tens of percent, or almost not at all. The changes tend to be greatest in high-chlorophyll waters, at large solar zenith angles, and at large off-nadir viewing angles.

Figure: Comparison of exact normalized water-leaving reflectances $ [\rho_w]_{\rm{N}}^{\rm{ex}}$ (solid lines) with unnormalized reflectances (dashed lines) for $ \theta_{\rm{s}} = 50, \theta_{\rm{v}} = 30, \phi = 90$ . The color identifies the chlorophyll values of 0.05, 0.5, and 5 $ \unit{mg~m^{-3}}$ .
Image 6b83da644e4767985fec5ecaf60dc358

When formulated in terms of reflectances, the partitioning of Eq. ([*]) becomes (e.g., Gordon and Wang (1994a))

$\displaystyle \rho_{\rm {t}} = \rho_{\rm {R}} + [\rho_{\rm {a}} + \rho_{\rm {Ra}} ] + T \rho_{\rm {g}} + t \rho_{\rm {wc}} + t \rho_{\rm {w}} \,,$ (13)

where the terms correspond to those of Eq. ([*]). A similar equation applies to the reflectance form of Eq. ([*]).

Figure [*] shows the radiances of Fig. [*] recast as normalized reflectances. It should be noted in Fig. [*] that the solar-spectrum features (most noticeable below 600 nm) seen in the TOA $ E_{\rm {d}}$ spectrum of Fig. [*] are removed by the normalization process. However, the atmospheric absorption features (most noticeable beyond 600 nm) remain in the TOA reflectances, but are not present in the surface reflectances. Thus the dotted curves for $ \rho_{\rm {t}}$ , $ \rho_{\rm {g}}$ , and $ \rho_{\rm {w}}$ , are very smooth functions of wavelength. The surface glint reflectance spectrum $ \rho_{\rm {g}}$ is almost independent of wavelength, whereas the surface glint radiance seen in Fig. [*] depends of wavelength in the same manner as the surface irradiance spectrum. (Note that the glint reflectance $ \rho_{\rm {g}}$ considered here is not the same at the surface radiance reflectance factor $ \rho$ of Mobley (1999) Eq. 4, and Mobley (2015), which is a ratio of reflected to incident sky radiances for the given wind speed, and Sun and viewing directions.)

Figure: Normalized reflectances corresponding to the radiances of Fig. [*]. These reflectances are $ [\rho]_{\rm{N}}$ , not $ [\rho]_{\rm{N}}^{\rm{ex}}$ .
Image 63b0aec416f69e7c447b2d661ff95f99

Summary

The interpretation of $ [\rho_{\rm {w}}]_{\rm {N}}^{\rm {ex}}$ as developed above can be summarized as follows:

  • Start with a measured $ L_{\rm {w}}(\theta_{\rm {s}}, \theta_{\rm {v}}, \phi)$ . Then
  • The division by $ \cos \theta_{\rm {s}}$ in Eq. ([*]) moves the Sun to the zenith.
  • The division by $ t(\theta_{\rm {s}})$ in Eq. ([*]) rescales the radiance magnitude to account for atmospheric attenuation, but the angular distribution of the radiance is unchanged. This diffuse transmission is for the actual atmosphere at the time of observation.
  • The $ (R/R_{\rm {o}})^2$ factor in Eq. ([*]) corrects for the Earth-Sun distance.
  • The BRDF factors in Eq. ([*]) normalize $ L_{\rm {w}}$ as measured for the actual atmospheric and in-water radiance distributions to what $ L_{\rm {w}}$ would be for a reference atmospheric radiance distribution and for Case 1 water with the given chlorophyll value.

The proposed atmospheric retrieval accuracy requirements for the PACE (Pre-Aerosol, Clouds, and ocean Ecosystem) mission are (Franz (2015))

  • for $ [\rho]_{\rm{N}}^{\rm{ex}}$ in the 350-400 nm range, a maximum error of 0.002 or 10%
  • for $ [\rho]_{\rm{N}}^{\rm{ex}}$ in the 400-710 nm range, a maximum error of 0.001 or 5%
In terms of $ R_{\rm {rs}} = [\rho]_{\rm {N}}^{\rm {ex}}/\pi$ , this gives
  • for $ R_{\rm {rs}}$ in the 350-400 nm range, a maximum error of $ 6 \times 10^{-4}~\unit{sr^{-1}}$
  • for $ R_{\rm {rs}}$ in the 400-710 nm range, a maximum error of $ 3\times10^{-4}~\unit{sr^{-1}}$
Given that the water-leaving radiance $ L_{\rm {w}}$ is at most 10% of the TOA radiance (e.g., Fig. [*]), a requirement for a 5% maximum error in $ L_{\rm {w}}$ (expressed as the same percentage error in $ [\rho]_{\rm{N}}^{\rm{ex}}$ or $ R_{\rm {rs}}$ ) implies a maximum error of roughly 0.5% in the measured TOA radiance $ L_{\rm {t}}$ . Such a small error in the TOA radiance cannot be achieved by pre-launch sensor radiometric calibration alone. Thus ocean color sensors require post-launch "vicarious calibration." This term refers to the process of (1) performing atmospheric correction to recover sea-level $ R_{\rm {rs}}$ spectra from the measured $ L_{\rm {t}}$ ; (2) propagating the corresponding sea-level $ L_{\rm {w}}$ back to the TOA (including the atmospheric path radiance contributions) to obtain a TOA radiance $ L_{\rm {t}}^{\rm {AC}}$ for comparison with the sensor-measured TOA radiance $ L_{\rm {t}}$ . The difference between the retrieved $ L_{\rm {t}}^{\rm {AC}}$ and the measured $ L_{\rm {t}}$ then (3) gives a final set of sensor calibration correction factors to be applied to the measured TOA radiances so that they lead, after atmospheric correction, to the proper sea-level water-leaving radiances and associated reflectances.

The $ Q$ part of the Morel et al. correction has been validated against empirical data for Case 1 waters by Voss et al. (2007) and found to give good agreement with measured radiance distributions. As concluded there, "...the bidirectional corrections based on the lookup tables generated from the model, and presently applied to ocean color imagery, are sound and amply validated for Case 1 waters..." However, it should be remembered that the above BRDF correction is based on a particular atmospheric model and on particular bio-optical models for Case 1 water. For different atmospheric conditions, the $ f/Q$ correction would be different, although the difference would probably be small. However, for different water IOPs, in particular for Case 2 waters with high mineral particle loads or high concentrations of CDOM, the differences in the water-column scattering and absorption properties could have a significant effect on the $ f/Q$ correction. That variability has not yet been studied, and the above correction remains the current state of the science and is implemented by the OBPG for operational ocean color image processing.