环境科学  2022, Vol. 43 Issue (1): 74-84    PDF    
运城市PM2.5时空分布特征和潜在源区季节分析
王妘涛1 , 张强1 , 温肖宇1, 窦乃超1, 赵文婷1, 罗淑贞1, 陈志1, 瞿程凯2     
1. 西北大学城市与环境学院, 西安 710127;
2. 中国地质大学(武汉)生物地质与环境地质国家重点实验室, 武汉 430074
摘要: 基于2019年3月~2020年2月环境空气质量监测数据,分析了运城市PM2.5污染的时空分布特征,并利用HYSPLIT后向轨迹模型和聚类分析等方法探讨不同季节运城市PM2.5污染的输送路径和潜在源区.结果表明,运城市ρ(PM2.5)冬季最高(111.24 μg·m-3),夏季最低(30.02 μg·m-3),PM2.5/PM10秋冬季均大于0.6,表明运城市秋冬两季颗粒物污染以细颗粒物为主;空间上ρ(PM2.5)年均值呈现北部和中部高、东部和西部低的分布特征,高值区PM2.5与SO2、NO2和CO呈显著强相关,表明本地排放对高值区ρ(PM2.5)影响较大,春季和冬季最高值分别位于河津市(58.50 μg·m-3)和稷山县(142.33 μg·m-3),夏季最高值位于南部的平陆县(36.92 μg·m-3),秋季则为中部的盐湖区(62.94 μg·m-3).输送分析表明,春季和冬季主要受西北向输送气团的影响,其占比可达58.69%和55.77%,夏季主要为偏东和偏南向气团,秋季来自河南西南部的短距离输送气团污染占比最大,为32.89%;潜在源区范围春夏最小,秋季开始变大,冬季范围最广且贡献最大,最主要源区春夏位于河南西部和中部,秋冬集中在陕西、河南和山西这3省交界区域以及陕西西部.
关键词: PM2.5      时空分布      后向轨迹      输送路径      潜在源区     
Spatiotemporal Distribution and Seasonal Characteristics of Regional Transport of PM2.5 in Yuncheng City
WANG Yun-tao1 , ZHANG Qiang1 , WEN Xiao-yu1 , DOU Nai-chao1 , ZHAO Wen-ting1 , LUO Shu-zhen1 , CHEN Zhi1 , QU Cheng-kai2     
1. College of Urban and Environmental Science, Northwestern University, Xi'an 710127, China;
2. State Key Laboratory of Biogeology and Environmental Geology, China University of Geosciences, Wuhan 430074, China
Abstract: The spatiotemporal distribution characteristics of PM2.5 pollution in Yuncheng City were analyzed based on the ambient air quality monitoring data from March 2019 to February 2020, and the hybrid single-particle Lagrangian integrated trajectory (HYSPLIT) and cluster analysis were used to discuss the contamination transportation pathway and the distribution of potential pollution sources affecting PM2.5 in Yuncheng City in different seasons. The results showed that the PM2.5 concentration in Yuncheng was the highest in winter (111.24 μg·m-3) and the lowest in summer (30.02 μg·m-3). PM2.5/PM10 was greater than 0.6 in autumn and winter, indicating that the particulate matter pollution in Yuncheng City in autumn and winter was mainly fine particulate matter; spatially, the annual mean value of PM2.5 presented distribution characteristics that were high in the northern and central regions and low in the eastern and western regions. The high-value areas of PM2.5 were significantly and strongly correlated with SO2, NO2, and CO, indicating that local emissions contributed to the high values. The highest values in spring and winter were located in Hejin City (58.50 μg·m-3) and Jishan County (142.33 μg·m-3), and the highest values in summer and autumn were located in Pinglu County (36.92 μg·m-3) and the Salt Lake area (62.94 μg·m-3), respectively. Transportation analysis showed that spring and winter were mainly affected by air masses transporting northwestward, accounting for 58.69% and 55.77%. In summer, air masses were mainly eastward and southward, and in autumn, pollution was from short-distance transport air masses from southwestern Henan. The largest ratio was 32.89%; the potential source area was the smallest in spring and summer; it then began to grow in autumn, and it was the widest in the winter and had the highest contribution. The main source area in spring and summer was located in western and central Henan, and in autumn and winter, source areas were concentrated at the junction of Shaanxi, Shanxi, and Henan and west of Shaanxi.
Key words: PM2.5      spatiotemporal distribution      backward trajectory      transport pathway      potential source contribution     

随着经济和城市化进程的飞速发展, 我国面临着严峻的大气污染问题, 尤其是细颗粒物(PM2.5)污染[ 1~ 3].有研究表明, PM2.5的时空分布特征存在显著的区域性差异, 其主要受本地排放强度[ 4]、地形起伏度[ 5, 6]、植被覆盖度[ 7]和气象条件[ 8, 9]等因素影响.不同区域PM2.5的首要污染来源也有所差异, 如北京为机动车排放源[ 10], 天津[ 11]和郑州[ 12]为二次污染源和机动车排放源、南京[ 13]为工业源和二次污染源、太原[ 14]为二次污染源和燃煤源等.因此, 研究PM2.5本地化的时空特征及其来源对于科学治理当地颗粒物污染起着重要作用.此外, PM2.5又因其具有空间溢出效应, 存在跨区域传输[ 15~ 17], 邵玄逸等[ 18]的研究发现, 北京市和唐山市PM2.5外来传输贡献可达48.74%和30.67%, 主要来自临近局地, 西北通道和西南通道的输送; 段时光等[ 19]利用后向轨迹模型分析了郑州市秋冬季PM2.5的传输路径和潜在源区, 发现京津冀传输通道城市群是郑州市PM2.5污染的首要贡献源; Kong等[ 20]的研究利用PSCF和CWT方法分析发现四川盆地的南部、东南部和东部是成都市PM2.5污染的主要源区.分析区域大气污染输送路径和潜在源区, 研究城市与周边区域的传输影响, 对于区域间大气污染联防联控工作有着重要意义.

