**Page updated:**
January 9, 2021 **Author:** Curtis Mobley

View PDF

# The Volume Scattering Function

We now take into account the angular distribution of the scattered power, with two assumptions. We ﬁrst assume that the medium is isotropic, i.e., its inﬂuence on light is the same in all directions at a given point. This is a reasonable assumption for natural waters in which the particles are randomly oriented by turbulence. We also assume that the incident light is unpolarized. If these two assumptions are true, then the scattering process is azimuthally symmetric. This means that the scattering depends only on the scattering angle $\psi $, which is measured from the direction of the unscattered beam as shown in Fig. 1. Clearly 0 $\le \psi \le \pi $.

### Conceptual Deﬁnition of the Volume Scattering Function

With these two assumptions, the fraction of incident power scattered out of the beam through an angle $\psi $ into a solid angle $\Delta \Omega $ centered on $\psi $, is ${\Delta}^{2}{\Phi}_{s}\left(\psi ,\lambda \right)\u2215{\Phi}_{i}\left(\lambda \right)$. The notation ${\Delta}^{2}{\Phi}_{s}\left(\psi ,\lambda \right)$ is adopted from Gordon (1994) to indicate that ${\Delta}^{2}{\Phi}_{s}\left(\psi ,\lambda \right)$ is a second-order small quantity, that is, it is small because both $\Delta r$ and $\Delta \Omega $ are small. The solid angle $\Delta \Omega $ now includes all directions within the two red rings shown in the ﬁgure, corresponding to all directions between scattering angles $\psi $ and $\psi +\Delta \psi $. The volume scattering function $\beta \left(\psi ,\lambda \right)$ is deﬁned as the limit of this fraction as $\Delta r\to 0$ and $\Delta \Omega \to 0$:

Recalling the deﬁnition of spectral radiant intensity as scattered power per unit solid angle, the corresponding intensity scattered into the given solid angle $\Delta \Omega $ is $\Delta {I}_{s}\left(\psi ,\lambda \right)={\Delta}^{2}{\Phi}_{s}\left(\psi ,\lambda \right)\u2215\Delta \Omega $. Moreover, if the incident power ${\Phi}_{i}\left(\lambda \right)$ falls on an area $\Delta A$, then the corresponding incident irradiance is ${E}_{i}\left(\lambda \right)$ = ${\Phi}_{i}\left(\lambda \right)\u2215\Delta A$. Noting that $\Delta V=\Delta r\Delta A$ is the volume of water that is illuminated by the incident beam allows Eq. (1) to be rewritten as

This form of $\beta \left(\psi ,\lambda \right)$ suggests the name volume scattering function (commonly abbreviated as VSF) and the physical interpretation of scattered intensity per unit incident irradiance per unit volume of water. In the language of a physicist, the VSF also can be interpreted as the diﬀerential scattering cross section per unit volume.

The forms of Eqs. (1) and (2) involving partial derivatives are conceptual deﬁnitions. Instruments that measure the VSF always have to deal with ﬁnite path lengths $\Delta r$ and ﬁnite solid angles $\Delta \Omega $. The transition from a deﬁnition as a derivative to an instrument design is discussed on the Measuring IOPs page.

### Other Measures of Scattering

Integrating $\beta \left(\psi ,\lambda \right)$ over all directions (solid angles) gives the total scattered power per unit incident irradiance and unit volume of water, in other words the scattering coeﬃcient:

The $2\pi $ in the last equation follows from our assumption that the scattering is azimuthally symmetric about the incident direction. This integration is often divided into forward scattering, 0 $\le \psi <\pi \u22152$, and backward scattering, $\pi \u22152\le \psi \le \pi $, parts. The corresponding forward and backward scattering coeﬃcients are, respectively,

$${b}_{f}\left(\lambda \right)\equiv 2\pi {\int}_{0}^{\pi \u22152}\phantom{\rule{0.3em}{0ex}}\beta \left(\psi ,\lambda \right)\phantom{\rule{0.3em}{0ex}}sin\psi \phantom{\rule{0.3em}{0ex}}d\psi $$ | (4) |

