环境科学  2023, Vol. 44 Issue (11): 5986-5996   PDF    
周口市2022年冬季重污染过程中细颗粒物污染特征及成因分析
马英歌1, 杨露1, 狄睿苗1, 马南2, 乔利平1, 吴宇航1, 周文鑫1, 赵新华2, 张越2, 孙志华2, 陈长虹1, 陈学军2, 楼晟荣1, 黄成1     
1. 上海市环境科学研究院国家环境保护城市大气复合污染成因与防治重点实验室, 上海 200233;
2. 河南省生态环境科学研究院, 郑州 450004
摘要: 为研究周口市冬季大气细颗粒物污染特征及导致其暴发增长的主要影响因素, 利用城区环境空气在线高时间分辨仪器对2022年1月周口市大气常规污染因子、细颗粒物中水溶性离子等进行观测分析. 结果表明, 二次无机气溶胶(SNA)、碳质气溶胶(CA, 包括有机碳OC和无机碳EC)以及重构后的地壳物质(CM, 如Al2O3、SiO2、CaO和Fe2O3等)是PM2.5中含量前三的组成, 占比分别为61.3%、24.3%和9.72%, SNA、CA、CM和二次有机气溶胶(SOA)浓度均随AQI升高而升高. 1月硫氧化率(SOR)和氮氧化率(NOR)分别为0.53和0.46, SO42-和NO3-增长速率[μg·(m3·h)-1]分别为0.027(-5.89~9.47, 下同)和0.051(-23.1~12.4).重污染时段SO42-和NO3-增长速率分别为0.13 μg·(m3·h)-1和0.24 μg·(m3·h)-1, 较月均值高4.8倍和4.7倍. 1月整体SOR大于NOR, 但由于气态前体物浓度和相对湿度变化等影响, NO3-增长速率约是SO42-的1.8倍, 重污染日NO3-增长速率显著高于SO42-. 高AQI和高湿度时段SOR、NOR、SNA和SOA浓度均高于低AQI和低湿时段, 大气氧化性Ox(NO2+O3)随湿度增加而下降. SOA夜间浓度高于白天, 夜间浓度随湿度增加幅度大于白天, SNA昼夜浓度差异不显著随湿度有缓慢增加.结合周口本地情况, 在低温、高湿和低风速条件下, 需重点关注SNA气态前体物排放, 加强对SO2和NO2主要排放源如移动源和燃煤源提前管控, 可降低冬季大气重污染峰值.
关键词: 细颗粒物(PM2.5)      二次无机气溶胶(SNA)      氮氧化率(NOR)      硫氧化率(SOR)      二次无机盐增长速率     
Characteristics and Cause of PM2.5 During Haze Pollution in Winter 2022 in Zhoukou, China
MA Ying-ge1 , YANG Lu1 , DI Rui-miao1 , MA Nan2 , QIAO Li-ping1 , WU Yu-hang1 , ZHOU Wen-xin1 , ZHAO Xin-hua2 , ZHANG Yue2 , SUN Zhi-hua2 , CHEN Chang-hong1 , CHEN Xue-jun2 , LOU Sheng-rong1 , HUANG Cheng1     
1. State Environmental Protection Key Laboratory of Formation and Prevention of Urban Air Pollution Complex, Shanghai Academy of Environmental Science, Shanghai 200233, China;
2. Henan Academy of Ecology and Environmental Science, Zhengzhou 450004, China
Abstract: The characteristics and main factors of causes of haze in Zhoukou in January 2022 were analyzed. Six air pollutants, water-soluble ions, elements, OC, EC, and other parameters in fine particulate matter were monitored and analyzed using a set of online high-time-resolution instruments in an urban area. The results showed that the secondary inorganic aerosols(SNA), carbonaceous aerosols(CA, including organic carbon OC and inorganic carbon EC), and reconstructed crustal materials(CM, such as Al2O3, SiO2, CaO, and Fe2O3, etc.) were the three main components, accounting for 61.3%, 24.3%, and 9.72% in PM2.5, respectively. The concentrations of SNA, CA, CM, and SOA were increased, accompanied with higher AQI. The sulfur oxidation rate(SOR) and nitrogen oxidation rate(NOR) in January were 0.53 and 0.46, respectively. The growth rates [μg·(m3·h)] of sulfate and nitrate were 0.027(-5.89-9.47, range) and 0.051(-23.1-12.4), respectively. During the haze period, the growth rates of sulfate and nitrate were 0.13 μg·(m3·h)-1and 0.24 μg·(m3·h)-1, which were 4.8 and 4.7 times higher than the average value of January, respectively. Although the sulfur oxidation rate was greater than the nitrogen oxidation rate, the growth rate of nitrate was approximately 1.8 times that of sulfate owing to the difference in the concentration of gaseous precursors and the influence of relative humidity. The growth rates of nitrate in SNA were significantly higher than those of sulfate on heavily polluted days. The values of SOR, NOR, and concentrations of SNA and SOA during higher AQI and humidity periods were higher than those in lower AQI and humidity periods. The Ox(NO2+O3) decreased with the increase in relative humidity. The SOA was higher at nighttime, increasing faster with the humidity than that in daytime. Under the situation of lower temperature, higher humidity, and lower wind speed, the emission of gaseous precursors of SNA requires further attention in Zhoukou in winter. Advanced control strategies of emissions of SO2 and NO2, such as mobile sources and coal-burning sources, could reduce the peak of haze in winter efficiently.
Key words: fine particulate matter(PM2.5)      secondary inorganic aerosol(SNA)      nitrogen oxidation rate(NOR)      sulfur oxidation rate(SOR)      SNA growth rate     

大气细颗粒物(PM2.5)中主要组分包括无机碳、有机碳、离子和元素等一次排放, 以及在大气中通过不同氧化途径生成的二次无机盐和有机物, 不同的物质组成代表了特定的污染来源和颗粒物氧化生成方式[1, 2].京津冀[3, 4]和苏皖鲁豫区域[4~6]秋冬季大气颗粒物污染频发, 由于地理位置和工业分布等差异, 以及本地排放和区域传输叠加气象因素[5~7], 导致造成PM2.5污染的主控因子和特征各异[8~10], 通过对城市PM2.5的中代表性污染物的污染特征及气象要素分析, 可为污染过程解析以及防控政策的制定提供科学依据[6, 11~13].

