三波长激光雷达拟合卷云的粒子形状

卷云在全球出现的概率可达30%,其散射特性在气候模式、光辐射传输和遥感领域都具有非常重要的意义。卷云的散射特性主要由冰晶粒子形状、尺度谱、折射率等因素所决定。利用355, 532和1 064 nm三个波长激光雷达数据反演卷云的后向散射系数颜色比,利用模拟计算获得不同形状冰晶粒子的卷云在上述三个波长上的后向散射系数颜色比,通过拟合得出被测卷云的冰晶粒子形状。拟合结果表明,合肥上空卷云中冰晶粒子大部分可能呈聚合物状。     

A simple new radiative transfer model for simulating the effect of cirrus clouds in the microwave spectral region

The upwelling atmospheric radiation in the millimeter wave spectral range is influenced by the presence of cirrus clouds. A plane parallel radiative transfer model which can take into account the effect of multiple scattering by ice particles in the cirrus has been developed and is used to simulate the brightness temperatures as they would be measured by a satellite instrument. The model uses an iterative procedure to solve the radiative transfer equation. The formulation of the model is such that it can easily be adapted to treat the full specific intensity vector instead of just the scalar total intensity. A convergence test for the model is explained and two cirrus cloud scenarios are simulated. The results illustrate the linearity of microwave radiative transfer for not too strong cirrus clouds in this frequency region.     

Sensitivity of depolarized lidar signals to cloud and aerosol particle properties

Measurements from depolarized lidars provide a promising method to retrieve both cloud and aerosol properties and a versatile complement to passive satellite-based sensors. For lidar observations of clouds and aerosols, multiple scattering plays an important role in the scattering process. Monte Carlo simulations are carried out to investigate the sensitivity of lidar backscattering depolarization to cloud and aerosol properties. Lidar parameters are chosen to be similar to those of the upcoming space-based CALIPSO lidar. Cases are considered that consist of a single cloud or aerosol layer, as well as a case in which cirrus clouds overlay different types of aerosols. It is demonstrated that besides thermodynamic cloud phase, the depolarized lidar signal may provide additional information on ice or aerosol particle shapes. However, our results show little sensitivity to ice or aerosol particle sizes. Additionally, for the case of multiple but overlapping layers involving both clouds and aerosols, the depolarized lidar contains information that can help identify the particle properties of each layer.     

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.     

激光在准球形边界卷云中传输的衰减和透过特性

高空卷云主要由各种不同形状的冰晶粒子组成,是地空链路上激光信号传输的重要影响因素。依据高空卷云中冰晶粒子的分布特征和散射特性,采用C版本的离散纵标法（CDISORT）,充分考虑地球球形曲率及云层冰晶粒子多次散射影响因素,研究准球形边界云层的激光透过率和衰减特性,并比较了太阳天顶角不同时平面平行模式和准球面模式下卷云大气激光透过率的差异,数值计算了三种激光波长（0.65,1.06和3.8 μm）在卷云中传输时的衰减和透过特性。计算结果表明：较小太阳天顶角（小于80°）入射时,两种模式下卷云大气激光透过率相对误差很小,其中0.65 μm激光波长入射时两种模式下的相对误差仅为1.72%,较大太阳天顶角（大于80°）入射时,两种模式下卷云大气激光透过率相对误差明显增大,0.65 μm激光波长入射时两种模式下的相对误差最大达到69%;卷云粒子单次散射时,激光在云层的衰减与卷云粒子有效半径、传输距离、光学厚度及激光波长等因素有关,随光学厚度的增加,云层的激光透过率减少,1.06 μm激光波长入射时透过率最大,3.8 μm激光波长入射时透过率最小;0.65和1.06 μm激光波长入射时,随云层粒子有效半径的增加激光透过率逐渐增加,而3.8 μm波长激光,随云层粒子有效半径的增加激光透过率逐渐减少,随相对方位角的增加,云层的激光透过率减少,且不同卷云传输模型对激光透过率也存在不同的影响。该研究工作将为开展地空链路星载、机载激光通信、激光雷达探测等工程系统中的激光信号云层传输特性的应用提供理论支持,同时也可进一步拓展为地空链路激光遥感、制导和预警等应用提供预先理论研究基础。     

