环境科学  2025, Vol. 46 Issue (2): 736-745   PDF    
成渝地区臭氧生成敏感性时空分布及影响因素
胥雯1, 杨皓1, 何敏1,2, 杨震1, 张洋1,2, 刘志红1,2, 何雨明3     
1. 成都信息工程大学资源环境学院, 成都 610225;
2. 成都平原城市气象与环境四川省野外科学观测研究站, 成都 610225;
3. 四川省环保科技工程有限责任公司, 成都 610041
摘要: 为研究成渝地区臭氧生成敏感性变化特征, 采用OMI卫星观测的对流层柱浓度FNR指示剂(HCHO/NO2), 分析夏季臭氧生成敏感性(OFS), 并基于地理探测器分析污染源排放、气象因子及其交互作用对OFS的影响. 结果表明:①成渝地区夏季OFS分区以协同控制区为主, 面积占比约为42.42%;VOCs控制区较少, 主要分布在成都市和重庆市主城区, 面积占比约为16.37%;NOx控制区主要分布在雅安市、乐山市、宜宾市、泸州市以及重庆市东部, 面积占比约为41.21%. ②从月份变化来看, 6~8月, VOCs控制区面积逐渐减少, 8月达到最低值, 占比仅有3.53%, 9月VOCs控制区面积增加显著. ③大气污染物(CO、PM2.5、NMVOC和NOx)排放是影响成渝地区夏季OFS的主要因素, 随着排放量的增加, FNR逐渐减小, OFS趋向于VOCs控制区. ④地表大气压力(SP)、2 m气温(T2), 大气柱总水量(TCW)是影响成渝地区夏季OFS变化的3个主要气象因素, 三者均与FNR呈负相关, SP、T2, TCW的值越大, OFS更趋向于VOCs控制. ⑤因子间的交互作用会加剧臭氧生成敏感性的变化, 污染物和气象因子之间的交互作用是成渝地区夏季OFS变化的主要驱动力.
关键词: 臭氧检测仪(OMI)      臭氧生成敏感性(OFS)      时空分布      影响因素      成渝地区     
Spatial and Temporal Characteristics of Ozone Formation Sensitivity and Its Influencing Factors in Chengdu-Chongqing Area
XU Wen1 , YANG Hao1 , HE Min1,2 , YANG Zhen1 , ZHANG Yang1,2 , LIU Zhi-hong1,2 , HE Yu-ming3     
1. College of Resources and Environment, Chengdu University of Information Technology, Chengdu 610225, China;
2. Chengdu Plain Urban Meteorology and Environment Observation and Research Station of Sichuan Province, Chengdu University of Information Technology, Chengdu 610225, China;
3. Sichuan Environmental Protection Technology Engineering Co., Ltd., Chengdu 610041, China
Abstract: To study the characteristics of ozone formation sensitivity changes in the Chengdu-Chongqing area, the summer ozone formation sensitivity (OFS) was analyzed using satellite observation of the OMI tropospheric column concentration of HCHO and NO2. Additionally, the effects of meteorological factors, source emissions, and their interactions on OFS were identified with the geodetic detector. The results showed that the summer OFS zoning in the Chengdu-Chongqing area was dominated by a transitional regime, with an area share of 42.42%. Less VOCs-limited regime was observed, which was mainly distributed in the urban areas of Chengdu City and Chongqing Municipality, with an area share of 16.37%. The NOx-limited regime was mainly located in Ya'an City, Leshan City, Yibin City, Luzhou City, and the eastern part of Chongqing Municipality, with an area share of 41.21%. The area of the VOCs-limited regime gradually decreased from June to August and reached the lowest value in August, accounting for only 3.53%, and the area of the VOCs-limited regime increased significantly in September. Pollutant emissions (CO, PM2.5, NMVOC, and NOx) were the main factors affecting summer OFS in the Chengdu-Chongqing area. With the increase in emissions, the FNR gradually decreased, and the OFS tended to be the VOCs-limited regime. SP, T2, and TCW were the three main meteorological factors that affected the change in summer OFS in the Chengdu-Chongqing area. All three factors were negatively correlated with FNR; the larger the values of SP, T2, and TCW, the more easily OFS changed to a VOCs-limited regime. Interactions between factors exacerbated OFS changes and were the main drivers of spatial and temporal variations for summer OFS in the Chengdu-Chongqing area.
Key words: ozone monitoring instrument (OMI)      ozone formation sensitivity (OFS)      spatial and temporal characteristics      influencing factors      Chengdu-Chongqing area     

挥发性有机物(VOCs)和氮氧化物(NOx)等前体物在太阳辐射下发生光化学反应[1~4], 导致近地面臭氧浓度超标的现象称为臭氧污染, 这已成为当今社会面临的严重环境问题之一, 高浓度臭氧可对人体健康、气候变化和农作物生长等造成严重危害[5~8]. 2022年, 全国339个地级及以上的城市中, 有126个城市的空气质量未达到二级标准(占比37.2%), 其中92个城市臭氧浓度超标[9].