硫酸盐(SO42-)、硝酸盐(NO3-)和铵盐(NH4+)是北方城市大气PM2.5中重要组成成分[4, 14, 15], 统称为二次无机气溶胶(secondary inorganic aerosol, SNA), SNA主要由煤炭等化石燃料燃烧、工业、机动车和农业畜牧业等排放的SO2、NOx和NH3等经大气氧化、吸附和沉降等物理化学过程二次转化生成, 统称二次无机气溶胶[9, 12, 16, 17], 在河南省濮阳市2018年1月污染时段SNA在PM2.5占比为66%[18], 在洛阳市2019年冬季占比为47.3%[13].冬季气团传输、煤炭和生物质燃烧、交通排放等均是重要的贡献源[2, 5, 17], 低温有利于SNA以颗粒态形式存在, 静稳天气使颗粒物滞留时间变久, 也有利于污染物的二次转化, 使SNA成为冬季PM2.5污染过程中主要贡献因子[19], 硫、氮氧化率以及SO42-和NO3-增长速率则是评价其对灰霾影响的重要指标[20].碳质气溶胶(carbonaceous aerosol, CA)也是局地和区域性灰霾时段的重要污染物, 一般分为无机碳(element carbon, EC)和有机碳(organic carbon, OC)[21].北方城市秋冬季污染多以SNA为主, 源解析结果显示, 来自煤炭燃烧、工业过程、工业锅炉和机动车排放的SO2和NO2是生成颗粒物中SO42-和NO3-的主要前体物[22~24].EC主要源于化石燃料或生物质燃料等含碳物质的不完全燃烧; OC来源复杂, 既有来自燃烧源的直接排放的一次有机碳(primary organic carbon, POC), 又有通过大气的氧化作用形成的二次有机碳(secondary organic carbon, SOC)[10, 25~27], 两者通常占PM2.5总质量的20%~40%. PM2.5中的K+是秸秆焚烧的常用指示物[5, 16]. NO2/SO2通常代表移动源和燃烧源的来源贡献比例, NO3-和SO42-为NO2和SO2经大气化学反应生成的二次无机盐, NO2/SO2或NO3-/SO42-越大, 除表示某类源贡献大[23, 28]之外, 洛阳[13]、北京、天津[29]和部分亚洲东北部城市的研究结果表明[30], 大气细颗粒物中NO3-的占比和生成速率[31, 32]超过了SO42-, 成为细颗粒物中的主要增加组分[13].

周口市位于黄淮平原腹地苏皖鲁豫区域, 东临安徽省亳州和阜阳市, 属于苏皖鲁豫区域管控范围, 整体空气质量在省内排名较前, 但冬季在气象不利和区域污染时, 颗粒物污染可达重度.受观测仪器等条件所限, 尚未有对其冬季细颗粒物污染主要推高因素的综合分析.本文基于周口边界站高时间分辨率的在线观测数据, 并基于2022年1月空气质量参数、细颗粒物中离子和元素等成分特征的统计, 对细颗粒物污染分布特征和真实大气环境下SO42-和NO3-快速增长趋势及其潜在影响因素进行分析, 研究结果可有助于进一步了解苏皖鲁豫和河南南部城市冬季PM2.5生成特征, 以期为当地大气污染管控提供科学依据.

1 材料与方法 1.1 观测时间和观测点位

周口市位于河南省东南部(33°03′~34°20′N, 114°05′~115°39′E), 东部与安徽省亳州市和阜阳市接壤, 属于暖温带半湿润季风性气候, 冬季寒冷少雨雪, 风向以东北偏北为主, 主要受蒙古冷高环流影响.由省监测平台统计结果表明, 2021年年均和冬防期间(2021-10-01~2022-03-31)ρ(PM2.5)均值分别为44 μg·m-3和66 μg·m-3, 分别超出我国《环境空气质量标准》(GB 3095-2012)二级浓度限值35 μg·m-3的1.3倍和1.9倍.其中2022年1月, 轻度~中度17 d, 重度8 d, 是冬防期间颗粒物污染天数最多的月份.本市工业企业以塑编、家具制造、水泥制品和烧结砖瓦制品等中小企业为主, 其中2018年成立的河南安钢周口钢铁有限责任公司位于市区东南约60 km的沈丘县, 一期高炉和电炉等生产系统于2020年4月开始试运行.

