环境科学  2023, Vol. 44 Issue (8): 4241-4249   PDF    
天津市"十三五"期间O3污染特征和驱动因子
李源1, 肖致美1, 毕晓辉2, 蔡子颖3, 徐虹1, 高璟赟1, 郑乃源1, 杨宁1     
1. 天津市生态环境监测中心, 天津 300191;
2. 南开大学环境科学与工程学院, 国家环境保护城市空气颗粒物污染防治重点实验室, 天津 300071;
3. 天津市环境气象中心, 天津 300074
摘要: 为了解"十三五"期间天津市O3污染特征和驱动因素, 基于2016~2020年高时空分辨率的在线监测数据, 利用空间自相关、空间热点分析和STIRPAT模型分析了O3污染空间分布、聚集特征和驱动因子.结果表明, 2016~2020年天津市O3浓度变化特征呈现污染发生时间点提前和污染范围扩大的趋势.6~10月O3污染分布具有显著聚集性, 高值-高值聚集区主要为市内六区、北辰区、津南区和静海区, O3浓度在西南部地区形成高值热点聚集区, 在东北部地区形成低值冷点聚集区.气温、小风百分率和日照时数等气象因子与NOx排放量、VOCs排放量和机动车保有量等社会因子对O3浓度有显著性影响, 综合驱动STIRPAT模型的回归拟合效果比单一气象因子或社会因子模型更好.为科学高效地开展"十四五"期间O3污染的防治, 在关注气象条件基础上, 在"双碳"目标的约束下, 天津市应进一步提升钢铁、石化、火电和建材等行业全过程排放的绩效水平, 引导企业清洁化提升改造、转型升级和绿色发展, 减少企业VOCs和NOx的排放量, 同时进一步控制燃油机动车保有量的增加, 大力推广新能源机动车, 减少机动车尾气的排放.
关键词: "十三五"期间      天津市      空间自相关      空间热点分析      STIRPAT模型     
Characteristics and Driving Factors of O3 Pollution During 13th Five-Year Period in Tianjin
LI Yuan1 , XIAO Zhi-mei1 , BI Xiao-hui2 , CAI Zi-ying3 , XU Hong1 , GAO Jing-yun1 , ZHENG Nai-yuan1 , YANG Ning1     
1. Tianjin Eco-Environmental Monitoring Center, Tianjin 300191, China;
2. State Environmental Protection Key Laboratory of Urban Ambient Air Particulate Matter Pollution Prevention and Control, College of Environmental Science and Engineering, Nankai University, Tianjin 300071, China;
3. Tianjin Environmental Meteorological Center, Tianjin 300074, China
Abstract: The spatial distribution, accumulation features, and driving factors of O3 pollution were analyzed using spatial autocorrelation and hotspot analysis and the STIRPAT model based on the high spatiotemporal resolution online monitoring data from 2016 to 2020 in Tianjin. The results showed that the variation characteristics of O3 concentration in Tianjin from 2016 to 2020 had the trend of pollution occurring in advance and the scope of the pollution expanding. The distribution of O3 pollution showed significant aggregation from June to October. High-high value clustering areas included six urban districts, Beichen District, Jinnan District, and Jinghai District. O3 concentration formed high value hot spots in the southwest and low value cold spots in the northeast. Meteorological factors such as temperature, breeze percentage, and sunshine duration, as well as social factors such as NOx emission, VOCs emission, and motor vehicle ownership had significant effects on O3 concentration. The regression fitting effect of the integrated drive STIRPAT model was better than that of the single meteorological factor or social factor models. In order to promote scientific and efficient prevention and control of ozone pollution during the 14th Five-Year Plan period, meteorological conditions require attention; under the goal of "peaking carbon dioxide emissions and achieving carbon neutrality, " it is necessary for Tianjin to further improve the emission performance of steel, petrochemicals, thermal power, building materials, and other industries, Additionally, clean upgrading, transformation, and green development should be guided for enterprises to reduce VOCs and NOx emissions. At same time, the increase in fuel vehicle numbers should be controlled, and new energy vehicles should be vigorously promoted to reduce vehicle emissions.
Key words: 13th Five-Year Period      Tianjin      spatial autocorrelation      spatial hotspot analysis      STIRPAT model     

随着《大气污染防治行动计划》和《打赢蓝天保卫战三年行动计划》深入实施, “十三五”期间, 我国城市PM2.5平均浓度持续下降[1, 2], 但城市臭氧(O3)污染日益凸显, 在全国多个区域都没能得到有效控制.O3浓度水平上升, 污染范围扩大, 已逐渐成为制约空气质量持续改善的主要因素[3~6].对于京津冀、长三角、珠三角和成渝地区等特大型城市群及其周边地区而言, O3污染现已成为困扰经济快速发展地区的重要环境问题[7, 8].

