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.     

On vector radiative transfer equation in curvilinear coordinate systems

The differential operator of polarized radiative transfer equation is examined in case of homogeneous medium in Euclidean three-dimensional space with arbitrary curvilinear coordinate system defined in it. This study shows that an apparent rotation of polarization plane along the light ray with respect to the chosen reference plane for Stokes parameters generally takes place, due to purely geometric reasons. Analytic expressions for the differential operator of transfer equation dependent on the components of metric tensor and their derivatives are found, and the derivation of differential operator of polarized radiative transfer equation has been made a standard procedure. Considerable simplifications take place if the coordinate system is orthogonal.     

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.     

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.     

Development of a GPU-based high-performance radiative transfer model for the Infrared Atmospheric Sounding Interferometer (IASI)

Satellite-observed radiance is a nonlinear functional of surface properties and atmospheric temperature and absorbing gas profiles as described by the radiative transfer equation (RTE). In the era of hyperspectral sounders with thousands of high-resolution channels, the computation of the radiative transfer model becomes more time-consuming. The radiative transfer model performance in operational numerical weather prediction systems still limits the number of channels we can use in hyperspectral sounders to only a few hundreds. To take the full advantage of such high-resolution infrared observations, a computationally efficient radiative transfer model is needed to facilitate satellite data assimilation. In recent years the programmable commodity graphics processing unit (GPU) has evolved into a highly parallel, multi-threaded, many-core processor with tremendous computational speed and very high memory bandwidth. The radiative transfer model is very suitable for the GPU implementation to take advantage of the hardware’s efficiency and parallelism where radiances of many channels can be calculated in parallel in GPUs.In this paper, we develop a GPU-based high-performance radiative transfer model for the Infrared Atmospheric Sounding Interferometer (IASI) launched in 2006 onboard the first European meteorological polar-orbiting satellites, METOP-A. Each IASI spectrum has 8461 spectral channels. The IASI radiative transfer model consists of three modules. The first module for computing the regression predictors takes less than 0.004% of CPU time, while the second module for transmittance computation and the third module for radiance computation take approximately 92.5% and 7.5%, respectively. Our GPU-based IASI radiative transfer model is developed to run on a low-cost personal supercomputer with four GPUs with total 960 compute cores, delivering near 4 TFlops theoretical peak performance. By massively parallelizing the second and third modules, we reached 364× speedup for 1 GPU and 1455× speedup for all 4 GPUs, both with respect to the original CPU-based single-threaded Fortran code with the –O₂ compiling optimization. The significant 1455× speedup using a computer with four GPUs means that the proposed GPU-based high-performance forward model is able to compute one day’s amount of 1,296,000 IASI spectra within nearly 10 min, whereas the original single CPU-based version will impractically take more than 10 days. This model runs over 80% of the theoretical memory bandwidth with asynchronous data transfer. A novel CPU–GPU pipeline implementation of the IASI radiative transfer model is proposed. The GPU-based high-performance IASI radiative transfer model is suitable for the assimilation of the IASI radiance observations into the operational numerical weather forecast model.     

Optical properties of ice particles in young contrails

The single-scattering properties of four types of ice crystals (pure ice crystals, ice crystals with an internal mixture of ice and black carbon, ice crystals coated with black carbon, and soot coated with ice) in young contrails are investigated at wavelengths 0.65 and 2.13 μm using Mie codes for coated spheres. The four types of ice crystals show differences in their single-scattering properties because of the embedded black carbon whose volume ratio is assumed to be 5%. The bulk-scattering properties of young contrails consisting of the four types of ice crystals are further investigated by averaging their single-scattering properties over a typical ice particle size distribution found in young contrails. The effect of the radiative properties of the four types of ice particles on the Stokes parameters I, Q, U, and V is also investigated for different viewing zenith angles and relative azimuth angles with a solar zenith angle of 30° using a vector radiative transfer model based on the adding-doubling technique. The Stokes parameters at a wavelength of 0.65 μm show pronounced differences for the four types of ice crystals, whereas the counterparts at a wavelength of 2.13 μm show similar variations with the viewing zenith angle and relative azimuth angle. However, the values of the results for the two wavelengths are noticeably different.     

灰气体加权和辐射模型综合评估及分析

