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.     

Discrete-ordinate method with matrix exponential for a pseudo-spherical atmosphere: Vector case

The paper is devoted to the extension of the matrix-exponential formalism for the scalar radiative transfer to the vector case. Using basic results of the theory of matrix-exponential functions we provide a compact and versatile formulation of the vector radiative transfer. As in the scalar case, we operate with the concept of the layer equation incorporating the level values of the Stokes vector. The matrix exponentials which enter in the expression of the layer equation are computed by using the matrix eigenvalue method and the Padé approximation. A discussion of the computational efficiency of the proposed method for both an aerosol-loaded atmosphere as well as a cloudy atmosphere is also provided.     

Matrix formulations of radiative transfer including the polarization effect in a coupled atmosphere-ocean system

A vector radiative transfer model has been developed for a coupled atmosphere-ocean system. The radiative transfer scheme is based on the discrete ordinate and matrix operator methods. The reflection/transmission matrices and source vectors are obtained for each atmospheric or oceanic layer through the discrete ordinate solution. The vertically inhomogeneous system is constructed using the matrix operator method, which combines the radiative interaction between the layers. This radiative transfer scheme is flexible for a vertically inhomogeneous system including the oceanic layers as well as the ocean surface. Compared with the benchmark results, the computational error attributable to the radiative transfer scheme has been less than 0.1% in the case of eight discrete ordinate directions. Furthermore, increasing the number of discrete ordinate directions has produced computations with higher accuracy. Based on our radiative transfer scheme, simulations of sun glint radiation have been presented for wavelengths of 670 nm and 1.6 μm. Results of simulations have shown reasonable characteristics of the sun glint radiation such as the strongly peaked, but slightly smoothed radiation by the rough ocean surface and depolarization through multiple scattering by the aerosol-loaded atmosphere. The radiative transfer scheme of this paper has been implemented to the numerical model named Pstar as one of the OpenCLASTR/STAR radiative transfer code systems, which are widely applied to many radiative transfer problems, including the polarization effect.     

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.     

Radiative transfer solutions for coupled atmosphere ocean systems using the matrix operator technique

Accurate radiative transfer models are the key tools for the understanding of radiative transfer processes in the atmosphere and ocean, and for the development of remote sensing algorithms. The widely used scalar approximation of radiative transfer can lead to errors in calculated top of atmosphere radiances. We show results with errors in the order of±8% for atmosphere ocean systems with case one waters. Variations in sea water salinity and temperature can lead to variations in the signal of similar magnitude. Therefore, we enhanced our scalar radiative transfer model MOMO, which is in use at Freie Universität Berlin, to treat these effects as accurately as possible. We describe our one-dimensional vector radiative transfer model for an atmosphere ocean system with a rough interface. We describe the matrix operator scheme and the bio-optical model for case one waters. We discuss some effects of neglecting polarization in radiative transfer calculations and effects of salinity changes for top of atmosphere radiances. Results are shown for the channels of the satellite instruments MERIS and OLCI from 412.5 nm to 900 nm.     

A direct method for the integration of the equation of radiative transfer in a turbid atmosphere

In this paper we apply a numerical method on the equation of radiative transfer in a turbid atmosphere. The solution is obtained by means of direct integration of the equation of radiative transfer without any circuitous series development.     

A vector radiative transfer model of coupled ocean-atmosphere system using matrix-operator method for rough sea-surface

A vector radiative transfer model termed PCOART has been developed for the coupled ocean-atmosphere system, using the matrix-operator (or adding-doubling) method, which considers the rough sea-surface. The theoretical formulations of the solution of the vector radiative transfer equation of the coupled ocean-atmosphere system, and the reflection-transmission matrices and internal radiation sources for rough sea surface are described. The model intercomparison is performed for several radiative transfer problems in the atmosphere and ocean, and the results show that PCOART can exactly predict the radiance fields for both flat and rough sea surface. Also, the polarizing remote sensing data from POLDER is used to test the capacity of PCOART to simulate the polarization radiance at the top-of-atmosphere, which shows that PCOART can perfectly reproduce the linear polarization reflectance measured by POLDER. PCOART can not only simulate the total radiance field in the coupled ocean-atmosphere system with wind-induced rough sea surface but also predict the polarization radiance field both in the atmosphere and in the ocean, which can serve as a good tool for the ocean optics and ocean color remote sensing communities.     

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.     

Exact solution of the standard transfer problem in a stellar atmosphere

We come back to the analytical solution of the standard transfer problem in a stellar atmosphere. It consists in solving the radiative transfer equation in a homogeneous and isothermal plane-parallel atmosphere, with light scattering taken as isotropic and monochromatic. The literature on the subject is reviewed and the existing solution in a finite slab is improved thanks to the introduction of non classical auxiliary functions. Eleven-figure tables of the solution are given for typical values of the input parameters currently met in stellar atmospheres.     

Attenuation of solar radiation by a water mist from the ultraviolet to the infrared range

A model is developed for the hemispherical transmittance of direct and scattered solar radiation from a cloudless atmosphere by a mist layer of water droplets in order to investigate the potential of water misting systems to serve as a protection from solar irradiation with particular emphasis on harmful UV radiation. The proposed model is based on published spectral experimental data for solar irradiation, Mie theory for interaction of the radiation with single spherical droplets, and radiative transfer theory. Known limiting solutions are employed to simplify the Mie calculations. The modified two-flux approximation is used to account for both direct and diffuse irradiation in lieu of a numerical solution for the full radiative transfer equation in anisotropically scattering media. The role of the governing parameters of a disperse water curtain of water droplets, water content, and droplet size for sample conditions is studied in some detail, particularly in the near-ultraviolet part of the spectrum where radiation can result in human tissue damage.     