近地面O3主要是人类活动排放的氮氧化物(NOx)和挥发性有机物(VOCs)等前体物在光照下发生复杂光化学反应的产物[9~16], 其浓度水平一方面受气象条件影响, 另一方面与社会因素密切相关[17, 18].一般认为, 影响O3的气象因子主要有气温、紫外辐射、相对湿度、边界层高度和风等[19~26], O3污染多出现在高温、低湿、强辐射的天气[27~29].王玫等[30]分析京津冀O3变化特征和与气象要素的关系表明, 当气温大于30℃、相对湿度为50% ~60%、风速为2.0~3.0 m·s-1时, O3超标概率最高.在社会因子上, Li等[31]利用GEOS-Chem模型发现NOx排放对O3浓度有一定的驱动作用; 柯碧钦等[32]指出第二产业占GDP的比重是华北地区O3浓度的主要社会经济驱动因素, 工业生产用电量和工业SO2排放量对O3浓度变化也有影响; 胡雪萍等[33]依据STIRPAT模型分析对大气污染产生影响的关键社会因素, 指出人口密度对大气污染未产生显著的影响, 能源效率对大气污染产生显著的负向抑制作用, 国民收入增长带来大气污染的增加.

现有研究成果对O3驱动因子的分析大多从单一方面入手, 在城市尺度上气象因子和社会因子综合驱动作用对O3影响的研究仍显缺乏, 且研究主要集中在四川盆地和华北地区等较大尺度范围区域[32, 34].天津市位于京津冀区域的中北部, 工业发达, 是京津冀区域大气污染较重的城市之一, 2020年天津市O3日最大8h平均浓度第90百分位数ρ(O3)达190 μg·m-3, 较2013年上升25.8%, 相关研究主要集中在O3和前体物污染特征和来源上[5, 35~37], 缺少O3驱动因子研究.本文以“十三五”期间天津市O3监测数据为基础, 运用空间自相关分析和空间热点分析, 分析O3污染的空间分布和聚集规律, 并利用STIRPAT模型对O3驱动因子进行研究, 比较不同驱动因子对O3的影响, 以期为“十四五”期间天津市O3污染的治理与防控提供理论依据.

1 材料与方法 1.1 采样和分析方法

O3浓度数据来源于天津市环境空气质量监测网络, 监测点位分布见图 1.采用Thermo Model 49i监测O3, 时间分辨率为1 h.仪器每天进行自动校准, 每周人工校准一次, 质量控制均严格按照《环境空气气态污染物(SO2、NO2、O3、CO)连续自动监测系统运行和质控技术规范》(HJ 818-2018)要求进行.本文使用的O3浓度数据为O3日最大8 h滑动平均值(O3-8h).

图 1 监测点位分布示意 Fig. 1 Location of the sampling station

1.2 驱动因子数据

本研究从气象因子和社会因子两方面选取驱动因子, 其中气象因子选取气温、相对湿度、小风百分率、紫外辐射、日照时数、混合层厚度、静稳指数和逆温百分率, 观测数据来自中国气象局天津大气边界层观测站(台站编号为54517); 社会因子选取SO2排放量、NOx排放量、烟粉尘排放量、VOCs排放量、人口密度、生产总值、第二产业比重和机动车保有量, VOCs排放量数据来自文献[38], 其他数据来源于天津市统计局2016~2020年《天津统计年鉴》和《天津市国民经济和社会发展统计公报》.

1.3 研究方法 1.3.1 空间自相关分析

采用全局Moran's I评价天津市O3浓度在空间上的关联程度, 判断是否有聚集、分异或随机等空间分布特征存在.全局Moran's I取值范围在±1之间, 大于0表示空间正相关, O3浓度高低值存在聚集性; 小于0表示空间负相关, O3浓度高低值存在分异性; 等于0表示整个区域O3浓度随机分布.I值越接近1, 表明在空间分布上呈现出的聚集效应越显著, I值越接近-1, 表明在空间分布上呈现出的分异性质越显著.其计算公式如下:

(1)

式中, n为样本数; xixj分别为监测站ij的O3-8h; x为所有监测站O3-8h的平均值; Wij为权重矩阵, 由于本研究的O3-8h是点文件, 故选择欧氏距离权重生成空间权重矩阵.

全局Moran's I只能判断整个区域O3浓度在空间是否存在自相关性, 需要采用局部Moran's I探测聚集或分异的具体位置, 评价某一局部区域与其周边临近区域的浓度高低值聚集情况.将基于全局空间自相关分析的统计方法分解到局部空间上[39], 得到局部Moran's I计算公式如下:

(2)

式中, 各符号含义与式(1)相同.本研究运用GeoDa软件来计算Ii值.

1.3.2 空间热点分析

为进一步分析天津市O3浓度高值或低值在空间上发生聚集的位置和范围, 采用空间热点分析(Getis-Ord Gi*)研究局部空间聚类分布特征, 更详细地表示出是高值聚集(热点), 还是低值聚集(冷点).热点分析不同于全局Moran's I, 其对于O3浓度在空间上随机分布的情况也可能探测出某种程度的聚集, 也不同于局部Moran's I只针对某一局部区域, 而是可以同时反映某几个区域的高低值(冷热点)分布.空间热点分析可在整个研究区域内探测出O3浓度显著高于或低于其他地方的区域, 将空间分布聚集的程度通过热点与冷点进行区分[40].公式如下:

