張 薔,李令軍*,趙文慧,劉保獻,3,鹿海峰,秦麗歡,陳韻如,李 琪

北京森林BVOCs排放現(xiàn)狀及動態(tài)變化特征分析

張 薔1,2,李令軍1,2*,趙文慧1,2,劉保獻1,2,3,鹿海峰1,2,秦麗歡1,2,陳韻如1,2,李 琪1,2

(1.北京市生態(tài)環(huán)境監(jiān)測中心,北京 100048；2.大氣顆粒物監(jiān)測技術(shù)北京市重點實驗室,北京 100048；3.清華大學(xué)環(huán)境學(xué)院,北京 100084)

以北京地區(qū)森林植被為研究對象,基于森林資源清查蓄積資料和逐小時氣象數(shù)據(jù),采用光溫影響模型對2000～2020年北京森林BVOCs排放量進行估算,并分析其對空氣質(zhì)量的影響.結(jié)果顯示,2020年北京森林BVOCs排放量為39.57×109g C,異戊二烯、單萜烯和OVOCs分別占72.19%、17.48%和10.32%,楊樹、櫟樹等闊葉樹是主要的異戊二烯排放源,油松等針葉樹是主要的單萜烯排放源.2000～2020年森林BVOCs排放量從20.30×109g C/a 增加到39.57×109g C/a,年平均增長率4.75%;BVOCs排放量的變化表現(xiàn)出明顯階段性特征,2000～2010年增長緩慢,2010～2020年出現(xiàn)大幅上升.20年間異戊二烯所占比重呈下降趨勢,單萜烯和OVOCs所占比重則呈上升趨勢;楊樹對BVOCs排放量的貢獻逐漸降低,櫟樹和其他闊葉樹的貢獻明顯增加,北京新增森林更加注重物種多樣化.2000～2020年,BVOCs的O3生成潛勢從181.76×109g增加到331.07×109g,異戊二烯占92.70%,是主要的貢獻者;SOA 生成潛勢從1.11×109g增加到2.65×109g,單萜烯和異戊二烯分別占75.40%和24.60%.O3生成潛勢最大的樹種是楊樹,SOA生成潛勢最大的樹種是油松.森林BVOCs排放在夏季對O3污染的貢獻最大,未來綠化中應(yīng)考慮優(yōu)化樹種組成.

蓄積量；異戊二烯；單萜烯；臭氧；二次有機氣溶膠；揮發(fā)性有機物

近年來,北京市大氣細(xì)顆粒物濃度顯著下降,但是臭氧(O3)污染逐步凸顯,成為影響空氣質(zhì)量的關(guān)鍵因素,開展細(xì)顆粒物和O3的協(xié)同控制是當(dāng)前大氣污染防治工作的當(dāng)務(wù)之急[1].揮發(fā)性有機化合物(VOCs)是O3和二次有機氣溶膠(SOA)的重要前體物,可以與大氣中的氮氧化物和羥基自由基等氧化劑發(fā)生化學(xué)反應(yīng),加劇環(huán)境空氣的污染[2].作為VOCs的主要來源之一,天然源揮發(fā)性有機化合物(BVOCs)排放量大、且化學(xué)活性高,占全球VOCs總排放量的90%以上,對O3和SOA的形成具有重要影響[2-4].因此,建立一個準(zhǔn)確的BVOCs排放清單,對評估其對大氣環(huán)境的影響、制定有效的空氣污染防治措施至關(guān)重要.

自20世紀(jì)90年代以來,科學(xué)家使用不斷更新改進的生物模型對不同空間尺度的BVOCs排放量進行了估算.Graedel等[5]和Guenther等[2,4]編制了全球陸地生態(tài)系統(tǒng)VOCs排放清單,研究人員[6-8]使用不同空間分辨率數(shù)據(jù)對中國的BVOCs排放量及分布特征進行了研究.在北京地區(qū),研究人員基于不同的排放模型、土地覆蓋和氣象數(shù)據(jù),對全市植被的BVOCs排放進行了測算,如張倩等[9]采用GLOBEIS模型和遙感解譯資料估算了北京BVOCs排放通量; Ren等[10]使用改進的G95算法分別計算了北京天然林和城市綠色空間的BVOCs排放量;謝揚飏等[11]則基于園林綠化普查資料建立了北京市園林綠地BVOCs排放清單.此外,研究人員進一步分析了BVOCs對空氣質(zhì)量的影響,估算了北京地區(qū)異戊二烯對O3生成的貢獻[12-14],對異戊二烯和單萜的O3和SOA生成潛勢進行了測算[10,15-16].但是,目前對北京地區(qū)的研究大多是基于植被面積來估算葉生物量,不同研究之間的結(jié)果具有較大差異(近10倍),現(xiàn)有研究較少以具體樹種為對象、基于蓄積量進行BVOCs排放測算.事實上,研究顯示樹木的蓄積量與生物量具有顯著的相關(guān)關(guān)系[17],通過生物量計算公式可實現(xiàn)二者間的換算[18].

此外,已有的研究主要集中于BVOCs的排放現(xiàn)狀,對于BVOCs排放的動態(tài)變化特征關(guān)注較少.植被覆蓋和氣候的改變會引起B(yǎng)VOCs排放發(fā)生顯著的年際變化[1,19].近年來,為改善生態(tài)環(huán)境,北京市實施了百萬畝造林、京津風(fēng)沙源治理、綠化隔離帶建設(shè)等一系列的生態(tài)恢復(fù)工程,北京地區(qū)的森林面積和蓄積量不斷增加.據(jù)統(tǒng)計,2006～2020年北京森林覆蓋率從35.9%增長到44.4%,森林蓄積量增加了1225萬m3[20].研究BVOCs排放的變化趨勢,將有助于了解植被覆蓋快速變化背景下BVOCs排放對大氣污染的影響.因此,本研究基于Guenther等[2,21]提出的光溫影響模型,利用北京市森林資源清查蓄積資料和逐小時氣象數(shù)據(jù),建立2000～2020年北京市森林BVOCs排放清單,分析其排放現(xiàn)狀及變化趨勢,并量化BVOCs對區(qū)域空氣質(zhì)量的影響,為制定有效的污染防治措施提供科學(xué)支撐.