本文以汾渭平原典型城市——运城市为研究对象, 汾渭平原产业结构偏重, 能源结构偏煤, 交通运输结构以公路为主, 加上不利于扩散的地形和气象条件, 导致污染形势严峻.目前, 针对汾渭平原PM2.5污染已经开展了较多研究[ 21~ 23].李雁宇等[ 24]的研究指出, 汾渭平原东部PM2.5浓度白天高晚上低, 周末比周内高0.2%, 空间上渭南市污染最为严重, 其次为运城市、洛阳市, 三门峡市污染最轻; 黄小刚等[ 25]研究了汾渭平原PM2.5污染的地形效应, 发现PM2.5污染与地形起伏度存在明显的负空间自相关性, 洛阳盆地和西安至临汾间的冲击平原污染最为严重; 张海东[ 26]的研究发现, 西安冬季PM2.5污染的主要潜在源区位于关中地区、山西运城以及河南三门峡、洛阳地区.运城市位于汾渭平原核心位置, 且处在西安至临汾冲击平原的污染传输带上, 其大气质量好坏对于周边城市有着重要影响, 运城市2019年PM2.5年均浓度为61 μg·m-3, 与2018年持平, 在汾渭平原城市群中排名从倒4降为倒3, 这表明运城市PM2.5污染问题仍然比较突出, 是汾渭平原颗粒物污染最为严重的城市之一[ 27].已有学者对运城市颗粒物的化学组成[ 28]和来源解析[ 29]做了相关研究, 但是针对运城市PM2.5污染的本地化时空分布特征和区域传输影响的研究还较为鲜见.

基于此, 本文结合2019年3月~2020年2月PM2.5浓度监测数据, 分析运城市PM2.5污染的时空分布特征, 同时利用TrajStat软件提供的聚类分析、潜在源贡献因子分析和浓度权重轨迹分析, 研究运城市不同季节PM2.5区域传输特征及其潜在源区, 以期为运城市PM2.5治理和减排政策制定提供科学依据, 同时为汾渭平原城市群大气污染联防联控治理提供参考依据.

1 材料与方法 1.1 研究区概况

运城市位于山西省西南部, 北依吕梁山, 东峙中条山, 西、南与陕西、河南隔黄河相望, 地形以低山、丘陵、平原和黄土地貌为主, 相对海拔差异较大, 气候属暖温带大陆性季风气候, 冬季寒冷干燥, 夏季高温多雨, 年均气温13℃, 年均降雨537 mm, 约50%降雨量集中在7~9月, 受地形和气候因素影响, 降雨量从东南向西北递减[ 30, 31].综合地形地貌、气候因素、第二产业布局和社会经济发展状况, 本研究将运城市划分为北部地区(新绛、稷山、河津、万荣、闻喜)、西部地区(临猗、永济)、东部地区(垣曲、绛县)、南部地区(夏县、芮城、平陆)和中心城区(盐湖区)五大片区( 图 1), 本文在此分区基础上分析PM2.5空间分布特征.

图 1 运城市五大片区分区示意 Fig. 1 Divisional schematic diagram of five regions in Yuncheng City

1.2 数据来源