(3)

式中, S为样本标准差, 其他符号含义与式(1)相同.本研究运用ArcGIS软件进行空间热点探测分析.

1.3.3 STIRPAT模型

采用STIRPAT模型分析O3驱动因素, 该模型由IPAT模型演变而来[41].为克服IPAT模型中各因素等比例影响的不足, York等基于IPAT模型构建了STIRPAT模型[42], 其表达式为:

(4)

为消除模型中可能存在的异方差影响, 本研究将所有变量进行对数化处理, 对数化之后的扩展STIRPAT模型如下:

(5)

式中, I为大气污染情况; P为人口规模; A为富裕程度; T为技术水平; a为常数; bcd分别为PAT的指数项; e为误差项.

在使用STIRPAT模型时, 由于气象因子或社会因子之间具有相关的共同趋势, 会产生多重共线问题, 导致估计失真或难以估计准确.本研究在构建模型之前采用方差膨胀因子(VIF)检验多重共线性, 通常以10作为判断边界, 当VIF<10时, 不存在多重共线性, 当VIF≥10时, 存在较强多重共线性[43].

2 结果与讨论 2.1 O3浓度时空变化特征

“十三五”期间, 天津市ρ(O3)由2016年的157 μg·m-3升至2020年的190 μg·m-3, 累计上升31.2%(图 2), 2017~2020年连续4 a超过国家二级浓度限值标准(160 μg·m-3), 其中2018年达到峰值(201 μg·m-3).从月浓度分布上看, 天津市O3浓度在每年的1~5月上升, 6月达到全年峰值, 7月后呈现下降趋势, 12月达到全年谷值, 与大多数城市变化特征相同[44]. 5~9月天津市气温较高, 日最高气温大于30℃的天数多[45], 大气光化学反应强, 该时期各年O3浓度均较高, 各月均超过160 μg·m-3.“十三五”后期(2018~2020年), 4月O3浓度明显上升, 各年ρ(O3)均超标, 分别为208、162和190 μg·m-3, 与“十三五”后期4月降水相对较少, 日最高气温大于30℃的天数明显增加有关[45].“十三五”期间, 天津市O3月浓度变化特征反映出O3污染发生时间节点提前, O3污染加重的特点.

图 2 2016~2020年天津市O3浓度 Fig. 2 O3 concentration from 2016 to 2020 in Tianjin

从空间分布上看(图 3), “十三五”初期(2016年), 天津市西部和中部地区O3浓度较高, 东南部区域O3浓度相对较低; 2017年开始, 天津市16个区O3浓度均超过国家二级标准限值, 高值区域呈现明显扩大趋势, 其中中部、西南部、西部和北部区域O3浓度均较高, 东南部区域O3浓度相对较低.天津市O3浓度空间分布特征与污染源排放和气象因素的分布有关, 中心城区机动车数量和密度高, 上下班高峰期交通堵塞严重, 机动车排放量大, 同时餐饮油烟、干洗、加油站和市政有机溶剂使用等面源排放影响较高, 受城市“热岛效应”影响, 温度相对较高, O3污染较重; 中部的环城区域钢铁、火电和建材等行业O3前体物排放量大, 西南部静海的黑色金属冶炼和压延加工业、医药、家具和运输设备制造业以及西部的武清区域的汽车制造、化学品制造、医药和家具行业等行业VOCs排放量大, O3污染较重; 北部的蓟州区O3污染较重与植物源排放量大有关.东南部的滨海新区石化行业排放量大, 但受天津港船舶和柴油车货运影响, NO排放量高, O3可以通过NO滴定效应去除[46], 较高的NO对O3生成有明显的抑制作用[47, 48], 同时东南部区域受海陆风影响, 东南风发生频率高[36], 大气扩散条件相对较好, O3污染相对较轻.

各区域对应行政区名称与图 1相同, 下同 图 3 2016~2020年天津市O3浓度空间分布 Fig. 3 Spatial distribution of O3 concentration from 2016 to 2020 in Tianjin

2.2 空间自相关分析

“十三五”期间天津市O3浓度全局空间自相关分析结果显示(表 1), 6~10月各月的全局Moran's I均为正值且均通过显著性检验(P<0.001), 说明6~10月天津市O3浓度呈现出高低值聚集的态势, 其分布明显受自身和邻近区域之间的影响; 其余各月全局Moran's I接近0且均不显著, 说明区域O3浓度之间相互独立.因此, 6~10月天津市O3污染的防控需要各区之间密切合作和协同治理.

表 1 2016~2020年天津市O3浓度全局Moran's I Table 1 Global Moran's I of O3 concentration in Tianjin from 2016 to 2020