1 材料與方法

1.1 研究區(qū)概況

北京(39°26′～41°04′N,115°25′～117°30′E)地處華北平原北部,背靠燕山,全市總面積1.64×104km2.北京屬于典型的北溫帶半濕潤大陸性季風(fēng)氣候,四季分明,夏季炎熱多雨,冬季寒冷干燥.多年平均氣溫12.3℃,1月份氣溫最低,平均-3.7℃,7月份氣溫最高,平均26.2℃.多年平均降水量571.9mm,全年降水量約80%集中在6～8月份[22].北京地區(qū)植被資源豐富,主要的地帶性植被為暖溫帶闊葉林和針闊混交林.根據(jù)北京市園林綠化局公布的2020年森林資源調(diào)查結(jié)果,全市林地面積11299.80km2,林木綠化率為62.50%,其中森林面積8483.14km2,森林蓄積量2520.67萬m3,森林覆蓋率達(dá)44.40%[20].

1.2 研究方法

1.2.1 BVOCs排放計算模型 植物排放的VOCs種類繁多,其中研究較多的是異戊二烯和單萜烯,這是由于二者排放量大且化學(xué)活性強[23].本研究將北京森林源BVOCs劃分為異戊二烯、單萜烯和其他VOCs(OVOCs)3類,分樹種分別進行估算.植物排放VOCs的速率主要受溫度和光照的影響,本研究采用Guenther等[2,21]提出的光溫影響模型進行排放量估算.其中,異戊二烯排放受溫度和太陽輻射的影響,采用葉綠素(CHL)算法;單萜烯和OVOCs排放主要受到溫度的影響,采用防御專屬組織(DST)算法.計算公式如下:

ISOP=[][][γPγTγS][] (1)

MON/OVOC=[][][γTγS][] (2)

式中:ISOP、MON和OVOC分別代表異戊二烯、單萜烯和OVOC的排放量, g;為標(biāo)準(zhǔn)狀態(tài)下該樹種的基本排放速率, μg C/(g·h);為葉生物量, g;P、T和S為環(huán)境校正因子,分別代表光合有效輻射、溫度和季節(jié)的影響;為逸出效率,表示排放總量中逸入樹冠上方大氣中的部分所占比例.

1.2.2 O3和SOA生成潛勢估算模型 通過估算BVOCs的O3和SOA生成潛勢,表征其對空氣質(zhì)量的潛在影響.O3生成潛勢(OFP)代表VOCs物種在最佳反應(yīng)條件下對O3生成的最大貢獻,采用“增量反應(yīng)”法進行計算,即每加入一單位的VOC所形成的O3的量.SOA生成潛勢使用“氣溶膠生成系數(shù)(FAC)”法進行量化.公式如下:

OFP= VOC×MIR(3)

SOA=VOC×FAC(4)

式中: OFP, VOCs組分的OFP值, g; VOC, VOCs組分的排放量, g; MIR, VOCs組分的最大增量反應(yīng)活性.本文MIR的取值采用Carter[24]的研究成果,其中異戊二烯為10.61g/g,單萜烯為4.04g/g;最新的研究表明,該MIR值的應(yīng)用可能會低估我國VOCs對臭氧生成的貢獻量, 且對臭氧生成關(guān)鍵物種的判定方面可能與實際會有一定的差異[25]. SOA, VOCs組分的二次有機氣溶膠生成量, g; FAC, VOCs組分的SOA生成系數(shù),以Grosjean[26]的研究數(shù)據(jù)為參考,異戊二烯FAC取0.02,單萜烯FAC取0.3.

基于上述計算結(jié)果,使用Ren等[10]和Ghirardo等[16]提出的簡化的“箱式”模型計算BVOCs排放對空氣質(zhì)量的影響.公式如下:

ozone(i)=OFP×BVOCi∕(×) (5)

SOA(i)=SOA×SOA∕(×) (6)

式中:Ozone(i)、SOA(i),分別代表VOCs組分導(dǎo)致的O3、SOA的濃度增量;為研究區(qū)域的面積, km2;為大氣混合層高度, 1200m[14,27];BVOCi、SOA,分別代表BVOCs和SOA在大氣中的壽命[10].

1.3 模型參數(shù)的確定

1.3.1 標(biāo)準(zhǔn)排放因子 本文中各樹種的異戊二烯和單萜烯的標(biāo)準(zhǔn)排放因子通過文獻調(diào)研獲取.根據(jù)森林資源調(diào)查統(tǒng)計結(jié)果,北京地區(qū)主要的優(yōu)勢樹種分為油松、落葉松、側(cè)柏、楊樹、櫟樹、樺木、刺槐和闊葉樹(包括榆樹、臭椿、銀杏、柳樹、懸鈴木、核桃楸等其他闊葉樹).不同研究之間的樹種排放潛力觀測結(jié)果具有較大差異.為使排放因子數(shù)據(jù)盡可能準(zhǔn)確,本文參考張莉等[28]的方法,優(yōu)先選擇北京及周邊地區(qū)的實測值,取異戊二烯和單萜烯排放速率的平均值作為標(biāo)準(zhǔn)排放因子,結(jié)果見表1[29].OVOCs由于排放速率低、研究較少,其標(biāo)準(zhǔn)排放因子采用Guenther等[30]的推薦值1.5μg C/(g·h).

表1 北京優(yōu)勢樹種BVOCs標(biāo)準(zhǔn)排放因子

1.3.2 葉生物量 根據(jù)北京市第5～9次森林資源調(diào)查數(shù)據(jù),統(tǒng)計優(yōu)勢樹種各齡級的蓄積量,推算葉生物量.具體計算公式如下:

=×T∕T×L

=×T/T×L

