快速检索
  气象   2023, Vol. 49 Issue (11): 1359-1370.  DOI: 10.7519/j.issn.1000-0526.2023.052401

论文

引用本文 [复制中英文]

周涛, 周青, 张勇, 等, 2023. 汾渭平原PM2.5-O3复合污染特征及气象成因分析[J]. 气象, 49(11): 1359-1370. DOI: 10.7519/j.issn.1000-0526.2023.052401.
[复制中文]
ZHOU Tao, ZHOU Qing, ZHANG Yong, et al, 2023. Characteristics of PM2.5-O3 Compound Pollution and Meteorological Impact in Fenwei Plain[J]. Meteorological Monthly, 49(11): 1359-1370. DOI: 10.7519/j.issn.1000-0526.2023.052401.
[复制英文]

资助项目

山西省气象局科学技术研究重点项目(SXKZDTC20226302)、中国环境科学研究院重点课题(DQGG2021101)和秦岭和黄土高原生态环境气象重点实验室开放研究基金重点课题(2021K-10)共同资助

第一作者

周涛, 主要从事天气预报和环境气象工作.E-mail: 805746177@qq.com

通信作者

周青, 主要从事大气成分与环境气象研究.E-mail: zhouqing@cma.gov.cn.

文章历史

2022年10月27日收稿
2023年5月10日收修定稿
汾渭平原PM2.5-O3复合污染特征及气象成因分析
周涛 1, 周青 2, 张勇 3, 吴昱树 4, 孙健 5    
1. 山西省运城市气象局, 运城 044000
2. 中国气象局气象探测中心, 北京 100081
3. 中国气象局气象发展与规划院, 北京 100081
4. 乌兰察布市气象局, 乌兰察布 012000
5. 莒县气象局, 莒县 276500
摘要:PM2.5和O3已经成为汾渭平原城市最主要的污染物,两者之间相互影响,在暖季经常同时出现构成污染,其污染程度与气象条件密切相关。利用2015—2021年汾渭平原12个城市逐日PM2.5和O3浓度、地面气象观测数据以及ERA5高空再分析数据等资料,分析了汾渭平原PM2.5和O3的时空变化特征以及复合污染发生时PM2.5和O3的关系,并研究了局地气象条件和天气形势对复合污染的影响。结果显示,该地区年均PM2.5和日最大8小时O3浓度分别在2017年和2018年开始持续下降,复合污染日数也在2019年后开始持续下降;复合污染主要发生在3—9月,在汾渭平原东部城市出现次数较多,多出现在高温、低湿的环境下;最后利用T-PCA算法(正交主成分分析)将复合污染的天气环流形势分为4种类型,主要呈现出以高空西北气流或偏西气流、低层为暖区偏南风或微风为主的天气特征。研究结果对汾渭平原的大气污染防治提供了参考依据,尤其对于PM2.5和O3污染协同治理具有重要意义。
关键词汾渭平原    PM2.5和O3复合污染    T-PCA分型    气象成因    
Characteristics of PM2.5-O3 Compound Pollution and Meteorological Impact in Fenwei Plain
ZHOU Tao1, ZHOU Qing2, ZHANG Yong3, WU Yushu4, SUN Jian5    
1. Yuncheng Meteorological Office of Shanxi Province, Yuncheng 044000;
2. CMA Meteorological Observation Centre, Beijing 100081;
3. Institute of Meteorological Development and Planning, CMA, Beijing 100081;
4. Wulanchabu Meteorological Office of Inner Mongolia, Wulanchabu 012000;
5. Juxian Meteorological Station of Shanxi Province, Juxian 276500
Abstract: PM2.5 and O3 which interact with each other have become the most important pollutants in the cities of Fenwei Plain. The pollution often occurs simultaneously in warm seasons and is closely related to meteorological conditions. Based on the daily PM2.5 and O3 concentration data, surface meteorological observation data and ERA5 high-altitude reanalysis data of 12 cities in the Fenwei Plain from 2015 to 2021, the spatio-temporal characteristics of PM2.5 and O3 in the Fenwei Plain, as well as the relationship between PM2.5 and O3 when compound pollution occurs are analyzed. Besides, the impact of local meteorological conditions and synoptic situation on compound pollution is also studied. The results show that the average annual PM2.5 and MDA8 O3 (maximum daily 8 h average ozone) concentrations in the Fenwei Plain began to decline continuously in 2017 and 2018, respectively, and the number of compound pollution days also began to decrease continuously after 2019. Compound pollution mainly occurs in March to September, frequently seen in the cities of eastern Fenwei Plain, mostly in condition of high temperature and low humidity. Finally, the synoptic circulation situation of compound pollution is divided into four types by principal component analysis with the T-PCA algorithm, featured with the northwest or westerly air flow in high-altitude and the southerly wind or breeze in the low-level warm zone. These findings could provide some basis for controlling pollution in the Fenwei Plain and are of great significance for cooperative governance of PM2.5 and O3 pollution.
Key words: Fenwei Plain    PM2.5 and O3 compound pollution    T-PCA typing    meteorological impact    
引言

汾渭平原是我国第四大平原,由汾河、渭河冲击形成,由于其大部分城市位于黄土高原,地形封闭,能源消费结构以煤炭为主,年消耗量达2.6亿t,占全国能源消费的80%以上(Feng et al,2016),加上繁重的道路运输和较高的人口密度,导致汾渭平原的污染排放量极高(Dao et al,2021),大气污染严重、生态环境脆弱(黄小刚等,2021; 王妘涛等,2022; 秦卓凡等,2021; 南洋等,2021张天航等,2020),并且可能加重邻近地区的污染(黄小刚等,2019)。2013年以来全国空气质量明显好转(Zhang et al,2018a2018b),但汾渭平原空气质量改善并不明显(Feng et al,2019),并且已经取代京津冀地区成为中国污染最严重的地区(Dao et al,2021; 杨乐超等,2018)。2018年后随着汾渭平原减排措施的强制实施(秦卓凡等,2021),汾渭平原的PM2.5浓度已经开始下降(郝永佩等,2022),但O3污染却越来越严重(关茜妍等,2021)。2021年汾渭平原污染物超标天数比例为29.8%,PM2.5污染和O3污染分别占总超标天数的39.3%和38.0%,污染形势依旧严峻(中华人民共和国生态环境部,2022)。

