气溶胶激光雷达的国内外研究进展与展望

黄忠伟，王雍恺，闭建荣，王天河，李武仁，李泽，周天

1.兰州大学西部生态安全省部共建协同创新中心,兰州730000;

2.兰州大学大气科学学院半干旱气候变化教育部重点实验室,兰州730000

1 引言

气溶胶是大气中的重要组分，对气候、环境、生态系统均有重要的影响（Andreae 和Crutzen，1997；毛节泰等，2002；石广玉等，2008；Bi等，2014；Huang 等，2014）。而要准确评估其气候与环境效应，必须获知它们物理化学光学特性的时空分布及其演变规律。不同地区、不同季节及在不同的高度上，气溶胶性质有着显著的差异，使得气溶胶研究存在很大不确定性（Aloyan 等，1997）。

由于具有探测距离远、相干性好、时空分辨率高等独特技术优势，激光雷达成为气溶胶研究领域的不可或缺手段。激光雷达结合了传统雷达技术和现代激光技术，是常规地面监测技术的有力补充，可探测气溶胶垂直分布及其快速变化（周军等，1998，2000；邹炳芳等，2008；赵一鸣等，2014）。1963年，世界首台红宝石米散射激光雷达诞生于美国（Fiocco 和Smullin，1963）。1966年，中国科学院大气物理研究所与上海光学精密机械研究所共同研制中国第一台大气探测激光雷达（王青梅和张以谟，2006）。随着探测业务需求与光电技术的快速发展，中国多家科研院所也开展了大量的大气探测激光雷达的仪器研发与应用研究，例如，中国科学院安徽光学精密机械研究所、上海光学精密机械研究所、武汉物理与数学研究所，西安理工大学、武汉大学、中国海洋大学、浙江大学、兰州大学等（Gong等，2015；吕达仁等，1977；刘金涛等，2003；程学武等，2007；毛建东等，2010；刘东等，2013；闫庆等，2013；周天等，2013；黄寒璐等，2014；冯长中等，2018；Huang 等，2020），有力推动了中国大气探测激光雷达的发展，并逐步引领国际激光雷达领域的技术前沿。

当前，激光雷达已被广泛应用于地基观测网与卫星观测，获得连续、可靠的区域或全球气溶胶光学物理特性观测数据资料。例如，美国宇航局（NASA）建设以全自动运行的微脉冲米散射激光雷达为主的MPLNET（Welton等，2018），欧洲航空局（ESA）主导建设的以拉曼激光雷达系统为主的欧洲激光雷达观测网EARLINET（Pappalardo等，2014），日本科研机构建设的以偏振激光雷达为主的东亚沙尘观测激光雷达网AD-net（Murayama等，2001；Sugimoto 等，2018），以及全球大气成分变化探测网NDACC（De Mazière 等，2018）、美国东部激光雷达观测网REALM（Hoff等，2005）、独联体激光雷达网CIS-LINET（Chaikovsky 等，2006）等。此外，近年来星载激光雷达的发展也越来越快。1994年，搭载发现号航天飞机升空的LITE 激光雷达是世界上第一台星载激光雷达（McCormick和Hostetler，1995）。2006-04，CALIPSO 卫星搭载全球首个星载云与气溶胶激光雷达（CALIOP）成功投入使用。CALIOP 采用532 nm 和1064 nm 两个波长，拥有偏振探测能力可更准确地区分云和气溶胶类型，为气溶胶—云—气候效应评估提供了全球气溶胶的三维时空分布（Vaughan 等，2004；Winker等，2007）。

激光雷达已成为国内外大气探测的重要手段，在探测气溶胶物理光学特性取得了长足的进步。本文主要总结国内外在气溶胶激光雷达硬件研发、算法开发以及大气遥感与环境监测应用等方面的研究进展。

2 激光雷达基本探测原理

气溶胶激光雷达的主要探测原理是激光器向大气中发射的激光被气溶胶、空气分子散射与吸收后，其中散射回波信号被天文望远镜接收，然后通过精细分光、提取再使用光电探测器转换成电信号，采用高速数据采集卡获得数据并保存，最后进行数据处理与反演就可以得到气溶胶特征信息，如图1所示。通常情况下，激光雷达系统包括激光发射、信号接收探测与数据采集处理等3个子系统。