式中:是葉生物量, g;是樹種蓄積量, m3;T是樹干密度, 106g/m3;T是樹干生物量占喬木層總生物量的比例, %;L是葉生物量占喬木層總生物量的比例, %.根據(jù)張莉等[28]、呂迪[35]提出的方法,將各樹種的樹干密度、主要森林類型各林齡的T和L值整理為表2.

表2 優(yōu)勢樹種的植物學(xué)參數(shù)

1.3.3 校正因子 (1)溫度校正因子

異戊二烯的溫度校正因子計算公式:

γT= exp[T1(?S)/S] / {1+exp[T2(?M)/

S]} (7)

式中:T1、T2、M均為經(jīng)驗系數(shù),值分別為95000J/ mol、230000J/mol和314K;為氣體平衡常數(shù),值為8.314J/(K·mol);為采樣時葉片實際溫度,單位為K,本文用氣溫代替;S為標(biāo)準(zhǔn)條件下的葉溫,303K.

單萜烯和OVOCs的溫度校正因子計算公式:

γT= exp[(?S)] (8)

式中:經(jīng)驗系數(shù)值為0.09K-1,其它變量與公式(7)相同.

(2)光照校正因子

光照校正因子計算公式如下:

P=(L)∕(1+22)0.5(9)

式中:、L為經(jīng)驗系數(shù),取值分別為0.0027和1.066,為光合量子密度[PPFD, μmol/(m2·s)],來源于太陽輻射.

(3)季節(jié)校正因子

植物的葉齡和葉生物量會隨季節(jié)發(fā)生變化,進而影響B(tài)VOCs排放.參照Staudt等[49]的研究,各月季節(jié)校正因子計算公式如下:

S=1?{1?exp[-(?0)2∕]} (10)

式中:指年排放速率變化幅度與最大排放速率之比,落葉樹種為1,常綠樹種為0.8[10];為月份;0為排放速率最大的月份,北京是7月;為BVOCs進行活躍排放的月數(shù),常綠樹種為12,北京地區(qū)落葉樹種為6[28].

1.3.4 氣象數(shù)據(jù) 本文使用的溫度、光合有效輻射逐小時數(shù)據(jù)來源于北京城市生態(tài)系統(tǒng)研究站、國家氣象科學(xué)數(shù)據(jù)中心及《1961—2014年中國光合有效輻射重構(gòu)數(shù)據(jù)集》[50].

2 結(jié)果與討論

2.1 北京森林BVOCs排放現(xiàn)狀

2020年北京森林BVOCs排放量為39.57×109g C,其中異戊二烯、單萜烯和OVOCs排放量分別為28.57×109, 6.92×109和4.08×109g C,分別占BVOCs總排放量的72.19%、17.48%和10.32%,異戊二烯是北京森林BVOCs的主要組分(圖1a).

根據(jù)樹種的BVOCs排放量對其進行排序(圖1b),結(jié)果顯示不同樹種的排放量及VOCs組成均存在較大差異.楊樹的BVOCs排放量最高,年排放量達(dá)到15.43×109g C,占總排放量的38.99%;櫟樹次之,年排放量12.09×109g C,占總量的30.55%;油松BVOCs年排放量4.05×109g C,對總量的貢獻率為10.25%.闊葉樹種是主要的異戊二烯排放源,其中楊樹和櫟樹對異戊二烯排放的貢獻率分別達(dá)到48.81%和32.01%,這主要是由于二者具有較大的蓄積量和異戊二烯排放速率;而油松等針葉樹種是主要的單萜烯排放源,油松的單萜烯排放量為3.17×109g C,占單萜烯總量的45.81%.

圖1 北京森林優(yōu)勢樹種BVOCs排放量及組成

與前人的研究結(jié)果進行對比,發(fā)現(xiàn)不同地區(qū)BVOCs各組分的占比具有較大差異,整體來看北方地區(qū)植物排放的VOCs中異戊二烯占比較高,由北向南異戊二烯占比呈下降趨勢.桂麗[37]對西安市BVOCs排放的研究顯示,異戊二烯、單萜烯和OVOCs分別占總排放量的44.91%、24.16%和30.93%;劉巖等[51]對長江三角洲地區(qū)的調(diào)查發(fā)現(xiàn),異戊二烯占BVOCs總排放量的37.3%;鄭君瑜等[52]在珠江三角洲地區(qū)的研究則顯示,異戊二烯僅占總排放量的24.7%.與上述研究相比,北京地區(qū)BVOCs中異戊二烯的貢獻率更高,這一差異主要是由于各地區(qū)植被類型的不同導(dǎo)致的.楊屬和櫟屬植物是北京市主要的植被類型,二者均具有較高的異戊二烯排放強度,且葉密度大、葉生物量高,因此北京森林排放的VOCs組分以異戊二烯為主.河北與北京地理位置相近、自然條件相似,樊沖[36]對河北省森林BVOCs排放的研究顯示,異戊二烯、單萜烯和OVOCs對總排放量的貢獻率分別為76.5%、13.6%和9.9%,與北京地區(qū)相似.

北京森林BVOCs的排放量和組分組成均表現(xiàn)出明顯的季節(jié)變化特征(圖2).夏季(6～8月)BVOCs排放量最高,達(dá)到30.74×109g C,占年排放總量的77.67%;春季(3～5月)和秋季(9～11月)排放量分別占總量的11.80%和10.27%;冬季(12～次年2月)排放最低,僅占0.26%.夏季是BVOCs排放的主要時段,此時氣溫和太陽輻射均較高,達(dá)到一年中的最大值,葉面積指數(shù)和酶活性也較高,可以促進植物VOCs的釋放.此外,BVOCs中不同組分的貢獻率也隨季節(jié)發(fā)生變化,其中異戊二烯的變化幅度最大,1～7月份占比呈上升趨勢,7～12月份呈下降趨勢;單萜烯貢獻率的變化趨勢則與之相反,呈冬季高、夏季低.北京地區(qū)的異戊二烯排放源主要是落葉闊葉樹,隨季節(jié)變化呈現(xiàn)夏盛冬衰的規(guī)律;單萜烯的排放源主要是針葉樹種,葉片四季常綠,在全年均有排放.因此,與其他BVOCs組分相比,異戊二烯排放量對季節(jié)變化更為敏感.