为进一步说明哪些区域出现聚集和其对天津市全局空间自相关的贡献, 对“十三五”期间6~10月的O3浓度进行局部自相关分析, 结果如图 4, 其中Moran's I散点图的第一、第二、第三和第四象限分别为高值-高值聚集区(H-H, 即该区域与其周边区域O3浓度均为高值)、低值-高值聚集区(L-H, 即该区域O3浓度较低, 但其周边O3浓度较高)、低值-低值聚集区(L-L, 即该区域与其周边区域O3浓度均为低值)和高值-低值聚集区(H-L, 即该区域O3浓度较高, 但其周边O3浓度较低). H-H和L-L聚集区表明O3污染呈正空间自相关, O3浓度呈现高低值聚集性; H-L和L-H聚集区表明O3污染呈负空间自相关, O3浓度呈现高低值分异性. “十三五”期间6~10月天津市O3浓度空间分布呈正相关聚集, H-H和L-L显著聚集区占62.5%, 其中H-H聚集区为市内六区、北辰区、津南区和静海区, 是O3污染急需集中治理的区域.

图 4 “十三五”期间6~10月天津市O3浓度局部Moran's I散点图和空间分布 Fig. 4 Local Moran's I scatter plot and spatial distribution of O3 concentration in Tianjin from June to October during 13th Five-Year Period

2.3 空间冷热点分析

基于空间自相关分析结果, 进行冷热点探测, 发现“十三五”期间6~10月天津市O3浓度整体表现出强聚集性, 其中市内六区、西青区和津南区呈现99%置信度的热点聚集, 北辰区和静海区呈现95%置信度的热点聚集, 这些区域的O3浓度显著高于其他地方; 宝坻区呈现99%置信度的冷点聚集, 蓟州区和宁河区呈现95%置信度的冷点聚集, 这些区域的O3浓度显著低于其他地方(图 5).由此可见, 天津市O3浓度冷热点空间分布非常明显, 即西南部地区形成浓度高值聚集区域, O3污染形势严重, 是防治的重点地区, 东北部地区形成浓度低值聚集区域, O3污染相对较轻.

图 5 “十三五”期间6~10月天津市O3浓度冷热点区域 Fig. 5 Cold and hot spot area of O3 concentration in Tianjin from June to October during the 13th Five-Year Period

2.4 O3驱动因子分析 2.4.1 相关性分析与多重共线性检验

本研究选取的8个气象因子均为与O3浓度相关的气象要素, 但在构建STIRPAT模型之前, 首先对气象因子和社会因子进行相关性分析, 以确保因子独立性.气象因子之间的相关性分析表明(表 2), 紫外辐射、混合层厚度和静稳指数与其他气象因子的相关性较高, 均存在4个相关系数绝对值大于0.7的情况, 另外5个气象因子之间的相关系数较低.在此基础上, 对气温、相对湿度、小风百分率、日照时数和逆温百分率进行多重共线性检验, VIF分别为6.63、6.05、2.75、4.79和1.85, VIF均小于10, 说明5个气象因子之间不存在多重共线性, 不会由于相互作用而影响模型的准确性和稳定性.

表 2 气象因子之间的相关系数 Table 2 Correlation coefficient between meteorological factors

社会因子之间的相关性分析表明(表 3), 烟粉尘排放量、人口密度和生产总值与其他社会因子的相关性较高, 均存在4个相关系数绝对值大于0.7的情况, 另外5个社会因子之间的相关系数较低.对SO2排放量、NOx排放量、VOCs排放量、第二产业占比和机动车保有量进行多重共线性检验, VIF分别为5.18、3.61、2.63、1.28和4.02, VIF均小于10, 说明5个社会因子之间共同趋势不大, 不存在多重共线性.

表 3 社会因子之间的相关系数 Table 3 Correlation coefficient between social factors

2.4.2 气象和社会单因子驱动分析

利用式(5)建立O3浓度与气温、相对湿度、小风百分率、日照时数和逆温百分率等气象因子STIRPAT模型, 模型调整后的R2为0.868、标准误差为0.188, 回归拟合效果好.模型结果可知(表 4), 气温、小风百分率和日照时数通过了显著性检验(P<0.05), 说明气温的上升、小风百分率的加大和日照时数的增多均对O3浓度的升高有明显驱动作用.气象因子中, 气温对O3浓度的影响最大, 较高的气温一方面有助于生物源VOCs等O3前体物的挥发, 另一方面能够提供光化学反应条件, 使反应速率大大提高, 驱动O3的生成[49].日照时数的增多一定程度上意味着光照时间的延长、紫外辐射的增加, 在为生成O3的光化学反应提供良好条件的同时, 也延长了反应时间, 有利于O3的生成与累积.小风百分率的加大会减弱大气湍流作用, 降低大气输送能力, 造成本地累积.

表 4 气象因子和社会因子STIRPAT模型结果 Table 4 STIRPAT model results of meteorological factors and social factors

建立O3浓度与SO2排放量、NOx排放量、VOCs排放量、第二产业占比和机动车保有量等社会因子STIRPAT模型, 模型调整后的R2为0.740, 标准误差为0.264, 回归拟合效果较好.根据模型结果(表 4), NOx排放量、VOCs排放量和机动车保有量通过了显著性检验(P<0.05), 这3个社会因子对O3有明显的驱动作用, 其中机动车保有量和VOCs排放量对O3的影响较大.“十三五”期间, 天津市机动车保有量由2016年的281.8万辆增加至2020年的338.5万辆[50], 机动车尾气中的NOx和VOCs是O3的主要前体物, 汽车保有量大的城市通常有更高的O3浓度[51~53].同时, 天津市O3的生成处于VOCs控制区[45, 54], “十三五”期间的VOCs排放量仍处于高位[55], 对O3的驱动作用明显.