Photon conservation in scattering by large ice crystals with the SASKTRAN radiative transfer model

The scattering of visible light by ice crystals and dust in radiative transfer models is challenging in part due to the large amount of scattering in the forward direction. We introduce a technique that ensures numerical conservation of photons in any radiative transfer model and that quantifies the integration error associated with highly asymmetric phase functions. When applied to a successive-orders of scatter model, the technique illustrates the high accuracy obtained in numerical integration of molecular and aerosol scattering. As well, a phase function truncation and renormalization technique is applied to scattering by ice crystals with very large size parameters, between 100 and 1000, and the scaled radiative transfer equation is solved with the spherical successive-orders model, SASKTRAN. Since computations shown this work are performed in a fully spherical model atmosphere, the computed radiances are not subject to the discontinuity at the horizon that is inherent in models using a plane–parallel assumption. The methods introduced in this work are of particular interest in modeling limb radiances in the presence of thin cirrus clouds.     

An algorithm to determine cirrus properties from analysis of multiple-scattering influence on lidar signals

A Monte Carlo-based evaluation of the multiple-scattering influence on ground-based Raman lidar measurements is presented. The lidar returns from cirrus clouds are analyzed in order to evaluate vertical profiles of the extinction and backscattering coefficients. Results show that for the typical cirrus cloud, the presence of the multiple scattering can lead to an underestimation of the extinction coefficient by as large as 200% whereas the backscattering coefficient is almost unaffected for the Raman lidar technique. An algorithm to select one or a set of phase functions which fit to the lidar data is also presented. It is an iterative procedure based on Monte Carlo scattering simulation. By comparison of the experimental value of the lidar ratio, corrected for the multiple scattering influence, and the phase function used in the Monte Carlo simulation, one can determine a suitable phase function. The validity and sensitivity of the algorithm have been demonstrated by applying it to simulated cases. The application to some real cases indicates that our procedure allows for the establishing of a practical scattering model for the cirrus clouds.     

A Monte Carlo method for 3D thermal infrared radiative transfer

A 3D Monte Carlo model for specific application to the broadband thermal radiative transfer has been developed in which the emissivities for gases and cloud particles are parameterized by using a single cubic element as the building block in 3D space. For spectral integration in the thermal infrared, the correlated k-distribution method has been used for the sorting of gaseous absorption lines in multiple-scattering atmospheres involving 3D clouds. To check the Monte-Carlo simulation, we compare a variety of 1D broadband atmospheric fluxes and heating rates to those computed from the conventional plane-parallel (PP) model and demonstrate excellent agreement between the two. Comparisons of the Monte Carlo results for broadband thermal cooling rates in 3D clouds to those computed from the delta-diffusion approximation for 3D radiative transfer and the independent pixel-by-pixel approximation are subsequently carried out to understand the relative merits of these approaches.     

复杂大气背景下机载通信终端与无人机目标之间的激光传输特性研究

云层、气溶胶和大气分子是大气环境的主要组成部分.本文基于逐次散射法求解辐射传输方程,建立了复杂大气背景下机载无线光通信终端与地空无人机目标之间的激光传输模型.考虑真实大气背景中卷云、大气分子和气溶胶存在的情况下,数值计算了1.55 mm激光经机载通信终端发出后通过大气背景的直接传输和一阶散射传输后接收功率随无人机目标高度的变化关系,分析了飞机在云上、云中和云下以及卷云冰晶粒子有效半径、飞机与无人机之间的水平距离对接收激光信号传输功率的影响.数值结果表明:激光通过卷云传输的功率很大程度上取决于飞机在云上、云下或云中的位置;飞机与无人机目标之间的水平距离和卷云冰晶粒子的有效半径对激光直接传输和一阶散射传输影响较大;与云上大气相比,云下的大气分子和气溶胶对激光有较大的衰减.本文工作可为进一步开展地空链路上复杂大气背景对机载与低空无人机目标激光通信实验、无人机编队、指挥和组网技术的研究提供理论支撑.     