圖2 不同月份BVOCs的排放量及組成

2.2 2000～2020年BVOCs排放變化趨勢

如表3所示, 2000～2020年北京市森林源BVOCs年排放量從20.30×109g C上升到39.57×109g C,增加了近1倍,年平均增長率4.75%.北京地區(qū)BVOCs年排放量的增長率高于全國平均水平,Li等[1,53]的研究顯示,1981～2003年和2008～2018年間,全國BVOCs排放年增長率分別為1.27%和2.03%.北京森林BVOCs排放量增加,主要是由于樹木生長以及植樹造林引起的生物量增加導(dǎo)致的.從圖3可以看出,二者隨時間的變化趨勢基本保持一致,呈逐年上升的態(tài)勢.北京森林BVOCs排放量的變化表現(xiàn)出明顯的階段性特征,2000～2010年間增長較為緩慢,2010～2020年間則出現(xiàn)大幅度上升,這與北京市百萬畝造林等綠化工程的實施密切相關(guān).

2000～2020年,異戊二烯所占比重呈下降趨勢,單萜烯和OVOCs所占比重則呈上升趨勢. 20年間異戊二烯、單萜烯和OVOCs的年排放量分別增加了12.44×109, 4.29×109和2.54×109g C(表3),分別占BVOCs總增量的64.54%、22.28%和13.16%.BVOCs各組分在不同調(diào)查階段所占比重不同,異戊二烯所占比重逐漸下降,2000年占總排放量的79.45%,2020年下降到72.19%;單萜烯和OVOCs所占比重則逐漸上升,單萜烯占比從12.93%增加到17.48%, OVOCs占比從7.63%增加到10.32%.楊樹作為異戊二烯的強排放源,由于生長速度快、可快速成林,在北京早期植樹造林中得到廣泛應(yīng)用,導(dǎo)致異戊二烯排放量占比較高;但楊樹容易招病蟲害、易倒伏、生長后勁不足,且產(chǎn)生的楊絮容易引發(fā)過敏性疾病等,近年來種植逐漸減少,而油松、側(cè)柏等樹種的種植逐漸增多,因此異戊二烯排放量占比下降、單萜烯占比增加.

表3 2000～2020年北京森林BVOCs排放量

近幾十年來,研究人員對北京市BVOCs排放特征進行了較多研究,由于使用的數(shù)據(jù)源和算法不同,導(dǎo)致不同研究結(jié)果之間存在差異(圖3).Wang等[54]對1998年北京市BVOCs排放量的估算結(jié)果為16.2×109g C;Ren等[10]估算了1990年和2015年北京市BVOCs的排放量,結(jié)果分別為10.9×109和33.0×109g C;張倩等[9]和井瀟溪[15]分別對2015年和2017年北京BVOCs排放量進行測算,結(jié)果分別為35.5×109和43.6×109g C.

圖3 不同研究對北京市BVOCs排放量的估算

從圖3可以看出,本文測算的北京市BVOCs排放量和年際變化趨勢與上述研究成果較為接近,但是低于Klinger等[6]和謝揚飏等[11]的計算結(jié)果.Klinger等[6]在2001年對全國生態(tài)系統(tǒng)的VOCs排放量進行了測算,其中北京市BVOCs釋放量為48.0×109g C;謝揚飏等[11]依據(jù)園林綠地普查結(jié)果,估算得到2000年北京綠地天然源VOCs排放量為38.5×109g C.不同研究所使用的土地覆蓋數(shù)據(jù)的不同是導(dǎo)致估算結(jié)果產(chǎn)生差異的主要原因之一,Klinger等[6]所使用的森林面積是Wang等[54]的約3倍,因此其估算結(jié)果較高.此外,各研究使用的估算模型不同也會導(dǎo)致結(jié)果產(chǎn)生較大差異.如謝揚飏等[11]使用Guenther等[21,55]提出的算法進行BVOCs測算,但未考慮葉齡隨季節(jié)變化的影響,導(dǎo)致估算結(jié)果偏大.

2.3 優(yōu)勢樹種的BVOCs排放變化趨勢

對2000～2020年間各優(yōu)勢樹種的BVOCs排放變化趨勢進行分析,首先從排放量來看,楊樹、櫟樹、闊葉樹、油松和側(cè)柏的排放量均呈逐年上升的趨勢;落葉松、樺木和刺槐的排放量整體上呈增加趨勢,在個別年份出現(xiàn)下降(表4).這與各樹種葉生物量的變化趨勢相一致(表5),表明其BVOCs排放量的變化主要是由葉生物量變化引起的.其中,櫟樹、闊葉樹、楊樹和油松的BVOCs年排放量分別增加7.17×109, 4.09×109, 3.93×109和2.69×109g C,占全市森林BVOCs排放增量的92.84%.特別是櫟樹,是BVOCs排放增量最多的森林類型,占總增加量的37.23%;同時其葉生物量增量也最高,占總增量的27.15%.在北京地區(qū),櫟樹主要分布在山區(qū),以蒙古櫟為主的櫟屬植物是北京山區(qū)面積最大、分布最廣的森林類型,也是北京山區(qū)地帶性的落葉闊葉林頂級群落之一.

表4 2000～2020年優(yōu)勢樹種BVOCs排放量及貢獻率

從貢獻率來看,在所有林分類型中,楊樹、櫟樹和油松對BVOCs排放量的貢獻在各調(diào)查階段均較高,占總排放量的79.79%～88.85%.其中,楊樹的貢獻率最高,2000年占森林BVOCs總排放量的56.61%,之后整體呈現(xiàn)下降趨勢,到2020年貢獻率減少到38.99%;與之相對應(yīng)的是櫟樹,2000年對森林BVOCs排放的貢獻率為24.22%,隨后保持穩(wěn)中有升的態(tài)勢,到2020年占比達(dá)到30.55%.此外,闊葉樹對森林BVOCs的貢獻率也表現(xiàn)出明顯的變化趨勢, 2000年占比為6.50%,到2020年占比增加到13.68%.特別是2015～2020年間,闊葉樹占比上升了5.97%,而同一時期楊樹占比下降了13.49%.據(jù)統(tǒng)計,林業(yè)資源調(diào)查中的闊葉樹包括榆樹、銀杏、柳樹、懸鈴木、核桃楸等30余種樹種,表明北京新增的森林植被更加注重物種多樣化.

