2. 中国海洋大学海洋大气相互作用与气候实验室, 物理海洋教育部重点实验室, 青岛 266100
2. Ocean-Atmosphere Interaction and Climate Laboratory, Key Laboratory of Physical Oceanography of Ministry of Education, Ocean University of China, Qingdao 266100, China
大气颗粒物是大气重要组成部分,粒子数浓度是大气颗粒物污染的重要特征之一[1]. 由于大气颗粒物污染,最近一段时间,我国中东部地区持续雾霾天气,给人们的生产和生活带来了巨大影响. 为了定量研究大气颗粒物对空气质量的影响,需要详细了解其微物理特征,包括颗粒物的数密度、 粒子谱分布等[2]. 研究表明大气颗粒物来源复杂,受气象条件的影响很大,具有很大的时空变化[3, 4, 5, 6, 7, 8, 9]. 我国大气颗粒物呈现大气复合型污染特征,不同地区不同季节大气颗粒物主要来源和相对贡献存在差异[10].
国内学者在多个城市研究了颗粒物数浓度的粒径谱分布,其与风速、 相对湿度的关系[11, 12, 13, 14],但是许多研究的采样时间较短, 除北京、 杭州外[15, 16],季节变化研究比较少. 另外关于颗粒物与气团来源、 混合层高度的相关关系研究比较少. 气象条件对大气气溶胶的影响是多种因素综合作用的结果,但在一定条件下某些特定因素起主要作用[17, 18, 19, 20, 21, 22 ,23, 24, 25, 26, 27]. 青岛是海滨城市, 内陆污染物输送和海洋排放以及局地污染都会对大气气溶胶浓度的变化产生影响,以前针对青岛地区大气颗粒物的研究主要集中在海洋因子对大气颗粒物的贡献、 气溶胶成分、 沉降通量,以及气象条件对气溶胶光学特性的影响等方面[28, 29, 30, 31, 32, 33, 34, 35],但是针对数浓度观测研究还不多. 本研究的目的是通过现场观测研究青岛地区大气颗粒物数浓度的变化特征并分析颗粒物来源和其对能见度的影响.
粒子浓度观测在中国海洋大学校园内的八关山(120°20.134′E,36°3.809′N)上进行,海拔高度大约为82.4 m. 观测点视野开阔,人为干扰因素少. 观测时间为2010年9月18日~2011年8月24日,约每周观测一次,每次观测时间固定在13:00~17:00左右. 雨雪天气时停止观测.
用便携式light house激光粒子计数器HANDHELD 5016测量最小粒径为0.5 μm的气溶胶粒子数浓度以及温湿度. 计数器具有0.5、 0.7、 1.0、 3.0、 5.0、 10.0 μm共6个粒子通道,流量为0.1 ft3 ·min-1. 每30 s采样一次,每次的采样容积为0.05 ft3. 采样同时观测和记录天气现象.
青岛气象台距离观测站点约1 km,海拔高度相近. 从气象台获取观测时间内每1 h间隔的风、 气温、 气压、 相对湿度、 水汽压和每10 min自动观测的能见度等气象要素.
利用Hysplit模式,计算大气颗粒物的后向轨迹,后推出大气颗粒物的源地. 利用观测当日14:00气温和08:00青岛站的探空资料,用干绝热法计算出14:00混合层高度. 运用统计分析方法,初步探讨大气颗粒物粒子数浓度变化与气象因子的关系.
表 1为粒子数浓度分布表. 颗粒物粒子数浓度随粒径呈指数减小,0.5-0.7 μm粒子数浓度量级达100~102个 ·cm-3,而10.0 μm以上粒子数浓度量级仅10-4~10-2 个 ·cm-3. 从多日平均来看,各粒径粒子数浓度在下午一段时间内没有显著变化,说明气溶胶粒子在较短的时间内浓度是较稳定的,但不同的观测日间相差较大. 图 1和图 2给出了颗粒物数浓度随季节的变化和月变化,结合图 1和图 2可以看出,颗粒物粒子数浓度季节变化较明显,除了0.5-0.7 μm的小粒子秋季较高以外,其他粒径都是冬春较高; 除了大于10 μm的大粒子夏季较高以外,其他粒径都是夏季最低; 0.5-3.0 μm颗粒物数浓度11月最高,9月最低; 3.0 μm以上颗粒物数浓度3月、 12月较高,6、 7、 8月较低.
![]() | 表 1 粒子数浓度统计表 Table 1 Statistic of the particle concentrations |
![]() | 图 1 颗粒物数浓度随季节的变化 Fig. 1 Seasonal variation of the particle number concentrations |
![]() | 图 2 颗粒物数浓度的月变化 Fig. 2 Monthly variation of the particle number concentrations |
与其他地区观测结果比较,春季各粒径粒子数浓度都低于广州地区[39],夏季粒子数浓度低于济南[36]、 南京[37]、 兰州[38],秋冬季0.5-1.0 μm颗粒物数浓度与上海地区[39]相当,高于北京地区[39]; 1.0-2.5 μm颗粒物数浓度较低. 可能是由于13:00以后海风发生频率较高[40],有利于大气污染物浓度的降低; 观测时段没有涵盖早晚高峰,也可能是造成颗粒物数浓度较低的原因之一.
冬春季节冷气压控制本地区,常有下沉逆温出现,大气层结较为稳定, 水平扩散减弱,因此颗粒物数浓度较高; 夏季受西太平洋副热带高压的影响,大气水平扩散增强[24],对大气颗粒数有稀释作用. 另外,冬季采暖、 春季沙尘等因素也会对大气颗粒物浓度增加有显著影响. 观测期间3月19日、 5月1日有沙尘,12月2日有雾. 沙尘一般为较粗粒径颗粒物,由外源输送而来,使颗粒物浓度增加; 雾是水汽附着在凝结核上形成,颗粒粒径较大,因此使观测期间粗颗粒浓度明显偏高. 张仁健等[32]研究了青岛大气气溶胶的质量浓度随季节的分布,发现青岛气溶胶质量浓度春季最高,冬秋季次之,夏季最低,这与本研究结果基本一致.
利用Hysplit模式对观测当天大气颗粒物进行72 h溯源分析(以观测中间时刻15:00开始往前溯源72 h,起始高度为500 m). 图 3给出了观测期间颗粒物的后向轨迹,图 4是各方向粒子平均值的统计结果,从北向开始顺时针各方向数据所占的比例分别为6.98%、 6.98%、 16.28%、 2.33%、 4.65%、 4.65%、 20.93%和37.21%. 结合图 3和图 4可以看出,从西向新疆、 甘肃一带地区和西南方向长江中下游一带输送的颗粒物数浓度较高,西北方向大粒径粒子数浓度偏高; 从东北方向及东南向海上输送的气团空气质量较好,大气颗粒物数浓度较低,其中东向大粒径粒子数浓度较低.
![]() | 图 3 观测期间气团轨迹 Fig. 3 Air mass trajectory during the observation period |
![]() | 图 4 各方向粒子数浓度平均值 Fig. 4 Mean value in every direction |
冬季颗粒物主要来源于西、 西北方向的长距离输送(西、 西北向气团来源占75%,地面风向为偏北的占66.7%),颗粒物数浓度较高,尤其是粗粒子数浓度较高; 夏季气团主要来源于偏东方向(占50%),颗粒物主要来源于海上及局地源的影响. 夏季观测的10 d中,有8 d为有雾天气,相对湿度较高,混合层高度较低,大气较稳定,暖湿气流经过冷海面形成海雾,且青岛站多处于高压入海后的东南气流影响下,有利于海雾传输至观测点. 为验证以上分析结果,选取了冬季12月2日和夏季7月10日的观测资料,利用Hysplit模式计算了当观测中间时刻(15:00)一个单位质量浓度的颗粒物扩散到青岛时前3 h 0-50 m近地层的颗粒物质量浓度(图 5). 从图 5可以看出,12月2日风向为西北风,污染物主要来自青岛西北方向,当天风向有利于颗粒物向青岛输送; 而7月10日污染物主要集中在青岛东南沿海,在偏东风的作用下颗粒物向青岛输送. 刘咸德等[41]的研究表明,冬季土壤扬尘和燃煤飞灰,硫酸盐二次颗粒物为大气颗粒物的主要来源,因此粗粒子数浓度较高; 而在夏季,海盐颗粒的浓度水平是冬季的4倍,说明颗粒物主要来自于海上. 这与本研究的分析结果吻合.
![]() | 图 5 观测前3小时扩散到青岛站的颗粒物质量浓度与青岛站观测期间的平均风速风向 Fig. 5 Particle mass concentrations 3 hours before observation and the average wind speed and direction at Qingdao station |
分析表明,大气颗粒物粒子数浓度变化与风速以及相对湿度的变化基本上呈现负相关关系,其中风速与较小粒径(0.5-3.0 μm)粒子数浓度负相关比较明显; 大气颗粒物数浓度与气温、 气压的相关关系不是很明显(图略).
将风速按照蒲福风力分级,与对应时刻各粒径段的粒子数浓度观测结果比较得到图 6; 将相对湿度按20%以下、 20%~30%、 30%~40%、 40%~50%、 50%~60%、 60%~70%、 70%~80%和80%以上分为8个等级,与对应时刻各粒径段的粒子数浓度观测结果比较得到图 7. 结合图 6和图 7可以看出,对于较小粒径粒子(0.5-0.7 μm),当风速较低(1-3级),相对湿度在40%~60%时,颗粒物数浓度最大; 当风速在5级以上时,0.5-3.0 μm粒子数浓度较高; 当相对湿度在30%以下和40%~60%时,0.5-3.0 μm粒子数浓度较高; 风速和相对湿度对较大粒径粒子(3.0 μm以上)数浓度影响不大.
![]() | 图 6 不同风速下大气颗粒物数浓度 Fig. 6 Atmospheric particle number concentrations under different wind speed |
![]() | 图 7 不同相对湿度下大气颗粒物数浓度 Fig. 7 Atmospheric particle concentrations under different relative humidity |
根据报道,当风速小于6 m ·s-1时,扩散作用占主导[33],一般情况下随着风速的增大, 单位时间内从污染源排放出来的污染物很快地扩散, 这时混入的大气量越多, 气溶胶浓度越小, 因此在其他条件不变的情况下,气溶胶浓度与风速成反比[22]; 相对湿度<60%时,水汽不是影响颗粒物浓度变化特征的第一要素,但相对湿度会促进超细颗粒物的聚积,使之生成较大粒径颗粒物[13]; 而对于更大粒径的大气颗粒物,由于聚积的作用会使其增大,质量增加,引起沉降,从而大气中的颗粒物数浓度降低. 风速与相对湿度的结论与林俊等[13]关于较粗粒径(0.2-2.4 μm)颗粒物的研究结论基本一致,但其并未对更大粒径(大于2.4 μm)颗粒物数浓度与气象因子的关系做讨论,这可能是不同地区大气颗粒物性质的差异和采样条件的限制所致.
为研究稳定度对大气颗粒物数浓度的影响,利用观测当日14:00地面温度和08:00青岛站的探空资料,用干绝热法计算出午后混合层高度,此混合层高度即为14:00混合层高度,一般为当日最大混合层高度,可以代表由于湍流混合作用而使污染物在铅直方向能充分稀释、 扩散的最大范围. 图 8绘出了大气颗粒物数浓度和混合层高度的变化,从中可以看出大气颗粒物数浓度基本上和混合层高度之间呈负相关关系,其中0.5-3.0 μm粒子数浓度与混合层高度之间的简单相关系数达到-0.27. 颗粒物数浓度较大的12月2日和2月27日,地面都存在明显的贴地逆温, 温度露点温度差较大,大气处于极稳定的状态. 这与王开燕等[23]对番禺地区秋季气溶胶浓度与气象条件的关系研究一致. 混合层高度越低,大气越稳定,使颗粒物在混合层内堆积,浓度升高.
![]() | 图 8 大气颗粒物数浓度和混合层高度的变化Fig. 8 Variation of atmospheric particle concentrations and mixed-layer height |
为定量讨论风、 温度、 压强、 相对湿度等气象要素对大气颗粒物粒子数浓度的影响,经过计算,得到了大气颗粒物粒子数浓度与各气象要素的偏相关系数(表 2). 选择风速、 气温、 相对湿度、 气压、 水汽压作为可能的因子,利用逐步回归,建立各粒径粒子数浓度与气象因子的回归方程. 回归系数列于表 3.
![]() | 表 2 偏相关系数计算1)Table 2 Partial correlation coefficient |
![]() | 表 3 回归系数计算1)Table 3 Regression coefficient |
结合表 2和表 3可以看出,风速和相对湿度与大气颗粒物的粒子数浓度基本上呈负相关,而气温与气压对大气颗粒物数浓度的相关性不太明显. 这与定性分析的结论基本一致.
温度主要对大气垂直对流产生影响. 温度较高时,大气垂直对流作用加剧,有利于大气扩散,因而一般与气溶胶浓度呈负相关[13]. 温度和各粒径粒子数浓度相关性不太明显,可能是测量时大多是晴好天气,垂直对流较弱,从而大气扩散也较弱造成的.
计算了各粒径粒子数浓度与能见度的简单相关系数分别为-0.705 0、 -0.657 1、 -0.376 7、 -0.137 1、 -0.139 5、 0.113 7,均通过了显著性水平为0.05的检验,可以看出粒子数浓度与能见度呈现较好的负相关关系,其中0.5-0.7 μm、 0.7-1.0 μm、 1.0-3.0 μm粒径段粒子数浓度与能见度的负相关关系最为明显. 图 9是不同粒径粒子数浓度与能见度的散点图,并对其做了线性和乘幂的拟合. 从中可以看出细粒子数浓度与能见度有较好的负相关关系.
![]() | 图 9 粒子数浓度与能见度散点图Fig. 9 Scatter plot of particulate number concentrations and visibility |
图 10显示了低于6 000 m能见度下气团来源与地面风向各方向日数. 总体来看,低能见度下粒子数浓度较高,气团主要来源于西向和西北向(总共占50%),地面风向为南到东南风(75%),地面风速较低,混合层高度较低. 由此可见,粒子数浓度是影响能见度的重要因子,当细粒子数浓度较高时,在西到西北向气团来源,较低风速和混合层高度条件下,容易出现低能见度现象.
![]() | 图 10 低于6 000 m低能见度下气团来源与地面风向玫瑰图Fig. 10 Rose diagram of air mass source and surface wind direction when the low visibility was below 6 000 m |
(1)青岛大气颗粒物浓度季节变化比较明显,冬春最高, 秋季次之,夏季最低.
(2)气团来源对大气颗粒物浓度有重要影响. 气团源自新疆、 甘肃及长江中下游一带的颗粒物数浓度较高,而来自于东北方向及海上的大气颗粒物数浓度较低.
(3)大气颗粒物粒子数浓度变化与风速、 相对湿度以及混合层高度的变化有较好的负相关关系,而与气温、 气压的相关关系不是很明显.
(4)在西到西北向气团来源,较低风速和混合层高度条件下,细粒子数浓度较高,容易出现低能见度现象.
[1] | 杨柳, 吴烨, 宋少洁, 等. 不同交通状况下道路边大气颗粒物数浓度粒径分布特征[J]. 环境科学, 2012, 33 (3): 694-700. |
[2] | 黄红莲, 黄印博, 饶瑞中. 合肥地区气溶胶数浓度和谱分布特征的观测研究[J]. 红外与激光工程, 2006, 35 (S3): 390-395. |
[3] | 银燕, 崔振雷, 张华, 等. 2006年中国地区大气气溶胶浓度分布特征的模拟研究[J]. 大气科学学报, 2009, 32 (5): 595-603. |
[4] | 刘毅, 周明煜. 中国近海大气气溶胶的时间和地理分布特征[J]. 海洋学报, 1999, 21 (1): 32-40. |
[5] | 何炳文, 胡振菊, 张东之, 等. 2006年9-10月常德市太阳山大气中PM10气溶胶变化特征[J]. 安徽农业科学, 2009, 37 (5): 2107-2108, 2125. |
[6] | 潘小乐, 王自发, 王喜全, 等. 2008北京奥运会前后城区黑碳气溶胶浓度的变化特征[J]. 气候与环境研究, 2010, 15 (5): 616-623. |
[7] | 周家茂, 曹军骥, 张仁健. 北京大气中PM2.5及其碳组分季节变化特征与来源[J]. 过程工程学报, 2009, 9 (2): 248-252. |
[8] | 谭建成, 林国杨, 陈政豪, 等. 利用激光雷达测量都市上空气溶胶的浓度分布[J]. 光散射学报, 2008, 20 (4): 375-378. |
[9] | Ji J H, Bae G N, Hwang J. Characteristics of aerosol charge neutralizers for highly charged particles[J]. Journal of Aerosol Science, 2004, 35 (11): 1347-1358. |
[10] | 胡敏, 唐倩, 彭剑飞, 等. 我国大气颗粒物来源及特征分析[J]. 环境与可持续发展, 2011, 36 (5): 15-19. |
[11] | 孙俊玲, 刘大锰, 扬雪. 北京市海淀区大气颗粒物污染水平及其影响因素[J]. 资源与产业, 2009, 11 (1): 96-100. |
[12] | 张涛, 陶俊, 王伯光, 等. 2008年1月广州颗粒物数浓度污染特征[J]. 中国环境监测, 2009, 25 (2): 31-34. |
[13] | 林俊, 刘卫, 李燕, 等. 大气气溶胶粒径分布特征与气象条件的相关性分析[J]. 气象与环境学报, 2009, 25 (1): 1-5. |
[14] | 徐鹏炜, 徐乔根, 蔡菊珍. 杭州城市可吸入颗粒物污染与气象条件关系的分级评价和BP神经网络评估[J]. 环境污染与防治, 2009, 31 (6): 92-95. |
[15] | 姜忠, 石广玉, 陈焕森. 北京市1998-2001年大气气溶胶粒子数浓度分析[J]. 气候与环境研究, 2003, 8 (4): 495-502. |
[16] | 谢小芳, 孙在, 付志民, 等. 杭州市大气超细颗粒数浓度谱季节性特征[J]. 环境科学, 2013, 34 (2): 434-440. |
[17] | 梁延刚, 胡文志, 杨敬基. 香港能见度、大气悬浮粒子浓度与气象条件的关系[J]. 气象学报, 2008, 66 (3): 261-269. |
[18] | 姚青, 孙玫玲, 张长春, 等. 潮州沿海大气气溶胶无机离子浓度分布与气象要素的相关分析[J]. 气象与环境学报, 2007, 23 (4): 39-42. |
[19] | 王开燕, 王雪梅, 唐建辉, 等. 番禺秋季气溶胶浓度的变化与气象条件的关系[J]. 广东气象, 2008, 30 (6): 9-12. |
[20] | 冯建军, 沈家芬, 梁任重, 等. 广州市PM10与气象要素的关系分析[J]. 中国环境监测, 2009, 25 (1): 78-82. |
[21] | 赵文慧, 宫辉力, 赵文吉, 等. 基于地统计学的北京市可吸入颗粒物时空变异性及气象因素分析[J]. 环境科学学报, 2010, 30 (11): 2154-2163. |
[22] | 申占营, 陆斌, 陈海波, 等. 降水过程中气象条件对郑州市区气溶胶浓度的影响[J]. 气象与环境科学, 2009, 32 (3): 55-58. |
[23] | 王开燕, 王雪梅, 张仁健. 气象条件对北京北郊冬季气溶胶污染的影响[J]. 过程工程学报, 2006, 6 (2): 70-74. |
[24] | 金维明. 主要气象因素对可吸入颗粒物浓度影响规律探讨[J]. 中国环境监测, 2009, 25 (4): 71-75. |
[25] | Wróbel A, Rokita E, Maenhaut W. The influence of meteorological parameters on short range transport of aerosols[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 1999, 150 (1-4): 403-408. |
[26] | Stulov E A, Plaude N O, Monakhova I A. Weather effects on the aerosol characteristics in the surface layer of the atmosphere[J]. Russian Meteorology and Hydrology, 2010, 35 (2): 94-100. |
[27] | Hussein T, Karppinen A, Kukkonen J, et al. Meteorological dependence of size-fractionated number concentrations of urban aerosol particles[J]. Atmospheric Environment, 2006, 40 (8): 1427-1440. |
[28] | Manton M J, Ayers G P. On the number concentration of aerosols in towns[J]. Boundary-Layer Meteorology, 1982, 22 (2): 171-181. |
[29] | 徐新华, 姚荣奎, 李金龙. 青岛地区大气气溶胶海洋因子贡献研究[J]. 海洋环境科学, 1997, 16 (2): 55-59. |
[30] | 陈静静, 盛立芳, 郑元鑫. 2006年春季青岛气溶胶光学厚度与气溶胶指数、能见度的关系分析[J]. 中国海洋大学学报 (自然科学版), 2010, 40 (1): 10-16. |
[31] | 邱明燕, 盛立芳, 房岩松, 等. 气象条件对青岛地区气溶胶光学特性的影响[J]. 中国海洋大学报 (自然科学版), 2004, 34 (6): 925-930. |
[32] | 张仁健, 浦一芬, 徐永福, 等. 青岛大气气溶胶的浓度分布和干沉降的观测研究[J]. 气候与环境研究, 2004, 9 (2): 390-395. |
[33] | 盛立芳, 耿敏, 王园香, 等. 2002年春季沙尘暴对青岛大气气溶胶的影响[J]. 环境科学研究, 2003, 16 (5): 11-13. |
[34] | Smirnov A, Villevalde Y, O'Neill N T, et al. Aerosol optical depth over the oceans: Analysis in terms of synoptic air mass types[J]. Journal of Geophysical Research: Atmospheres, 1995, 100 (D8): 16639-16650. |
[35] | 邱明燕. 青岛地区气溶胶特性的初步研究[D]. 青岛: 中国海洋大学, 2004.12-35. |
[36] | 高健, 王进, 程淑会, 等. 济南夏季大气颗粒物粒径分布特征及来源机理分析[J]. 中国科学院研究生院学报, 2007, 24 (5): 680-687. |
[37] | 王红磊, 朱彬, 沈利娟, 等. 南京市夏季大气气溶胶新粒子生成事件分析[J]. 环境科学, 2012, 33 (3): 701-710. |
[38] | 赵素平, 余晔, 陈晋北, 等. 兰州市夏秋季颗粒物谱分布特征研究[J]. 环境科学, 2012, 33 (3): 687-693. |
[39] | 张涛, 陶俊, 王伯光, 等. 广州市春季大气颗粒物的粒径分布及能见度研究[J]. 中国科学院研究生院学报, 2010, 27 (3): 331-337. |
[40] | 吴增茂, 龙宝森. 青岛局地风特征的分析[J]. 海洋湖沼通报, 1993, 1 (1): 16-22. |
[41] | 刘咸德, 封跃鹏, 贾红, 等. 青岛市大气颗粒物来源的定量解析——化学质量平衡方法[J]. 环境科学研究, 1998, 11 (5): 51-54. |