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

# The Level Sea Surface

This chapter discusses how light is reﬂected by and transmitted through air-water surfaces, and how it is reﬂected by opaque surfaces such as the ocean bottom. The Level 1 discussion begins with the basics of reﬂection and transmission of unpolarized light by a level or ﬂat water surface. Although the ocean is rarely glassy smooth, reﬂection and transmission by rough, wind-blown surfaces are modeled using the equations for a ﬂat surface applied to each small patch of sea surface, which although tilted from the normal to the mean sea surface can be assumed to be locally ﬂat. Other Level 1 material introduces the bidirectional reﬂectance distribution function or BRDF. The BRDF is the fundamental quantity for specifying how an opaque surface reﬂects light.

The Level 2 material ﬁrst considers reﬂection and transmission for polarized light. Then several pages show how wind-blown sea surfaces can be described in terms of energy spectra and, conversely, how random sea surfaces can be generated starting with energy spectra. These techniques are widely used both for generation of sea surfaces for quantitative modeling of reﬂection and transmission, and for computer animation of sea surfaces as used in many movies.

#### Geometric Relations

The wavelength of visible light is much, much less than the millimeter and larger spatial wavelengths of the waves on wind-blown surfaces. Therefore, the laws of geometrical optics and the idealization of a narrow ray of collimated light give a good description of the relevant physical processes.

Figure 1 illustrates a level surface with light incident onto the surface from the air side (panel a), and from the water side (panel b). The real index of refraction of the air is ${n}_{a}$, which is taken to be one. ${n}_{w}$ is the real index of refraction of the water, which is approximately 1.34 at visible wavelengths. $\widehat{n}$ is a unit vector normal to the surface. $\widehat{\xi}i$ is a unit vector in the incident direction; $\widehat{\xi}r$ and $\widehat{\xi}t$ are respectively the directions of the reﬂected and transmitted rays. ${\mathit{\theta}}_{i}={cos}^{-1}\left(|\widehat{\xi}i\cdot \widehat{n}|\right)$ is the acute angle between the incident direction and the normal, and ${\mathit{\theta}}_{r}$ and ${\mathit{\theta}}_{t}$ are the angles of the reﬂected and transmitted rays relative to the normal.

The incident, reﬂected, and refracted directions all lie in the plane deﬁned by $\widehat{\xi}i$ and $\widehat{n}$. The reﬂected angle is always equal to the incident angle: ${\mathit{\theta}}_{r}={\mathit{\theta}}_{i}$, which is known as the Law of Reﬂection. The incident and transmitted angles are related by

$${n}_{1}sin{\mathit{\theta}}_{1}={n}_{2}sin{\mathit{\theta}}_{2}\phantom{\rule{0.3em}{0ex}},$$ | (1) |

where subscripts 1 and 2 refer to any two media. This equation is usually called Snell’s law, although more properly it should be Snel’s law. It was rediscovered in the west by Willebrord Snel van Royen (1580-1626). In the days when European scientists published in Latin, Snel’s name was Latinized to Snellius, which became Snell as Latin was replaced by German and then English as the common language of physical science. However, this law (in a diﬀerent but equivalent form) can be traced back to the treatise On Burning Mirrors and Lenses published by the Persian Abu ibn Sahl in Bagdad in 984.

Figure 2 gives a visual representation of the relations between the various unit vectors, angles, and indices of refraction.

For air-incident light, ${n}_{a}=1$ and Snel’s law reads $sin{\mathit{\theta}}_{i}={n}_{w}sin{\mathit{\theta}}_{t}$. Then the angle of transmission is given by

$${\mathit{\theta}}_{t}={sin}^{-1}\left(\frac{1}{{n}_{w}}sin{\mathit{\theta}}_{i}\right)\phantom{\rule{0.3em}{0ex}}.$$ | (2) |

The relations between the unit vectors are given by the following equations (with ${n}_{a}=1$):