2.4.3 多因子综合驱动分析

选取对O3有显著性影响的6个因子进行多重共线性检验, 得到气温、小风百分率、日照时数、NOx排放量、VOCs排放量和机动车保有量的VIF分别为5.02、2.40、2.42、3.55、4.54和2.17, VIF均小于10, 说明因子之间不存在多重共线性, 可以确保模型的准确性和稳定性.建立的综合驱动STIRPAT模型, 调整后的R2为0.891, 标准误差为0.169, 回归拟合效果比单一建立气象因子或社会因子模型更好(表 5).结果表明, 气温、日照时数、NOx排放量、VOCs排放量和机动车保有量均通过了显著性检验(P<0.05), 其中气温与O3浓度呈极显著的正相关关系, 是最重要的O3驱动因子.为科学高效地开展“十四五”期间O3污染的防治, 在关注气象条件的基础上, 在“双碳”目标的约束下, 天津市应进一步提升钢铁、石化、火电和建材等行业全过程排放的绩效水平, 引导企业清洁化提升改造、转型升级和绿色发展, 减少企业VOCs和NOx的排放量; 同时, 要进一步控制燃油机动车保有量的增加, 大力推广新能源机动车, 减少机动车尾气的排放.

表 5 综合驱动STIRPAT模型结果 Table 5 STIRPAT model results of integrated driving

3 结论

(1)“十三五”期间天津市O3浓度整体呈上升趋势, O3污染主要集中在每年的5~9月, 中部、西南部、西部和北部区域O3浓度相对较高, 变化特征呈现出O3污染发生时间点提前、污染范围扩大的趋势.

(2) 空间自相关分析表明, 6~10月天津市O3浓度全局Moran's I为正值且通过显著性检验(P<0.001), O3污染分布具有明显的聚集性, H-H聚集区主要为市内六区、北辰区、津南区和静海区, 是O3污染急需集中治理的区域.

(3) 空间热点分析表明, 6~10月天津市西南部地区形成O3浓度高值热点聚集区域, O3污染形势严重, 是防治的重点地区, 东北部地区形成O3浓度低值冷点聚集区域, O3污染相对较轻.

(4) 气象因子中, 气温的上升、小风百分率的加大和日照时数的增多对O3浓度的升高有显著的促进作用; 社会因子中, NOx排放量、VOCs排放量和机动车保有量对O3有显著的驱动作用.综合驱动STIRPAT模型回归拟合效果比单一气象因子或社会因子模型更好, 气温、日照时数、NOx排放量、VOCs排放量和机动车保有量是驱动O3浓度变化的重要因子.

