深圳大氣中VOCs的二次有機氣溶膠生成潛勢

王扶潘,朱 喬,馮 凝,劉芮伶,黃曉鋒,何凌燕 (北京大學深圳研究生院,環境與能源學院,城市人居環境科學與技術重點實驗室,廣東 深圳 518055)

深圳大氣中VOCs的二次有機氣溶膠生成潛勢

王扶潘,朱 喬,馮 凝,劉芮伶,黃曉鋒,何凌燕*(北京大學深圳研究生院,環境與能源學院,城市人居環境科學與技術重點實驗室,廣東 深圳 518055)

使用氣溶膠生成系數法和產率法,對深圳四季5種芳香烴和異戊二烯所生成的二次有機氣溶膠(SOA)量分季節計算并比較得出:四季中除春季外,其余三季產率法計算得到的 SOA 量均高于生成系數法計算的結果,四季所生成 SOA 的平均值分別為(2.48±2.02)μg/m3和(2.10±1.21)μg/m3. 2種方法計算結果均為夏季SOA的生成量最大,秋季和冬季次之,春季最小.生成系數法計算得出人為源和天然源的貢獻分別為96%和4%,而產率法得到兩種源的貢獻分別為86%和14%.將生成系數法和產率法計算得到的夏季SOA值與實測值進行比較,發現計算值均低于實測值,分別占實測值的21%和31%.最后計算芳香烴和異戊二烯的自由基反應速率得出,與OH自由基反應是其生成SOA最主要的途徑,比例為75%,與NO3自由基和O3的反應比例分別為22%和3%.在生成SOA速度上,苯乙烯的速度最快,苯的速度最慢.

氣溶膠生成系數法；產率法；芳香烴；異戊二烯；二次有機氣溶膠

大氣中的揮發性有機物 (VOCs)如高碳烷烴、烯烴、芳香烴等由于其氧化態蒸汽壓比還原態蒸汽壓低,具有被氧化成二次有機氣溶膠(顆粒物)SOA的潛力,可與大氣中的 OH自由基、NO3自由基和O3等發生氧化反應而生成SOA[1].在北京、上海、珠江三角洲等地區的研究中均發現SOA是PM2.5的重要成分[2-4],而SOA中的苯環類(包括硝基苯、硝基酚、苯酚類)、脂肪族酸類、脂肪族硝酸酯類和羰基化合物類對人體健康和環境有重要影響[5],因此對SOA與前體物之間關系進行研究具有重要意義.

SOA的前體物包括揮發性有機物VOCs和半揮發性/中等揮發性有機物 S/IVOCs,而關于VOCs生成SOA的理論主要有氣相氧化理論、氣/固分配理論和顆粒相反應理論[6-9].研究者通過實測的VOCs濃度和煙霧箱實驗來研究SOA的生成,在實驗中將 VOCs引入煙霧箱并與其中的氧化劑反應,并加入一些顆粒物如(NH4)2SO4作為非均相成核表面,測量所生成的顆粒物,建立SOA 生成量與前體物之間的關系[10-13].Grosjean[5,14]假設VOCs生成SOA的產率是固定的,不隨環境條件而變化,根據煙霧箱實驗提出氣溶膠生成系數FAC,來反映SOA與初始VOCs濃度間的關系,并對多種 VOCs的氣溶膠生成系數進行了測定;Odum 等[15]對 FAC法進行了發展,考慮溫度、顆粒物濃度等因素的影響,在煙霧箱實驗中使用有機吸附分配的概念推測有機顆粒物的產率,對 VOCs的 SOA生成能力進行量化.國外的研究發現在計算得到的SOA中,芳香烴生成的SOA所占比例最大,如希臘雅典計算得到的SOA中 60%～90%的來自芳香烴的生成[16];美國南加州地區61%的SOA來自芳香烴氧化生成[13],在加拿大英屬哥倫比亞地區,這個比例約為80%[17];異戊二烯來主要自于天然源的排放,且從全球尺度上來看是排放量最大的揮發性有機物,對 SOA有重要的貢獻[18],因而本文主要研究芳香烴和異戊二烯的SOA生成潛勢.