臭氧污染不仅与VOCs和NOx等前体物有着密切关联[10], 还受到气温[11, 12]、辐射强度[13]、湿度和风速[14]等多种气象因素的影响[6, 15, 16], 其浓度水平对前体物排放变化有显著的非线性响应特征[17~19]. 臭氧生成敏感性(ozone formation sensitivity, OFS)的研究方法包括源示踪法、敏感性测试法和指示剂法[20]. 源示踪法利用气象数据、空气质量数据以及数学模型来监测和分析大气中的不同化学物质的浓度, 进而确定不同污染源的贡献. 敏感性测试法利用空气质量模型模拟臭氧浓度在各类排放场景下的变化, 以评估不同前体物对臭氧敏感性的影响[21]. 指示剂法使用特定的化学指标或指示剂来测量O3生成的敏感性, 常用的指示剂包括NOx、O3/NOy、HCHO/NO2、过氧化氢产率/硝酸产率[P(H2O2)/P(HNO3)]和HCHO/NOy[22~24]. 其中, 采用卫星观测的柱浓度产品FNR指示剂[甲醛(HCHO)/二氧化氮(NO2), ratio of the tropospheric columns of formaldehyde to nitrogen dioxide]判断臭氧控制区[20, 25~29]的方法具有时空连续以及受人为影响较小等优势得到了广泛应用[20]. Wang [30]等使用FNR分析了2016~2019年中国上空臭氧形成敏感性的时空变异性, 结果表明, 在中国的城市地区, O3生成敏感性倾向于受挥发性有机物控制, 而在农村和偏远地区, 则倾向于受氮氧化物控制. 李毅[31]采用了WRF-CMAQ模型, 并与OMI卫星遥感数据中的HCHO/NO2(柱浓度比值)相结合, 从而计算出本地化的FNR指标阈值, 并进一步分析了珠三角地区在臭氧生成敏感性方面的差异性. 白杨[27]等采用卫星观测OMI FNR指示剂分析河南省夏季臭氧生成敏感性以及驱动因素.

成渝地区是中国西部的核心经济区, 包括重庆和成都两大城市以及其周边地区. 近年来, 伴随着成渝地区经济的加速发展, 工业活动、交通排放和城市扩张等人类活动使得该区域内臭氧污染形势愈发严峻, 已成为我国除京津冀、长三角和珠三角地区之外的另一污染高发区[32, 33]. 掌握成渝地区大气臭氧的生成敏感性对于臭氧污染预测和管理至关重要, 而目前还鲜见成渝地区臭氧敏感性时空变化及影响因素的研究. 本研究利用OMI卫星观测的对流层NO2和HCHO柱浓度, 采用FNR指示剂研究成渝地区O3生成敏感性的时空分布, 并探讨污染源排放和气象条件等因素的影响, 以更好地理解区域臭氧生成的机制, 以期为臭氧污染防控提供科学依据.

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

研究区域为成渝地区(见图 1), 包括成都平原地区(成都市、德阳市、绵阳市、眉山市、资阳市、乐山市、雅安市和遂宁市)、川南地区(内江市、自贡市、宜宾市和泸州市)和川东北地区(广安市、南充市和达州市)以及重庆地区.

图 1 2016~2021年成渝地区臭氧污染概况 Fig. 1 Ozone pollution in Chengdu-Chongqing area from 2016 to 2021

采用国控站点监测数据, 对2016~2021年成渝地区各城市第90百分位日最大8 h O3浓度进行统计(见图 1), 结果显示, O3浓度高值主要集中在四川的成都平原地区以及重庆的西南部, 其中沙坪坝区、江津区以及璧山区ρ(O3)超过160 μg·m-3.

1.2 数据来源

