A linearized vector radiative transfer model for atmospheric trace gas retrieval

We present a plane parallel radiative transfer model for polarized light, that provides the intensity vector as well as the derivatives of the four Stokes parameters with respect to atmospheric trace gas profiles. These derivatives are essential for retrieval of height resolved trace gas information from satellite measurements of backscattered sunlight. The model uses the Gauss-Seidel iteration technique for solving the radiative transfer equation. For the first time, the forward-adjoint radiative perturbation theory is applied for the linearization of a radiative transfer model including polarization. The accuracy of the model is better than 0.025% for all four Stokes parameters and better than 0.03% for the derivatives.     

Extension of the discrete-ordinate algorithm and efficient radiative transfer calculation

As an accurate and efficient algorithm, the discrete-ordinate method (DOM) has been used to solve the radiative transfer problem of plane-parallel scattering atmosphere illuminated by a parallel beam, an idealized case of the sun, from above the atmosphere. In this paper, we extend this algorithm so that radiative problems of more general sources, such as parallel surface sources that illuminate with a parallel beam in any direction from any vertical position, and general surface sources that illuminate continuously in a hemisphere, can be solved. For a problem where intensity distributions are sought for a number of different sources within the same atmosphere-surface system, the intrinsic properties of DOM are used so that the time required for the solution for extra sources is reduced to a substantially small amount. In the case of parallel surface sources, numerical testing has shown that the amount can be reduced to as little as 15% of a full solution. Examples of applications are presented.     

Robust and accurate filtered spherical harmonics expansions for radiative transfer

We present a novel application of filters to the spherical harmonics (P_N) expansion for radiative transfer problems in the high-energy-density regime. The filter we use is based on non-oscillatory spherical splines and a filter strength chosen to (i) preserve the equilibrium diffusion limit and (ii) vanish as the expansion order tends to infinity. Our implementation is based on modified equations that are derived by applying the filter after every time step in a simple first-order time integration scheme. The method is readily applied to existing codes that solve the P_N equations. Numerical results demonstrate that the solution to the filtered P_N equations are (i) more robust and less oscillatory than standard P_N solutions and (ii) more accurate than discrete ordinates solutions of comparable order. In particular, the filtered P₇ solution demonstrates comparable accuracy to an implicit Monte Carlo solution for a benchmark hohlraum problem in 2D Cartesian geometry.     

Linearization of radiative transfer with respect to surface properties

The linearization of radiative transfer with respect to surface properties in the UV and visible part of the solar spectrum is presented. The proposed method is a rigorous extension of the radiative perturbation theory with respect to surface properties. Given the forward and adjoint intensity field, analytical expressions are provided for the linearization of any observable related to the radiation field with respect to surface properties characterized by Minnaert's and Lambertian bidirectional reflection distribution function. For the considered surface reflection characteristics, we also discuss an extension of the reduction approach of Chandrasekhar as an alternative linearization method. The suitability of both approaches for the combined retrieval of trace gas and surface properties from the backscattered sunlight in the UV and visible part of the spectrum is discussed. The authors come to the conclusion that the perturbation theory, for this purpose, represents the superior method because of its general applicability to any parameter characterizing the optical properties of the atmosphere and the underlying surface.     

Generalized form of the direct and adjoint radiative transfer equations

The main goal of this paper is to give a rigorous derivation of the generalized form of the direct (also referenced as forward) and adjoint radiative transfer equations. The obtained expressions coincide with expressions derived by Ustinov [Adjoint sensitivity analysis of radiative transfer equation: temperature and gas mixing ratio weighting functions for remote sensing of scattering atmospheres in thermal IR. JQSRT 2001;68:195-211]. However, in contrast to [Ustinov EA. Adjoint sensitivity analysis of radiative transfer equation: temperature and gas mixing ratio weighting functions for remote sensing of scattering atmospheres in thermal IR. JQSRT 2001;68:195-211] we formulate the generalized form of the direct radiative transfer operator fully independent from its adjoint. To illustrate the application of the derived adjoint radiative transfer operator we consider the angular interpolation problem in the framework of the discrete ordinate method widely used to solve the radiative transfer equation. It is shown that under certain conditions the usage of the solution of the adjoint radiative transfer equation for the angular interpolation of the intensity can be computationally more efficient than the commonly used source function integration technique.     

Sparse tensor spherical harmonics approximation in radiative transfer

The stationary monochromatic radiative transfer equation is a partial differential transport equation stated on a five-dimensional phase space. To obtain a well-posed problem, boundary conditions have to be prescribed on the inflow part of the domain boundary.     

ARTS, the atmospheric radiative transfer simulator

