2. 大气颗粒物监测技术北京市重点实验室, 北京 100048
2. Beijing Key Laboratory of Atmospheric Particulate Monitoring Technology, Beijing 100048, China
O3是O2的同素异形体, 是天然大气的重要微量成分, 约90%的O3存在于平流层, 仅有10%左右的O3分布在对流层中[1].平流层中的O3能够吸收紫外线, 保护人类与生态环境, 但是如果对流层大气中O3浓度增高, 则会对人体的健康造成危害, 如加重呼吸疾病、损害肺功能等[2].另外, 作为一种强氧化剂, O3还在对流层许多化学过程中起着重要作用, 是光化学烟雾的主要标识物[3].城市中O3的大部分是由NOx、CO和VOCs等前体物在合适的气象条件下反应生成的[4], 近年来随着经济的发展、人口的增长以及机动车保有量的增加, 北京市O3污染问题也愈发严重[5].
近年来, 部分学者对北京市O3历史长期变化趋势开展过研究. Ding等[6]对北京地区O3的探空数据进行了分析, 发现1995~2005年间对流层下部的O3浓度以每年2%的速率增长. Tang等[7]在北京地区的监测发现, 2001~2006年间NOx浓度以每年(3.9±0.5)×10-9体积分数的速率下降, 但O3体积分数以每年(1.1±0.5)×10-9的速率上升. Wang等[8]在北京地区的研究发现, 2002~2010年间对流层O3柱浓度以每年4.6%的速率增长. Zhang等[9]在北京的监测发现, 2005~2011年间NOx和VOCs浓度分别以每年1.4×10-9体积分数和1.6×10-9体积分数的速率下降, 但白天平均O3浓度却以每年2.6×10-9体积分数的速率增长.
虽然自2013年新空气质量标准发布以来, O3浓度才纳入国家考核, 但北京市自2004年即开始对O3浓度进行了系统性监测, 本研究对2006~2015年间北京市不同站点O3浓度进行分析, 摸清其浓度水平及变化趋势, 以期为相关领域研究以及北京市O3污染防治提供科学依据.
1 材料与方法目前, 北京市环保局对外实时发布35个监测站点的空气质量, 按照监测功能的不同分为城区环境评价点、郊区环境评价点、对照点及区域点以及交通污染监控点.选择其中的定陵站(北部清洁点)、官园(城区环境点)、琉璃河(南部边界点)和前门(交通点)进行分析, 气象数据来自于北京观象台, 如图 1. O3分析仪采用Thermo Fisher 49C紫外光度法分析仪, 原理为O3分子吸收波长为254 nm的紫外光, 该波段紫外光被吸收的程度直接与O3的体积分数相关, 根据检测样品通过紫外光时被吸收的程度来计算出O3体积分数.分析仪最低检测限:1×10-9(体积分数); 精度:1×10-9(体积分数); 零漂:0.4%/24 h; 跨漂:±1%/24 h, ±2%/7 d. CO分析仪采用Thermo Fisher 48C气体过滤相关法分析仪, 最低检测限:40×10-9(体积分数); 零漂:小于100×10-9/24 h; 跨漂:±1%/24 h.在监测的质保质控方面, O3标准每年与中国计量科学研究院O3一级标准NIST SRP41进行比对和标准溯源, 并对监测站O3进行标准传递. NOx监测采用Thermo Fisher 42C化学发光NO-NO2-NOx分析仪, 其原理是NO与O3发生化学反应时产生激发态的NO2分子, 当激发态的NO2分子返回基态时发出一定波长的光, 所发出光的强度与NO的体积分数呈线性正相关.检测NO2时先将NO2通过钼转换器转换成NO, 然后再通过化学发光反应进行定量分析.该分析仪最低检测限:0.05×10-9(体积分数); 零漂:小于0.025×10-9/24 h; 跨漂:±1%/24 h.每2 d进行1次零点和跨点校准, 每1个月进行1次多点校准, 每3个月进行一次精密度审核.各操作流程严格按照文献[10]进行.
![]() |
图 1 监测站点位置示意 Fig. 1 Locations of the monitoring sites |
在进行数据处理时, 若每天中缺失数据多于4个, 则该天作为无效天数处理; 若每月中缺失天数超过6 d, 则该月作为无效月处理, 若每年中缺失月份超过3个, 则该年作为无效年处理.
2 结果与讨论 2.1 O3年际变化趋势表 1显示了2006~2015年各站点O3监测数据的有效率, 可以看出, 10年间各监测站点O3小时数据有效率基本保持在90%以上, 数据质量较高, 可以满足本研究的需求. 图 2显示了2006~2015年各站点O3平均浓度、小时标准差及拟合趋势线.从10年平均浓度来看, 清洁点定陵(65.2 μg·m-3)最高, 交通点前门浓度最低(40.4 μg·m-3), 官园(49.6 μg·m-3)跟琉璃河(53.4 μg·m-3)较为接近.从小时标准差来看, 则是官园(56.4 μg·m-3)和琉璃河(55.5 μg·m-3)最高, 其次为定陵(52.1 μg·m-3)和前门(46.5 μg·m-3), 即官园和琉璃河虽然O3浓度水平低, 但O3变化幅度较大.从年际变化来看, 各站点也表现出不同的规律, 定陵和琉璃河呈较为明显的波浪型曲线, 定陵2010~2013年间浓度较低, 为曲线的波谷, 而琉璃河在该段时间内浓度较高, 为曲线的波峰.官园和前门呈锯齿型波动变化.从各站点的10年线性拟合趋势来看, 定陵整体O3浓度呈下降趋势[0.5 μg·(m3·a)-1], 而官园[0.9 μg·(m3·a)-1]、琉璃河[0.3 μg·(m3·a)-1]和前门[0.3 μg·(m3·a)-1]均呈上升趋势.定陵站O3浓度的下降可能与局地气象因素影响以及近年来周边机动车保有量的增加有关, 较高的NO浓度滴定反应了O3[11].
![]() |
表 1 各站点O3监测数据个数及有效率 Table 1 Values and the effective rates of O3 concentrations at the monitoring sites |
![]() |
(a) O3平均浓度、小时标准差; (b) O3趋势线 图 2 2006~2015年各站点O3平均浓度、小时标准差及趋势线 Fig. 2 Average concentrations, hourly standard deviations, and trend lines for O3 at different sites from 2006-2015 |
近地面O3主要通过前体物的复杂的化学反应生成, 包括自由基的生成、传递和终止反应等[12~16], 反应速率受温度、光照和太阳辐射等气象条件影响.夏季的高温度、强光照和太阳辐射有利于O3的生成, 因此O3浓度最高. 图 3显示了2006~2015年各站点O3月均浓度最高值和最低值, 可以看出, 浓度最高值均出现在6~8月, 出现频次最高的为7月(17次), 平均月均浓度为99.8 μg·m-3; 最低值均出现在11、12月和1、2月, 出现频次最高的为1月(14次), 平均月均浓度为16.6 μg·m-3, 是最高月均浓度的16.6%.分不同站点来看, 琉璃河和定陵最高月均值分布较为均匀, 6~8月均有; 前门和官园最高月均值基本分布在7、8月.各站点最低月均值分布相对较为均匀.
![]() |
图 3 2006~2015年各站点O3月均浓度最高值和最低值 Fig. 3 Highest and lowest monthly average concentrations of O3 at different sites from 2006-2015 |
O3日变化曲线呈现规律的单峰型曲线, 图 4显示了2006~2015年各站点O3日变化最高值和最低值.可以看出, 各站点日变化峰值基本集中在15:00~18:00, 谷值基本集中在07:00~09:00.分站点来看, 处于夏季上风向的琉璃河站峰值出现的时间早于下风向的定陵和官园站, 体现了夏季O3区域传输的特征.从谷值来看, 官园、琉璃河和前门基本维持20 μg·m-3左右或以下, 而定陵谷值浓度明显偏高, 为40 μg·m-3左右, 表明北部清洁地区NO浓度相比于城区和南部地区仍较低, 夜间对O3的滴定作用较弱.另外, 各站点共有的特征是, 近年与前几年相比, O3浓度出现的时间明显提前, 近3年峰值均在15:00~16:00出现, 提前了1~2 h.朱彬等[17]、汪伟峰等[18]在观测中发现的现象一致, 也是近年来北京市大气氧化效率逐渐增强的验证.
![]() |
图 4 2006~2015年各站点O3日变化最高值和最低值 Fig. 4 Highest and lowest diurnal variation in the concentrations of O3 at different sites from 2006-2015 |
2012年国家颁布新空气质量标准[19]以后, 增加了O3-8h浓度限值, 本研究根据《环境空气质量指数(AQI)技术规定》(HJ 633-2012)中的分级办法, 统计了各站点2013~2015年各级别天数, 如表 2.从超标率来看, 交通站前门最低, 城区站官园最高, 从重污染率来看, 清洁点定陵最高, 前门站近3年未出现O3重污染现象.从年际变化来看, 前门站和官园站O3超标天数呈逐年增加趋势, 琉璃河和定陵站则呈波动规律.近3年琉璃河和官园站的重污染天数均稳定在1~3 d之间.值得注意的是, 2015年定陵站O3重污染天数达到11 d, 比2013年增加了10 d, 表明近年来夏季北京下风向山区的O3重污染状况愈发严重.
![]() |
表 2 2013~2015年各站点O3浓度各级别天数/d Table 2 Number of days reaching each level of O3 concentration at the sites from 2013-2015/d |
图 5显示了2008~2015年官园站超标日和非超标日平均O3-8h浓度, 可以看出, 由于冬季O3浓度较低, 所以非超标日O3-8h浓度明显低于国家二级标准限值(160 μg·m-3), 而超标日O3-8h浓度明显高于二级标准, 表明夏季北京市O3超标现象严重. 2011年开始, 超标日和非超标日O3-8h浓度基本呈现上升的趋势.
![]() |
图 5 2008~2015年官园站超标日和非超标日平均O3-8h浓度 Fig. 5 Eight-hour average O3 concentrations at the Guanyuan site from 2008-2015 |
图 6显示了2013~2015年各站点NO2浓度和O3日均浓度的散点图, 可以看出, NO2浓度水平前门站最高, 其次为官园站和琉璃河站, 定陵站浓度水平最低.从相关性来看, 定陵站O3浓度与NO2浓度呈正相关, 其余站点两者浓度均呈负相关, 暗示定陵站O3生成的前体物控制区可能为NO2控制区, 而其他站点为VOCs控制区, 这与徐峻等[20]和聂滕等[21]的研究结论较为一致.
![]() |
图 6 2013~2015年夏季各站点NO2和O3日均值散点图和拟合直线 Fig. 6 Scatter plots and best fit lines for the daily average concentrations of NO2 and O3 at the sites in summer from 2013-2015 |
为进行相关性分析, 首先对各组数据分别进行K-S(Kolmogorov-Smirnov)统计检验, 得出所有数据均不服从正态分布, 因此对各组数据进行Spearman相关分析. 表 3显示了2008~2015年官园站O3浓度与各气象要素的相关性系数, 可以看出, O3浓度与温度呈正相关关系, 与湿度和气压呈负相关关系.仅从相关性分析来看, 温度对O3浓度的影响最大, 其次是气压和湿度.从年际变化来看, 近3年O3浓度与各气象要素的相关性呈加强趋势, 表明北京市O3污染受气象条件的影响更为显著.
![]() |
表 3 2006~2015年间O3浓度与各气象要素的相关性系数1) Table 3 Correlation coefficients among O3 concentrations and meteorological factors from 2006-2015 |
图 7显示了2008~2015年官园站超标日和达标日气象要素水平, 可以看出, 超标日风速水平整体低于达标日, 8年平均风速低于5.4%.超标日湿度水平整体高于达标日, 8年平均湿度高出14.6%.超标日的温度水平明显高于达标日, 8年平均温度高出172.1%.超标日气压水平整体低于达标日, 8年平均气压低1.2%.
![]() |
图 7 2008~2015年官园站超标日和达标日气象要素水平 Fig. 7 Meteorological factors at the Guanyuan site from 2008-2015 |
图 8显示了2008~2015年5~8月份官园站O3日最大8 h滑动平均浓度与日最高气温和平均相对湿度的关系, 同样可以看出, O3浓度与日最高温度呈现显著的正相关关系, Pearson相关系数为0.50, 线性拟合方程为y=6.3x-25.0, O3浓度与相对湿度则整体呈现负相关关系.当日最高温度超过30℃, 相对湿度介于30%~70%之间时, 北京市O3日最大8h滑动平均浓度超过200 μg·m-3的概率较高, 空气质量级别会达到轻度至中度污染的级别.
![]() |
图 8 2008~2015年的5~8月O3日最大8 h滑动平均浓度与日最高气温和平均相对湿度的散点图 Fig. 8 Scatter plot of the 8 h O3 concentrations and daily maximum temperatureand average relative humidity during May and August, 2008-2015 |
(1) 2006~2015年北京市定陵站点O3平均浓度为65.2 μg·m-3, 为4个站点最高; 其他依次为琉璃河(53.4 μg·m-3)、官园(49.6 μg·m-3)和交通站前门(40.4 μg·m-3).从变化趋势来看, 定陵整体O3浓度呈下降趋势[0.5 μg·(m3·a)-1], 而官园[0.9 μg·(m3·a)-1]、琉璃河[0.3 μg·(m3·a)-1]和前门[0.3 μg·(m3·a)-1]均呈上升趋势.
(2) 从月变化来看, 浓度最高值均出现在6~8月, 出现频次最高的为7月(17次), 平均月均浓度为99.8 μg·m-3; 最低值均出现在11、12月和1、2月, 出现频次最高的为1月(14次), 平均月均浓度为16.6 μg·m-3, 是最高月均浓度的16.6%.日变化规律来看, 各站点日变化峰值基本集中在15:00~18:00, 谷值基本集中在07:00~09:00.近年与前几年相比, O3浓度出现的时间明显提前, 近3年峰值均在15:00~16:00出现, 提前了1~2 h.
(3) 近3年超标情况来看, 交通站前门最低, 城区站官园最高, 从重污染率来看, 清洁点定陵最高, 前门站近3年未出现O3重污染现象. 2015年定陵站O3重污染天数达到11d, 比2013年增加了10d, 表明近年来夏季北京下风向山区的O3重污染状况愈发严重.
(4) 与前体物的相关性分析表明, 定陵站O3浓度与NO2浓度呈正相关, 其余站点两者浓度均呈负相关, 暗示定陵站O3生成的前体物控制区可能为NO2控制区, 而其他站点为VOCs控制区.与气象要素的相关性分析表明, O3浓度与温度呈正相关关系, 与湿度和气压呈负相关关系, 温度对O3浓度的影响最大, 其次是气压和湿度.当日最高温度超过30℃, 相对湿度介于30%至70%之间时, 北京市O3日最大8 h滑动平均浓度超过200 μg·m-3的概率较高, 空气质量级别会达到轻度至中度污染的级别.
致谢: 感谢监测中心子站维护工作人员的辛苦和努力工作.[1] | Pan X, Kanaya Y, Tanimoto H, et al. Examining the major contributors of ozone pollution in a rural area of the Yangtze river delta region during harvest season[J]. Atmospheric Chemistry and Physics, 2015, 15(11): 6101-6111. DOI:10.5194/acp-15-6101-2015 |
[2] | Lehman J, Swinton K, Bortnick S, et al. Spatio-temporal characterization of tropospheric ozone across the eastern United States[J]. Atmospheric Environment, 2004, 38(26): 4357-4369. DOI:10.1016/j.atmosenv.2004.03.069 |
[3] | Jiang X Y, Wiedinmyer C, Carlton A G. Aerosols from fires:an examination of the effects on ozone photochemistry in the western United States[J]. EnvironmentalScience &Technology, 2012, 46(21): 11878-11886. |
[4] | Wang Y G, Hopke P K, Xia X Y, et al. Source apportionment of airborne particulate matter using inorganic and organic species as tracers[J]. Atmospheric Environment, 2012, 55: 525-532. DOI:10.1016/j.atmosenv.2012.03.073 |
[5] |
王占山, 李云婷, 陈添, 等. 北京城区臭氧日变化特征及与前体物的相关性分析[J]. 中国环境科学, 2014, 34(12): 3001-3008. Wang Z S, Li Y T, Chen T, et al. Analysis on diurnal variation characteristics of ozone and correlations with its precursors in urban atmosphere of Beijing[J]. China Environmental Science, 2014, 34(12): 3001-3008. |
[6] | Ding A J, Wang T, Thouret V, et al. Tropospheric ozone climatology over Beijing:analysis of aircraft data from the MOZAIC program[J]. Atmospheric Chemistry and Physics, 2008, 8(1): 1-13. DOI:10.5194/acp-8-1-2008 |
[7] | Tang G, Li X, Wang Y, et al. Surface ozone trend details and interpretations in Beijing, 2001-2006[J]. Atmospheric Chemistry and Physics, 2009, 9(22): 8813-8823. DOI:10.5194/acp-9-8813-2009 |
[8] | Wang Y, Konopka P, Liu Y, et al. Tropospheric ozone trend over Beijing from 2002-2010:ozonesonde measurements and modeling analysis[J]. Atmospheric Chemistry and Physics, 2012, 12(18): 8389-8399. DOI:10.5194/acp-12-8389-2012 |
[9] | Zhang Q, Yuan B, Shao M, et al. Variations of ground-level O3 and its precursors in Beijing in summertime between 2005 and 2011[J]. Atmospheric Chemistry and Physics, 2014, 14(12): 6089-6101. DOI:10.5194/acp-14-6089-2014 |
[10] |
HJ/T 193-2005环境空气质量自动监测技术规范[S]. HJ/T 193-2005Automated method for ambient air quality monitoring[S]. |
[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] |
陆克定, 张远航, 苏杭, 等. 珠江三角洲夏季臭氧区域污染及其控制因素分析[J]. 中国科学(化学), 2010, 53(4): 651-663. Lu K D, Zhang Y H, Su H, et al. Regional ozone pollution and key controlling factors of photochemical ozone production in Pearl river delta during summer time[J]. Science China Chemistry, 2010, 53(3): 651-663. |
[13] |
王占山, 张大伟, 李云婷, 等. 北京市夏季不同O3和PM2.5污染状况研究[J]. 环境科学, 2016, 37(3): 807-815. Wang Z S, Zhang D W, Li Y T, et al. Different air pollution situations of O3 and PM2.5 during summer in Beijing[J]. Environmental Science, 2016, 37(3): 807-815. |
[14] |
邵平, 安俊琳, 杨辉, 等. 南京北郊夏季近地层臭氧及其前体物体积分数变化特征[J]. 环境科学, 2014, 35(11): 4031-4043. Shao P, An J L, Yang H, et al. Variation characteristics of surface ozone and its precursors during summertime in Nanjing northern suburb[J]. Environmental Science, 2014, 35(11): 4031-4043. |
[15] |
杨俊益, 辛金元, 吉东生, 等. 2008~2011年夏季京津冀区域背景大气污染变化分析[J]. 环境科学, 2012, 33(11): 3693-3704. Yang J Y, Xin J Y, Ji D S, et al. Variation analysis of background atmospheric pollutants in North China during the summer of 2008 to 2011[J]. Environmental Science, 2012, 33(11): 3693-3704. |
[16] |
王占山, 李云婷, 陈添, 等. 2013年北京市PM2.5的时空分布[J]. 地理学报, 2015, 70(1): 110-120. Wang Z S, Li Y T, Chen T, et al. Spatial-temporal characteristics of PM2.5 in Beijing in 2013[J]. Acta GeographicaSinica, 2015, 70(1): 110-120. |
[17] |
朱彬, 安俊岭, 王自发, 等. 光化学臭氧日变化特征与其前体物关系的探讨[J]. 南京气象学院学报, 2006, 29(6): 744-749. Zhu B, An J L, Wang Z F, et al. Relations of diurnal variations of photochemical ozone to its precursors[J]. Journal of Nanjing Institute of Meteorology, 2006, 29(6): 744-749. |
[18] |
汪伟峰, 朱彬, 王跃思, 等. 一次光化学污染过程中O3和NOx的垂直梯度观测研究[J]. 环境工程学报, 2010, 4(2): 411-416. Wang W F, Zhu B, Wang Y S, et al. Studies on gradient observation of O3 and NOx in an atmospheric photochemical pollution[J]. Chinese Journal of Environmental Engineering, 2010, 4(2): 411-416. |
[19] |
GB 3095-2012环境空气质量标准[S]. GB 3095-2012Ambient air quality standards[S]. |
[20] |
徐峻, 张远航. 北京市区夏季O3生成过程分析[J]. 环境科学学报, 2006, 26(6): 973-980. Xu J, Zhang Y H. Process analysis of O3 formation in summer at Beijing[J]. Acta Scientiae Circumstantiae, 2006, 26(6): 973-980. |
[21] |
聂滕, 李璇, 王雪松, 等. 北京市夏季臭氧前体物控制区的分布特征[J]. 北京大学学报(自然科学版), 2014, 50(3): 557-564. Nie T, Li X, Wang X S, et al. Characteristics of the spatial distributions of ozone-precursor sensitivity regimes in summer over Beijing[J]. Acta ScientiarumNaturaliumUniversitatisPekinensis, 2014, 50(3): 557-564. |