本文所研究的空气质量污染因子数据来自位于川汇区政府站, 该国控站点位于市中心区域, 周边主要为高新科技企业园区、居民住宅及城市主干道中州大道和太清路, 丰泉环保电力有限公司距离该点位约2 km.气象数据来自中国气象数据网(http://data.cma.cn)发布的城市小时数据, 包括气温、相对湿度、风速和降水量等.空气质量AQI的分类, 根据河南省省辖市预报系统中采用的《环境空气质量评价计算规范(试行)》(HJ 663-2013)中对于小时和日评价时段和对应的统计方法, 分为优(0~50)、良(51~100)、轻度(101~150)、中度(151~200)、重度(201~300)和严重污染(301~500)等不同污染时段和污染日.本研究时段为2022年1月1~31日, 其中包含了4次PM2.5中度到重度污染过程, 污染日采用日评价方法, AQI分级统计采用小时评价方法.

1.2 观测仪器

周口市川汇区政府站点位空气质量数据均来自在线仪器, 其中水溶性离子如硫酸盐(SO42-)、硝酸盐(NO3-)、氯离子(Cl-)、铵根离子(NH4+)和钾离子(K+)等的测定, 采用离子色谱分析仪(URG 9000, Thermo Fisher Scientific, 美国); 元素如Al、Si、Ca、Fe、Ti和K等采用无机元素分析仪(AMMS-100, 聚光科技); EC和OC采用在线EC/OC分析仪(RT-4, Sunset, 美国); 常规空气6参数如PM2.5、PM10、NO2、SO2、CO和O3分析仪采用安徽蓝盾光电子股份有限公司LGH系列仪器.所用仪器检测方法和质控措施, 遵从环境空气颗粒物(PM10、PM2.5)连续自动监测系统技术要求及检测方法(HJ 653-2013)和环境空气气态污染物(SO2、NO2、O3、CO)连续自动监测系统技术要求及检测方法(HJ 654-2013).

1.3 分析方法

已有研究中对PM2.5中阴阳离子组成、硫氧化率、氮氧化率和二次生成有机碳的计算, 阴离子电荷当量用AE表示(anion equation, AE), 阳离子电荷当量用CE(cation equation, CE)表示[33], 如式(1)和式(2)所示.

(1)
(2)

通常用硫氧化率(SOR)和氮氧化率(NOR)分别表示大气中气态的SO2和NO2在大气中经氧化反应后生成二次无机离子SO42-和NO3-转换效率[4, 18, 23, 33], 如式(3)和式(4)所示.

(3)
(4)

式中, SO42-、NO3-、SO2和NO2均需换算为mol·m-3单位.

在大气PM2.5中, 通常存在未被监测到的有机物含量, 为计算包含C、H、O、N和S等的有机物, 常进行有机物(organic matter, OM)重构[6], OM采用有机碳OC来计算, 如式(5):

(5)

式中, f为折算系数.

由于大气氧化程度各异, 生成OM的f值范围通常在1.2~2.6, 数值从小到大分别代表了新鲜到逐步老化的气溶胶.根据统计, 我国北方城市冬季平均f值通常取1.59±0.18, 本文中采用[OM]=1.6×[OC]来计算[19].二次有机气溶胶(secondary organic aerosol, SOA)[25, 33, 34]代表在大气氧化过程中二次生成的有机物, 通常用公式SOA=1.4×SOC计算, 二次有机碳(SOC)采用OC/EC最小比值法[式(6)]来估算:

(6)
2 结果与讨论 2.1 1月空气质量及PM2.5组成特征

周口市2022年1月空气质量和主要气象参数(日均值)如表 1所示, 从中可知, 1月ρ(PM2.5)均值为113 μg·m-3, 风速(WS)均值为1.1 m·s-1, 相对湿度(RH)均值为77%, 风速较低, 降水量较少, 整体不利于污染物的扩散和清除.观测值统计结果表明, PM2.5占PM10的质量分数为83.8%(范围16.7%~99.5%), 颗粒物中细颗粒物质量分数较高, NO2/SO2为4.57(1.0~16.7), 除SO2有短时高值外, 多数情况下NO2浓度高于SO2, 表明本地以移动源、锅炉等燃烧排放的NO2对大气的贡献排放大于煤炭燃烧排放的SO2, 但偶发性燃烧源贡献也不容忽视.

表 1 周口市2022年1月空气质量与主要气象参数1) Table 1 Air quality and main meteorological parameters of Zhoukou in January 2022

周口市1月小时AQI范围为19~275, 以100为单位将数据分为3个污染段, 不同污染小时指标按AQI分: 0~100(优良)、101~200(轻度~中度)和201~300(重度), 不同污染等级下细颗粒物和主要组成变化如表 2所示, 其为主要污染因子和离子、元素(小时值)在不同污染段的分布特征.由表 2可见, 随着污染程度的增加, 粗、细颗粒物浓度增加, 细颗粒物在粗颗粒物中所占质量分数在轻~中、重度污染分别是优良时段的1.10倍和1.16倍, 大气氧化性在中度污染时是清洁时段的0.99倍, 重度污染时是清洁时段的1.10倍).SOR在重污染时段比轻~中度和优良时段增加了1.23倍和1.46倍, NOR则分别为1.33倍和1.45倍, NO3-和SO42-比值在重度污染时段是优良时段的1.16倍, 轻~中度污染时段比值与清洁时段相似.在重度污染时段, 大气氧化性升高, NOR和SOR的升高导致了NO3-、SO42-氧化率和增长速率的升高, 来自机动车和锅炉等燃烧排放的氮氧化物不仅增强了大气氧化性, 也为大气氧化反应生成NO3-提供了气态前体物.

表 2 周口市1月不同污染程度下PM2.5组成分布特征1) Table 2 Distribution of atmosphere parameters during different AQI level periods

周口市PM2.5中NO3-/SO42-较高, 比北京清洁和污染时段的1.12和1.17均高出约2~3倍[19], 且在重度污染时段NO3-增加量超过SO42-, NO3-/SO42-在3个不同的污染段均值分别为2.66、2.66和3.07.1月NO3-/SO42-均值为2.75(范围0.82~6.64), 也略高于本省工业城市洛阳市2019年10月至2020年1月的均值2.4[13, 22], 说明机动车和锅炉等排放对周口市细颗粒物的影响大于燃煤, 氮氧化物排放对PM2.5的贡献要高于北京、洛阳等北方城市.SNA浓度在重度污染时段是优良和轻-中度污染时段的2.89倍和1.75倍, SNA中NO3-的快速增长是造成1月PM2.5重污染的主要原因.CA和SOA随污染程度加重也有所增加, 仅CM呈下降趋势.

对细颗粒物中组分结果进行统计分析, 1月周口PM2.5中AE与CE相关性R2为0.985, 斜率1.00, 说明在线观测数据质量可靠, 细颗粒物中SNA来源相对稳定, 无机盐组成接近分子形态.在优良时段, AE与CE相关性R2为0.911, 斜率1.05, 轻~重度污染时段, R2为0.982, 斜率1.01, 为污染时段阴离子如SO42-、NO3-占比升高所致.1月PM2.5中AE与CE均值分别为0.93(0.20~1.84)和0.94(0.25~1.90), CE略高于AE, 阴阳离子分布特征相似.

1月SOR和NOR分别为0.53(0.13~0.87)和0.46(0.17~0.73), SOR约为NOR的1.15倍, 分别高于濮阳2017~2018年秋冬季的2.4倍和2.1倍[18].随着污染等级的增加, SOR和NOR随之增加, NO3-和SO42-同步升高, 该趋势与2018年和2019年河南北部城市SNA分布趋势相似[6, 13]. PM2.5在PM10中的占比在优良和轻中度污染时段分别为76%(18%~98%)和85%(33%~98%), 重度污染时段为89%(80%~97%).同时, SOA随AQI升高而升高, 但在PM2.5中占比下降, 由优良时段的4.28%(0.82%~9.63%)到轻-中度的2.29%(0.09%~4.55%)和重度的1.90%(0.31%~3.03%).可能由于冬季污染时段SOA气相氧化通道作用变弱, 导致了污染程度增加, SOA浓度上升但在PM2.5中的贡献率下降.