$$\widehat{\xi}r=\widehat{\xi}i-2\left(\widehat{\xi}i\cdot \widehat{n}\right)\widehat{n}\phantom{\rule{0.3em}{0ex}},$$ |

and

$$\widehat{\xi}t=\frac{1}{{n}_{w}}\left(\widehat{\xi}i-c\widehat{n}\right)\phantom{\rule{0.3em}{0ex}},$$ |

where

$$c=\widehat{\xi}i\cdot \widehat{n}+\sqrt{{\left(\widehat{\xi}i\cdot \widehat{n}\right)}^{2}+{n}_{w}^{2}-1}\phantom{\rule{0.3em}{0ex}}.$$ |

For the water-incident case, Snel’s law reads ${n}_{w}sin{\mathit{\theta}}_{i}=sin{\mathit{\theta}}_{t}$, in which case the angle of transmission is given by

$${\mathit{\theta}}_{t}={sin}^{-1}\left({n}_{w}sin{\mathit{\theta}}_{i}\right)\phantom{\rule{0.3em}{0ex}}.$$ | (3) |

The relations between the unit vectors are then

$$\widehat{\xi}r=\widehat{\xi}i-2\left(\widehat{\xi}i\cdot \widehat{n}\right)\widehat{n}\phantom{\rule{0.3em}{0ex}},$$ |

and

$$\widehat{\xi}t={n}_{w}\widehat{\xi}i-c\widehat{n}\phantom{\rule{0.3em}{0ex}},$$ |

where

$$c={n}_{w}\widehat{\xi}i\cdot \widehat{n}-\sqrt{{\left({n}_{w}\widehat{\xi}i\cdot \widehat{n}\right)}^{2}-{n}_{w}^{2}+1}\phantom{\rule{0.3em}{0ex}}.$$ |

If ${n}_{w}sin{\mathit{\theta}}_{i}<1$, Eq. (3) gives a real value for ${\mathit{\theta}}_{t}$ and light is transmitted from the water to the air. However, if ${\mathit{\theta}}_{i}$ is greater than the critical angle for total internal reﬂection

$${\mathit{\theta}}_{c}={sin}^{-1}\left(1\u2215{n}_{w}\right)\phantom{\rule{0.3em}{0ex}},$$ | (4) |

then there is no real solution for the inverse sine. In this case, all light incident onto the water side of the air-water surface is reﬂected back into the water. This is called total internal reﬂection. The dotted line in the right panel of Fig. 1 represents the critical angle. The red unit vectors illustrate the case ${\mathit{\theta}}_{i}<{\mathit{\theta}}_{c}$ with both reﬂected and transmitted light, and the yellow vectors represent the case of ${\mathit{\theta}}_{i}>{\mathit{\theta}}_{c}$ and total internal reﬂection.

It should be noted that air-incident light with a grazing incident angle of ${\mathit{\theta}}_{i}=90\phantom{\rule{1em}{0ex}}deg$ is transmitted into the water at the critial angle: ${\mathit{\theta}}_{t}={\mathit{\theta}}_{c}$. Thus light from the entire sky is transmitted through the surface into a cone of half angle ${\mathit{\theta}}_{c}$, which is known as Snel’s cone.

#### Fresnel’s Equations for Unpolarized Light

The equations of the previous section show the relations between the angles and directions of the incident and the reﬂected and transmitted light. However, they do not show how much light is reﬂected or transmitted. That information is given by Fresnel’s equations.

Consider a collimated beam of unpolarized incident light, which has some irradiance measured on a surface normal to the direction $\widehat{\xi}i$ of propagation. The fraction of this incident irradiance that is reﬂected by the air-water surface is

which holds for ${\mathit{\theta}}_{i}\ne 0$. For normally incident light, ${\mathit{\theta}}_{i}=0$, the reﬂectance is

$${R}_{F}\left({\mathit{\theta}}_{i}=0\right)={\left(\frac{{n}_{w}-1}{{n}_{w}+1}\right)}^{2}\phantom{\rule{0.3em}{0ex}}.$$ | (6) |

