A semi-analytic technique to speed up successive order of scattering model for optically thick media

A semi-analytic technique has been developed to speed up the integration of radiative transfer over optically thick media for the successive order of scattering (SOS) method. Based on the characteristics of the internal distribution of scattering intensity, this technique uses piecewise analytic eigenfunctions to fit internal scattering intensities and integrates them analytically over optical depth. This semi-analytic approach greatly reduces the number of sub-layers required for accurate radiative transfer calculation based on the SOS method. Results show that an accuracy of 1% for both flux and radiance (polar angle less than 67^°) can be achieved with a significantly small number of layers. This technique is accurate and efficient and makes the SOS method applicable for optically thick scattering media.     

Multiple-scaling methods for Monte Carlo simulations of radiative transfer in cloudy atmosphere

Two multiple-scaling methods for Monte Carlo simulations were derived from integral radiative transfer equation for calculating radiance in cloudy atmosphere accurately and rapidly. The first one is to truncate sharp forward peaks of phase functions for each order of scattering adaptively. The truncated functions for forward peaks are approximated as quadratic functions; only one prescribed parameter is used to set maximum truncation fraction for various phase functions. The second one is to increase extinction coefficients in optically thin regions for each order scattering adaptively, which could enhance the collision chance adaptively in the regions where samples are rare. Several one-dimensional and three-dimensional cloud fields were selected to validate the methods. The numerical results demonstrate that the bias errors were below 0.2% for almost all directions except for glory direction (less than 0.4%) and the higher numerical efficiency could be achieved when quadratic functions were used. The second method could decrease radiance noise to 0.60% for cumulus and accelerate convergence in optically thin regions. In general, the main advantage of the proposed methods is that we could modify the atmospheric optical quantities adaptively for each order of scattering and sample important contribution according to the specific atmospheric conditions.     

Fast and accurate 4 and 6 stream linearized discrete ordinate radiative transfer models for ozone profile retrieval

The global and long-term measurement of ozone vertical and horizontal distributions is one of the most important tasks in the monitoring of the earth's atmosphere. A number of satellite instruments are capable of delivering ozone profile distributions from UV nadir backscatter measurements. Retrieval algorithms should be efficient enough to deliver profiles in real-time without compromising accuracy. Such algorithms require a radiative transfer model that can generate quickly and accurately both simulated radiances and Jacobian matrices of weighting functions. We develop fast and analytic 4 stream and 6 stream linearized discrete ordinate models designed to satisfy performance and accuracy requirements for such an algorithm. The models have the pseudo-spherical treatment of the direct beam attenuation. For anisotropic scattering we use the delta-M scaling method to deal with strong forward scattering peaks. We demonstrate that the accuracy of the models is improved greatly upon application of a single scatter correction based on an exact specification of the phase function. For wide-angle off-nadir viewing, a sphericity correction is developed to deal more precisely with attenuation in a curved atmosphere. Radiances and weighting functions for the 4 and 6 stream models are compared with 20 stream output from the LIDORT model. We show that for the UV range pertinent to ozone profile retrieval from space, the 4 stream model generates backscatter radiances to an accuracy >1.25% for all viewing situations in a clear sky Rayleigh and background aerosol reference atmosphere, and up to 1.75% for a number of special scenarios with optically thick particulate layers. Six stream radiances are accurate to the 0.25% level for clear sky situations, and 0.65% for the special cases; weighting functions for the 6 stream output are accurate to ±2% in all cases. We discuss the implications of these comparisons regarding the performance and accuracy of the radiative transfer forward model in the ozone profile retrieval context.     

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.     

SASKTRAN: A spherical geometry radiative transfer code for efficient estimation of limb scattered sunlight

The inversion of satellite-based observations of limb scattered sunlight for the retrieval of constituent species requires an efficient and accurate modelling of the measurement. We present the development of the SASKTRAN radiative transfer model for the prediction of limb scatter measurements at optical wavelengths by method of successive orders along rays traced in a spherical atmosphere. The component of the signal due to the first two scattering events of the solar beam is accounted for directly along rays traced in the three-dimensional geometry. Simplifying assumptions in successive scattering orders provide computational optimizations without severely compromising the accuracy of the solution. SASKTRAN is designed for the analysis of measurements from the OSIRIS instrument and the implementation of the algorithm is efficient such that the code is suitable for the inversion of OSIRIS profiles on desktop computers. SASKTRAN total limb radiance profiles generally compare better with Monte-Carlo reference models over a large range of solar conditions than the approximate spherical and plane-parallel models typically used for inversions.     

Evaluation method for radiative heat transfer in polydisperse water droplets