表5 2000～2020年優(yōu)勢樹種葉生物量(109g)

2.4 BVOCs排放對空氣質(zhì)量的影響

如表6所示, O3生成潛勢從2000年的181.76×109g增加到2020年的331.07×109g.其中,異戊二烯對O3生成潛勢的貢獻率從94.17%下降到91.56%,平均貢獻率為92.70%,是主要的貢獻者;單萜烯對O3生成潛勢的貢獻率從5.83%增加到8.44%,平均貢獻率為7.30%.SOA 生成潛勢從2000年的1.11×109g增加到2020年的2.65×109g,其中異戊二烯對SOA生成潛勢的貢獻率從29.07%下降到21.59%,平均貢獻率為24.60%;單萜烯的貢獻率從70.93%增加到78.41%,平均貢獻率為75.40%,是主要的貢獻者.O3生成潛勢最大的樹種是楊樹,SOA生成潛勢最大的樹種是油松.

表6 2000～2020年北京森林BVOCs的O3和SOA生成潛勢

為量化北京森林源BVOCs對區(qū)域空氣質(zhì)量的影響,采用簡化的“箱式”模型進行估算.結(jié)果顯示,近20年來森林BVOCs排放對大氣環(huán)境中O3和SOA濃度的貢獻量逐年上升(圖4).截止到2020年,BVOCs排放產(chǎn)生的O3生成潛勢導(dǎo)致全市范圍內(nèi)大氣年均O3濃度上升20.50μg/m3,全市大氣年均SOA濃度增加2.36μg/m3.

圖4 BVOCs排放對空氣質(zhì)量的影響

圖5 2020年BVOCs排放對環(huán)境O3污染的貢獻

環(huán)境O3數(shù)據(jù)來源于北京市生態(tài)環(huán)境局《2016-2020北京市生態(tài)環(huán)境質(zhì)量報告書》

森林源BVOCs排放具有明顯的季節(jié)變化特征,由其引起的大氣O3濃度增加主要集中在夏季.從圖5可以看出,在排放量最高的6～8月份,森林BVOCs導(dǎo)致大氣O3濃度平均上升65.23μg/m3.在其他相關(guān)研究中,也發(fā)現(xiàn)了類似的規(guī)律.Ren等[10]的研究顯示,2015年BVOCs排放導(dǎo)致北京城市核心區(qū)7月份臭氧濃度增加71.5μg/m3,占環(huán)境濃度的30%.Li等[56]在長江三角洲地區(qū)的研究發(fā)現(xiàn),BVOCs排放對日間O3濃度的貢獻在18.6%～24.9%之間.此外,雖然北京地區(qū)植物源VOCs排放量低于人為源,但是植物源VOCs組分具有更高的反應(yīng)活性和O3形成潛力[14-15],因此其對大氣O3污染的影響不容忽視.

為改善環(huán)境空氣質(zhì)量,未來進行城市綠化時,應(yīng)優(yōu)先考慮VOCs排放較低的樹種,建議在現(xiàn)有基礎(chǔ)上進一步減少楊樹等O3生成潛勢較大的樹種;此外,雖然BVOCs對PM2.5影響較小,但目前貢獻占比最大的油松為北京城市綠化主種的樹種之一,應(yīng)考慮總量控制.

3 結(jié)論

3.1 2020年,北京市森林BVOCs排放量為39.57× 109g C,其中異戊二烯、單萜烯和OVOCs分別占總排放量的72.19%、17.48%和10.32%.北京主要的優(yōu)勢樹種中,對BVOCs排放貢獻較高的是楊樹、櫟樹和油松,闊葉樹種是主要的異戊二烯排放源,油松等針葉樹種是主要的單萜烯排放源.

3.2 2000～2020年,北京市森林BVOCs年排放量從20.30×109g C上升到39.57×109g C,年平均增長率4.75%.BVOCs排放量的變化表現(xiàn)出明顯階段性特征,2000～2010年增長緩慢,2010～2020年則出現(xiàn)大幅上升,這與北京市百萬畝造林等綠化工程的實施密切相關(guān).研究期內(nèi),異戊二烯所占比重呈下降趨勢,單萜烯和OVOCs所占比重呈上升趨勢;異戊二烯的強排放源楊樹的占比減少,櫟樹和其他闊葉樹的占比增加,北京新增的森林植被更加注重物種多樣化.

3.3 2000～2020年,北京森林BVOCs的O3生成潛勢從181.76×109g增加到331.07×109g,其中異戊二烯是主要的貢獻者,占92.70%;SOA 生成潛勢從1.11× 109g增加到2.65×109g,單萜烯是主要貢獻者,占75.40%.O3生成潛勢最大的樹種是楊樹,SOA生成潛勢最大的樹種是油松.

3.4 2020年,BVOCs排放導(dǎo)致全市大氣年均O3濃度上升20.50μg/m3, SOA濃度增加2.36μg/m3,森林BVOCs排放在夏季對O3污染的影響最為明顯.未來城市綠化中應(yīng)考慮優(yōu)化樹種組成,以減少BVOCs對環(huán)境空氣質(zhì)量的影響.

[1] Li L Y, Yang W Z, Xie S D, et al. Estimations and uncertainty of biogenic volatile organic compound emission inventory in China for 2008～2018 [J]. Science of the Total Environment, 2020,733:139301.

[2] Guenther A, Hewitt C, Erickson D, et al. A global-model of natural volatile organic-compound emission [J]. Journal of Geophysical Research Atmospheres, 1995,100(5):8873-8892.

[3] Fehsenfeld F, Calvert J, Fall R, et al. Emissions of volatile organic compounds from vegetation and the implications for atmospheric chemistry [J]. Global Biogeochemical Cycles, 1992,6:389-430.

