环境科学  2022, Vol. 43 Issue (5): 2284-2293   PDF    
海-陆大气交汇作用下青岛冬季大气PM2.5污染特征与来源解析
庹雄, 杨凌霄, 张婉, 齐安安, 王浥铭, 王鹏程, 黄琦, 赵彤, 张雄飞, 徐鹏, 张天琪, 王文兴     
山东大学环境研究院, 青岛 266237
摘要: 为了研究海-陆大气交汇作用对沿海城市大气污染物的传输与扩散的影响, 于2019年11月18日至12月23日在青岛观测站采集大气PM2.5样品, 对PM2.5中的水溶性离子、无机元素和碳质组分特征进行了分析, 并结合后向轨迹聚类分析模型和PMF模型等方法对青岛市冬季大气污染来源进行分析.结果表明, 青岛冬季ρ(PM2.5)平均值为61.0 μg ·m-3, 其中, ρ(水溶性离子)、ρ(无机元素)、ρ(OC)和ρ(EC)平均值分别为29.9、5.46、10.2和3.82 μg ·m-3; 二次离子(SO42-、NO3-、NH4+)和地壳元素(Si、K、Ca、Fe、Ti)是主要的离子和元素成分, 分别占水溶性离子和无机元素的89.3%和61.1%.青岛市大气主要受局部海陆风气流影响(43.4%), 其次是季风气流(36.2%), 冷空气气流对青岛影响较低(20.4%).PMF模型结果表明, 青岛冬季PM2.5主要来自机动车源(22.0%)、海盐源(21.0%)和煤燃烧源(19.8%).局部海陆风气流中的污染物主要来自煤燃烧(25.4%)和机动车(18.9%), 该气流会加剧本地污染向市内聚集; 冷空气气流则以海盐源为主(26.8%), 有利于市内大气污染扩散; 季风气流对青岛无机元素贡献较高, 也会将外地的机动车污染输送至青岛.
关键词: 海-陆大气交汇      PM2.5      水溶性离子      无机元素      碳质组分      正定矩阵因子分解模型(PMF)     
Characteristics and Source Analysis of PM2.5 in Qingdao in Winter Under the Action of Sea-Land-Atmosphere Convergence
TUO Xiong , YANG Ling-xiao , ZHANG Wan , QI An-an , WANG Yi-ming , WANG Peng-cheng , HUANG Qi , ZHAO Tong , ZHANG Xiong-fei , XU Peng , ZHANG Tian-qi , WANG Wen-xing     
Environment Research Institute, Shandong University, Qingdao 266237, China
Abstract: In order to explore the impact of sea-land-atmosphere convergence on the transmission and diffusion of atmospheric pollutants in coastal cities, PM2.5 samples were collected from November 18 to December 23, 2019, at the Qingdao Observatory.The characteristics of water-soluble inorganic ions, inorganic elements, and carbon components in PM2.5 were analyzed, and the sources of air pollution in Qingdao in winter were analyzed via the backward trajectory cluster analysis model and PMF. The results showed that the average concentration of PM2.5 in Qingdao in winter was 61.0 μg ·m-3. The average concentrations of water-soluble ions, inorganic elements, OC, and EC were 29.9, 5.46, 10.2, and 3.82 μg ·m-3, respectively. The secondary ions (SO42-, NO3, and NH4+) and crustal elements (Si, K, Ca, Fe, and Ti) were the main ions and elemental components, accounting for 89.3% and 61.1% of water-soluble ions and inorganic elements, respectively. Atmosphere of Qingdao was mainly affected by the local sea-land airflow (43.4%), followed by the temperate monsoon airflow (36.2%). The impact of cold air airflow on Qingdao was relatively low (20.4%) The results of the PMF model indicated that the main sources of PM2.5 in Qingdao in winter were motor vehicle sources (22.0%), sea salt sources (21.0%), and coal combustion sources (19.8%). The pollutants in the local sea-land airflow mainly came from coal combustion (25.4%) and motor vehicles (18.9%), and this airflow aggravated local pollution in the city. The cold air airflow was dominated by sea salt sources (26.8%), which was conducive to the spread of air pollution in the city. The monsoon airflow had a high contribution to Qingdao's inorganic elements, and it transported the pollution of motor vehicles from other places to Qingdao.
Key words: sea-land-atmosphere convergence      PM2.5      water-soluble inorganic ions      inorganic elements      carbon components      positive matrix factorization model(PMF)     

随着我国大气污染防控政策的深入实施, 大气PM2.5污染整体减轻, 但是仍有较多城市的PM2.5浓度超过《环境空气质量标准》中的二级标准, 并且形成PM2.5与臭氧的复合污染[1]. PM2.5化学组成多样, 来源和成因复杂, 对环境质量和人体健康都有着较高的危害[2~4].水溶性无机离子(WSIIs)作为PM2.5的主要成分, 通常占比达到30% ~80%[5, 6], 其中二次离子(SNA, 即NH4+、NO3-和SO42-)不仅是WSIIs的主要成分[7], 其与臭氧也有着较强的依存关系, 更是PM2.5和臭氧协同控制的重点[1].无机元素虽然在PM2.5中占比较低, 但其易附着于细颗粒物中, 且金属元素对人体的危害性较大[2]; 无机元素还具有独特的源特征, 常被用作颗粒物的示踪成分[8]. PM2.5中的碳组分对光有散射和吸收作用, 能明显影响大气能见度[9].其中无机碳(EC)主要来自化石或生物质燃料的不完全燃烧, 其稳定性较强, 可作为一次污染源的标志物[10]; 有机碳(OC)成分和来源都较复杂, 其对人体的危害性也较高[11, 12], 其中二次有机碳(SOC)更是被认定为引发重污染的关键组分[13].因此研究大气PM2.5的化学组分与污染特征, 不但能有效探究其来源, 也有利于控制PM2.5污染, 降低对环境和人体的危害.