2.2 重污染过程分析

2022年1月共经历了4次细颗粒物由中度到重度的污染过程(如图 1图 2中阴影部分所示), 其中1月2、10、12和24日为中度污染, ρ(PM2.5)日均值为116~140 μg·m-3, 1月3~4、11、16~19和26日为重度污染, ρ(PM2.5)日均值为156~196 μg·m-3, 1月4日00:00达1月小时最高值ρ(PM2.5 236 μg·m-3).如图 1图 2所示, 在中~重度污染前期, CO、NO2和SO2浓度有不同程度上升.与优良时段相比, 污染状态下AE和CE均同步升高, 但AE/CE基本稳定, SOR、NOR和SOA呈不同程度升高, 二次有机气溶胶SOA和碳质气溶胶CA对颗粒物浓度升高有一定贡献.在1月的4次污染过程中, 随着污染过程加重, 由于污染过程中气态前体物浓度不同程度地升高, SOR和NOR升高, 二次离子生成速率加快, 导致颗粒物中二次无机盐SNA含量增加.同时, 污染时段相对湿度增加, 风速变小, 气象条件也有利于细颗粒物的二次生成和聚集.高排放和高的二次转化率, 叠加高湿、低风速的静稳天气, 是周口市1月的4次污染过程的主要原因.污染过程后期, 随着前体物浓度降低, 风速增大, 二次离子生成速率呈下降趋势(表现为负值, 负值绝对值越小, 表明逐小时下降速度越快), SNA浓度下降, 颗粒物浓度也迅速下降至污染解除.

阴影部分Ⅰ~Ⅳ为4次PM2.5污染过程 图 1 周口市2022年1月空气质量与主要气象参数时间变化 Fig. 1 Time series of air quality and main meteorological parameters of Zhoukou in January 2022

阴影部分Ⅰ~Ⅳ为4次细颗粒物污染过程; (SOA/PM2.5)×100表示曲线为(SOA/PM2.5)原数值乘相应倍数后数值 图 2 周口市2022年1月细颗粒物中主要组分及组成变化 Fig. 2 Time series of variation in main components in PM2.5 of Zhoukou in January 2022

除SO2的短时高值外, NO2和SO2均低于文献报道中濮阳相应数值, 可能由于位于市区的观测站点周边分布了较多的固定源和移动源, 汽车尾气排放的NOx、碳质气溶胶和钢铁厂或散煤燃烧排放的SO2以及与喷涂等有关的工厂VOCs排放, 在冬季大气环境中发生氧化反应, 造成了较高的SOR和NOR, 因此在较低的气态前体物情况下导致了较高浓度的SO42-和NO3-, 对PM2.5中的SNA生成起到了促进作用.

2.3 重污染过程成因分析

硫/氮氧化率数值与大气氧化性、温度和相对湿度等相关, 直接影响大气中SO2和NO2的二次无机盐生成效率和增长速率.当SOR和NOR大于0.25和0.1时, 表明大气中发生了较强的光氧化反应[4].由表 1图 2可知, 1月周口市的SOR和NOR均在0.25和0.1以上, 表明大气中存在着较多的SO2和NO2氧化反应.如图 2所示, 在1月的Ⅰ~Ⅳ次污染过程中, SOR和NOR呈上升趋势, 多数时段SOR高于NOR.

二次无机盐中SO42-和NO3-增长速率的计算, 是用当前小时SO42-或NO3-浓度值减去前N个小时相应浓度值, 分母即为时间间隔N个小时.本文采用1h时间间隔, 分母为1, 该值可代表细颗粒物中SO42-和NO3-小时变化趋势, 比如在颗粒物污染增长期间, 该值为正, 数值越高代表PM2.5中的SNA增加越快, 在污染消退阶段, 该值为负, 负值绝对值越大代表了PM2.5污染消退越快, 数值大小与前体物排放浓度变化强度、SOR/NOR大小以及风速和相对湿度等有关.如在前体物高排放时段叠加SOR和NOR增加, 增长速率呈正值, 表现在颗粒物中为二次无机盐浓度增加, 如前体物排放降低、SOR和NOR增长趋缓或风速加大或降雨时段, 二次无机盐的增长速率则为负值, 表明此阶段细颗粒物中SO42-和NO3-浓度呈下降趋势.

在第一次污染过程中, 空气质量日AQI值从1月1日轻度到1月2日的中度、3~4日的重度, 再降低到5日的良, ρ(SNA)由1月1日轻度污染时的均值38.2 μg·m-3增加到了3~4日重污染日的96.2 μg·m-3, 升高了近2.5倍, SO42-浓度从1月1日的4.86 μg·m-3, 逐渐增加到2~4日的8.14、10.3和16.9 μg·m-3, 与此同时, NO3-浓度也从24.6 μg·m-3增加到39.5、54.7和70.2 μg·m-3.SOA则由1月1日的4.89 μg·m-3降到了3~4日的3.44 μg·m-3, 下降了约70%, 此次污染过程中二次组分的增加以SNA为主, SOA对本市冬季细颗粒物污染贡献不显著.

在此过程中, SOR日均值在1月1~5日分别为: 0.21、0.20、0.30、0.47和0.55, NOR分别为: 0.26、0.36、0.46、0.54和0.46, 可见随着污染程度的升高, SOR在逐步升高, 在污染消退时仍处于升高趋势, NOR日均值则从升高到5日有所下降, 说明在污染抬升过程中, SOR和NOR同时促进了SNA的快速生成, NOR略高于SOR, 但在污染消退阶段, NOR下降, SOR仍持续上升, 重度污染时段约是优良时段的1.4倍.