ARTS is a modular program that simulates atmospheric radiative transfer. The paper describes ARTS version 1.0, which is applicable in the absence of scattering. An overview over all major parts of the model is given: calculation of absorption coefficients, the radiative transfer itself, and the calculation of Jacobians. ARTS can be freely used under a GNU general public license.Unique features of the program are its scalability and modularity, the ability to work with different sources of spectroscopic parameters, the availability of several self-consistent water continuum and line absorption models, and the analytical calculation of Jacobians.     

Adjoint problem of radiative transfer for a pseudo-spherical atmosphere and general viewing geometries

We formulate the adjoint radiative transfer for a pseudo-spherical atmosphere and various retrieval scenarios. The single scattering radiance is computed in a spherical atmosphere by using the source integration technique, while for the multiple scattering radiance we formulate an one-dimensional adjoint radiative transfer equation in a plane-parallel atmosphere. The adjoint solution of the radiative transfer equation is obtained by employing the discrete ordinate method with matrix exponential. We provide an abbreviated derivation of our formalism as well as a discussion of the numerical implementation of the theory.     

Accelerated higher order radiative perturbation computation: an application of GDOM

In a previous paper by the authors, the computational expressions for the higher order terms of the radiative perturbation series have been developed, and numerical results have been obtained. In this paper, new computational expressions are developed that use the analytic Green's function and the GDOM code, which were developed recently by Qin and Box. Our analysis and numerical computations indicate that the new scheme dramatically improves the efficiency of the computation (by reducing the CPU-time to less than 2% of that required by the previous scheme), and the huge demand on memory usage has been completely removed. This new scheme allows us to expand the high order perturbation computations to more general cases, for example, to include azimuth dependence, and to use the perturbation theory in more potential applications.     

One-dimensional solar radiative transfer: Perturbation approach and its application to independent-pixel calculations for realistic cloud fields

The radiative transfer perturbation theory (RTPT), which has already been introduced in atmospheric radiative transfer several years ago, is applied to cloud related problems. The RTPT requires the solution of the radiative transfer equation in the forward and the adjoint mode. The basic principles of this technique are presented as well as its extensions to isotropic surface reflection and its conjunction with the Hermite interpolation. This set of methods is applied to different atmospheric conditions including realistic cloud scenes. The results are compared with the usual (forward) independent-pixel calculations with respect to errors of individual pixels and domain-averaged values. The RTPT turns out to be sufficiently accurate in the case the clouds’ internal vertical variations remain moderate. It is also shown that, depending on the specific radiative transfer problem, the RTPT can offer some advantages on computational speed. However, the limitations of the RTPT with regard to realistic clouds are addressed as well.     

Simultaneous derivation of intensities and weighting functions in a general pseudo-spherical discrete ordinate radiative transfer treatment

The retrieval of atmospheric constituents from measurements of backscattered light requires a radiative transfer forward model that can simulate both intensities and weighting functions (partial derivatives of intensity with respect to atmospheric parameters being retrieved). The radiative transfer equation is solved in a multi-layer multiply-scattering atmosphere using the discrete ordinate method. In an earlier paper dealing with the upwelling top-of-the-atmosphere radiation field, it was shown that a full internal perturbation analysis of the plane-parallel discrete ordinate solution leads in a natural way to the simultaneous generation of analytically-derived weighting functions with respect to a wide range of atmospheric variables. In the present paper, a more direct approach is used to evaluate explicitly all partial derivatives of the intensity field. A generalization of the post-processing function is developed for the derivation of weighting functions at arbitrary optical depth and stream angles for both upwelling and downwelling directions. Further, a complete treatment is given for the pseudo-spherical approximation of the direct beam attenuation; this is an important extension to the range of viewing geometries encountered in practical radiative transfer applications. The numerical model LIDORT developed for this work is able to generate intensities and weighting functions for a wide range of retrieval scenarios, in addition to the passive remote sensing application from space. We present a number of examples in an atmosphere with O₃ absorption in the UV, for satellite (upwelling radiation) and ground-based (downwelling radiation) applications. In particular, we examine the effect of various pseudo-spherical parameterizations on backscatter intensities and weighting functions with respect to O₃ volume mixing ratio. In addition, the use of layer-integrated multiple scatter output from the model is shown to be important for satellite instruments with wide-angle off-nadir viewing geometries.     

Adjoint sensitivity analysis of radiative transfer equation: 2. Applications to retrievals of temperature in scattering atmospheres in thermal IR