在O_2/CO_2气氛下,参与性介质的非灰气体辐射特性表现出不同于空气气氛下的特性,因此,非灰气体辐射模型的选择和应用在换热过程中将变得十分重要.基于统计窄谱带模型,本文综合评估近年发展应用较广的灰气体加权和(WSGG)模型.结果表明,几种WSGG模型的预测值总体趋势正确,但仍存在着一定的差别.对于发射率,Dorigon等(2013 Int.J.Heat Mass Transfer 64 863)和Bordbar等(2014 Combust.Flame 161 2435)的WSGG模型与基准模型符合较好,相对误差小于20%.与离散坐标法结合,本文求解了PH_2O/PCO_2=1,2时的一维平行平板间辐射换热问题.结果显示,由Dorigon等和Bordbar等的WSGG模型得到的辐射热源和热流密度分布的相对误差均较小(10%左右).Johansson等(2011 Combust.Flame 158893)和Bordbar等的WSGG模型具有更广的适用范围.     

An improved model to calculate radiative heat transfer in hot combustion gases

Radiative heat transfer plays an important role in the chemical reactions in the combustor. The widely used WSGG model proposed by Smith is established for normal pressure, which shows inevitable computational errors when dealing with radiative heat transfer problems at reduced or elevated pressures. In this paper, an improved global model is established to calculate the radiant energy exchanges between combustion gases and combustor chamber walls. Compared with the Smith model, the new model shows better performance in a wide range of pressure regions. The model accuracy is examined by computing the emissivity, radiative heat flux as well as the radiative source of H₂O–CO₂ gas mixtures at different pressure values. Finally, the radiative heat transfer inside a 3D TBCC(turbine-based combined cycle) engine exhaust system where strong gradients of pressure and temperature exist, is also addressed. The computational results show that the developed model provides approximate results at much less computational costs than the high-precision MSMGFSK-c8 model, which makes it competitive in complicated combustion systems.     

Fresnel boundary and interface conditions for polarized radiative transfer in a multilayer medium

In many applications of the theory of radiative transfer, it is important to consider the changes in the index of refraction that occur when the physical domain being studied consists of material regions with distinct optical properties. When polarization effects are taken into account, the radiation field is described by a vector of four components known as the Stokes vector. At an interface between two different material regions, the reflected and transmitted Stokes vectors are related to the incident Stokes vector by means of reflection and transmission matrices, which are derived from the Fresnel formulas for the amplitude coefficients of reflection and transmission. Having seen that many works on polarized radiative transfer that allow for changes in the index of refraction exhibit discrepancies in their expressions for the transmission matrix, we present in this work a careful derivation of the relations between the reflected and transmitted Stokes vectors and the Stokes vector incident on an interface. We obtain a general form of a transmission factor that is required to ensure conservation of energy and we show that most of the discrepancies encountered in existing works are associated with the use of improper forms of this factor. In addition, we derive explicit and compact expressions for the Fresnel boundary and interface conditions appropriate to the study of polarized radiative transfer in a multilayer medium.     

Sensitivity of trace gas abundances retrievals from infrared limb emission spectra to simplifying approximations in radiative transfer modelling

The accuracy of trace gas abundances retrieved from spectrally resolved infrared limb emission measurements relies, among other things, on the appropriate modelling of radiative transfer through the actual atmosphere. We quantify the mapping of several commonly applied simplifications in radiative transfer modelling on the trace gas abundances retrieval error at the example of the Michelson interferometer for passive atmospheric sounding (MIPAS)/environmental satellite (ENVISAT) space borne Fourier transform infra-red limb emission experiment. The Karlsruhe optimized and precise radiative transfer algorithm (KOPRA) which was used as a tool for this study will be introduced. KOPRA supports accurate modelling of the particular instrument requirements of MIPAS and the observation scenarios during the ENVISAT mission, in particular with respect to its viewing direction and its altitude coverage of the atmosphere. We show that disregarding of horizontal temperature inhomogeneities and non-local thermodynamic equilibrium effects, insufficient accuracy in modelling of field-of-view and apodisation effects, and disregarding individual profiles of isotopic species play the key roles in radiative transfer modelling and lead to systematic retrieval errors which can by far exceed the expected random retrieval error caused by measurement noise.     

Passive remote sensing of planetary atmospheres and retrievals of atmospheric macro- and microphysical parameters