图1 气溶胶激光雷达系统基本探测原理示意图Fig.1 Schematic diagram of basic detection principle of aerosol lidar system

激光发射子系统包括激光器和扩束器。按照工作物质类型不同，激光器可分为固体、半导体、染料和气体激光器；扩束器可扩大激光光斑，同时使得发散角变小。信号接收子系统包括望远镜、光阑、分光镜、探测器等。常用望远镜有牛顿反射式望远镜、卡塞格林反射式望远镜；光阑放置在望远镜的焦点，可减少天空背景噪声和多次散射效应影响；目前探测器主要使用雪崩二极管（APD），光电倍增管（PMT）以及电荷耦合器件（CCD）等。数据采集子系统包括高速数据采集卡、电脑等。有两种数据采集方式：一是模数转换技术，具有高线性度的特点，常用于采集强回波信号；二是光子计数技术，具有高灵敏度的特点，适合微弱回波信号。通常情况下，同时使用两种采集技术，提高探测信号动态范围，有效提升激光雷达探测能力。常用的气溶胶激光雷达类型，如表1所示。

表1 气溶胶激光雷达主要类型介绍Table 1 Main types of aerosol lidar

2.1 米散射激光雷达

米散射激光雷达是激光雷达中发展的最成熟、应用最广泛的类型。当散射粒子尺寸与激光波长相当的时候发生米散射现象，其特点是散射前后波长没有变化，大气中与激光发生米散射的主要是气溶胶与云。由于其散射截面较大，通常米散射激光雷达的回波信号较强。

在使用激光雷达的数据反演气溶胶消光系数或后向散射系数之前，需要先对原始数据进行一系列订正，依次包括背景噪声、距离平方及几何重叠因子等订正（黄忠伟，2012）。一般选取高空无气溶胶的区域信号平均值作为背景信号，将回波信号减去该平均值即完成背景噪声订正，然后把观测数据乘以对应探测距离的平方即可完成距离平方订正。由于发射场与接收场的交迭并不完全重合，导致在近距离只有部分信号被探测，需对此进行信号修正，因而准确地订正几何重迭因子极为重要。目前常用的几何重叠因子获取方法有解析法、水平实验法及拉曼观测法等。把经过距离订正的信号除以几何重叠因子，就可完成该订正步骤。

利用米散射激光雷达数据反演气溶胶消光系数，需要先假定雷达比（LR），即气溶胶消光系数与后向散射系数之比，由气溶胶组分、大小、形状与湿度等因子决定。不同类型气溶胶的雷达比，如表2所示。

米散射激光雷达方程为

式中，P为激光雷达回波信号强度，C是雷达常数，E是发射激光脉冲能量，O（z）是几何重叠因子，σ与β分别是大气总消光系数和总后向散射系数，Z是探测距离，PB是接收背景信号。目前求解米散射激光雷达方程的方法主要有Fernald 和Klett法。

（1）Fernald 法。假设大气由空气分子和气溶胶两部分组成，于是大气消光系数（后向散射系数）也是由大气分子和气溶胶两部分的消光系数（后向射射系数）共同组成。因此，气溶胶后向散射系数反演方程为（Fernald，1984；Welton和Campbell，2002）

（2）Klett 法。假设气溶胶消光系数σ与后向散射系数β存在如下关系（Collis和Russell，1976；Klett，1981）：

式中，C为常数，k受气溶胶物理化学特性和入射激光波长影响，通常情况下k取值范围为0.67—1.0。

将式（3）代入激光雷达方程，解得：

式中，X为已订正处理后的激光雷达数据，Zm往往是最大探测距离或者气溶胶浓度很低的高度。

k值的确定主要通过两个步骤，一是采用不同k值反演气溶胶消光系数随高度的变化值及其垂直分布趋势。一般认为只要有气溶胶粒子存在，消光系数大于零，所以若反演的消光系数存在负值或者随高度升高出现明显振荡，便可缩小k的取值范围。二是通过太阳光度计结果和激光雷达反演光学厚度比较差异性，当差值近似为零时，可精确k值如表3所示（王治华等，2006；张文煜等，2008；张京朋等，2016；孙海波，2017）。