Equations (5) and (6) hold for both air- and water-incident light. Given the incident angle ${\mathit{\theta}}_{i}$, the transmitted angle ${\mathit{\theta}}_{t}$ is computed using either Eq. (2) or (3), and then Eq. (5) (or 6) can be evaluated. For water-incident light and ${\mathit{\theta}}_{i}\ge {\mathit{\theta}}_{c}$, ${R}_{F}=1$. Figure 3 shows the Fresnel reﬂectance for the range of water indices of refraction at visible wavelengths.

To be completely general, the Fresnel equations should use the complex index of refraction $m=n+ik$, where $n$ is the real index of refraction seen above and $k\left(\lambda \right)=\lambda a\left(\lambda \right)\u22152\pi $ is the complex part ($a$ is the absorption coeﬃcient). Thus Eq. (6) should be

$${R}_{F}\left({\mathit{\theta}}_{i}=0\right)={\left|\frac{m-1}{m+1}\right|}^{2}=\frac{{\left({n}_{w}-1\right)}^{2}+{k}^{2}}{{\left({n}_{w}+1\right)}^{2}+{k}^{2}}\phantom{\rule{0.3em}{0ex}}.$$ |

However, for water at near-UV to near-IR wavelengths, $k<1{0}^{-6}$ and the diﬀerence is negligible. However $k$ can be of order 0.1 to 1 at some UV and far-IR wavelengths, in which case the complex index of refraction must be used.

Conservation of energy requires that the sum of the reﬂected and transmitted energy equal the incident energy. Thus the fraction of the incident energy that is transmitted is ${T}_{F}=1-{R}_{F}$. It can be confusing to see that energy is conserved when diﬀerent quantities such as plane irradiance, scalar irradiance, or radiance are used to describe the light, or when the incident light is not a single collimated beam. This is discussed in detail on the page on energy conservation.

Reﬂection and transmission are much more complicated when the incident light is polarized. The Fresnel reﬂectance and transmittance equations for polarized light are given in the Level 2 material for this chapter.

#### The ${n}^{2}$ Law for Radiance

Snel’s law yields an important result governing how unpolarized radiance changes when going from one medium to another, e.g., when crossing an air-water surface. Figure 4 shows two beams of radiance, one incident onto an interface and one transmitted. Let ${L}_{1}$ be the incident radiance in medium 1 deﬁned by power ${\Phi}_{1}$ passing through an area $\Delta {A}_{1}$ normal to the direction of photon travel and contained in a solid angle $\Delta {\Omega}_{1}=sin{\mathit{\theta}}_{1}\Delta {\mathit{\theta}}_{1}\Delta {\varphi}_{1}$, where ${\mathit{\theta}}_{1}$ is polar angle measured relative to the normal to the surface and $\Delta {\varphi}_{1}$ is the width of the solid angle in the azimuthal direction. Likewise, ${L}_{2}$ is the transmitted radiance in medium 2 deﬁned by the corresponding quantities as illustrated. The azimuthal angle does not change when crossing the surface, so $\Delta {\Omega}_{2}=sin{\mathit{\theta}}_{2}\Delta {\mathit{\theta}}_{2}\Delta {\varphi}_{1}$. The incident and transmitted power passes through a common area $\Delta A$ at the interface.

The indices of refraction ${n}_{1}$ and ${n}_{2}$ are ﬁxed, but the polar angle $\mathit{\theta}$ changes when crossing the interface. Squaring Eq. (1) and diﬀerentiating gives

$${n}_{1}^{2}sin{\mathit{\theta}}_{1}cos{\mathit{\theta}}_{1}\Delta {\mathit{\theta}}_{1}={n}_{2}^{2}sin{\mathit{\theta}}_{2}cos{\mathit{\theta}}_{2}\Delta {\mathit{\theta}}_{2}\phantom{\rule{0.3em}{0ex}}.$$ |