Jacobson等[19]的研究結果認為 VOCs中含碳數大于6的烷烴、烯烴、芳香烴和羰基化合物等才可以形成SOA,低碳數VOCs的氧化產物飽和蒸汽壓太高而不能凝結成顆粒相,無法生成SOA,但Volkamer等[20]的研究表明乙二醛和乙炔等低碳的VOCs也具有生成SOA的能力.如前所述芳香烴和異戊二烯對SOA生成有重要的貢獻,結合國內外關于生成系數和SOA生成產率的報道,以及實驗室PTR-MS所監測的物種,本文所研究的VOCs物種包含有苯、甲苯、苯乙烯、C-8芳香烴、C-9芳香烴和異戊二烯.

本文根據 PTR-MS所測得的深圳四季VOCS的濃度,分別使用氣溶膠生成系數法和氣溶膠生成產率法對四季芳香烴和異戊二烯所生成的SOA進行計算比較分析,并將夏季計算得到的SOA值與實測值進行比較;最后計算芳香烴和異戊二烯的OH、NO3和O3自由基反應速率,探討SOA的生成方式和速度.

1 采樣與計算方法

1.1 采樣

為了對深圳四季芳香烴和異戊二烯的 SOA生成潛勢進行研究,作者在2011年8月10日至9月15日,2011年12月25日至2012年1月30日,2012年10月8日至11月11日和2013年3月1日至4月10日對VOCs濃度進行監測,分別代表深圳夏、冬、秋、春季的濃度水平,在每季采樣中剔除異常值,均保證有30d的有效數據,采樣地點位于北京大學深圳研究生院(22°35′28″N,113°58′30″E),校園及周邊多植被,無明顯人為源,該區域南面靠近留仙大道,東南面為南坪快速,附近以居民住宅區為主,兼有零星的工廠企業和商業設施.

在四季采樣中均使用 PTR-MS—質子轉移反應質譜[21](奧地利 Ionicon公司)在線對 VOCs的濃度進行測定,所監測的VOCs有19種,其中芳香烴物種有苯、甲苯、苯乙烯、C-8芳香烴和C-9芳香烴,以及異戊二烯.PTR-MS 基于VOCs與H3O+的質子轉移反應,使用質譜檢測器能快速測量生成的質子化VOCs離子(RH+),具有不需要前處理、時間分辨率高、靈敏度高等優點,可以在線測量質子親和性高于水 VOCs物種,包括烯烴、醇、醛、酮、芳香烴和有機酸等.劉芮伶[22]在2011年夏季觀測中對同一臺PTR- MS與在線GC-FID/PID和DOAS所測結果進行了比對.發現三臺儀器對甲苯的測量比對結果較好,變化趨勢一致;苯、甲苯和 C8芳香烴的線性較好(r>0.95);GC-FID/PID的測量結果稍高于 PTR-MS測量結果,差異在30%之內.因而使用PTR-MS能很好的測定芳香烴物種的濃度.

1.2 SOA的計算

由于VOCs生成SOA的具體反應過程非常復雜,研究者在進行煙霧箱實驗時,不考慮具體的反應過程,建立實測的VOCs濃度與SOA生成量之間的關系,使用參數化的方法計算不同條件下SOA的生成量[23-27].國內外研究中使用較多的方法有氣溶膠生成系數法和氣溶膠生成產率法.

1.2.1 氣溶膠生成系數法 氣溶膠生成系數FAC是SOA生成潛勢的一種表達方式,對于可以生成SOA的VOCs組分i,所生成的SOA用式(1)進行計算:

VOCs的初始濃度(VOCsi.initial)乘以相應的FAC系數就可以得到生成的SOA量.Grosjean在研究中認為VOCs生成SOA的生成系數FAC是固定的,本文使用的FAC值是Grosjean等[5,14]在大量煙霧箱實驗數據和大氣化學動力學數據的基礎上提出的.對于VOCs初始濃度的計算,考慮到受體點測得的 VOCs往往是經過氧化后的VOCst,它與排放源排出的初始濃度之間的關系可通過式(2)來表示:

式中:VOCst是實測的濃度;FVOCr是VOCs物種參與反應的分數,%,該值使用Grosjean在煙霧箱實驗中得到的數據.

對式(1)和式(2)進行整理綜合,最終得到如下的計算式(3):

