基于组合拟合法的冰晶粒子的光散射计算

利用组合拟合法计算了冰晶粒子的单次散射特性。给出了消光效率因子、单次散射反照度及非对称因子的拟合公式,利用拟合公式对有效粒子尺度为20μm和120μm的六种冰晶粒子的消光效率因子、单次散射反照度及非对称因子进行了计算。结果表明,粒子的消光效率因子、单次散射反照率和非对称因子随着入射波长的增加有着较大的起伏,后两者随着波长的增加而变化趋势基本一致;对于单次散射反照率来说,在可见光波段,反照率非常接近于1;在短波段,粒子的非对称因子变化较小,并且随着波长的增加,非对称因子会逐渐增大。     

随机取向非球形沙尘气溶胶粒子的光学特性

利用离散偶极子近似法分析了一种随机取向旋转椭球体沙尘气溶胶粒子模型在尺度参数变化范围为0.1~23时(波长0.55!m对应有效半径为0.01~2!m)的光学特性,研究了沙尘粒子非球形性程度对其光学特性的影响,并考察了非球形粒子的随机取向能否用等体积球体来代替。就随机取向单分散和多分散旋转椭球体沙尘气溶胶而言,粒子非球形特征越明显,消光效率因子、不对称因子和单次散射反照率基本上偏离其等体积球体越大;对于相同的非球形,不对称因子偏离其等体积球体的相对偏差要比消光效率因子和单次散射反照率要大。非球形粒子的随机取向并不能使其光学特性严格等效为其等体积球体的光学特性。如果粒子形状偏离球体较小,则非球形粒子的随机取向的平均效果能使其消光效率因子、不对称因子和单次散射反照率近似用等体积球体的对应光学参量来等效;而如果粒子形状偏离球形较大,仅有单次散射反照率可以近似用等体积球体的单次散射反照率来等效,例如,轴半径比为16的旋转椭球体沙尘粒子的单次散射反照率偏离其等体积球体仅在3%以内。     

随机取向双层椭球粒子偏振散射特性研究

基于T矩阵方法,给出了随机取向、轴对称、含核椭球粒子的散射计算方法.散射体的核和外壳均可为非球形粒子,整个粒子具有轴对称性.以含核椭球粒子为模型,计算了含有吸收性内核(黑炭,black carbon)的水凝物气溶胶的散射特性,分析了核的大小、形状对消光系数、散射系数、吸收系数、不对称因子、单次散射反照率以及Muller矩阵等的影响. 关键词： 光散射 T矩阵 Muller矩阵 椭球粒子     

激光信号在沙尘天气下的脉冲时延和展宽

将修正的TTHG(two term Henyey-Greenstein)散射相位函数运用到电子散射理论中,并利用小角度近似条件推导出沙尘天气下激光信号脉冲时延和展宽的表达式。研究了脉冲时延、展宽与沙尘能见度之间的关系,以及散射反照率、散射系数和不对称因子等因素对脉冲时延和展宽的影响;分析了不同波长条件下脉冲时延和展宽随传输距离的变化规律。结果表明,在能见度较低的沙尘天气下,脉冲的时延量和展宽量均随传输距离、反照率和散射系数的增大而增大,随大气能见度的增大而减小。尤其是当大气能见度增加至3.5km时,脉冲时延和展宽逐渐趋于一个与激光波长有关的稳定值。随着激光信号波长和沙粒不对称因子的增大,脉冲的时延和展宽呈负指数趋势减小。     

深圳城市气溶胶物理光学特性的观测研究