随着经济发展和城市化进程, 城市污染现象逐渐增加, 较多的研究以城市局部环境气象要素为主, 进行PM2.5化学组分与来源解析分析[14~16].有研究发现, 海陆风能缓解京津冀地区的大气污染, 而山谷风的输送和汇聚作用会增加北京及北京以南地区和太行山东部地区的大气重污染[17].邵玄逸等[18]的研究发现, 西北和西南通道的大气传输明显增加了京津冀地区的污染, 外来源对北京和唐山的贡献达到48.7%和30.7%.长江三角洲地区各城市之间的污染传输比较明显, 城市群污染来源复杂, 影响因素多样[19], 余钟奇等[20]的研究发现秋冬季的长江三角洲内部排放贡献达到84.1%; 二次转化对长江三角洲地区PM2.5有着显著贡献[19].而青岛位于山东半岛的东南部, 东、南皆与黄海相接, 是一座典型的滨海城市, 也是我国华东地区大气和海洋交汇的一个重要通道, 其气象条件不仅受区域性温带季风气候的影响, 还受局部海陆风环流的影响; 同时, 屈文军等[21]的研究发现局部海陆风中的陆风输送也会增加市内PM2.5的质量浓度, 这种现象在冬季影响更为明显.因此分析青岛市大气PM2.5污染情况, 有利于研究海-陆大气交汇作用对沿海城市PM2.5的影响.

受体模型是常见的源解析方法, 其中化学质量平衡模型(chemical mass balance model, CMB)需要使用具体排放源成分谱协助分析, 主成分分析法(principal components analysis, PCA)无法直接定量解析[4]; 而正定矩阵因子分解模型(positive matrix factorization model, PMF)仅需受体点成分谱信息, 就可定量解析出来源贡献, 因此近年研究多以PMF模型进行PM2.5的来源定量分析[2, 4, 22].混合型单粒子拉格朗日综合轨迹(hybrid single particle lagrangian integrated trajectory, HYSPLIT)能模拟大气后向轨迹, 是研究气流轨迹的主要工具[3, 23], 结合PMF模型可以分析出每条轨迹的污染源, 从而判断出大气交汇作用对PM2.5污染的传输影响.为了深入研究青岛冬季PM2.5的污染特征, 本文使用HYSPLIT模型和PMF分析青岛冬季PM2.5的主要来源, 重点分析海-陆大气交汇作用下不同气流轨迹对青岛冬季大气PM2.5的污染特征与来源的影响, 以期为确定沿海区域大气污染的基本特征及其潜在来源提供参考.

1 材料与方法 1.1 样品采集

本研究于2019年11月18日至12月23日在青岛即墨区超级观测站(120.67°E, 36.35°N, 距地面高15 m), 使用2台武汉天虹仪器仪表有限公司生产的智能中流量悬浮颗粒物采样器(型号:TH-150A, 流量100 L ·min-1)连续采集23.5 h的PM2.5样品, 一台使用特氟龙膜(What man, UK)收集样品, 并用于分析水溶性离子和无机元素; 另一台则用石英膜(Pall Gelman, USA)收集PM2.5并用于碳质组分分析.样品采集之前, 石英膜需在600℃的马弗炉中煅烧6 h来去除膜上残留的有机物对实验的干扰; 将采样后的石英膜保存好并放置在-20℃的冰箱中, 以防有机物挥发.采样前的特氟龙膜一直保存于室温下的干燥皿中, 采样前后用北京赛多利斯科技仪器有限公司的电子天平(型号SOP, 精确度0.01 mg)称重.每个膜需平衡24 h后再称重, 并且每次称重要保证连续3次重量之差小于0.05 mg.为防止样品交叉污染, 每次采样后对切割头用酒精棉球进行擦拭.根据已有研究[24~26], 本研究采样期间的11月18日和19日有明显的沙尘天气现象.

1.2 分析方法

本研究使用离子色谱(ICS-600, 美国, 阴离子淋洗液:20 mmol ·L-1 KOH; 阳离子淋洗液:20 mmol ·L-1甲基磺酸)、Epsilon 4能量色散X射线荧光光谱仪(ED-XRF, PANalytical, 荷兰)分析特氟龙膜中的水溶性离子和无机元素; 使用OC/EC分析仪(Sunset Laboratory Inc, 美国)分析石英膜的EC与OC.具体仪器原理、分析方法和质控见文献[27].水溶性离子的仪器检出限在0.01~0.022 μg ·m-3之间; 碳质组分的仪器检出限均为0.20 μg ·m-3; Si、Mg和Al的仪器检出限在0.2~0.4 μg ·m-3之间, 其余无机元素的均在0.001~0.015 μg ·m-3.

所有样品测试每测10个样品后随机进一个重复样, 离子和碳质组分保证样品两次测量的相对标准偏差低于5%, 无机元素的标准偏差低于10%.所有测试均有实验空白和野外空白, 空白样品处理与样品处理完全相同, 所有离子、元素和碳质组分的浓度数据均减去空白样品的浓度, 以排除野外空白和空白膜对实验数据的干扰.

