**Page updated:**
April 12, 2020 **Author:** Curtis Mobley

View PDF

# Energy Conservation

When making radiance or irradiance measurements in the ocean, it is sometimes hard to verify that energy is being conserved. This is equally true when looking at the output from a numerical model such as HydroLight. Confusion most often arises when comparing measurements made just above and below the sea surface and observing that the measured (or predicted) radiances do not seem to ”add up” correctly. To resolve such confusion, ﬁrst note that it’s the Law of Conservation of Energy, not the law of conservation of radiance or irradiance. We must therefore express the law of conservation of energy in terms of radiance or irradiance as appropriate to the measurements being made.

#### Within the Water

Within the water column, if there are no internal sources such as bioluminescence or inelastic scatter contributing to the wavelength of interest, conservation of energy is expressed in terms of irradiances by Gershun’s law:

Here $a$ is the absorption coeﬃcient, ${E}_{o}$ is the scalar irradiance, and ${E}_{d}$ and ${E}_{u}$ are the downwelling and upwelling plane irradiances. This equation can be illustrated using numerical results generated by HydroLight. A run was made with the sun placed at a 60 deg zenith angle in a clear sky; the wind speed was $5\phantom{\rule{1em}{0ex}}m\phantom{\rule{1em}{0ex}}{s}^{-1}$, and the water was inﬁnitely deep. The absorption coeﬃcient was set to $a=0.1\phantom{\rule{1em}{0ex}}{m}^{-1}$, the scattering coeﬃcient to $b=0.4\phantom{\rule{1em}{0ex}}{m}^{-1}$, and an ”average-particle” phase function was chosen. (The run used HydroLight’s ”abconst” IOP model for these inputs.) No internal sources were included in the run. Table 1 shows the computed irradiances at depths of 5.00 and 5.01 m. Closely spaced depths are used to illustrate Eq. (1) so that the depth derivative can be accurately approximated as a ﬁnite diﬀerence and the scalar irradiance can be taken as the average of the values at 5.00 and 5.01 m.

Using these values in the ﬁnite diﬀerence approximation of Eq. (1) gives

$$\begin{array}{llll}\hfill a=& -\phantom{\rule{0.3em}{0ex}}\frac{1}{0.5\left(0.51833+0.51749\right)}\phantom{\rule{0.3em}{0ex}}\frac{\left(0.32559-0.013129\right)-\left(0.32613-0.013151\right)}{5.01-5.00}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill \phantom{\rule{1em}{0ex}}=& \phantom{\rule{1em}{0ex}}0.1000\phantom{\rule{1em}{0ex}}{m}^{-1}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

The absorption coeﬃcient as computed from the irradiances via Gershun’s law is the same as the one used as input to the run, so HydroLight is conserving energy in its in-water solution of the radiative transfer equation. (HydroLight also conserves energy if internal sources are present. Demonstration of that requires adding an appropriate source term to Gershun’s equation; see Light and Water (1994) Eq. 5.34.) Such accuracy cannot be expected with real data because of imperfect instrument calibration and other sources of error such as Raman scatter (an internal source not accounted for in Eq. 1), which will always be present to some degree during daylight measurements. Nevertheless, Gershun’s law can be a useful way to retrieve absorption coeﬃcients from in-water radiance or irradiance measurements, as seen in papers by Voss (1989) and Maﬃone et al. (1993).

#### Across the Air-Water Surface

The situation is more complicated when considering energy conservation across the air-water surface. Consider a collimated beam of light incident onto a level sea surface. Part will be reﬂected and part will be transmitted into the water. Figure 1 shows this geometry and various quantities needed to express conservation of energy at the surface.

Conservation of energy (or equivalently, power) across the surface requires that

$${P}_{i}={P}_{r}+{P}_{t}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ | (2) |

where ${P}_{i}$ is the power in an incident collimated beam of cross section ${A}_{i}$, ${P}_{r}$ is the power in the reﬂected beam of cross section ${A}_{r}$, and ${P}_{t}$ is the power in transmitted beam of cross section ${A}_{t}$. These beam cross sectional areas are shown in green in Fig. 1. The beam of incident power illuminates a horizontal area ${A}_{s}$ of the level sea surface (shown in blue).

##### Conservation in terms of plane irradiances

Conservation of energy can be expressed in terms of incident (downwelling), reﬂected (upwelling), and transmitted (downwelling) plane irradiances using the horizontal projections of the beam cross section areas. Recalling that plane irradiance is power per unit horizontal area, Eq. (2) becomes

$$\begin{array}{llll}\hfill {A}_{p}{E}_{d}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {A}_{p}{E}_{u}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{p}{E}_{d}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

where the horizontally projected areas ${A}_{p}$ corresponding to each beam are shown in yellow in Fig. 1. Since each beam has the same horizontally projected area, this equation reduces to just

$$\begin{array}{llll}\hfill {E}_{d}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {E}_{u}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {E}_{d}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(3)}\end{array}$$

Equation (3) is conservation of energy across the surface expressed in terms of plane irradiances for a collimated incident beam.

##### Conservation in terms of scalar irradiances

To express Eq. (2) in terms of the scalar irradiances, ﬁrst recall how ${E}_{od}$ is obtained from the radiance (Eq 5 of the Geometrical Radiometry page):

