It looks like you're using Internet Explorer 6: Features on this site are not supported by that browser version. Please upgrade to the latest version of Internet Explorer.

You are viewing a version of the Ocean Optics Web Book site specially-formatted for printing.

# level 2

## The Vector Radiative Transfer Equation

 Page updated: Jan 3, 2014 Principal author: Curtis Mobley

When describing polarized light at the phenomenological level of radiative transfer theory, the radiance is specified by the 4-component Stokes column vector , where the superscript T indicates the transpose (on this page, underlined bold font indicates 4 component vectors or 4 x 4 matrices). Vector component I is the total radiance, as is measured by a detector that is not sensitive to polarization. (L is usually used for radiance in oceanography, but I is customary in polarization studies involving Stokes vectors.) Q and U describe the state of linear polarization in a coordinate system chosen for the problem at hand, and V describes the state of circular polarization. Unpolarized radiance has the Stokes vector , with I being the radiance in .

In oceanic optics, it is convenient to pick a coordinate system with the and axes lying in the mean sea surface and the axis normal to the mean sea surface (bold font not underlined indicates unit vectors specifying directions in Euclidean 3-space). In HydroLight, for example, the axis points in the downwind direction and points downward into the water. Figure 1 shows such a coordinate system. The red arrow in the figure indicates the direction of propagation of an unscattered beam of light; the beam direction is specified by polar and azimuthal angles in the coordinate system. The shaded plane defined by and is called the meridian plane. is a unit vector that is normal to the meridian plane and thus parallel to the mean sea surface, i.e., "horizontal." Direction lies in the meridian place and thus is "vertical" to the mean sea surface when . Stokes vectors are usually specified with reference to the meridian plane and the "horizontal-vertical" coordinate system of Fig. 1. If the light is linearly polarized with the electric field vector E vibrating in the meridian plane as illustrated by the blue arrows in Fig. 1, then the light is said to be "vertically" polarized. The corresponding Stokes vector is , if the magnitude of the total radiance is .

Because the electric field describing a propagating wave at the level of Maxwell's equations is a vector, and the corresponding Stokes representation of the radiance is a 4-component column vector, care must be taken when referring these vectors to different coordinate systems when light is scattered. When light is scattered from incident direction to a final direction , the meridian planes are different. Thus the initial and final Stokes vectors are specified in different coordinate systems. In addition, the 4x4 Mueller matrix M, which describes scattering in the Stokes-vector representation, is specified in the scattering plane defined the by incident and final directions and , respectively.

Figure 2 shows the coordinate systems needed to define scattering of polarized radiance. The coordinate system is the same as in Fig. 1. The coordinate system is similarly defined for the meridian plane of the light scattered into direction . The scattering plane is shown shaded. The scattering angle can be computed from

 (1)

The Stokes vector must be specified in a coordinate system defined by the scattering plane before the Mueller matrix can be applied. The coordinate system of the incident light in the scattering plane is defined with being perpendicular (senkrecht, in German) to the scattering plane, and being parallel to the scattering plane, as shown in Fig. 2. The Stokes vector of the scattered light obtained from application of the Mueller matrix to the rotated incident vector is specified in the scattering plane system, but is measured in the meridian plane system. The steps necessary to describe the scattering process are thus

• Specify the Stokes vector of the incident radiance in the coordinate system defined by the incident meridian plane, where is measured.
• Rotate the coordinate system into the coordinate system of the scattering plane, to obtain the components of in the scattering plane, where the Mueller matrix is specified.
• Apply the Mueller matrix to scatter through scattering angle into in the scattering plane system .
• Rotate the system into the final meridian plane , where the scattered Stokes vector is measured.

Define a rotation through a positive angle as being a counterclockwise rotation about the beam direction, when looking into the beam. Thus the initial rotation from into is a rotation through a positive angle , and the final rotation from into is a rotation through a negative angle , as illustrated in Fig. 2. Spherical trigonometry gives these angles as (van de Hulst, 1980), p 499)

 (2) (3)

With this definition of a positive angle rotation, the 4x4 rotation matrix that transforms from the initial meridian plane into the scattering plane is (e.g., Mischenko et al., 2002, Eq. 1.97)

 (4)

[Other rotation angle conventions are common in the literature. Mischenko et al. (2002) define a positive rotation as being clockwise looking along the beam, which is equivalent to the present counterclockwise looking into the beam. However, Bohren and Huffman, (1983, their Eq. 2.83), Lenoble, (1993, her Eq. 14.24), and Kattawar and Adams (1989, their Eq. 8) define a positive rotation as being counterclockwise looking along the beam, or clockwise looking into the beam. Those authors thus derive a Stokes vector rotation matrix that is the transpose of Eq. (4).]

The matrix multiplications involved in scattering into are thus

 (5)

These preliminaries allow us to write (in analogy with the scalar RTE), the 1D, time-independent, vector (or polarized) radiative transfer equation (VRTE) as

 (6)

Here c is the beam attenuation coefficient. is any internal source of polarized or unpolarized radiance in direction , with units of radiance per unit distance, i.e. . The Stokes vectors, , c, and generally depend on depth and wavelength, but those arguments are not shown here for brevity. It should be remembered that (just as for the volume scattering function) each element M(i,j) of a Mueller matrix is also a function of the scattering angle , and that the functional form of these elements is determined by the type of scattering particle. In the integration over all incident directions, the angles are obtained from Eqs. (1-3) for the given .

The M(1,1) element of is the volume scattering function for unpolarized radiance, . Mueller matrices are often presented as reduced Mueller matrices in which each element of has been divided by the M(1,1) element. This amounts to "factoring out" the volume scattering function for unpolarized radiance. Thus, for example, the Mueller matrix for Rayleigh scattering can be written as

 (7)

or

 (8)

Here

is the reduced Mueller matrix for Rayleigh scattering seen in Eq. (7).

is the Rayleigh volume scattering function for unpolarized radiance [with units of ].

is the Rayleigh scattering coefficient [units of ], which has a wavelength dependence.

is the Rayleigh phase function for unpolarized radiance [with units of ]. This phase function satisfies the normalization condition .

The VRTE can thus be written as

 (9)

The product of the rotation and reduced Mueller matrices is often written as

 (10)

where is called the phase matrix in analogy to the phase function of scalar theory. The VRTE is then

 (11)

Note that if and all other elements are 0, then Eq. (11) reduces to the scalar RTE.

Mueller matrices for spherical particles larger than those of the Rayleigh domain (particle size ) can be computed by Mie theory. However, there are almost no actual measurements of Mueller matrices for ocean waters. Figure 3 shows one of the few actual data sets, acquired by Voss and Fry (1984) in Atlantic waters. It should be noted that the reduced Mueller matrix for ocean waters is similar, but not identical, to that for Rayleigh scattering.

The similarities seen in Fig. 3 between the reduced Muller matrices for Rayleigh scattering and for ocean waters suggests that for some (but probably not all) purposes it may be reasonable to use a Mueller matrix that combines a VSF for unpolarized scattering by oceanic particles and the reduced Mueller matrix for Rayleigh scattering, which has a convenient analytic form. With this approximation, in Eq. (11) would be replaced by . Whether or not this is an acceptable simplification for general oceanographic calculations awaits measurements of in-water Mueller matrices for a variety of oceanic particles and subsequent comparison of numerical computations based on those measurements with those based on the simplification.