1.3 后向轨迹聚类分析

本文采用已广泛用于大气污染传输轨迹计算的HYSPLIT模型[28], 以观测站(120.67°E, 36.35°N)为中心, 起始高度设置为50 m(当地平均海拔与距地面距离之和), 以一日4次(00:00、06:00、12:00和18:00)进行后向时间尺度为48h的轨迹计算, 并进行聚类分析.

1.4 PMF分析

正定矩阵因子分解法(PMF)最早是由Paatero等[29]提出的基于最小二乘法的定量受体源解析模型.PMF模型将样品数据分解为源成分谱矩阵F(k×j)和源贡献矩阵G(i×k), 以及一个残差矩阵E(i×j).

(1)

式中, Xij是第i个样品的第j个化学成分的浓度, Fkj是源k中第j个化学成分浓度, Gik是源k对第i个样品的贡献率, p为源的数量.

PMF模型残差矩阵与输入的不确定度矩阵比值Q的最小化为模型最优解.

(2)

式中, eijuij分别为残差矩阵与不确定度矩阵.

当样品浓度小于等于相应方法检出限(MDL)时, uij为:

(3)

当大于MDL时, uij为:

(4)

式中, σ为样品误差, 取经验值15%; c为样品浓度.

本文利用美国环保署(USEPA)推荐的PMF 5.0软件对PM2.5的相关成分进行源解析.

1.5 气象数据来源

此超级观测站使用美国赛默飞世尔公司的5014i、49i、42i、43i和48i型号仪器在线观测PM2.5、O3、NOx、SO2和CO; 同时利用气象站(WS-6P, 智翔宇仪器)在线监测风向和风速等数据, 其观测能力达到1 min的分辨率.再利用美国气象环境预报中心(NCEP)和美国国家大气研究中心(NCAR)联合制作NCEP再分析数据集, 其采用了当今最先进的全球资料同化系统和完善的数据库, 对各种资料来源(地面、船舶、无线电探空和卫星等)的观测资料进行质量控制和同化处理, 它不仅包含的要素多, 范围广, 而且延伸的时段长, 能够准确分析出大气气流轨迹.

2 结果与讨论 2.1 后向轨迹聚类分析

为了有效揭示海-陆大气交汇作用对青岛冬季PM2.5的影响, 本文利用HYSPLIT模型并结合NCEP再分析数据集, 对采样期间青岛的气流进行48 h后推轨迹计算并聚类分析, 如图 1(a)所示; 并通过超站观测到的风速风向等在线数据, 绘制采样期间风玫瑰图[图 1(b)].

图 1 研究期间青岛市气流后向轨迹聚类与风玫瑰图 Fig. 1 Backward trajectory clustering and wind rose diagram of Qingdao during the study period

图 1中气流后向轨迹聚类可知, 采样期间的48 h气流轨迹可分为3类, 表现为轨迹1(短轨迹)和轨迹2、3(长轨迹), 其气流占比顺序为:轨迹1>轨迹3>轨迹2.根据风玫瑰图可知, 采样期间盛行西北风, 其次是南风, 而采样点南部靠海, 南风主要表现为海风, 西北风表现为陆地风, 这与颜梅[30]根据百年青岛气场和风场的数据统计的结果相似.其中短轨迹1呈现环形, 并且以采样点南部为主, 且气流传输高度也较低, 符合局部海陆风的低气层环流特征, 所以推测短轨迹1为局部海陆风形成的气流轨迹.轨迹1是本次采样期间的后向聚类轨迹中占比最高, 达到43.4%, 说明局部海陆风气流对青岛的大气有着较大的影响, 可能会促进青岛市大气污染物向市内聚集, 使其不易扩散[31].

轨迹2、3为长轨迹, 均来自采样点西北方向, 表现为温带季风气候影响的气流轨迹, 但轨迹2、3在气流起源和途径上都有着一定的差异.轨迹2起源于俄罗斯, 途经蒙古国传入我国内蒙古自治区, 并经过我国辽宁省跨渤海传到青岛, 该轨迹在前期主要为西北风, 传至辽宁省时, 由于受西伯利亚冷空气影响[32, 33], 该气流轨迹转而南下, 途径渤海传至青岛.该轨迹传入青岛前主要受冷空气影响, 并以北风形式传入青岛, 故将该轨迹视为冷空气气流轨迹, 而其又经过渤海的净化作用, 可能有助于青岛市内污染物的扩散.

长轨迹3起源于我国内蒙古自治区和蒙古国边境地区, 途经河北省, 也是跨渤海传至山东, 但该轨迹途经河北省唐山和山东省淄博、潍坊等重工业区, 容易携带较多的大气污染物, 并以西北风的形式传至青岛, 因此将该轨迹视为季风气流轨迹, 并且轨迹3可能会增加外来污染源对青岛市大气污染的贡献.

2.2 各轨迹PM2.5组分分析