該方法的優點是計算簡單,能從VOCs的排放清單或實測的環境濃度直接估算出 SOA的量,并可以反映各SOA前體物的貢獻.根據Grosjean的假設:SOA 的生成只在白天(08:00～17:00)發生,且VOCs只與OH自由基發生反應,但實際VOCs生成SOA的途徑還包括了與NO3自由基和O3的反應,因而根據FAC值結算的結果將低于實際值.

1.2.2 氣溶膠生成產率法 Pankow[7-8]和Odum等[15]研究者認為 SOA 生成產率并非固定的,其值要受到溫度、前體物、光照條件和 NOx濃度等因素的影響,結合氣/固分配理論,最后發展成了SOA生成產率理論.VOCs在大氣中的消耗是SOA生成的重要來源,因此以VOCs的消耗量作為基礎,使用式(4)[26]對SOA的生成量進行計算:

VOCs在大氣中的消耗主要是與OH自由基的反應,de Gouw[28]提出了關于VOCs消耗量的計算式(5):

式中:VOCsi,consumed為 VOCsi的消耗量;ERi為VOCsi相對于CO的排放比;[OH]Δt為OH自由基暴露量;Yi為SOA產率[25].

Shilling 等[25]、Ng 等[26]、Carlton 等[27]對很多重要的 VOCs的 SOA生成產率進行了測定,發現在低 NOx狀況下芳香烴的 SOA產率在0.3～0.4之間,高于在高 NOx狀況下對應的產率,而且與顆粒物濃度沒有關系,因而認為芳香烴氧化生成了揮發性很低的單一產物.本研究中使用的是 Ng[26]所測定的產率,且均為低 NOx濃度下VOCs的SOA產率,因而計算得到的SOA值是VOCs所能生成SOA的上限值,可以反映VOCs能生成SOA的最大潛勢.

在計算光化學齡時傾向于選擇2個 KOH值相差較大的 VOCs物種,如正丁烷/丙烷、甲苯/苯、二甲苯/乙苯、乙炔/CO,本研究中選擇甲苯/苯來計算光化學齡[29-30].將OH自由基濃度[OH]與光化學齡Δt的乘積定義為 OH自由基的暴露量[31],使用式(6)進行計算:

在式(5)計算 VOCsi的消耗量時需要知道各VOCsi相對于CO的排放比ERi、CO與OH自由基反應常數 KCO,二者可以利用式(7)進行多元非線性回歸得到[28]:

式中:[VOCsi]是測到的 VOCsi濃度;[CO]為 CO濃度與CO背景濃度的差值;ERi為VOCsi相對于 CO的排放比;Ki、KCO分別為 VOCsi和 OH自由基與CO的反應常數;[OH]Δt為OH自由基暴露量.

圖1 春季甲苯與苯濃度比值的日變化和苯與甲苯的相關性(黑色虛線為二者的初始排放比)Fig.1 The diurnal variation of the ratio of toluene and benzene in the spring and their correlation. (The dotted line shows their initial emission ratio)

2 結果與討論

2.1 芳香烴和異戊二烯的四季濃度水平

在春季觀測期間6種VOCs的平均濃度為(7.67±5.34)×10-9夏季平均濃度為(11.13±8.59)×10-9,秋季平均濃度為(7.98±5.15)×10-9,冬季濃度為(8.30±8.78)×10-9,可以看出在 2012～2013 年深圳四季之中,6種VOCs物種濃度在夏季最高,冬季和秋季次之,春季最低,而四季濃度的平均值為(8.77±6.97)×10-9.

各VOCs物種的季節變化規律如圖2所示:苯的濃度在春季和夏季相當,然后呈逐漸升高的趨勢,冬季濃度最大,為(1.14±1.01)×10-9;甲苯是四季觀測中濃度最高的物種且在夏季濃度最高,為(4.89±4.57)×10-9,秋季最低,為(3.46±3.19)×10-9.苯乙烯、C-9芳香烴和異戊二烯春季到夏季濃度增加,并在夏季達到四季濃度最大值,然后從夏季到冬季濃度逐漸降低;苯乙烯是濃度最低的VOCs,其最高濃度出現在夏季為(0.36±0.30)×10-9,其余三季相當;C-9芳香烴的最高濃度為(0.68±0.53)×10-9;異戊二烯主要來自于天然源的排放,因而在植物排放強烈的夏季出現濃度最大值為(0.75±0.59)×10-9;C-8芳香烴從春季到秋季濃度呈逐漸升高的趨勢,然后降低并在冬季出現最低濃度,是排放量第二大的物種,其最高濃度出現在秋季,為(3.24±3.14)×10-9.

