激光雷达探测整层大气昼夜气溶胶光学厚度

汪惜今 , 徐青山 , 范传宇 ,4, 程晨 , 戚鹏 ,5, 徐赤东

( 1 中国科学院合肥物质科学研究院安徽光学精密机械研究所基础科学研究中心, 安徽 合肥 230031;2 中国科学技术大学, 安徽 合肥 230026;3 中国科学院合肥物质科学研究院安徽光学精密机械研究所, 中国科学院通用光学定标与表征技术重点实验室,安徽 合肥 230031;4 皖西学院电气与光电工程学院, 安徽 六安 237012;5 安徽建筑大学电子信息工程学院, 安徽 合肥 230601)

0 引 言

大气气溶胶是液态或固态微粒在空气中的悬浮体系, 是地球大气的主要成分之一[1]。气溶胶通过对太阳辐射能的吸收和散射影响气候变化[2]。气溶胶对太阳辐射的影响通常采用气溶胶光学参数来表征。针对大气气溶胶光学参数的探测, 一直是大气科学领域学者的密切关注对象。气溶胶光学厚度 (AOD) 表征了气溶胶对太阳辐射的衰减特性, 是气溶胶光学参数中的关键参数。

AOD定义为沿着传播路径, 由大气柱中的气溶胶引起的消光量。探测AOD的技术手段主要分为被动探测技术和主动光学探测技术。被动探测技术包括太阳光度计探测、月光和星光探测, 太阳光度计应用广泛、稳定性强, 但是无法在夜间进行观测; 月光和星光探测能够在夜晚探测AOD, 是夜间气溶胶被动遥感手段的有效补充, 但其工程应用仍在探索之中。激光雷达作为一种主动探测工具, 具备昼夜观测能力, 为夜间观测气溶胶光学参数提供了一种有效手段。王珍珠等[3]使用米散射激光雷达在北京地区进行了探测, 指出大气边界层内的气溶胶粒子浓度较高, 大气边界层高度在夏季相对稳定。周碧等[4]利用连续5年的激光雷达资料,反演得到兰州地区的气溶胶消光系数垂直廓线, 指出消光系数的反演会受到沙尘等特殊天气的影响。

在气溶胶消光系数廓线遵循指数型分布的情形下, Penndorf[5]指出消光系数减少至地面1/e时的高度称为气溶胶标高, 结合激光雷达探测数据, 可以获得AOD。Qiu等[6]分析了1994—2001年中国11个站点的季节性气溶胶标高和消光系数廓线的特征, 得出气溶胶标高变化与地域以及季节密切相关, 多日数据平均值反演得出的消光系数更符合指数型分布。范伟等[7]根据大气AOD和近地面水平能见度的测量结果, 得出了气溶胶随时间、季节以及地理位置变化的规律: 秋季内陆地区标高最大, 冬季次之, 夏季沿海地区气溶胶标高最小。李成才等[8]应用MODIS遥感产品对北京及周边地区的气溶胶进行了分析, 发现受沙尘暴天气的影响, 春季气溶胶标高较高, 夏季由于混合层高度较高, 气溶胶标高明显高于冬季; 此外, 将AOD和能见度拟合得到的AOD产品, 在同一季节内进行对比, 不同季节内进行校准, 可用于研究地面污染等有关问题。韩永等[9]利用能见度仪和太阳辐射计获得了气溶胶标高水平方向上的变化特征, 指出气溶胶标高与波长成反比, 夏季气溶胶标高小于冬季, 但是由于观测条件限制无法得出夜间的气溶胶标高。

为了验证激光雷达用于整层大气AOD测量的可行性, 胡顺星等[10]将雷达所得结果与太阳辐射计的测量数据进行了分析对比, 发现两者具有较好的一致性, 证明将激光雷达用于测量整层大气AOD 的方法可行。陈莎莎等[11]利用Fernald后向积分法反演得出气溶胶垂直消光系数廓线, 并将消光系数廓线分为四个不同的类型, 从而拟合得出气溶胶标高, 计算出整层大气AOD, 与太阳高度计所测整层大气AOD进行对比, 平均相对误差小于6.7%, 一致性较好。

本文利用激光雷达进行昼夜连续观测气溶胶光学参数的特性, 得出气溶胶垂直方向上的消光系数廓线,并在此基础上得出气溶胶标高; 进而基于地面能见度数据推算出近地面消光系数, 将二者相结合, 计算得出整层大气AOD; 最后以合肥地区一个月的观测数据为例, 初步尝试分析了整层大气AOD的昼夜变化特征。

1 仪器和基本原理

1.1 仪器