In two recent papers, referred to as Paper 1 and Paper 2 in the main text, we have pursued an idea that atmospheric weighting functions for any geophysical parameters essentially consist of two kinds of entities: weighting functions for the atmospheric radiative parameters, which directly enter the radiative transfer equation, and partial derivatives of the atmospheric radiative parameters with respect to atmospheric geophysical parameters. Corresponding analysis was performed for non-scattering and scattering planetary atmospheres for nadir-viewing geometry and thermal spectral region. In the present paper we conduct a study of the second group of above entities, the partial derivatives of atmospheric radiative parameters with respect to atmospheric geophysical parameters. Along the way, we analyze the role of radiative parameters of individual atmospheric constituents, which serve as intermediaries between the total radiative parameters of the atmosphere and its geophysical parameters. The obtained expressions, combined with radiative weighting functions, are used in the end-to-end expressions for the weighting functions obtained for various specific geophysical parameters.     

Radiative gas dynamics of the Fire-II superorbital space vehicle

The rates of convective and radiative heating of the Fire-II reentry vehicle are calculated, and the results are compared with experimental flight data. The computational model is based on solving a complete set of equations for (i) the radiative gas dynamics of a physically and chemically nonequilibrium viscous heatconducting gas and (ii) radiative transfer in 2D axisymmetric statement. The spectral optical parameters of high-temperature gases are calculated using ab initio quasi-classical and quantum-mechanical methods. The transfer of selective thermal radiation in terms of atomic lines is calculated using the line-by-line method on a specially generated computational grid that is nonuniform in radiation wavelength.     

被动多轴差分吸收光谱技术监测大气NO_2垂直廓线研究

痕量气体垂直廓线的监测,对大气污染研究具有重要意义。介绍了被动多轴差分吸收光谱(MAX-DOAS)技术监测痕量气体垂直廓线的光学遥感方法。研究中MAX-DOAS测量多个角度的斜柱浓度、结合大气辐射传输模型,利用最优估算法反演出痕量气体垂直廓线,并对最优估算法参数选取和反演评估进行了详细描述。将该技术应用于合肥地区NO2垂直廓线的监测:通过与长光程差分吸收光谱仪的测量结果对比,相关系数达到0.80。该技术为大气环境立体监测提供了一种简便的方法。     

On the derivation of vector radiative transfer equation for polarized radiative transport in graded index media

Light transport in graded index media follows a curved trajectory determined by Fermat's principle. Besides the effect of variation of the refractive index on the transport of radiative intensity, the curved ray trajectory will induce geometrical effects on the transport of polarization ellipse. This paper presents a complete derivation of vector radiative transfer equation for polarized radiation transport in absorption, emission and scattering graded index media. The derivation is based on the analysis of the conserved quantities for polarized light transport along curved trajectory and a novel approach. The obtained transfer equation can be considered as a generalization of the classic vector radiative transfer equation that is only valid for uniform refractive index media. Several variant forms of the transport equation are also presented, which include the form for Stokes parameters defined with a fixed reference and the Eulerian forms in the ray coordinate and in several common orthogonal coordinate systems.     

气粒混合物非灰辐射特性合并宽窄谱带K分布模型

气粒混合物辐射问题具有全场性、非灰性、耦合性等特点,准确预估高温燃气/粒子非灰辐射特性是非常重要的。本文将合并宽窄谱带K分布模础(CWNBCK)与离散坐标法(DOM)结合,开展了非灰气粒混合物辐射换热问题的模拟工作,分别验证了一维和三维情况下应用该模型的准确性,给出不同工况下的热流源项、壁面热流或辐射热流等。结果表明:该模型能够给出与SNB模型精度基本相同的结果,考虑其计算效率的提高,可以在工程实际中应用该模型计算非灰气粒混合物辐射换热。     

Conversion issues between microwave radiance and brightness temperature

Microwave radiances are usually converted into brightness temperatures for data assimilation and retrievals. The Rayleigh-Jeans approximation has been believed to be a good approximation for the conversion at low frequencies, but inaccurate at high frequencies. However, the simplified radiative transfer models under the Rayleigh-Jeans approximation (hereafter referred as BT-RTE) have been successfully applied in radiance simulations for frequencies below 183 GHz, which has somewhat puzzled the radiative transfer community. This paper clarifies the confusion. In addition, the conversion formula for the third and the fourth Stokes components are derived.Simulations for a polarized sensor, the Special Sensor Microwave Imager and Sounder, show that the BT-RTE is generally accurate. Results for a polarimetric sensor, WINDSAT, show that the third and the fourth Stokes radiances should be converted using the exact conversion formula given in this study rather than using a direct Planck function conversion.     