Multiplying each side of this equation by the common value of $\Delta \varphi $ and rewriting in terms of solid angles gives

$${n}_{1}^{2}cos{\mathit{\theta}}_{1}\Delta {\Omega}_{1}={n}_{2}^{2}cos{\mathit{\theta}}_{2}\Delta {\Omega}_{2}\phantom{\rule{0.3em}{0ex}},$$ |

which is known as Straubel’s invariant.

The radiances are deﬁned by

$${L}_{1}=\frac{\Delta {\Phi}_{1}}{\Delta {A}_{1}\Delta {\Omega}_{1}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}and\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}{L}_{2}=\frac{\Delta {\Phi}_{2}}{\Delta {A}_{2}\Delta {\Omega}_{2}}\phantom{\rule{0.3em}{0ex}}.$$ |

Fresnel’s equation gives the transmitted power as $\Delta {\Phi}_{2}=\left[1-{R}_{F}\left({\mathit{\theta}}_{1}\right)\right]\Delta {\Phi}_{1}={T}_{F}\Delta {\Phi}_{1}$. The areas are related by $\Delta {A}_{1}=\Delta Acos{\mathit{\theta}}_{1}$ and $\Delta {A}_{2}=\Delta Acos{\mathit{\theta}}_{2}$. Thus the ratios of the incident and transmitted radiances can be written as

$$\begin{array}{llll}\hfill \frac{{L}_{2}}{{L}_{1}}=& \phantom{\rule{1em}{0ex}}\frac{\Delta {\Phi}_{2}}{\Delta {\Phi}_{1}}\frac{\Delta {A}_{1}\Delta {\Omega}_{1}}{\Delta {A}_{2}\Delta {\Omega}_{2}}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill =& \phantom{\rule{1em}{0ex}}{T}_{F}\frac{cos{\mathit{\theta}}_{1}\Delta {\Omega}_{1}}{cos{\mathit{\theta}}_{2}\Delta {\Omega}_{2}}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill =& \phantom{\rule{1em}{0ex}}{T}_{F}\frac{{n}_{2}^{2}}{{n}_{1}^{2}}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$or

$$\frac{{L}_{2}}{{n}_{2}^{2}}=\phantom{\rule{1em}{0ex}}{T}_{F}\frac{{L}_{1}}{{n}_{1}^{2}}\phantom{\rule{0.3em}{0ex}}.$$ |

This result is called the n-squared law for radiance. The quantity $L\u2215{n}^{2}$ is sometimes called the reduced radiance or the basic radiance.

Although energy is conserved when crossing a boundary, the radiance changes by a factor proportional to the change in the index of refraction squared. This is a simple consequence of the change in solid angle resulting from the change in $\mathit{\theta}$ when crossing the boundary. Note that for normal incidence and ${n}_{w}=1.34$, ${T}_{F}\approx 0.979$ and the radiance just below a water surface is $0.979{\left(1.34\right)}^{2}\approx 1.76$ times the radiance in the air. Conversely, when going from water to air, the in-water radiance is reduced by a factor of 1.76.

To the extent that losses to absorption and scattering out of the beam can be ignored (sometimes a good approximation for atmospheric transmission, but almost never the case in water), the radiance divided by the square of the index of refraction is constant along any path. This result has even been called the fundamental theorem of radiometry, which is perhaps a bit grandiose given that real beams always lose radiance due to absorption and can lose or gain radiance due to scattering.

Finally, note that the ${n}^{2}$ law applies only to radiance transmission. When tracing photons in a Monte Carlo simulation, from which the radiance can be estimated by appropriate binning of the transmitted photons, no ${n}^{2}$ factor is applied to the energy of the transmitted photons or to the radiance estimated from the detected photons. This is because the ${n}^{2}$ eﬀect is automatically built into the radiance estimate photon by photon as the directions of the individual photons are computed by Snel’s law.