表3 不同天气和k值下消光系数平均值（王治华等，2006）Table 3 Average extinction coefficient under different weather and k value（Wang et al.，2006）

2.2 偏振激光雷达

偏振激光雷达探测原理是发射线偏振激光与气溶胶发生散射，非球形粒子会改变激光偏振方向。当回波信号经过偏振晶体（PBS）后，被分为与发射激光偏振方向平行的P 通道和垂直的S 通道。其中，气溶胶退偏比δ可定义为S通道信号P⊥与P通道信号P||之比，其方程如下：式中，K为偏振订正系数，由两个偏振探测通道的光学和电子学增益差别造成（刘东等，2003）。一般来说，球形粒子的退偏比为零，而非球形粒子的退偏比大于零。

图2 偏振雷达的分光装置图Fig.2 Polarization lidar splitter device diagram

由于两个偏振通道光电效率不一定相等，即K往往不等于1，需要进行有效定标。目前常用的偏振订正方法（李保卫等，2019）有：一是大气分子法，利用大气分子的退偏比作为已知量，再获取某一高度退偏比实际值。该值与大气分子退偏比的比值即为两通道订正系数；二是半波片法，在偏振晶体之前插入半波片，通过半波片使P、S通道发生旋转，比较前后偏振信号的差异进而完成定标（Freudenthaler等，2009）。

2.3 拉曼激光雷达

拉曼散射是激光与大气分子之间的一种非弹性散射过程，在相互作用中散射光因为分子内部结构的振动发生频率漂移，导致散射光波长与入射光波长不一致。目前一般使用短波飘移部分，常用波长为355 nm，532 nm。大气中各种主要气体的拉曼频移量为：氧气（1556 cm-1），氮气（2331 cm-1），水汽（3652 cm-1）。拉曼激光雷达探测气溶胶利用大气中高度分布稳定的氮气分子振动拉曼散射信号只受气溶胶消光影响的特点，在无需激光雷达比、边界条件假设条件下实现对气溶胶光学特性的精细探测，从而减小误差（高飞等，2017）。由于拉曼信号比较弱，比米散射信号要小3—4 个量级，且散射截面小，受背景噪音影响大，不适合白天探测（刘东等，2015）。所以必须对噪音进行有效抑制，一般使用高灵敏度的PMT和光子计数进行信号探测。拉曼激光雷达方程为如下所示：

式中，P(λR，z)为拉曼通道回波信号，λR为经过氮分子拉曼散射后的波长，λL为激光波长，CR为拉曼通道的雷达常数，Y(z)为雷达接收视场角内的重叠系数，N(z)为高度z处的氮气分子数密度，σN2为氮气分子后向散射截面，αm(λL，z)和αm(λR，z)分别为高度z处的大气分子在λL、λR的消光系数，αa(λL，z)和αa(λR，z)分别为高度z处的大气气溶胶分子在λL、λR的消光系数。

拉曼方法求解气溶胶消光系数的公式为（Ansman等，1992）

式中，N2N(z)是高度为z处的氮气分子数密度，一般氮气在大气中含量比较稳定常用作参考值，可以利用标准大气模式获得。αm(z)、αa(z)分别是高度z处的大气分子和气溶胶消光系数。λL、λR分别为激发波长和拉曼散射波长。气溶胶消光波长指数V一般取1。

此外，可选择石英作为探测矿物气溶胶的示踪剂，利用拉曼激光雷达探测二氧化硅的拉曼散射，估计气溶胶中石英或沙尘的质量浓度（Tatarov 和Sugimoto，2005；Tatarov 等，2008）。由于石英的拉曼后向散射系数几乎是恒定的。在探测矿物气溶胶，特别是在球形和非球形气溶胶的混合情况下，这种方法相较于偏振探测方法具有显著的优势，且有助于表征大气气溶胶和矿物粉尘的研究，对区分沙尘和冰云特别有用。

2.4 高光谱分辨率激光雷达

由于米散射激光雷达反演气溶胶参数的精度较低，而拉曼激光雷达又由于信号较弱，白天探测能力受限。Fiocco 和DeWolf 于1968年提出了高光谱分辨率激光雷达（HSRL）（Dewolf 和DeWolf，1968）。其基本原理是利用大气气溶胶粒子和大气分子光谱宽度的不同，使用窄带滤光器将回波信号的大气气分子瑞利散射和气溶胶米散射分开，从而实现气溶胶的高精度精细探测（刘东等，2015），如图3所示。