Retrieval of tropospheric NO2 columns from satellite measurements in presence of cirrus: A theoretical sensitivity study using SCIATRAN and prospect application for the A-Train

A theoretical sensitivity study of the influence of cirrus cloud properties on tropospheric NO₂ columns retrieved from the spaceborne Ozone Monitoring Instrument (OMI) measurements is performed. It is conducted within the framework of the synergetic use of A-Train sensors to derive more representative trace gas products. We aim to study the potential effects of cirrus clouds on tropospheric NO₂ retrievals using a retrieval algorithm that, unlike the OMI Standard and DOMINO algorithms, does not correct for the effects of clouds. The sensitivity study is based on the radiative transfer code SCIATRAN that performs both simulations of top of atmosphere (TOA) reflectances as measured by an OMI-like band and tropospheric NO₂ column retrievals based on the differential optical absorption spectroscopy (DOAS) method. The results of the sensitivity study show that if a correction for cirrus clouds is not included in our simple retrieval that does not account for clouds in the first place, the tropospheric column can be underestimated by 55%. This underestimation depends strongly on cirrus parameters as, in order of importance, cloud fraction, cloud optical depth, asymmetry factor of cirrus cloud phase function and cloud top height. The perspective of the synergy between OMI and cloud information obtained from cloud-derived products of the A-Train is evaluated in two parts by applying a simple cloud correction scheme based on the independent pixel approximation (IPA). Firstly, we evaluated the tropospheric NO₂ column retrievals error caused by uncertainties in cirrus cloud properties. Secondly we studied the influence of subpixel cloud optical depth variability on NO₂ retrievals. From our simulations, it is demonstrated that the error will be reduced significantly if the cloud fraction is lower or equal to 0.5. In this case, the cloud fraction and the cloud optical depth must be known within accuracy less than 0.05% and 50%, respectively. The cloud top height and the asymmetry factor must be known within uncertainty of at least 1 km and less than 0.05, respectively. The latter result shows that the uncertainty of the asymmetry factor is a major source of error in the cloud correction for tropospheric NO₂ retrieval in the presence of cirrus.     

高斯光束在双层云中传输的蒙特卡罗模拟

基于辐射传输理论, 利用蒙特卡罗方法模拟了无限窄(冲击函数)准直光束入射到典型水云以及冰水双层云时的后向散射特性, 进而将得到的冲击响应与高斯光束卷积, 得到高斯光束在云层中传输的多次散射特性. 文中给出了两种波束入射时水云以及冰水双层云的反射函数随径向r和天顶角α的变化关系, 并给出了光强在云层内部的二维分布图. 计算结果表明, 高斯光束入射时, 云层反射函数的特点与无限窄准直光束入射时有较大区别. 因此在利用激光雷达进行云层探测时需要考虑激光的散斑, 文中的方法可以为此提供理论依据.     

Validation of the community radiative transfer model