图 2所示, 不同轨迹上的PM2.5浓度平均值大小顺序为:轨迹1>轨迹3>轨迹2; 同时, 根据表 1可知, 轨迹1的PM2.5浓度平均值高于《环境空气质量标准》中的二级标准, 说明轨迹1存在较高的污染, 尤其在12月7~10日, 青岛出现高浓度的PM2.5污染, 而此时青岛市气流均为轨迹1.主要原因是此时冷空气较弱, 青岛以局部海陆风气流为主, 形成了低空低温层结[33], 导致市内污染不易扩散, 从而促使青岛市内PM2.5污染加剧.而轨迹2的PM2.5浓度明显低于青岛冬季的平均值, 表明轨迹2对青岛有一定的净化作用.轨迹1对青岛市冬季PM2.5、水溶性离子、OC和EC都有着较高的贡献, 轨迹3对无机元素有着较高的贡献, 而轨迹2对青岛市PM2.5及其各组分的贡献均最低.综上所述, 轨迹1主要表现为污染轨迹, 说明海陆风气流对青岛市PM2.5及其各组分有着较高的影响; 轨迹2为清洁轨迹, 说明冷空气能有效缓解青岛市的污染情况; 轨迹3对青岛市的无机元素贡献较高, 可能与沙尘天气有关.

图 2 不同轨迹PM2.5的组分分析 Fig. 2 Composition analysis of PM2.5 in different trajectories

表 1 各轨迹PM2.5成分/μg ·m-3 Table 1 PM2.5 components of each trajectory/μg ·m-3

2.2.1 各轨迹PM2.5中OC和EC特征分析

PM2.5中的EC主要来自一次污染排放[2], 并且在大气中相对稳定; 而OC的来源分为直接排放的一次有机碳(POC)和在大气中二次反应生产的二次有机碳(SOC)[34].有研究表明[35], 当OC/EC>2时, 有SOC生产, 而在本研究中青岛冬季OC/EC的平均值为4.01, 说明青岛冬季OC中有明显的SOC生成.本研究采用常用的环境样本中的OC和EC最小比值法来估算PM2.5中的POC和SOC的量[13, 36], 其计算公式如下:

(5)
(6)

表 1可知, 轨迹1的OC和EC含量都是各轨迹中最高的, 但其SOC含量却是最低, 也仅占轨迹1中OC的26.7%;而轨迹2和轨迹3的SOC分别占OC的56.6%和52.0%.表明轨迹2和轨迹3的传输距离长, 污染物容易二次转化, 因此轨迹2和轨迹3中的OC主要来自二次污染源; 轨迹1传输距离短, 污染物不易二次转化, OC主要受一次污染源影响.

根据相关研究, OC/EC可以用来判断颗粒物的来源, 当OC/EC为1.0~4.2时, OC和EC主要来源是机动车的排放[37]; 当OC/EC为2.5~10.5时, OC和EC就是典型的煤燃烧所排放的[38]; 而生物质燃烧排放的OC和EC, 其OC/EC为4.1~14.5[39].在本研究期间,轨迹1~3的OC/EC平均值分别为3.86、4.42和3.75, 表明轨迹1和轨迹3都主要来自煤燃烧和机动车尾气; 而轨迹2主要来自煤燃烧和生物质燃烧.

2.2.2 各轨迹PM2.5中水溶性离子特征分析

轨迹1~3的ρ(WSII)平均值分别为54.4、11.6和22.1 μg ·m-3, 其中SNA质量分数分别达到92.0%、81.6%和86.2%(图 3).通过分析NO3-和SO42-的比值, 可以初步判断移动源和固定源对颗粒物的影响[2].轨迹1~3中两种离子比值均值分别为2.37、1.73和2.48, 均大于1, 说明各轨迹上的移动源对颗粒物影响均要大于固定源, 也表明青岛市冬季大气污染主要受机动车尾气影响.

图 3 各轨迹WSII的组分质量分数 Fig. 3 Mass fraction of WSII components in each trajectory

大气中的SO42-和NO3-主要来自大气中SO2和NO3的气相氧化和云内反应[7], 而硫氧化率(SOR)和氮氧化率(NOR)可以表示这种转化程度, 也有利于研究青岛市冬季大气气溶胶的二次转化, 具体计算公式如下:

(7)
(8)

当SOR和NOR大于0.1时, 说明大气中有明显的SO2和NO3向SO42-和NO3-转化, 而其值越大, 则说明转化程度越高[40]. 根据超站在线观测到的SO2和NO2数据, 计算出轨迹1~3的SOR分别为0.402、0.214和0.228, NOR分别为0.324、0.133和0.190; 说明各轨迹均有明显的转化, 其中轨迹1虽然为短轨迹, 但其SOR和NOR均为最高, 说明轨迹1的前体物转化程度最高, 可能与轨迹1中前体物浓度较高有关.由表 2可知, 虽然NO2浓度明显高于SO2, 但各轨迹中的SOR均高于NOR, 说明在同样的大气环境下, 冬季SO2更容易向二次污染物转化.

表 2 各轨迹SO42-和NO3-分析 Table 2 Analysis of SO42- and NO3- for each trajectory

2.2.3 各轨迹PM2.5中无机元素特征分析

轨迹1~3的无机元素浓度平均值分别为4.98、4.77和6.57 μg ·m-3; 均以S、K、Fe、Cl、Si和Ca等10种元素为主(图 4), 分别占到总无机元素的93.3%、95.2%和95.6%.其中轨迹2和轨迹3以Si和Ca等地壳元素为主[4], 表明长距离传输更容易输入地壳元素.轨迹1中S的占比最高, 而S作为煤燃烧标志[41], 说明轨迹1主要通过局部海陆风气流将本地煤燃烧产生的污染物传入青岛.

图 4 各轨迹中无机元素质量分数 Fig. 4 Mass fraction of inorganic elements in each trajectory