圖2 深圳四季大氣VOCs濃度變化Fig.2 Seasonal variations of atmospheric VOCs in Shenzhen

2.2 基于FAC法計算得到的SOA生成量

根據 Grosjean煙霧箱實驗得到的氣溶膠生成系數FAC和FVOCr值,結合式(3)計算得到深圳四季大氣中各芳香烴和異戊二烯生成的SOA量,如表1所示.

從表1可以得出:四季中夏季生成的SOA量最大,為(2.70±1.55)μg/m3,冬季和秋季次之,分別為(1.92±1.26)μg/m3和(1.91±1.10)μg/m3,春季生成的SOA 量最少,為(1.85±0.93)μg/m3,全年 SOA 計算平均值為(2.10±1.21)μg/m3;5種芳香烴和異戊二烯所生成的 SOA與濃度呈正相關關系,即濃度越大則所生成SOA量也越大.甲苯由于濃度最高,SOA生成量最大,因而是對 SOA生成貢獻最大的物種,其次是C-8芳香烴,濃度最低的苯,對SOA的生成貢獻也最小;從全年看,6種物質對 SOA生成貢獻的平均值分別為苯3%、甲苯46%、苯乙烯4%、C-8芳香烴33%、C-9芳香烴10%和異戊二烯4%;在春季和冬季異戊二烯對SOA生成貢獻為3%,秋季為4%,夏季為7%,四季異戊二烯對SOA生成貢獻的平均值為 4%.夏季異戊二烯的貢獻最大,這與異戊二烯濃度夏季最高一致,可見天然源對 SOA的生成貢獻受季節影響較大.以芳香烴為代表的人為源排放的VOCs對SOA生成貢獻為96%,雖然沒有對a蒎烯和β蒎烯生成的SOA進行計算,會低估天然源對 SOA的貢獻,但仍然可以看出在深圳地區人為源的貢獻遠大于天然源的貢獻.

表1 FAC法計算得到的四季SOA生成量Table 1 The calculated SOA values of the four seasons by the FAC method

SOA的生成受到溫度、濕度、有機氣溶膠質量等很多因素的影響,因而用一個固定的 FAC值來表示SOA的產率會出現較大的偏差,但是通過FAC計算得到的SOA生成潛勢仍能給出一些SOA的重要信息,比如SOA生成的大致數量級,以及各前體物的相對貢獻等.

2.3 基于產率法計算得到的SOA生成量

使用Ng等[26]研究得到SOA生成產率,對深圳四季芳香烴和異戊二烯生成 SOA進行計算,結果如表2所示.

表2 產率法計算得到的四季SOA生成量Table 2 The calculated SOA values of the four seasons by the Yields method

根據產率法計算得到:四季中春季生成的SOA 最小,為(1.74±1.30)μg/m3,夏季生成的 SOA最高,為(4.00±3.56)μg/m3,秋季生成的 SOA 為(2.24±1.67)μg/m3,冬季生成 的 SOA 為 (1.94±1.54)μg/m3,雖然冬季 VOCs濃度高于秋季,但所生成的SOA較秋季低,全年SOA的計算平均值為(2.48±2.02)μg/m3;四季中,苯乙烯生成的 SOA值均為最大,是貢獻最大的物種,而苯生成 SOA量最小,是貢獻最小的物種;(3)從全年看,6種物質對SOA生成貢獻的平均值分別為:苯小于1%、甲苯 6%、苯乙烯 48%、C-8芳香烴 7%、C-9芳香烴19%和異戊二烯14%;(4)春季天然源排放的異戊二烯對SOA生成的貢獻為14%,夏季貢獻最大且為20%,秋季貢獻為17%,冬季貢獻最小且為 5%,四季貢獻的平均值為 14%,因而人為源排放的芳香烴對SOA生成平均貢獻為86%,可以看出人為源的貢獻遠大于天然源的貢獻.

在計算過程中發現,SOA計算值與濃度間并不存在一致的正相關,這是因為產率法的計算不僅要受到VOCs自身濃度的影響,還要受OH自由基暴露量、CO濃度的影響.影響因素的增多,會加大計算結果的不確定性,使得計算值和實測值間會存在較大的差距.