[4] Guenther A, Karl T, Harley P, et al. Estimates of global terrestrial isoprene emissions using MEGAN (Model of Emissions of Gases and Aerosols from Nature) [J]. Atmospheric Chemistry and Physics, 2006,6:3181-3210.

[5] Graedel T E, Bates T S, Bouwman A F, et al. A compilation of inventories of emissions to the atmosphere [J]. Global Biogeochemical Cycles, 1993,7:1–26.

[6] Klinger L F, Li Q J, Guenther A B, et al. Assessment of volatile organic compound emissions from ecosystems of China [J]. Journal of Geophysical Research, 2002,107(D21):1-21.

[7] Li L Y, Chen Y, Xie S D. Spatio-temporal variation of biogenic volatile organic compounds emissions in China [J]. Environmental Pollution, 2013,182:157-168.

[8] Yin L, Xu Z, Liu M, et al. Estimation of biogenic volatile organic compound (BVOC) emissions in China using WRF–CLM–MEGAN coupled model, Biogeosciences Discuss. [preprint], https://doi.org/ 10.5194/bg-2019-458, 2020.

[9] 張 倩,李洪遠(yuǎn).城市綠地對大氣環(huán)境的潛在負(fù)面效應(yīng)[J]. 景觀設(shè)計, 2019,18(2):4-11.

Zhang Q, Li H Y. Potential negative effects of urban green space on the atmospheric environment [J]. Landscape Design, 2019,18(2):4-11.

[10] Ren Y, Qu Z, Du Y, et al. Air quality and health effects of biogenic volatile organic compounds emissions from urban green spaces and the mitigation strategies [J]. Environmental Pollution, 2017,230:849- 861.

[11] 謝揚飏,邵 敏,陸思華,等.北京市園林綠地植被揮發(fā)性有機物排放的估算[J]. 中國環(huán)境科學(xué), 2007,27(4):498-502.

Xie Y Y, Shao M, Lu S H, et al. The estimation of volatile organic compounds emission from landscape plants in Beijing [J]. China Environmental Science, 2007,27(4):498-502.

[12] Xie X, Shao M, Liu Y, et al. Estimate of initial isoprene contribution to ozone formation potential in Beijing, China [J]. Atmospheric Environment, 2008,42(24):6000-6010.

[13] Pang X, Mu Y, Zhang Y, et al. Contribution of isoprene to formaldehyde and ozone formation based on its oxidation products measurement in Beijing, China [J]. Atmospheric Environment, 2009, 43(13):2142-2147.

[14] Mo Z, Shao M, Wang W, et al. Evaluation of biogenic isoprene emissions and their contribution to ozone formation by ground-based measurements in Beijing, China [J]. Science of the Total Environment, 2018,627:1485-1494.

[15] 井瀟溪.北京市森林植物揮發(fā)性有機物排放研究[D]. 北京:北京林業(yè)大學(xué), 2020.

Jing X X. Study on biogenic volatile organic compounds emission from forest plants in Beijing [D]. Beijing: Beijing Forestry University, 2020.

[16] Ghirardo A, Xie J, Zheng X, et al. Urban stress-induced biogenic VOC emissions and SOA-forming potentials in Beijing [J]. Atmospheric Chemistry and Physics, 2016,16(5):2901-2920.

[17] 方精云,陳安平.中國森林植被碳庫的動態(tài)變化及其意義[J]. 植物學(xué)報, 2001,43(9):967-973.

Fang J Y, Chen A P. Dynamic forest biomass carbon pools in China and their significance [J]. Bulletin of Botany, 2001,43(9):967-973.

[18] 張 莉.中國森林生態(tài)系統(tǒng)異戊二烯排放研究[D]. 南京:南京氣象學(xué)院, 2002.

Zhang L. Isoprene emission from forest ecosystem in China [D]. Nanjing: Nanjing Institute of Meteorology, 2002.

[19] Li L Y, Xie S D. Historical variations of biogenic volatile organic compound emission inventories in China, 1981～2003 [J]. Atmospheric Environment, 2014,95:185-196.

[20] 北京市園林綠化局.2020年森林資源情況[EB/OL]. http://yllhj. beijing.gov.cn/zwgk/tjxx/.

Beijing Municipal Forestry and Parks Bureau. Forest resources in 2020 [EB/OL]. http://yllhj.beijing.gov.cn/zwgk/tjxx/.

[21] Guenther A, Baugh B, Brasseur G, et al. Isoprene emission estimates and uncertainties for the central African EXPRESSO study domain [J]. Journal of Geophysical Research: Atmospheres, 1999,104(D23): 30625-30639.

[22] 馮惠生.北京生態(tài)環(huán)境保護[M]. 北京:中國環(huán)境出版集團, 2019.

Feng H S. Beijing ecological and environmental protection [M]. Beijing: China Environment Publishing Group, 2019.

[23] Benjamin M T, Winer A M. Estimating the ozone-forming potential of urban trees and shrubs [J]. Atmospheric Environment, 1998,32:53-68.

[24] Carter W P L. Reactivity estimates for selected consumer product compounds [R]. California: Center for Environmental Research and Technology, College of Engineering, University of California, February 19,2008.

[25] 邱婉怡,劉禹含,譚照峰,等.基于中國四大城市群計算的最大增量反應(yīng)活性[J]. 科學(xué)通報, 2020,65(7):610-621.

Qiu W Y, Liu Y H, Tan Z F, et al. Calculation of maximum incremental reactivity scales based on typical megacities in China [J]. Science Bulletin, 2020,65(7):610-621.

[26] Grosjean D. In situ organic aerosol formation during a smog episode: Estimated production and chemical functionality [J]. Atmospheric Environment, Part A, General Topics, 1992,26(6):953-963.

[27] 杜吳鵬,房小怡,黃宏濤,等.北京近年地表風(fēng)速和大氣混合層厚度變化特征研究[J]. 環(huán)境科學(xué)與技術(shù), 2017,40(6):149-156.