污染物局地排放是大气污染的主要因素,而不利的气象因素会加剧空气污染(Wang et al,2018; 祁海霞等, 2019; 陈镭等, 2020; 胡春梅等, 2020; 宋佳琨等, 2021; 刘兆东等, 2020; 杨颖川等, 2022),2014—2018年汾渭平原的PM2.5质量浓度(以下简称浓度)呈4.7 μg·m-3·a-1的趋势下降,气象因素的贡献为0.7 μg·m-3·a-1,贡献率达15%,O3浓度呈6.6 μg·m-3·a-1的趋势增长,气象因素贡献率占12% (Chen et al,2020),有研究显示汾渭平原2015—2019年各年度臭氧单因子所导致污染超标所占比例分别为1.4%、5.4%、13.0%、11.1%、14.4%(解淑艳等,2021)。O3污染事件通常与高温低湿条件有关(何慧根等, 2021; 杨镇江等, 2023),高温和低湿环境导致大量生物挥发性有机物(BVOC)的排放促进光化学反应(Lu et al,2019),进而加重地面O3污染;而PM2.5污染通常发生在弱风、高湿和低边界层的条件下(Xu et al,2018; 张琴和姚秀萍, 2021),两者条件虽有不同,但高浓度的O3会导致二次PM2.5粒子产生,使得空气中的PM2.5浓度增大(Ding et al,2013; Wang et al,2016),而高浓度PM2.5会通过吸收太阳辐射、改变光解速率从而降低O3浓度(Unger et al,2009; Menon et al,2008; Li et al,2018),两者之间会相互影响,因此在许多城市都出现了PM2.5和O3同时污染的现象(He et al,2022; Schnell and Prather, 2017; Luo et al,2022; Yang et al,2021徐栋夫等,2020),我国的大气污染已经呈现出以PM2.5和O3为主的复合污染特征(Zhang et al,2008),这种复合污染的现象会极大增加人类健康风险(Anenberg et al,2010)。

综上所述,以前的学者对汾渭平原污染的气象成因研究大多是针对单一PM2.5或O3污染,对汾渭平原PM2.5和O3复合污染时的气象成因分析较少。因此我们将从气温、湿度和风速等气象要素及天气形势方面分析PM2.5和O3复合污染成因,以期为汾渭平原大气复合污染防治以及环境气象预报服务等提供参考依据。

1 资料与方法 1.1 研究区域

此次研究区域汾渭平原城市群共有12个城市(图 1),包括陕西省西安市、渭南市、铜川市、咸阳市和宝鸡市,山西省太原市、晋中市、吕梁市、临汾市和运城市,河南省三门峡市和洛阳市。2021年这12个城市中有6个位于全国污染倒数20名之内(中华人民共和国生态环境部,2022)。

图 1 汾渭平原各城市地理位置分布及海拔高度信息 注:陕西省:蓝色,山西省:黑色,河南省:红色。 Fig. 1 Geographic location distribution and elevation information of cities in Fenwei Plain
1.2 污染物和气象数据