2.4 計算結果分析

對兩種方法計算結果進行比較可以發現:四季中除春季外,其余三季產率法計算得到的SOA值均高于FAC法計算的結果,兩種方法計算得到四季所生成SOA平均值分別為(2.10±1.21)μg/m3和(2.48±2.02)μg/m3,產率法計算值高于FAC法計算結果;兩種方法計算結果均表明夏季芳香烴和異戊二烯生成的 SOA量最大,秋季和冬季次之,春季生成SOA量最小;產率法中四季異戊二烯對SOA生成貢獻的平均值為 4%,產率法中異戊二烯貢獻的平均值為 14%,因而人為源的貢獻分別為96%和86%,兩種方法均表現為人為源的貢獻大于天然源的貢獻;兩者在評估物種對SOA生成貢獻大小上存在差異,FAC法中對SOA生成貢獻最大的物種是甲苯,而在產率法中貢獻最大的為苯乙烯,貢獻最小的均為苯;產率法對苯乙烯、C-9芳香烴和異戊二烯三種化學活性較強的VOCs生成的SOA計算值較FAC計算值高,因而產率法能更好地反映化學活性較強的 VOCs所能生成的SOA的水平.FAC法計算得到苯生成的SOA值約為產率法計算結果的5～8倍,可見產率法會低估活性較弱的VOCs生成SOA的能力.

在2011年夏季觀測中,宮照恒[33]使用 AMS(氣溶膠質譜儀)測得大氣中 SOA的量為(12.98±10.06)μg/m3,則 FAC法計算值占實測值的 18%～30%,平均值為21%,產率計算值占實測值的15%～46%,平均為 31%,FAC法的計算過程較產率法簡單,使用的是固定生成系數,誤差也更大,產率法是對 FAC法的繼承和發展,其計算結果更接近實測值.呂子峰[34]使用 FAC法計算得北京地區 VOCs對 SOA生成潛勢的貢獻 30%;Weitkamp[35]和Grieshop[36]通過煙霧箱實驗結合產率法估算出的SOA生成量僅占實測值的2%～8%和5%～30%.在國內外的研究結中計算值與實測值的比在一個較大的區間波動,約為 0.5%～50%,由于計算過程中只考慮了與 OH自由基的反應,同時缺失其他VOCs生成SOA的數據,因而絕大多數的估算結果僅能解釋不足一半的實測值.

不同的研究結果說明了參數化估算方法對SOA生成的量化有很大程度的不確定性,并且均表現為估算值比實測值低.本研究中產率法計算得到SOA值是芳香烴和異戊二烯生成SOA的上限值,但仍低于 AMS測得值,其主要原因是許多未測量的半揮發性有機物可能是SOA的重要前體物,并且很多VOCs的SOA生成系數和產率都還沒有報道,使得計算值比實際值低.

2.5 自由基反應率對SOA生成的影響

在不考慮具體的反應過程的情況下,可以通過計算得到 VOCs所生成的 SOA,從反應機理上講VOCs主要與OH自由基、NO3自由基和O3等發生反應而生成SOA,但3種自由基在SOA生成中的重要程度不盡相同,下面將使用深圳春季芳香烴和異戊二烯濃度數據計算自由基反應速率,來研究其生成SOA的主要方式和生成速度.VOCs的3種自由基反應速率可以用如下式進行計算:

式中:LOH,i、LNO3,i、LO3,i分別為 VOCsi的 OH、NO3、O3自由基反應速率;KOH,i、KNO3,i、KO3,i分別為 VOCsi與 OH、NO3和 O3自由基反應率常數[37];[VOCsi]為 VOCsi的濃度;[OH]、[NO3]和[O3][32]為 3種自由基濃度,在觀測過程中沒有對OH自由基的濃度進行測定,因而取文獻中深圳地區OH自由基24h平均濃度2.5×106molecules/cm[22].

圖3 不同VOCs的3種自由基反應速率比例Fig.3 The relative fractions of the reaction rates with the three radicals of different VOCs