2.2.4 PMF源解析

为了更准确地判断出青岛冬季PM2.5的主要来源, 本文利用美国环保署(USEPA)推荐的PMF 5.0软件对PM2.5中化学组分进行源解析.由于无机元素的种类较多, 故本文选取10种具有特征的代表元素, 与水溶性离子和OC、EC进行PMF模型分析.按3~7个因子数分别进行20次的基础运行, 对各因子数下运行的结果进行质量平衡约束.以上运行均通过了误差及DISP检验.最后, 根据每个因子的成分图谱, 对照已有研究中排放源指纹谱或标志物, 结合采样点周边实际排放源信息确定了6个因子, 分别为扬尘源、海盐源、二次生成源、生物质燃烧源、机动车源和煤燃烧源.

图 5所示, 因子1被认定为扬尘源, 因为典型地壳元素Ti、Si、Al和Ca对因子1有着显著贡献[4, 42]; 扬尘标志物Ca2+对因子1也有着较高的贡献[43, 44].因子2主要特征组分为Na+、Cl-和Mg2+, Na+和Cl-是海盐源的标志物[43]; 海盐源也是Mg2+的重要来源之一[27], 故将因子2判定为海盐源.因子3中以NH4+和NO3-为主, SO42-对其贡献也较高, 具有明显的二次源特征[4], 故将因子3解释为二次生成源.因子4中EC的负荷最高, 达到59.9%, EC主要来自生物质燃烧[45], 因此将因子4视为生物质燃烧源.因子5中, 贡献率较高的物种为Cu、Zn、Fe和Pb, 而机动车润滑油、发动机磨损、刹车片和轮胎磨损等都会增加这些元素的浓度[3, 8, 15, 42], 因子5可被认定为机动车源.因子6主要以SO42-和S为主, S主要来自含硫煤的燃烧[15, 42], 煤燃烧对SO42-也有着较高的贡献[44], 故可将因子6可认定为煤燃烧源.

图 5 青岛冬季细颗粒物PMF解析结果 Fig. 5 PMF analysis results of fine particulate matter in Qingdao in winter

通过PMF计算, 采样期间的PM2.5主要来自机动车源(22.0%)、海盐源(21.0%)和煤燃烧源(19.8%); 为了探究大气作用对青岛市冬季的影响, 本研究分析了各簇气流轨迹中的来源贡献率, 结果如图 6所示.轨迹1以煤燃烧源为主, 二次生成源和机动车源对其也有着较高的贡献, 说明局部海陆风气流中的污染主要来自青岛市内, 通过海风将煤燃烧和机动车等产生的污染物输送回市内.海盐源和机动车源对轨迹2贡献显著, 其中海盐源可能是该轨迹途经渤海时携带而至; 机动车源可能是青岛市本地产生的.轨迹3和轨迹2的污染来源相似, 以机动车和海盐源为主, 其海盐源也主要来自渤海; 而机动车源不仅与青岛市有关, 也可能与淄博等拥有大量机动车的工业城市有关, 其产生的污染物也可通过轨迹3传至青岛.相比于轨迹1和轨迹2, 轨迹3中的扬尘源有着明显的增加, 可能与沙尘天有关, 说明沙尘主要通过轨迹3传至青岛.

图 6 PMF解析结果对PM2.5的贡献率 Fig. 6 Contribution rate of PMF analysis results to PM2.5

3 结论

(1) 青岛冬季ρ(PM2.5)平均值为61.0 μg ·m-3, 其中, ρ(水溶性离子)、ρ(无机元素)、ρ(OC)和ρ(EC)平均值分别为29.9、5.46、10.2和3.82 μg ·m-3; 青岛市大气主要受局部海陆风气流影响(43.4%), 其次是季风气流(36.2%); 冷空气气流对青岛影响较低(20.4%), 但能改善青岛市的大气质量.

(2) 组分分析结果表明, 局部海陆风气流对青岛市冬季PM2.5、水溶性离子、OC和EC都有着较高的贡献; 局部海陆风气流、冷空气气流和季风气流中的SOR分别达到0.402、0.214和0.228, NOR分别为0.324、0.133和0.190; SNA占比为92.0%、81.6%和86.2%, 都有明显的二次转化存在.

(3) 源解析结果表明, 局部海陆风气流中的污染物主要来自煤燃烧(25.4%)和机动车(18.9%), 该气流会加剧本地污染向市内聚集, 从而导致青岛市内雾-霾天的形成.冷空气气流则以海盐源为主(26.8%), 有利于市内大气污染扩散.季风气流对青岛市大气环境影响较低, 主要因沙尘因素而增加对青岛市无机元素的贡献, 也会将外地的机动车污染输送至青岛.

(4) 青岛冬季受移动源影响较大, 各轨迹中的NO3-和SO42-的比值均大于1; 机动车尾气对青岛市冬季大气污染的贡献率达到22.0%, 在局部海陆风气流、冷空气气流和季风气流中的贡献分别达到18.9%、25.8%和23.2%.故本文建议青岛市应该加强对冬季机动车的控制, 限制私家车, 鼓励公共出行, 能有效减少青岛市冬季的大气污染并改善青岛市内大气质量.