The backscatter fraction, deﬁned by

$${B}_{b}\left(\lambda \right)\equiv \frac{{b}_{b}\left(\lambda \right)}{b\left(\lambda \right)}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ | (6) |

gives the fraction of scattered light that is deﬂected through scattering angles greater than 90 deg. This quantity is fundamental to remote sensing because most of the light leaving the ocean in an upward direction comes from sunlight that was originally heading downward, but which was backscattered into upward directions.

The volume scattering phase function, $\stackrel{\u0303}{\beta}\left(\psi ,\lambda \right)$, is deﬁned by

Writing the VSF $\beta \left(\psi ,\lambda \right)$ as the product of the scattering coeﬃcient $b\left(\lambda \right)$ times the phase function $\stackrel{\u0303}{\beta}\left(\psi ,\lambda \right)$ partitions $\beta \left(\psi ,\lambda \right)$ into a factor giving the magnitude of the total scattering, $b\left(\lambda \right)$ with units of ${m}^{-1}$, and a factor giving the angular distribution of the scattered light, $\stackrel{\u0303}{\beta}\left(\psi ,\lambda \right)$, with units of ${sr}^{-1}$. Combining Eqs. (3) and (7) gives the normalization condition for the phase function:

This normalization implies that the backscatter fraction can be computed from

Comment: The term “phase function” has its historical origins in astronomy, where the phase angle is the angle between the directions of the light incident onto and reﬂected from an object. In the case of sunlight being reﬂected by the Moon and seen on Earth, this angle depends on the phase of the Moon—hence the name phase function for the function that gives the phase angle as a function of time. In our terminology, the lunar phase angle is $\pi $ minus the scattering angle. The term phase function is rather misleading because the phase function $\stackrel{\u0303}{\beta}\left(\psi ,\lambda \right)$ deﬁned here has nothing whatsoever to do with the phase of an electromagnetic wave or with its historical origin tied to the phase of the Moon. The Russian literature calls $\stackrel{\u0303}{\beta}$ the “scattering diagram,” which is more descriptive of what it really is. You also see the term “indicatrix.” Neither of these terms is commonly used in the English-language literature.

The asymmetry parameter g or mean cosine of the phase function is average over all scattering directions of the cosine of the scattering angle $\psi $, namely

The asymmetry parameter is a convenient measure of the “shape” of the phase function. For example, if $\stackrel{\u0303}{\beta}\left(\psi \right)$ is very large for small $\psi $, then $g$ is near one. If $\stackrel{\u0303}{\beta}\left(\psi \right)$ is symmetric about $\psi $ = 90 degrees , then $g=0$. Typical ocean waters have $g$ values in the range of 0.8 to 0.95.

Comment: We have assumed above that the scattering is azimuthally symmetric, so that the directional pattern or angular shape of the VSF depends only on the scattering angle $\psi $. This is not the case for polarized incident light, even if the medium is isotropic. Thus a linearly polarized laser beam will scatter diﬀerently for diﬀerent azimuthal directions (measured relative to the plane of polarization), even if the medium is isotropic. Description of polarized light requires the Stokes Vector formalism and a more complicated description of scattering; see Polarization: Scattering Geometry and the The General Vector RTE.

If the medium contains non-spherical particles that are not randomly oriented, even unpolarized light will scatter diﬀerently for diﬀerent azimuthal directions. This is sometimes the case, for example, in cirrus clouds and ice fogs, which can have non-spherical ice crystals that become oriented in a particular way as the crystals fall through a calm atmosphere. The atmosphere is then an optically anisotropic medium, and scattering is not azimuthally symmetric. Scattering of sunlight in such an atmosphere gives rise to phenomena such as “sun dogs” or parhelia, which are bright spots to either side of the sun.