從圖3可以看出:芳香烴和異戊二烯與3種自由基反應中,OH自由基反應速率最大,NO3自由基反應速率次之,O3自由基的反應速率最小,說明大氣中對其的去除主要是通過與 OH自由基反應,比例為 75%,NO3的去除作用占 22%,O3的去除作用占 3%;具體到每種物質,苯乙烯和異戊二烯由于有烯烴雙鍵的存在,使得與NO3自由基和O3的反應消耗速率增加,其余芳香烴的OH自由基反應速率均大于 97%,NO3的貢獻占極小的一部分,而基本不存在O3自由基消耗.

圖4 芳香烴的3種自由基反應速率Fig.4 The reaction rates with the three radicals of different aromatics

如圖 4所示,苯乙烯具有較高的 LOH值和LNO3值,同時也與 O3自由基發生反應,三自由基反應速率合計為 0.27×10-9/h,使得其在大氣中的化學活性最為活躍,生成SOA的速度最快;甲苯、C-8芳香烴和C-9芳香烴的LOH值較高,分別為0.18×10-9/h、0.14×10-9/h 和 0.16×10-9/h,使得三者生成 SOA 速度較快,但是受其濃度的影響,所生成的SOA量與自由基反應速率不是正相關;苯的OH自由基反應消耗速率在3種自由基反應速率中最大,但也僅約 0.01×10-9/h,可見苯的大氣化學反應活性很低,生成 SOA 的速度很慢,同時其濃度也最低,使得相同時間內所能生成的SOA量最低,對SOA生成貢獻也最小.

綜上,從反應方式上講,芳香烴和異戊二烯與OH自由基反應是其生成二次有機氣溶膠 SOA最主要的方式,如果是烯烴或含有烯烴雙鍵,則與NO3自由基和O3的反應也是其生成SOA重要的途徑;從生成速度上講,苯乙烯生成 SOA 的速度最快,甲苯、C-8芳香烴和C-9芳香烴次之,苯生成速度最慢.

3 結論

3.1 四季除春季外,其余三季產率法計算得到的SOA值均高于FAC法計算的結果,兩種方法計算得到四季所生成 SOA平均值分別為(2.10±1.21)μg/m3和(2.48±2.02)μg/m3,產率法計算值高于FAC法計算結果;夏季生成的SOA最大,秋季和冬季次之,春季最小;兩種方法計算得到異戊二烯對SOA生成貢獻的平均值分別為4%和14%,人為源的貢獻分別為96%和86%,均表現為人為源的貢獻大于天然源的貢獻.

3.2 FAC法中對 SOA生成貢獻最大的物種是甲苯,而在產率法中貢獻最大的為苯乙烯,貢獻最小的均為苯;產率法能更好的反映化學活性較強的VOCs所能生成的SOA的水平,但會低估活性較弱的VOCs物種生成SOA的能力.

3.3 FAC法計算值占實測值的21%,而產率計算值占實測值的 31%.產率法使用低 NOx濃度下SOA產率計算得到的SOA值是VOCs所能生成SOA 的上限,但仍較實際值低,說明使用的參數化方法估算SOA生成潛勢有很大程度的不確定性.

3.4 芳香烴和異戊二烯主要通過與OH自由基反應而生成SOA,但苯乙烯和異戊二烯生成SOA的途徑還包括了與O3和NO3自由基反應;在生成SOA速度上,苯乙烯的速度最快,甲苯、C-8芳香烴和C-9芳香烴次之,苯的速度最慢.

[1]唐孝炎,張遠航,邵 敏.大氣環境化學. [M]. 2版.北京:高等教育出版社, 2006

[2]Feng Y, Chen Y, Guo H, et al. Characteristics of organic and elemental carbon in PM2.5samples in Shanghai, China [J].Atmospheric Research, 2009,92(4):434-442.

[3]Lin P, Hu M, Deng Z, et al. Seasonal and diurnal variations of organic carbon in PM2.5 in Beijing and the estimation of secondary organic carbon [J]. Journal of Geophysical Research,2009,114(D2):D00G11.

[4]Hu W W, Hu M, Deng Z Q, et al. The characteristics and origins of carbonaceous aerosol at a rural site of PRD in summer of 2006[J]. Atmospheric Chemistry and Physics, 2012,12(4):1811-1822.

[5]Grosjean D. In situ organic aerosol formation during a smog episode estimated production and chemical functionality [J].Atmospheric Environment, 1992,26A:953-963.

