基于多轴差分吸收光谱技术的合肥地区气溶胶观测研究

李晓梅，谢品华，3*，徐 晋，李 昂，田 鑫，任 博，胡肇焜，吴子扬

1. 中国科学院安徽光学精密机械研究所，中国科学院环境光学与技术重点实验室，安徽 合肥 230031 2. 中国科学技术大学科学岛分院，安徽 合肥 230026 3. 中国科学院区域大气环境研究卓越创新中心，福建 厦门 361021

引 言

大气气溶胶是大气成分的重要组成部分，是大气液态或固态悬浮微粒的总称。气溶胶在大气物理中的作用是影响地球辐射平衡，作为云凝结核，影响云的特性；另外，它在大气化学转化中也起到了重要的作用，参与多种大气化学过程，为气体转化提供场所，加速了污染气体的生成[1]。近年来在我国京津冀地区、华北平原、长江中下游地区等工业发达、人口稠密的特大城市和城市群周围频繁发生严重的雾霾事件，严重危害人类的健康。有研究发现[2]，在中国不同地区的大城市中，大多数气溶胶颗粒是次生的，即在霾污染事件期间通过前体气体的光化学反应形成。因此，研究气溶胶的光学特性非常必要。

气溶胶的光学遥感探测技术主要有太阳光度计、激光雷达等，太阳光度计可以提供较为准确的气溶胶光学厚度(AOD)，目前已经建立全球气溶胶地基观测网络(AERONET)，但是无法提供气溶胶垂直消光系数廓线；激光雷达技术是一种有效的大气探测技术，能够实现对大气颗粒物的高精度、高时空分辨率的连续观测，但是在气溶胶较为活跃的边界层存在盲区[3], 且价格昂贵。

多轴差分吸收光谱仪(MAX-DOAS)通过多个仰角观测太阳散射光，再利用差分吸收光谱技术反演出对流层的痕量气体和气溶胶的垂直分布[4]，在国际上很多研究机构已经开展了基于该技术的气溶胶观测研究。Zieger等2009年在荷兰Cabauw利用MAX-DOAS研究了气溶胶消光系数对湿度(RH)的依赖性[5]；Irie等在日本筑波市进行了MAX-DOAS，腔衰荡光谱仪，激光雷达和天空辐射计的协同气溶胶观测，评估MAX-DOAS在467 nm反演的气溶胶消光系数(AESC)和AOD[6]；Davis等2013年8月至9月在麦克默里堡北部的艾伯塔省开展MAX-DOAS观测，并将气溶胶结果与激光雷达和GEM-MACH空气质量模型作对比[7]；Zhang等2017年4月至12月份在上海郊区开展MAX-DOAX观测，研究O4校正因子的优化调整，以及冬季气溶胶垂直分布特征[8]。

利用地基MAX-DOAS在2017年12月至2018年1月在合肥地区开展两个月的连续观测，反演了气溶胶光学厚度信息，并将观测的AOD结果与安装在同一观测点的太阳光度计CE318的结果进行对比，来验证MAX-DOAS仪器性能以及研究合肥地区的气溶胶光学特性。将MAX-DOAS反演出的AOD和近地面气溶胶消光系数分别与太阳光度计CE318和地面站点的PM2.5作对比；选取了一次雾霾过程，分析了观测期间气溶胶AOD和廓线的变化特征。

1 原 理

1.1 O4DSCD计算原理

差分光学吸收光谱(DOAS)技术，基于朗伯比尔定律[9]，如式(1)

(1)

其中I(λ)是接收的散射太阳光，即经过各种气体吸收、瑞利散射、米散射等衰减后的光强，εR为瑞利散射系数，εM为米散射系数，l为光程长度，σi(λ，P，T)为吸收截面，c为气体浓度，i为气体种类。

散射的影响可以通过多项式拟合等方法去除，进而得到痕量气体吸收的部分，再与标准气体截面拟合，最后得到痕量气体的斜柱浓度。

(2)

MAX-DOAS采用将差分光学吸收光谱(DOAS)技术，通过多仰角的测量获取多角度的O4DSCD。

O4DSCDα=O4SCDα-O4SCD90°

(3)

其中α为不同仰角。

1.2 气溶胶反演原理