To validate the Community Radiative Transfer Model (CRTM) developed by the U.S. Joint Center for Satellite Data Assimilation (JCSDA), the discrete ordinate radiative transfer (DISORT) model and the line-by-line radiative transfer model (LBLRTM) are combined in order to provide a reference benchmark. Compared with the benchmark, the CRTM appears quite accurate for both clear sky and ice cloud radiance simulations with RMS errors below 0.2 K, except for clouds with small ice particles. In a computer CPU run time comparison, the CRTM is faster than DISORT by approximately two orders of magnitude. Using the operational MODIS cloud products and the European Center for Medium-range Weather Forecasting (ECMWF) atmospheric profiles as an input, the CRTM is employed to simulate the Atmospheric Infrared Sounder (AIRS) radiances. The CRTM simulations are shown to be in reasonably close agreement with the AIRS measurements (the discrepancies are within 2 K in terms of brightness temperature difference). Furthermore, the impact of uncertainties in the input cloud properties and atmospheric profiles on the CRTM simulations has been assessed. The CRTM-based brightness temperatures (BTs) at the top of the atmosphere (TOA), for both thin (τ<5) and thick (τ>30) clouds, are highly sensitive to uncertainties in atmospheric temperature and cloud top pressure. However, for an optically thick cloud, the CRTM-based BTs are not sensitive to the uncertainties of cloud optical thickness, effective particle size, and atmospheric humidity profiles. On the contrary, the uncertainties of the CRTM-based TOA BTs resulting from effective particle size and optical thickness are not negligible in an optically thin cloud.     

Scattering in thick multifractal clouds, Part II: Multiple scattering

In Part I of this paper, we developed asymptotic approximations for single photon scattering in thick, highly heterogeneous, “Log-Lévy” multifractal clouds. In Part II, theoretical multiple scattering predictions are numerically tested using Monte Carlo techniques, which show that, due to long range correlations, the photon paths are “subdiffusive” with the corresponding fractal dimensions tending to increase slowly with mean optical thickness. We develop reasonably accurate statistical relations between N scatter statistics in thick clouds and single scatter statistics in thin clouds. This is explored further using discrete angle radiative transfer (DART) approach in which the radiances decouple into non-interacting families with only four (for 2-D clouds) radiance directions each. Sparse matrix techniques allow for rapid and extremely accurate solutions for the transfer; the accuracy is only limited by the spatial discretization.By “renormalizing” the cloud density, we relate the mean transmission statistics to those of an equivalent homogeneous cloud. This simple idea is remarkably effective because two complicating effects act in contrary directions: the “holes” which lead to increased single scatter transmission and the tendency for multiply scattered photons to become “trapped” in optically dense regions, thus decreasing the overall transmission.     

太赫兹波被动遥感卷云微物理参数的敏感性试验分析

太赫兹波长和典型卷云的冰晶粒子尺度处于同一量级,其在遥感卷云微物理参数(粒子尺度和冰水路径)方面具有广阔的应用前景.为了评估卷云微物理参数对太赫兹波传输特性的影响及其在太赫兹波段的敏感性,基于大气辐射传输模式分别模拟计算了晴空和有云条件下大气层顶的太赫兹辐射光谱特征,分析了这两种条件下辐射亮温差值的特点,研究了卷云冰晶粒子形状、粒子尺度及冰水路径对太赫兹辐射传输特性的影响,并定量计算了相关敏感系数.结果表明:卷云冰晶粒子形状、粒子尺度、冰水路径等对太赫兹波传输特性均有不同程度的影响,卷云效应也因通道频率而异,太赫兹波对卷云的粒子尺度和冰水路径有较高的敏感性,是理论上被动遥感卷云微物理特性的最佳波段.研究结果对于进一步发展太赫兹波被动遥感卷云技术、提高卷云参数的反演精度具有重要意义.     

A strong ice cloud event as seen by a microwave satellite sensor: Simulations and observations

In this article, brightness temperatures observed by channels of the Advanced Microwave Sounding Unit-B (AMSU-B) instrument are compared to those simulated by a radiative transfer model, which can take into account the multiple scattering due to ice particles by using a discrete ordinate iterative solution method. The input fields, namely, the pressure, temperature, humidity, and cloud water content are taken from the short range forecast from the Met Office mesoscale model (UKMES). The comparison was made for a case study on the 25 January 2002 when a frontal system associated with significant cloud was present over the UK. It is demonstrated that liquid clouds have maximum impact on channel 16 of AMSU whereas ice clouds have maximum impact on channel 20. The main uncertainty for simulating microwave radiances is the assumptions about microphysical properties, such as size distribution, shape and orientation of the cloud particles, which are not known in the mesoscale model. The article examines the impact of these parameters on the cloud signal. The polarisation signal due to oriented ice particles at these frequencies is also discussed.     

