Discrete-ordinated finite-element method for atmospheric radiative transfer and remote sensing

Advantages and disadvantages of modern discrete-ordinates finite-element methods for the solution of radiative transfer problems in meteorology, climatology, and remote sensing applications are evaluated. After the common basis of the formulation of radiative transfer problems in the fields of neutron transport and atmospheric optics is established, the essential features of the discrete-ordinates finite-element method are described including the limitations of the method and their remedies. Numerical results are presented for 1-D and 2-D atmospheric radiative transfer problems where integral as well as angular dependent quantities are compared with published results from other calculations and with measured data. These comparisons provide a verification of the discrete-ordinates results for a wide spectrum of cases with varying degrees of absorption, scattering, and anisotropic phase functions. Accuracy and computational speed are also discussed. Since practically all discrete-ordinates codes offer a builtin adjoint capability, the general concept of the adjoint method is described and illustrated by sample problems. Our general conclusion is that the strengths of the discrete-ordinates finite-element method outweight its weaknesses. We demonstrate that existing general-purpose discrete-ordinates codes can provide a powerful tool to analyze radiative transfer problems through the atmosphere, especially when 2-D geometries must be considered.