图3 大气回波散射信号的光谱图（刘东，2015）Fig.3 Spectrum of atmospheric echo scattering signal（Liu et al.，2015）

分离后的空气分子信号和气溶胶信号光子数为（刘东等，2015）

式中，Kmol和Kaer包含了所有与测量高度无关的系统常数；fam和fmm分别为分子通道中气溶胶信号和分子信号的通过率；faa和fma分别为气溶胶通道中气溶胶信号和分子信号的通过率，这些参数均可通过实验前的校准过程获得。又因为βmol可以根据大气密度直接计算得到，所以通过式（8）和式（9）即可反演得到气溶胶后向散射系数βaer以及气溶胶消光系数αaer的值。

HSRL 是一种将以精细分光技术为核心的激光雷达，具有较高的光谱分辨能力。而目前为了实现激光系统的小型化，促进星载HSRL的发展，将多纵模激光器应用于HRSL 也是一大热点（成中涛等，2017），这一技术不需要对谐振腔实施精密控制，极大地提高了激光器的稳定性和环境适应性。鉴频技术是HSRL实现对大气参数精确探测的关键所在，其核心是滤波器，主要有F-P 干涉仪滤光器、原子/分子吸收滤光器以及视场展宽迈克尔逊干涉滤光器（FWMI）等。1983年Wisconsin大学的第1 台基于F-P 干涉仪的HSRL（Shipley 等，1983）。Shimizu 等（1983）提出了在HSRL 中使用窄带原子吸收滤波器测量大气参数的概念，利用原子或分子的特征吸收峰实现对米散射和瑞利散射的分离。之后钡原子吸收滤光器、碘分子吸收滤光器相继被应用；Hair 等（2001，2008）分别利用碘分子吸收滤光器搭建了HSRL 对大气参数进行测量。目前碘吸收滤波器是应用最多、性能最好的滤波器，但使用要求高和适用范围较小，所以HSRL的发展还是受限于硬件稳定性。

碘吸收滤波器所测到的两个通道的后向散射信号强度为

式中，T(z)=TR(z)Ta(z)；fa、fR分别为激光、气溶胶与碘分子光谱的卷积结果和激光、大气分子散射和碘分子光谱卷积的结果；CM、CF分别为两路的系统常数，其中包括激光能量、光学效率、探测器的量子效率和激光雷达的距离因子等。对式（10）、（11）两式进行相除，就可以得到气溶胶后向散射比R：

2.5 荧光激光雷达

荧光激光雷达常被用于探测生物气溶胶。与拉曼散射相似，荧光效应是激光与大气分子或原子相互作用的能级跃迁造成的，并且荧光效应产生的是连续较宽光谱。紫外激光诱导荧光技术可通过不同种类生物气溶胶的吸收光谱和发射光谱，将有机粒子从无机气溶胶粒子区分出来，再经过比较不同固有荧光谱带的光谱强度，进行反演计算，获得生物气溶胶的相关信息。多数生物气溶胶主要的组成成分为芳香族氨基酸和辅酶，荧光光谱范围多在200—600 nm（饶志敏等，2017）。为了保证回波信号质量要求尽可能的避免杂散光的影响，提高荧光回波强度，一般激光发射波长选择为266 nm或355 nm。

荧光光谱强度受温度、pH 值、激光能量、气溶胶浓度的影响。在实际应用中，荧光雷达回波信号的强度较弱，探测距离有限，一般采用光子计数法探测，通常辅以氮气拉曼信号，有效保证精度的同时避免繁琐校准。

荧光激光雷达方程为（Measures，1984）

Stowers 等（2006）成功利用四倍频266 nm 紫外激光测量了生物气溶胶颗粒的荧光光谱。德国CBRN 研究中心实现了1064 nm 波长探测弹性散射信号，532 nm 波长测量偏振，266/355 nm 作为生物气溶胶荧光激发波长的多波长雷达系统（Joshi 等，2013）。Wojtanowski 等（2015）构建了一台266/355 nm作为激发波长的激光诱导生物气溶胶荧光雷达。目前，在激光诱导荧光探测生物气溶胶的研究方面，国内仍处于起步阶段。2014年兰州大学黄建平教授团队成功搭建了一台多波段拉曼荧光激光雷达系统，355 nm 作为激发光源，接收38 个波段，对大气状态实现了高精细探测（2014）。西安理工大学饶志敏等构建了一台用于大气生物气溶胶预警和探测的双波长激光诱导荧光雷达系统，覆盖紫外到红外波段，能够实现对生物气溶胶云团的预警跟踪和辨别（饶志敏等，2017）。

