Three-dimensional vector radiative transfer in a semi-infinite, Rayleigh scattering medium exposed to spatially varying polarized radiation: generalized reflection matrix

Three-dimensional vector radiative transfer in a semi-infinite medium exposed to spatially varying, polarized radiation is studied. The problem is to determine the generalized reflection matrix for a multiple scattering medium characterized by a 4×4 scattering matrix. A double integral transform is used to convert the three-dimensional vector radiative transfer equation to a one-dimensional form, and a modified Ambarzumian's method is then applied to derive a nonlinear integral equation for the generalized reflection matrix. The spatially varying backscattered radiation for an arbitrarily polarized incident beam can be found from the generalized reflection matrix. For Rayleigh scattering and normal incidence and emergence, the generalized reflection matrix is shown to have five non-zero elements. Benchmark results for these five elements are presented and compared to asymptotic results. When the incident radiation is polarized, the vector approach used in this study correctly predicts three-dimensional behavior, while the scalar approach does not. When the incident radiation is unpolarized, both the vector and scalar approaches predict a two-dimensional distribution of the intensity, but the error in the scalar prediction can be as high as 20%.     

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

The focus of this study is the generalized reflection function of multidimensional radiative transfer. The physical situation considered is spatially varying, collimated radiation incident on the upper boundary of an isotropically scattering, semi-infinite medium. 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 resulting equation is said to be in double-integral form because the integration is over both angular variables. Computational issues associated with this generalized reflection function formulation are investigated. The source function and reflection function formulations are compared, and the relative merits of the two approaches are discussed.     

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.     

Three-dimensional radiative transfer in an anisotropically scattering, plane-parallel medium: generalized reflection and transmission functions

The problem of spatially varying, collimated radiation incident on an anisotropically scattering, plane-parallel medium is considered. A very general phase function is allowed. An integral transform is used to reduce the three-dimensional radiative transport equation to a one-dimensional form, and a modified Ambarzumian's method is applied to derive nonlinear integral and integro-differential equations for the generalized reflection and transmission functions. The integration is over the polar and azimuthal angles—this formulation is referred to as the double-integral formulation. The integral equations are used to illustrate symmetry relationships and to obtain single- and double-scattering approximations. The generalized reflection and transmission functions are important in the construction of the solutions to many multidimensional problems. Coupled integral equations for the interior and emergent intensities are developed and, for the case of two identical homogeneous layers, used to formulate a doubling procedure. Results for an isotropic and Rayleigh scattering medium are presented to illustrate the computational characteristics of the formulation.     

Emergent intensity from a plane-parallel medium exposed to a Gaussian laser beam: A comparison of Rayleigh and isotropic scattering

The focus of this two-dimensional study is the radially varying intensity emergent from a plane-parallel scattering medium exposed to a collimated, Gaussian laser beam directed perpendicular to the upper surface. The method of analysis is the integral transform technique. Specifically, this work uses the generalized reflection and transmission functions from a previous study to construct the emergent intensity with the use of an inverse Hankel transform. Radially varying backscattered and transmitted intensities are calculated for media with isotropic and Rayleigh scattering phase functions and optical thicknesses that range from 0.125 to 8.0. The behavior of the emergent radiation inside and outside the beam is investigated for both narrow and wide beams. A new integration method is implemented to compute the emergent intensity at the beam center. The emergent intensity at the beam center is used to quantify when a one-dimensional model may be used. As expected, for small optical thicknesses and near the beam the phase function has significant influence, while far from the beam multiple scattering reduces the influence of the Rayleigh phase function. Results from this study will be useful in understanding and interpreting more complicated situations, such as those that include polarization.     

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.     

Spectral element method for vector radiative transfer equation

A spectral element method (SEM) is developed to solve polarized radiative transfer in multidimensional participating medium. The angular discretization is based on the discrete-ordinates approach, and the spatial discretization is conducted by spectral element approach. Chebyshev polynomial is used to build basis function on each element. Four various test problems are taken as examples to verify the performance of the SEM. The effectiveness of the SEM is demonstrated. The h and the p convergence characteristics of the SEM are studied. The convergence rate of p-refinement follows the exponential decay trend and is superior to that of h-refinement. The accuracy and efficiency of the higher order approximation in the SEM is well demonstrated for the solution of the VRTE. The predicted angular distribution of brightness temperature and Stokes vector by the SEM agree very well with the benchmark solutions in references. Numerical results show that the SEM is accurate, flexible and effective to solve multidimensional polarized radiative transfer problems.     

Multi-coupled single scattering method of solving vector radiative transfer equations

