基于米氏散射的气溶胶激光雷达比研究.docxVIP

基于米氏散射的气溶胶激光雷达比研究 1 气溶胶分布的反演方法 气溶胶垂直分布的特点是该局的气候和气候条件非常重要。同时，气溶胶的垂直分布变化往往是许多科学问题的定量化研究，包括气候、云物理、卫星气候、地球化学循环和相对湿度。空气辐射强制性、土壤辐射预算和分层中心的研究包括气溶胶的垂直分布信息。因此,对大气气溶胶消光系数垂直分布的高精度探测,对于全球气候变化的研究有着重要意义。 微脉冲激光雷达是测探气溶胶垂直分布的主要手段,目前已建成全球微脉冲激光雷达探测网(MPLNET)。但在激光雷达探测气溶胶垂直分布反演方法方面,尤其是白天高精度反演方法方面的研究还不够成熟。主要表现在一个激光雷达探测方程中存在两个未知量,即消光系数和后向散射系数,定量反演这两个量必须有其他信息作为辅助[3~6]。若没有辅助信息而使用激光雷达观测数据精确地反演气溶胶消光系数垂直分布,气溶胶的消光系数和后向散射系数的比值,消光后向散射比[也称为激光雷达比(LR)]是一个必需的参数。一般假设该参数为常数,将消光系数和后向散射系数这两个未知参数转化为一个未知量,从而达到求解的目的。 由于LR是时间和空间的复杂函数,它不仅与气溶胶粒子的尺度谱有关,还跟气溶胶的组成有紧密关系,因此,可以通过反演不同区域气溶胶的LR来分析气溶胶的类型。有研究指出在气溶胶消光系数的反演过程中,LR的误差对结果的精确性有很大的影响,在气溶胶的分布不均匀的情况下,这种影响更大。在不同的地理位置和气象条件下反演消光系数时,必须准确地给出对应的LR。另一方面,由于LR与气溶胶的分布类型有着密切的关系,可以通过获得不同区域气溶胶的LR并配合气溶胶的Angstrom波长指数来分析气溶胶的类型。 目前得到LR的方法主要有两大类:独立测量法和联合反演法。独立测量分为多角度测量法、高光谱分辨率雷达法和拉曼雷达法,目前广泛使用拉曼雷达法。联合反演包括用卫星测得的光学厚度限制激光雷达测得的气溶胶消光系数的垂直积分法和利用太阳光度计与激光雷达反演LR。考虑白天情况下拉曼雷达探测精度的限制,本文采取利用太阳光度计反演的方法对中国部分典型地区的LR进行反演与分析。 2 反演方法、数据和误差分析 2.1 气溶胶谱分布的反演算法 激光雷达向大气中发射的激光脉冲由于大气分子和气溶胶粒子的作用向空间各个方向散射,其后向散射部分被激光雷达接收系统接收用于分析气溶胶粒子的物理和光学性质。设β是气溶胶后向散射系数、σe是气溶胶消光系数,则消光系数和后向散射系数的比值就是LR(单位:sr),用S表示,即 单个气溶胶粒子对光的散射与消光作用由其尺度参数x和复折射指数m决定,其中x=2πr/λ,r为粒子半径,λ为光的波长。粒子系统整体光学特性由单个粒子光学特性按粒子谱分布集合决定: 式中积分限rmin、rmax分别为粒子系统的最小和最大特征半径,分别取0.01μm和15μm;Qe、Qπ分别为消光和后向散射效率因子,可利用复折射指数和粒子半径根据米氏散射理论得到;为气溶胶粒子数谱分布,N(r)为与气溶胶半径有关的气溶胶个数。太阳光度计可以反演得到气溶胶体积谱分布,可以表示为 式中下标C,F分别代表气溶胶监测网(AERONET)反演出的气溶胶粗粒子模态和细粒子模态,δ为标准偏差,Cv为体积浓度,rv为体积中值半径,V(r)为与气溶胶半径有关的气溶胶体积。 由气溶胶体积谱分布ddV(lnrr)来反演气溶胶数谱分布,两者之间的关系可以表示为 由(1)~(5)式可以算出LR。选择计算550nm波段的LR,主要原因是国际气象组织选择550nm波段作为研究气溶胶标准,且气溶胶的Angstrom波长指数在550nm与532nm相差很少,LR的误差可以忽略。Angstrom给出了气溶胶光学厚度(AOD)与波长之间的关系式: 式中τ(λ)为对应波段的AOD,k为Angstrom混浊系数,可以代表大气中气溶胶的浓度;α为Angstrom波长指数。(6)式代入不同波长且相除后取对数,即可计算出Angstrom波长指数α: 对于计算Angstrom波长指数,选取的波段很重要,雷达常见波段532nm和1064nm反映出的信息量比较少。本文采用440nm和870nm波段来得到Angstrom波长指数。光学厚度是消光系数的积分,可得 2.2 aerenet站位数据的应用 AERONET是全球布站的气溶胶特性地基观测网,它以法国CIMEL公司生产的全自动CE318太阳光度计为观测仪器,该仪器扫描获得的太阳直射辐射数据可以用于反演气溶胶光学厚度、大气可降水量等。当其进行天空扫描时,所获得的数据可用于反演气溶胶粒子谱分布、复折射指数及气溶胶散射相函数。 利用的AERONET站点数据是榆林(35.94°,104.13°)、北京(39.97°,116.38°)、寿县(116