实验采用波长为 532 nm 的微脉冲激光雷达 (MPL) 进行探测, 用于研究合肥地区大气气溶胶的光学特性。MPL精度高、体积小、移动性好, 对颗粒物、大气能见度等可开展连续观测工作, 获得其时空分布特征和变化规律[12]。

激光进入大气后, 会受到传输路径上气溶胶粒子以及空气分子的衰减, 发生弹性散射, 最终被探测器接收。激光雷达方程可表示为[13]

1.2 基本原理

已有的研究表明, 气溶胶浓度随高度呈指数型下降, 从而气溶胶消光系数σz随高度变化也遵循指数型下降[14], 其表达式为

式中σ0为近地面消光系数,H为气溶胶标高。

光学厚度τ是消光系数在垂直方向上的积分, 可以表示为

消光系数的反演一般有斜率法[15]和 Fernald法[16]。通常情况下, 斜率法常用于水平探测, 前提是大气均匀分布, 而Fernald法常用于垂直探测。吕立慧等[17]对这两种激光雷达反演算法进行了对比, 指出Fernald法更适用于非理想状态下消光系数的反演, 相对误差较小。

近地面消光系数可以利用激光雷达进行水平测量获得, 但由于不可能频繁交替测量垂直和水平数据, 因此可以使用地面能见度数据推算近地面消光系数。已有研究结果表明, 利用能见度计算得出的近地面消光系数和激光雷达水平方向测得的消光系数具有较好的一致性[18,19]。根据Koschmieder 定律, 能见度V可表示为

式中ε定义为人眼视角阈值, 取值为 0.02[20], 由此可得

根据安徽省气象观测站 (58321) 提供的合肥城区的近地面能见度资料, 利用式 (6) 可计算得出近地面消光系数。

基于垂直方向上的气溶胶消光系数廓线, 可以求解得出气溶胶标高。利用Fernald后向积分法可计算得出气溶胶消光系数垂直廓线, 其计算公式为

式中σa(r) 和σm(r) 分别代表气溶胶粒子、空气分子在高度r处的消光系数;Sa是气溶胶消光后向散射比, 本实验的观测点为合肥城郊, 该地区的气溶胶消光后向散射比Sa通常取值为 50 sr 较为合理[21];Sm是空气分子消光后向散射比, 通常取 8π/3 。rc为标定高度, 一般选择气溶胶含量微乎其微的大气层所在高度, 并且可以认为这一高度上 [P(r)r2]/βm(r) 有最小值, 本研究选取 4～6 km 为标定高度;σa(rc) 和σm(rc) 分别为标定高度处气溶胶粒子和空气分子的消光系数的初值, 标定高度处气溶胶消光系数的初值由该处的气溶胶散射比R决定 [R= 1+βa(rc)/βm(rc) = 1.01], 空气分子的消光系数则取决于大气模式[22]。

2 数据处理及结果分析

利用MPL 2020年1月的观测数据开展数据处理和分析验证工作。观测数据的时间分辨率为 1 min, 能见度的观测数据同步为 1 min 数据。

2.1 气溶胶消光系数廓线的计算

利用MPL接收到的回波信号, 通过Fernald法求解出气溶胶消光系数垂直廓线。图 1 (a)、(b) 分别给出了北京时间2020年1月2日夜间01:42和白天08:30垂直方向上的气溶胶消光系数廓线, 其中所选两个时刻的有效数据均满足晴天无云条件。对两组观测数据分别进行指数型拟合, 发现观测期间气溶胶消光系数垂直分布基本符合指数型分布特征。对比实测的昼夜两组气溶胶消光系数廓线, 可以看出, 白天气溶胶消光系数廓线层结构较为丰富, 没有夜间平滑, 这可能是由于白天对流层的垂直活动较为剧烈, 这种气溶胶消光系数廓线的昼夜差异将会影响昼夜气溶胶标高的变化。

图 1 气溶胶消光系数廓线。(a) 夜间; (b)白天Fig. 1 Aerosol extinction coefficient profiles. (a) Night; (b) day

2.2 气溶胶标高的计算

气溶胶标高的计算采用陈莎莎等[11]提出的方法, 将气溶胶消光系数廓线分为四种类型。

1) 如果气溶胶消光系数廓线在垂直方向上大致服从指数型变化, 那么可以根据式 (3) 直接采取指数拟合的方法得出标高。

2) 如果气溶胶消光系数廓线在边界层上方呈现指数型变化, 边界层内气溶胶充分混合, 则有

标高H=H′a+H1。

3) 如果边界层内气溶胶分布不均匀, 在某一高度层内有积累, 则有

针对上述四种方法拟合计算出的气溶胶标高与太阳光度计的结果进行对比, 发现两者具有较好的相关性, 证明上述方法可行。