一层随机非球形粒子高阶散射矢量辐射传输迭代解的参数求逆

根据一层随机非球形粒子矢量辐射传输(VRT)方程的一阶与高阶Mueller矩阵解的差异,提出了在已知一些参数的条件下,用水平与垂直同极化和交叉极化后向散射测量反演非球形粒子复介电常数和单位面积粒子数的迭代方法,并用极化后向散射强度的模拟测量数据及机载合成孔径雷达(Air SAR)实测数据进行了反演试验,并对反演值与各参数的关系以及测量噪声对反演的影响进行了讨论.     

Thermal analysis for transient radiative cooling of a conducting semitransparent layer of ceramic in high-temperature applications

A method is developed for obtaining transient temperature distribution in a cooling semitransparent layer of ceramic. The layer is emitting, absorbing, isotropically scattering and heat conducting with a refractive index ranging from 1 to 2. The solution involves solving simultaneously the energy equation and the integral equation for the radiative flux gradient. The energy equation is solved using an implicit finite volume scheme and the integral equation of radiative heat transfer is solved using the singularity technique and Gaussian integration. The effects of scattering are investigated. It is shown that scattering has a significant effect on the transient temperature distribution and the transient mean temperature of the layer.     

Iterative approach of high-order scattering solution for vector radiative transfer of inhomogeneous random media

By stratifying a random scatter media into multiple thin layers in the vertical z direction, the first-order scattering solution of each thin layer is employed to derive high-order scattering solution of whole random media. Using the Fourier transform and Mueller matrices in discrete ordinates, an iterative approach to solve high-order scattering solution of vector radiative transfer (VRT) equation is newly developed. Numerical results are well compared with the Mueller matrix solutions of the first order for a single layer medium, second order for a half-space, and the results of the discrete ordinate and eigen analysis method. It demonstrates our approach as feasible, effective and especially applicable to high-order solution of VRT for both bistatic scattering and thermal emission of inhomogeneous non-spherical scatter media.     

The Monte Carlo atmospheric radiative transfer model McArtim: Introduction and validation of Jacobians and 3D features

A new Monte Carlo atmospheric radiative transfer model is presented which is designed to support the interpretation of UV/vis/near-IR spectroscopic measurements of scattered Sun light in the atmosphere. The integro differential equation describing the underlying transport process and its formal solution are discussed. A stochastic approach to solve the differential equation, the Monte Carlo method, is deduced and its application to the formal solution is demonstrated. It is shown how model photon trajectories of the resulting ray tracing algorithm are used to estimate functionals of the radiation field such as radiances, actinic fluxes and light path integrals. In addition, Jacobians of the former quantities with respect to optical parameters of the atmosphere are analyzed. Model output quantities are validated against measurements, by self-consistency tests and through inter comparisons with other radiative transfer models.     

An exact solution to the multi-dimensional line transfer equation

The theory of generalized analytic functions is used to obtain an exact closed form analytical solution to a transfer problem for spectral line radiation in a multi-dimensional atmosphere. The multi-dimensional full-space and half-space Green's functions so obtained are quite general and may be used, along with the corresponding orthogonality relationships, to obtain solutions to any general multi-dimensional radiative transfer problem involving model two-level atoms. An application of the method using perturbation techniques is illustrated.     

非均匀散射层矢量辐射传输(VRT)方程高阶散射解的迭代法

将散射介质层在z轴方向划分成薄层，用薄层的一阶散射强度、Fourier变换和迭代方法求解散射介质整层的矢量辐射传输(VRT)方程的高阶散射解.该方法将一阶散射与高阶散射迭代结合起来，计算公式简明，可计算高阶迭代解，计算时间少.计算结果与一层均匀散射介质的VRT方程一阶Mueller矩阵解、半空间均匀散射介质二阶Mueller矩阵解、以及离散坐标-特征值特征矢量法的VRT热辐射的数值解作了全面的比较.提出并讨论了非均匀散射层主动与被动VRT方程的高阶解.本计算程序可以通用于非球形粒子多层结构及非均匀介质的散射和热辐射计算. 关键词： VRT方程 分层 迭代解     

A perturbation solution of the radiative transfer equation in a differentially moving atmosphere

We present a new technique for solving the radiative transfer equation in a differentially moving atmosphere. The method is based on a pertubation of the solution of the transfer problem in a static atmosphere. The perturbation technique may be applied with any method for solving the static atmosphere problem and leads to significant reductions in computer time and storage requirements.The method is flexible and may be used to solve problems involving depth dependence in any of the parameters of the transfer equation.     

The temporal Mueller matrix solution for polarimetric scattering from a layer of random non-spherical scatterers

The Mueller matrix solution of the vector radiative transfer equation with time dependence is derived in this paper. It is applied to simulation of polarimetric bistatic scattering from a layer of non-uniformly oriented, random non-spherical scatterers when a Gaussian plane pulse is incident upon. Co-polarized and cross-polarized bistatic scattering are numerically calculated. The pulse echoes are compared with the incidence, and demonstrate its functional dependence on the physical parameters of random medium, such as spatial orientation and fractional volume of scatterers, incidence angle and polarization, the layer depth and others.