You are viewing a version of the Ocean Optics Web Book site specially-formatted for printing.
We now further examine the diffuse attenuation functions, which are often called "K functions." Under typical oceanic conditions, for which the incident lighting is provided by the sun and sky, the various radiances and irradiances all decrease approximately exponentially with depth in homogeneous water, when far enough below the surface (and far enough above the bottom, in optically shallow water) to be free of boundary effects. It is therefore convenient to write the depth dependence of , for example, as
where is the diffuse attenuation function for spectral downwelling plane irradiance. Solving for gives the equation previously used for the definition of :
Similar equations can be written for other radiometric variables and their corresponding K functions.
As explained in the previous discussion of AOPs, in order to be useful for relating light measurements to water properties, the K functions should depend strongly on water IOPs but only weakly on external environmental conditions like sun location, sky condition, or surface waves. To illustrate these dependencies, various K functions were numerically computed using the HydroLight radiative transfer numerical model. In most situations it is preferable to work with real data. However, use of this model gives us the ability to simulate different environmental conditions at will and to control things that cannot be controlled in nature, such as turning chlorophyll fluorescence on or off. This can be very useful for understanding the interdependence of various quantities.
HydroLight simulations were first performed at one wavelength for homogeneous idealized water bodies dominated by either scattering or absorption. For the "highly scattering" water, the absorption coefficient was set to , and the scattering coefficient was , so that the albedo of single scattering was . These values correspond roughly to what might be found in Case 1 water with a chlorophyll concentration of at blue or green wavelengths. An "average particle" scattering phase function was used, the sun was placed in a clear sky, and the water was infinitely deep. Note that since , the optical depth is numerically equal to the geometric depth z in meters.
Figure 1 shows various K functions as a function of depth for the highly scattering water when the sun was placed at a zenith angle of 40 degrees and the surface was level (wind speed of ). As is conventional in radiative transfer theory, the are for radiance propagating in the direction. HydroLight measures depth positive downward from the mean sea surface, and polar angle is measured from the or nadir direction. Thus refers to light heading straight down into the water (corresponding to ), refers to light heading straight up (corresponding to , and refers to light traveling horizontally through the water. Azimuthal angle refers to light heading towards the sun, which was placed at ; thus refers to light heading away from the sun. In this simulation, corresponds to looking into the sun's refracted beam underwater, which is light heading downward and away from the sun.
There are several important features to note in Fig. 1:
Figure 2 shows the K functions corresponding to the same conditions as Fig. 1, except that the wind speed was . We see that there is very little difference between Figs. 1 and 2. Thus, as hoped, the K functions are almost unaffected by the surface waves.
Figure 3 shows the K functions for a level surface and the sun at the zenith, rather than at 40 deg. Again, the irradiance K functions are almost unchanged. However, the radiance function now corresponds to looking straight upward into the sun's direct beam. Thus now looks very similar to in the previous figures. Similarly, now looks much like in the previous figures. This is because moving the sun from 40 deg in air (28 deg in water) to the zenith gives almost the same scattering angle relation to the sun's direct beam as had in the previous figures.
Figure 4 shows K functions for "highly absorbing" water: the absorption coefficient was , and the scattering coefficient was , so that the albedo of single scattering was . These values correspond roughly to what might be found at red wavelengths, where absorption by the water itself begins to dominate the IOPs. Other conditions were the same as for Fig. 1.
Comparing Figs. 1 and 4, we note that
Comparing Figs. 1 and 4 also shows that the K's have changed greatly because of the change in IOPs, which is what is desired in any AOP. For the highly absorbing water, .
The distinction between beam and diffuse attenuation is important. The beam attenuation coefficient c is defined in terms of the radiant power lost from a collimated beam of photons. The downwelling diffuse attenuation function , for example, is defined in terms of the decrease with depth of the ambient downwelling irradiance , which comprises photons heading in all downward directions (a diffuse, or uncollimated, radiance distribution). In the above simulations, at all depths, but all K functions except the near-surface in the high absorption case are less than c. Radiative transfer theory shows (e.g., Light and Water (1994) Eq. 5.71), for example, that in general , where is the mean cosine of the downwelling radiance. These inequalities are seen to hold true in the above simulations.
HydroLight was next used to simulate a homogeneous Case 1 water body with a chlorophyll concentration of . As for Figs. 1 and 4, the sun was at 40 deg and the surface was level. Figure 5 plots several quantities as a function of wavelength at 10 m depth for this simulation. We see that between 300 and about 600 nm, the various functions are very similar and proportional to the total (including water) absorption coefficient . However, beyond 600 nm the K functions differ from each other, and they are all much different from a. The reason for this behavior is inelastic scattering.
At the near-UV to blue to green wavelengths below 600 nm, most of the radiance at 10 m depth (for these IOPs) comes from sunlight being transmitted through the upper 10 m of the water column. Above about 600 nm, absorption by the water itself has removed most of the sunlight. For example, at 700 nm where , we expect roughly of the surface light to reach 10 m. However, Raman scatter and CDOM fluorescence inelastically scatter light from shorter wavelengths, where sunlight is present, into the red wavelengths, and thus create additional red light at 10 m. Thus, beyond 600 nm, the various radiances and irradiances no longer decrease with depth in the simple exponential fashion expected. Note that tracks a longer than do and . This is because continues to collect whatever downwelling sunlight remains, and thus the inelastic contribution to is not noticeable until the chlorophyll fluorescence contribution begins near 670 nm. and on the other hand have only a small amount of backscattered sunlight, so that the inelastic contribution becomes significant sooner, at around 600 nm.
It is easy to verify that the peculiar behavior of the K functions beyond 600 nm is due to inelastic scatter. The HydroLight run was repeated with Raman scatter and CDOM and chlorophyll fluorescence "turned off." Figure 6 shows the results. Now, the K functions all track the absorption nicely at all wavelengths.
As seen in these figures, radiative transfer theory shows that K functions are very "absorption like," meaning that the K functions are strongly correlated with the total absorption coefficient when inelastic scatter effects are negligible. For , the approximate relation gives close agreement between the exact (computed by HydroLight) and the value estimated from the absorption coefficient and the downwelling mean cosine of the radiance distribution, which was also obtained from the HydroLight simulation.
These few simulations are enough to establish the salient features of diffuse attenuation functions. Their use has a venerable history in optical oceanography. Smith and Baker (1978) listed some of their virtues:
Gordon's normalization of Kd
Gordon (1989) developed a simple way to normalize measured values. His normalization for all practical purposes removes the effects of the sea state and incident sky radiance distribution from , so that the normalized can be regarded as an IOP. The theory behind the normalization is given in his paper; the mechanics of the normalizing process are as follows.
Let be the irradiance incident onto the sea surface due to the sun's direct beam, and let be the irradiance due to diffuse sky radiance. Then the fraction f of the direct sunlight in the incident irradiance that is transmitted through the surface into the water is
Here t(sun) and t(sky) are respectively the fractions of the direct beam and of the diffuse irradiance transmitted through the surface; these quantities can be computed using methods described in Light and Water (1994) Chapter 4 [where they are denoted by t(a,w)]. However, if the solar zenith angle in air, , is less than 45 degrees , then . If the sky radiance distribution is nearly uniform (as it is for a clear sky), then . In this case, we can accurately estimate f from measurements made just above the sea surface:
The sun and sky irradiances are easily obtained from an instrument on the deck of a ship. When both direct and diffuse light fall onto the instrument, it records . When the direct solar beam is blocked, the instrument records . (Advanced technology is not required here: just hold your hat so that its shadow falls on the instrument.)
Next compute the nadir angle of the transmitted solar beam in water, , using Snell's law:
Finally, compute the quantity
This value of is valid for flat or rough sea surfaces as long as . For larger values of , or for an overcast sky, a correction must be applied to to account for surface wave effects on the transmitted light; the correction factors are given in Gordon (1989, his Fig. 6). Gordon's normalization then consists simply of dividing the measured by :
Physically, is a function (essentially ) that reduces values to the values that would be measured if the sun were at the zenith, if the sea surface were level, and if the sky were black (i.e., if there were no atmosphere). The zenith-sun, level-surface, black-sky case is the only physical situation for which . In other words, normalization by removes the influence of incident lighting and sea state on . The same normalization can be applied to depth-averaged values of .
We recommend that experimentalists routinely make the simple measurements necessary to determine , because normalization of enhances its value in the recovery of IOP's from irradiance measurements.