**Page updated:**
May 19, 2021 **Author:** Curtis Mobley

View PDF

# Gershun's Law

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

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

$$\begin{array}{llll}\hfill cos\mathit{\theta}\frac{dL\left(z,\mathit{\theta},\varphi ,\lambda \right)}{dz}=& -c\left(z,\lambda \right)\phantom{\rule{0.3em}{0ex}}L\left(z,\mathit{\theta},\varphi ,\lambda \right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {\int}_{0}^{2\pi}{\int}_{0}^{\pi}L\left(z,{\mathit{\theta}}^{\prime},{\varphi}^{\prime},\lambda \right)\phantom{\rule{0.3em}{0ex}}\beta \left(z,{\mathit{\theta}}^{\prime},{\varphi}^{\prime}\to \mathit{\theta},\varphi ,\lambda \right)sin{\mathit{\theta}}^{\prime}d{\mathit{\theta}}^{\prime}d{\varphi}^{\prime}\phantom{\rule{2em}{0ex}}& \hfill \text{(1)}\end{array}$$ over all directions. Dropping the wavelength argument for brevity and writing the diﬀerential element of solid angle $sin\mathit{\theta}d\mathit{\theta}d\varphi $ as $d\Omega \left(\mathit{\theta},\varphi \right)$, the left hand side of Eq. (1) yields $$\begin{array}{llll}\hfill {\int}_{0}^{2\pi}{\int}_{0}^{\pi}cos\mathit{\theta}\frac{dL\left(z,\mathit{\theta},\varphi \right)}{dz}d\Omega \left(\mathit{\theta},\varphi \right)=& \frac{d}{dz}{\int}_{0}^{2\pi}{\int}_{0}^{\pi}L\left(z,\mathit{\theta},\varphi \right)cos\mathit{\theta}d\Omega \left(\mathit{\theta},\varphi \right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill =& \frac{\phantom{\rule{1em}{0ex}}d}{dz}\left[{E}_{d}\left(z\right)-{E}_{u}\left(z\right)\right]\phantom{\rule{2em}{0ex}}& \hfill \text{(2)}\end{array}$$

after noting that $cos\mathit{\theta}<0$ for $\pi \u22152<\mathit{\theta}\le \pi $ and recalling the deﬁnitions of the upwelling and downwelling plane irradiances as integrals of the radiance. The $-cL$ term becomes

$$\begin{array}{llll}\hfill \iint -c\left(z\right)L\left(z,\mathit{\theta},\varphi \right)d\Omega \left(\mathit{\theta},\varphi \right)=& -c\left(z\right)\iint L\left(z,\mathit{\theta},\varphi \right)d\Omega \left(\mathit{\theta},\varphi \right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill =& -c\left(z\right){E}_{o}\left(z\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill \text{(3)}\end{array}$$ where the double integration over all directions is the same as shown in Eq. (2), and ${E}_{o}\left(z\right)$ is the scalar irradiance. The elastic scatter path function gives

$$\begin{array}{llll}\hfill .& \iint \left[\iint L\left(z,{\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)\phantom{\rule{0.3em}{0ex}}\beta \left(z,{\mathit{\theta}}^{\prime},{\varphi}^{\prime}\to \mathit{\theta},\varphi \right)d\Omega \left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)\right]d\Omega \left(\mathit{\theta},\varphi \right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill =& \iint L\left(z,{\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)\left[\iint \beta \left(z,{\mathit{\theta}}^{\prime},{\varphi}^{\prime}\to \mathit{\theta},\varphi \right)d\Omega \left(\mathit{\theta},\varphi \right)\right]d\Omega \left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill =& b\left(z\right)\iint L\left(z,{\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)d\Omega \left({\mathit{\theta}}^{\prime},{\varphi}^{\prime}\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill =& b\left(z\right){E}_{o}\left(z\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(4)}\end{array}$$ Here we recall that the integral of the volume scattering function over all directions is the scattering coeﬃcient.

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

$$\frac{d}{dz}\phantom{\rule{0.3em}{0ex}}\left[{E}_{d}-{E}_{u}\phantom{\rule{0.3em}{0ex}}\right]=-c{E}_{o}+b{E}_{o}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

or

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

The physical signiﬁcance of Eq. (5) is that it relates the depth rate of change of the net irradiance ${E}_{d}-{E}_{u}$ to the absorption coeﬃcient $a$ and the scalar irradiance ${E}_{o}$. If inelastic scattering (ﬂuorescence 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 coeﬃcient $a$ from in situ measurements of the irradiance triplet ${E}_{d},{E}_{u}$, and ${E}_{o}$. This is an example of an explicit inverse model—a model that retrieves an inherent optical property from measurements of the light ﬁeld.

Voss (1989) used Gershun’s law (5) to recover $a$ values to within an estimated error of order 20%. Inelastic scattering and internal source eﬀects 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 ﬁelds, as shown in Light and Water section 5.10. The result is known as the divergence law for irradiance:

$$\frac{1}{v}\frac{\partial {E}_{o}}{\partial t}+\nabla \cdot \overrightarrow{E}=-a{E}_{o}+{E}_{o}^{S}\phantom{\rule{0.3em}{0ex}}.$$ | (6) |

Here $v$ is the speed of light in the water, $\overrightarrow{E}$ is the vector irradiance, and ${E}_{o}^{S}$ is a source term. For time-independent, 1D, source-free water, this equation reduces to Eq. (5).

Maﬃone 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, artiﬁcial, isotropic light source. The artiﬁcial 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 ﬂuorescence are signiﬁcant. For this reason, and because of calibration diﬃculties if diﬀerent instruments are used to measure ${E}_{d},{E}_{u}$, and ${E}_{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.

### Heating Rates in the Upper Ocean

Gershun’s law has a very important application in the computation of heating rates in the upper ocean. The rate of heating of water depends on how much scalar irradiance is available and on the total absorption coeﬃcient of the water (pure water plus all other constituents). Combining the First Law of Thermodynamics (i.e., conservation of energy) with Eq. (5) gives

$$\frac{\partial T}{\partial t}\left(z,\text{at}\phantom{\rule{2.6108pt}{0ex}}\lambda \right)=\frac{1}{{c}_{v}\rho}\phantom{\rule{0.3em}{0ex}}a\left(z,\lambda \right){E}_{o}\left(z,\lambda \right)=-\frac{1}{{c}_{v}\rho}\frac{\partial \left[{E}_{d}\left(z,\lambda \right)-{E}_{u}\left(z,\lambda \right)\right]}{\partial z}\phantom{\rule{1em}{0ex}}\left[\frac{deg\phantom{\rule{2.6108pt}{0ex}}C}{sec}\right]\phantom{\rule{0.3em}{0ex}}.$$ |

Here

- $T$ is the temperature in deg C
- $t$ is the time in seconds
- $\frac{\partial T}{\partial t}\left(z,\text{at}\phantom{\rule{2.6108pt}{0ex}}\lambda \right)$ denotes the rate of change of temperature at depth $z$ due to energy absorbed at wavelength $\lambda $. It does not mean that temperature is a function of wavelength.
- ${c}_{v}=3900\phantom{\rule{2.6108pt}{0ex}}J\phantom{\rule{2.6108pt}{0ex}}{\left(kg\phantom{\rule{2.6108pt}{0ex}}deg\phantom{\rule{2.6108pt}{0ex}}C\right)}^{-1}$ is the speciﬁc heat of sea water at constant volume
- $\rho =1025\phantom{\rule{2.6108pt}{0ex}}kg\phantom{\rule{2.6108pt}{0ex}}{m}^{-3}$ is the density of sea water

This equation is usually applied with wavelength-integrated irradiances. In ocean modeling, “short-wave” radiation is usually taken to be the range of 400-1000 nm. Letting

$${E}_{d,u}\left(z\right)={\int}_{400}^{1000}{E}_{d,u}\left(z,\lambda \right)d\lambda $$ |

gives

In optically deep water, ${E}_{u}\left(z\right)<<{E}_{d}\left(z\right)$, and the upwelling irradiance is usually dropped. The downwelling irradiance is then modeled, e.g. in terms of the chlorophyll concentration in Case 1 water. The EcoLight-S radiative transfer code (Mobley (2011)) was developed to provide extremely fast solutions of the radiative transfer equation to obtain ${E}_{d}\left(z\right)$ and ${E}_{u}\left(z\right)$ for any water body (deep or shallow, Case 1 or Case 2) for use in Eq. (7) in coupled physical-biological-optical ecosystem models. (EcoLight-S also computes other quantities such as PAR and remote-sensing reﬂectance, which are inputs to photosynthesis calculations or are useful for model validation.) For an example of such a model, see Mobley et al. (2015).