由于氧二聚体O4在大气中的浓度基本稳定不变，影响O4差分斜柱总量(O4DSCD)的主要因素就是光的传输路径。在晴朗无云的条件下，光的传输路径主要受气溶胶廓线的影响，因此可以利用O4DSCD反演出大气中气溶胶的信息。研究中采用中科院安徽光学精密机械研究所(AIOFM)与德国马普化学所(MPIC)合作开发的气溶胶消光和痕量气体浓度廓线反演算法(PriAM)[10]结合多角度的O4DSCD进一步反演出气溶胶廓线信息。

首先，测量值是y(α，φ，λ)，即在测量序列中具有相同的扫描方位角φ和O4吸收波长带的不同仰角处的O4DSCD观测值，假设测量值y可以通过正演模型重现，而正演模型结果取决于气溶胶廓线x和气溶胶光学特性。假设每层中的气溶胶均匀分布，使得观测向量(y)可以描述为

y+ε=f(x)+δ

(4)

其中ε和δ分别是测量误差和模型误差，f(x)为大气辐射传输模型(RTM)。

其次，观测向量y中包含的信息很可能不足以反演独特的气溶胶消光特征，因此使用最优估计方法进行气溶胶反演[11]。

大气中的辐射传输是非线性的，因此气溶胶消光的反演通过使用Gauss-Newton方法迭代求解[12]，定义为

(5)

为加速迭代过程，提高反演精度，结合Levenberg-Marquardt算法对Gauss-Newton方法进行修正后的迭代过程变为

(6)

根据式(6)，初始气溶胶廓线xa通过不断迭代得到最接近真实情况的气溶胶廓线。

2 实验部分

2.1 二维MAX-DOAS

二维MAX-DOAS仪器安装在合肥(117.16°E，31.91°N)西北角的科学岛安光所综合实验楼楼顶，位置如图1所示。仪器结构图如图2(a,b)所示，望远镜收集太阳散射光，并通过石英光纤束引入光谱仪单元, 光谱仪光谱分辨率为0.47 nm, 然后光谱仪将范围为293～414 nm的光谱映射到具有2 048像元的面阵CCD探测器，最后传输到计算机中进行数据处理。光谱仪单元需要冷却到+19 ℃的稳定温度，以便使光谱仪的光学性质变化最小并减少探测器暗电流[13]。室外的二维转动平台可以实现0～360°方位角和0～90°仰角的扫描，并且装有电子倾角仪，可现实反馈控制，从而提高仰角指向精度。仪器的正北方设定为方位角0°，垂直扫描选取1°, 2°, 3°, 4°, 5°, 6°, 8°, 15°, 30°, 90°(天顶角)10个仰角。

图1 测量点的位置

图2 二维MAX-DOAS装置结构图和实物图

2.2 太阳光度计

在二维MAX-DOAS仪器北7 m处，运行着一台CE318标准型太阳光度计，如图3所示，用于测AOD精度较高的仪器，从紫外到近红外波段设有8 个通道，其中通道中心波长为340和380 nm所测得的AOD用于与MAX-DOAS的结果做对比。我们使用云屏蔽和去除了突变值的质量有保证的level1.5数据。

图3 太阳光度计CE318实物图

2.3 数据处理

利用DOAS算法，将采集到的光谱在337～370 nm波段反演O4DSCD，取当天中午的天顶光谱作为Fraunhofer参考光谱(FRS)。利用QDOAS软件，将FRS、若干痕量气体吸收截面、Ring谱的三阶多项式和二阶偏移多项式的对数拟合到测量光谱的对数，并针对暗电流和偏移进行了修正，QDOAS参数设置如表1所示。图4(a,b,c,d)分别给出了2018年6月2日13:08在方位角为0°，仰角30°时记录的DOAS拟合的实例。对于每个测量周期，将相应的天顶光谱(仰角90°)作为离轴仰角处的光谱的FRS。在很大程度上消除了对DSCD的平流层贡献。但是，O4DSCD仅受平流层吸收的轻微影响因为O4主要存在于对流层。工作中使用高分辨率吸收截面，通过仪器狭缝函数进行卷积以匹配仪器的分辨率。

表1 QDOAS参数设置Table 1 QDOAS parameter settings