The NIMO Monte Carlo model for box-air-mass factor and radiance calculations

A new fully spherical multiple scattering Monte Carlo radiative transfer model named NIMO (NIWA Monte Carlo model) is presented. The ray tracing algorithm is described in detail along with the treatment of scattering and absorption, and the simulation of backward adjoint trajectories. The primary application of NIMO is the calculation of box-air-mass factors (box-AMFs), which are used to convert slant column densities (SCDs) of trace gases, derived from UV–visible multiple axis Differential Optical Absorption Spectroscopy (MAX-DOAS) measurements, into vertical column densities (VCDs). Box-AMFs are also employed as weighting functions for optimal estimation retrievals of vertical trace gas profiles from SCDs. Monte Carlo models are well suited to AMF calculations at high solar zenith angles (SZA) and at low viewing elevation angles where multiple scattering is important. Additionally, the object-oriented structure of NIMO makes it easily extensible to new applications by plugging in objects for new absorbing or scattering species. Box-AMFs and radiances, calculated for various wavelengths, SZAs, viewing elevation and azimuth angles and aerosol scenarios, are compared with results from nine other models using a set of exercises from a recent radiative transfer model intercomparison. NIMO results for these simulations are well within the range of variability of the other models.     

大气折射对可见光波段辐射传输特性的影响

折射是影响辐射传输的重要因素. 为分析大气折射对辐射传输的影响, 基于Monte Carlo方法, 给出了考虑大气折射的矢量辐射传输模型, 实现了均匀气层和耦合面处光子随机运动过程的模拟, 实现了直射光及漫射光Stokes矢量、偏振度和辐射通量等参数的计算. 在考虑和不考虑大气折射两种条件下, 验证了模型的准确性; 在纯瑞利散射条件下, 讨论了大气折射对不同方向漫射光Stokes矢量的影响; 在不同太阳天顶角、大气廓线、气溶胶及含云大气条件下, 分析了大气折射对辐射传输过程的影响. 结果表明: 大气折射对漫射光Stokes矢量的影响主要体现在天顶角70°–110°区间, 且随着太阳入射角增大, 其影响更为显著; 不同大气廓线情形下, 大气折射对Stokes矢量的影响不一致, 其原因是不同大气廓线对应的折射率廓线存在差异. 含云及含气溶胶大气条件下, 大气折射对辐射传输的影响变弱, 沙尘型及海盐型气溶胶条件下, 折射对辐射传输的影响强于可溶型气溶胶情形; 不同形状气溶胶条件下, 大气折射对辐射传输的影响也存在显著差异; 不同云高条件下, 大气折射对漫射光Stokes矢量的影响无显著差异, 但随着云光学厚度增大, 大气折射的影响减弱.     

Coupled radiation and turbulent multiphase flow in an aluminised solid propellant rocket engine

We present in this paper numerical simulations of coupled radiative transfer and turbulent flows at high temperature and pressure, typical of multiphase flows encountered in aluminised solid propellant rocket engines. The radiating medium is constituted of gases and of liquid or solid particles of oxidised aluminum. The turbulent flow of the gaseous phase is treated by using a four equation, low Reynolds number, boundary-layer-type turbulence model. The distributions of concentrations, temperatures, and temperature fluctuation variances of particles are calculated from a Lagrangian approach and a turbulence dispersion model. Thermal and mechanical non-equilibrium between the gas and different classes of particles is allowed. A locally one dimensional, iteratively based, radiative transfer solver is developed to compute wall fluxes and radiative source terms. It is shown that the thermal boundary layer attenuates significantly the radiative fluxes coming from the outer regions. Particle radiation is found to be much more important than gas radiation. Turbulent dispersion of particles in the boundary layer induces a decrease of particle concentration in the region of maximum turbulent kinetic energy, and then, decreases the attenuation effect of wall fluxes due to the boundary layer. The effects of turbulent temperature fluctuations are found to be small in the problem under consideration.