Simplifications of the model for nongray radiative heat transfer analysis in participating media comprised of polydisperse water droplets are presented. Databases of the radiative properties for a water droplet over a wide range of wavelengths and diameters are constructed using rigorous Mie theory. The accuracy of the radiative properties obtained from the database interpolation is validated by comparing them with those obtained from the Mie calculations. The radiative properties of polydisperse water droplets are compared with those of monodisperse water droplets with equivalent mean diameters. Nongray radiative heat transfer in the anisotropic scattering fog layer, including direct and diffuse solar irradiations and infrared sky flux, is analyzed using REM². The radiative heat fluxes within the fog layer containing polydisperse water droplets are compared with those in the layer containing monodisperse water droplets. Through numerical simulation of the radiative heat transfer, polydisperse water droplets can be approximated by using the Sauter diameter, a technique that can be useful in several research fields, such as engineering and atmospheric science. Although this approximation is valid in the case of pure radiative transfer problems, the Sauter diameter is reconfirmed to be the appropriate diameter for approximating problems in radiative heat transfer, although volume-length mean diameter shows better accordance in some cases. The CPU time for nongray radiative heat transfer analysis with a fog model is evaluated. It is proved that the CPU time is decreased by using the databases and the approximation method for polydisperse particulate media.     

Comparison of single and multiple scattering approaches for the simulation of limb-emission observations in the mid-IR

The validity of single scattering radiative transfer calculations for simulation of limb-emission measurements of clouds in the mid-infrared spectral region was investigated by comparison with a multiple scattering model. For in limb direction optically thin clouds, like polar stratospheric clouds, errors of the single scattering scheme range below 3%. For optically thick clouds deviations are below 3% in case of low single scattering albedo (ω₀=0.24) increasing up to 10-30% for ω₀=0.84. Clouds which are optically thick in limb, but thin in nadir direction, can cause limb radiances which are by a factor of 1.7 higher than the blackbody radiance at cloud altitude.     

Radiative transfer with internal reflection via the converged discrete ordinates method

In response to the challenge of establishing highly accurate solutions to the plane layer radiative transfer equation with the simplest of methods, the converged discrete ordinates method is presented. With this algorithm of only finite difference, quadrature and acceleration, we show how to obtain highly accurate intensities for radiative transfer in a finite layer with internal surface reflection. The method features angular smoothing and angular interpolation through “faux” quadrature. In addition, a manufactured solution demonstrates the high accuracy of the method for forward peaked scattering. We consider scattering in a heterogeneous medium as a final demonstration.     

On the accuracy of double scattering approximation for atmospheric polarization computations

Interpretation of multi-angle spectro-polarimetric data in remote sensing of atmospheric aerosols requires fast and accurate methods of solving the vector radiative transfer equation (VRTE). The single and double scattering approximations could provide an analytical framework for the inversion algorithms and are relatively fast; however accuracy assessments of these approximations for the aerosol atmospheres in the atmospheric window channels have been missing. This paper provides such analysis for a vertically homogeneous aerosol atmosphere with weak and strong asymmetry of scattering. In both cases, the double scattering approximation gives a high accuracy result (relative error ～0.2%) only for the low optical path ～10^?2. As the error rapidly grows with optical thickness, a full VRTE solution is required for the practical remote sensing analysis. It is shown that the scattering anisotropy is not important at low optical thicknesses neither for reflected nor for transmitted polarization components of radiation.     

Monte Carlo ray-tracing simulation for radiative heat transfer in turbulent fluctuating media under the optically thin fluctuation approximation

The Monte Carlo ray-tracing method (MCRT) based on the concept of radiation distribution factor is extended to solve radiative heat transfer problem in turbulent fluctuating media under the optically thin fluctuation approximation. A one-dimensional non-scattering turbulent fluctuating media is considered, in which the mean temperature and absorption coefficient distribution are assumed and the shape of probability density function is given. The distribution of the time-averaged volume radiation heat source is solved by MCRT and direct integration method. It is shown that the results of MCRT based on the concept of radiation distribution factor agree with these of integration solution very well, but results of MCRT based on the concept of radiative transfer coefficient do not agree with these of integration solution. The solution of time-averaged radiative transfer equation by the concept of radiative transfer coefficient should be treated with caution.     

Tau method approximation for radiative transfer problems in a slab medium

An approximate method for solving the radiative transfer equation in a slab medium with an isotropic scattering is proposed. The method is based upon constructing the double Legendre series to approximate the required solution using Legendre tau method. The differential and integral expressions which arise in the radiative transfer equation are converted into a system of linear algebraic equations which can be solved for the unknown coefficients. Numerical examples are included to demonstrate the validity and applicability of the method and a comparison is made with existing results.     

Solution of the equation of transfer for coherent anisotropic scattering by double interval spherical harmonic method