3 结果与讨论

3.1 MAX-DOAS与CE318测量AOD相关性对比

利用QDOAS软件，将测量光谱计算出O4DSCD，然后再利用PriAM算法反演出AOD。由MAX-DOAS和太阳光度计测量的气溶胶光学厚度(AOD)的时间序列如图5(a)所示，AOD的日均值时间序列如图5(b)所示，MAX-DOAS和太阳光度计数据随时间显示出相似的变化趋势。另外，在测量时间段，MAX-DOAS测得的AOD均值为0.71，日均值在0.26～1.90之间。由于MAX-DOAS和太阳光度计的采样时间不同，因此进行小时均值和日均值处理，MAX-DOAS和太阳光度计AOD测量之间的相关性如图5(c)和(d)所示，MAX-DOAS测量结果与AERONET观测值的小时均值Pearson相关系数rhourly=0.92，日均值相关系数rdaily=0.91，均表现出良好的一致性，截距为0.14，且太阳光度计测得的AOD比MAX-DOAS所测结果高出20%。由于MAX-DOAS

图4 QDOAS反演示例(a)：观测光谱与参考光谱；(b)：残差；(c)：NO2DSCD=3.22×1016 (molec·cm-2)；(d): O4DSCD=3.48×1043(molec2·cm-5)

图5 MAX-DOAS和CE318测量的AOD的时间序列图(a)AOD小时均值, (b)AOD日均值MAX-DOAS和太阳光度计，AOD测量之间的相关性(c) AOD小时均值, (d)AOD日均值

Fig.5TimeseriesdiagramofAODdailyandhourlydatameasuredbyMAX-DOASandCE318areshownin(a)and(b),respectively.Linearregressionplotsofdailyandhourlydataareshownin(c), (d)respectively

测量对对流层低层的气溶胶最敏感，而太阳光度计则通过直接太阳测量得到AOD，对整个大气中的气溶胶很敏感[16]。因此，当气溶胶在对流层上层运动时，两次测量之间的差异可能很大[17]。气溶胶光学特性的假设也导致反演到的MAX-DOAS AOD的不确定性。

3.2 近地面气溶胶消光系数

由于在观测地点地表和垂直方向缺乏气溶胶消光系数测量，另一种半定量验证MAX-DOAS反演气溶胶消光系数的方法是将其与颗粒质量浓度进行比较。将地面站点的地表PM2.5浓度与MAX-DOAS反演近地面气溶胶消光系数进行比较，采用气溶胶消光剖面50 m以下的底层作为代表性的近地面。图6(a)和(b)显示了MAX-DOAS反演的AOD和近地面气溶胶消光系数与用地面站点测量的地表PM2.5浓度的时间序列图，其中图中橘色阴影部分的AOD值较高，而气溶胶近地面消光系数和PM2.5浓度均较低。具体分析橘色阴影时段(12月3日至6日)的气溶胶垂直廓线，将3 km以下的气溶胶消光系数分成6段进行统计分析，如图6(c)所示，气溶胶主要分布在0.8～1.4 km处，表明这一时段气溶胶主要位于高层，因此导致AOD与近地面仪器测量的PM2.5浓度相差较大。

如图6(d)(e)所示，MAX-DOAS反演的气溶胶近地面消光系数与地面PM2.5浓度之间存在良好的一致性，这意味着MAX-DOAS反演获得的近地面消光系数真实地反映了地面附近的颗粒物含量。气溶胶消光系数和颗粒质量浓度的日均值和小时均值线性拟合相关系数r分别为0.83和0.62。通过对两种气溶胶类型进行回归分析，气溶胶消光系数基本随颗粒质量浓度变化; 然而，它也受到粒子组成、吸湿性和气象条件等多种因素的影响；其次，原位PM2.5数据直接在地表面处或靠近表面测量，而气溶胶近地面消光系数是从MAX-DOAS反演到的从地面到50 m的剖面中提取的，在大的水平距离上取平均值。因此，颗粒的分布不均匀在垂直和水平方向上，降低了这两个数据系列的一致性[18]。同时，整个反演的数据集没有考虑云的影响。