2.6 其他类型激光雷达

随着光电技术快速发展，近年来也涌现出一些新型激光雷达用于探测气溶胶，如沙式激光雷达、CCD探测激光雷达和LED激光雷达等。

沙式雷达技术以沙式成像原理为基础，可实现大气后向散射信号的距离分辨探测。沙氏成像原理的优点在于激光雷达回波信号不随距离的平方衰减，不需要运用单光子计数技术，降低了系统对测量动态范围的要求；且近场探测盲区小于100 m，有效解决了传统脉冲式大气激光雷达近距离探测盲区大的问题（梅亮，2018）。2015年，瑞典隆德大学联合挪威纳斯克公司最早利用SLidar技术开展了大气气溶胶遥感探测实验研究工作，初步论证了SLidar技术在大气气溶胶探测方面的可行性（Brydegaard，2014；Mei和Brydegaard，2015）。

沙式激光雷达方程可表示为

式中，λ是波长；z是测量距离；P0(λ)是激光器的输出功率；K是系统常数；O(r)是几何重叠因子；β(λ，z)是大气的后向散射系数；α(λ，z)是大气消光系数；P(λ，z)是SLidar 系统探测到的距离分辨的大气后向散射信号。

CCD 雷达，也称侧向激光雷达。该系统将发射装置与接收装置分两处放置，激光器发射出激光，利用大张角的CCD 相机对激光光柱进行拍照，进而分析被探测区域的大气性质（陶宗明等，2014）。侧向激光雷达在近距离段的探测精度很高，可以实现无盲区探测，没有几何因子的问题，特别适合于近距离段的大气气溶胶探测，缺点在于远距离段空间分辨率不高、信噪比较低。目前常把CCD 探测与后向散射激光雷达联合在一起进行工作。Barnes 等（2003）利用CCD 相机侧向激光雷达系统在夏威夷的莫纳罗亚山上对大气边界层的气溶胶空间分布进行了探测。国内相关单位利用CCD 雷达对近地面的气溶胶进行了探测研究（史博等，2015；麻晓敏等，2018；陶宗明等，2019）。

LED光源雷达是近年来出现的一种新型主动式遥感雷达。LED雷达的优点在于选用LED作为探测光源，波长丰富，在进行多波长探测时容易实现气溶胶高分辨率的粒子谱分布探测（Thurmond 等，2016）。但是目前来说LED 雷达光源的发散角太大，发射能量偏低，在实际应用中需要复杂的发射光源整形系统。日本千叶大学Tatsuo等学者研制了第一台LED 光源大气遥感雷达，该雷达可以实现夜间200 m范围内的气溶胶探测。由于LED雷达体积小、坚固耐用，Tatsuo 和Sonoko 等在2016年设计了一款可用于火星车探测烟雾和灰尘的LED雷达（Shiina，2013；Shiina 等，2016）。钟文婷等（2016，2018）设计并研发了新型单波长LED 光源雷达系统，探测西安城区上空气溶胶分布。

3 探测气溶胶性质的进展

3.1 光学特性

气溶胶散射信号与粒径谱、波长、大气折射率、浓度等都有关系。通过测量不同波段的回波信号，可以获得气溶胶光学特性，如消光系数、后向散射系数、雷达比（LR）、气溶胶光学厚度（AOD）、Ångström 波长指数等。气溶胶消光系数和后向散射系数是激光雷达方程的两个未知参量。后向散射系数越大表明大气中颗粒物的散射能力越强，反之越弱。影响后向散射系数的因素较多，例如颗粒物浓度，粒子半径，相互之间的关系较复杂。气溶胶消光系数反映了气溶胶消光能力。Ångström 波长指数是度量米散射粒子尺度分布的重要参数。Ångström波长指数的值一般在0—4.0的范围内，与气溶胶的大小密切相关，揭示了消光对入射光波长的依赖关系。当气溶胶粒子半径与入射激光波长相比很大时，Ångström 指数接近于0，并随着粒子尺寸的减小而变大（狄慧鸽等，2014）。