NO2对流层柱浓度数据来自NASA(https://disc.gsfc.nasa.gov)发布的OMI对流层二氧化氮柱浓度L3日均值级产品(OMNO2d), 使用IDL(interactive data language)语言对其进行处理, 筛选出云覆盖比例小于30%且地表反射率小于0.3的像元, 数据产品的空间分辨率为0.25°×0.25°.

HCHO对流层柱浓度数据来自于NASA发布的OMI对流层甲醛垂直柱浓度L3级产品(OMHCHOd), 使用IDL语言对该数据进行处理, 排除负值、云量大于30%以及误差较大的数据, 数据产品的空间分辨率为0.1°×0.1°.

污染物排放数据来源于中国多尺度排放清单模型MEIC(http://www.meicmodel.org/), MEIC清单由清华大学开发, 被认为是最具代表性的区域排放源清单, 空间分辨率为0.25°×0.25°, 包含电力、工业、交通、民用和农业这5个排放部门[34]. 本研究使用逐月的NOx、NMVOC、CO、PM2.5产品网格化数据研究OFS与人为源污染物排放的关系.

气象数据来自欧洲中期天气预报中心ECMWF[35]发布的月均尺度的气象再分析产品ERA5, 选择UTC时间为05:00和06:00(北京时间13:00和14:00)的均值代表OMI过境时刻. 主要考虑的气象要素包括边界层高度(BLH)、10 m风速(WS10)、地表大气压力(SP)、地表净太阳辐射(SSR)、2 m气温(T2)以及大气柱总水量(TCW). 气象数据的空间分辨率为0.25°×0.25°.

1.3 研究方法 1.3.1 臭氧敏感性分区方案

提取OMI卫星观测的成渝地区2016~2021年6~9月的网格化HCHO和NO2柱浓度数据, 为了实现不同数据间的空间分辨率匹配, 对分辨率为0.25°的NO2数据重采样至0.1°. 随后, 对研究区域内各网格的FNR(HCHO/NO2)值进行计算, 其公式如下:

(1)

式中, c(HCHO)为HCHO柱浓度, c(NO2)为NO2柱浓度, FNR为二者比值. 采用Wang等[30]的分区方案进行臭氧敏感性研究, 如表 1所示.

表 1 臭氧敏感性分区方案 Table 1 Ozone sensitive zoning program

1.3.2 地理探测器

采用地理探测器[36]解析成渝地区夏季OFS的空间分异性和影响因素. 其中, 分异及因子探测器能够揭示因变量在空间上的差异, 并评估因子在多大程度上解释了因变量的空间分异, 数学模型如下式所示:

(2)
(3)
(4)

式中, h=1, 2, 3, …, L为O3浓度或驱动因子的分层, 即分类或分区;NhN分别为层h和川渝地区的地级市数;σh2σ2分别为层h和选取驱动因子的方差. q统计量用来描述因子X对OFS空间异质性的解释程度. 因子X涵盖了6种气象因子(BLH、WS10、SP、SSR、T2和TCW)和4种污染物(NOx、NMVOC、CO和PM2.5)的排放量. q值的值域为[0, 1], 值越大说明驱动因子对OFS变化的影响越显著.

交互作用探测器用于检测不同因子之间的相互作用效应, 评估这些因子共同作用时是否会改变对因变量的解释力. 首先分别探测两个因子对臭氧生成敏感性的解释能力qX1)和qX2), 再计算因子交互时的q值:qX1X2), 并对qX1)、qX2)和qX1X2)进行比较. 两个因子交互作用类型可分为以下5类(见表 2).

表 2 两个单因子交互作用的类型 Table 2 Types of two single-factor interactions

2 结果与讨论 2.1 O3前体物时空分布

在阳光照射下, NOx和VOCs经过一系列的光化学链式反应后生成臭氧[18, 37]. NO2是NOx的关键组分[38], HCHO由VOCs氧化生成, 在大气中的氧化和光解寿命很短[31, 39], 且对自由基化学过程具有显著影响[40], 因此可通过分析NO2和HCHO柱浓度, 探讨O3前体物的时空分布特征.

图 2(a)所示, 各区域HCHO柱浓度在一定范围内波动, 16市的HCHO柱浓度均值在2018年达到峰值, 为197×1014 molec·cm-2. 2021年, 川南地区浓度最高, 均值为199×1014 molec·cm-2, 其次是重庆市和成都平原, 均值分别为189×1014 molec·cm-2和187×1014 molec·cm-2, 川东北浓度最低, 均值为179×1014 molec·cm-2. 根据图 3(a)所示, HCHO柱浓度高值地区呈点状分布, 主要位于成都市中心城区, 绵阳市和德阳市西边的部分地区以及重庆市的部分区域. 与图 2(a)相对应, 川东北地区HCHO柱浓度普遍较低.

图 2 2016~2021年的6~9月HCHO和NO2柱浓度 Fig. 2 Column concentrations of HCHO and NO2 for June to September from 2016 to 2021

图 3 2021年的6~9月HCHO和NO2柱浓度 Fig. 3 Column concentrations of HCHO and NO2 from June to September in 2021

图 2(b)所示, 各区域NO2柱浓度变化趋势基本一致, 16市的NO2柱浓度平均值变化平稳, 基本保持在48×1014 molec·cm-2左右. 2021年, 川东北地区的NO2柱浓度最高, 其次是重庆市和成都平原, 川南地区最低. 根据图 3(b)可以看出, 2021年的NO2柱浓度呈现片状分布, 以成都市、德阳市以及重庆市主城区为中心, 向外逐渐降低. NO2柱浓度高值区域主要分布于成都平原、重庆市主城区、主城以西、主城以东和川东北部分城市, 主要是一些煤源型大型工业企业分布较多和人口、道路密集的城市.

2.2 O3-NOx-VOCs敏感性的空间变化特征 2.2.1 O3控制区年际变化