Du W P, Fang X Y, Huang H T, et al. Variation characteristics of surface wind speed and atmospheric mixing layer height in recent years in Beijing [J]. Environmental Science and Technology, 2017, 40(6):149-156.

[28] 張 莉,王效科,歐陽志云,等.中國森林生態(tài)系統(tǒng)的異戊二烯排放研究[J]. 環(huán)境科學(xué), 2003,24(1):8-15.

Zhang L, Wang X K, Ooyang Z Y, et al. Estimation of isoprene emission from forest ecosystems in China [J]. Environmental Science, 2003,24(1):8-15.

[29] 張 薔,李令軍,趙文慧,等.北京森林BVOCs排放特征及對區(qū)域空氣質(zhì)量的影響[J]. 中國環(huán)境科學(xué), 2021,41(2):622-632.

Zhang Q, Li L J, Zhao W H, et al. Emission characteristics of VOCs from forests and its impact on regional air quality in Beijing [J]. China Environmental Science, 2021,41(2):622-632.

[30] Guenther A, Zimmerman P, Wildermuth M. Natural volatile organic compound emission rate estimates for U.S. woodland and landscapes [J]. Atmospheric Environment, 1994,28(6):1197-1210.

[31] 王志輝,張樹宇,陸思華,等.北京地區(qū)植物VOCs排放速率的測定[J]. 環(huán)境科學(xué), 2003,24(2):7-12.

Wang Z H, Zhang S Y, Lu S H, et al. Screenings of 23plant species in Beijing for volatile organic compound emissions [J]. Environmental Science, 2003,24(2):7-12.

[32] 趙 靜,白郁華,王志輝,等.我國植物VOCs排放速率的研究[J]. 中國環(huán)境科學(xué), 2004,24(6):654-657.

Zhao J, Bai Y H, Wang Z H, et al. Studies on the emission rates of plants VOCs in China [J]. China Environmental Science, 2004,24(6): 654-657.

[33] 趙美萍,邵 敏,白郁華,等.我國幾種典型樹種非甲烷烴類的排放特征[J]. 環(huán)境化學(xué), 1996,15(1):69-75.

Zhao M P, Shao M, Bai Y H, et al. Study on NMHC emission characteristics of several typical trees in China [J]. Environmental Chemistry, 1996,15(1):69-75.

[34] 陳 穎,李德文,史 奕,等.沈陽地區(qū)典型綠化樹種生物源揮發(fā)性有機物的排放速率[J]. 東北林業(yè)大學(xué)學(xué)報, 2009,37(3):49-51.

Chen Y, Li D W, Shi Y, et al. Emission rate of biogenic volatile organic compounds from urban trees in Shenyang, China [J]. Journal of Northeast Forestry University, 2009,37(3):49-51.

[35] 呂 迪.陜西森林植被BVOCs排放特征研究[D]. 楊凌:西北農(nóng)林科技大學(xué), 2016.

Yang D. Estimation of BVOCs emission from forests in Shaanxi [D]. Yangling: Northwest A & F University, 2016.

[36] 樊 沖.河北省森林植物源揮發(fā)性有機物排放量研究[D]. 北京:北京林業(yè)大學(xué), 2019.

Fan C. Study on biogenic volatile organic compounds emission from forest plants in Hebei province [D]. Beijing: Beijing Forestry University, 2019.

[37] 桂 麗.西安市揮發(fā)性有機物排放量估算及其臭氧生成潛勢研究[D]. 西安:西北大學(xué), 2019.

Gui L. Estimation of volatile organic organic compounds emissions and ozone formation potential in Xi’an city [D]. Xi’an: Northwest University, 2019.

[38] Tani A, Kawawata Y. Isoprene emission from the major native Quercus spp. in Japan [J]. Atmospheric Environment, 2008,42(19): 4540-4550.

[39] Yaman B, Aydin Y M, Koca H, et al. Biogenic volatile organic compound (BVOC) emissions from various endemic tree species in Turkey [J]. Aerosol and Air Quality Research, 2015,15(1):341-356.

[40] 張福珠,苗 鴻,魯 純.落葉闊葉林釋放異戊二烯的研究[J]. 環(huán)境科學(xué), 1994,19(1):1-5.

Zhang F Z, Miao H, Lu C. Emission of isoprene from deciduous forest [J]. Environmental Science, 1994,19(1):1-5.

[41] Tiwary A, Namdeo A, Fuentes J, et al. Systems scale assessment of the sustainability implications of emerging green initiatives [J]. Environmental Pollution, 2013,183:213-223.

[42] 王效科,牟玉靜,歐陽志云,等.太湖流域主要植物異戊二烯排放研究[J]. 植物學(xué)通報, 2002,20(2):224-230.

Wang X K, Mu Y J, Ouyang Z Y, et al. Study on emission of isoprene from major plants living in Taihu basin [J]. Chinese Bulletin of Botany, 2002,20(2):224-230.

[43] Aydin Y M, Yaman B, Koca H, et al. Biogenic volatile organic compound (BVOC) emissions from forested areas in Turkey: Determination of specific emission rates for thirty-one tree species [J]. Science of the Total Environment, 2014,490:239-253.

[44] Prendez M, Carvajal V, Corada K, et al. Biogenic volatile organic compounds from the urban forest of the Metropolitan Region, Chile [J]. Environmental Pollution, 2013,183:143-150.

[45] Noe S M, Penuelas J, Niinemets U. Monoterpene emissions from ornamental trees in urban areas: a case study of Barcelona, Spain [J]. Plant Biology, 2008,10(1):163-169.

[46] 牟玉靜,宋文質(zhì),張曉山,等.落葉闊葉樹異戊二烯排放研究[J]. 環(huán)境化學(xué), 1999,18(1):21-27.

Mu Y J, Song W Z, Zhang X S, et al. Study on emissions of isoprene from deciduous and broadleaf trees [J]. Environmental Chemistry, 1999,18(1):21-27.

[47] 白郁華,李金龍,夏 羽,等.楊樹排放碳?xì)浠衔锏南嚓P(guān)因素[J]. 環(huán)境化學(xué), 1995,14(2):118-123.