The double interval spherical harmonic method introduced effectively by Wilson and Sen has already been used by Ghosh and Karanjai to solve the equation of radiative transfer in coherent isotropic scattering atmosphere, originally developed by Woolley and Stibbs. The same method has been successfully used in this paper to solve the equation of transfer for coherent anisotropic scattering.     

Curved ray tracing method for one-dimensional radiative transfer in the linear-anisotropic scattering medium with graded index

The curved ray tracing method (CRT) is extended to radiative transfer in the linear-anisotropic scattering medium with graded index from non-scattering medium. In this paper, the CRT is presented to solve one-dimensional radiative transfer in the linear-anisotropic scattering gray medium with a linear refractive index and two black boundaries. The predicted temperature distributions and radiative heat flux at radiative equilibrium are determined by the proposed method, and numerical results are compared with the data in references. The results show that the CRT has a good accuracy for radiative transfer in the linear-anisotropic scattering medium with graded index and the dimensionless emissive power and dimensionless radiative heat flux depend on the dimensionless refractive index gradient. It can also be seen that the dimensionless refractive index gradient has important effects on the temperature discontinuity at the boundaries.     

Three-dimensional radiative transfer in an anisotropically scattering, semi-infinite medium: generalized reflection function

Three-dimensional radiative transfer in an anisotropic scattering medium exposed to spatially varying, collimated radiation is studied. The generalized reflection function for a semi-infinite medium with a very general scattering phase function is the focus of this investigation. An integral transform is used to reduce the three-dimensional transport equation to a one-dimensional form, and a modified Ambarzumian's method is applied to formulate a nonlinear integral equation for the generalized reflection function. The integration is over both the polar and azimuthal angles; hence, the integral equation is said to be in the double-integral form. The double-integral, reflection function formulation can handle a variety of anisotropic phase functions and does not require an expansion of the phase function in a Legendre polynomial series. Complicated kernel transformations of previous single-integral studies are eliminated. Single and double scattering approximations are developed. Numerical results are presented for a Rayleigh phase function to illustrate the computational characteristics of the method and are compared to results obtained with the single-integral method. Agreement between the two approaches is excellent; however, as the transform variable increases beyond five the number of quadrature points required for the double-integral method to produce accurate solutions significantly increases. A new interpolation scheme produces accurate results when the transform variable is large.     

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.     

A full spectrum k-distribution based non-gray radiative property model for unburnt char

By using the concept of weighted sum of four gray particles and spectrum k-distribution (WSGP-SK), a non-gray radiative property model for unburnt char particles is developed. Based on the carbon burnout kinetic model for structure during oxidation, and the linear mixed approximation theory for complex index of refraction, spectral radiative properties of unburnt char particles are first calculated as function of the burnout ratio by Mie theory. Referring to the full spectrum k-distribution model, k-distribution is applied to reorder absorption and scattering efficiencies of particles. Then, weighting factors and efficiency factors of the non-gray radiative property model are directly obtained from Gaussian integral points of k-distribution. The model is validated against the benchmark solutions of line-by-line (LBL) model. Maximum relative errors of this model are 3% and 15% for radiative heat fluxes and source terms in non-isothermal inhomogeneous particulate media, respectively. The assumption of linearly varying radiative properties with burnout ratio (Lockwood et al. 1986) will result in a predicted deviation of 53% for radiative source terms. Results also show that this non-gray model is remarkably better than the Planck mean method. Moreover, a satisfactory comparison with LBL solutions is achieved in the gas and particle mixture by combining the non-gray WSGG-SK model (Guo et al. 2015). As a radiation sub-model, this non-gray radiative property model can significantly improve prediction accuracy of radiative heat transfer in oxy-fuel combustion.     

Finite element method for radiation heat transfer in multi-dimensional graded index medium

In graded index medium, ray goes along a curved path determined by Fermat principle, and curved ray-tracing is very difficult and complex. To avoid the complicated and time-consuming computation of curved ray trajectories, a finite element method based on discrete ordinate equation is developed to solve the radiative transfer problem in a multi-dimensional semitransparent graded index medium. Two particular test problems of radiative transfer are taken as examples to verify this finite element method. The predicted dimensionless net radiative heat fluxes are determined by the proposed method and compared with the results obtained by finite volume method. The results show that the finite element method presented in this paper has a good accuracy in solving the multi-dimensional radiative transfer problem in semitransparent graded index medium.     

Radiative perturbation theory for polarized radiances: 1. Theory

Radiative perturbation theory has proven to be a useful tool in radiative transfer calculations, especially in situations where repeated solution of the radiative transfer equation is required. So far however, its use has been restricted to non-polarized situations, including such applications as surface fluxes, UV indices, and the inversion of satellite radiance observations. Here, we extend the structure of radiative perturbation theory to incorporate the full Stokes formalism of polarization, to obtain the relevant equations for the first order term. This formalism will be applied to fluxes in a follow-up paper, and eventually to satellite observations.