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.

Ocean Optics Web Book

Skip to main content
Ocean Optics Web Book

Radiative Transfer Theory

Solving the Radiative Transfer Equation

Page updated: Oct 31, 2013
Principal author: Curtis Mobley

As seen on the previous page, the RTE is easy to derive. However, it is quite difficult to solve. The 1D, time-independent RTE is a linear two-point boundary value problem. That is to say, there are boundary conditions describing the radiance at the top and bottom of the atmosphere or ocean (the two spatial points), and the propagation of radiance within the medium--between the boundaries--is governed by the linear integro-differential RTE. In the oceanographic setting, the upper boundary condition specifies the sky radiance incident onto the sea surface. The lower boundary condition specifies how the sea bottom reflects the downwelling radiance. (The "sea bottom" can be a physical bottom, or the deepest depth in the water column at which the RTE needs to be solved to obtain the radiances above that depth.)

Exact Analytical Solutions

Exact analytical (i.e., pencil and paper) solutions of the RTE can be obtained only for very simple situations, such as no scattering. There is no function (that anyone has ever found) that gives

$\displaystyle L(z,\theta,\phi,\lambda) = f(a,~\rm {VSF,~sun~angle,~bottom~reflectance,~etc.}),

where $ f$ is some function in which you can "plug in" the absorption coefficient, VSF, and other parameters and get back the radiance. This is true even for very simple situations such as homogeneous water with isotropic scattering. Even the conceptually simple geometry of an isotropically emitting point light source in an infinite homogeneous ocean is unsolved. (A very complicated solution for the the scalar irradiance $ E_o(r)$ around an isotropically emitting point source with isotropic scattering does exist; see Davison and Sykes (1957), eq. 5.25 or Mobley (1996), Section 9.2.) This may seem surprising because other point-source problems, such as the electric filed around a point charge or the gravitational field around a point mass, are often very simple. The difference with optics lies in the complications caused by scattering within the medium surrounding the point source (which do not exist for problems like the gravitational field around a point mass).

If there is no scattering, the RTE of Eq. (9) on the previous page reduces to just

$\displaystyle \frac{d L (r, \theta , \phi , \lambda )}{d r} ~=~- a(r,\lambda ) L(r, \theta , \phi , \lambda ) 
 +~S(r, \theta , \phi ,\lambda )\,\,.$ (1)

This is a linear, first order, ordinary differential equation, which is easily solved (see any text on differential equations, e.g., Rainville (1964), page 36). If the medium is homogeneous, so that the absorption coefficient and source function do not depend on distance $ r$ , the solution of Eq. (1) is (dropping the wavelength and direction arguments for brevity)

$\displaystyle L(r)~=~L(0)e^{-ar} + \frac{S}{a} [ 1 - e^{-ar}]\,\, ,$ (2)

where $ L(0)$ is the initial radiance at distance $ r = 0$ . In source-free water, $ S = 0$ and the solution is a simple exponential decay of the initial radiance with distance. This result was derived in another way in the discussion on measuring IOPs; see Eq. (3) that page.

Note that for great distances, $ r \to \infty$ , the radiance depends only on the source function and the absorption coefficient: $ L \to S/a$ . This result for the asymptotic behavior of the radiance with distance also holds when scattering is present, as shown on the page for an analytical asymptotic solution for internal sources.

Approximate Analytical Solutions

A number of approximate analytical solutions to the RTE can be derived after simplifying the RTE in various ways. One of these approximate solutions is the single-scattering approximation (SSA). This solution is developed by assuming that the water in homogeneous, the sea surface is level, the sun is a point source in a black sky, there are no internal sources, and only single-scattering of photons is considered. The quasi-single-scattering approximation (QSSA) is a further simplification of the SSA. These are discussed in detail in Level 2 and in Gordon (1994). A much more complicated solution including second-order scattering (i.e., including photons that have been scattered twice) is developed in Walker (1994), Section 2-6.

These approximate solutions of the RTE are useful for isolating the main factors influencing underwater radiances. However, the solutions depend on various simplifying assumptions and the predicted radiances are typically accurate to a few tens of percent at best, and can be off by an order of magnitude.

Numerical Solutions

If accurate solutions of the RTE are to be obtained for realistic oceanic conditions, numerical methods must be used to solve the RTE. Many such methods have been developed. Some of these solution techniques have been tailored to specific environments, such as stellar or planetary atmospheres, and are not commonly used in oceanography. The numerical methods most commonly employed in oceanographic radiative transfer, and their salient characteristics, can be summarized as follows:

Monte Carlo
  • based on conceptually simple physics that mimics how nature absorbs and scatters photons
  • completely general; can solve time-dependent and 3D problems with arbitrary geometry
  • mathematical "tricks" to speed up the calculations can be sophisticated
  • easy to program
  • computed radiances have statistical errors, which can be reduced by tracing more photons (requiring longer computer times)
  • computer run times can be extremely slow for some problems (e.g., solving the RTE to large optical depths; run times increase exponentially with optical depth)

Discrete Ordinates
  • highly mathematical
  • difficult to program
  • does not handle highly peaked scattering phase functions well
  • models the medium as a stack of homogeneous layers
  • is fast for irradiance calculations and homogeneous water, but can be slow for radiances or if many layers are needed to resolve depth-dependent IOPs

Invariant Imbedding
  • highly mathematical
  • difficult to program
  • can solve only 1D problems (the one dimension being the depth in optical oceanography)
  • includes all orders of multiple scattering
  • computed radiances do not have statistical errors
  • is extremely fast (run times increase linearly with optical depth)

Because of their simplicity and generality, Monte Carlo methods are widely used to solve the RTE in fields as diverse as oceanography, atmospheric sciences, astronomy, medical physics, and nuclear engineering. The trade-off for their simplicity and generality is long computer run times for many problems. Nevertheless, Monte Carlo methods warrant a chapter of their own. The scattering phase functions for atmospheric aerosols are not as highly peaked at very small scattering angles as are those for oceanic particles. Discrete ordinates can handle aerosol phase functions well and is often used in atmospheric optics, but is not much used for underwater calculations because of the need to resolve highly peaked phase functions and to have many layers if the IOPs vary greatly with depth. Invariant imbedding is the solution technique used by the HydroLight numerical model, which is widely used in oceanography. Each of these techniques gives the same answer for the same inputs and boundary conditions for the RTE, for problems where all three techniques are applicable, as can be seen in the model comparison study of Mobley et al. (1993). They differ only in their internal mathematics and the resulting computer run times, and well-debugged and validated computer programs exist for each. In this sense, solving the RTE in the oceanographic setting can be considered a "solved problem."