运城市PM2.5和PM10逐小时浓度数据来源于国家空气质量联网监测管理平台( http://60.222.230.181:4017), 对应时段为2019年3月1日至2020年2月29日.后向轨迹模型(HYSPLIT)所用气象场资料来源于NCEP(美国国家环境预报中心)提供的GDAS(全球资料同化系统)数据( ftp://arlftp.arlhq.noaa.gov/pub/archives/gdas1), 该数据每日4个时次, 分别为00:00, 06:00, 12:00, 18:00(UTC, 世界时), 空间分辨率为1°×1°.本研究季节划分春季为3~5月, 夏季为6~8月, 秋季为9~11月, 冬季为12月~次年2月.

1.3 HYSPLIT模型

HYSPLIT模型是由美国国家海洋大气管理中心(NOAA)和澳大利亚气象局共同研发的一种用以计算和分析气流运动、沉降以及扩散的专业模型, 可用来分析多种气象要素输入场、多种物理过程以及不同类型污染物的排放源[ 32, 33], 该模型平流和扩散计算采用拉格朗日方法, 能够直观地表现出气流或粒子的运动轨迹[ 34], 被广泛应用于研究大气污染物的传输特征与来源分析.

本研究利用TrajStat软件[ 35]对运城市进行后向轨迹计算、聚类统计、PSCF和CWT分析, 将运城市市中心(35.02°N, 110.99°E)作为轨迹模拟起始点, 起始高度设置为300m, 该高度既能够反映区域间气流的流动特征, 又可减少近地面摩擦力对其的影响[ 36].每日计算4个时次(00:00, 06:00, 12:00, 18:00, UTC), 模拟时长为24h, 依据Angle Distance算法[ 37]和总空间方差(total spatial variance, TSV)[ 38]的突变点进行聚类, 并结合各季节气流所对应污染物的浓度进行统计分析.

1.4 潜在源贡献因子分析(PSCF)

PSCF是一种基于HYSPLIT后向轨迹终点的受体模型, 可用来识别大气污染物的潜在来源[ 39~ 41].PSCF函数基于空间网格进行概率计算, 本文将研究区域按照0.25°×0.25°划分网格, 按文献[ 42]设定ρ(PM2.5)=75 μg·m-3为秋冬季阈值, 春夏两季由于污染较轻, 阈值设定为当季日均浓度的75分位数[ 43], 当轨迹对应的浓度高于此阈值时, 认为该轨迹为污染轨迹.PSCF值为经过网格ij的污染轨迹端点数mij与该网格内所有轨迹端点数nij的比值, 即:

(1)

由于PSCF是一种条件概率, 当网格内轨迹端点数(nij)比较小时, 可能会使计算结果存在较大不确定性, 为此引入权重系数Wij来增强其统计平衡以降低误差[ 44], 当某一网格内端点数小于研究区内所有网格平均端点数的2倍时, 就使用Wij乘以PSCFij, 即WPSCFij=Wij×PSCFij, 权重系数Wij定义为[ 45]:

(2)
1.5 浓度权重轨迹分析(CWT)

PSCF只能反映网格污染轨迹的占比, 并不能体现污染轨迹对于受点的污染贡献.因此, 本文采用浓度权重轨迹(CWT)分析法[ 46, 47]来计算源区格网的浓度权重, 分析各个轨迹源区对于受点污染的贡献大小, CWT计算公式为:

(3)

式中, CWTij是网格ij对应的PM2.5平均权重浓度, l是轨迹, Cl是轨迹l经过网格ij时所对应PM2.5的浓度, τijl为轨迹l在网格ij内的停留时间(可用网格内轨迹端点数来代替), 为减少误差, CWT同样适用PSCF引用的权重系数Wij, 即:

2 结果与讨论 2.1 颗粒物浓度时间分布特征

2019年3月~2020年2月, 运城市颗粒物污染较为严重, ρ(PM2.5)年均值为58.95 μg·m-3, ρ(PM10)年均值为95.86 μg·m-3, 分别比文献[ 42]中年均值二级浓度限值高68%和36.94%, 分别高出汾渭平原年均值10.55%和4.72%, 高出山西省年均值27.32%和8.06%, PM2.5污染较PM10更为严重.由 图 2 图 3可知, 不同季节颗粒物浓度差异较大. ρ(PM10)表现为冬季高, 夏季低, 春秋居中且相差不大的季节特征, 四季大小为: 冬季>春季>秋季>夏季, 浓度依次为141.85、101.55、86.26和54.18 μg·m-3. ρ(PM2.5)表现为冬季最高, 夏季最低, 秋季大于春季的季节特征, 四季大小为: 冬季>秋季>春季>夏季, 浓度依次为111.24、51.36、43.67和30.02 μg·m-3.从月份来看, ρ(PM2.5)1月最高, 可达135.92 μg·m-3, 7月最低, 为28.93 μg·m-3, ρ(PM10)1月最高, 为161.29 μg·m-3, 8月最低, 为52.07 μg·m-3. ρ(PM2.5)最高值出现在12月22日, 为269.20 μg·m-3, 当日PM2.5/PM10达到0.90, 细颗粒物污染严重, ρ(PM10)最高值出现在5月13日, 为320.50 μg·m-3, 当日PM2.5/PM10为0.28, 主要为粗颗粒物污染.此外, 10月至次年2月之间颗粒物污染持续时间长, 波动剧烈, 这是由于采暖期化石燃料的大量使用增大了污染物的排放量, 加之运城市冬季高湿静稳天气频发, 不利于污染物的扩散, 容易形成区域持续性污染.

上下短横线分别代表最大和最小浓度, 矩形代表位于25%下四分位数和75%上四分位数之间的浓度, 矩形中的×号和横线分别代表均值和中位数 图 2 运城市PM2.5和PM10季节变化特征 Fig. 2 Seasonal changes in PM2.5 and PM10 in Yuncheng City

图(a)和(b)中虚线为环境空气质量标准日均值二级浓度限值, (c)中虚线为分界线 图 3 运城市PM2.5, PM10和PM2.5/PM10逐日变化特征 Fig. 3 Daily variations of PM2.5, PM10 and PM2.5/PM10 in Yuncheng City

PM2.5与PM10的比值可以反映区域内大气细颗粒物与粗颗粒物的占比, 运城市PM2.5/PM10季节均值变化特征为: 冬季(0.76)>秋季(0.61)>夏季(0.56)>春季(0.45), 秋冬季PM2.5/PM10均大于0.6, 说明这两个季节细颗粒物污染更为严重.从 图 3(c)可以明显看出, 运城市PM2.5/PM10小于0.4的主要集中在3~6月之间, 说明运城市春夏两季粗颗粒物污染事件较多, 而且从春到冬, PM2.5/PM10呈逐渐上升趋势, 细颗粒物污染随季节逐渐严重.

2.2 PM2.5空间分布特征

为分析运城市PM2.5浓度空间分布特征, 本研究对2019年运城市区及所属区县共29个站点的PM2.5年均值和季节均值进行了空间插值, 结果见 图 4, 插值采用Arcgis克里金插值法(Kriging)[ 48], 该方法可较好地表征大气污染物的空间分布特征[ 49].同时考虑到各个污染物之间的协同性, 也分别做了SO2、NO2、CO和PM2.5/PM10年均值的空间分布, 结果见 图 5.

图 4 运城市年均PM2.5空间分布特征 Fig. 4 Spatial distribution of average annual mass concentration of PM2.5 in Yuncheng City

图 5 运城市年均CO、SO2、NO2、PM2.5/PM10空间分布特征 Fig. 5 Spatial distribution of average annual mass concentration of CO, SO2, NO2, and PM2.5/PM10in Yuncheng City

图 4可知, 运城市PM2.5年均值呈现北部和中部高, 东部和西部低的污染特征, 北部新稷河一带ρ(PM2.5)最高, 区间可达62.58~72.66 μg·m-3, 西部垣曲县PM2.5污染最轻, ρ(PM2.5)为37.11 μg·m-3.

Pearson相关性分析可知( 表 1), 运城市PM2.5与SO2、NO2、CO和PM10均为显著正相关关系(P < 0.01), 尤其是北部地区和中心城区, 其PM2.5与SO2、NO2和CO相关性均在0.6以上, 属于强相关, PM2.5和PM10相关性高达0.918和0.878, 属极强相关, 这表明该区域PM2.5与SO2、NO2、CO和PM10具有强同源性.SO2和CO来自于工业排放和化石燃料的燃烧, NO2主要来源于汽车尾气排放[ 50], 如 图 5(a)~ 5(c)所示, SO2、NO2和CO高值区均集中在北部地区, 这是由于运城北部属于工业区, 分布着大量冶金、发电、建材和石油化工等高排放企业, 化石燃料使用量大, 而且原辅料、燃料和产品都是通过重型柴油车公路运输, 污染物排放量大, 加之北部地区和中心城区属于盆地区域, 地势较低, 年平均风速较小, 不利于污染物扩散, 导致PM2.5污染严重.而东部和西南部属于山区地形, 产业主要以农、林业为主, 污染物排放少, 加之山脉阻挡作用, 区域传输较少, 所以PM2.5污染较轻.

表 1 运城市PM2.5与SO2、NO2、CO、PM10的Pearson相关系数1) Table 1 Pearson's correlation coefficient between PM2.5, SO2, NO2, CO, and PM10 in Yuncheng City