同时, 1月1~5日SO42-日均增长速率分别为: -0.07、0.11、0.40、-0.53和-0.16 μg·(m3·h)-1, NO3-日均增长速率分别为: -0.25、1.03、1.15、-1.24和-2.14 μg·(m3·h)-1, 在污染增长初期阶段, PM2.5中的SO42-和NO3-有短时的增长和消减, 日均增长速率为负, 随着污染加重, 1月2~3日两种离子均呈现出增长趋势, NO3-增长速率是SO42-的近10和3倍, 虽然1月4日仍是重污染, 但由图 1可见, 由于风速增大, SO42-和NO3-日均增长速率均为负值, 表现为下降趋势, 到1月5日AQI转良时下降速度更快.1月2~4日SO42-和NO3-增长速率分别同比1日同时段升高了0.25 μg·(m3·h)-1和1.0 μg·(m3·h)-1.结合图 1可见, 在1月3日, 温度下降, RH%增加, WS变小, 低温和高湿促进了SO42-和NO3-的转化和生成, 低风速加剧了SNA本地聚集.

Ⅱ~Ⅳ次污染过程呈现类似情形, 快速升高的SOR和NOR使得SO42-和NO3-更快地生成, 在PM2.5中的浓度快速增加, 表现为1月污染时段较高的SO42-和NO3-增长速率和被迅速推高的PM2.5浓度, 说明在冬季污染过程中, 污染时段更高的SOR和NOR是导致细颗粒物浓度增加的重要因素.

在1月Ⅰ~Ⅳ(图 2阴影部分)4次细颗粒物污染过程中, SNA(NH4+、NO3-、SO42-)和有机物OM是细颗粒物中含量较高且增加较显著的组分, OC/EC在2.54~6.89范围之间, 污染日OC/EC增加不明显, OC和SOA在污染时段有所升高, 但SNA是污染期间升高最显著组分, 且以NO3-增加为主.

在细颗粒物污染消退阶段, 虽然SOR和NOR仍呈增加趋势, 但由于气态前体物浓度下降, 风速增加, 两种无机离子浓度在同步下降, 表现为SO42-和NO3-增长速率呈负值.1月NO3-下降的极值出现在1月19日06:00, SO42-和NO3-的增长速率分别为-5.59μg·(m3·h)-1和-23.1 μg·(m3·h)-1, 两者分别由05:00的15.1 μg·m-3和55.6 μg·m-3下降到了06:00的9.47 μg·m-3和32.5 μg·m-3, 此时风速由0.3 m·s-1增为1.4 m·s-1, ρ(PM2.5)也由199 μg·m-3降到155 μg·m-3.

在真实的大气环境下, 1月NO3-的平均增长速率[0.05 μg·(m3·h)-1]略高于SO42-平均增长速率[0.03 μg·(m3·h)-1], NO3-浓度在细颗粒物污染过程阶段较SO42-上升更快, 在污染消退阶段下降速率也大于SO42-.由于尚未证实气态前体物与SNA线性对应关系[35], 因此气象不利情况下, 需对NO3-的气态前体物氮氧化物进行提前管控, 将有助于污染峰值削减.

基于在线小时时间分辨率的观测数据, 对不同污染等级时段分为优良、轻度及以上和中~重度污染这3个不同污染范围, 分析不同污染等级范围中主要污染因子日变化规律, 结果如图 3所示.由图 3可见, 在优良天PM2.5呈现出早晚交通高峰升高的趋势, 在轻度及以上污染时夜间PM2.5略有增加, 在中~重度污染时段下午开始有明显增加, 与O3的高值时段吻合, 增加趋势持续到夜间和凌晨.同时, 污染时段相对湿度高出优良时段约10%, 高相对湿度也促进了SNA的二次生成和向颗粒相转化.

绿色系点线表示优良时段,橙色系点线表示轻度及以上污染,紫色系点线表示中~重度污染时段 图 3 不同污染等级时段下颗粒物及其组成元素和离子分布特征 Fig. 3 Distribution characteristics of PM2.5 and its components during three AQI level periods

气态前体物SO2在中~重度污染时段浓度显著高于清洁时段, 也高于轻度及以上污染时段, SO2高值时段仅在部分中~重度污染时段出现, 显著升高时段出现在夜间.SOR也随着污染加重而升高, 中~重度污染的高值时段接近SO2高值时段, 在轻度及以上时段SOR显著高于优良时段, SO42-则呈现明显分层, 中~重度污染范围下SO42-浓度高于轻度以上和优良时段, 日变化呈现出下午、夜间和凌晨高值, 与SOR和SO2高值时段吻合.不同污染程度下NO2的高峰在夜间有明显抬升, NOR峰值出现在下午和凌晨时段, 与O3高值时段重合.SOR在优良、轻度、中度和重度时段均值分别为0.41、0.54、0.54和0.62, NOR在这些污染等级下的均值分别为0.36、0.41、0.49和0.54, 重污染时段SOR与NOR高于优良时段约1.5倍, 这与北京冬季(SOR霾天和清洁天分别为0.20和0.15, NOR分别为0.23和0.10)[33]、郑州年均(SOR为0.35, NOR为0.19)[4]、郑州冬季(SOR霾天和清洁天分别为0.28和0.15, NOR分别为0.29和0.09)[19]、濮阳(重污染时段SOR为0.32, NOR为0.36)[18]和川南地区(NOR在清洁天为0.11~0.17, 污染天为0.24~0.32; SOR在清洁天为0.22~0.51, 污染天为0.30~0.53)[11]等研究的结果趋势一致, SOR和NOR与濮阳接近, 高于其他地区, 说明周口市冬季大气氧化性和SNA转化率较高.同时, 一次污染物普遍存在的夜间抬升现象, 与夜间随着地表辐射冷却出现逆温层结、大气边界层高度下降有关.

随着PM2.5污染程度增加, 在高RH、AE/CE趋近为1、pH接近中性条件下, 轻度及以上污染时段O3浓度略有下降, 但NO2整体升高, 因此总氧化性Ox上升, SOR和NOR也随之升高, 促进了SO2和NO2向SO42-和NO3-的转化.清洁时段NO3-峰值出现在10:00、12:00和18:00左右, 与同污染等级下SO42-峰值时段接近.但在污染时段, 高值时段出现在夜间和凌晨.1月SOR整体略高于NOR, 但NO3-对PM2.5的贡献显著超过SO42-, 污染时段有明显抬升, 证实了文献中所论证的夜间富氨、高RH%条件下, 大气的非均相反应促进了NO3-的显著增长[36].夜间和上午时段RH高于下午, 污染天RH大于优良天, 高湿度有利于SNA转化, 高氧化性和高湿状态均有利于SNA的生成.Ge团队研究了华北的冬季雾-霾事件[14], 结合WRF-Chem模型得出, SO42-以气相氧化为主, 在液态水中发生的是H2O2氧化反应.Xiao等[29]采用贝叶斯混合模型定量分析天津大气PM2.5中NO3-的形成途径和来源结果表明, 白天高浓度NOx下HONO光解生成NO, NO2与·OH反应生成HNO3, HNO3与大气中的NH3结合, 生成二次无机盐NH4NO3, 夜间NO2表面水解生成HONO后再与NO2氧化, 和高HNO3生成NO3-后与NH3结合生成NO3-[37], 因此, 高浓度NO和NO2均有利于NO3-的生成.