[6]Kroll J H, Seinfeld J H. Chemistry of secondary organic aerosol:Formation and evolution of low-volatility organics in the atmosphere [J].Atmospheric Environment, 2008,42(16):3593-3624.

[7]Pankow J F. An absorption-model of gas-particle partitioning of organic-compounds in the atmosphere [J]. Atmospheric Environment, 1994,28(2):185-188.

[8]Pankow J F. An absorption-model of the gas aerosol partitioning involved in the formation of secondary organic aerosol [J].Atmospheric Environment, 1994,28(2):189-193.

[9]陳文泰,邵 敏,袁 斌,等.大氣中揮發性有機物(VOCs)對二次有機氣溶膠(SOA)生成貢獻的參數化估算 [J]. 環境科學學報, 2013,33(1): 163-172.

[10]Alfarra M R, Paulsen D, Gysel M, et al. A mass spectrometric study of secondary organic aerosols formed from the photooxidation of anthropogenic and biogenic precursors in a reaction chamber [J]. Atmospheric Chemistry and Physics, 2006,6(12):5279-5293.

[11]Hao L Q, Wang Z Y, Huang M Q, et al. Effects of seed aerosols on the growth of secondary organic aerosols from the photooxidation of toluene [J]. Journal of Environmental Sciences,2007,19(6):704-708.

[12]Robinson A L, Donahue N M, Shrivastava M K, et al. Rethinking organic aerosols: semivolatile emissions and photochemical Aging [J]. Science, 2007,315(5816):1259-1262.

[13]George I J, Abbatt J P D. Chemical evolution of secondary organic aerosol from OH-initiated heterogeneous oxidation [J].Atmospheric Chemistry and Physic, 2010,10(12):5551-5563.

[14]Grosjean D, Seinfeld JH. Parameterization of the formation potential of secondary organic aerosols [J]. Atmospheric Environment, 1989,23:1733-1747.

[15]Odum J R, Jungkamp T P W, Griffin R J ,et al.The atmospheric aerosol-forming potential of whole gasoline vapor [J]. Science,1997,276(5309):96-99.

[16]Kourtidis K, Ziomas I. Estimation of secondary organic aerosol( SOA)production from traffic emissions in the city of Athens [J]. Global Nest, 1999,1(1):33-39.

[17]Barthelmie R J, Pryor S C. Secondary organic aerosols: formation potential and ambient data [J]. Sci. Total.Environ, 1997,205(2)-3:167-178.

[18]Claeys M, Graham B, Vas G, et al. Formation of secondary organic aerosols through photo oxidation of isoprene [J]. Science,2004,303:1173–1176.

[19]Jacobson M C, Hansson H C, Noone K J, et al. Organic atmospheric aerosols: review and state of the science [J]. Rev Eophys, 2000,38(2):267-294.

[20]Volkamer R, Martini F S, Molina L T D, et al. A missing sink for gas-phase glyoxal in Mexico City: Formation of secondary organic aerosol [J]. GeophysicalResearchLetters, 2007,34(L19807),doi:1029/2007GL030752.

[21]Hansel A, Jordan A. Proton-transfer reaction massspectrometry—online trace gas-analysis at the ppb level [J].International Journal of Mass Spectrometry and Ion Processes 1995,149:609-619.

[22]劉芮伶.基于質子轉移反應譜的大氣含氧揮發性有機物氣在線測量和來源研究 [D]. 北京:北京大學碩士學位論文, 2012.

[23]DeCarlo P F, Kimmel J R, Trimborn A, et al. Field-deployable,high-resolution, time-of-flight aerosol mass spectrometer [J].Anal Chem, 2006,78(24):8281-8289

[24]Kimmel J R, Farmer D K, Cubison M J, et al. Real-time aerosol mass spectrometry with millisecond resolution [J]. Int J Mass Spectrom, 2011, 303(1):15–26.

[25]Shilling J E, Chen Q, King S M, et al. Particle mass yield in secondary organic aerosol formed by the dark ozonolysis of alpha-pinene [J]. Atmospheric Chemistry and Physics, 2008,8(7):2073-2088.

[26]Ng N L, Kroll A. Chan W, et al.Secondary organic aerosol formation from m-xylene, toluene, and benzene [I], Atmospheric Chemistry and Physics, 2007,7(14):3909-3922.