参考文献
[1] 李红, 鲍捷萌, 毕方, 等. PM2.5与臭氧污染协同控制: 挑战与应对[J]. 世界环境, 2020(5): 24-29.
Li H, Bao J M, Bi F, et al. Synergistic control of PM2.5 and ozone pollution: Challenges and responses[J]. World Environment, 2020(5): 24-29.
[2] 赵清, 李杏茹, 王国选, 等. 运城秋冬季大气细粒子化学组成特征及来源解析[J]. 环境科学, 2021, 42(4): 1626-1635.
Zhao Q, Li X R, Wang G X, et al. Chemical composition and source analysis of PM2.5 in Yuncheng, Shanxi Province in autumn and winter[J]. Environmental Science, 2021, 42(4): 1626-1635.
[3] 张剑飞, 姜楠, 段时光, 等. 郑州市PM2.5化学组分的季节变化特征及来源解析[J]. 环境科学, 2020, 41(11): 4813-4824.
Zhang J F, Jiang N, Duan S G, et al. Seasonal chemical composition characteristics and source apportionment of PM2.5 in Zhengzhou[J]. Environmental Science, 2020, 41(11): 4813-4824. DOI:10.3969/j.issn.1000-6923.2020.11.021
[4] 王成, 闫雨龙, 谢凯, 等. 阳泉市秋冬季PM2.5化学组分及来源分析[J]. 环境科学, 2020, 41(3): 1036-1044.
Wang C, Yan Y L, Xie K, et al. Analysis of chemical components and sources of PM2.5 during autumn and winter in Yangquan city[J]. Environmental Science, 2020, 41(3): 1036-1044.
[5] Wang H L, Zhu B, Shen L J, et al. Water-soluble ions in atmospheric aerosols measured in five sites in the Yangtze River Delta, China: Size-fractionated, seasonal variations and sources[J]. Atmospheric Environment, 2015, 123: 370-379. DOI:10.1016/j.atmosenv.2015.05.070
[6] Xu L L, Duan F K, He K B, et al. Characteristics of the secondary water-soluble ions in a typical autumn haze in Beijing[J]. Environmental Pollution, 2017, 227: 296-305. DOI:10.1016/j.envpol.2017.04.076
[7] 冯炎鹏, 张军科, 黄小娟, 等. 成都夏冬季PM2.5中水溶性无机离子污染特征[J]. 环境科学, 2020, 41(7): 3012-3020.
Feng Y P, Zhang J K, Huang X J, et al. Pollution characteristics of water-soluble inorganic ions in Chengdu in summer and winter[J]. Environmental Science, 2020, 41(7): 3012-3020.
[8] 周安琪, 刘建伟, 周旭, 等. 北京大气PM2.5载带金属浓度、来源及健康风险的城郊差异[J]. 环境科学, 2021, 42(6): 2595-2603.
Zhou A Q, Liu J W, Zhou X, et al. Concentrations, sources, and health risks of PM2.5 carrier metals in the Beijing urban area and suburbs[J]. Environmental Science, 2021, 42(6): 2595-2603. DOI:10.3969/j.issn.1000-6923.2021.06.012
[9] 李安娜, 温天雪, 华维, 等. 鼎湖山大气颗粒物中OC与EC的浓度特征及粒径分布[J]. 环境科学, 2020, 41(9): 3908-3917.
Li A N, Wen T X, Hua W, et al. Characterization and size distribution of carbonaceous aerosols at Mountain Dinghu[J]. Environmental Science, 2020, 41(9): 3908-3917.
[10] Chen P F, Kang S C, Tripathee L, et al. Light absorption properties of elemental carbon (EC) and water-soluble brown carbon (WS-BrC) in the Kathmandu Valley, Nepal: A 5-year study[J]. Environmental Pollution, 2020, 261. DOI:10.1016/j.envpol.2020.114239
[11] Wu X F, Chen C R, Vu T V, et al. Source apportionment of fine organic carbon (OC) using receptor modelling at a rural site of Beijing: Insight into seasonal and diurnal variation of source contributions[J]. Environmental Pollution, 2020, 266. DOI:10.1016/j.envpol.2020.115078
[12] 王成, 曹靖原, 段小琳, 等. 山西省四城市冬季PM2.5中碳质组分特征及来源分析[J]. 环境工程, 2021, 39(6): 114-121.
Wang C, Cao J Y, Duan X L, et al. Characteristics and sources analysis of carbonaceous components in PM2.5 in winter in four cities of Shanxi Province[J]. Environmental Engineering, 2021, 39(6): 114-121.
[13] 董贵明, 唐贵谦, 张军科, 等. 北京南部城区PM2.5中碳质组分特征[J]. 环境科学, 2020, 41(10): 4374-4381.
Dong G M, Tang G Q, Zhang J K, et al. Characteristics of carbonaceous species in PM2.5 in southern Beijing[J]. Environmental Science, 2020, 41(10): 4374-4381.
[14] 王德羿, 王体健, 韩军彩, 等. "2+26"城市大气重污染下PM2.5来源解析[J]. 中国环境科学, 2020, 40(1): 92-99.
Wang D Y, Wang T J, Han J C, et al. Source apportionment of PM2.5 under heavy air pollution conditions in "2+26" cities[J]. China Environmental Science, 2020, 40(1): 92-99. DOI:10.3969/j.issn.1000-6923.2020.01.010
[15] 刘可可, 张红, 刘桂建. 合肥市PM2.5和PM10中元素组成特征及重污染成因分析[J]. 环境科学, 2019, 40(8): 3415-3420.
Liu K K, Zhang H, Liu G J, et al. Elemental composition characteristics of PM2.5 and PM10, and heavy pollution analysis in Hefei[J]. Environmental Science, 2019, 40(8): 3415-3420.
[16] 赵孝囡, 王申博, 杨洁茹, 等. 郑州市PM2.5组分、来源及其演变特征[J]. 环境科学, 2021, 42(8): 3633-3643.
Zhao X N, Wang S B, Yang J R, et al. Chemical components and sources of PM2.5 and their evolutive characteristics in Zhengzhou[J]. Environmental Science, 2021, 42(8): 3633-3643.
[17] 李青春, 李炬, 郑祚芳, 等. 冬季山谷风和海陆风对京津冀地区大气污染分布的影响[J]. 环境科学, 2019, 40(2): 513-524.
Li Q C, Li J, Zheng Z F, et al. Influence of mountain valley breeze and sea land breeze in winter on distribution of air pollutants in Beijing-Tianjin-Hebei region[J]. Environmental Science, 2019, 40(2): 513-524.
[18] 邵玄逸, 王晓琦, 钟嶷盛, 等. 京津冀典型城市冬季人为源减排与气象条件对PM2.5污染影响[J]. 环境科学, 2021, 42(9): 4095-4103.
Shao X Y, Wang X Q, Zhong Y S, et al. Impacts of anthropogenic emission reduction and meteorological conditions on PM2.5 pollution in typical cities of Beijing-Tianjin-Hebei in Winter[J]. Environmental Science, 2021, 42(9): 4095-4103.
[19] 郭振东, 朱彬, 王红磊, 等. 长江三角洲霾天气PM2.5中水溶性离子特征及来源解析[J]. 中国环境科学, 2019, 39(3): 928-938.
Guo Z D, Zhu B, Wang H L, et al. Characteristics and source analysis of water-soluble ions in PM2.5 in the haze weather over in Yangtze River Delta[J]. China Environmental Science, 2019, 39(3): 928-938. DOI:10.3969/j.issn.1000-6923.2019.03.004
[20] 余钟奇, 瞿元昊, 周广强, 等. 2018年秋冬季长江三角洲区域PM2.5污染来源数值研究[J]. 中国环境科学, 2020, 40(10): 4237-4246.
Yu Z Q, Qu Y H, Zhou G Q, et al. Numerical simulations of PM2.5 pollution source in the Yangtze River Delta region in fall and winter in 2018[J]. China Environmental Science, 2020, 40(10): 4237-4246. DOI:10.3969/j.issn.1000-6923.2020.10.006
[21] 屈文军, 张小曳, 盛立芳, 等. 青岛海滨大气PM2.5、PM10质量浓度与海陆风输送的关系[A]. 见: 中国气象学会2008年年会大气环境监测、预报与污染物控制分会场论文集[C]. 北京: 中国气象学会, 2008.
[22] 王国祯, 任万辉, 于兴娜, 等. 沈阳市冬季大气PM2.5中水溶性离子污染特征及来源解析[J]. 环境科学, 2021, 42(1): 30-37.
Wang G Z, Ren W H, Yu X N, et al. Characteristics and sources of water-soluble ion pollution in PM2.5 in winter in Shenyang[J]. Environmental Science, 2021, 42(1): 30-37.
[23] 孙佳傧, 董喆, 李利萍, 等. 洛阳市大气细颗粒物化学组分特征及溯源分析[J]. 环境科学, 2021, 42(12): 5624-5632.
Sun J B, Dong Z, Li L P, et al. Characteristics of chemical composition and source apportionment of atmospheric fine particulate matter in Luoyang[J]. Environmental Science, 2021, 42(12): 5624-5632.
[24] 中华人民共和国环境保护部. 关于印发《受沙尘天气过程影响城市空气质量评价补充规定》的通知[EB/OL]. https://www.mee.gov.cn/gkml/hbb/bgt/201701/t20170106_394054.htm, 2017-01-04.
[25] 胡元洁, 蒋楠. 沙尘天气对西安市环境空气质量的影响[J]. 中国沙漠, 2020, 40(6): 53-60.
[26] 孙韧, 陈冠益, 董海燕. 天津市春季颗粒物污染特征及典型沙尘过程分析[J]. 环境科学与技术, 2018, 41(6): 8-14.
Sun R, Chen G Y, Dong H Y. Pollution characteristics of particulate matters and a dust storm episode in the spring of Tianjin[J]. Environmental Science & Technology, 2018, 41(6): 8-14.
[27] 张艳. 滨海地区背景点和城区大气PM1. 和PM2.5的化学组成及对人体健康的影响[D]. 青岛: 山东大学, 2020.
[28] Wang G, Zhu Z Y, Zhao N, et al. Variations in characteristics and transport pathways of PM2.5 during heavy pollution episodes in 2013-2019 in Jinan, a central city in the north China Plain[J]. Environmental Pollution, 2021, 284. DOI:10.1016/j.envpol.2021.117450
[29] Paatero P, Tapper U. Analysis of different modes of factor analysis as least squares fit problems[J]. Chemometrics and Intelligent Laboratory Systems, 1993, 18(2): 183-194. DOI:10.1016/0169-7439(93)80055-M
[30] 颜梅, 徐勤爱, 邵滋和. 青岛百年气压与风场统计特征分析[J]. 海洋预报, 2003, 20(3): 12-17. DOI:10.3969/j.issn.1003-0239.2003.03.003
[31] 王宏, 郑秋萍, 洪有为, 等. 2017年漳州海陆风特征与冬春季污染物浓度关系[J]. 气象与环境学报, 2020, 36(3): 33-40.
Wang H, Zheng Q P, Hong Y W, et al. Characteristics of sea-land breeze and its relationship with pollutant concentrations during winter and spring in 2017 in Zhangzhou[J]. Journal of Meteorology and Environment, 2020, 36(3): 33-40.
[32] 张天航, 徐冉, 桂海林. 2019年11月大气环流和天气分析[J]. 气象, 2020, 46(2): 283-288.
Zhang T H, Xu R, Gui H L. Analysis of the November 2019 atmospheric circulation and weather[J]. Meteorological Monthly, 2020, 46(2): 283-288.
[33] 徐冉, 桂海林, 江琪, 等. 2019年12月大气环流和天气分析[J]. 气象, 2020, 46(3): 441-448.
Xu R, Gui H L, Jiang Q, et al. Analysis of the December 2019 atmospheric circulation and weather[J]. Meteorological Monthly, 2020, 46(3): 441-448.
[34] 康晖, 朱彬, 王红磊, 等. 长三角典型站点冬季大气PM2.5中OC、EC污染特征[J]. 环境科学, 2018, 39(3): 961-971.
Kang H, Zhu B, Wang H L, et al. Characterization and variation of organic carbon (OC) and elemental carbon (EC) in PM2.5 during the winter in the Yangtze River Delta Region, China[J]. Environmental Science, 2018, 39(3): 961-971.
[35] Chow J C, Watson J G, Lu Z Q, et al. Descriptive analysis of PM2.5 and PM10 at regionally representative locations during SJVAQS/AUSPEX[J]. Atmospheric Environment, 1996, 30(12): 2079-2112. DOI:10.1016/1352-2310(95)00402-5
[36] 武高峰, 王丽丽, 武志宏, 等. 石家庄市采暖季PM2.5碳组分昼夜污染特征及来源分析[J]. 环境科学学报, 2020, 40(7): 2356-2364.
Wu G F, Wang L L, Wu Z H, et al. Pollution diurnal variation and source analysis of PM2.5 carbon components in heating season in Shijiazhuang City[J]. Acta Scientiae Circumstantiae, 2020, 40(7): 2356-2364.
[37] Ding X X, Kong L D, Du C T, et al. Characteristics of size-resolved atmospheric inorganic and carbonaceous aerosols in urban Shanghai[J]. Atmospheric Environment, 2017, 167: 625-641.
[38] Zhang H Y, Cheng S Y, Li J B, et al. Investigating the aerosol mass and chemical components characteristics and feedback effects on the meteorological factors in the Beijing-Tianjin-Hebei region, China[J]. Environmental Pollution, 2019, 244: 495-502.
[39] Liang X X, Huang T, Lin S Y, et al. Chemical composition and source apportionment of PM1 and PM2.5 in a national coal chemical industrial base of the Golden Energy Triangle, Northwest China[J]. Science of the Total Environment, 2019, 659: 188-199.
[40] 孙有昌, 姜楠, 王申博, 等. 安阳市大气PM2.5中水溶性离子季节特征及来源解析[J]. 环境科学, 2020, 41(1): 75-81.
Sun Y C, Jiang N, Wang S B, et al. Seasonal characteristics and source analysis of water-soluble ions in PM2.5 of Anyang city[J]. Environmental Science, 2020, 41(1): 75-81.
[41] 刘威杰, 石明明, 程铖, 等. 夏季大气PM2.5中元素特征及源解析: 以华中地区平顶山-随州-武汉为例[J]. 环境科学, 2020, 41(1): 23-30.
Liu W J, Shi M M, Cheng C, et al. Characteristics and sources of elements in PM2.5 during summer for three typical cities in Pingdingshan-Suizhou-Wuhan, Central China[J]. Environmental Science, 2020, 41(1): 23-30.
[42] 李雪梅, 牟玲, 田妹, 等. 山西大学城PM2.5中元素特征、来源及健康风险评估[J]. 环境科学, 2020, 41(11): 4825-4831.
Li X M, Mu L, Tian M, et al. Characteristics, sources, and health risks of elements in PM2.5 in Shanxi University town[J]. Environmental Science, 2020, 41(11): 4825-4831.
[43] 李明燕, 杨文, 魏敏, 等. 典型沿海城市采暖期细颗粒物组分特征及来源解析[J]. 环境科学, 2020, 41(4): 1550-1560.
Li M Y, Yang W, Wei M, et al. Characteristics and sources apportionment of fine particulate matter in a typical coastal city during the heating period[J]. Environmental Science, 2020, 41(4): 1550-1560.
[44] 马红璐, 赵欣, 陆建刚, 等. 宿迁市PM2.5中水溶性无机离子的季节特征和来源分析[J]. 环境科学, 2020, 41(9): 3899-3907.
Ma H L, Zhao Xin, Lu J G, et al. Seasonal characteristics and source analysis of water-soluble inorganic ions in PM2.5 in Suqian city[J]. Environmental Science, 2020, 41(9): 3899-3907.
[45] 史芳天, 罗彬, 张巍, 等. 成都平原PM2.5中碳质组分时空分布特征与来源[J]. 环境科学, 2020, 41(2): 520-528.
Shi F T, Luo B, Zhang W, et al. Spatio-temporal variations and source apportionment of carbonaceous species in PM2.5 across multiple sampling locations in the Chengdu Plain[J]. Environmental Science, 2020, 41(2): 520-528.