基于臭氧生成敏感区划分方案, 成渝地区不同年份6~9月O3生成敏感性分类结果和O3控制区面积占比变化分别如图 4图 5所示. 可以看出, 成渝地区臭氧生成以协同控制区为主, 主要分布在成都市、绵阳市、德阳市、广安市、遂宁市和南充市等成渝地区中部城市. 2017年后, 随着NOx控制区的增加, 协同控制区范围呈下降趋势, 在2020年后有所上升, 2021年占比最高(约为46.13%). VOCs控制区主要分布在成都市和重庆市主城区, 绵阳市、德阳市和南充市等部分城市有少量分布. 2016~2021年间, VOCs控制区面积占比呈“M”形变化趋势, 其中2017年的VOCs控制区比例最高(约为21.66%). 2019年后, 呈现出由VOCs控制区转变为协同控制区的趋势. NOx控制区主要分布在成渝区域边缘, 如雅安市和乐山市. 2018年NOx控制区面积最大(约占49.12%), 除雅安市、自贡市和乐山市等地区之外, 眉山市、成都市、内江市、资阳市以及重庆东部也出现较大范围的NOx控制区, 相比2017年增长了45%.

图 4 不同年份的6~9月O3生成敏感性分类结果 Fig. 4 Classification of ozone formation sensitivity for different years from June to September

图 5 不同年份的6~9月O3控制区占比变化 Fig. 5 Changes in the percentage of O3 control area for different years from June to September

在臭氧生成敏感性分区结果的基础上, 基于MEIC清单, 进一步分析不同控制区内VOCs和NOx的排放特征. 从图 6可以看出, 各敏感性分区中NOx和VOCs排放来源较为一致, 其中NOx主要来源于工业源及交通源, VOCs主要来源于居民源及工业源. 不同敏感性分区的源排放构成相比而言, VOCs控制区中工业NOx和VOCs排放贡献占比略高于协同控制区及NOx控制区, 而NOx控制区中居民污染物排放贡献则高于其它两个区域, 表明不同污染源排放对臭氧生成敏感性有重要影响.

臭氧控制区类型:1.VOCs控制区, 2.协同控制区, 3.NOx控制区 图 6 不同臭氧控制区内人为源NOx和NMVOC排放特征 Fig. 6 NOxand NMVOC emissions from anthropogenic sources in different ozone control areas

2.2.2 O3控制区月际变化

图 7图 8给出了2016~2021年不同月份O3敏感性空间分布及面积占比变化, 可以看出, 6~8月期间VOCs控制区逐渐减少, 8月达到最低值, 占比仅有4.02%. 9月VOCs控制区呈明显上升趋势, 可能受天然源VOCs排放减少的影响[20, 41]. 图 9显示, 9月HCHO柱浓度与6~8月相比显著下降, 主要受到植被挥发性有机化合物排放变化的影响[20]. 毛红梅等[42]研究结果表明, 在夏季, 由于高温、强烈的太阳辐射和长时间的日照等影响因素, 植被挥发性有机化合物的排放在6~8月达到峰值, 进入9月, 随着气温的下降, 排放量也随之减少.

图 7 2016~2021年的不同月份O3生成敏感性分类结果 Fig. 7 Classification of ozone formation sensitivity for different months from 2016 to 2021

图 8 2016~2021年不同月份O3控制区占比变化 Fig. 8 Changes in the percentage of O3 control area for different months from 2016 to 2021

图 9 各城市2016~2021年的不同月份HCHO柱浓度 Fig. 9 Column concentrations of HCHO for different months from 2016 to 2021 in each city

2.3 O3生成敏感性变化的影响因素分析

以2020年为例, 基于MEIC清单、ERA5气象数据以及各网格的FNR, 采用地理探测器分析4个污染物排放因子(CO、NOx、PM2.5和NMVOC)、6个气象因子(SSR、T2、TCW、WS10、BLH和SP)以及它们之间的交互作用对成渝地区夏季臭氧生成敏感性变化的影响, 并研究了污染物排放、气象要素与FNR之间的相关性.

图 10所示, 污染物排放是成渝地区夏季臭氧生成敏感性的主要影响因素;同时, 通过对图 11(a)~11(d)的观察, 随着4种污染物排放的增加, FNR逐渐减小, 臭氧生成敏感性趋向于VOCs控制. 相比于NMVOC(q=0.22)和NOxq=0.21)这2种驱动因子, CO(q=0.33)对臭氧生成敏感性变化的解释力更强, 这可能是因为CO的排放量显著高于VOCs和NOx, 更容易与·OH和O2发生反应, 生成HO2·, 加速光化学反应速率, 从而促使臭氧生成敏感性更易改变[27]. PM2.5的增加(q=0.30)导致臭氧生成敏感区趋向于VOCs控制, 这是因为PM2.5的化学反应涉及HOx和NOx·的活性吸收, 其化学作用导致臭氧的产生比其他情况下更受VOCs的限制, 即对于NOx减排的敏感性降低[43].