<正>A new method of multi-coupled single scattering(MCSS) for solving a vector radiative transfer equation is developed and made public on Internet.Recent solutions from Chandrasekhar’s X-Y method is used to validate the MCSS’s result,which shows high precision.The MCSS method is theoretically simple and clear,so it can be easily and credibly extended to the simulation of aerosol/cloud atmosphere’s radiative properties,which provides effective support for research into polarized remote sensing.     

Discrete-ordinate radiative transfer in a stratified medium with first-order rotational Raman scattering

Rotational Raman scattering (RRS) by air molecules in the Earth's atmosphere is predominantly responsible for the Ring effect: Fraunhofer and absorption-feature filling-in observed in UV/visible backscatter spectra. Accurate determination of RRS effects requires detailed radiative transfer (RT) treatment. In this paper, we demonstrate that the discrete-ordinate RT equations may be solved analytically in a multi-layer multiple scattering atmosphere in the presence of RRS treated as a first-order perturbation. Based on this solution, we develop a generic pseudo-spherical RT model LIDORT-RRS for the determination of backscatter radiances with RRS included; the model will generate output at arbitrary viewing geometry and optical thickness. Model comparisons with measured RRS filling-in effects from OMI observations show very good agreement. We examine telluric RRS filling-in effects for satellite-view backscatter radiances in a spectral range covering the ozone Huggins absorption bands. The model is also used to investigate calcium H and K Fraunhofer filling-in through cloud layers in the atmosphere.     

Time-dependent radiation transfer in a semi-infinite stochastic medium with Rayleigh scattering

The time-dependent radiation transfer in a semi-infinite stochastic medium of binary Markovian mixture with Rayleigh scattering is presented. A formalism, developed to treat radiation transfer in statistical mixtures, is used to obtain the ensemble-averaged solution. The average reflectivity, radiant energy and net flux are computed for specular-reflecting boundary. For the sake of comparison, we use two different weight functions in our calculations.     

Two-dimensional rayleigh scattering in a semi-infinite cylindrical medium exposed to a laser beam

A modification of Ambarzumian's method is used to develop the integro-differential equations for the source function, flux, and intensity at the boundary of a two-dimensional, semi-infinite cylindrical medium with second order Legendre phase function scattering. The incident radiation is collimated, normal to the top surface, and is dependent only on the radial coordinate. Boundary conditions which vary as a Bessel function and as a Gaussian distribution are investigated. The Gaussian distribution approximates a laser beam. Numerical results are presented in graphical and tabular forms for a Rayleigh scattering medium. The results are compared with those of isotropic scattering.     

Time-dependent radiation transfer with Rayleigh scattering in finite slab media

The time-dependent radiation transfer equation in a finite plane geometry with Rayleigh scattering is studied. The traveling wave transformation is used to obtain the corresponding stationary-like equation. Pomraning-Eddington approximation is then used to find the solution. Numerical results for reflectivity at the left boundary and transmissivity from the right boundary are presented at different times. The medium is assumed to have specular-reflecting boundaries with angular-dependent externally incident flux. Two different weight functions are introduced to force the boundary conditions to fulfill.     

Two-dimensional linearly anisotropic scattering in a semi-infinite cylindrical medium exposed to a laser beam

A modification of Ambarzumian's method is used to develop the integro-differential equations for the source function, flux, and intensity at the boundary of a two-dimensional, semi-infinite cylindrical medium which scatters linearly. The incident radiation is collimated, normal to the top surface of the medium, and is dependent only on the radial coordinate. The radial variation is assumed to be a Bessel function or a Gaussian distribution. The Gaussian boundary condition is used to simulate a laser beam. Numerical results are presented in graphical and tabular forms for both boundary conditions. Results for forward and backward scattering phase functions are compared with those for isotropic scattering. A method is presented for extending these results to the problem of a strongly anisotropic phase function which is made up of a spike in the forward direction superimposed on a linear phase function.     

Three-dimensional radiative transfer in an isotropically scattering, plane-parallel medium: generalized X- and Y-functions

The topic of this work is the generalized X- and Y-functions of multidimensional radiative transfer. The physical problem considered is spatially varying, collimated radiation incident on the upper boundary of an isotropically scattering, plane-parallel medium. An integral transform is used to reduce the three-dimensional transport equation to a one-dimensional form, and a modified Ambarzumian's method is used to derive coupled, integro-differential equations for the source functions at the boundaries of the medium. The resulting equations are said to be in double-integral form because the integration is over both angular variables. Numerical results are presented to illustrate the computational characteristics of the formulation.     

Infrared radiative transfer modelling in a 3D scattering cloudy atmosphere: Application to limb sounding measurements of cirrus