利用Fernald后向积分法得到2020年1月2日气溶胶消光系数廓线, 通过对消光系数进行处理得到气溶胶标高。图 2为利用MPL测得的全天数据进行处理得到的当日气溶胶标高, 数据按小时平均。从图中可以看出, 气溶胶标高值呈现早晚低、中午高的趋势, 凌晨气溶胶标高值较小, 随着时间的推移, 气溶胶标高值白天逐渐升高, 夜晚又逐渐减小, 昼夜气溶胶标高值存在差异。

图 2 气溶胶标高的日变化Fig. 2 Diurnal variation of aerosol scale height

2.3 近地面消光系数的计算

根据合肥地区2020 年1 月份的能见度数据, 利用式 (6) 计算得出1 月2 日的近地面消光系数, 如图 3 所示。可以看出, 近地面消光系数从00:00开始由0.5缓慢减小, 然后不断升高, 在12:00达到当日最高值0.637,之后又逐渐减小, 到了夜晚又呈现增大趋势。将能见度得到的消光系数和激光雷达MPL得到的结果进行对比, 结果如图 4所示。由图可知, 相关系数为 0.8821, 证明利用能见度获得消光系数是可靠的。

图 3 近地面消光系数分布Fig. 3 Distribution of surface extinction coefficient

图 4 近地面消光系数结果 (能见度仪和MPL) 对比Fig. 4 Comparison of surface extinction coefficients from MPL and visibility

2.4 AOD日变化的计算

通过前面计算所得气溶胶近地面消光系数和标高, 根据式 (4) 可以得出1月2日合肥地区整层大气AOD分布, 结果如图 5 所示。由图可知, 1 月2 日的AOD 呈现单峰型, 凌晨AOD 值较小, 在05:00 为当天最小值0.191; 随着太阳的升起, AOD值逐渐增大, 在14:00达到当日最大值1.38, 随后又不断降低。当日AOD昼夜分布也存在差异, 白天AOD平均值为 0.783, 大于夜晚AOD平均值 0.315。

图 5 AOD日变化Fig. 5 Diurnal variation of AOD

2.5 气溶胶光学参数的月变化及影响因素分析

图 6给出了2020年1月连续一个月的气溶胶标高和AOD的昼夜分布计算结果。从中可以看出, 1月白天和夜晚的气溶胶标高和AOD均存在差异, 都通过了显著水平为0.05的检验, 白天气溶胶标高和AOD明显大于夜晚。

图 6 2020年1月气溶胶标高 (a) 和AOD (b) 昼夜对比Fig. 6 Day-night comparion of aerosol scale height (a) and AOD (b) in January 2020

图 7 气溶胶标高和相对湿度的关系Fig. 7 Relationship between aerosol scale height and relative humidity

进一步对于影响气溶胶标高和AOD的气象因素进行了分析。图 7表示了2020年1月的气溶胶标高和相对湿度之间的相关性, 其中相对湿度数据为逐日平均值, 来源于安徽省气象观测站 (58321), 由图可知, 相对湿度和气溶胶标高呈现显著的负相关关系, 通过了99%的显著性水平检验, 即气溶胶标高随着相对湿度的减少而增大。造成这种现象的原因, 可能是由于相对湿度的增加, 气溶胶粒子对水汽的吸收能力增强, 使气溶胶粒子本身的质量不断增大, 受重力的影响出现沉降作用, 导致气溶胶标高减小[23]。1月气溶胶标高和相对湿度的相关系数为 -0.881, 证明两者之间存在明显的负相关。

此外, 还对2020年1月AOD和能见度的关系进行分析, 结果如图 8所示。总体来看, 1月昼间、夜间和总的AOD 与能见度均呈现负相关关系。值得注意的是, 1 月白天AOD 与能见度的相关性要强于夜晚, 白天AOD与能见度的相关系数为 -0.701, 而夜间的相关数据为 -0.612。

3 结 论

利用MPL接收到的回波信号, 通过Fernald法求解出气溶胶消光系数垂直廓线, 计算得出气溶胶标高; 同时利用能见度和消光系数的关系得到近地面水平方向的消光系数, 由此计算出AOD。对合肥地区的2020年1月连续一个月的昼夜气溶胶标高以及AOD进行了计算分析, 结果显示, 两者昼夜均存在差异。进一步对影响的气象因素进行分析, 得出气溶胶标高和相对湿度呈现负相关关系, 相对湿度大则气溶胶标高值小, AOD和能见度也存在负相关关系。提出的方法可为连续获得整层大气AOD 昼夜分布特征、以及气溶胶标高和AOD的季节性差异提供技术参考。