目前通过多通道激光雷达系统，可实现气溶胶消光系数、后向散射系数和激光雷达比等光学特性垂直分布的。太阳光度计提供的准确AOD 值作为约束条件为精确求解米散射激光雷达方程给予了极大帮助，拉曼或者高光谱分辨率激光雷达的使用使反演总消光系数误差降低到10%，甚至后向散射系数的误差降到5%（Bösenberg 和Hoff，2017）。但是要得到较为准确的气溶胶光学特性，难点在于反演算法的不稳定。国内许多学者一直致力于反演方法与激光多次散射问题的研究，发展了不同的多次散射激光雷达方程（吕达仁等，1977；周诗健等，1981；邱金桓等，1992，2003）。

激光雷达在探测气溶胶光学特性方面具有明显的技术优势，例如，黄忠伟等人通过多部激光雷达追踪2010年春季发生在东亚地区的一次大尺度沙尘暴，发现亚洲沙尘气溶胶在5 d 内从中国西北部迅速传输到北极。其中SACOL、仙台两地的激光雷达探测沙尘层垂直分布结果如图4所示（Huang等，2015b）。

图4 位于SACOL、仙台的2部激光雷达追踪同一沙尘事件的观测结果（Huang等，2015）Fig.4 Observation results of two lidar systems in SACOL and Sendai tracking the same dust event（Huang et al.，2015）

3.2 大小

气溶胶尺度大小常用色比CR（Color Ratio）来表示，往往是长波变量（如衰减后向散射系数或消光系数）除以短波变量的比值，即较大的色比值表示颗粒较粗，反之则反。下面给出以532 nm 波长和355 nm波长之间的色比：

实际应用中，很难通过色比大小来定量确定粒子尺寸。一般云的色比往往大于气溶胶的，而沙尘气溶胶的色比大于城市污染气溶胶的色比。

气溶胶粒子大小也会随环境因素发生改变。例如，相对湿度的增加会使气溶胶粒径变大，且相对湿度越大，粒径变化越明显。2008年中美联合沙尘天气试验在中国西北进行，利用3个多脉冲激光雷达（MPL）系统以及其他地面仪器和星载遥感技术，测量了沙尘暴的沙尘气溶胶垂直廓线及其远距离输送。图5 显示为CALIPSO 观测到高浓度气溶胶层，该气溶胶层出现在北京时间03：50左右，其色比介于0.6—0.9（Huang等，2010）。

图5 CALIPSO星载激光雷达探测2008-05-02中国西北地区沙尘暴期间的沙尘色比大小（Huang等，2010）Fig.5 The size of dust color ratio detected by CALIPSO during the sandstorm in Northwest China on May 2，2008（Huang et al.，2010）

3.3 形状

气溶胶退偏比可反映其形状的不规则性，理论上球形粒子的退偏比为零，非球形粒子的退偏比大于零。退偏比被认为是鉴别沙尘气溶胶的重要参数，532 nm波长的气溶胶退偏振比达到30%—35%时，可以表明是纯沙尘（Sugimoto 和Lee，2006；Freudenthaler 等，2009），冰晶粒子也有较高的退偏比，例如卷云（Mishchenko和Sassen，1998）。

Scholand 等（1971）首先利用偏振激光雷达对云进行了探测。经过几十年的发展，目前偏振激光雷达已被广泛应用于探测气溶胶和云的垂直分布（Ansmann 等，1992；Cairo 等，1999；McGill 等，2015）。偏振激光雷达在区分云相位具有独特的探测能力（Sassen，1991）。Liu 等（2002）对亚洲沙尘暴沙尘气溶胶的退偏比和激光雷达比进行了统计分析，更准确地区分云和气溶胶粒子层。Sassen等（2001）根据激光雷达区分了各种大气目标，如冰水云、幡状云、降水和气溶胶等。Liu等（2008）利用CALIPSO激光雷达观测的层平均退偏比研究了对流层沙尘粒子的特征。周天等（2013）利用层积分后向散射衰减系数和层积分退偏比之间的关系识别塔克拉玛干沙漠上空的沙尘。Burton 等（2014）提出了一种参数化方案，整合了激光雷达比、后向散射比和退偏比三者用以识别气溶胶类型。