以往研究表明[25], NO3-主要通过光化学氧化和气粒分配方式进入颗粒相, SO42-在冬季则主要通过液相氧化方式进入颗粒相, 相对湿度对二次无机颗粒物生成和气-固分配有较大影响[38].不同相对湿度下颗粒物浓度、气温、风速和离子浓度等昼夜分布如图 4所示, 从中可知, 观测期间低湿度时细颗粒物浓度白天大于夜间, 当湿度增加到70%以上时, 夜间浓度超过白天, Ox(O3+NO2)的浓度随RH增大而呈近线性下降趋势, RH大于70%后夜间Ox浓度超过白天.观测期间昼夜风速近似, 气温白天高于夜间, 高湿伴随低温出现, 叠加高湿时夜间高Ox, 共同推高了夜间高湿条件下SO42-和NO3-的升高, SO42-浓度随RH升高而升高, 夜间SO42-升高更显著.NO3-浓度随RH升高无明显变化, 相对湿度大于90%后NO3-浓度呈下降趋势, 低RH时段NO3-与SO42-的比值较高RH时大, 在RH<50%时NO3-/SO42-为4.02, RH>50%时该比值为2.63, 说明周口1月低湿度时NO3-生成效率大于SO42-, 与北京2014~2015年冬天低湿度时(RH<50%)比值2.1下降到高湿时(RH>50%)的0.40趋势一致, 数值则高北京数值近2倍[31, 33].

K+×5和SOA×10分别表示曲线为K+和SOA原数值乘相应倍数后数值 图 4 不同相对湿度下颗粒物浓度、气象要素和离子等参数变化 Fig. 4 Distribution trends of PM, meteorological parameters, and ions under different RH

结合SOR和NOR与RH相关性图可知, 随相对湿度增加, SOR增加幅度大于NOR, 在湿度由45%增加到68%时多数NOR大于SOR, 但在湿度69%之后, NOR增加趋缓, SOR则较快增加, 证明了高湿有利于SOR转化, 真实大气观测结果证明了高湿条件下SO42-生成速率加快.NO3-的急剧增加主要是由于硝酸与氨的均相反应, SO42-则可能是NO2、SO2和NH3的非均相反应[16].阳离子NH4+随相对湿度增加有略微增加, 高湿有利于NH4NO3和NH4(SO4)2生成.

代表生物质燃烧的K+对RH影响不显著, 昼夜均有高值出现, 二次无机盐SNA与RH变化趋势与铵盐近似, 可能与NH4+主导的SNA液相化学生成有关.夜间SOA随RH增加有显著增加, 相对湿度对夜间SOA的形成有促进作用, 但对白天SOA影响较小, 可能由于昼夜SOA生成路径不同, 夜间以NO3液相化学反应为主, 白天则以O3的气相氧化为主[1], 符合此前研究结果随湿度增加大气氧化性呈下降的趋势.冬季SOA生成气相氧化通道作用较弱, SNA随湿度增加的同时, 白天SOA随湿度变化并不明显, 该结果与2017~2018年河南北部城市濮阳研究结果相似[18].在冬季高湿条件下, 白天和夜间的气相氧化和颗粒相表面的非均相氧化分别促进昼夜SOA和SNA的生成[39].

3 结论

(1) 周口市2022年1月PM2.5中占比最高前3的组分为水溶性离子(均值为61.3%, 范围为24.6%~93.2%, 下同)、CA(24.3%, 11.0%~52.0%)和CM(9.72%, 1.16%~62.3%).在1月的Ⅰ~Ⅳ共4次污染过程中, SNA分别占PM2.5的50.6%、65.3%、52.2%和58.0%, 是PM2.5中占比最高的一类组成.其次为CA, 占比分别为24.2%、19.3%、25.1%和20.9%, CM占比分别为9.97%、7.10%、8.24%和6.81%.

(2) 随着AQI增加, PM2.5在PM10中的占比升高, 污染过程中NO3-增长速率大于SO42-增长速率, NO3-/CO和NOR升高幅度大于SO42-/CO和SOR, PM2.5中SNA的增加以NO3-为主.SOR和NOR在不同污染等级下昼夜变化与SO42-和NO3-的浓度变化吻合, 虽然在第Ⅲ和Ⅳ次污染过程中SOR大于NOR, 但由于气态前体物浓度的差异, NO3-增长速度仍大于SO42-, PM2.5中仍以NO3-为主.因此, 控制NO和NO2排放有助于气象不利条件下降低PM2.5污染峰值.

(3) 随相对湿度的增加PM2.5浓度升高但Ox下降, SOR和NOR同步增加, SOR和NOR污染时段大于清洁时段且随RH增加而升高, 随RH增加昼夜折线式上升, SOR昼夜随湿度增加较NOR显著.SOA夜间浓度和增加幅度均大于昼间, SNA在RH<60%时昼间浓度大于夜间, RH>60%时夜间浓度升高幅度增大, 部分时段大于夜间.