污染物数据包括2015—2021年12个城市的PM2.5和MDA8 O3(日最大8小时平均O3浓度)观测资料,来源于中国环境监测总站(https://quotsoft.net/air/)(Wang et al,2021; Jiang et al,2021; Zhang et al,2022);地面气象数据来自中国气象局国家气象数据观测中心(http://data.cma.cn/),包括同时期逐小时的地面风速(WS)、风向(WD)、气温(T)、降水(RAIN)、相对湿度(RH);此外还使用了欧洲中期天气预报中心第5代全球大气再分析数据(ERA5),包括位势高度、风向、风速,时间分辨率为1 h,空间分辨率为0.25°×0.25°。

1.3 污染日定义

本研究重点关注PM2.5和O3复合污染时的气象条件。世界卫生组织在综合考虑过去15年间环境健康研究进展的基础上,制定了2021版空气质量指南(朱彤等,2022),针对PM2.5和O3提出了更为严格的目标。根据中国环境空气质量标准(生态环境部,2012),将同时满足日均PM2.5>35 μg·m-3 (一级浓度限值)和MDA8 O3>100 μg·m-3(一级浓度限值)两个条件时的日期定义为复合污染日,同时满足PM2.5>35 μg·m-3和MDA8 O3 < 100 μg·m-3定义为PM2.5污染日,同时满足PM2.5 < 35 μg·m-3和MDA8 O3>100 μg·m-3定义为O3污染日,同时满足PM2.5 < 35 μg·m-3和MDA8 O3 < 100 μg· m-3定义为无污染日。

1.4 研究方法

自2013年政府发布《大气污染防治行动计划》以来,中国的空气污染物浓度呈下降趋势(Zhang et al,2018a),首要污染物以PM2.5和O3为主。要研究每日PM2.5和O3的关系,必须消除PM2.5和O3的长期变化趋势所产生的影响,因此我们采用去趋势化和去季节化浓度的处理方式,Tai et al(2010)Dai et al(2021)均用此方法进行研究,去趋势化算法如下:

$ \operatorname{Dev}_{m i}=C_{m i}-\sum\limits_{i=1}^n \frac{C_{m i}}{n} $ (1)

式中: Devmi为去趋势化后的污染物浓度,Cmim月份第i天的污染物浓度,nm月份的天数,$\sum\limits_{i=1}^n C_{m i} / n$m月污染物的平均浓度,因此当某日的Devmi为正时,说明其浓度高于当月的平均浓度。

在天气分型方面,使用正交主成分分析方法(T-PCA)来对汾渭平原的天气形势进行分类。T-PCA作为一种可靠的分类方法(Ye et al,2016; Miao et al,2018; Huth et al,2008),可以准确地反映原始环流场的特征,不会因分型对象的调整而有太大变化,得到的时空场也更加稳定(Huth et al,2008)。使用cost733 class分型软件(http://cost733.met.no)中的T-PCA算法对研究区域850 hPa位势高度场进行客观分型,进而识别出易受污染的天气形势。此外,通过计算不同天气分型数量下的累积解释方差(ECV)以及累积解释方差的变率(ΔECV),从而完成对天气分型性能的评估,其中ΔECV=ECVk+1-ECV k,其中k为分类数量,ΔECV值越大表示在天气模式分类中的性能提升幅度越大,并且趋于稳定(Ning et al, 2019)。

2 结果和讨论 2.1 2015—2021年PM2.5和O3浓度变化趋势

图 2为2015—2021年汾渭平原各城市PM2.5和MDA8 O3年均浓度空间分布,其中O3年评价值采用日最大8小时平均浓度的第90百分位数浓度值计算得到。整体来看汾渭平原PM2.5年均浓度于2015—2017年呈上升趋势,2017年后开始下降,由2017年的66.8 μg·m-3下降至2021年的44.5 μg· m-3,下降幅度达33.4%,西安、渭南、咸阳、临汾、运城和洛阳等城市浓度下降明显,但浓度仍远高于世界卫生组织推荐的年均标准(朱彤等,2022);虽然2013年大气污染防治行动计划要求控制生成O3的氮氧化物和挥发性有机物等的排放,但MDA8 O3年评价值仍由2015年的131.1 μg·m-3升至2017年的183.4 μg·m-3,上升幅度达39.9%,说明氮氧化物和挥发性有机物尚未得到有效控制,2018年加大整治力度后缓慢下降。通过统计汾渭平原PM2.5和MDA8 O3逐年月平均浓度变化趋势(图略),发现PM2.5浓度主要在11月至次年2月呈现高值,而MDA8 O3则在4—9月呈现高值,这些特征与京津冀等区域的污染情况类似(Ou et al,2022)。

图 2 2015—2021年汾渭平原各城市(a)PM2.5和(b)MDA8 O3年均浓度空间分布 注:括号内数字为相应的年评价值。 Fig. 2 Spatial distributions of the annual average concentrations of (a) PM2.5 and (b) the 90th quantile in MDA8 O3 in cities of Fenwei Plain from 2015 to 2021

图 3为2015—2021年汾渭平原各城市PM2.5污染、O3污染和复合污染天数的空间分布,可以看出汾渭平原7年来各城市PM2.5污染天数(图 3a)介于852~1132 d,其中污染天数最多的地区是咸阳市;O3污染天数(图 3b)则介于328~671 d,O3污染天数最多的地区是铜川市,主要是由于铜川市第二产业占比较高,O3前体物排放较其他区域突出(郑小华等,2021);从复合污染情况(图 3c)看,各城市污染天数介于332~737 d,运城和临汾的污染天数较多,占总天数28%,东部城市太原、晋城、临汾、运城和三门峡普遍在600 d以上,西部西安、宝鸡、咸阳、渭南和铜川则均在600 d以下,东部的运城、临汾和三门峡、洛阳复合污染比西部城市严重,主要是由于这几个城市处于两省交界处,也可能与政策缺乏协同规划有关(郑小华等,2021)。从各城市逐年复合污染次数二维分布来看(图 3d),2017年和2018年为复合污染天数最多的两年,2019—2021年复合污染天数开始逐年下降,表明O3和PM2.5协同治理已经发挥了作用,2021年西部城市西安、宝鸡、咸阳、渭南和铜川的复合污染天数已经降至40 d以下,洛阳和三门峡分别为52 d和57 d,而太原、晋中、临汾和运城的天数仍在70 d以上,因此这几个城市仍需加大协同治理力度。

图 3 2015—2021年汾渭平原(a)PM2.5, (b)MDA8 O3, (c)复合污染总天数及(d)各城市逐年复合污染次数二维分布 Fig. 3 Total number of days of (a) PM2.5, (b) MDA8 O3, (c) compound pollution and (d) two-dimensional distribution of times of yearly compound pollution days in Fenwei Plain from 2015 to 2021
2.2 汾渭平原地区PM2.5和MDA8 O3的关系

统计了汾渭平原2015—2021年各城市逐月累计复合污染天数二维分布情况(图略),各城市复合污染主要集中在3—9月,10月至次年2月PM2.5污染严重,但冬季MDA8 O3浓度偏低。使用式(1)计算了每日的Dev_PM2.5(去趋势后的PM2.5)和Dev_MDA8 O3(去趋势后的MDA8 O3),从二者2015—2021年不同季节的相关系数箱线图(图 4)看,各城市春夏秋三季多以正相关为主,其中夏季相关性明显高于春秋两季,春季在2015—2016年后呈现波动趋势,夏季和秋季均在2019—2021年呈现上升的趋势,说明夏秋两季O3对PM2.5浓度的影响呈现加剧的态势;而两者在冬季多呈现负相关,且自2015年开始负相关逐渐减小,到2021年各城市的相关性均值接近于0,说明冬季PM2.5对O3的影响在持续减弱。

图 4 2015—2021年汾渭平原不同季节各城市Dev_PM2.5和Dev_MDA8 O3相关系数箱线图 注:箱体上下边界分别表示第75、第25百分位数,黑色圆点表示平均值,中间横线表示中位数,空心圆表示异常值。 Fig. 4 Boxplots of DEV_PM2.5 and DEV_MDA8 O3 correlation coefficients of each city in different seasons in Fenwei Plain from 2015 to 2021

由于复合污染主要集中在3—9月,因此重点分析3—9月每日PM2.5和MDA8 O3的关系。统计发现,2015—2021年,汾渭平原PM2.5呈下降趋势,而MDA8 O3呈上升趋势,且去趋势化后的PM2.5与MDA8 O3具有弱正相关性(相关系数为0.21)。

图 5a为各城市PM2.5浓度偏大时(即Dev_PM2.5>0 μg·m-3),Dev_PM2.5和Dev_MDA8 O3的相关系数。可以看出4、5、6和9月基本呈负相关,表明当PM2.5浓度偏大时,会制约影响O3的生成,但7月和8月太阳辐射较强、气温偏高、大气扩散条件好,因此PM2.5等颗粒污染物对太阳辐射的削弱作用有限,因此7—8月大部分城市的Dev_PM2.5和Dev_MDA8 O3为正相关。而当MDA8 O3浓度偏大(即Dev_MDA8 O3>0 μg·m-3)时(图 5b), 大气氧化作用强,使得一次气态污染物SO2、NOx和NH3转换成硫酸盐和硝酸盐气溶胶粒子,加重PM2.5污染(王芃等, 2022),因此在3—9月所有城市基本都表现为正相关,尤其是高海拔的太原、晋中(在4月和5月相关系数高达0.5)。

图 5 2015—2021年3—9月汾渭平原(a)Dev_PM2.5浓度偏大时和(b)Dev_MDA8 O3浓度偏大时各城市Dev_PM2.5和Dev_MDA8 O3的相关系数 Fig. 5 Correlation coefficients between Dev_PM2.5 and Dev_MDA8 O3 for each city in Fenwei Plain when (a) Dev_PM2.5>0 μg·m-3 and (b) Dev_MDA8 O3>0 μg·m-3 from March to September from 2015 to 2021

图 6为复合污染时Dev_MDA8 O3与Dev_PM2.5的散点关系,图中展示了当Dev_MDA8 O3在不同取值区间时Dev_PM2.5>0 μg·m-3的样本所占比例。当Dev_MDA8 O3>0 μg·m-3时,Dev_PM2.5>0 μg·m-3的比例均能超过69%,并且随着Dev_MDA8 O3的增大呈上升趋势,在Dev_MDA8 O3>40 μg·m-3时,汾渭平原12个城市的Dev_PM2.5>0 μg·m-3的比例超过80%,表明了大气氧化对PM2.5浓度增加起到积极影响。然而大气强氧化性并不能无限制导致PM2.5浓度的增长,一方面受制于硫酸盐、硝酸盐等粒子的数量,另一方面由于O3的生成特性要求大气环境为高温和强辐射,而高温条件下会导致硫酸盐、硝酸盐等粒子的挥发和扩散,进而导致PM2.5浓度的下降(张远远等, 2022)。

图 6 2015—2021年3—9月汾渭平原Dev_MDA8 O3与Dev_PM2.5的散点图 注:不同颜色代表Dev_MDA8 O3不同取值区间,上横坐标各刻度值表示当Dev_MDA8 O3在各取值区间时Dev_PM2.5>0 μg·m-3的样本所占比例。 Fig. 6 Scatter plot of Dev_MDA8 O3 and Dev_PM2.5in Fenwei Plain from March to September from 2015 to 2021
2.3 气象条件对污染的影响

除排放因素外,PM2.5和O3浓度也受到气象条件的影响(何国文等, 2022; 肖致美等, 2022; 花丛等, 2022),图 7为汾渭平原12个城市3—9月中发生不同污染时的气温(T)、相对湿度(RH)和风速(W)的箱线图,可以看出O3污染和复合污染发生在高温、低湿的环境条件下,日平均气温主要介于18.0~26.5℃,单独O3污染时较复合污染更高,日平均相对湿度主要介于46.3%~59.0%;而PM2.5污染和无污染发生在气温相对较低、湿度偏高的条件下,PM2.5污染时气温介于8.6~18.6℃,相对湿度在46.5%~80.3%,相对湿度范围大于O3和复合污染,而无污染时气温介于13.4~21.5℃,高于PM2.5污染,但相对湿度介于53%~89%,略高于污染时的水平。在风速上4种类型差别不明显,均介于1.2~2.4 m·s-1。高温利于大气O3的光化学反应和二次污染物的生成,而高的相对湿度则利于PM2.5的吸湿增长(闫小利等, 2021),低的风速不利于污染物的扩散(肖致美等, 2022)。无污染时表现的低温、高湿、弱风气象特征则可能由于PM2.5和O3排放较少或降雨的影响。降雨可以有效地清除PM2.5污染,并且降雨、降温和较厚云层阻挡太阳辐射不利于O3产生,统计发现无污染天气(2385 d)中汾渭平原降雨(超过6个城市出现降雨)的天数为1264 d,占比为53%。由于PM2.5和无污染样本在3种气象要素下不易区分,剔除无污染中的降水样本后,无污染和PM2.5污染样本主要分布在3、4和9月,对比两种情况在3种气象要素下的一阶差分的分布状况(后一日要素值减去前一日的要素值,图略),发现无污染时气温与前一日相比会略有降低,而PM2.5污染时则温度均值会略有上升,这主要是由于春秋季冷空气活动利于污染物扩散;而相对湿度方面,无污染时相对湿度则会明显降低,而PM2.5污染时相对湿度则不明显,这主要是由于低的相对湿度不利于污染物吸湿增长;风速方面两者呈现出不明显的区别,这主要因为在不同风向的传输影响下可能呈现出污染加重或减轻两种现象。

图 7 2015—2021年3—9月汾渭平原各污染类型下(a)2 m气温、(b)2 m相对湿度、(c)10 m风速的箱线分布 注:横坐标标注和不同颜色分别代表单独的PM2.5污染,单独O3污染,复合污染,无污染。 Fig. 7 Boxplots showing distribution of (a) 2 m air temperature, (b) 2 m relative humidity, and (c) 10 m wind speed related to each pollution type in Fenwei Plain from March to September from 2015 to 2021

为了更好地了解复合污染天气发生时的天气形势,提取了2015—2021年3—9月汾渭平原复合污染时的日期作为统计样本,这些污染日中有一半以上的城市都出现了复合污染的现象,使用T-PCA算法将复合污染发生时的850 hPa位势高度进行天气分型,同时计算了2~12种天气分型结果的ECV以及ΔECV,结果显示当天气分型数量为4时ΔECV取最大值,因此这里根据统计样本分为4种天气形势,绘制出对应的500 hPa形势图(图 8)。图 9为对应的复合污染时每种天气类型在各月的分布比例情况,分型方法详见1.4节。可以看到,类型1为最主要的污染天气形势,占总天数的62.8%,平均PM2.5和MDA8 O3浓度分别为55 μg·m-3和138 μg·m-3,此种类型下PM2.5浓度最高,O3浓度最低,这是由于该类型多出现在3—5月,气温较其余类型偏低,此时500 hPa副热带高压偏南,汾渭平原受槽后脊前的西北风影响,导致天气晴朗,850 hPa受暖脊控制,以弱偏南风为主,利于污染物的输送,但温度偏低在一定程度上影响了O3的产生。类型2占比为15.2%,主要发生在7—8月,平均PM2.5和MDA8 O3浓度分别为47 μg·m-3和144 μg·m-3,PM2.5污染程度较类型1偏弱,O3浓度则相对偏高,500 hPa副热带高压较类型1偏北,退居海上,脊线位于25°N附近,汾渭平原上空气流平直,而850 hPa与类型1类似,以东南风为主,同样存在暖湿气流输送,但气压梯度偏大,风力大于类型1,加之7—8月温度较高,利于污染物湍流扩散,导致PM2.5浓度低于类型1而O3浓度高于类型1。类型3占比为15.2%,主要发生在6月和7月,平均PM2.5和MDA8 O3浓度分别为45 μg·m-3和152 μg·m-3,此类型O3污染最强,这是一年中O3浓度最高的月份,500 hPa汾渭平原上空为槽后西北风,850 hPa汾渭平原同样为东南风,虽然6月和7月污染物排放少,但高的O3浓度会使得大气氧化性增强,提高PM2.5浓度。类型4占比为6.8%,主要在7—9月,平均PM2.5和MDA8 O3浓度分别为44 μg·m-3和142 μg·m-3,副热带高压偏东,汾渭平原上空为槽后弱西北气流控制,地面气压场弱,风速小,利于PM2.5和O3前体物的累积。综合分析发现,汾渭平原污染时500 hPa以槽后西北风和偏西风为主,而850 hPa多受暖脊的影响,以偏南风或静风为主,槽后西北风和偏西风的动力作用使天气以晴朗活动云为主,利于太阳辐射,低层偏南暖湿气流输送或者静风使得污染物累积,导致发生PM2.5和O3复合污染。

图 8 2015—2021年3—9月汾渭平原50%以上城市出现复合污染时所对应的4种天气形势(a)500 hPa位势高度分型图和风场(风矢),(b)每种类型所对应的850 hPa位势高度形势图和风场(风矢) 注:黑框代表汾渭平原地区。 Fig. 8 Four synoptic patterns corresponding to compound pollution at the time for compound pollution to occur in more than half of cities in Fenwei Plain from March to September from 2015 to 2021 (a) 500 hPa geopotential height (colored) superposed with the wind (wind vector), (b) 850 hPa geopotential height corresponding to each type with the wind (vector)

图 9 2015—2021年3—9月汾渭平原复合污染时不同类型天气形势的分布 Fig. 9 Distribution of different synoptic patterns during compound pollutionin Fenwei Plain from March to September from 2015 to 2021
3 结论

本文基于2015—2021年间汾渭平原12个城市逐日PM2.5和O3浓度、地面气象观测数据以及ERA5高空再分析数据等资料,对汾渭平原7年来PM2.5和O3的时空变化特征尤其是复合污染特征开展了统计分析,并重点研究了复合污染发生时PM2.5和O3二者的相关性,局地气象条件以及天气形势对复合污染的影响,主要结论如下:

汾渭平原的PM2.5和MDA8 O3年均浓度均在2017年达到近年来的峰值,得益于大气污染防治行动计划和蓝天保卫行动计划的实施,其年评价值(PM2.5为年均浓度,O3为日最大8小时滑动平均值第90百分位数) 分别呈5.6、3.7 μg·m-3·a-1的速率下降。汾渭平原各城市复合污染天数介于332~ 737 d,东部城市污染日数多于西部,运城复合污染日数占比达30.7%;复合污染日数在2016年开始增加,2018年达到峰值,2019—2021年呈下降趋势,体现了政府在污染协同治理方面的成效。

去趋势化后的PM2.5和MDA8 O3存在弱正相关性,Dev_MDA8 O3>0 μg·m-3时,Dev_PM2.5>0 μg·m-3的比例均能超过69%,Dev_MDA8 O3>40 μg·m-3时Dev_PM2.5>0 μg·m-3的比例超过80%,表明了大气氧化对增加颗粒物浓度所起到的正贡献影响。

O3污染和复合污染发生在高温、低湿的环境条件下,PM2.5污染和无污染发生在气温相对较低、湿度偏高和风速较小的条件下,无污染时出现降雨的概率为53%,体现了降雨对污染物的湿清除作用。使用T-PCA算法将复合污染时的天气形势分为4种类型,这些形势的共同特点是500 hPa位势高度场以西北或偏西气流为主,低层850 hPa位势高度场受暖脊影响,以偏南风或静风为主,利于O3的生成和PM2.5的传输和累积。

参考文献
陈镭, 周广强, 毛卓成, 等, 2020. 上海地区2017年三次短时重度污染过程变化特征及其气象影响因素[J]. 气象, 46(5): 675-686. Chen L, Zhou G Q, Mao Z C, et al, 2020. Variation characteristics and meteorological impact factors of three short-time severe air pollutions in Shanghai in 2017[J]. Meteor Mon, 46(5): 675-686 (in Chinese).
关茜妍, 陆克定, 张宁宁, 等, 2021. 西安市大气臭氧污染光化学特征与敏感性分析[J]. 科学通报, 66(35): 4561-4573. Guan X Y, Lu K D, Zhang N N, et al, 2021. Analysis of the photochemical characteristics and sensitivity of ozone pollution in Xi'an[J]. Chin Sci Bull, 66(35): 4561-4573 (in Chinese).
郝永佩, 宋晓伟, 赵文珺, 等, 2022. 汾渭平原大气污染时空分布及相关因子分析[J]. 生态环境学报, 31(3): 512-523. Hao Y P, Song X W, Zhao W J, et al, 2022. Spatiotemporal distribution of air pollution and correlation factors in Fenwei Plain[J]. Ecol Environ Sci, 31(3): 512-523 (in Chinese). DOI:10.16258/j.cnki.1674-5906.2022.03.010
何国文, 邓涛, 欧阳珊珊, 等, 2022. 广州地区秋季PM2.5和臭氧复合污染的观测研究[J]. 环境科学学报, 42(6): 250-259. He G W, Deng T, Ouyang S S, et al, 2022. Observation studies on the PM2.5 and O3 complex episodes during autumn in Guangzhou[J]. Acta Sci Circumstantiae, 42(6): 250-259 (in Chinese).
何慧根, 唐红玉, 李永华, 等, 2021. 2014—2018年重庆主城区大气污染的特征及其与大气环流之间的关系[J]. 气象, 47(10): 1233-1245. He H G, Tang H Y, Li Y H, et al, 2021. Characteristics of air pollution and its relationship with atmospheric circulation in Chongqing City from 2014 to 2018[J]. Meteor Mon, 47(10): 1233-1245 (in Chinese). DOI:10.7519/j.issn.1000-0526.2021.10.006
胡春梅, 陈道劲, 周国兵, 等, 2020. 基于自组织神经网络算法的重庆秋冬季空气污染与天气分型的关系[J]. 气象, 46(9): 1222-1234. Hu C M, Chen D J, Zhou G B, et al, 2020. Relationship between air pollution events in autumn and winter in Chongqing and the classification of synoptic situation based on self-organizing maps[J]. Meteor Mon, 46(9): 1222-1234 (in Chinese).
花丛, 江琪, 迟茜元, 等, 2022. 我国中东部地区2015—2020年夏半年PM2.5和臭氧复合污染气象特征分析[J]. 环境科学研究, 35(3): 650-658. Hua C, Jiang Q, Chi X Y, et al, 2022. Meteorological characteristics of PM2.5-O3 air combined pollution in central and eastern China in the summer half years of 2015-2020[J]. Res Environ Sci, 35(3): 650-658 (in Chinese).
黄小刚, 邵天杰, 赵景波, 等, 2019. 汾渭平原PM2.5浓度的影响因素及空间溢出效应[J]. 中国环境科学, 39(8): 3539-3548. Huang X G, Shao T J, Zhao J B, et al, 2019. Influence factors and spillover effect of PM2.5 concentration on Fen-Wei Plain[J]. China Environ Sci, 39(8): 3539-3548 (in Chinese). DOI:10.3969/j.issn.1000-6923.2019.08.049
黄小刚, 赵景波, 孙从建, 等, 2021. 汾渭平原PM2.5空间分布的地形效应[J]. 环境科学, 42(10): 4582-4592. Huang X G, Zhao J B, Sun C J, et al, 2021. Orographic influences on the spatial distribution of PM2.5 on the Fen-Wei Plain[J]. Environ Sci, 42(10): 4582-4592 (in Chinese).
刘兆东, 王宏, 沈新勇, 等, 2020. 京津冀及周边地区冬季能见度与PM2.5浓度和环境湿度的多元回归分析[J]. 气象学报, 78(4): 679-690. Liu Z D, Wang H, Shen X Y, et al, 2020. Multiple regression analysis of winter visibility, PM2.5 concentration and humidity in Beijing-Tianjin-Hebei and its surrounding regions[J]. Acta Meteor Sin, 78(4): 679-690 (in Chinese).
南洋, 饶晓琴, 尤媛, 等, 2021. 2020年11月大气环流和天气分析[J]. 气象, 47(2): 253-260. Nan Y, Rao X Q, You Y, et al, 2021. Analysis of the November 2020 atmospheric circulation and weather[J]. Meteor Mon, 47(2): 253-260 (in Chinese).
祁海霞, 崔春光, 赵天良, 等, 2019. 2015年冬季湖北省PM2.5重污染传输特征及影响天气系统的数值模拟[J]. 气象, 45(8): 1113-1122. Qi H X, Cui C G, Zhao T L, et al, 2019. Numerical simulation on the characteristics of PM2.5 heavy pollution and the influence of weather system in Hubei Province in winter 2015[J]. Meteor Mon, 45(8): 1113-1122 (in Chinese).
秦卓凡, 廖宏, 陈磊, 等, 2021. 汾渭平原空气质量及气象要素对其日变化和年际变化的影响[J]. 大气科学, 45(6): 1273-1291. Qin Z F, Liao H, Chen L, et al, 2021. Fenwei Plain air quality and the dominant meteorological parameters for its daily and interannual variations[J]. Chin J Atmos Sci, 45(6): 1273-1291 (in Chinese).
生态环境部, 2012. 环境空气质量标准(GB 3095—2012)[S]. 北京: 中国标准出版社. Ministry of Ecological Environment, 2012. Environmental air quality standard (GB 3095-2012)[S]. Beijing: Standards Press of China(in Chinese).
宋佳琨, 陈耀登, 陈丹, 2021. 气象-气溶胶资料联合同化对秋季PM2.5浓度模拟的影响研究[J]. 气象学报, 79(3): 477-491. Song J K, Chen Y D, Chen D, 2021. A study of meteorology-aerosol joint data assimilation on autumn PM2.5 concentration simulation[J]. Acta Meteor Sin, 79(3): 477-491 (in Chinese).
王芃, 朱盛强, 张梦媛, 等, 2022. 大气氧化性及其对二次污染物形成的贡献[J]. 科学通报, 67(18): 2069. Wang P, Zhu S Q, Zhang M Y, et al, 2022. Atmospheric oxidation capacity and its contribution tosecondary pollutants formation[J]. Chin Sci Bull, 67(18): 2069-2078 (in Chinese).
王妘涛, 张强, 温肖宇, 等, 2022. 运城市PM2.5时空分布特征和潜在源区季节分析[J]. 环境科学, 43(1): 74-84. Wang Y T, Zhang Q, Wen X Y, et al, 2022. Spatiotemporal distribution and seasonal characteristics of regional transport of PM2.5 in Yuncheng City[J]. Environ Sci, 43(1): 74-84 (in Chinese).
肖致美, 李源, 孔君, 等, 2022. 天津市PM2.5-O3复合污染特征及气象影响分析[J]. 环境科学, 43(6): 2928-2936. Xiao Z M, Li Y, Kong J, et al, 2022. Characteristics and meteorological factors of PM2.5-O3 compound pollution in Tianjin[J]. Environ Sci, 43(6): 2928-2936 (in Chinese).
解淑艳, 霍晓芹, 曾凡刚, 等, 2021. 2015—2019年汾渭平原臭氧污染状况分析[J]. 中国环境监测, 37(1): 49-57. Xie S Y, Huo X Q, Zeng F G, et al, 2021. Analysis of ozone pollution in Fenwei Plain from 2015 to 2019[J]. Environ Monit China, 37(1): 49-57 (in Chinese). DOI:10.19316/j.issn.1002-6002.2021.01.08
徐栋夫, 曹萍萍, 王源程, 2020. 成都一次重污染过程的气溶胶光学特性垂直分布[J]. 气象, 46(7): 948-958. Xu D F, Cao P P, Wang Y C, 2020. Study of the vertical distribution of aerosol optical properties during a heavy pollution event in Chengdu[J]. Meteor Mon, 46(7): 948-958 (in Chinese).
闫小利, 叶东, 牛广山, 等, 2021. 2017年冬焦作市一次持续严重污染天气过程分析[J]. 气象与环境科学, 44(4): 33-42. Yan X L, Ye D, Niu G S, et al, 2021. Analysis of a continuous heavy pollution weather process in Jiaozuo City in winter 2017[J]. Meteor Environ Sci, 44(4): 33-42 (in Chinese). DOI:10.16765/j.cnki.1673-7148.2021.04.005
杨乐超, 董雪丽, 徐波, 2018. 汾渭平原雾霾时空变化特征及其溢出效应[J]. 环境经济研究, 3(3): 75-87. Yang L C, Dong X L, Xu B, 2018. Spatial distribution and spillover effects of haze pollution in the Fen-Wei Plain[J]. J Environ Econ, 3(3): 75-87 (in Chinese). DOI:10.19511/j.cnki.jee.2018.03.006
杨颖川, 叶倩, 魏颖, 等, 2022. 内蒙古呼包鄂区域冬季大气细颗粒物污染成因解析[J]. 大气科学, 46(6): 1332-1348. Yang Y C, Ye Q, Wei Y, et al, 2022. Analysis of the causes of atmospheric fine particle pollution in winter in the Hohhot-Baotou-Ordos Area of Inner Mongolia[J]. Chin J Atmos Sci, 46(6): 1332-1348 (in Chinese).
杨镇江, 李柯, 廖宏, 等, 2023. 2022年夏季历史极端高温下我国近地表臭氧污染及气象成因分析[J]. 大气科学. Yang Z J, Li K, Liao H, et al, 2023. Analysis of surface ozone pollution in China amid the record summertime extreme heat of 2022[J]. Chin J Atmos Sci (in Chinese). DOI:10.3878/j.issn.1006-9895.2302.22211(inChinese)
张琴, 姚秀萍, 2021. 鲁中地区霾期间气溶胶的垂直分布及其与气象条件的关系[J]. 气象, 47(9): 1099-1112. Zhang Q, Yao X P, 2021. Vertical distribution of aerosol and its relationship with meteorological conditions during haze in central part of Shandong Province[J]. Meteor Mon, 47(9): 1099-1112 (in Chinese).
张天航, 徐冉, 桂海林, 2020. 2019年11月大气环流和天气分析[J]. 气象, 46(2): 283-288. Zhang T H, Xu R, Gui H L, 2020. Analysis of the November 2019 atmospheric circulation and weather[J]. Meteor Mon, 46(2): 283-288 (in Chinese).
张远远, 戴维, 华楠, 等, 2022. 基于不同周期PM2.5组成高时间分辨观测的PMF源解析研究[J]. 环境科学学报, 42(2): 308-317. Zhang Y Y, Dai W, Hua N, et al, 2022. PMF source apportionment based on high time-resolved measurements of PM2.5 components during different observation periods[J]. Acta Sci Circumstantiae, 42(2): 308-317 (in Chinese).
郑小华, 李明星, 娄盼星, 2021. 不同时间尺度下汾渭平原臭氧浓度变化及气象环境影响[J]. 高原气象, 40(4): 954-964. Zheng X H, Li M X, Lou P X, 2021. Different-scale changes in ozone concentration and meteorological environment in Fenwei Plain[J]. Plateau Meteor, 40(4): 954-964 (in Chinese).
中华人民共和国生态环境部, 2022.2021中国生态环境状况公报[R]. 北京: 中华人民共和国生态环境部. Ministry of Ecology and Environment of the People's Republic of China, 2022.2021 State of China's ecological environment bulletin[R]. Beijing: Ministry of Ecology and Environment of the People's Republic of China(in Chinese).
朱彤, 万薇, 刘俊, 等, 2022. 世界卫生组织《全球空气质量指南》修订解读[J]. 科学通报, 67(8): 697-706. Zhu T, Wan W, Liu J, et al, 2022. Insights into the new WHO global air quality guidelines[J]. Chin Sci Bull, 67(8): 697-706 (in Chinese).
Anenberg S C, Horowitz L W, Tong D Q, et al, 2010. An estimate of the global burden of anthropogenic ozone and fine particulate matter on premature human mortality using atmospheric modeling[J]. Environ Health Perspect, 118(9): 1189-1195. DOI:10.1289/ehp.0901220
Chen L, Zhu J, Liao H, et al, 2020. Meteorological influences on PM2.5 and O3 trends and associated health burden since China's clean air actions[J]. Sci Total Environ, 744: 140837. DOI:10.1016/j.scitotenv.2020.140837
Dai H B, Zhu J, Liao H, et al, 2021. Co-occurrence of ozone and PM2.5 pollution in the Yangtze River Delta over 2013-2019:spatiotemporal distribution and meteorological conditions[J]. Atmos Res, 249: 105363. DOI:10.1016/j.atmosres.2020.105363
Dao X, Ji D S, Zhang X, et al, 2021. Characteristics, sources and health risk assessment of PM2.5 in China's coal and coking heartland: insights gained from the regional observations during the heating season[J]. Atmos Pollutn Res, 12(12): 101237. DOI:10.1016/j.apr.2021.101237
Ding A J, Fu C B, Yang X Q, et al, 2013. Ozone and fine particle in the western Yangtze River Delta: an overview of 1 yr data at the SORPES station[J]. Atmos Chem Phys, 13(11): 5813-5830. DOI:10.5194/acp-13-5813-2013
Feng S L, Gao D, Liao F, et al, 2016. The health effects of ambient PM2.5 and potential mechanisms[J]. Ecotoxicol Environ Saf, 128: 67-74. DOI:10.1016/j.ecoenv.2016.01.030
Feng Y Y, Ning M, Lei Y, et al, 2019. Defending blue sky in China: effectiveness of the "air pollution prevention and control action plan" on air quality improvements from 2013 to 2017[J]. J Environ Manage, 252: 109603. DOI:10.1016/j.jenvman.2019.109603
He Y P, Li L, Wang H L, et al, 2022. A cold front induced co-occurrence of O3 and PM2.5 pollution in a Pearl River Delta city: temporal variation, vertical structure, and mechanism[J]. Environ Pollut, 306: 119464. DOI:10.1016/j.envpol.2022.119464
Huth R, Beck C, Philipp A, et al, 2008. Classifications of atmospheric circulation patterns: recent advances and applications[J]. Ann NY Acad Sci, 1146(1): 105-152. DOI:10.1196/annals.1446.019
Jiang Z J, Li J, Lu X, et al, 2021. Impact of western Pacific subtropical high on ozone pollution over eastern China[J]. Atmos Chem Phys, 21(4): 2601-2613. DOI:10.5194/acp-21-2601-2021
Li M M, Wang T J, Xie M, et al, 2018. Agricultural fire impacts on ozone photochemistry over the Yangtze River Delta region, East China[J]. J Geophys Res Atmos, 123(12): 6605-6623. DOI:10.1029/2018JD028582
Lu X, Zhang L, Chen Y F, et al, 2019. Exploring 2016-2017 surface ozone pollution over China: source contributions and meteorological influences[J]. Atmos Chem Phys, 19(12): 8339-8361. DOI:10.5194/acp-19-8339-2019
Luo Y H, Zhao T L, Yang Y J, et al, 2022. Seasonal changes in the recent decline of combined high PM2.5 and O3 pollution and associated chemical and meteorological drivers in the Beijing-Tianjin-Hebei region, China[J]. Sci Total Environ, 838: 156312. DOI:10.1016/j.scitotenv.2022.156312
Menon S, Unger N, Koch D, et al, 2008. Aerosol climate effects and air quality impacts from 1980 to 2030[J]. Environ Res Lett, 3(2): 024004. DOI:10.1088/1748-9326/3/2/024004
Miao Y C, Guo J P, Liu S H, et al, 2018. Impacts of synoptic condition and planetary boundary layer structure on the trans-boundary aerosol transport from Beijing-Tianjin-Hebei Region to Northeast China[J]. Atmos Environ, 181: 1-11.
Ning G C, Yim S H L, Wang S G, et al, 2019. Synergistic effects of synoptic weather patterns and topography on air quality: a case of the Sichuan Basin of China[J]. Climate Dyn, 53(11): 6729-6744.
Ou S J, Wei W, Cai B, et al, 2022. Exploring the causes for co-pollution of O3 and PM2.5 in summer over North China[J]. Environ Monit Assess, 194(4): 1-14.
Schnell J L, Prather M J, 2017. Co-occurrence of extremes in surface ozone, particulate matter, and temperature over eastern North America[J]. Proc Natl Acad Sci USA, 114(11): 2854-2859.
Tai A P K, Mickley L J, Jacob D J, 2010. Correlations between fine particulate matter(PM2.5) and meteorological variables in the United States: implications for the sensitivity of PM2.5 to climate change[J]. Atmos Environ, 44(32): 3976-3984.
Unger N, Menon S, Koch D M, et al, 2009. Impacts of aerosol-cloud interactions on past and future changes in tropospheric composition[J]. Atmos Chem Phys, 9(12): 4115-4129.
Wang D F, Zhou B, Fu Q Y, et al, 2016. Intense secondary aerosol formation due to strong atmospheric photochemical reactions in summer: observations at a rural site in eastern Yangtze River Delta of China[J]. Sci Total Environ, 571: 1454-1466.
Wang G, Zhu Z Y, Zhao N, et al, 2021. 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]. Environ Pollut, 284: 117450.
Wang X Y, Dickinson R E, Su L Y, et al, 2018. PM2.5 pollution in China and how it has been exacerbated by terrain and meteorological conditions[J]. Bull Amer Meteor Soc, 99(1): 105-119.
Xu Y L, Xue W B, Lei Y, et al, 2018. Impact of meteorological conditions on PM2.5 pollution in China during winter[J]. Atmosphere, 9(11): 429.
Yang K J, Kong L D, Tong S Y, et al, 2021. Double high-level ozone and PM2.5 co-pollution episodes in Shanghai, China: pollution characteristics and significant role of daytime HONO[J]. Atmosphere, 12(5): 557.
Ye X X, Song Y, Cai X H, et al, 2016. Study on the synoptic flow patterns and boundary layer process of the severe haze events over the North China Plain in January 2013[J]. Atmos Environ, 124: 129-145.
Zhang N N, Ma F, Qin C B, et al, 2018a. Spatiotemporal trends in PM2.5 levels from 2013 to 2017 and regional demarcations for joint prevention and control of atmospheric pollution in China[J]. Chemosphere, 210: 1176-1184.
Zhang Q Q, Ma Q, Zhao B, et al, 2018b. Winter haze over North China Plain from 2009 to 2016:influence of emission and meteorology[J]. Environ Pollut, 242: 1308-1318.
Zhang X X, Cheng C X, Zhao H, 2022. A health impact and economic loss assessment of O3 and PM2.5 exposure in China from 2015 to 2020[J]. GeoHealth, 6(3): e2021GH000531.
Zhang Y H, Hu M, Zhong L J, et al, 2008. Regional integrated experiments on air quality over Pearl River Delta 2004(PRIDE-PRD2004): overview[J]. Atmos Environ, 42(25): 6157-6173.