应指出的是，研究还发现退偏比不仅取决于颗粒形状，与颗粒的大小和浓度也有关系。青藏高原夏季的沙尘气溶胶在不同海拔高度上的体退偏比发现具有显著的差异（Huang 等，2007）；大气中相对湿度的增加会使非球形气溶胶转化为球形气溶胶，且非球形气溶胶与球形气溶胶随着海拔升高而逐渐增加（孙冉等，2017）；沙尘气溶胶随着远距离传输和抬升高度增加，受自身重力沉降作用以及其他污染物的吸附和包裹影响，造成不规则程度降低（茹建波，2017）；Huang 等（2020）发现沙尘气溶胶在532 nm 和355 nm 处的DRs比值与气溶胶的吸收系数有很好的关系，证实了沙尘气溶胶的吸收特性导致了紫外和可见波长的退偏比差异。如图6所示，气溶胶和云在不同波段退偏比具有显著的差异。

图6 兰州大学自主研发的双波段偏振激光雷达于2014-03-11在兰州市探测的云和气溶胶垂直结构（齐思齐，2020）Fig.6 Vertical structure of cloud and aerosol detected by dual-band polarization lidar independently developed by Lanzhou University in Lanzhou city on March 11，2014（Qi et al.，2020）

3.4 组分

大气气溶胶的种类很多，按其来源可以分为自然源和人为源，自然源主要是源于海洋、土壤、火山爆发、森林燃烧的粒子，人为源主要是源于人为活动，如化石燃烧等。由于不同气溶胶类型具有不同光学物理特性，为区分气溶胶提供了物理基础。并且气溶胶排放源的空间、时间非均匀性、以及气溶胶老化和二次生成引起的气溶胶变化，也让通过激光雷达来开展气溶胶类型探测研究显得非常必要。目前随着地基激光雷达观测站的增加、星载激光雷达的全球覆盖和观测时间的累积，结合地基和卫星资料，发挥各自所长对气溶胶进行分类，使划分的气溶胶类型更为准确合理，有望极大地提高对气溶胶在全球分布及其气候效应的认识。利用不同波长的偏振测量可以显著改善气溶胶分类（Burton 等，2012，2014；Huang等，2020）。目前，CALIOP 激光雷达观测数据已被广泛用于评估多种气溶胶类型的气溶胶模型模拟。CALIPSO 卫星气溶胶和云的官方分类识别算法经历了多个版本的改进，2016年11月发布的版本4.10（V4）2级气溶胶数据产品，对气溶胶分类和激光雷达比值选择算法进行了重大改进。重新定义气溶胶分类规则，更新所有类型气溶胶的激光雷达比率，使气溶胶的分类更加完善合理，可靠性大大提高。图7 显示了CALIPSO 卫星经过东欧到北非的气溶胶类型分布，通过气溶胶分类算法发现6 种气溶胶类型都出现在轨道上（Omar等，2008）。

图7 CALIPSO星载激光雷达在2006-08-12观测的云和气溶胶分类结果（Omar等，2008）Fig.7 Cloud and aerosol classification results observed by CALIPSO on August 12，2006（Omar et al.，2008）

基于CALIPSO 星载激光雷达观测，Huang 等（2015b）融合了MODIS 地表类型、土地利用等多元数据，根据自然沙尘与人为沙尘的后向散射以及退偏比等光学特性的差异，提出了一种人为沙尘气溶胶识别的新方法。

3.5 浓度

一般来说，激光雷达回波信号强度越大，说明气溶胶浓度较高，光学厚度也较大。于是气溶胶光学厚度也被广泛用于确定大气中颗粒物的浓度。

现阶段针对全球尺度的气溶胶质量浓度分布的研究主要还是利用被动遥感，从AOD 数据中提供区域或全球范围内白天颗粒物浓度的分布。但是被动遥感缺少夜间信息与高垂直分辨率信息，且在高反照率地区反演不确定性较高、受云的影响严重等，难以满足高精度气溶胶反演研究的需求（Stubenrauch等，2013）。