X1X2分别表示第一个和第二个影响因子, X1X2表示两个影响因子之间的交互作用:1.CO∩NOx, 2.CO∩PM2.5, 3.CO∩NMVOC, 4.CO∩SSR, 5.CO∩T2, 6.CO∩TCW, 7.CO∩WS10, 8.CO∩BLH, 9.CO∩SP, 10.NOx∩PM2.5, 11.NOx∩NMVOC, 12.NOx∩SSR, 13.NOx∩T2, 14. NOx∩TCW, 15. NOx∩WS10, 16.NOx∩BLH, 17.NOx∩SP, 18.PM2.5∩NMVOC, 19.PM2.5∩SSR, 20.PM2.5∩T2, 21.PM2.5∩TCW, 22.PM2.5∩WS10, 23.PM2.5∩BLH, 24.PM2.5∩SP, 25.NMVOC∩SSR, 26.NMVOC∩T2, 27.NMVOC∩TCW, 28.NMVOC∩WS10, 29.NMVOC∩BLH, 30.NMVOC∩SP, 31.SSR∩T2, 32.SSR∩TCW, 33.SSR∩WS10, 34.SSR∩BLH, 35.SSR∩SP, 36.T2∩TCW, 37.T2∩WS10, 38.T2∩BLH, 39.T2∩SP, 40.TCW∩WS10, 41.TCW∩BLH, 42.TCW∩SP, 43.WS10∩BLH, 44.WS10∩SP, 45.BLH∩SP 图 10 人为源排放和气象因子的q Fig. 10 The q values of anthropogenic source emissions and meteorological factors

图 11 成渝地区FNR与各影响因子间的相关关系 Fig. 11 Correlation relationship between FNR and each influencing factor in Chengdu-Chongqing area

SP(q=0.268)、T2(q=0.271)和TCW(q=0.259)这3个因子的q值均大于0.25, 是影响臭氧生成敏感性的主导气象因子. BLH(q=0.098)、WS10(q=0.13)、SSR(q=0.08)、对成渝地区夏季臭氧生成敏感性的变化影响较小[图 11(e)11(f)11(h)], 且与FNR的相关性不显著. FNR与SP呈负相关, 随着SP的增加, 臭氧对VOCs更加敏感. 如图 11(g), 当SP < 90 000 Pa时, 臭氧生成敏感性趋向于NOx控制;当SP > 95 000 Pa时, 臭氧生成敏感性趋向于VOCs控制. 这可能是因为高气压通常与晴朗的天气相联系, 日照更强, 从而促进NO2的光化学产生, 使得FNR减小[27]. FNR与T2呈显著负相关, T2 > 302 K时, 臭氧生成敏感性趋近于VOCs控制;T2 < 300 K时, 臭氧生成敏感性趋向于NOx控制. 主要因为温度升高, VOCs的氧化反应速率增加, 产生的氧自由基(RO2·和HO2·)将NO氧化为NO2[44], 从而使得FNR减小, 臭氧生成敏感性趋向于VOCs控制. TCW对成渝地区臭氧生成敏感性变化也有一定影响, 随着TCW的增加, FNR减小, 当TCW < 45 kg·m-2时, 臭氧生成敏感性趋向于NOx控制;当TCW > 55 kg·m-2时, 臭氧生成敏感性趋向于VOCs控制. 研究表明在污染较重的城市地区, 水汽的增加会促进NO2转换为硝酸[45], NO2与VOCs生成臭氧的反应受到抑制, 更趋向于VOCs控制[图 11(j)].

从交互作用探测结果来看[图 10(b)], 两个因子间的交互作用对臭氧生成敏感性空间变化的影响大于单一因子的独立因子, 这是因为气象条件的改变会影响大气中的光解反应, 进而改变前体物排放生成臭氧的速率, 对成渝地区夏季臭氧生成的敏感性产生影响[32]. CO、PM2.5与气象因素之间的交互作用以及各气象因素之间的交互作用是成渝地区夏季臭氧生成敏感性的主要影响因素. 例如, TCW与WS10的交互作用导致臭氧生成敏感性易发生改变, 在低风速低湿度的环境下, FNR较大, 臭氧生成敏感性趋向于NOx控制;而在高风速高湿度的环境下, FNR减小, 臭氧生成敏感性趋向于VOCs控制.

3 结论

(1)2016~2021年臭氧生成敏感性的年际变化特征表明, 成渝地区OFS分区以协同控制区为主, 主要分布在成都市、绵阳市、德阳市、广安市、遂宁市和南充市等成渝地区的中部城市, 2021年占比最高(约为46.13%);VOCs控制区较少, 主要分布在成都市和重庆市主城区, 其面积占比呈现“M”型变化趋势;NOx控制区主要分布在川渝边缘(如雅安市和乐山市), 2018年NOx控制区(约占49.12%)最高, 相比前一年增长约45%.

