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

level 2

Gershun's Law

Page updated: Jun 20, 2017
Principal author: Curtis Mobley

The radiative transfer equation is a statement of energy conservation in the sense that it accounts for all the losses and gains to a swarm of photons moving through the water along a path in a fixed direction. We now derive a useful conservation statement that holds at a fixed point in the water, through which photons are moving in all directions.

The desired result is obtained by integrating the 1-D, time-independent, source-free RTE

$\displaystyle \cos \theta \frac{d L (z, \theta , \phi , \lambda )}{d z} ~=$ $\displaystyle ~- c(z,\lambda )\,L(z, \theta , \phi , \lambda )$    
$\displaystyle +$ $\displaystyle ~\int_{0}^{2 \pi } \int_{0}^{\pi} L(z,\theta' ,\phi' , \lambda ) \, \beta (z,\theta' ,\phi' \rightarrow \theta , \phi ,\lambda ) \sin \theta' d \theta' d \phi'$ (1)

over all directions. Dropping the wavelength argument for brevity and writing the differential element of solid angle $ \sin \theta d \theta d \phi$ as $ d\Omega(\theta, \phi)$, the left hand side of Eq. (1) yields

$\displaystyle \int_{0}^{2 \pi } \int_{0}^{\pi} \cos \theta {\frac{d L(z,\theta ,\phi )}{d z}} d\Omega(\theta ,\phi)~=$ $\displaystyle ~
 {\frac{d}{d z}} \int_{0}^{2 \pi } \int_{0}^{\pi} L(z,\theta ,\phi ) \cos \theta d\Omega (\theta ,\phi )$    
$\displaystyle =$ $\displaystyle ~{\frac{~d}{d z}} [E _{\rm d}(z) - E _{\rm u}(z) ]$ (2)

after noting that $ \cos \theta < 0$ for $ \pi /2 < \theta \leq \pi$ and recalling the definitions of the upwelling and downwelling plane irradiances as integrals of the radiance. The $ {-cL}$ term becomes

$\displaystyle \iint - c(z) L(z,\theta ,\phi ) d\Omega (\theta ,\phi )~=$ $\displaystyle ~- c(z) \iint L(z,\theta ,\phi) d\Omega (\theta ,\phi)$    
$\displaystyle =$ $\displaystyle ~-c(z) E _{\rm o} (z)\,\, ,$ (3)

where the double integration over all directions is the same as shown in Eq. (2), and $ E_{\rm o}(z)$ is the scalar irradiance. The elastic scatter path function gives

$\displaystyle .$ $\displaystyle ~\iint \left [ \iint L(z,\theta ',\phi') \, \beta (z,\theta ', \phi ' \rightarrow \theta, \phi ) d\Omega (\theta ',\phi ' ) \right ] d\Omega(\theta, \phi)$    
$\displaystyle ~=$ $\displaystyle ~\iint L(z,\theta ',\phi') \left [ \iint \beta (z, \theta ', \phi ' \rightarrow \theta, \phi ) d\Omega (\theta, \phi ) \right ] d\Omega (\theta ', \phi ' )$    
$\displaystyle =$ $\displaystyle ~b(z) \iint L(z, \theta ', \phi ') d \Omega (\theta ',\phi ' )$    
$\displaystyle =$ $\displaystyle ~b(z) E _{\rm o} (z)\,\, .$ (4)

Here we recall that the integral of the volume scattering function over all directions is the scattering coefficient.

Collecting terms (2)-(4) resulting from the directional integration of the RTE, we have

$\displaystyle {\frac{d}{d z}} \, [E _{\rm d} - E _{\rm u} \, ]~=~- c E _{\rm o} + b E _{\rm o} \,\, ,


$\displaystyle {\frac{d}{d z}} [E _{\rm d}(z, \lambda) - E _{\rm u}(z, \lambda) \, ]~=~- a(z, \lambda) E _{\rm o}(z, \lambda)~~~~~{\rm {(W~m^{-3}~nm ^{-1})}\,\, ,}$ (5)

which is the desired result. This equation is known as Gershun's law (Gershun (1936), Gershun (1939)).

The physical significance of Eq. (5) is that it relates the depth rate of change of the net irradiance $ E_{{\rm d}} - E_{{\rm u}}$ to the absorption coefficient a and the scalar irradiance $ E_{{\rm o}}$. If inelastic scattering (fluorescence and Raman scattering) and internal sources (such as bioluminescence) are negligible at the wavelength of interest, then Eq. (5) can be used to obtain the absorption coefficient a from in situ measurements of the irradiances $ E_{{\rm d}}$, $ E_{{\rm u}}$, and $ E_{{\rm o}}$. This is an example of an inverse model - a model that retrieves an inherent optical property from measurements of the light field.

Voss (1989) used Gershun's law (5) to recover a values to within an estimated error of order 20%. Inelastic scattering and internal source effects were reasonably assumed to be negligible in his study. The needed irradiances were all computed from a measured radiance distribution, so that no intercalibration of instruments was required.

A more general development can be made to account for internal sources or inelastic scatter and for 3-D and time-dependent light fields, as shown in Light and Water section 5.10. The result is known as the divergence law for irradiance.

Maffione et al. (1993) determined absorption values by writing the source-free form of the 3-D divergence law in spherical coordinates and applying the result to irradiance measurements made using an underwater, artificial, isotropic light source. The artificial light source allowed measurements to be made at night, thus there was no inelastic scattering from other wavelengths. Their instrument did not require absolute radiometric calibration.

Note, however, Gershun's law will give incorrect absorption values if naively applied to waters and wavelengths where inelastic processes such as Raman scattering or fluorescence are significant. For this reason, and because of calibration difficulties if different instruments are used to measure $ E_{\rm d}$, $ E_{\rm u}$, and $ E_{\rm o}$, Gershun's law is seldom used as a way to measure absorption. Nevertheless, it is sometimes a useful check on the internal consistency of numerical models or measured data, and it leads to a convenient way of calculating radiant heating rates.