此外, 由 图 5(d)可知, 北部地区和中心城区PM2.5/PM10相对较小, 这说明该区域粗颗粒物污染更为严重, PM10不仅受人为源排放影响, 还与沙尘传输有关[ 16], 该区域位于西安-临汾冲击平原的中心位置, 易受到来自西安、咸阳、渭南和临汾等地沙尘的传输.

运城市PM2.5空间分布季节变化如 图 6所示.春季运城市ρ(PM2.5)主要分布在34.03~58.50 μg·m-3之间, 高值区集中在北部地区, 与年均值分布相比, 盐湖区污染相对较轻, 这可能是由于采暖期结束, 本地排放量降低, 再加上春季来自西北方位的气团传输距离较长, 风力损耗严重, 对盐湖区传输较少.夏季PM2.5空间分布与春季大体一致, 但污染较轻, 除去平陆县(36.92 μg·m-3)之外, ρ(PM2.5)均低于国家一级限值, 这是由于夏季温度上升, 大气热力条件增强, 加剧大气湍流的交换, 有利于污染物的扩散[ 51].而平陆县ρ(PM2.5)较高可能是受到了三门峡和洛阳等地的传输影响, 由于其地处中条山脉间的峡口处, 从后文气流传输结果来看, 来自东南部的气团都是经此峡口进入运城.秋季ρ(PM2.5)主要分布在32.72~62.94 μg·m-3之间, 较夏季污染加重, 空间上与春夏相比, 盐湖区转为高值区, 这可能是由于秋季来自洛阳和三门峡的污染气团与本地排放相互影响, 加之盐湖区地处运城盆地海拔最低处, 周边山地环绕, 扩散条件差, 导致污染物累积.冬季PM2.5高值区分布与年均相似, 主要为北部地区和中心城区, ρ(PM2.5)分布在57.44~136.00 μg·m-3之间, 污染最为严重, 主要是由于冬季气温较低, 相对湿度较大, PM2.5二次转化和累积严重, 而且风力较弱, 不利于污染物扩散.

图 6 运城市PM2.5空间分布季节特征 Fig. 6 Seasonal characteristics of PM2.5 spatial distribution in Yuncheng City

2.3 输送路径特征

为研究运城市大气PM2.5各季节输送路径特征, 本研究采用HYSPLIT模型计算了运城市(35.02°N, 110.99°E)2019年3月1日至2020年2月29日每6h后向轨迹, 通过总空间方差法(TSV)将后向轨迹气团按季节分别聚类为4、4、5和4类( 图 7), 为进一步量化表征各气流轨迹的污染物浓度特征, 将运城市PM2.5逐日小时浓度与气团轨迹相结合进行统计分析( 表 2).

图 7 运城市不同季节后向轨迹聚类分布 Fig. 7 Backward trajectory clusters in different seasons in Yuncheng City

表 2 四季各类轨迹路径及出现概率和对应的PM2.5浓度统计分析结果 Table 2 Statistical analysis results of all types of trajectory paths, occurrence probability, and PM2.5 concentration in all tracks in the four seasons

分析 图 7 表 2可知, 不同季节聚类轨迹出现概率、途经区域、污染浓度和季节占比均存在明显差异.春季来自宁夏东部, 途经甘肃东北部、陕西北部、铜川和渭南到达运城的轨迹2出现概率和污染轨迹占比均最大, 分别为34.51%和28.57%, 但是其平均ρ(PM2.5)和污染轨迹浓度均最低, 分别为38.48 μg·m-3和67.84 μg·m-3, 来自陕西南部, 途经三门峡到达运城的轨迹3为短距离输送气团, 移动速度慢, 不利于污染物扩散, 导致其污染轨迹占比和浓度均较高, 分别为28.57%和75.41 μg·m-3, 污染轨迹浓度最高的是来自内蒙古中部, 经陕西北部、临汾到达运城的轨迹4, 浓度可达82.72 μg·m-3, 其出现概率和污染轨迹季节占比分别为24.11%和24.18%.夏季主要以短距离输送为主, 来自河南西南部, 经洛阳和三门峡进入运城的轨迹3为主要聚类, 出现概率为36.68%, 污染轨迹占比35.11%, 污染轨迹浓度最高的是起自山西南部, 途经河南西北部、洛阳和三门峡到达运城的轨迹2, 浓度为53.15 μg·m-3, 其污染轨迹季节占比也最高, 为43.62%.与其他季节相比, 夏季到达运城的各气团携带的ρ(PM2.5)低, 气团较为清洁.秋季来自河南的轨迹占比最大, 总共可达51.10%, 其中来自河南西南部, 途经洛阳和三门峡的轨迹1污染轨迹占比和浓度最高, 分别为32.89%和120.22 μg·m-3, 来自铜川和渭南的轨迹2污染轨迹占比及其污染浓度分别为22.37%和108.02 μg·m-3, 污染程度仅次于轨迹1.冬季来自内蒙古南部, 途经陕西北部、渭南到达运城的轨迹1为主要聚类, 出现概率高达33.24%, 污染轨迹占比28.03%, 其污染轨迹浓度较高, 为146.95 μg·m-3, 来自甘肃东北部, 途经铜川和渭南的轨迹4平均ρ(PM2.5)和轨迹浓度最高, 分别为127.44 μg·m-3和153.14 μg·m-3, 其污染轨迹占比为26.36%, 仅次于轨迹1.

