Towards linearization of atmospheric radiative transfer in spherical geometry

We present a general approach for the linearization of radiative transfer in a spherical planetary atmosphere. The approach is based on the forward-adjoint perturbation theory. In the first part we develop the theoretical background for a linearization of radiative transfer in spherical geometry. Using an operator formulation of radiative transfer allows one to derive the linearization principles in a universally valid notation. The application of the derived principles is demonstrated for a radiative transfer problem in simplified spherical geometry in the second part of this paper. Here, we calculate the derivatives of the radiance at the top of the atmosphere with respect to the absorption properties of a trace gas species in the case of a nadir-viewing satellite instrument.     

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.     

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.     

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.     

Computation of Green's function for radiative transfer

In an accompanying paper, we develop the computational expressions for the higher order perturbation of the radiative transfer equation, and present some numerical results for typical cases. In this article, we discuss a number of issues regarding the implementation of the HOP computation: obtaining the Green's function, its expansion as a double series of Legendre polynomials, and obtaining the adjoint radiance of more general sources such as those for the fluxes at arbitrary altitudes. Examples of Green's function and its expansion coefficients are presented.     

Numerical test of an inverse polarized radiative transfer algorithm

A procedure is tested with which to determine the single-scattering albedo from polarization measurements of the angle-dependent intensity at two locations within, or on the boundaries of, a homogeneous finite or infinite atmosphere that scatters radiation according to the Rayleigh law with true absorption.     

Scaling algorithms for the calculation of solar radiative fluxes

We derived new scaling formulae based on the method of successive orders of scattering to calculate solar radiative flux. In this report, we demonstrate a multiple scaling method, in which we introduce scaling factors for each scattering order independently. The formula of radiative transfer by the method of successive orders of scattering cannot be solved rapidly except in the case of optically thin atmospheres. Then we further derived a double scaling method, which scales the ordinary radiative transfer equation by two scaling factors. We applied the double scaling method to two-stream and four-stream approximations of the discrete ordinates method. Comparing the results of the double scaling method with those of the delta-M method, we found that the double scaling method improved the accuracy of radiative fluxes at large solar zenith angles, especially in the optically thin region, and that in the region where multiple scattering dominates, its accuracy was comparable to that of the delta-M method. Once we determined the scaling factors appropriately, the double scaling method calculated radiative fluxes as rapidly as the delta-M method in the two-stream and four-stream approximations. This method, therefore, is useful for accurate computation of solar radiative fluxes in general circulation models.     

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.     

Vector Green's function algorithm for radiative transfer in plane-parallel atmosphere

Green's function is a widely used approach for boundary value problems. In problems related to radiative transfer, Green's function has been found to be useful in land, ocean and atmosphere remote sensing. It is also a key element in higher order perturbation theory. This paper presents an explicit expression of the Green's function, in terms of the source and radiation field variables, for a plane-parallel atmosphere with either vacuum boundaries or a reflecting (BRDF) surface. Full polarization state is considered but the algorithm has been developed in such way that it can be easily reduced to solve scalar radiative transfer problems, which makes it possible to implement a single set of code for computing both the scalar and the vector Green's function.     

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.     

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.     

Higher-order radiative perturbation theory

As a computationally effective tool, the first-order term of the radiative perturbation theory has been computed successfully, and has been applied in a number of areas. In this article, we develop the computational expressions for the higher-order terms of the perturbation expansion in a plane parallel atmosphere. These expressions are then implemented, and numerical results for some typical cases are presented. These results indicate that the computation is successful and that the higher-order terms are essential in cases where the first-order term alone cannot predict the perturbation with sufficient accuracy.     

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.     

Effect of the improvement of the HITRAN database on the radiative transfer calculation

The line parameters of the HITRAN 2004 have been updated, as compared with the older editions (the 2000 edition and the 1996 edition). In order to know the effect of the modifications on radiative transfer calculation with high spectral resolution, comparison in optical depth and radiance spectrum have been given between different editions. Four infrared spectral regions are selected, and they cover the three bands of atmospheric infrared sounder (AIRS) and one of geosynchronous imaging fourier transform spectrometer (GIFTS). The comparison has shown that the relative difference between HITRAN 2000 and 2004 and that between HITRAN 1996 and 2004 is decreasing. But the maximal discrepancy between the latest two editions in some spectral intervals is over 1%. It is important to estimate the error of calculation with the line parameters correctly or one has to use the new edition of HITRAN.     

Stability of explicit radiation-material coupling in radiative transfer calculations

We show that explicit radiation-material coupling, which is essentially always stable for infrared radiative transfer is conditionally stable in the high energy density regime. A linearized stability analysis is first performed for a simple infinite-medium problem that yields both a criterion for unconditional stability, a time-step restriction that applies for conditional stability, and a time-step criterion that always applies for non-oscillatory solutions. This analysis is then extended to include space dependence with the result that the system is always conditionally stable, but with a time step restriction somewhat different from the infinite-medium case. Nonetheless, the time step restriction for non-oscillatory solutions remains the same. Computations are presented that confirm the predictions of our analysis, and conclusions are given.     

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.     

Approximate polarization calculations with scalar radiative transfer theory

The Pomraning phase function can be used to perform approximate polarized Rayleigh transfer calculations with a scalar radiative transfer equation. The approximation is numerically tested for the albedo problem consisting of azimuthally independent radiation incident on a homogeneous semi-infinite atmosphere. The numerical tests were carried out with the same approach used by Viik (JQSRT 68 (2000) 319-326) to numerically test the approximate phase function for solving the Milne problem. Away from the surface the Pomraning phase function gives marginally better results for the diffuse radiation than the usual scalar Rayleigh phase function because it was derived from an asymptotic limit more appropriate for deeper locations in an atmosphere. For optical depths less than unity, though, the scalar Rayleigh approximation is better than the Pomraning approximation.     

Variational method for solving radiative transfer problems

The variational principle is used to solve two problems of radiative transfer. The first one is the temperature distribution and radiative heat flux for a plane layer of ceramic material. While the second is the calculation of the integral blankness degree of a sphere filled with dust of an arc steel-melting furnace. Numerical results obtained shows good agreement with the published data.     

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.     

Application of the spatial averaging theorem to radiative heat transfer in two-phase media

The spatial averaging theorem is applied to rigorously derive continuum-scale equations of radiative transfer in two-phase media consisting of arbitrary-type phases in the limit of geometrical optics. The derivations are based on the equations of radiative transfer and the corresponding boundary conditions applied at the discrete-scale to each phase, and on the discrete-scale radiative properties of each phase and the interface between the phases. The derivations confirm that radiative transfer in two-phase media consisting of arbitrary-type phases in the range of geometrical optics can be modeled by a set of two continuum-scale equations of radiative transfer describing the variation of the average intensities associated with each phase. Finally, a Monte Carlo based methodology for the determination of average radiative properties is discussed in the light of previous pertinent studies.