Page updated: April 12, 2020
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 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

cos𝜃dL(z,𝜃,ϕ,λ) dz = c(z,λ)L(z,𝜃,ϕ,λ) + 02π0πL(z,𝜃,ϕ,λ)β(z,𝜃,ϕ 𝜃,ϕ,λ)sin𝜃d𝜃dϕ (1)  over all directions. Dropping the wavelength argument for brevity and writing the differential element of solid angle sin𝜃d𝜃dϕ as dΩ(𝜃,ϕ), the left hand side of Eq. (1) yields

02π0π cos𝜃dL(z,𝜃,ϕ) dz dΩ(𝜃,ϕ) = d dz02π0πL(z,𝜃,ϕ)cos𝜃dΩ(𝜃,ϕ) = d dz[Ed(z) Eu(z)] (2) after noting that cos𝜃 < 0 for π2 < 𝜃 π and recalling the definitions of the upwelling and downwelling plane irradiances as integrals of the radiance. The cL term becomes

c(z)L(z,𝜃,ϕ)dΩ(𝜃,ϕ) = c(z)L(z,𝜃,ϕ)dΩ(𝜃,ϕ) = c(z)Eo(z), (3) where the double integration over all directions is the same as shown in Eq. (2), and Eo(z) is the scalar irradiance. The elastic scatter path function gives

. L(z,𝜃,ϕ)β(z,𝜃,ϕ 𝜃,ϕ)dΩ(𝜃,ϕ)dΩ(𝜃,ϕ) = L(z,𝜃,ϕ) β(z,𝜃,ϕ 𝜃,ϕ)dΩ(𝜃,ϕ)dΩ(𝜃,ϕ) = b(z)L(z,𝜃,ϕ)dΩ(𝜃,ϕ) = b(z)Eo(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

d dz[Ed Eu] = cEo + bEo,


d dz[Ed(z,λ) Eu(z,λ)] = a(z,λ)Eo(z,λ)(Wm3nm1), (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 Ed Eu to the absorption coefficient a and the scalar irradiance Eo. 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 Ed, Eu, and Eo. 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 Ed, Eu, and Eo, 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.

Comments for Gershun's Law:

Loading Conversation