因此，通过高时空分辨率的主动遥感长期连续观测获得可靠的全球颗粒物浓度非常重要。国内外学者提出了利用激光雷达观测反演颗粒物质量浓度垂直分布的方法（Koelemeijer 等，2006；Lin 等，2015；Tao 等，2016；Toth 等，2019）。近些年，深度学习技术也被引入反演气溶胶质量浓度的研究中（Chen等，2018；Shen等，2018），但这些研究没有充分考虑气溶胶类型对反演结果的影响，难以满足支持全球气候变化研究的需求。华雯丽等人结合气溶胶退偏比进行粗细粒子分离的反演方法，在实时反演沙尘质量消光系数的基础上，反演了敦煌地区沙尘的质量浓度垂直分布特征（华雯丽等，2018）。Ma 等（2020）根据CALIPSO 星载激光雷达观测，开发了一种反演全球PM10&PM2.5质量浓度分布的方法，为全球范围内的颗粒物质量浓度提供了一个可靠的高分辨率数据库。2007年—2016年全球平均PM10和PM2.5质量浓度分布昼夜差异结果。

目前通过传统的算法反演得到气溶胶消光系数廓线和基于粒子谱算法反演得到颗粒物粒子谱分布，结合由米散射理论可以求得颗粒物质量消光效率（MEE）（饶志敏，2018）：

式中，r气溶胶粒子半径；ρ为粒子标准密度，ρ＝2（g · cm-3）；Qext为球形粒子消光效率；n为粒子谱分布情况；m为负折射率；λ为激光波长。

最后，通过气溶胶消光系数与质量消光效率相结合获得颗粒物的质量浓度：

但是由于现有卫星遥感技术的局限性，直接获取MEE 很难。目前常用的方法是地基偏振激光雷达和太阳光度计相结合获得MEE。首先分离沙尘和污染物，获得了沙尘MEE，王天河等（2021）证实了基于微脉冲激光雷达MPL 观测可有效获得沙尘质量浓度（DMC）剖面，并评估上述不确定性的影响。2012-04-24 敦煌站激光雷达观测到沙尘过程的结果，以及沙尘质量浓度反演结果，如图8所示（Wang等，2021）。

图8 2012年4月24日敦煌站激光雷达观测沙尘过程（Wang等，2021）Fig.8 Dust process observed by liDAR at Dunhuang Station on April 24，2012（Wang et al.，2021）

4 结语

激光雷达是大气遥感与环境监测不可或缺的手段，是跨领域多学科交叉的应用系统，需要各技术领域互相协调。目前国内外已经进行了大量有关气溶胶特性的激光雷达观测与研究，积累了大量宝贵的观测资料，取得一批批有价值的成果。随着光电技术、人工智能、精密加工等方面的进一步发展，激光雷达发展必将再上一个新的台阶。今后，气溶胶激光雷达的主要发展趋势有以下几方面：

（1）硬件技术方面：将朝着小型化、智能化、易维护、高稳定性等方面发展，甚至在不久的将来实现可搭载在无人机平台上的气溶胶激光雷达；

（2）反演算法方面：随着探测参量越来越多，且人工智能算法的加入，将获得更为准确、更丰富的气溶胶参数，例如粒子谱、复折射指数、吸收系数、单次散射反照率SSA等。

（3）地基雷达观测网：目前，全球已经有多个地基激光雷达网，今后将会有更多覆盖更大范围的气溶胶观测网络陆续建立或完善，如兰州大学正在建设的“一带一路”激光雷达网。这些观测网有望发展成为各具特色的、高度融合的全球性数据共享激光雷达网络；

（4）星载激光雷达观测：CALIPSO 卫星已提供连续15 a 的全球气溶胶和云三维时空分布观测资料，2022年4月中国成功发射了大气环境监测卫星，这是全球首颗激光雷达二氧化碳探测卫星，也是中国第一颗可探测全球气溶胶和云三维分布的星载激光雷达卫星。今后还会发射更多的可探测气溶胶的星载激光雷达，例如欧空局和日本宇航局合作的EarthCARE 卫星，以及中国即将发射的另一颗星载激光雷达。这些丰富的星载激光雷达观测数据，将极大地推动气溶胶环境与气候效应相关研究。