图6 (a)MAX-DOAS反演的AOD与地面站点测得的PM2.5浓度的时间序列图；(b)MAX-DOAS反演气溶胶近地面消光系数与地面站点测得的PM2.5浓度的时间序列图；(c)橘色阴影时段3 km以下气溶胶消光系数分段统计图；(d)气溶胶近地面消光系数与PM2.5浓度的日均值相关性; (e)小时均值相关性

Fig.6(a)TimeseriesoftheAODsat360nmretrievedfromMAXDOASandcomparisonwithinsituPM2.5concentration；(b)TimeseriesofsurfaceaerosolextinctioncoefficientretrievedfromMAX-DOASandcomparisonwithinsituPM2.5concentration；(c)Segmentationchartofaerosolextinctioncoefficientbelow3kminorangeshadedperiod;Linearregressionplotsofdailyandhourlydataareshownin(d)and(e),respectively

3.3 一次雾霾过程分析

图6(a)和(b)的蓝色阴影处是2017年12月3日至6日的一次明显雾霾过程，图7(a)是MAX-DOAS测得的AOD与地面站点测得的PM2.5浓度在此次污染过程的时间序列图，4日的气溶胶含量达到最大值，5日稍有下降。图7(b)和(c)展现了这一段时间的温湿度及风速风向等气象参数，4日和5日的相对湿度基本在80%以上，这种高湿度环境更利于气溶胶吸湿增长，使污染加重；4日凌晨左右风向从东南风变成了北风，而且污染严重的两天风向均是来自西北方向，表明北边有污染气团传输过来，5日夜间风向又变为东南方向，气溶胶含量下降，污染减弱。图7(d)是基于后向轨迹分析(Hysplit模型，http：//www.ready.noaa.gov)的4日12点的48小时风场后向轨迹图，500 m高度风场来自东北方向，1 500和3 000 m高度风场来自西北方向，三者占比相近，可见西北方向的风场影响更大，这也与图7(c)结果吻合。从图8气溶胶消光系数垂直分布可以看出，气溶胶主要聚集在地面上1 km附近，且在雾霾期间(4日至5日)，气溶胶均是上午聚集在1公里高度，到了中午便开始沉降至近地面。另外，结合图6(b)这一时段气溶胶近地面消光系数与PM2.5趋势一致，更加说明MAX-DOAS测量对对流层低层的气溶胶最敏感。

图7 (a)2017年12月3日至6日MAX-DOAS反演的AOD与地面站点测量的PM2.5浓度时间序列图；(b)温湿度；(c)风速风向；(d)2017年12月4日12点的48 h风场后向轨迹图

Fig.7(a)TimeseriesoftheAODsat360nmretrievedfromMAXDOASandcomparisonwithinsituPM2.5concentrationon3thto6thDec2017；Diurnalvariationofmeteorologicalparametersareshowedin(b)and(c)；(d)49hback-trajectoriesarrivinginHefeiat12:004Dec2017

图8 气溶胶消光系数的垂直分布

4 结 论

利用MAX-DOAS从2017年12月至2018年1月在合肥市科学岛观测并反演了气溶胶AOD和消光系数。通过反演O4斜柱浓度，利用廓线反演算法PriAM反演对流层4 km以下的气溶胶消光系数。将MAX-DOAS在360 nm反演的AOD与太阳光度计CE318的结果进行比较，二者均表现出相似的时间变化规律，小时均值和日均值的相关系数r均达到0.90以上，但太阳光度计测得的AOD比MAX-DOAS测得的高出20%。MAX-DOAS测量对对流层低层的气溶胶最敏感，而太阳光度计通过直射太阳光测量得到AOD，对整个大气中的气溶胶很敏感。将MAX-DOAS反演的近地面气溶胶消光系数与地面站点的PM2.5浓度进行了对比，两种数据之间的对比呈现出良好一致性。选择2017年12月3日至6日的雾霾期间MAX-DOAS反演气溶胶，与地面站点PM2.5对比，并通过气团后向轨迹分析辅助MAX-DOAS反演识别污染来源，结果表明，除近地面二者表现出相同的变化趋势外，MAX-DOAS还获得了高空气溶胶的分布变化信息，说明MAX-DOAS技术是获取对流层的气溶胶信息的有效工具，其结果可以为气象、环境管理部门提供重要的数据支持。