Bai Y H, Li J L, Xia Y, et al. Relative factor of hydrocarbons emission from popular tree [J]. Environmental Chemistry, 1995,14(2):118-123.

[48] 張鋼鋒,謝紹東.基于樹種蓄積量的中國森林VOC排放估算 [J]. 環(huán)境科學(xué), 2009,30(10):2816-2822.

Zhang G F, Xie S D. Estimation of VOC emission from forests in China based on the volume of tree species [J]. Environmental Science, 2009,30(10):2816-2822.

[49] Staudt M, Bertin N, Frenzel B, et al. Seasonal variation in amount and composition of monoterpenes emitted by young Pinus pinea trees -Implications for emission modeling [J]. Journal of Atmospheric Chemistry, 2000,35(1):77-99.

[50] 唐利琴,劉 慧,胡 波,等.1961～2014年中國光合有效輻射重構(gòu)數(shù)據(jù)集[J/OL]. 中國科學(xué)數(shù)據(jù), 2017,2(3).DOI: 10.11922/csdata.170. 2017.0135.

Tang L Q, Liu H, Hu B, et al. A dataset of reconstructed photosynthetically active radiation in China (1961 – 2014) [J/OL]. China Scientific Data, 2017,2(3).DOI:10.11922/csdata.170.2017.0135.

[51] 劉 巖,李 莉,安靜宇,等.長江三角洲2014年天然源BVOCs排放、組成及時空分布[J]. 環(huán)境科學(xué), 2018,39(2):608-617.

Liu Y, Li L, An J Y, et al. Emissions, chemical composition, and spatial and temporal allocation of the BVOCs in the Yangtze River Delta Region in 2014 [J]. Environmental Science, 2018,39(2):608-617.

[52] 鄭君瑜,鄭卓云,王兆禮,等.珠江三角洲天然源VOCs排放量估算及時空分布特征[J]. 中國環(huán)境科學(xué), 2009,29(4):345-350.

Zheng J Y, Zheng Z Y, Wang Z L, et al. Biogenic VOCs emission inventory and its temporal and spatial characteristics in the Pearl River Delta area [J]. China Environmental Science, 2009,29(4):345-350.

[53] Li L Y, Xie S D. Historical variations of biogenic volatile organic compound emission inventories in China, 1981～2003 [J]. Atmospheric Environment, 2014,95:185-196.

[54] Wang Z H, Bai Y H, Zhang S Y. A biogenic volatile organic compounds emission inventory for Beijing [J]. Atmospheric Environment, 2003,37(27):3771-3782.

[55] Guenther A, Geron C, Pierce T, et al. Natural emissions of non- methane volatile organic compounds, carbon monoxide, and oxides of nitrogen from North America [J]. Atmospheric Environment, 2000, 34(12-14):2205-2230.

[56] Li L, An J Y, Shi Y Y, et al. Source apportionment of surface ozone in the Yangtze River Delta, China in the summer of 2013 [J]. Atmospheric Environment, 2016,144:194-207.

The current status and characteristics of variations in forest BVOCs emissions in Beijing.

ZHANG Qiang1,2, LI Ling-jun1,2*, ZHAO Wen-hui1,2, LIU Bao-xian1,2,3, LU Hai-feng1,2, QIN Li-huan1,2, CHEN Yun-ru1,2, LI Qi1,2

(1.Beijing Municipal Ecological and Environment Monitoring Center, Beijing 100048, China；2.Beijing Key Laboratory of Airborne Particulate Matter Monitoring Technology, Beijing 100048, China；3.School of Environment, Tsinghua University, Beijing 100084, China)., 2022,42(9)：3991～4000

Accurate estimation of biogenic volatile organic compounds (BVOCs) emissions is necessary for developing effective air pollution control measures. To understand the current status and characteristics of variations in forest BVOCs emissions in Beijing, historical BVOCs emission inventories for the period from 2000 to 2020 were developed based on the volume data of forest resources survey, hourly meteorological data and light-temperature model. The results show that the forest BVOCs emissions in 2020 were 39.57×109g C, of which isoprene, monoterpene and other VOCs accounted for 72.19%, 17.48% and 10.32%, respectively. Broadleaf trees emitted the largest amount of isoprene and total BVOC emissions, while conifers made the largest monoterpene emission. Results show that Beijing’s annual forest BVOCs emissions had increased from 20.30×109g C in 2000 to 39.57×109g Cin 2020, with an annual average increase rate of 4.75%. Isoprene, monoterpene and other VOCs accounted for 64.54%, 22.28% and 13.16% of the total increment, respectively. In the past 20 years, the proportion of isoprene decreased, while the proportion of monoterpene and OVOCs increased. The contribution to the BVOCs emissions tended to decrease for poplar while to increase significantly for quercus and other broadleaf trees, indicating that the new afforestation in Beijing has led to a higher species diversity. From 2000 to 2020, the ozone formation potential of BVOCs had increased from 181.76×109to 331.07×109g C, of which isoprene was the main contributor, accounting for 92.70%. The secondary organic aerosol formation potential of BVOCs had increased from 1.11×109to 2.65×109g C, and the contribution of monoterpene and isoprene was 75.40% and 24.60%, respectively. Poplar exhibits the highest ozone formation potential, andshows the highest secondary organic aerosol formation potential. The forest BVOCs emissions made relatively higher contribution to ozone pollution in summer. Therefore, tree species optimization should be considered to mitigate BVOCs emissions.

timber volume；isoprene；monoterpene；ozone；secondary organic aerosol；volatile organic compounds

X511

A

1000-6923(2022)09-3991-10

2022-02-21

國家重點研發(fā)計劃資助項目(2021YFC1809000);北京市科技計劃(Z211100004321006)

*責(zé)任作者, 教授級高級工程師, lilj2000@126.com

張 薔(1992-),女,山東聊城人,高級工程師,博士,主要從事生態(tài)遙感、植物化學(xué)計量等方面的研究.發(fā)表論文8篇.