参考文献
[1] Jia C H, Tong S R, Zhang X R, et al. Atmospheric oxidizing capacity in autumn Beijing: analysis of the O3 and PM2.5 episodes based on observation-based model[J]. Journal of Environmental Sciences, 2023, 124: 557-569. DOI:10.1016/j.jes.2021.11.020
[2] Lv L L, Wei P, Hu J N, et al. Source apportionment and regional transport of PM2.5 during haze episodes in Beijing combined with multiple models[J]. Atmospheric Research, 2022, 266. DOI:10.1016/j.atmosres.2021.105957
[3] 徐双喜, 张众志, 杜晓惠, 等. 京津冀及周边民用散煤燃烧控制对北京市PM2.5的影响[J]. 环境科学研究, 2021, 34(12): 2876-2886.
Xu S X, Zhang Z Z, Du X H, et al. Impact of residential coal combustion control in Beijing-Tianjin-Hebei and surrounding region on PM2.5 in Beijing[J]. Research of Environmental Sciences, 2021, 34(12): 2876-2886.
[4] 张俊美, 陈仕霖, 王乾恒, 等. 郑州市大气PM2.5中水溶性离子的污染特征及来源解析[J]. 环境科学, 2023, 44(2): 602-610.
Zhang J M, Chen S L, Wang Q H, et al. Pollution characteristics and sources of water-soluble Ions in PM2.5 in Zhengzhou City[J]. Environmental Science, 2023, 44(2): 602-610.
[5] Li Y L, Geng Y P, Hu X M, et al. Seasonal differences in sources and formation processes of PM2.5 nitrate in an urban environment of north China[J]. Journal of Environmental Sciences, 2022, 120: 94-104. DOI:10.1016/j.jes.2021.08.020
[6] 王申博, 王玲玲, 范相阁, 等. 河南省北部区域霾污染过程中城市和农村点位PM2.5组分差异[J]. 环境科学, 2023, 44(1): 11-19.
Wang S B, Wang L L, Fan X G, et al. Differences in PM2.5 components between urban and rural sites during heavy haze event in northern Henan Province[J]. Environmental Science, 2023, 44(1): 11-19.
[7] 熊一帆, 丁秋冀, 舒卓智, 等. 基于数值模拟与资料同化探究长三角地区冬季PM2.5污染过程的气象影响[J]. 环境科学学报, 2022, 42(4): 293-303.
Xiong Y F, Ding Q J, Shu Z Z, et al. The influence of meteorological parameters on particulate matter in the Yangtze River Delta Region based on numerical simulation and data assimilation[J]. Acta Scientiae Circumstantiae, 2022, 42(4): 293-303.
[8] Fu Z H, Cheng L B, Ye X N, et al. Characteristics of aerosol chemistry and acidity in Shanghai after PM2.5 satisfied national guideline: Insight into future emission control[J]. Science of the Total Environment, 2022, 827. DOI:10.1016/j.scitotenv.2022.154319
[9] Hong J, Tang M, Wang Q Q, et al. Measurement report: wintertime new particle formation in the rural area of north China plain: influencing factors and possible formation mechanism[J]. Atmospheric Chemistry and Physics, 2022. DOI:10.5194/acp-2022-784
[10] Li M Y, Yu S C, Chen X, et al. Impacts of condensable particulate matter on atmospheric organic aerosols and fine particulate matter(PM2.5) in China[J]. Atmospheric Chemistry and Physics, 2022, 22(17): 11845-11866. DOI:10.5194/acp-22-11845-2022
[11] 吴安南, 黄小娟, 何仁江, 等. "大气十条"实施结束川南城市群秋季霾污染过程中水溶性离子特征[J]. 环境科学, 2022, 43(3): 1170-1179.
Wu A N, Huang X J, He R J, et al. Characteristics of water-soluble Ions in an autumn haze process in the southern Sichuan urban agglomeration after the implementation of China's Air Pollution Prevention and Control Action Plan[J]. Environmental Science, 2022, 43(3): 1170-1179.
[12] 路娜, 李治国, 周静博, 等. 2015年石家庄市采暖期一次重污染过程细颗粒物在线来源解析[J]. 环境科学, 2017, 38(3): 884-893.
Lu N, Li Z G, Zhou J B, et al. Online source analysis of particulate matter(PM2.5) in a heavy pollution process of Shijiazhuang city during heating period in 2015[J]. Environmental Science, 2017, 38(3): 884-893. DOI:10.13227/j.hjkx.201605058
[13] Xu M, Liu Z R, Hu B, et al. Chemical characterization and source identification of PM2.5 in Luoyang after the clean air actions[J]. Journal of Environmental Sciences, 2022, 115: 265-276. DOI:10.1016/j.jes.2021.06.021
[14] Wang T T, Liu M X, Liu M Y, et al. Sulfate formation apportionment during winter haze events in north China[J]. Environmental Science & Technology, 2022, 56(12): 7771-7778.
[15] Wang J Z, Yang Y Q, Jiang X F, et al. Observational study of the PM2.5 and O3 superposition-composite pollution event during spring 2020 in Beijing associated with the water vapor conveyor belt in the northern hemisphere[J]. Atmospheric Environment, 2022, 272. DOI:10.1016/j.atmosenv.2022.118966
[16] Bao Z E, Chen L H, Li K W, et al. Meteorological and chemical impacts on PM2.5 during a haze episode in a heavily polluted basin city of eastern China[J]. Environmental Pollution, 2019, 250: 520-529. DOI:10.1016/j.envpol.2019.04.045
[17] Lin Y C, Yu M Y, Xie F, et al. Anthropogenic emission sources of sulfate aerosols in Hangzhou, east China: insights from isotope techniques with consideration of fractionation effects between gas-to-particle transformations[J]. Environmental Science & Technology, 2022, 56(7): 3905-3914.
[18] 陈楚, 王体健, 李源昊, 等. 濮阳市秋冬季大气细颗粒物污染特征及来源解析[J]. 环境科学, 2019, 40(8): 3421-3430.
Chen C, Wang T J, Li Y H, et al. Pollution characteristics and source apportionment of fine particulate matter in autumn and winter in Puyang, China[J]. Environmental Science, 2019, 40(8): 3421-3430. DOI:10.13227/j.hjkx.201901119
[19] Dong Z S, Su F C, Zhang Z Y, et al. Observation of chemical components of PM2.5 and secondary inorganic aerosol formation during haze and sandy haze days in Zhengzhou, China[J]. Journal of Environmental Sciences, 2020, 88: 316-325. DOI:10.1016/j.jes.2019.09.016
[20] Wang D X, Hu J L, Xu Y, et al. Source contributions to primary and secondary inorganic particulate matter during a severe wintertime PM2.5 pollution episode in Xi'an, China[J]. Atmospheric Environment, 2014, 97: 182-194. DOI:10.1016/j.atmosenv.2014.08.020
[21] Li Y C, Yu J Z, Ho S S H, et al. Chemical characteristics of PM2.5 and organic aerosol source analysis during cold front episodes in Hong Kong, China[J]. Atmospheric Research, 2012, 118: 41-51. DOI:10.1016/j.atmosres.2012.05.026
[22] 孙佳傧, 董喆, 李利萍, 等. 洛阳市大气细颗粒物化学组分特征及溯源分析[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. DOI:10.13227/j.hjkx.202104256
[23] 陈巧, 谷超, 徐涛, 等. 伊犁河谷夏季PM2.5和PM10中水溶性无机离子浓度特征和形成机制[J]. 环境科学, 2022, 43(11): 5009-5017.
Chen Q, Gu C, Xu T, et al. Characterization and formation mechanism of water-soluble inorganic ions in PM2.5 and PM10 in summer in the urban agglomeration of the Ili River valley[J]. Environmental Science, 2022, 43(11): 5009-5017.
[24] 赵倩彪, 胡鸣, 伏晴艳. 2016-2020年上海市大气细颗粒物化学组成特征和来源解析研究[J]. 中国环境科学, 2022, 42(11): 5036-5046.
Zhao Q B, Hu M, Fu Q Y. Chemical characterization and source apportionment of fine particulate matter in Shanghai during 2016-2020[J]. China Environmental Science, 2022, 42(11): 5036-5046. DOI:10.3969/j.issn.1000-6923.2022.11.009
[25] 孙玺, 周声圳, 高敏, 等. 珠海2019年秋季一次大气污染过程的特征、成因及来源分析[J]. 环境科学学报, 2022, 42(8): 64-75.
Sun X, Zhou S Z, Gao M, et al. Analysis of the characteristics, formation mechanisms and sources of an atmospheric pollution in Zhuhai during the autumn of 2019[J]. Acta Scientiae Circumstantiae, 2022, 42(8): 64-75.
[26] Tan J H, Duan J C, He K B, et al. Chemical characteristics of PM2.5 during a typical haze episode in Guangzhou[J]. Journal of Environmental Sciences, 2009, 21(6): 774-781. DOI:10.1016/S1001-0742(08)62340-2
[27] 尉鹏, 李娟, 赵森, 等. 西安市"阶梯型"持续污染过程演变及PM2.5区域输送分析[J]. 环境科学学报, 2022, 42(2): 298-307.
Wei P, Li J, Zhao S, et al. Analysis of the evolution of the "stepped" continuous pollution process and the regional transport of PM2.5 in Xi'an[J]. Acta Scientiae Circumstantiae, 2022, 42(2): 298-307.
[28] Du H H, Kong L D, Cheng T T, et al. Insights into summertime haze pollution events over Shanghai based on online water-soluble ionic composition of aerosols[J]. Atmospheric Environment, 2011, 45(29): 5131-5137. DOI:10.1016/j.atmosenv.2011.06.027
[29] Xiao H, Ding S Y, Ji C W, et al. Strict control of biomass burning inhibited particulate matter nitrate pollution over Tianjin: perspective from dual isotopes of nitrate[J]. Atmospheric Environment, 2023, 293. DOI:10.1016/j.atmosenv.2022.119460
[30] Ha Y, Kim J, Lee S, et al. Spatiotemporal differences on the real-time physicochemical characteristics of PM2.5 particles in four Northeast Asian countries during Winter and Summer 2020-2021[J]. Atmospheric Research, 2023, 283. DOI:10.1016/j.atmosres.2022.106581
[31] Huang R J, He Y, Duan J, et al. Contrasting sources and processes of particulate species in haze days with low and high relative humidity in wintertime Beijing[J]. Atmospheric Chemistry and Physics, 2020, 20(14): 9101-9114. DOI:10.5194/acp-20-9101-2020
[32] Qiao B Q, Chen Y, Tian M, et al. Characterization of water soluble inorganic ions and their evolution processes during PM2.5 pollution episodes in a small city in southwest China[J]. Science of the Total Environment, 2019, 650(Pt 2): 2605-2613.
[33] 李杏茹, 白羽, 陈曦, 等. 北京冬季重污染过程大气细颗粒物化学组成特征及来源分析[J]. 环境化学, 2018, 37(11): 2397-2409.
Li X R, Bai Y, Chen X, et al. Chemical composition and source apportionment of PM2.5 during winter in Beijing[J]. Environmental Chemistry, 2018, 37(11): 2397-2409. DOI:10.7524/j.issn.0254-6108.2018011401
[34] 钱骏, 冯小琼, 陈军辉, 等. 四川盆地典型城市PM2.5污染过程组分特征和来源解析[J]. 环境科学学报, 2021, 41(11): 4366-4376.
Qian J, Feng X Q, Chen J H, et al. Composition characteristics and source apportionment of PM2.5 pollution process in typical cities in the Sichuan Basin[J]. Acta Scientiae Circumstantiae, 2021, 41(11): 4366-4376.
[35] Li J W, Zhang Z S, Wu Y F, et al. Effects of chemical compositions in fine particles and their identified sources on hygroscopic growth factor during dry season in urban Guangzhou of South China[J]. Science of the Total Environment, 2021, 801. DOI:10.1016/j.scitotenv.2021.149749
[36] Li Q H, Wu B G, Liu J L, et al. Characteristics of the atmospheric boundary layer and its relation with PM2.5 during haze episodes in winter in the North China Plain[J]. Atmospheric Environment, 2020, 223. DOI:10.1016/j.atmosenv.2020.117265
[37] Fan M Y, Zhang Y L, Lin Y C, et al. Important role of NO3 radical to nitrate formation aloft in urban Beijing: Insights from triple oxygen isotopes measured at the tower[J]. Environmental Science & Technology, 2022, 56(11): 6870-6879.
[38] Prasad P, Basha G, Ratnam M V. Impact of Relative Humidity on the vertical distribution of aerosols over India[J]. Atmospheric Research, 2023, 281. DOI:10.1016/j.atmosres.2022.106468
[39] Huang X J, Zhang J K, Luo B, et al. Water-soluble ions in PM2.5 during spring haze and dust periods in Chengdu, China: Variations, nitrate formation and potential source areas[J]. Environmental Pollution, 2018, 243: 1740-1749.