综上, 运城市春冬两季源自内蒙古和宁夏等地的西北向气团占比最高, 分别为58.69%和55.77%, 且轨迹污染最为严重, ρ(PM2.5)最高可达82.72 μg·m-3和153.14 μg·m-3. 夏季主要受偏东和偏南气团的影响, 气团较为清洁, ρ(PM2.5)处于26.57~35.58 μg·m-3之间, 相对污染程度较高的气团来自山西南部和河南西南部, 秋季来自河南西南部的气团污染占比最大, 为32.89%, ρ(PM2.5)最高, 可达120.22 μg·m-3.

2.4 污染潜在源区分析

在输送路径分析的基础上, 利用TrajStat软件对运城市PM2.5进行潜在源分析(WPSCF), 计算结果见 图 8.将WPSCF大于0.7的区域作为受点污染最主要的潜在源区.

图 8 运城市PM2.5潜在源贡献四季分布 Fig. 8 Potential source contribution function maps of PM2.5 in four seasons in Yuncheng City

图 8可知, 运城市PM2.5的污染潜在源区季节差异明显.春季最主要的污染源区分布在湖北西北部和河南中部, 其次, 在洛阳市西南部、临汾市南部也存在着零星高值区.夏季WPSCF高值主要分布在运城本地、河南的中部和东部.秋季WPSCF高值区主要分布在陕西的北部、西部以及东南部与三门峡市交界区域、河南的西南部.冬季WPSCF大于0.7的范围明显增大, 高值区主要集中在陕西、河南和山西这3省交界区域、陕西西部、甘肃东北部、宁夏中部和河南东部.

为了进一步探讨源区对运城市PM2.5的污染贡献程度并验证WPSCF分析结果, 对其进行了浓度权重轨迹计算(WCWT), 结果见 图 9, WCWT值越大代表着对受点的污染贡献越大, 本文将WCWT值大于75 μg·m-3的区域定义为主要贡献区.

图 9 运城市PM2.5浓度权重轨迹四季分布 Fig. 9 Concentration weighted trajectory analysis maps of PM2.5 in four seasons in Yuncheng City

图 9可知, WCWT分析结果与WPSCF分析结果相似.春季和夏季WCWT值基本都小于75 μg·m-3, 春季相对高值区域主要分布在陕西东南部、河南南部和山西西南部等区域, 内蒙古南部和宁夏西北部也存在零星高值区, 夏季相对高值区主要分布在河南中部、运城本地和洛阳市东北部.秋季WCWT值大于75 μg·m-3主要集中在运城本地、三门峡市、陕西的西部和东南部, 在湖北的北部和河南的南部也存在着零星的WCWT大于90 μg·m-3的区域.在冬季, WCWT值普遍超过其它季节且大于120 μg·m-3的范围最大, 主要集中在陕西、河南和山西这3省交界区域, 此外, 湖北的西北部、河南的中部和北部、陕西的西北部、甘肃东北部和宁夏中部等零星区域WCWT值也大于100 μg·m-3, 对于运城市冬季PM2.5有着显著影响.

3 结论

(1) 时间分布特征上, 运城市ρ(PM2.5)和ρ(PM10)冬季最高, 分别为111.24μg·m-3和141.85 μg·m-3, 夏季最低, 春秋居中且相差不大. PM2.5/PM10季节均值冬季(0.76)>秋季(0.61)>夏季(0.56)>春季(0.45).秋冬季颗粒物污染主要以细颗粒物为主, PM2.5/PM10日均值小于0.4的主要集中在3~6月, 春夏季粗颗粒物污染事件较多.

(2) 空间分布特征上, 运城市年均PM2.5呈现北部和中部高, 东部和西部低的污染特征, 新稷河一带污染最重, ρ(PM2.5)为72.66 μg·m-3, 西部垣曲县污染最轻, ρ(PM2.5)为37.11 μg·m-3, 高值区PM2.5与SO2、NO2和CO呈显著强相关, 本地排放对高值区PM2.5影响较大; 季节分布特征与年分布有所差异, 春季中心城区污染较轻, 夏季ρ(PM2.5)均值最低, 为28.95 μg·m-3, 南部平陆县污染相对最重, 秋季中心城区转为高值区且ρ(PM2.5)相对最高, 可达62.94 μg·m-3, 冬季污染最为严重, ρ(PM2.5)分布在57.44~136.00 μg·m-3之间, 高值区与年分布一致, 集中在北部地区和中心城区.

(3) 输送路径分析结果表明, 西北向气团是影响运城市春季和冬季污染的主要气团, 其季节占比分别为58.69%和55.77%, ρ(PM2.5)最高可达82.72 μg·m-3和153.14 μg·m-3, 夏季主要受偏东和偏南气团的影响, 气团较为清洁, 相对污染程度较高的气团来自山西南部和河南西南部, 秋季来自河南西南部的气团污染占比最大, 为32.89%, ρ(PM2.5)最高, 可达120.22 μg·m-3.