[27]Carlton A G, Wiedinmyer C, Kroll J H. A review of Secondary Organic Aerosol(SOA)formation from isoprene [J]. Atmospheric Chemistry and Physics, 1999,9(14),4987-5005.

[28]de Gouw J A. Budget of organic carbon in a polluted atmosphere:Results from the New England Air Quality Study in 2002 [J],Journal of Geophysical Research Atmospheres,2005,110(D16305),doi:10.1029/2004JD005623.

[29]McKeenS A, Liu S C. Hydrocarbon ratios an photochemical of air mass [J]. Geophysical Research Letters, 1993,20(21):2363-2366.

[30]Grosjean E, Rasmussen R A. Toxic air contaminants in Porto Alegre, Brazil [J]. Environmental Science and Technology,1999,33(12): 1970-1978.

[31]Jimenez J L. Evolution of organic aerosols in the atmosphere [J],science, 2009,326(5959),1525-1529.

[32]袁 斌.揮發性有機物(VOCs)化學轉化的量化表征及其應用研究 [b]. 北京:北京大學博士學位論文.

[33]宮照恒.基于高分辨質譜在線觀測的2011深圳大運會前后PM1化學組成與粒徑分布 [J]. 中國科學:化學, 2013,43(3):363-372.

[34]呂子峰,郝吉明,段菁春,等.北京市夏季二次有機氣溶膠生成潛勢的估算 [J]. 環境科學, 2009,30(4):969-975.

[35]Weitkamp E A, Sage A M, Pierce J R, et al. Organic aerosol formation from photochemical oxidation of diesel exhaust in a smog chamber [J]. Environmental Science and Technology, 2007,41(20):6969-6975.

[36]Grieshop A P, Logue J M, Donahue N M, et al. Laboratory Investigation of photochemical oxidation of organic aerosol from wood fires 1: measurement and simulation of organic aerosol evolution [J]. Atmospheric Chemistry and Physics, 2009,9(4):1263-1277.

[37]Atkinson R, Baulch D L, Cox R A, et al. Evaluated kinetic and photochemical data for atmospheric chemistry: Volume II; gas phase reactions of organic species [J]. Atmospheric Chemistry and Physics, 2006,6:3625-4055.

The generation potential of secondary organic aerosol of atmospheric VOCs in Shenzhen.

WANG Fu-pan, ZHU Qiao,FENG Ning, LIU Rui-ling, HUANG Xiao-feng, HE Ling-yan*(Key Laboratory for Urban Habitat Environmental Science and Technology, School of Environment and Energy, Shenzhen Graduate School, Peking University, Shenzhen 518055,China). China Environmental Science, 2014,34(10)：2449～2457

Both the fractional aerosol coefficients (FAC) and yields (Y) methods was wsed to calculate the secondary organic aerosol (SOA)production from atmospheric volatile aromatics and isoprene in ShenZhen. Throughant the year results from the Yields method were all larger than the results by the FAC method except in spring, and the total average SOA production values were (2.48±2.02)μg/m3and (2.10±1.21)μg/m3. The calculated SOA value in the summer was the largest, followed by in the autumn, winter, and spring. By the FAC method, the contribution from anthropogenic sources was 96%, while the contribution from natural sources was 4%. By the Yields method, contributions of anthropogenic and natural sources were 86% and 14%, respectively. During summer time, the calculated SOA values by the FAC and Yields methods only accounted for 21% and 31% of the measured SOA value. The OH radical reaction activities of aromatic hydrocarbons and isoprene were also calculated, and the results showed that the aromatic hydrocarbons' reaction with OH radicals was the main source of SOA with a percentage of 75%, while the percentages of NO3and O3were 22% and 3%.In terms of generating speed of SOA, styrene was found to have the fastest speed, while the speed of benzene was the slowest.

FAC method；Yields method；aromatics；isoprene；secondary organic aerosol

X513

A

1000-6923(2014)10-2449-09

2013-11-27

國家自然科學基金(U1301234;21177001);國家“973”項目(2013CB228503);深圳市科技計劃

* 責任作者, 教授, hely@pku.edu.cn

王扶潘(1987-),男,四川綿陽人,北京大學環境與能源學院碩士研究生,主要從事大氣揮發性有機物研究.