Simple relation between lidar multiple scattering and depolarization for water clouds

An empirical relationship is derived between the multiple-scattering fraction and the linear depolarization ratio by using Monte Carlo simulations of water clouds measured by backscatter lidar. This relationship is shown to hold for clouds having a wide range of extinction coefficients, mean droplet sizes, and droplet size distribution widths. The relationship is also shown to persist for various instrument fields of view and for measurements made within broken cloud fields. The results obtained from the Monte Carlo simulations are verified by using multiple-field-of-view lidar measurements. For space-based lidars equipped to measure linear depolarization ratios, this new relationship can be used to accurately assess signal perturbations due to multiple scattering within nonprecipitating water clouds.     

Temperature dependence of ice optical constants: Implications for simulating the single-scattering properties of cold ice clouds

For a spectrum from ultraviolet to microwave and a temperature range from 160 to 270 K, the optical constants of water ice are compiled on the basis of the Kramers-Kronig relation in conjunction with existing datasets reported in literature. Significant temperature dependence is observed in both the mid-infrared and longer wavelengths. A sensitivity study at wavelengths in the infrared split window region indicates that the temperature dependence of the single-scattering properties of ice crystals is not negligible. Thus, it is necessary to take into account the temperature dependence of ice optical constants when simulating the radiative properties of cirrus clouds for various applications to remote sensing under cirrus cloud conditions.     

冰水混合云对可见光的吸收和散射特性

根据Mie理论，分别计算了由纯水、纯冰和冰-水同心球形粒子构成的云层在可见光波段的单次散射特性. 根据辐射传输理论，利用叠加法数值计算了这三种不同构成的云层的反射函数，以及它们的平面反照率、透过率和吸收率. 结果表明，冰云和冰_水云的反射函数和平面反照率在大部分散射角下要略小于水云，而透过率却比水云的大. 通过具体的数值计算结果和理论分析，研究了水云的异常吸收现象. 关键词： 冰_水粒子 Mie理论 叠加法 光散射     

An efficient diffusion approximation for 3D radiative transfer parameterization: application to cloudy atmospheres

The three-dimensional (3D) diffusion radiative transfer equation, which utilizes a four-term spherical harmonics expansion for the scattering phase function and intensity, has been efficiently solved by using the full multigrid numerical method. This approach can simulate the transfer of solar and thermal infrared radiation in inhomogeneous cloudy conditions with different boundary conditions and sharp boundary discontinuity. The correlated k-distribution method is used in this model for incorporation of the gaseous absorption in multiple-scattering atmospheres for the calculation of broadband fluxes and heating rates in the solar and infrared spectra. Comparison of the results computed from this approach with those computed from plane-parallel and 3D Monte Carlo models shows excellent agreement. This 3D radiative transfer approach is well suited for radiation parameterization involving 3D and inhomogeneous clouds in climate models.     

The role of cloud-scale resolution on radiative properties of oceanic cumulus clouds

Both individual and combined effects of the horizontal and vertical variability of cumulus clouds on solar radiative transfer are investigated using a two-dimensional (x- and z-directions) cloud radar dataset. This high-resolution dataset of typical fair-weather marine cumulus is derived from ground-based cloud radar observations. The domain-averaged (along x-direction) radiative properties are computed by a Monte Carlo method. It is shown that (i) different cloud-scale resolutions can be used for accurate calculations of the mean absorption, upward and downward fluxes; (ii) the resolution effects can depend strongly on the solar zenith angle; and (iii) a few cloud statistics can be successfully applied for calculating the averaged radiative properties.