(4) PSCF和CWT分析结果表明, 运城市PM2.5潜在源区存在明显的季节特征.主要贡献源区范围春夏最小, 秋季开始变大, 冬季范围最广且贡献最大.春季主要源区位于陕西东南部、河南中部和西南部与湖北接壤区域, 夏季主要位于河南中部和东部区域, 秋季主要位于河南西南部、陕西西部和东南部, 冬季则主要集中在陕西、河南和山西这3省交界区域以及河南中部和陕西西部.

参考文献
[1] Zhang Y, Shen L Y, Shuai C Y, et al. How is the environmental efficiency in the process of dramatic economic development in the Chinese cities?[J]. Ecological Indicators, 2019, 98: 349-362. DOI:10.1016/j.ecolind.2018.11.006
[2] Cao S S, Zhao W J, Guan H L, et al. Comparison of remotely sensed PM2.5 concentrations between developed and developing countries: Results from the US, Europe, China, and India[J]. Journal of Cleaner Production, 2018, 182: 672-681. DOI:10.1016/j.jclepro.2018.02.096
[3] Li Y, Liao Q, Zhao X G, et al. Premature mortality attributable to PM2.5 pollution in China during 2008-2016: Underlying causes and responses to emission reductions[J]. Chemosphere, 2021, 263. DOI:10.1016/J.chemosphere.2020.127925
[4] 杨文涛, 谯鹏, 刘贤赵, 等. 2011~2017年中国PM2.5多尺度时空分异特征分析[J]. 环境科学, 2020, 41(12): 5236-5244.
Yang W T, Qiao P, Liu X Z, et al. Analysis of multi-scale spatio-temporal differentiation characteristics of PM2.5 in China from 2011 to 2017[J]. Environmental Science, 2020, 41(12): 5236-5244.
[5] 李梦真, 张廷斌, 易桂花, 等. 四川盆地PM2.5浓度时空变化特征遥感监测与影响因子分析[J]. 环境科学, 2021, 42(7): 3136-3146.
Li M Z, Zhang T B, Yi G H, et al. Spatio-temporal variation characteristics monitored by remotely sensed technique of PM2.5 concentration and its influencing factor analysis in Sichuan basin, China[J]. Environmental Science, 2021, 42(7): 3136-3146.
[6] 张波, 宋国君, 周芳. 基于PM2.5监测点空间聚类的关中五市空气污染区域识别[J]. 环境科学学报, 2021, 41(3): 797-805.
Zhang B, Song G J, Zhou F. Regional delimitation of PM2.5 pollution: A case study of five cities in Guanzhong Plain[J]. Acta Scientiae Circumstantiae, 2021, 41(3): 797-805.
[7] 刘海猛, 方创琳, 黄解军, 等. 京津冀城市群大气污染的时空特征与影响因素解析[J]. 地理学报, 2018, 73(1): 177-191.
Liu H M, Fang C L, Huang J J, et al. The spatial-temporal characteristics and influencing factors of air pollution in Beijing-Tianjin-Hebei urban agglomeration[J]. Acta Geographica Sinica, 2018, 73(1): 177-191.
[8] Liu Y L, Shi G M, Zhan Y, et al. Characteristics of PM2.5 spatial distribution and influencing meteorological conditions in Sichuan Basin, southwestern China[J]. Atmospheric Environment, 2021, 253(12). DOI:10.1016/j.atmosenv.2021.118364
[9] 丁净, 唐颖潇, 郝天依, 等. 天津市冬季空气湿度对PM2.5和能见度的影响[J]. 环境科学, 2021, 42(11): 5143-5151.
Ding J, Tang Y X, Hao T Y, et al. Impact of air humidity on PM2.5 mass concentration and visibility in winter of Tianjin[J]. Environmental Science, 2021, 42(11): 5143-5151.
[10] Gao J J, Wang K, Wang Y, et al. Temporal-spatial characteristics and source apportionment of PM2.5 as well as its associated chemical species in the Beijing-Tianjin-Hebei region of China[J]. Environmental Pollution, 2018, 233: 714-724. DOI:10.1016/j.envpol.2017.10.123
[11] 肖致美, 徐虹, 李立伟, 等. 基于在线观测的天津市PM2.5污染特征及来源解析[J]. 环境科学, 2020, 41(10): 4355-4363.
Xiao Z M, Xu H, Li L W, et al. Characterization and source apportionment of PM2.5 based on the online observation in Tianjin[J]. Environmental Science, 2020, 41(10): 4355-4363.
[12] 赵孝囡, 王申博, 杨洁茹, 等. 郑州市PM2.5组分、来源及其演变特征[J]. 环境科学, 2021, 42(8): 3633-3643.
Zhao X N, Wang S B, Yang J R, et al. Chemical component and source of PM2.5 and their evolutive characteristics in Zhengzhou[J]. Environmental Science, 2021, 42(8): 3633-3643.
[13] 贺瑶, 韩秀秀, 黄晓虎, 等. 南京市不同功能区冬季大气PM2.5分布特征及其来源解析[J]. 环境科学学报, 2021, 41(3): 830-841.
He Y, Han X X, Huang X H, et al. Distribution characteristics and source apportionment of atmospheric PM2.5 in winter season from different functional areas of Nanjing[J]. Acta Scientiae Circumstantiae, 2021, 41(3): 830-841.
[14] 任娇, 尹诗杰, 郭淑芬. 太原市大气PM2.5中水溶性离子的季节污染特征及来源分析[J]. 环境科学学报, 2020, 40(9): 3120-3130.
Ren J, Yin S J, Guo S F. Seasonal variation and source analysis of water-soluble ions in PM2.5 in Taiyuan[J]. Acta Scientiae Circumstantiae, 2020, 40(9): 3120-3130.
[15] 李颜君, 安兴琴, 范广洲. 北京地区大气颗粒物输送路径及潜在源分析[J]. 中国环境科学, 2019, 39(3): 915-927.
Li Y J, An X Q, Fan G Z. Transport pathway and potential source area of atmospheric particulates in Beijing[J]. China Environmental Science, 2019, 39(3): 915-927. DOI:10.3969/j.issn.1000-6923.2019.03.003
[16] 刘娜, 余晔, 马学谦. 西宁市大气污染来源和输送季节特征[J]. 环境科学, 2021, 42(3): 1268-1279.
Liu N, Yu Y, Ma X Q. Seasonal Characteristics of air pollutant sources and transport pathways in Xining City[J]. Environmental Science, 2021, 42(3): 1268-1279.
[17] 雷雨, 张小玲, 康平, 等. 川南自贡市大气颗粒物污染特征及传输路径与潜在源分析[J]. 环境科学, 2020, 41(7): 3021-3030.
Lei Y, Zhang X L, Kang P, et al. Analysis of transport pathways and potential sources of atmospheric particulate matter in Zigong, in south of Sichuan Province[J]. Environmental Science, 2020, 41(7): 3021-3030.
[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, 40(1): 86-93.
Duan S G, Jiang N, Yang L M, et al. Transport pathways and potential sources of PM2.5 during the winter in Zhengzhou[J]. Environmental Science, 2019, 40(1): 86-93.
[20] Kong L W, Tan Q W, Feng M, et al. Investigating the characteristics and source analyses of PM2.5 seasonal variations in Chengdu, Southwest China[J]. Chemosphere, 2019, 243: 125267.
[21] 徐丹妮, 王瑾婷, 袁自冰, 等. 汾渭平原复杂地形影响下冬季PM2.5污染分布特征、来源及成因分析[J]. 环境科学学报, 2021, 41(4): 1184-1198.
Xu D N, Wang J T, Yuan Z B, et al. Temporal-spatial variations, source apportionment, and formation mechanisms of PM2.5 pollution over Fenwei Plain, China[J]. Acta Scientiae Circumstantiae, 2021, 41(4): 1184-1198.
[22] 郭庆元, 杨晓春, 吴其重, 等. 基于数值模拟对西安地区冬季一次重污染过程PM2.5区域来源解析[J]. 环境科学学报, 2020, 40(9): 3103-3111.
Guo Q Y, Yang X C, Wu Q Z, et al. An analysis of PM2.5 regional sources from heavy pollution in the winter of Xi'an based on numerical simulation[J]. Acta Scientiae Circumstantiae, 2020, 40(9): 3103-3111.
[23] 曹宁, 黄学敏, 祝颖, 等. 西安冬季重污染过程PM2.5理化特征及来源解析[J]. 中国环境科学, 2019, 39(1): 32-39.
Cao N, Huang X M, Zhu Y, et al. Pollution characteristics and source apportionment of fine particles during a heavy pollution in winter in Xi'an City[J]. China Environmental Science, 2019, 39(1): 32-39. DOI:10.3969/j.issn.1000-6923.2019.01.004
[24] 李雁宇, 李杰, 曾胜兰, 等. 2017年汾渭平原东部大气颗粒物污染特征分析[J]. 环境科学研究, 2020, 33(1): 63-72.
Li Y Y, Li J, Zeng S L, et al. Analysis of atmospheric particulates in the Eastern Fenwei Plain in 2017[J]. Research of Environmental Sciences, 2020, 33(1): 63-72.
[25] 黄小刚, 赵景波, 孙从建, 等. 汾渭平原PM2.5空间分布的地形效应[J]. 环境科学, 2021, 42(10): 4582-4592.
Huang X G, Zhao J B, Sun C J, et al. Orographic influences on the spatial distribution of PM2.5 on Fen-Wei Plain[J]. Environmental Science, 2021, 42(10): 4582-4592.
[26] 张海东. 西安市大气污染特征及颗粒物来源分析[D]. 西安: 西安建筑科技大学, 2020.
Zhang H D. Characteristics of air pollution and source analysis of particulate matter in Xi'an[D]. Xi'an: Xi'an University of Architecture and Technology, 2020.
[27] 黄小刚, 邵天杰, 赵景波, 等. 汾渭平原PM2.5浓度的影响因素及空间溢出效应[J]. 中国环境科学, 2019, 39(8): 3539-3548.
Huang X G, Shao T J, Zhao J B, et al. Influence factors and spillover effect of PM2.5 concentration on Fen-wei Plain[J]. China Environmental Science, 2019, 39(8): 3539-3548. DOI:10.3969/j.issn.1000-6923.2019.08.049
[28] 武媛媛, 李如梅, 彭林, 等. 运城市道路扬尘化学组成特征及来源分析[J]. 环境科学, 2017, 38(5): 1799-1806.
Wu Y Y, Li R M, Peng L, et al. Chemical compositions and source apportionment of road dust in Yuncheng[J]. Environmental Science, 2017, 38(5): 1799-1806.
[29] 赵清, 李杏茹, 王国选, 等. 运城秋冬季大气细粒子化学组成特征及来源解析[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.
[30] 李斌, 贾文毓. 运城市的自然村名与地形[J]. 山西师范大学学报(自然科学版), 2015, 29(S1): 59-60.
[31] 曹小虎. 运城市过去50a和未来气候变化趋势研究[J]. 人民黄河, 2014, 36(11): 13-15.
Cao X H. Analysis of characteristics and rules of climate change of Yuncheng City over the past 50 years and in the future[J]. Yellow River, 2014, 36(11): 13-15. DOI:10.3969/j.issn.1000-1379.2014.11.004
[32] Draxler R R, Hess G D. An overview of the HYSPLIT_4 modelling system for trajectories, dispersion, and deposition[J]. Australian Meteorological Magazine, 1998, 47(4): 295-308.
[33] Stein A F, Draxler R R, Rolph G D, et al. NOAA's HYSPLIT atmospheric transport and dispersion modeling system[J]. Bulletin of the American Meteorological Society, 2015, 96(12): 2059-2077. DOI:10.1175/BAMS-D-14-00110.1
[34] Draxler R R, Hess G D. Description of the HYSPLIT_4 modeling system[EB/OL]. https://www.arl.noaa.gov/wp_arl/wp-content/uploads/documents/reports/arl-224.pdf, 2021-04-16.
[35] Wang Y Q, Zhang X Y, Draxler R R. TrajStat: GIS-based software that uses various trajectory statistical analysis methods to identify potential sources from long-term air pollution measurement data[J]. Environmental Modelling & Software, 2009, 24(8): 938-939.
[36] 王芳, 陈东升, 程水源, 等. 基于气流轨迹聚类的大气污染输送影响[J]. 环境科学研究, 2009, 22(6): 637-642.
Wang F, Chen D S, Cheng S Y, et al. Impacts of air pollutant transport based on air trajectory clustering[J]. Research of Environmental Sciences, 2009, 22(6): 637-642.
[37] Sirois A, Bottenheim J W. Use of backward trajectories to interpret the 5-year record of PAN and O3 ambient air concentrations at Kejimkujik National Park, Nova Scotia[J]. Journal of Geophysical Research: Atmospheres, 1995, 100(2): 2867-2881.
[38] Draxler R R, Stunder B, Rolph G D, et al. HYSPLIT User's Guide[EB/OL]. http://www.arl.noaa.gov/documents/reports/hysplit_user_guide.pdf, 2021-04-16.
[39] Jiang N, Li L P, Wang S S, et al. Variation tendency of pollution characterization, sources, and health risks of PM2.5-bound polycyclic aromatic hydrocarbons in an emerging megacity in China: Based on three-year data[J]. Atmospheric Research, 2019, 217: 81-92. DOI:10.1016/j.atmosres.2018.10.023
[40] Dimitriou K, Kassomenos P. A meteorological analysis of PM10 episodes at a high altitude city and a low altitude city in central Greece-The impact of wood burning heating devices[J]. Atmospheric Research, 2018, 214: 329-337. DOI:10.1016/j.atmosres.2018.08.014
[41] Karaca F, Anil I, Alagha O. Long-range potential source contributions of episodic aerosol events to PM10 profile of a megacity[J]. Atmospheric Environment, 2009, 43(36): 5713-5722. DOI:10.1016/j.atmosenv.2009.08.005
[42] GB 3095-2012, 环境空气质量标准[S].
GB3095-2012, Ambient air quality standard[S].
[43] Zhang L, Shen F Z, Gao J L, et al. Characteristics and potential sources of black carbon particles in suburban Nanjing, China[J]. Atmospheric Pollution Research, 2020, 11(5): 981-991. DOI:10.1016/j.apr.2020.02.011
[44] Polissar A V, Hopke P K, Harris J M. Source regions for atmospheric aerosol measured at Barrow, Alaska[J]. Environmental Science & Technology, 2001, 38(21): 4214-4226.
[45] Xu X, Akhtar U S. Identification of potential regional sources of atmospheric total gaseous mercury in Windsor, Ontario, Canada using hybrid receptor modeling[J]. Atmospheric Chemistry and Physics, 2010, 10(15): 7073-7083. DOI:10.5194/acp-10-7073-2010
[46] Dvorská A, Lammel G, Holoubek I. Recent trends of persistent organic pollutants in air in central Europe-Air monitoring in combination with air mass trajectory statistics as a tool to study the effectivity of regional chemical policy[J]. Atmospheric Environment, 2009, 43(6): 1280-1287. DOI:10.1016/j.atmosenv.2008.11.028
[47] Hsu Y K, Holsen T M, Hopke P K. Comparison of hybrid receptor models to locate PCB sources in Chicago[J]. Atmospheric Environment, 2003, 37(4): 545-562. DOI:10.1016/S1352-2310(02)00886-5
[48] 李杰, 翟亮, 桑会勇, 等. PM2.5浓度插值中不同空间插值方法对比[J]. 测绘科学, 2016, 41(4): 50-54, 101.
Li J, Zhai L, Sang H Y, et al. Comparison of different spatial interpolation methods for PM2.5[J]. Science of Surveying and Mapping, 2016, 41(4): 50-54, 101.
[49] 郭雯雯, 陈永金, 刘阁, 等. 2016~2019年长江中游城市群空气质量时空变化特征及影响因素分析[J]. 生态环境学报, 2020, 29(10): 2034-2044.
Guo W W, Chen Y J, Liu G, et al. Analysis on the characteristics and influencing factors of air quality of urban agglomeration in the middle reaches of the Yangtze River in 2016 to 2019[J]. Ecology and Environment Sciences, 2020, 29(10): 2034-2044.
[50] 任浦慧, 解静芳, 姜洪进, 等. 太原市大气PM2.5季节传输路径和潜在源分析[J]. 中国环境科学, 2019, 39(8): 3144-3151.
Ren P H, Xie J F, Jiang H J, et al. Transport pathways and potential sources of PM2.5 in different seasons in Taiyuan city[J]. China Environmental Science, 2019, 39(8): 3144-3151. DOI:10.3969/j.issn.1000-6923.2019.08.002
[51] 沈楠驰, 周丙锋, 李珊珊, 等. 2015~2019年天津市大气污染物时空变化特征及成因分析[J]. 生态环境学报, 2020, 29(9): 1862-1873.
Shen N C, Zhou B F, Li S S, et al. Temporal and spatial variation characteristics and origin analysis of air pollutants in Tianjin from 2015 to 2019[J]. Ecology and Environment Sciences, 2020, 29(9): 1862-1873.