参考文献
[1] 肖致美, 徐虹, 李立伟, 等. 基于在线观测的天津市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.
[2] 李源, 陈魁, 孔君, 等. 天津市PM2.5污染特征与来源解析[J]. 环境工程, 2019, 37(11): 132-137.
Li Y, Chen K, Kong J, et al. Pollution characteristics and sources analysis of PM2.5 in Tianjin[J]. Environmental Engineering, 2019, 37(11): 132-137.
[3] Ma Z Q, Xu J, Quan W J, et al. Significant increase of surface ozone at a rural site, North of eastern China[J]. Atmospheric Chemistry and Physics, 2016, 16(6): 3969-3977. DOI:10.5194/acp-16-3969-2016
[4] 王雪涵, 张文慧, 毕晓辉, 等. 2001—2020年天津市大气污染特征的演变与防治历程[J]. 环境科学研究, 2022, 35(4): 945-955.
Wang X H, Zhang W H, Bi X H, et al. Characteristics evolution and prevention development of ambient pollution in Tianjin, China[J]. Research of Environmental Sciences, 2022, 35(4): 945-955. DOI:10.13198/j.issn.1001-6929.2021.12.21
[5] 李源, 孔君, 徐虹, 等. 天津市臭氧污染特征及来源解析研究[J]. 环境污染与防治, 2019, 41(6): 647-651.
Li Y, Kong J, Xu H, et al. Analysis of ozone pollution characteristic in Tianjin and its source apportionment[J]. Environmental Pollution & Control, 2019, 41(6): 647-651. DOI:10.15985/j.cnki.1001-3865.2019.06.006
[6] 姜华, 常宏咪. 我国臭氧污染形势分析及成因初探[J]. 环境科学研究, 2021, 34(7): 1576-1582.
Jiang H, Chang H M. Analysis of China's ozone pollution situation, preliminary investigation of causes and prevention and control recommendations[J]. Research of Environmental Sciences, 2021, 34(7): 1576-1582. DOI:10.13198/j.issn.1001-6929.2021.05.20
[7] 李霄阳, 李思杰, 刘鹏飞, 等. 2016年中国城市臭氧浓度的时空变化规律[J]. 环境科学学报, 2018, 38(4): 1263-1274.
Li X Y, Li S J, Liu P F, et al. Spatial and temporal variations of ozone concentrations in China in 2016[J]. Acta Scientiae Circumstantiae, 2018, 38(4): 1263-1274. DOI:10.13671/j.hjkxxb.2017.0399
[8] 杜云松, 罗彬, 陈建文, 等. 气温在成都地区臭氧预报的运用研究[J]. 环境科学与技术, 2017, 40(S1): 329-334.
Du Y S, Luo B, Chen J W, et al. Study on the application of air temperature in ozone forecast in Chengdu area[J]. Environmental Science & Technology, 2017, 40(S1): 329-334.
[9] Doherty R M. Atmospheric chemistry: ozone pollution from near and far[J]. Nature Geoscience, 2015, 8(9): 664-665. DOI:10.1038/ngeo2497
[10] Fishman J, Crutzen P J. The origin of ozone in the troposphere[J]. Nature, 1978, 274(5674): 855-858. DOI:10.1038/274855a0
[11] Wang Z S, Li Y T, Chen T, et al. Ground-level ozone in urban Beijing over a 1-year period: temporal variations and relationship to atmospheric oxidation[J]. Atmospheric Research, 2015, 164-165: 110-117. DOI:10.1016/j.atmosres.2015.05.005
[12] Von Schneidemesser E, Monks P S, Allan J D, et al. Chemistry and the linkages between air quality and climate change[J]. Chemical Reviews, 2015, 115(10): 3856-3897. DOI:10.1021/acs.chemrev.5b00089
[13] Zhang Y N, Xiang Y R, Chan L Y, et al. Procuring the regional urbanization and industrialization effect on ozone pollution in Pearl River Delta of Guangdong, China[J]. Atmospheric Environment, 2011, 45(28): 4898-4906. DOI:10.1016/j.atmosenv.2011.06.013
[14] Liu T, Li T T, Zhang Y H, et al. The short-term effect of ambient ozone on mortality is modified by temperature in Guangzhou, China[J]. Atmospheric Environment, 2013, 76: 59-67. DOI:10.1016/j.atmosenv.2012.07.011
[15] Iqbal M A, Kim K H, Shon Z H, et al. Comparison of ozone pollution levels at various sites in Seoul, a megacity in Northeast Asia[J]. Atmospheric Research, 2014, 138: 330-345. DOI:10.1016/j.atmosres.2013.12.003
[16] Shi C Z, Wang S S, Liu R, et al. A study of aerosol optical properties during ozone pollution episodes in 2013 over Shanghai, China[J]. Atmospheric Research, 2015, 153: 235-249. DOI:10.1016/j.atmosres.2014.09.002
[17] Li Y, Wang J G, Li L W, et al. Research on the influence of weather patterns on ozone concentration: a case study in Tianjin[J]. Atmosphere, 2022, 13(8). DOI:10.3390/atmos13081312
[18] Jaén C, Udina M, Bech J. Analysis of two heat wave driven ozone episodes in Barcelona and surrounding region: Meteorological and photochemical modeling[J]. Atmospheric Environment, 2021, 246. DOI:10.1016/j.atmosenv.2020.118037
[19] 白宇, 李源, 郭瑞, 等. 气象因子对天津市臭氧浓度的影响分析[J]. 环境影响评价, 2022, 44(2): 14-19.
Bai Y, Li Y, Guo R, et al. Influence analysis of meteorological factors on ozone concentration in Tianjin[J]. Environmental Impact Assessment, 2022, 44(2): 14-19. DOI:10.14068/j.ceia.2022.02.003
[20] 齐冰, 牛彧文, 杜荣光, 等. 杭州市近地面大气臭氧浓度变化特征分析[J]. 中国环境科学, 2017, 37(2): 443-451.
Qi B, Niu Y W, Du R G, et al. Characteristics of surface ozone concentration in urban site of Hangzhou[J]. China Environmental Science, 2017, 37(2): 443-451.
[21] 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
[22] Xu X, Lin W, Wang T, et al. Long-term trend of surface ozone at a regional background station in eastern China 1991-2006: enhanced variability[J]. Atmospheric Chemistry and Physics, 2008, 8(10): 2595-2607. DOI:10.5194/acp-8-2595-2008
[23] 周贺玲, 周玉都, 闻静. 河北廊坊近地面层O3特征及其影响因素[J]. 干旱气象, 2017, 35(3): 405-411.
Zhou H L, Zhou Y D, Wen J. Characteristics of ozone over surface layer of Langfang City and its influencing meteorological factors[J]. Journal of Arid Meteorology, 2017, 35(3): 405-411.
[24] 崔梦瑞, 白林燕, 冯建中, 等. 京津唐地区臭氧时空分布特征与气象因子的关联性研究[J]. 环境科学学报, 2020, 40(12): 1-13.
Cui M R, Bai L Y, Feng J Z, et al. Analysis of temporal and spatial variations of ozone coupling with dynamics of meteorological factors in the Beijing-Tianjin-Tangshan region[J]. Acta Scientiae Circumstantiae, 2020, 40(12): 1-13.
[25] Yang L F, Xie D P, Yuan Z B, et al. Quantification of regional ozone pollution characteristics and its temporal evolution: Insights from identification of the impacts of meteorological conditions and emissions[J]. Atmosphere, 2021, 12(2). DOI:10.3390/atmos12020279
[26] Ma M L, Yao G B, Guo J P, et al. Distinct spatiotemporal variation patterns of surface ozone in China due to diverse influential factors[J]. Journal of Environmental Management, 2021, 288. DOI:10.1016/j.jenvman.2021.112368
[27] Solomon P, Cowling E, Hidy G, et al. Comparison of scientific findings from major ozone field studies in North America and Europe[J]. Atmospheric Environment, 2000, 34(12-14): 1885-1920. DOI:10.1016/S1352-2310(99)00453-7
[28] Colette A, Andersson C, Baklanov A, et al. Is the ozone climate penalty robust in Europe?[J]. Environmental Research Letters, 2015, 10(8). DOI:10.1088/1748-9326/10/8/084015
[29] Ejimofor C S, Okoro E C, Sivla W T. Effects of elevated humidity on stratospheric ozone content in the tropics[J]. International Journal of Physical Sciences, 2021, 15(4): 182-193.
[30] 王玫, 郑有飞, 柳艳菊, 等. 京津冀臭氧变化特征及与气象要素的关系[J]. 中国环境科学, 2019, 39(7): 2689-2698.
Wang M, Zheng Y F, Liu Y J, et al. Characteristics of ozone and its relationship with meteorological factors in Beijing-Tianjin-Hebei Region[J]. China Environmental Science, 2019, 39(7): 2689-2698. DOI:10.19674/j.cnki.issn1000-6923.2019.0318
[31] Li K, Jacob D J, Liao H, et al. Anthropogenic drivers of 2013-2017 trends in summer surface ozone in China[J]. Proceedings of the National Academy of Science of the United States of America, 2019, 116(2): 422-427. DOI:10.1073/pnas.1812168116
[32] 柯碧钦, 何超, 杨璐, 等. 华北地区地表臭氧时空分布特征及驱动因子[J]. 中国环境科学, 2022, 42(4): 1562-1574.
Ke B Q, He C, Yang L, et al. The spatiotemporal variation of surface ozone and the main driving factors in North China[J]. China Environmental Science, 2022, 42(4): 1562-1574. DOI:10.3969/j.issn.1000-6923.2022.04.009
[33] 胡雪萍, 方永丽. 中国大气污染的影响因素及防治措施研究——基于STIRPAT模型和固定效应面板模型[J]. 工业技术经济, 2018, 37(2): 107-113.
Hu X P, Fang Y L. Study on influencing factors and preventive measures of atmospheric pollution in China-based on STIRPAT model and fixed effect panel model[J]. Journal of Industrial Technological Economics, 2018, 37(2): 107-113.
[34] 汪可可, 康平, 周明卫, 等. 四川盆地臭氧浓度空间分异及驱动因子研究[J]. 中国环境科学, 2020, 40(6): 2361-2370.
Wang K K, Kang P, Zhou M W, et al. Spatial differentiation and driving factors of ozone concentration in Sichuan Basin[J]. China Environmental Science, 2020, 40(6): 2361-2370.
[35] 王文美, 高璟赟, 肖致美, 等. 天津市夏季不同臭氧浓度级别VOCs特征及来源[J]. 环境科学, 2021, 42(8): 3585-3594.
Wang W M, Gao J Y, Xiao Z M, et al. Characteristics and sources of VOCs at different ozone concentration levels in Tianjin[J]. Environmental Science, 2021, 42(8): 3585-3594.
[36] 关玉春, 肖致美, 陈魁, 等. 天津市臭氧污染特征与影响因素分析[J]. 中国环境监测, 2017, 33(4): 40-49.
Guan Y C, Xiao Z M, Chen K, et al. Characteristics of ozone pollution and its influencing factors in Tianjin[J]. Environmental Monitoring in China, 2017, 33(4): 40-49.
[37] 高璟赟, 肖致美, 徐虹, 等. 2019年天津市挥发性有机物污染特征及来源[J]. 环境科学, 2021, 42(1): 55-64.
Gao J Y, Xiao Z M, Xu H, et al. Characterization and source apportionment of atmospheric VOCs in Tianjin in 2019[J]. Environmental Science, 2021, 42(1): 55-64.
[38] 肖致美, 徐虹, 蔡子颖, 等. 天津市"十三五"期间PM2.5减排效果评估[J]. 环境科学, 2023, 44(6): 3054-3062.
Xiao Z M, Xu H, Cai Z Y, et al. Assessment of emission reduction effect of major air pollution control measures on PM2.5 concentrations during 13th Five-Year Period in Tianjin[J]. Environmental Science, 2023, 44(6): 3054-3062.
[39] Anselin L. Local indicators of Spatial association—LISA[J]. Geographical Analysis, 1995, 27(2): 93-115.
[40] 周磊, 武建军, 贾瑞静, 等. 京津冀PM2.5时空分布特征及其污染风险因素[J]. 环境科学研究, 2016, 29(4): 483-493.
Zhou L, Wu J J, Jia R J, et al. Investigation of temporal-spatial characteristics and underlying risk factors of PM2.5 pollution in Beijing-Tianjin-Hebei Area[J]. Research of Environmental Sciences, 2016, 29(4): 483-493.
[41] Ehrlich P R, Holdren J P. Impact of population growth[J]. Science, 1971, 171(3977): 1212-1217.
[42] York R, Rosa E A, Dietz T. STIRPAT, IPAT and ImPACT: analytic tools for unpacking the driving forces of environmental impacts[J]. Ecological Economics, 2003, 46(3): 351-365.
[43] Marquaridt D W. Generalized inverses, ridge regression, biased linear estimation, and nonlinear estimation[J]. Technometrics, 1970, 12(3): 591-612.
[44] 王笑哲, 赵莎, 郭灵辉, 等. 京津冀及周边地区"2+26"城市臭氧的季节性变化规律[J]. 环境科学研究, 2022, 35(8): 1786-1797.
Wang X Z, Zhao S, Guo L H, et al. Seasonal variation of ozone in '2+26' cities in Beijing-Tianjin-Hebei region and surrounding areas[J]. Research of Environmental Sciences, 2022, 35(8): 1786-1797.
[45] 肖致美, 李源, 孔君, 等. 天津市PM2.5-O3复合污染特征及气象影响分析[J]. 环境科学, 2022, 43(6): 2928-2936.
Xiao Z M, Li Y, Kong J, et al. Characteristics and meteorological factors of PM2.5-O3 compound pollution in Tianjin[J]. Environmental Science, 2022, 43(6): 2928-2936.
[46] Seguel R J, Morales S R G E, Leiva G M A. Ozone weekend effect in Santiago, Chile[J]. Environmental Pollution, 2012, 162: 72-79.
[47] 安俊琳, 王跃思, 李昕, 等. 北京大气中NO、NO2和O3浓度变化的相关性分析[J]. 环境科学, 2007, 28(4): 706-711.
An J L, Wang Y S, Li X, et al. Analysis of the relationship between NO, NO2 and O3 concentrations in Beijing[J]. Environmental Science, 2007, 28(4): 706-711.
[48] 毛敏娟, 刘厚通, 杜荣光. 不同时间尺度下杭州市O3污染特征及控制因素[J]. 环境科学研究, 2019, 32(11): 1844-1851.
Mao M J, Liu H T, Du R G. Characteristics and control factors of ozone pollution at different time scales in Hangzhou city[J]. Research of Environmental Sciences, 2019, 32(11): 1844-1851.
[49] Gong X, Hong S, Jaffe D A. Ozone in China: spatial distribution and leading meteorological factors controlling O3 in 16 Chinese cities[J]. Aerosol and Air Quality Research, 2018, 18(9): 2287-2300.
[50] 天津市统计局, 国家统计局天津调查总队. 2017-2021年天津统计年鉴[EB/OL]. https://stats.tj.gov.cn/tjsj_52032/tjnj/, 2022-08-30.
[51] 陈天增, 葛艳丽, 刘永春, 等. 我国机动车排放VOCs及其大气环境影响[J]. 环境科学, 2018, 39(2): 478-492.
Chen T Z, Ge Y L, Liu Y C, et al. VOCs Emission from motor vehicles in China and its impact on the atmospheric environment[J]. Environmental Science, 2018, 39(2): 478-492.
[52] Zhang Y L, Yang W Q, Simpson I, et al. Decadal changes in emissions of volatile organic compounds (VOCs) from on-road vehicles with intensified automobile pollution control: case study in a busy urban tunnel in South China[J]. Environmental Pollution, 2018, 233: 806-819.
[53] 王鑫龙, 赵文吉, 李令军, 等. 中国臭氧时空分布特征及与社会经济因素影响分析[J]. 地球与环境, 2020, 48(1): 66-75.
Wang X L, Zhao W J, Li L J, et al. Characteristics of spatiotemporal distribution of O3 in China and impact analysis of socio-economic factors[J]. Earth and Environment, 2020, 48(1): 66-75.
[54] Jin X M, Holloway T. Spatial and temporal variability of ozone sensitivity over China observed from the ozone monitoring instrument[J]. Journal of Geophysical Research Atmospheres, 2015, 120(14): 7229-7246.
[55] 肖致美, 徐虹, 高璟贇, 等. 天津市PM2.5-O3复合污染特征及来源分析[J]. 环境科学, 2022, 43(3): 1140-1150.
Xiao Z M, Xu H, Gao J Y, et al. Characteristics and Sources of PM2.5-O3 compound pollution in Tianjin[J]. Environmental Science, 2022, 43(3): 1140-1150.
天津市"十三五"期间O3污染特征和驱动因子
李源, 肖致美, 毕晓辉, 蔡子颖, 徐虹, 高璟赟, 郑乃源, 杨宁