气溶腔是影响气候变化和空气污染的重要因子,在深圳地区展开气溶胶观测实验,可以获得可靠的光学物理特征,进而有助于准确评估气溶胶在新型超极城市区域的气像和环境效应,本文利用2010年12月至2011年8月太阳光度计、黑碳和浊度计等气溶胶观测资料,分析了新型超级城市深圳地区的气溶胶物理光学特性,深圳地区气溶胶呈明显季节变化,冬、春季由于城市污染性气溶胶的影响,气溶胶光学厚度和Angstrom波长指数都较大,夏季受海盐气溶胶的影响,光学厚度较小Angstrom波长指数也较小,光学厚度与Angstrom波长指数对比表明城市综合性污染是引起深圳气溶胶高光学厚度的主要原因,深圳地区气溶胶的散射系数、吸收系数的平均值(标准偏差)分别为178.7×10-6 m-1 (126.6×10-6m-1和32.5×1O-6m 1 (18.1×10-6m1),均低于珠三角腹地多年观测平均值的二分之一和国内其他大型城市观测值,而单次散射反照率为0.81,与珠三角其他地区得到的结果接近.此外,气溶胶吸收、散射和单次散射反照率呈明显日变化,可能主要受大气边界层变化的影响.     

激光在云雾中传输的消光特性（英文）

研究了大气分子和云雾粒子对入射激光的吸收与散射效应.采用Kruse模型、单次散射近似法和蒙特卡洛法计算了谱分布参量不同的三种云雾对0.86μm波长激光的消光特性,结果表明:能见度小于0.5km时,Kruse模型不再适用于云雾消光系数的计算;能见度相同时,消光系数的单次散射近似结果大于蒙特卡洛的仿真结果,导致传输距离相同时,云雾对0.86μm波长激光透过率的单次散射近似结果小于蒙特卡洛的仿真结果;同一能见度下谱分布参量不同的三种云雾,消光特性的差别较小,云雾粒子尺度分布对消光特性的影响不明显.     

相态对冰水混合云微观特性反演的影响研究

当云层的温度在－40℃~0℃之间时, 云层中会存在冰和水两种相态的云滴, 其散射特性与纯水云以及纯冰云特性有较大差异, 因此遥感反演混合相云层的微观和宏观物理特性具有重要的意义。本文采用冰水双层球模型模拟了冰水混合云中的云滴, 利用Mie理论计算了纯水、纯冰和冰水颗粒的单次散射特性, 分析了单次散射相函数, 不对称因子, 单次散射反照率等随着有效半径、相态、内外半径比等的变化特性。利用离散纵标法(DISORT)计算了水云和冰云对0.75 μm、2.16 μm和3.3 μm的双向反射函数, 讨论了利用纯水滴和纯冰滴反演冰水混合云滴的误差。分析结果表明, 利用0.75 μm和2.16 μm的太阳光反演冰水混合云的光学厚度和有效半径时, 光学厚度误差较大, 有效半径误差较小; 结合0.75 μm和3.3 μm的太阳光反演冰水混合云的光学厚度和有效半径时, 光学厚度误差较小, 有效半径误差较大, 其会高估其有效半径; 另外结合0.75 μm和3.3 μm这两个波长的反射函数反演冰水云的冰水混合比更为有效。     

煤粉颗粒的后向散射特性研究

利用Mie理论研究了单分散煤粉颗粒对光波的散射特性,给出了单次散射相函数、消光效率因子、散射效率因子、不对称因子和单次散射反照率随粒子尺寸的变化关系;当煤粉浓度较大时,需要考虑颗粒系的多次散射特性,根据辐射传输理论,利用蒙特卡洛方法研究了煤粉颗粒的多次散射特性,并给出了其后向散射随观测角度、粒子半径和光学厚度的变化关系,计算结果对煤粉颗粒反演具有一定指导意义。     

基于Monte Carlo方法的雾红外传输仿真及分析

为了研究多次散射效果下红外辐射在雾中的传输衰减特性,基于Mie散射理论,计算了1.064μm、3.8μm和10.6μm激光在平流雾和辐射雾中的消光参量和多粒子散射相位函数.利用Monte Carlo法建立了辐射传输模型,分析了多次散射效应下接收屏、能见度和传输距离对透过率的影响,并与Lambert-Beer定律计算结果进行了比较.结果表明:当能见度和传输距离均为1km时,1.064μm激光在平流雾中的接收屏粒子数随接收屏边长的增加而显著增加;相同的能见度和传输距离下,激光在辐射雾中的衰减小于在平流雾中的衰减;10.6μm激光在平流雾中具有很好的传输性;当粒子散射能力越强、前向散射概率越大时,多次散射对透过率的贡献越明显;激光在雾中的传输衰减特性不仅与消光系数有关,还与散射系数密切相关;能见度和传输距离存在最佳组合,使得此传输条件下多次散射对透过率的贡献最大.     

沙尘气溶胶粒子群的散射和偏振特性

根据Mie散射理论,以对数正态分布函数描述沙尘气溶胶粒子群的粒径尺度分布,计算了沙尘气溶胶粒子群在0.2～40μm波段间对太阳短波辐射和地球大气长波辐射的单次散射反照率、散射相矩阵函数,揭示了不同相对湿度时,沙尘粒子群对入射辐射的散射和偏振的特征。结果表明,沙尘粒子群的单次散射反照率随着入射波长的增加有较大起伏,不同相对湿度条件下,变化趋势基本一致;在可见光、近红外波段单次散射反照率随湿度增加而变大,湿度95%时非常接近于1;大于10μm的热红外波段单次散射反照率随相对湿度增加而减小,具有较强的吸收辐射能力。散射辐射强度受湿度影响较小,随散射角的增加呈现先减小后增大的趋势,且增大的趋势随着波长的增加而减弱;不同波段上,线偏振和圆偏振随散射角和相对湿度变化存在差异;在前向和后向仅对入射辐射为圆偏振辐射产生圆偏振散射;散射光的偏振特性及其湿度差异主要表现在后向散射区,多以拱形形式体现。拱顶峰值散射角位置存在差异,且峰值散射角随相对湿度的降低向后向漂移。     

Feasibility of a Lidar Utilizing the Glory for Measuring Particle Size of Water Clouds

A new lidar method for measuring water cloud particle size is proposed, and the feasibility of the measurement is discussed. The method utilizes the phenomenon known as the glory which is observed in open air. The proposed lidar consists of a multicolor laser transmitter and two receiver systems looking at the scattering from the target cloud with different scattering angles. Results of the theoretical study show that a system with five laser wavelengths (355, 532, 750, 1064 and 1500 nm) and two receivers located at scattering angles of 180 and 177.5–179 deg is useful for measuring particle size (mode radius of the size distribution) in a range of 4 to 12μm.     

Two-Color Dual-Polarization Pulsed Bistatic Lidar for Measuring Water Cloud Droplet Size

The concept of a pulsed bistatic lidar for measuring water cloud particle size is presented. The method uses a two-color laser and a receiver with a polarization analyzer located at a suitable scattering angle. The dependence of Mie scattering on scattering angle, wavelength, and polarization is used to derive water cloud droplet size. The measurement was simulated for the C₁ and C₂ clouds, and the technique for determining mode radius was studied. The result shows the lidar system with a two-wavelength laser (1064 nm and 532 nm) and a dual-polarization receiver fixed at a scattering angle of around 178 deg can be used to measure a cloud particle size (mode radius) of 4 to 12 μm. Evaluation of the effect of multiple scattering showed that the method can be applied not only for the measurement at the cloud base but also in the cloud where multiple scattering is not negligible.     

Research on bimodal particle extinction coefficient during Brownian coagulation and condensation for the entire particle size regime

The extinction coefficient of atmospheric aerosol particles influences the earth’s radiation balance directly or indirectly, and it can be determined by the scattering and absorption characteristics of aerosol particles. The problem of estimating the change of extinction coefficient due to time evolution of bimodal particle size distribution is studied, and two improved methods for calculating the Brownian coagulation coefficient and the condensation growth rate are proposed, respectively. Through the improved method based on Otto kernel, the Brownian coagulation coefficient can be expressed simply in powers of particle volume for the entire particle size regime based on the fitted polynomials of the mean enhancement function. Meanwhile, the improved method based on Fuchs–Sutugin kernel is developed to obtain the condensation growth rate for the entire particle size regime. And then, the change of the overall extinction coefficient of bimodal distributions undergoing Brownian coagulation and condensation can be estimated comprehensively for the entire particle size regime. Simulation experiments indicate that the extinction coefficients obtained with the improved methods coincide fairly well with the true values, which provide a simple, reliable, and general method to estimate the change of extinction coefficient for the entire particle size regime during the bimodal particle dynamic processes.     

双光路DOAS系统获取大气颗粒物粒谱分布方法研究

当前大气复合污染日趋严重,造成大气氧化性增强,气体向颗粒物的转化加快。大气颗粒物粒径大小及谱分布决定其在大气中的行为,以差分吸收光谱法（DOAS）为基础,结合双光路设计技术,开展实时、在线、获取近地面大气气溶胶颗粒物的粒谱分布的光谱方法研究。首先构建低噪声性能稳定的宽带氙弧灯为光源的双光路差分吸收光谱系统,基于干净天气条件下大气的能见度数据对系统进行校准,通过两个不同光路获得的光谱信号强度之比获取近地面紫外-近红外波段的大气总宽波段消光系数。基于宽波段消光系数,在去除瑞利散射以及气体吸收对消光系数的影响后,解析出气溶胶颗粒物的消光系数。基于核函数准则,利用均匀球型粒子的电磁场Mie理论来反演气溶胶物理特性,获得气溶胶粒子在该测量谱段的体积谱分布,利用体积谱与数密度谱的关系,反演出气溶胶粒子的数密度谱分布。开展利用直方图方法来表现颗粒物的粒谱分布方法研究,首先将DOAS测量波段近似等分为若干谱段,利用谱段处平均值,获取气溶胶粒谱直方分布图。最后把该系统和方法应用于外场实验,获得了气溶胶颗粒物在300～650 nm范围内的消光系数,将测量波段等分为11个谱段,反演了颗粒物的在0.1～1.25 μm粒径范围的数密度谱分布。该研究为整治我国灰霾天气,研究大气气相/粒子非均相化学反应提供科学依据。同时将推动DOAS技术的进一步发展和应用。     

Radiative properties of cirrus clouds in the infrared (8–13 μm) spectral region

Atmospheric radiation in the infrared (IR) 8–13 μm spectral region contains a wealth of information that is very useful for the retrieval of ice cloud properties from aircraft or space-borne measurements. To provide the scattering and absorption properties of nonspherical ice crystals that are fundamental to the IR retrieval implementation, we use the finite-difference time-domain (FDTD) method to solve for the extinction efficiency, single-scattering albedo, and the asymmetry parameter of the phase function for ice crystals smaller than 40 μm. For particles larger than this size, the improved geometric optics method (IGOM) can be employed to calculate the asymmetry parameter with an acceptable accuracy, provided that we properly account for the inhomogeneity of the refracted wave due to strong absorption inside the ice particle. A combination of the results computed from the two methods provides the asymmetry parameter for the entire practical range of particle sizes between 1 and 10,000 μm over the wavelengths ranging from 8 to 13 μm. For the extinction and absorption efficiency calculations, several methods including the IGOM, Mie solution for equivalent spheres (MSFES), and the anomalous diffraction theory (ADT) can lead to a substantial discontinuity in comparison with the FDTD solutions for particle sizes on the order of 40 μm. To overcome this difficulty, we have developed a novel approach called the stretched scattering potential method (SSPM). For the IR 8–13 μm spectral region, we show that SSPM is a more accurate approximation than ADT, MSFES, and IGOM. The SSPM solution can be further refined numerically. Through a combination of the FDTD and SSPM, the extinction and absorption efficiencies are computed for hexagonal ice crystals with sizes ranging from 1 to 10,000 μm at 12 wavelengths between 8 and 13 μm.

Calculations of the cirrus bulk scattering and absorption properties are performed for 30 size distributions obtained from various field campaigns for midlatitude and tropical cirrus cloud systems. Ice crystals are assumed to be hexagonal columns randomly oriented in space. The bulk scattering properties are parameterized through the use of second-order polynomial functions for the extinction efficiency and the single-scattering albedo and a power-law expression for the asymmetry parameter. We note that the volume-normalized extinction coefficient can be separated into two parts: one is inversely proportional to effective size and is independent of wavelength, and the other is the wavelength-dependent effective extinction efficiency. Unlike conventional parameterization efforts, the present parameterization scheme is more accurate because only the latter part of the volume-normalized extinction coefficient is approximated in terms of an analytical expression. After averaging over size distribution, the single-scattering albedo is shown to decrease with an increase in effective size for wavelengths shorter than 10.0 μm whereas the opposite behavior is observed for longer wavelengths. The variation of the asymmetry parameter as a function of effective size is substantial when the effective size is smaller than 50 μm. For effective sizes larger than 100 μm, the asymmetry parameter approaches its asymptotic value. The results derived in this study can be useful to remote sensing studies of ice clouds involving IR window bands.     

近红外波段气溶胶的衰减特性研究

在已知大气气溶胶折射率和气溶胶谱分布的基础上, 对近红外波段的气溶胶衰减特性进行了研究。利用Mie散射理论计算并讨论了气溶胶的消光、散射、吸收效率因子随尺度参数的变化和消光系数随半径和波长的变化, 并且在MATLAB中对各种变化情况进行了仿真。结果表明, 三种气溶胶粒子的消光和散射能力依次为沙尘性粒子, 水溶性粒子, 烟煤。消光系数在粒子半径和入射波长相近时达到最大, 并且粒子半径对消光、散射、吸收系数的影响比入射波长更明显。这些结论可以为红外辐射在大气中的衰减计算和分析提供依据。     

颗粒系的可见消光光谱分析及最佳波长的选择

在光全散射法颗粒粒径测量中,被测颗粒系的消光光谱包含有颗粒粒径、折射率等信息。对常用的单峰及双峰R-R分布的可见消光光谱进行模拟,分析了可见消光光谱随粒径和相对折射率变化的规律。选取二阶微分可见消光光谱不连续点对应的波长作为测量波长,并将可见光边界波长同时作为测量波长,在非独立模式下,采用遗传优化算法反演粒径分布。仿真计算与实验验证结果表明,通过消光光谱分析,可以预先确定被测颗粒系的分布状况,缩小优化算法中反演参数的搜索范围。结合最佳波长选择方法,使测量的精度和可靠性明显提高。对测量光谱加入1%随机噪声时,单峰及多峰分布粒径反演均能得到满意的结果。     

Optical depth and multiple scattering depolarization in liquid clouds

In this paper, we calculate multiply scattered lidar signals with Monte Carlo method for measuring optical depth (extinction coefficient), effective size of water droplets, and liquid water content of clouds, and present algorithms that implement our method. We calculated multiply scattered lidar signals for various water droplet sizes and liquid water contents using a Monte Carlo method. A simple correspondence between water droplet optical depth and the degree of polarization in a modified gamma size distribution (C₁ cloud) is found. We also calculated the degree of polarization of a lidar signal for a given liquid water content, finding that the degree of polarization is only dependent on optical depth. Since the Raman lidar signal of liquid water depends on the total volume of the water droplet, the effective radius of the water droplet can thus be recovered from the degree of polarization of the lidar signal and the Raman signal of the liquid water.     

作用于水云粒子的太阳辐射压力

根据球形粒子的Mie散射理论和平行平面大气中的辐射传输理论分析了作用于水云粒子的太阳辐射压力。结果表明,一般情况下当云滴半径不大于10μm时,太阳辐射压力可达云滴重力的百分之几;云滴半径在0.3μm附近时,太阳辐射压力达到极大值,约为云滴重力的百分之几十。因此在云物理研究中应考虑作用于云粒子的太阳辐射压力。     

The impact of ice cloud particle microphysics on the uncertainty of ice water content retrievals

Ice water content (IWC) is a standard product of cloud radar measurements. In this work, cloud radar cross-sections of various ice clouds are modeled to examine the relationship between the radar signal and the IWC. We report that using backscatter signal at cloud radar wavelength to retrieve IWC results in large uncertainties. Particle size distribution is the primary cause for the uncertainty in the retrieved IWC at radar wavelengths, though particle shape and orientation also play significant roles. Particularly in this study, we demonstrate that using both transmitted waves through the clouds (extinction) and backscattered waves from the clouds to retrieve the mean particle size and then using the mean particle size for IWC retrieval reduces the uncertainty. IWC retrieval can be improved with size distribution derived from dual wavelength cloud radar.