(2)臭氧生成敏感性的月际变化特征表明, 6~8月期间VOCs控制区逐渐减少, 8月达到最低值, 占比仅有4.02%. 9月由于温度降低, 降水量减少等原因, 天然源挥发性有机化合物排放量显著降低, 导致NOx控制区向协同控制区转化, 协同控制区向VOCs控制区转化, VOCs控制区增加显著.

(3)污染物排放(CO、PM2.5、NMVOC和NOx)是成渝地区夏季OFS的主导因素;随着这4种污染物排放的增加, FNR值逐渐降低, OFS呈现向VOCs控制区变化的趋势, 对NOx减排的敏感性减弱. 同时, SP、T2, TCW这3个气象因素在夏季臭氧生成敏感性的变化中扮演关键角色, 三者与FNR呈现负相关关系. 各因子之间的交互作用加剧了臭氧生成敏感性的变化, 污染物与气象因子之间的交互作用是成渝地区夏季臭氧生成敏感性变化的主要驱动力.

参考文献
[1] Wang T, Xue L K, Brimblecombe P, et al. Ozone pollution in China: a review of concentrations, meteorological influences, chemical precursors, and effects[J]. Science of the Total Environment, 2017, 575: 1582-1596. DOI:10.1016/j.scitotenv.2016.10.081
[2] Lelieveld J, Evans J S, Fnais M, et al. The contribution of outdoor air pollution sources to premature mortality on a global scale[J]. Nature, 2015, 525(7569): 367-371. DOI:10.1038/nature15371
[3] Zhang K, Li L, Huang L, et al. The impact of volatile organic compounds on ozone formation in the suburban area of Shanghai[J]. Atmospheric Environment, 2020, 232. DOI:10.1016/j.atmosenv.2020.117511
[4] Zhang Y J, Gao J, Zhu Y J, et al. Evolution of ozone formation sensitivity during a persistent regional ozone episode in northeastern China and its implication for a control strategy[J]. Environmental Science & Technology, 2023, 58(1): 617-627.
[5] 程育恺, 戴海夏, 张蕴晖, 等. 长三角地区2017-2020年臭氧浓度时空分布与人群健康效益评估[J]. 环境科学, 2023, 44(2): 719-729.
Cheng Y K, Dai H X, Zhang Y H, et al. Spatial and temporal distribution characteristics of ozone concentration and population health benefit assessment in the Yangtze River Delta Region from 2017 to 2020[J]. Environmental Science, 2023, 44(2): 719-729.
[6] Xie Y, Dai H C, Zhang Y X, et al. Comparison of health and economic impacts of PM2.5 and ozone pollution in China[J]. Environment International, 2019, 130. DOI:10.1016/j.envint.2019.05.075
[7] Wang F Y, Qiu X H, Cao J Y, et al. Policy-driven changes in the health risk of PM2.5 and O3 exposure in China during 2013-2018[J]. Science of the Total Environment, 2021, 757. DOI:10.1016/j.scitotenv.2020.143775
[8] 周骥, 付世华, 彭丽, 等. 臭氧和PM2.5对慢阻肺死亡影响及气温修饰效应[J]. 中国环境科学, 2021, 41(12): 5904-5911.
Zhou J, Fu S H, Peng L, et al. Effect of ozone and PM2.5 on chronic obstructive pulmonary disease deaths and the effect modification by the air temperature[J]. China Environmental Science, 2021, 41(12): 5904-5911.
[9] 中华人民共和国生态环境部. 中国生态环境状况公报[EB/OL]. http://www.cnemc.cn/jcbg/zghjzkgb/202305/W020230531511478063909.pdf, 2023-05-31.
[10] Yu R L, Lin Y L, Zou J H, et al. Review on atmospheric ozone pollution in China: formation, spatiotemporal distribution, precursors and affecting factors[J]. Atmosphere, 2021, 12(12). DOI:10.3390/atmos12121675
[11] Pusede S E, Gentner D R, Wooldridge P J, et al. On the temperature dependence of organic reactivity, nitrogen oxides, ozone production, and the impact of emission controls in San Joaquin Valley, California[J]. Atmospheric Chemistry and Physics, 2014, 14(7): 3373-3395. DOI:10.5194/acp-14-3373-2014
[12] 肖致美, 李鹏, 孔君, 等. 天津市持续高温强光照天气下臭氧污染差异性[J]. 中国环境科学, 2023, 43(7): 3322-3330.
Xiao Z M, Li P, Kong J, et al. Difference of ozone pollution under the continuous high temperature and strong sunlight weather in Tianjin[J]. China Environmental Science, 2023, 43(7): 3322-3330.
[13] Fix M J, Cooley D, Hodzic A, et al. Observed and predicted sensitivities of extreme surface ozone to meteorological drivers in three US cities[J]. Atmospheric Environment, 2018, 176: 292-300. DOI:10.1016/j.atmosenv.2017.12.036
[14] 赵伟, 高博, 刘明, 等. 气象因素对香港地区臭氧污染的影响[J]. 环境科学, 2019, 40(1): 55-66.
Zhao W, Gao B, Liu M, et al. Impact of meteorological factors on the ozone pollution in Hong Kong[J]. Environmental Science, 2019, 40(1): 55-66.
[15] 李飞, 孔少杰, 屈志光, 等. 中国臭氧时空分异及热点城市群气象关联特征[J]. 中国环境科学, 2023, 43(4): 1539-1549.
Li F, Kong S J, Qu Z G, et al. Spatiotemporal differentiation of ozone pollution in Chinese cities and meteorological correlation among the identified hot urban agglomeration during 2018 to 2020[J]. China Environmental Science, 2023, 43(4): 1539-1549.
[16] 窦筱艳, 高海鹏, 陈珂, 等. 青海东北部一次臭氧污染过程机制研究[J]. 中国环境科学, 2023, 43(S1): 46-54.
Dou X Y, Gao H P, Chen K, et al. Study on the formation mechanism of an ozone pollution event in Northeast Qinghai[J]. China Environmental Science, 2023, 43(S1): 46-54.
[17] 李光耀, 陈强, 郭文凯, 等. 基于正交试验的臭氧及其前体物的非线性响应及控制方案[J]. 环境科学, 2021, 42(2): 616-623.
Li G Y, Chen Q, Guo W K, et al. Nonlinear response characteristics and control scheme for ozone and its precursors based on orthogonal experimental methods[J]. Environmental Science, 2021, 42(2): 616-623.
[18] Liu Y M, Wang T. Worsening urban ozone pollution in China from 2013 to 2017–Part 2: the effects of emission changes and implications for multi-pollutant control[J]. Atmospheric Chemistry and Physics, 2020, 20(11): 6323-6337.
[19] Soleimanian E, Wang Y X, Li W, et al. Understanding ozone episodes during the TRACER-AQ campaign in Houston, Texas: the role of transport and ozone production sensitivity to precursors[J]. Science of the Total Environment, 2023, 900. DOI:10.1016/j.scitotenv.2023.165881
[20] 武卫玲, 薛文博, 雷宇, 等. 基于OMI数据的京津冀及周边地区O3生成敏感性[J]. 中国环境科学, 2018, 38(4): 1201-1208.
Wu W L, Xue W B, Lei Y, et al. Sensitivity analysis of ozone in Beijing-Tianjin-Hebei (BTH) and its surrounding area using OMI satellite remote sensing data[J]. China Environmental Science, 2018, 38(4): 1201-1208.
[21] 朱家贤, 王晓琦, 欧盛菊, 等. 2019年7月石家庄市O3生成敏感性及控制策略解析[J]. 环境科学, 2022, 43(7): 3473-3482.
Zhu J X, Wang X Q, Ou S J, et al. Ozone sensitivity analysis and control strategy in Shijiazhuang City in July 2019[J]. Environmental Science, 2022, 43(7): 3473-3482.
[22] Tonnesen G S, Dennis R L. Analysis of radical propagation efficiency to assess ozone sensitivity to hydrocarbons and NOx: 2. Long‐lived species as indicators of ozone concentration sensitivity[J]. Journal of Geophysical Research: Atmospheres, 2000, 105(D7): 9227-9241.
[23] Sillman S, He D Y. Some theoretical results concerning O3‐NOx‐VOC chemistry and NOx‐VOC indicators[J]. Journal of Geophysical Research: Atmospheres, 2002, 107(D22). DOI:10.1029/2001JD001123
[24] Du X H, Tang W, Cheng M M, et al. Modeling of spatial and temporal variations of ozone-NOx-VOC sensitivity based on photochemical indicators in China[J]. Journal of Environmental Sciences, 2022, 114: 454-464.
[25] 赵伟, 高博, 卢清, 等. 2006~2019年珠三角地区臭氧污染趋势[J]. 环境科学, 2021, 42(1): 97-105.
Zhao W, Gao B, Lu Q, et al. Ozone pollution trend in the Pearl River Delta Region during 2006-2019[J]. Environmental Science, 2021, 42(1): 97-105.
[26] Jin X M, Fiore A, Boersma K F, et al. Inferring changes in summertime surface ozone-NOx-VOC chemistry over U.S. urban areas from two decades of satellite and ground-based observations[J]. Environmental Science & Technology, 2020, 54(11): 6518-6529.
[27] 白杨, 王盼, 赵鹏飞, 等. 河南省夏季臭氧生成敏感性及其驱动因素分析[J]. 遥感学报, 2022, 26(5): 988-1001.
Bai Y, Wang P, Zhao P F, et al. Summertime ozone formation sensitivity and driving factors in Henan Province[J]. National Remote Sensing Bulletin, 2022, 26(5): 988-1001.
[28] Chu W H, Li H, Ji Y Y, et al. Research on ozone formation sensitivity based on observational methods: development history, methodology, and application and prospects in China[J]. Journal of Environmental Sciences, 2024, 138: 543-560.
[29] Song H Y, Zhao W J, Yang X C, et al. Ozone sensitivity analysis and ozone formation regimes division in the Beijing–Tianjin–Hebei Region Based on satellite remote sensing data[J]. Atmosphere, 2023, 14(11). DOI:10.3390/atmos14111637
[30] Wang W N, van der A R, Ding J Y, et al. Spatial and temporal changes of the ozone sensitivity in China based on satellite and ground-based observations[J]. Atmospheric Chemistry and Physics, 2021, 21(9): 7253-7269.
[31] 李毅. 基于卫星遥感与空气质量模型的珠江三角洲臭氧生成敏感性研究[D]. 广州: 华南理工大学, 2021.
Li Y. Study of ozone formation sensitivity in the Pearl River Delta based on satellite remote sensing and air quality model[D]. Guangzhou: South China University of Technology, 2021.
[32] Zhao S P, Yu Y, Qin D H, et al. Analyses of regional pollution and transportation of PM2.5 and ozone in the city clusters of Sichuan Basin, China[J]. Atmospheric Pollution Research, 2019, 10(2): 374-385.
[33] Zhang X Y, Wang Y Q, Niu T, et al. Atmospheric aerosol compositions in China: spatial/temporal variability, chemical signature, regional haze distribution and comparisons with global aerosols[J]. Atmospheric Chemistry and Physics, 2012, 12(2): 779-799.
[34] Li M, Liu H, Geng G N, et al. Anthropogenic emission inventories in China: a review[J]. National Science Review, 2017, 4(6): 834-866.
[35] Owens R. ECMWF forecast user guide[EB/OL]. https://www.ecmwf.int/en/elibrary/81307-ecmwf-forecast-user-guide, 2018-05-01.
[36] 王劲峰, 徐成东. 地理探测器: 原理与展望[J]. 地理学报, 2017, 72(1): 116-134.
Wang J F, Xu C D. Geodetector: principle and prospective[J]. Acta Geographica Sinica, 2017, 72(1): 116-134.
[37] Martin R V, Fiore A M, Van Donkelaar A. Space-based diagnosis of surface ozone sensitivity to anthropogenic emissions[J]. Geophysical Research Letters, 2004, 31(6). DOI:10.1029/2004GL019416
[38] 高晋徽, 朱彬, 王言哲, 等. 2005~2013年中国地区对流层二氧化氮分布及变化趋势[J]. 中国环境科学, 2015, 35(8): 2307-2318.
Gao J H, Zhu B, Wang Y Z, et al. Distribution and long-term variation of tropospheric NO2 over China during 2005 to 2013[J]. China Environmental Science, 2015, 35(8): 2307-2318.
[39] Liu J H, Strode S A, Liang Q, et al. Change in tropospheric ozone in the recent decades and its contribution to global total ozone[J]. Journal of Geophysical Research: Atmospheres, 2022, 127(22). DOI:10.1029/2022JD037170
[40] Palmer P I, Jacob D J, Fiore A M, et al. Mapping isoprene emissions over North America using formaldehyde column observations from space[J]. Journal of Geophysical Research: Atmospheres, 2003, 108(D6). DOI:10.1029/2002JD002153
[41] Gong C, Liao H, Yue X, et al. Impacts of ozone-vegetation interactions on ozone pollution episodes in north China and the Yangtze River Delta[J]. Geophysical Research Letters, 2021, 48(12). DOI:10.1029/2021GL093814
[42] 毛红梅, 张凯山, 第宝锋. 四川省天然源VOCs排放量的估算和时空分布[J]. 中国环境科学, 2016, 36(5): 1289-1296.
Mao H M, Zhang K S, Di B F. Studies on estimates of biogenic VOC emission and its temporal and spatial distribution in Sichuan[J]. China Environmental Science, 2016, 36(5): 1289-1296.
[43] Li K, Jacob D J, Liao H, et al. A two-pollutant strategy for improving ozone and particulate air quality in China[J]. Nature Geoscience, 2019, 12(11): 906-910.
[44] 张琴, 莫慧偲, 周国治, 等. 基于转置EKMA曲线的湖南省臭氧生成敏感性分析[J]. 中国环境监测, 2023, 39(5): 105-115.
Zhang Q, Mo H S, Zhou G Z, et al. Study of ozone formation sensitivity in Hunan province based on the transformed-EKMA curve approach[J]. Environmental Monitoring in China, 2023, 39(5): 105-115.
[45] Jacob D J, Winner D A. Effect of climate change on air quality[J]. Atmospheric Environment, 2009, 43(1): 51-63.