$${E}_{od}\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{\iint}_{2{\pi}_{d}}L\left(\widehat{{\xi}^{\prime}}\right)d\Omega \left(\widehat{{\xi}^{\prime}}\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.$$ |

For a collimated incident beam in direction $\widehat{{\xi}_{o}}$, we can write the radiance in terms of the beam irradiance and a Dirac delta function to get

$${E}_{od}\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{\iint}_{2{\pi}_{d}}{E}_{\perp}\delta \left(\widehat{{\xi}^{\prime}}-\widehat{{\xi}_{o}}\right)d\Omega \left(\widehat{{\xi}^{\prime}}\right)\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{E}_{\perp}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

where ${E}_{\perp}$ is the incident irradiance on the beam cross sectional area; ${E}_{\perp}={P}_{i}\u2215{A}_{i}$. Similar equations hold for the reﬂected and transmitted scalar irradiances. These beam cross sectional areas are shown in green in Fig. 1. Thus Eq. (2) becomes

$$\begin{array}{llll}\hfill {A}_{i}{E}_{od}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {A}_{r}{E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{t}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

or

$$\begin{array}{llll}\hfill {E}_{od}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{t}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\frac{{A}_{t}}{{A}_{i}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill \text{(4)}\end{array}$$after noting that ${A}_{r}={A}_{i}$ by the law of reﬂection. By geometry, the ratio of areas ${A}_{t}\u2215{A}_{i}$ is just $cos{\mathit{\theta}}_{t}\u2215cos{\mathit{\theta}}_{i}$. ${\mathit{\theta}}_{t}$ can be written in terms of ${\mathit{\theta}}_{i}$ using Snell’s law, $sin{\mathit{\theta}}_{i}=nsin{\mathit{\theta}}_{t}$, where $n$ is the real index of refraction of the water. Equation (4) then becomes

$$\begin{array}{llll}\hfill {E}_{od}\left(\phantom{\rule{1em}{0ex}}incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{1em}{0ex}}=& \phantom{\rule{1em}{0ex}}{E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& \phantom{\rule{1em}{0ex}}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\frac{\phantom{\rule{1em}{0ex}}1}{cos{\mathit{\theta}}_{i}}\sqrt{1-\frac{{sin}^{2}{\mathit{\theta}}_{i}}{{n}^{2}}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(5)}\end{array}$$

Equation (5) is conservation of energy expressed in terms of scalar irradiance for a collimated incident beam.

To illustrate Equations (3) and (5) with HydroLight output, a run was made with the sun at a zenith angle of ${\mathit{\theta}}_{i}=60\phantom{\rule{1em}{0ex}}deg$ in a black sky, in order to obtain a collimated incident beam. The incident ${E}_{d}$ onto the surface was set to 1. The sea surface was level (zero wind speed), and the water index of refraction was $n=1.34$. The scattering coeﬃcient within in the water was set to zero, so that upwelling light from the water column would not confuse the illustration of energy conservation across the sea surface itself. The absorption coeﬃcient was set to $a=0.1\phantom{\rule{1em}{0ex}}{m}^{-1}$. (These runs used HydroLight IOP model ”abconst” and the ”idealized sky model”).

The computed irradiances just above and below the surface are shown in Table 2. The upwelling irradiances in the water are zero because there was no scattering by the water and the bottom was inﬁnitely deep. These computed irradiances satisfy Equations (3) and (5), showing that HydroLight conserves energy across the surface.

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.12690 | 2.0077 | 0.063209 | 1.0000 |

in water at depth 0 | 0.00000 | 1.2277 | 0.0000 | 0.93693 |

A run was then done with a heavily overcast sky to generate a diﬀuse sky radiance pattern with the sun’s location not discernable. The wind speed was $10\phantom{\rule{1em}{0ex}}m\phantom{\rule{1em}{0ex}}{s}^{-1}$ to generate a non-level sea surface. The resulting irradiances are shown in Table 3.

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.18863 | 1.7764 | 0.050456 | 1.0000 |

in water at depth 0 | 0.00000 | 1.1430 | 0.0000 | 0.94954 |

As before, these plane irradiances satisfy Eq. (3). However, the scalar irradiances do not satisfy Eq. (5). The reason is that Eq. (5) must be applied for each particular incident direction ${\mathit{\theta}}_{i}$ of the diﬀuse sky radiance, in order to correctly compute the angle factor seen in Eq. (5). If this is done one direction at a time (with the sun in a black sky), each individual incident direction satisﬁes Eq. (5) just as in the previous example. Thus energy in conserved across the surface ”beam by beam,” hence in toto, even if it is not obvious from the total scalar irradiances.

A further complication occurs when the water body has scattering. To illustrate this, the scattering coeﬃcient was set to $b=0.4\phantom{\rule{1em}{0ex}}{m}^{-1}$ and an ”average-particle” phase function was used in the previous simulation. The computed irradiance were then

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.21264 | 1.7764 | 0.062372 | 1.0000 |

in water at depth 0 | 0.07922 | 1.2046 | 0.031115 | 0.96873 |

Now even the plane irradiance appears to violate conservation of energy, Eq. (3). The reason is that ${E}_{ou}$ and ${E}_{u}$ in air contain both surface-reﬂected radiance (as before) and radiance transmitted upward through the sea surface. Likewise, ${E}_{od}$ and ${E}_{d}$ in the water contain both downward transmitted irradiance (as before) as well as upwelling radiance that has been reﬂected back downward by the surface. This makes it appear at ﬁrst glance that the irradiances do not add up correctly, even though energy actually is conserved across the air-water surface. The eﬀects of both the surface itself and the water body can be accounted for simply by tallying the plane irradiances incident onto and leaving the surface from above and below:

$$\begin{array}{llll}\hfill & {E}_{d}\left(\phantom{\rule{1em}{0ex}}downwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)-{E}_{u}\left(\phantom{\rule{1em}{0ex}}upwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & {E}_{d}\left(\phantom{\rule{1em}{0ex}}downwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)-{E}_{u}\left(\phantom{\rule{1em}{0ex}}upwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(6)}\end{array}$$

Equation (6) expresses conservation of energy for the surface plus water body, for any conditions of incident lighting and water IOPs. The irradiances seen in Table 4 satisfy this equation.

##### Conservation in terms of radiance

Another complication occurs when expressing conservation of energy across the surface in terms of radiance. Now it is necessary to account for changes in both the beam cross sectional area and in solid angle when going across the surface. When interpreting HydroLight output, it must be remembered that the computed radiances are averages over ”quads” or solid angles of ﬁnite size, as discussed on the HydroLight page. Such quads are represented by the green areas in Fig. 2.

Remembering that radiance is power per unit area per unit solid angle, Eq. (2) can be written as

$${L}_{i}{A}_{i}\Delta {\Omega}_{i}={L}_{r}{A}_{r}\Delta {\Omega}_{r}+{L}_{t}{A}_{t}\Delta {\Omega}_{t}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

or

$${L}_{i}={L}_{r}+{L}_{t}\phantom{\rule{0.3em}{0ex}}\frac{{A}_{t}\phantom{\rule{0.3em}{0ex}}\Delta {\Omega}_{t}}{{A}_{i},\Delta {\Omega}_{i}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

after noting that ${A}_{r}={A}_{i}$ and $\Delta {\Omega}_{r}=\Delta {\Omega}_{i}$ by the law of reﬂection. As noted before, ${A}_{t}\u2215{A}_{i}=cos{\mathit{\theta}}_{t}\u2215cos{\mathit{\theta}}_{i}$. The solid angles are given by $\Delta \Omega =\Delta \left(cos\mathit{\theta}\right)\Delta \varphi $, where the angle ranges are determined by the $\mathit{\theta}$ and $\varphi $ boundaries of the quads. Thus the last equation becomes

since $\Delta {\varphi}_{t}=\Delta {\varphi}_{i}$ for the HydroLight quad layout. These angle factors must be evaluated for the particular quads associated with the incident and transmitted beams.

Consider again the previous case of the sun in a black sky at ${\mathit{\theta}}_{i}=60\phantom{\rule{1em}{0ex}}deg$, a level water surface, and no scattering. The incident and reﬂected HydroLight quads are centered at 60 deg and extend from 55 to 65 deg. Snell’s law gives the transmitted direction as ${\mathit{\theta}}_{t}={sin}^{-1}\left[\left(sin{\mathit{\theta}}_{i}\right)\u2215n\right]=40.26\phantom{\rule{1em}{0ex}}deg$. Thus the transmitted radiance is contained in the HydroLight quad centered at 40 deg and extending from 35 to 45 deg. For the transmitted quad, for example, $\Delta \left(cos{\mathit{\theta}}_{t}\right)=cos35-cos45=0.11205$. Evaluating the angle factors in Eq. (7) for these quads gives

The corresponding computed quad-averaged radiances for this run were ${L}_{i}=50.7996$, ${L}_{r}=3.21094$, and ${L}_{t}=41.8544\phantom{\rule{1em}{0ex}}W\phantom{\rule{1em}{0ex}}{m}^{-2}\phantom{\rule{1em}{0ex}}s{r}^{-1}\phantom{\rule{1em}{0ex}}n{m}^{-1}$. These radiances satisfy Eq. (8). Thus energy expressed as radiance is conserved beam-by-beam across the water surface.

We have seen that conservation of energy across the sea surface can be expressed in various ways in terms of radiance and irradiance. However, those expressions are often written for collimated incident beams and consider only the surface itself (Eqs. 3, 5, and 7). It is therefore not always obvious that energy is conserved for the entire surface plus water column when looking at data or simulations for realistic conditions of sky radiances and water IOPs. In general, Eq. (6) should be used to check conservation of energy across the sea surface with measured or modeled plane irradiances.

When making radiance or irradiance measurements in the ocean, it is sometimes hard to verify that energy is being conserved. This is equally true when looking at the output from a numerical model such as HydroLight. Confusion most often arises when comparing measurements made just above and below the sea surface and observing that the measured (or predicted) radiances do not seem to ”add up” correctly. To resolve such confusion, ﬁrst note that it’s the Law of Conservation of Energy, not the law of conservation of radiance or irradiance. We must therefore express the law of conservation of energy in terms of radiance or irradiance as appropriate to the measurements being made.

#### Within the Water

Within the water column, if there are no internal sources such as bioluminescence or inelastic scatter contributing to the wavelength of interest, conservation of energy is expressed in terms of irradiances by Gershun’s law:

Here $a$ is the absorption coeﬃcient, ${E}_{o}$ is the scalar irradiance, and ${E}_{d}$ and ${E}_{u}$ are the downwelling and upwelling plane irradiances. This equation can be illustrated using numerical results generated by HydroLight. A run was made with the sun placed at a 60 deg zenith angle in a clear sky; the wind speed was $5\phantom{\rule{1em}{0ex}}m\phantom{\rule{1em}{0ex}}{s}^{-1}$, and the water was inﬁnitely deep. The absorption coeﬃcient was set to $a=0.1\phantom{\rule{1em}{0ex}}{m}^{-1}$, the scattering coeﬃcient to $b=0.4\phantom{\rule{1em}{0ex}}{m}^{-1}$, and an ”average-particle” phase function was chosen. (The run used HydroLight’s ”abconst” IOP model for these inputs.) No internal sources were included in the run. Table 1 shows the computed irradiances at depths of 5.00 and 5.01 m. Closely spaced depths are used to illustrate Eq. (1) so that the depth derivative can be accurately approximated as a ﬁnite diﬀerence and the scalar irradiance can be taken as the average of the values at 5.00 and 5.01 m.

Using these values in the ﬁnite diﬀerence approximation of Eq. (1) gives

$$\begin{array}{llll}\hfill a=& -\phantom{\rule{0.3em}{0ex}}\frac{1}{0.5\left(0.51833+0.51749\right)}\phantom{\rule{0.3em}{0ex}}\frac{\left(0.32559-0.013129\right)-\left(0.32613-0.013151\right)}{5.01-5.00}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill \phantom{\rule{1em}{0ex}}=& \phantom{\rule{1em}{0ex}}0.1000\phantom{\rule{1em}{0ex}}{m}^{-1}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

The absorption coeﬃcient as computed from the irradiances via Gershun’s law is the same as the one used as input to the run, so HydroLight is conserving energy in its in-water solution of the radiative transfer equation. (HydroLight also conserves energy if internal sources are present. Demonstration of that requires adding an appropriate source term to Gershun’s equation; see Light and Water (1994) Eq. 5.34.) Such accuracy cannot be expected with real data because of imperfect instrument calibration and other sources of error such as Raman scatter (an internal source not accounted for in Eq. 1), which will always be present to some degree during daylight measurements. Nevertheless, Gershun’s law can be a useful way to retrieve absorption coeﬃcients from in-water radiance or irradiance measurements, as seen in papers by Voss (1989) and Maﬃone et al. (1993).

#### Across the Air-Water Surface

The situation is more complicated when considering energy conservation across the air-water surface. Consider a collimated beam of light incident onto a level sea surface. Part will be reﬂected and part will be transmitted into the water. Figure 1 shows this geometry and various quantities needed to express conservation of energy at the surface.

Conservation of energy (or equivalently, power) across the surface requires that

$${P}_{i}={P}_{r}+{P}_{t}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ | (2) |

where ${P}_{i}$ is the power in an incident collimated beam of cross section ${A}_{i}$, ${P}_{r}$ is the power in the reﬂected beam of cross section ${A}_{r}$, and ${P}_{t}$ is the power in transmitted beam of cross section ${A}_{t}$. These beam cross sectional areas are shown in green in Fig. 1. The beam of incident power illuminates a horizontal area ${A}_{s}$ of the level sea surface (shown in blue).

##### Conservation in terms of plane irradiances

Conservation of energy can be expressed in terms of incident (downwelling), reﬂected (upwelling), and transmitted (downwelling) plane irradiances using the horizontal projections of the beam cross section areas. Recalling that plane irradiance is power per unit horizontal area, Eq. (2) becomes

$$\begin{array}{llll}\hfill {A}_{p}{E}_{d}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {A}_{p}{E}_{u}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{p}{E}_{d}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

where the horizontally projected areas ${A}_{p}$ corresponding to each beam are shown in yellow in Fig. 1. Since each beam has the same horizontally projected area, this equation reduces to just

$$\begin{array}{llll}\hfill {E}_{d}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {E}_{u}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {E}_{d}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(3)}\end{array}$$

Equation (3) is conservation of energy across the surface expressed in terms of plane irradiances for a collimated incident beam.

##### Conservation in terms of scalar irradiances

To express Eq. (2) in terms of the scalar irradiances, ﬁrst recall how ${E}_{od}$ is obtained from the radiance (Eq 5 of the Geometrical Radiometry page):

$${E}_{od}\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{\iint}_{2{\pi}_{d}}L\left(\widehat{{\xi}^{\prime}}\right)d\Omega \left(\widehat{{\xi}^{\prime}}\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.$$ |

For a collimated incident beam in direction $\widehat{{\xi}_{o}}$, we can write the radiance in terms of the beam irradiance and a Dirac delta function to get

$${E}_{od}\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{\iint}_{2{\pi}_{d}}{E}_{\perp}\delta \left(\widehat{{\xi}^{\prime}}-\widehat{{\xi}_{o}}\right)d\Omega \left(\widehat{{\xi}^{\prime}}\right)\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{E}_{\perp}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

where ${E}_{\perp}$ is the incident irradiance on the beam cross sectional area; ${E}_{\perp}={P}_{i}\u2215{A}_{i}$. Similar equations hold for the reﬂected and transmitted scalar irradiances. These beam cross sectional areas are shown in green in Fig. 1. Thus Eq. (2) becomes

$$\begin{array}{llll}\hfill {A}_{i}{E}_{od}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {A}_{r}{E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{t}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

or

$$\begin{array}{llll}\hfill {E}_{od}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{t}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\frac{{A}_{t}}{{A}_{i}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill \text{(4)}\end{array}$$after noting that ${A}_{r}={A}_{i}$ by the law of reﬂection. By geometry, the ratio of areas ${A}_{t}\u2215{A}_{i}$ is just $cos{\mathit{\theta}}_{t}\u2215cos{\mathit{\theta}}_{i}$. ${\mathit{\theta}}_{t}$ can be written in terms of ${\mathit{\theta}}_{i}$ using Snell’s law, $sin{\mathit{\theta}}_{i}=nsin{\mathit{\theta}}_{t}$, where $n$ is the real index of refraction of the water. Equation (4) then becomes

$$\begin{array}{llll}\hfill {E}_{od}\left(\phantom{\rule{1em}{0ex}}incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{1em}{0ex}}=& \phantom{\rule{1em}{0ex}}{E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& \phantom{\rule{1em}{0ex}}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\frac{\phantom{\rule{1em}{0ex}}1}{cos{\mathit{\theta}}_{i}}\sqrt{1-\frac{{sin}^{2}{\mathit{\theta}}_{i}}{{n}^{2}}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(5)}\end{array}$$

Equation (5) is conservation of energy expressed in terms of scalar irradiance for a collimated incident beam.

To illustrate Equations (3) and (5) with HydroLight output, a run was made with the sun at a zenith angle of ${\mathit{\theta}}_{i}=60\phantom{\rule{1em}{0ex}}deg$ in a black sky, in order to obtain a collimated incident beam. The incident ${E}_{d}$ onto the surface was set to 1. The sea surface was level (zero wind speed), and the water index of refraction was $n=1.34$. The scattering coeﬃcient within in the water was set to zero, so that upwelling light from the water column would not confuse the illustration of energy conservation across the sea surface itself. The absorption coeﬃcient was set to $a=0.1\phantom{\rule{1em}{0ex}}{m}^{-1}$. (These runs used HydroLight IOP model ”abconst” and the ”idealized sky model”).

The computed irradiances just above and below the surface are shown in Table 2. The upwelling irradiances in the water are zero because there was no scattering by the water and the bottom was inﬁnitely deep. These computed irradiances satisfy Equations (3) and (5), showing that HydroLight conserves energy across the surface.

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.12690 | 2.0077 | 0.063209 | 1.0000 |

in water at depth 0 | 0.00000 | 1.2277 | 0.0000 | 0.93693 |

A run was then done with a heavily overcast sky to generate a diﬀuse sky radiance pattern with the sun’s location not discernable. The wind speed was $10\phantom{\rule{1em}{0ex}}m\phantom{\rule{1em}{0ex}}{s}^{-1}$ to generate a non-level sea surface. The resulting irradiances are shown in Table 3.

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.18863 | 1.7764 | 0.050456 | 1.0000 |

in water at depth 0 | 0.00000 | 1.1430 | 0.0000 | 0.94954 |

As before, these plane irradiances satisfy Eq. (3). However, the scalar irradiances do not satisfy Eq. (5). The reason is that Eq. (5) must be applied for each particular incident direction ${\mathit{\theta}}_{i}$ of the diﬀuse sky radiance, in order to correctly compute the angle factor seen in Eq. (5). If this is done one direction at a time (with the sun in a black sky), each individual incident direction satisﬁes Eq. (5) just as in the previous example. Thus energy in conserved across the surface ”beam by beam,” hence in toto, even if it is not obvious from the total scalar irradiances.

A further complication occurs when the water body has scattering. To illustrate this, the scattering coeﬃcient was set to $b=0.4\phantom{\rule{1em}{0ex}}{m}^{-1}$ and an ”average-particle” phase function was used in the previous simulation. The computed irradiance were then

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.21264 | 1.7764 | 0.062372 | 1.0000 |

in water at depth 0 | 0.07922 | 1.2046 | 0.031115 | 0.96873 |

Now even the plane irradiance appears to violate conservation of energy, Eq. (3). The reason is that ${E}_{ou}$ and ${E}_{u}$ in air contain both surface-reﬂected radiance (as before) and radiance transmitted upward through the sea surface. Likewise, ${E}_{od}$ and ${E}_{d}$ in the water contain both downward transmitted irradiance (as before) as well as upwelling radiance that has been reﬂected back downward by the surface. This makes it appear at ﬁrst glance that the irradiances do not add up correctly, even though energy actually is conserved across the air-water surface. The eﬀects of both the surface itself and the water body can be accounted for simply by tallying the plane irradiances incident onto and leaving the surface from above and below:

$$\begin{array}{llll}\hfill & {E}_{d}\left(\phantom{\rule{1em}{0ex}}downwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)-{E}_{u}\left(\phantom{\rule{1em}{0ex}}upwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & {E}_{d}\left(\phantom{\rule{1em}{0ex}}downwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)-{E}_{u}\left(\phantom{\rule{1em}{0ex}}upwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(6)}\end{array}$$

Equation (6) expresses conservation of energy for the surface plus water body, for any conditions of incident lighting and water IOPs. The irradiances seen in Table 4 satisfy this equation.

##### Conservation in terms of radiance

Another complication occurs when expressing conservation of energy across the surface in terms of radiance. Now it is necessary to account for changes in both the beam cross sectional area and in solid angle when going across the surface. When interpreting HydroLight output, it must be remembered that the computed radiances are averages over ”quads” or solid angles of ﬁnite size, as discussed on the HydroLight page. Such quads are represented by the green areas in Fig. 2.

Remembering that radiance is power per unit area per unit solid angle, Eq. (2) can be written as

$${L}_{i}{A}_{i}\Delta {\Omega}_{i}={L}_{r}{A}_{r}\Delta {\Omega}_{r}+{L}_{t}{A}_{t}\Delta {\Omega}_{t}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

or

$${L}_{i}={L}_{r}+{L}_{t}\phantom{\rule{0.3em}{0ex}}\frac{{A}_{t}\phantom{\rule{0.3em}{0ex}}\Delta {\Omega}_{t}}{{A}_{i},\Delta {\Omega}_{i}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

after noting that ${A}_{r}={A}_{i}$ and $\Delta {\Omega}_{r}=\Delta {\Omega}_{i}$ by the law of reﬂection. As noted before, ${A}_{t}\u2215{A}_{i}=cos{\mathit{\theta}}_{t}\u2215cos{\mathit{\theta}}_{i}$. The solid angles are given by $\Delta \Omega =\Delta \left(cos\mathit{\theta}\right)\Delta \varphi $, where the angle ranges are determined by the $\mathit{\theta}$ and $\varphi $ boundaries of the quads. Thus the last equation becomes

since $\Delta {\varphi}_{t}=\Delta {\varphi}_{i}$ for the HydroLight quad layout. These angle factors must be evaluated for the particular quads associated with the incident and transmitted beams.

Consider again the previous case of the sun in a black sky at ${\mathit{\theta}}_{i}=60\phantom{\rule{1em}{0ex}}deg$, a level water surface, and no scattering. The incident and reﬂected HydroLight quads are centered at 60 deg and extend from 55 to 65 deg. Snell’s law gives the transmitted direction as ${\mathit{\theta}}_{t}={sin}^{-1}\left[\left(sin{\mathit{\theta}}_{i}\right)\u2215n\right]=40.26\phantom{\rule{1em}{0ex}}deg$. Thus the transmitted radiance is contained in the HydroLight quad centered at 40 deg and extending from 35 to 45 deg. For the transmitted quad, for example, $\Delta \left(cos{\mathit{\theta}}_{t}\right)=cos35-cos45=0.11205$. Evaluating the angle factors in Eq. (7) for these quads gives

The corresponding computed quad-averaged radiances for this run were ${L}_{i}=50.7996$, ${L}_{r}=3.21094$, and ${L}_{t}=41.8544\phantom{\rule{1em}{0ex}}W\phantom{\rule{1em}{0ex}}{m}^{-2}\phantom{\rule{1em}{0ex}}s{r}^{-1}\phantom{\rule{1em}{0ex}}n{m}^{-1}$. These radiances satisfy Eq. (8). Thus energy expressed as radiance is conserved beam-by-beam across the water surface.

We have seen that conservation of energy across the sea surface can be expressed in various ways in terms of radiance and irradiance. However, those expressions are often written for collimated incident beams and consider only the surface itself (Eqs. 3, 5, and 7). It is therefore not always obvious that energy is conserved for the entire surface plus water column when looking at data or simulations for realistic conditions of sky radiances and water IOPs. In general, Eq. (6) should be used to check conservation of energy across the sea surface with measured or modeled plane irradiances.

When making radiance or irradiance measurements in the ocean, it is sometimes hard to verify that energy is being conserved. This is equally true when looking at the output from a numerical model such as HydroLight. Confusion most often arises when comparing measurements made just above and below the sea surface and observing that the measured (or predicted) radiances do not seem to ”add up” correctly. To resolve such confusion, ﬁrst note that it’s the Law of Conservation of Energy, not the law of conservation of radiance or irradiance. We must therefore express the law of conservation of energy in terms of radiance or irradiance as appropriate to the measurements being made.

#### Within the Water

Within the water column, if there are no internal sources such as bioluminescence or inelastic scatter contributing to the wavelength of interest, conservation of energy is expressed in terms of irradiances by Gershun’s law:

Here $a$ is the absorption coeﬃcient, ${E}_{o}$ is the scalar irradiance, and ${E}_{d}$ and ${E}_{u}$ are the downwelling and upwelling plane irradiances. This equation can be illustrated using numerical results generated by HydroLight. A run was made with the sun placed at a 60 deg zenith angle in a clear sky; the wind speed was $5\phantom{\rule{1em}{0ex}}m\phantom{\rule{1em}{0ex}}{s}^{-1}$, and the water was inﬁnitely deep. The absorption coeﬃcient was set to $a=0.1\phantom{\rule{1em}{0ex}}{m}^{-1}$, the scattering coeﬃcient to $b=0.4\phantom{\rule{1em}{0ex}}{m}^{-1}$, and an ”average-particle” phase function was chosen. (The run used HydroLight’s ”abconst” IOP model for these inputs.) No internal sources were included in the run. Table 1 shows the computed irradiances at depths of 5.00 and 5.01 m. Closely spaced depths are used to illustrate Eq. (1) so that the depth derivative can be accurately approximated as a ﬁnite diﬀerence and the scalar irradiance can be taken as the average of the values at 5.00 and 5.01 m.

Using these values in the ﬁnite diﬀerence approximation of Eq. (1) gives

$$\begin{array}{llll}\hfill a=& -\phantom{\rule{0.3em}{0ex}}\frac{1}{0.5\left(0.51833+0.51749\right)}\phantom{\rule{0.3em}{0ex}}\frac{\left(0.32559-0.013129\right)-\left(0.32613-0.013151\right)}{5.01-5.00}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill \phantom{\rule{1em}{0ex}}=& \phantom{\rule{1em}{0ex}}0.1000\phantom{\rule{1em}{0ex}}{m}^{-1}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

The absorption coeﬃcient as computed from the irradiances via Gershun’s law is the same as the one used as input to the run, so HydroLight is conserving energy in its in-water solution of the radiative transfer equation. (HydroLight also conserves energy if internal sources are present. Demonstration of that requires adding an appropriate source term to Gershun’s equation; see Light and Water (1994) Eq. 5.34.) Such accuracy cannot be expected with real data because of imperfect instrument calibration and other sources of error such as Raman scatter (an internal source not accounted for in Eq. 1), which will always be present to some degree during daylight measurements. Nevertheless, Gershun’s law can be a useful way to retrieve absorption coeﬃcients from in-water radiance or irradiance measurements, as seen in papers by Voss (1989) and Maﬃone et al. (1993).

#### Across the Air-Water Surface

The situation is more complicated when considering energy conservation across the air-water surface. Consider a collimated beam of light incident onto a level sea surface. Part will be reﬂected and part will be transmitted into the water. Figure 1 shows this geometry and various quantities needed to express conservation of energy at the surface.

Conservation of energy (or equivalently, power) across the surface requires that

$${P}_{i}={P}_{r}+{P}_{t}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ | (2) |

where ${P}_{i}$ is the power in an incident collimated beam of cross section ${A}_{i}$, ${P}_{r}$ is the power in the reﬂected beam of cross section ${A}_{r}$, and ${P}_{t}$ is the power in transmitted beam of cross section ${A}_{t}$. These beam cross sectional areas are shown in green in Fig. 1. The beam of incident power illuminates a horizontal area ${A}_{s}$ of the level sea surface (shown in blue).

##### Conservation in terms of plane irradiances

Conservation of energy can be expressed in terms of incident (downwelling), reﬂected (upwelling), and transmitted (downwelling) plane irradiances using the horizontal projections of the beam cross section areas. Recalling that plane irradiance is power per unit horizontal area, Eq. (2) becomes

$$\begin{array}{llll}\hfill {A}_{p}{E}_{d}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {A}_{p}{E}_{u}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{p}{E}_{d}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

where the horizontally projected areas ${A}_{p}$ corresponding to each beam are shown in yellow in Fig. 1. Since each beam has the same horizontally projected area, this equation reduces to just

$$\begin{array}{llll}\hfill {E}_{d}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {E}_{u}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {E}_{d}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(3)}\end{array}$$

Equation (3) is conservation of energy across the surface expressed in terms of plane irradiances for a collimated incident beam.

##### Conservation in terms of scalar irradiances

To express Eq. (2) in terms of the scalar irradiances, ﬁrst recall how ${E}_{od}$ is obtained from the radiance (Eq 5 of the Geometrical Radiometry page):

$${E}_{od}\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{\iint}_{2{\pi}_{d}}L\left(\widehat{{\xi}^{\prime}}\right)d\Omega \left(\widehat{{\xi}^{\prime}}\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.$$ |

For a collimated incident beam in direction $\widehat{{\xi}_{o}}$, we can write the radiance in terms of the beam irradiance and a Dirac delta function to get

$${E}_{od}\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{\iint}_{2{\pi}_{d}}{E}_{\perp}\delta \left(\widehat{{\xi}^{\prime}}-\widehat{{\xi}_{o}}\right)d\Omega \left(\widehat{{\xi}^{\prime}}\right)\phantom{\rule{1em}{0ex}}=\phantom{\rule{1em}{0ex}}{E}_{\perp}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

where ${E}_{\perp}$ is the incident irradiance on the beam cross sectional area; ${E}_{\perp}={P}_{i}\u2215{A}_{i}$. Similar equations hold for the reﬂected and transmitted scalar irradiances. These beam cross sectional areas are shown in green in Fig. 1. Thus Eq. (2) becomes

$$\begin{array}{llll}\hfill {A}_{i}{E}_{od}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {A}_{r}{E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{t}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$

or

$$\begin{array}{llll}\hfill {E}_{od}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{t}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\frac{{A}_{t}}{{A}_{i}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill \text{(4)}\end{array}$$after noting that ${A}_{r}={A}_{i}$ by the law of reﬂection. By geometry, the ratio of areas ${A}_{t}\u2215{A}_{i}$ is just $cos{\mathit{\theta}}_{t}\u2215cos{\mathit{\theta}}_{i}$. ${\mathit{\theta}}_{t}$ can be written in terms of ${\mathit{\theta}}_{i}$ using Snell’s law, $sin{\mathit{\theta}}_{i}=nsin{\mathit{\theta}}_{t}$, where $n$ is the real index of refraction of the water. Equation (4) then becomes

$$\begin{array}{llll}\hfill {E}_{od}\left(\phantom{\rule{1em}{0ex}}incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{1em}{0ex}}=& \phantom{\rule{1em}{0ex}}{E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& \phantom{\rule{1em}{0ex}}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\frac{\phantom{\rule{1em}{0ex}}1}{cos{\mathit{\theta}}_{i}}\sqrt{1-\frac{{sin}^{2}{\mathit{\theta}}_{i}}{{n}^{2}}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(5)}\end{array}$$

Equation (5) is conservation of energy expressed in terms of scalar irradiance for a collimated incident beam.

To illustrate Equations (3) and (5) with HydroLight output, a run was made with the sun at a zenith angle of ${\mathit{\theta}}_{i}=60\phantom{\rule{1em}{0ex}}deg$ in a black sky, in order to obtain a collimated incident beam. The incident ${E}_{d}$ onto the surface was set to 1. The sea surface was level (zero wind speed), and the water index of refraction was $n=1.34$. The scattering coeﬃcient within in the water was set to zero, so that upwelling light from the water column would not confuse the illustration of energy conservation across the sea surface itself. The absorption coeﬃcient was set to $a=0.1\phantom{\rule{1em}{0ex}}{m}^{-1}$. (These runs used HydroLight IOP model ”abconst” and the ”idealized sky model”).

The computed irradiances just above and below the surface are shown in Table 2. The upwelling irradiances in the water are zero because there was no scattering by the water and the bottom was inﬁnitely deep. These computed irradiances satisfy Equations (3) and (5), showing that HydroLight conserves energy across the surface.

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.12690 | 2.0077 | 0.063209 | 1.0000 |

in water at depth 0 | 0.00000 | 1.2277 | 0.0000 | 0.93693 |

A run was then done with a heavily overcast sky to generate a diﬀuse sky radiance pattern with the sun’s location not discernable. The wind speed was $10\phantom{\rule{1em}{0ex}}m\phantom{\rule{1em}{0ex}}{s}^{-1}$ to generate a non-level sea surface. The resulting irradiances are shown in Table 3.

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.18863 | 1.7764 | 0.050456 | 1.0000 |

in water at depth 0 | 0.00000 | 1.1430 | 0.0000 | 0.94954 |

As before, these plane irradiances satisfy Eq. (3). However, the scalar irradiances do not satisfy Eq. (5). The reason is that Eq. (5) must be applied for each particular incident direction ${\mathit{\theta}}_{i}$ of the diﬀuse sky radiance, in order to correctly compute the angle factor seen in Eq. (5). If this is done one direction at a time (with the sun in a black sky), each individual incident direction satisﬁes Eq. (5) just as in the previous example. Thus energy in conserved across the surface ”beam by beam,” hence in toto, even if it is not obvious from the total scalar irradiances.

A further complication occurs when the water body has scattering. To illustrate this, the scattering coeﬃcient was set to $b=0.4\phantom{\rule{1em}{0ex}}{m}^{-1}$ and an ”average-particle” phase function was used in the previous simulation. The computed irradiance were then

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.21264 | 1.7764 | 0.062372 | 1.0000 |

in water at depth 0 | 0.07922 | 1.2046 | 0.031115 | 0.96873 |

Now even the plane irradiance appears to violate conservation of energy, Eq. (3). The reason is that ${E}_{ou}$ and ${E}_{u}$ in air contain both surface-reﬂected radiance (as before) and radiance transmitted upward through the sea surface. Likewise, ${E}_{od}$ and ${E}_{d}$ in the water contain both downward transmitted irradiance (as before) as well as upwelling radiance that has been reﬂected back downward by the surface. This makes it appear at ﬁrst glance that the irradiances do not add up correctly, even though energy actually is conserved across the air-water surface. The eﬀects of both the surface itself and the water body can be accounted for simply by tallying the plane irradiances incident onto and leaving the surface from above and below:

$$\begin{array}{llll}\hfill & {E}_{d}\left(\phantom{\rule{1em}{0ex}}downwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)-{E}_{u}\left(\phantom{\rule{1em}{0ex}}upwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & {E}_{d}\left(\phantom{\rule{1em}{0ex}}downwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)-{E}_{u}\left(\phantom{\rule{1em}{0ex}}upwelling,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}.\phantom{\rule{2em}{0ex}}& \hfill \text{(6)}\end{array}$$

Equation (6) expresses conservation of energy for the surface plus water body, for any conditions of incident lighting and water IOPs. The irradiances seen in Table 4 satisfy this equation.

##### Conservation in terms of radiance

Another complication occurs when expressing conservation of energy across the surface in terms of radiance. Now it is necessary to account for changes in both the beam cross sectional area and in solid angle when going across the surface. When interpreting HydroLight output, it must be remembered that the computed radiances are averages over ”quads” or solid angles of ﬁnite size, as discussed on the HydroLight page. Such quads are represented by the green areas in Fig. 2.

Remembering that radiance is power per unit area per unit solid angle, Eq. (2) can be written as

$${L}_{i}{A}_{i}\Delta {\Omega}_{i}={L}_{r}{A}_{r}\Delta {\Omega}_{r}+{L}_{t}{A}_{t}\Delta {\Omega}_{t}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

or

$${L}_{i}={L}_{r}+{L}_{t}\phantom{\rule{0.3em}{0ex}}\frac{{A}_{t}\phantom{\rule{0.3em}{0ex}}\Delta {\Omega}_{t}}{{A}_{i},\Delta {\Omega}_{i}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ |

after noting that ${A}_{r}={A}_{i}$ and $\Delta {\Omega}_{r}=\Delta {\Omega}_{i}$ by the law of reﬂection. As noted before, ${A}_{t}\u2215{A}_{i}=cos{\mathit{\theta}}_{t}\u2215cos{\mathit{\theta}}_{i}$. The solid angles are given by $\Delta \Omega =\Delta \left(cos\mathit{\theta}\right)\Delta \varphi $, where the angle ranges are determined by the $\mathit{\theta}$ and $\varphi $ boundaries of the quads. Thus the last equation becomes

since $\Delta {\varphi}_{t}=\Delta {\varphi}_{i}$ for the HydroLight quad layout. These angle factors must be evaluated for the particular quads associated with the incident and transmitted beams.

Consider again the previous case of the sun in a black sky at ${\mathit{\theta}}_{i}=60\phantom{\rule{1em}{0ex}}deg$, a level water surface, and no scattering. The incident and reﬂected HydroLight quads are centered at 60 deg and extend from 55 to 65 deg. Snell’s law gives the transmitted direction as ${\mathit{\theta}}_{t}={sin}^{-1}\left[\left(sin{\mathit{\theta}}_{i}\right)\u2215n\right]=40.26\phantom{\rule{1em}{0ex}}deg$. Thus the transmitted radiance is contained in the HydroLight quad centered at 40 deg and extending from 35 to 45 deg. For the transmitted quad, for example, $\Delta \left(cos{\mathit{\theta}}_{t}\right)=cos35-cos45=0.11205$. Evaluating the angle factors in Eq. (7) for these quads gives

The corresponding computed quad-averaged radiances for this run were ${L}_{i}=50.7996$, ${L}_{r}=3.21094$, and ${L}_{t}=41.8544\phantom{\rule{1em}{0ex}}W\phantom{\rule{1em}{0ex}}{m}^{-2}\phantom{\rule{1em}{0ex}}s{r}^{-1}\phantom{\rule{1em}{0ex}}n{m}^{-1}$. These radiances satisfy Eq. (8). Thus energy expressed as radiance is conserved beam-by-beam across the water surface.

We have seen that conservation of energy across the sea surface can be expressed in various ways in terms of radiance and irradiance. However, those expressions are often written for collimated incident beams and consider only the surface itself (Eqs. 3, 5, and 7). It is therefore not always obvious that energy is conserved for the entire surface plus water column when looking at data or simulations for realistic conditions of sky radiances and water IOPs. In general, Eq. (6) should be used to check conservation of energy across the sea surface with measured or modeled plane irradiances.

#### Within the Water

Using these values in the ﬁnite diﬀerence approximation of Eq. (1) gives

#### Across the Air-Water Surface

Conservation of energy (or equivalently, power) across the surface requires that

$${P}_{i}={P}_{r}+{P}_{t}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ | (2) |

##### Conservation in terms of plane irradiances

##### Conservation in terms of scalar irradiances

or

$$\begin{array}{llll}\hfill {E}_{od}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{t}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\frac{{A}_{t}}{{A}_{i}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill \text{(4)}\end{array}$$

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.12690 | 2.0077 | 0.063209 | 1.0000 |

in water at depth 0 | 0.00000 | 1.2277 | 0.0000 | 0.93693 |

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.18863 | 1.7764 | 0.050456 | 1.0000 |

in water at depth 0 | 0.00000 | 1.1430 | 0.0000 | 0.94954 |

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.21264 | 1.7764 | 0.062372 | 1.0000 |

in water at depth 0 | 0.07922 | 1.2046 | 0.031115 | 0.96873 |

##### Conservation in terms of radiance

Remembering that radiance is power per unit area per unit solid angle, Eq. (2) can be written as

or

#### Within the Water

Using these values in the ﬁnite diﬀerence approximation of Eq. (1) gives

#### Across the Air-Water Surface

Conservation of energy (or equivalently, power) across the surface requires that

$${P}_{i}={P}_{r}+{P}_{t}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},$$ | (2) |

##### Conservation in terms of plane irradiances

##### Conservation in terms of scalar irradiances

or

$$\begin{array}{llll}\hfill {E}_{od}\left(incident,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)=& {E}_{ou}\left(reflected,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}air\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill +& {A}_{t}{E}_{od}\left(transmitted,\phantom{\rule{1em}{0ex}}in\phantom{\rule{1em}{0ex}}water\right)\frac{{A}_{t}}{{A}_{i}}\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}},\phantom{\rule{2em}{0ex}}& \hfill \text{(4)}\end{array}$$

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.12690 | 2.0077 | 0.063209 | 1.0000 |

in water at depth 0 | 0.00000 | 1.2277 | 0.0000 | 0.93693 |

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.18863 | 1.7764 | 0.050456 | 1.0000 |

in water at depth 0 | 0.00000 | 1.1430 | 0.0000 | 0.94954 |

location | ${E}_{ou}$ | ${E}_{od}$ | ${E}_{u}$ | ${E}_{d}$ |

in air | 0.21264 | 1.7764 | 0.062372 | 1.0000 |

in water at depth 0 | 0.07922 | 1.2046 | 0.031115 | 0.96873 |

##### Conservation in terms of radiance

Remembering that radiance is power per unit area per unit solid angle, Eq. (2) can be written as

or