The Monte Carlo cloud scattering forward model (McClouds_FM) has been developed to simulate limb radiative transfer in the presence of cirrus clouds, for the purposes of simulating cloud contaminated measurements made by an infrared limb sounding instrument, e.g. the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). A reverse method three-dimensional Monte Carlo transfer model is combined with a line-by-line model for radiative transfer through the non-cloudy atmosphere to explicitly account for the effects of multiple scattering by the clouds. The ice cloud microphysics are characterised by a size distribution of randomly oriented ice crystals, with the single scattering properties of the distribution determined by accurate calculations accounting for non-spherical habit.A comparison of McClouds_FM simulations and real MIPAS spectra of cirrus shows good agreement. Of particular interest are several noticeable spectral features (i.e. H₂O absorption lines) in the data that are replicated in the simulations: these can only be explained by upwelling tropospheric radiation scattered into the line-of-sight by the cloud ice particles.     

Three-dimensional polarized Monte Carlo atmospheric radiative transfer model (3DMCPOL): 3D effects on polarized visible reflectances of a cirrus cloud

A polarized atmospheric radiative transfer model for the computation of radiative transfer inside three-dimensional inhomogeneous mediums is described. This code is based on Monte Carlo methods and takes into account the polarization state of the light. Specificities introduced by such consideration are presented. After validation of the model by comparisons with adding-doubling computations, examples of reflectances simulated from a synthetic inhomogeneous cirrus cloud are analyzed and compared with reflectances obtained with the classical assumption of a plane parallel homogeneous cloud (1D approximation). As polarized reflectance is known to saturate for optical thickness of about 3, one could think that they should be less sensitive to 3D effects than total reflectances. However, at high spatial resolution (80 m), values of polarized reflectances much higher than the ones predicted by the 1D theory can be reached. The study of the reflectances of a step cloud shows that these large values are the results of illumination and shadowing effects similar to those often observed on total reflectances. In addition, we show that for larger spatial resolution (10 km), the so-called plane-parallel bias leads to a non-negligible overestimation of the polarized reflectances of about 7–8%.     

Reverse electric field Monte Carlo simulation for vector radiative transfer in the atmosphere

In this paper, a reverse electric field Monte Carlo （REMC） method is proposed to study the vector radiation transfer in the atmosphere. The REMC is based on tracing the multiply scattered electric field to simulate the vector transmitted radiance. The reflected intensities with different total optical depth values are obtained, which accord well with the results in the previous research. Stokes vector and the degree of polarization are numerically investigated. The simulation result shows that when the solar zenith angle is determined, the zenith angle of detector has two points, of which the degree of polarization does not change with the ground albedo and the optical depth. The two points change regularly with the solar zenith angle. Moreover, our REMC method can be applied to the vector radiative transfer in the atmosphere-ocean system.     

The Iterative-DRESOR method to solve radiative transfer in a plane-parallel, anisotropic scattering medium with specular-diffuse boundaries

Using the intensity with high directional resolution obtained by the Basic-DRESOR method as an initial guess, which is substituted into the integrated radiative transfer equation (IRTE), an iterative algorithm is proposed, called the Iterative-DRESOR method. This method can reduce the error levels of the intensity from several percent using the Basic-DRESOR method to a level of less than 1.0×10⁻⁶ with acceptable computation costs. The method is also validated against the exact heat flux in literature in some cases. It further clarifies some uncertain results for the reflectance in a pure, linearly anisotropic scattering medium with specular-diffuse boundaries. The directional distributions of intensity are obviously influenced by the reflecting modes of the boundary, especially in the zone near the boundary. The reflecting mode of an emitting boundary has little effect on the transmittance or reflectance. The reflecting mode of a non-emitting boundary also has little effect on the transmittance, but it obviously influences the reflectance. The difference between the reflectance for specular and diffuse boundaries increases at first, and then decreases, as the optical thickness of the medium increases. The difference will decrease as the scattering albedo of the medium increases, and it is negligible when the medium is pure scattering. The effect of the scattering phase function of the medium on the difference can also not be ignored. The Iterative-DRESOR method is expected to strengthen the capability of the Monte Carlo method to produce accurate results and to validate the results of other methods to solve RTE.     

On the solution of a vectorial radiative transfer equation in an arbitrary three-dimensional turbid medium with anisotropic scattering

The authors developed a numerical method of the boundary-value problem solution in the vectorial radiative transfer theory applicable to the turbid media with an arbitrary three-dimensional geometry. The method is based on the solution representation as the sum of an anisotropic part that contains all the singularities of the exact solution and a smooth regular part. The regular part of the solution could be found numerically by the finite element method that enables to extend the approach to the arbitrary medium geometry. The anisotropic part of the solution is determined analytically by the special form of the small-angle approximation. The method development is performed by the examples of the boundary-value problems for the plane unidirectional and point isotropic sources in a turbid medium slab.