An approach to formulation of inversion algorithms for thermal sounding in the case of scattering atmosphere based on the adjoint equation of radiative transfer (Ustinov, JQSRT 68 (2001) 195, referred to as Paper 1 in the main text) is applied to temperature retrievals in the scattering atmosphere for the nadir viewing geometry. Analytical expressions for the weighting functions involving the integration of the source function are derived. Temperature weighting functions for a simple model of the atmosphere with scattering are evaluated and convergence to the case of pure atmospheric absorption is demonstrated. The numerical experiments on temperature retrievals are carried out to demonstrate the validity of the expressions obtained.     

VLIDORT: A linearized pseudo-spherical vector discrete ordinate radiative transfer code for forward model and retrieval studies in multilayer multiple scattering media

We describe a new vector discrete ordinate radiative transfer model with a full linearization facility. The VLIDORT model is designed to generate simultaneous output of Stokes vector light fields and their derivatives with respect to any atmospheric or surface property. We develop new implementations for the linearization of the vector radiative transfer solutions, and go on to show that the complete vector discrete ordinate solution is analytically differentiable for a stratified multilayer multiply scattering atmospheric medium. VLIDORT will generate all output at arbitrary viewing geometry and optical depth. The model has the ability to deal with attenuation of solar and line-of-sight paths in a curved atmosphere, and includes an exact treatment of the single scatter computation. VLIDORT also contains a linearized treatment for non-Lambertian surfaces. A number of performance enhancements have been implemented, including a facility for multiple solar zenith angle output. The model has been benchmarked against established results in the literature.     

Refractive index imaging from radiative transfer equation-based reconstruction algorithm: Fundamentals

We present the first reconstruction algorithm for refractive index imaging, which is based on the radiative transfer equation (RTE). An objective function is iteratively minimized to find a solution to the problem of inversion of the refractive index field. The function describes the discrepancies of the emerging light measurements on the surface of the sample to be probed with predicted data from the corresponding numerical model. The unknown refractive index field is updated within each reconstruction iteration according to a search direction on the index distribution given by the adjoint model to the RTE. In this paper, emphasis is placed on the theoretical aspects. Preliminary tests are demonstrated on generic phantoms.     

Fast multilevel radiative transfer

The vast majority of recent advances in the field of numerical radiative transfer relies on approximate operator methods better known in astrophysics as Accelerated Lambda-Iteration (ALI). A superior class of iterative schemes, in term of rates of convergence, such as Gauss-Seidel and successive overrelaxation methods were therefore quite naturally introduced in the field of radiative transfer by Trujillo Bueno and Fabiani Bendicho [A novel iterative scheme for the very fast and accurate solution of non-LTE radiative transfer problems. Astrophys J 1995;455:646]; it was thoroughly described for the non-LTE two-level atom case. We describe hereafter in details how such methods can be generalized when dealing with non-LTE unpolarised radiation transfer with multilevel atomic models, in monodimensional geometry.     

Layer-edge conditions for multislab atmospheric radiative transfer problems

We derive nonstandard layer-edge conditions for efficient solution of multislab atmospheric radiative transfer problems. We begin by defining a local radiative transfer problem on the lowermost layer of a multislab model atmosphere and we consider a standard discrete ordinates version of this local problem. We then make use of a recently developed computational method in order to derive layer-edge conditions involving incident, reflected and transmitted radiation. These layer-edge conditions for the lowermost layer are given in terms of inherent optical properties of the layer, the solar zenith angle and the quadrature set used in the discrete ordinates approach. They can be used to increase the efficiency of our computational method in solving practical problems in atmospheric radiative transfer. Moreover, they are amenable to incorporation into other discrete ordinates methods. To illustrate, we report numerical results for two atmospheric model problems.     

Many polarized radiative transfer models

This note is an introduction to the reprint of the 1991 JQSRT article “A new polarized atmospheric radiative transfer model” by K.F. Evans and G.L. Stephens. We discuss the significance of the article, how our two plane-parallel polarized radiative transfer codes came about, how our codes have been used, and more recent developments in polarized radiative transfer modeling.     

An analytic and numerical solution with spectral Green's function method for transport equation in spherical geometry

Since in many cases curvilinear geometry is more appropriate than cartesian geometry for precise modeling of the complex systems for reactor calculation, we have developed the spectral Green's function (SGF) method which is employed to obtain angular and scalar flux distributions in heterogeneous sphere geometry with isotropic scattering. In this study, we showed that the neutron transport problems of homogeneous spheres could be reduced to the solution of plane geometry equation.Finally, some results are discussed and compared with those already obtained by diamond difference scheme to test the accuracy of the results. The agreement is satisfactory. SGF method is very suitable for the numerical solution of the neutron transport equation with isotropic scattering.