北京森林BVOCs 排放特征及對區域空氣質量的影響

張 薔,李令軍,趙文慧,王新輝,姜 磊,劉保獻,李 翔,鹿海峰 (1.北京市生態環境監測中心,北京 100048；2.大氣顆粒物監測技術北京市重點實驗室,北京 100048；3.北京市生態環境局,北京 100044；4.清華大學環境學院,北京 100084)

近年來,北京市大氣污染治理取得顯著成效,顆粒物特別是PM2.5濃度逐年降低,但夏季臭氧(O3)污染居高不下,逐漸成為制約空氣質量改善的關鍵因素[1].O3的產生是揮發性有機物(VOCs)和多種氮氧化物(NOx)發生光化學反應的結果,同時VOCs 還是二次有機氣溶膠(SOA)的重要前體物,對大氣中顆粒物的形成具有重要貢獻[2-4].VOCs 按照來源可劃分為天然源VOCs(BVOCs)和人為源VOCs(AVOCs),研究發現全球每年BVOCs 的排放量約為1150Tg,占VOCs 年排放總量的90%以上,遠大于AVOCs 的排放量[2].在北京,由于近年來實施了嚴格的減排措施,如限制車輛使用、實施清潔生產等,AVOCs 排放得到有效控制[5-6];另一方面,調查顯示北京市森林覆蓋率逐年上升,從2006 年的35.9%增加到2018 年的43.5%[7],BVOCs 的排放量也隨之不斷增加[8-9].因此,研究森林植被BVOCs 的排放特征對制定更加有效的O3及顆粒物污染控制措施至關重要.

迄今為止,國內外學者針對VOCs 的天然源清單已經進行了大量研究,從采樣方法、分析方法到區域BVOCs 排放通量的計算都取得很大進展,并建立了不同空間尺度下的BVOCs 排放模型,如BEIS、BEIS2、GLOBEIS、G95 和MEGAN 模型等[2,10-16].對于北京地區BVOCs 的研究也取得了一些成果,如王志輝等[17]和趙靜等[18]測定了北京市多種典型植物的VOCs 排放速率,謝揚飏等[19]和Ghirardo 等[9]利用園林普查資料估算了北京市園林綠地植被的VOCs 排放量,Wang 等[20]和Ren 等[21]則分別對全北京植被的VOCs 排放量進行了估算.但是目前對北京市BVOCs 排放情況的研究多是通過植被的分布面積資料推算出葉生物量,進而估算植被的排放量,且不同研究之間的結果差異較大(相差3～10 倍)[22].已有的研究很少基于樹種蓄積量、以具體樹種為對象對北京地區森林源BVOCs 的排放特征進行分析,且在BVOCs 對區域空氣質量的具體影響方面分析不足.因此,本文以北京市森林優勢樹種蓄積量為基礎,通過Guenther 等提出的光溫影響模型[2,23],模擬各樹種的 BVOCs 排放特征,建立北京森林源BVOCs 排放清單,并探討BVOCs 排放的時間變化規律和在不同齡級林分中的分配規律;此外,進一步估算BVOCs 的O3和SOA 生成潛勢,并嘗試將其對環境空氣質量的影響進行量化,為未來評估植被變化可能對空氣質量產生的影響以及制定有效的管理措施提供數據支撐.

1 材料與方法

1.1 研究區概況

北京位于華北平原西北部邊緣,地域范圍39°26′～41°04′N,115°25′～117°30′E,面積1.64×104km2.北京地勢西北高、東南低,整體由西北東三面環山區、東南平原區和中部城市核心區3 部分組成.北京地區屬溫帶濕潤大陸性季風氣候,夏季高溫多雨,冬季寒冷干燥,多年平均氣溫12.3℃ ,多年平均降水量571.9mm 左右.地帶性植被為暖溫帶落葉闊葉林和溫帶針葉林,主要是以櫟屬(Quercus)和楊屬(Populus)等落葉闊葉樹種占優勢的落葉闊葉林,及側柏(Platycladus orientalis)和油松(Pinus tabulaeformis)占優勢的溫性針葉林.截止到2015 年,北京市森林面積 7449.6km2,森林覆蓋率達 41.6%,森林蓄積量1701.1×104m3.

1.2 研究方法

1.2.1 BVOCs 排放量計算模型 BVOCs 種類繁多,其中異戊二烯和單萜烯由于排放量大、活性強而受到廣泛關注[19,24].本研究將BVOCs 劃分為異戊二烯、單萜烯和其他VOCs(OVOCs) 3 類,OVOCs 物種主要包括醇、醛、酮、有機酸、低碳烷烴和烯烴等.參照Guenther 等[2,23]提出的光溫影響模型,按樹種分別進行排放量估算.其中,異戊二烯的計算采用葉綠素(CHL)算法,其排放主要受溫度和太陽輻射的影響;單萜烯和OVOCs 的計算采用防御專屬組織(DST)算法,其排放主要受溫度影響.

CHL 算法:

式中:EISOP為異戊二烯排放量;ε 為標準狀態下[溫度為303K,光合有效輻射(PAR)為1000μmol photon/(m2·s)]該樹種的基本排放速率, μg C/(g·h); B 為葉生物量; γP、γT和γS為環境校正因子,分別代表光合有效輻射、溫度和季節的影響; ρ 為逸出效率,表示排放總量中逸入樹冠上方大氣中的部分所占比例.

γT的計算公式為:

式中:T 為當前葉表面溫度,K,簡單用氣溫代替; TS為標準條件下的葉溫,303K; R 為氣體平衡常數,8.314J/(K·mol); CT1、CT2、TM為經驗常數,CT1=95000J/mol; CT2= 230000J/mol,TM= 314K.

γP的計算公式為:

式中: a、CL為經驗參數,取值分別為0.0027 和1.066,Q 為當前光合量子密度[PPFD, μmol/(m2·s)],來源于太陽輻射.

葉齡和葉生物量會隨季節發生變化,從而影響BVOCs 的排放.根據Staudt 等[25]的研究,各月季節校正因子γS的計算公式為:

式中: M 指月份; M0指具有最大排放速率的月份; n指一年中BVOCs 進行活躍排放的月數; λ 為年排放速率變化幅度與最大排放速率之比[26].DST 算法:

式中: ETMT和EOVOC分別為單萜烯和OVOC 排放量,其它變量與前述異戊二烯排放公式相同.由于單萜烯和其他VOCs 排放量僅與溫度相關,故僅考慮溫度校正因子γT,計算公式為γT= exp [ β ( T - TS)] ,經驗系數β 值為0.09K-1.

1.2.2 O3和SOA 生成潛勢估算模型 本研究通過估算VOCs 的O3和SOA 生成潛勢來表征其對大氣環境的潛在影響.BVOCs 種類繁多,不同組分間的反應活性差異極大,本文僅對在BVOCs 中占比較大的異戊二烯和單萜烯進行估算. O3生成潛勢(OFP)代表VOCs在最佳反應條件下對O3生成的最大貢獻,采用“增量反應”方法進行計算,即每加入一單位的VOCs 所形成的O3量[27-28].SOA 生成潛勢采用“氣溶膠生成系數(FAC)”法進行計算[29].計算公式為:

式中: OFPi為VOC 組分i 生成O3的最大值; Ei為VOCs 組分i 的排放量或濃度; MIRi為組分i 的最大增量反應活性,MIR 的取值采用Carter[30]的研究成果,其中異戊二烯為 10.61g/g,單萜烯為 4.04g/g;SOAi為VOC 組分i 的二次有機氣溶膠生成量; FACi為VOC 組分i 的SOA 生成系數.參考Grosjean[31]的煙霧箱試驗,假設BVOCs僅與OH-反應生成SOA,異戊二烯FAC 取0.02,單萜烯FAC 取0.3.

在上述計算結果的基礎上,利用的一個簡化的“箱式”模型[9,21]來估算BVOCs 排放對大氣質量的影響:

式中: COzone(i)、CSOA(i)分別為VOC 組分i 導致的O3、SOA 的濃度增量; S 為研究區面積, km2; h 為逆溫層的一般高度, 2km2[9];τBVOCi、τSOA分別為BVOCs 和SOA 在大氣中的壽命[21].

1.3 參數的確定

1.3.1 標準排放因子 對各物種的排放潛力進行準確測定,是正確估算排放量的重要基礎.本文中植物異戊二烯、單萜烯的標準排放因子數據通過文獻獲取.目前文獻報道的排放因子觀測結果多為一定域值,且不同研究之間的觀測值差異較大[26].為使排放因子數據盡可能準確,本文廣泛收集北京市及周邊地區的大量觀測結果,并參考張莉等[26]的方法,取異戊二烯和單萜烯排放速率的最小值、最大值和平均值,分別計算各樹種異戊二烯和單萜烯的最小、最大可能排放量和平均排放量.其他VOC 的標準排放因子采用Guenther等[32]提出的推薦值1.5μg C/(g·h).各樹種的標準排放因子見表1.

表1 優勢樹種異戊二烯和單萜烯排放因子Table 1 Emission factors of isoprene and monoterpenes for dominant species in Beijing

1.3.2 葉生物量 本文參考方精云等[49]研究森林植被碳庫時采用的由植被蓄積資料推算生物量的方法,利用各類樹種蓄積量對其葉生物量進行估算.具體計算公式如下[50]:

式中: B為葉生物量, g; V為樹種蓄積量, m3; DT為樹干密度,106g/m3; PT為樹干生物量占喬木層總生物量的比例,%; PL為葉生物量占喬木層總生物量的比例, %.

森林蓄積資料來源于北京市第八次森林資源二類調查統計結果,內容包括北京市主要森林類型各林齡級的面積和蓄積量.依照張莉等[26]提出的方法,將北京市優勢樹種的樹干密度、主要森林類型各林齡級的PT和PL值分別分別整理為表2 和表3.

表2 主要樹種樹干密度(106g/m3)Table 2 Trunk densities of dominant species

表3 各森林類型不同林齡級的PL 和PT 值(%)Table 3 The PL and PT values of different forest types

1.3.3 氣象數據 本文使用的溫度和光合有效輻射的逐小時數據來源于北京城市生態系統研究站.

2 結果與討論

2.1 北京市BVOCs 排放清單

表4 按貢獻率排名給出了2015 年北京市森林優勢樹種BVOCs排放清單.可以看出,北京市森林源BVOCs 排放總量平均值為27.97×109g C/a,變化域值范圍為9.46～76.45×109g C/a.其中異戊二烯、單萜烯和OVOCs 平均排放量分別為21.00, 4.37 和2.60×109g C/a,對森林源BVOCs 總排放量的貢獻率分別為75.09%、15.62%和9.29%.

森林源BVOCs 排放存在非常大的樹種差異(表4).楊樹的BVOCs 排放量最大,為14.68×109g C/a,占BVOCs 排放總量的52.47%;櫟樹次之,排放量為7.20×109g C/a,貢獻率為25.73%.同時,楊樹和櫟樹也是主要的異戊二烯排放源,對異戊二烯排放的貢獻率分別為 63.16% (13.27×109g C/a)和25.92% (5.44×109g C/a).這一方面與北京地區楊樹和櫟樹擁有較大的蓄積量有關;另一方面,楊樹與櫟樹也具有較高的異戊二烯排放速率.對于單萜烯,排放量最大的樹種是油松(1.79×109g C/a),貢獻率為40.90%;其次為櫟樹(1.16×109g C/a),貢獻率為26.45%.

通過與前人的研究結果對比發現,不同地區之間BVOCs 各組分所占比例具有較大差異: Guenther等[51]在北美地區的研究顯示,異戊二烯排放量占BVOCs 總排放量的35%;宋媛媛等[52]對中國東部地區的調查發現異戊二烯、單萜烯和OVOC 的排放量所占比例分別為44.9%、31.5%和23.6%;劉巖等[53]對長江三角洲BVOCs 的研究顯示,異戊二烯、單萜烯和OVOC 分別占總排放量的37.3%、16.1%和46.6%.與上述研究相比,北京地區異戊二烯占BVOCs 總排放量的比例明顯更高,這主要是由于北京地區分布有較多的楊屬和櫟屬植物,二者均屬于高異戊二烯排放源,而且葉生物量較大.謝揚飏等[19]和Ren 等[21]對北京地區綠地植被BVOCs 的研究結果也表明,異戊二烯占總排量的比例分別為80.3%和74.0%,這與本文結果相似.

在過去的十幾年間,研究人員使用不同的數據源和算法對北京地區的BVOCs 排放量進行了估算,結果如圖1 所示.將本文基于樹種蓄積量的估算結果與以往的研究進行對比,可以看出本文估算值與Ren 等[21]對北京BVOCs 排放量的估算結果(33Gg C/a)較為接近.但是要高于Wang 等[20]和Ghirardo等[9]對1998, 2005 和2010 年BVOCs 排放量的估算值,低于Klinger等[34]、謝揚飏等[19]和Li等[54]對2000,2001 和2003 年排放量的估算值.

表4 北京市森林優勢樹種BVOCs 排放清單Table 4 BVOCs emission inventory of dominant species in Beijing

圖1 不同研究北京市BVOCs 排放量估算值比較Fig.1 Comparison of BVOCs emissions estimated by this study and previous studies

上述各研究之間估算結果相差較大的原因來自多個方面,如標準排放因子、葉生物量、植被組成、溫度和輻射強度等,任何一個模型輸入參數的變化都會導致估算結果的差異[17].已有的研究表明,不同研究中所使用的土地利用數據的不同是導致差異產生的一個主要原因,如Wang 等[20]對1998 年BVOCs 排放量進行估算所使用的森林面積為987km2,而Klinger 等[34]對2001 年BVOCs 估算使用的森林面積為2671km2,二者相差2.7 倍,這使得兩項研究得到的BVOCs 年排放值相差3 倍(分別為16 和48Gg C).事實上,從1990～2015 年北京市的森林覆蓋率增加了20.64%,從2005～2015 年森林蓄積量上升了405.8 萬m3,因而森林源BVOCs 排放量也隨之增加[21].此外,排放因子和生物量的差異也會在很大程度上影響估算結果. Li 等[55]、Li 等[54]的研究使用的排放因子和生物量為全國平均值,這使得二者對北京地區BVOCs 的估算結果要遠高于其它研究[55];謝揚飏等[19]和Ghirardo 等[9]同樣僅計算了城市綠地BVOCs 的排放,但謝揚飏等的估算結果遠高于Ghirardo 等,這主要也是因為二者所使用的排放因子和葉生物量數據的差異.Ren 等[21]的研究通過異速生長方程計算葉生物量,估算出北京市2015 年天然林BVOCs 排放量為28Gg,與本文基于森林蓄積量的估算結果基本持平;張倩等[56]使用傳統的基于葉面積指數的方法計算葉生物量,得到2015 年度北京市BVOCs 的排放總量為35Gg,略高于本文結果.

2.2 各齡級優勢樹種的BVOCs 貢獻

根據樹木生長年份可將其劃分為幼齡林、中齡林、近熟林、成熟林和過熟林.本文對不同齡級森林的BVOCs 排放情況進行分析,并按排放量高低給出北京市主要森林類型BVOCs 排放的齡級分配(圖2).可以看出,不同樹種各齡級貢獻情況有所差異,其中側柏和刺槐對BVOCs 排放量貢獻最大的是幼齡林;楊樹、櫟樹、樺木、落葉松和闊葉樹對BVOCs 排放量貢獻最大的是中齡林;油松貢獻率較大的是成熟林和中齡林.此外多數森林類型5 個齡級均有BVOCs 排放,但是部分森林因為生長年限短導致過熟林所占比例非常低,如櫟樹林、油松林、和樺木林等.

整體來看,北京地區中齡林對BVOCs 排放的貢獻率最大,占排放總量的39.14%;其次為幼齡林,貢獻率為23.05%;近熟林、成熟林和過熟林的貢獻率較小,分別為15.84%、14.59%和7.37%(圖3(a)).森林BVOCs 排放量與葉生物量在不同齡級林分中的分配情況較為一致(圖3),其中近熟林、成熟林和過熟林在BVOCs 排放量中所占比重略高于其在葉生物量中所占比重,可能是由于這3 個齡級森林中排放潛力高的樹種(本文中為楊樹)所占比例相較其它齡級更大導致[57].

圖2 北京市主要森林類型BVOCs 排放量的齡級分配Fig.2 Distribution of BVOCs emission among different age-class of main forest types in Beijing

圖3 不同齡級森林BVOCs 排放量和葉生物量分配情況Fig.3 Distribution of BVOCs emission and leaf biomass among different age-class forests

2.3 北京市BVOCs 排放的季節變化

由圖4 可以看出,北京森林源BVOCs 排放呈現出明顯的季節變化特征.在春季,地面溫度和光照強度逐漸上升,植被生長,BVOCs 排放量逐漸增加.但由于樹葉仍處于發芽生長階段,葉面積較小、酶活性也較低,春季排放量僅3.51×109g C,占全年的12.55%(表5).夏季是BVOCs 排放的主要季節,這一時期溫度和太陽輻射達到1a 中的最大值,樹葉也達到成熟期,葉生物量大、酶活性高,各樹種的BVOCs 排放速率也因此增加,排放量達到21.67×109g C,占全年總排放量的77.48%.秋季氣溫逐漸降低,樹葉由成熟變為衰老,落葉樹種葉片相繼脫落,因此秋季BVOCs 排放量驟減,僅2.73×109g C,占全年排放總量的9.76%.冬季氣溫及光照強度均降至全年最低,僅有部分常綠樹種以較低的速率排放VOCs[20],排放量0.06×109g C,只占全年排放總量的0.21%.

圖4 不同月份北京森林源BVOCs 排放量Fig.4 Monthly BVOCs emissions of forest in Beijing

本研究分析的3 類BVOCs 中,在不同季節排放量相差最大的是異戊二烯.異戊二烯對BVOCs 總排放量的貢獻從1～7 月呈上升趨勢,從7～12 月呈下降趨勢,夏季占總排放量的77.30%,而冬季僅占5.27%(表5),說明異戊二烯排放受季節變化的影響要高于其他2 類.原因主要有兩方面,一是根據Guenther等[2,23]的研究,異戊二烯屬于葉綠體排放類型,其排放量受溫度和太陽輻射的雙重影響;而單萜烯和OVOCs 為防御專屬組織排放類型,其排放主要受溫度影響,因此隨季節變化的波動要小于異戊二烯.二是異戊二烯主要來源于落葉闊葉林,其排放量隨著樹葉的春生、夏長、秋衰、冬落而表現出巨大的季節差異;而單萜烯主要來源于針葉林,樹葉一年四季保持常綠,全年均有排放.

表5 北京森林源BVOCs 各季節排放量及比重Table 5 Seasonal variation of BVOCs emissions and proportions from forest in Beijing

2.4 北京市BVOCs 排放的日變化

圖5 北京地區1 月和7 月氣溫、光合有效輻射和BVOCs排放日變化規律Fig.5 Diurnal variation of air temperature, photosynthetic active radiation and BVOCs emissions in January and July

如圖5 所示,異戊二烯的排放量隨溫度和太陽輻射強度的變化而改變,由于夜晚無光照,因此排放僅出現在白天,表現為明顯的單峰曲線(圖5(b),(c)).1月份異戊二烯的排放時間在07:00～17:00,上午隨著溫度和輻射強度的增加,排放量持續上升,在中午13:00 達到排放峰值,之后開始下降; 7 月份光照時間長,異戊二烯的排放時間在04:00～20:00,排放峰值出現在14:00 點前后.

單萜烯和OVOC 的排放主要受溫度的影響,全天均有排放(圖5(b),(c)).單萜烯和OVOC 呈現出相似的日變化規律:1 月份的最低排放值出現在07:00,因為此時溫度最低,之后隨著白天溫度的回升排放量逐漸增加,并在15:00左右達到排放峰值; 7月份的最低排放值出現在05:00,最大排放值出現在14:00～16:00.冬季和夏季排放的日變化趨勢整體相似,但最低排放值出現的時間不同,這主要是由于冬夏季日最低溫出現的時間不同導致.此外,本文顯示異戊二烯排放的日變化幅度遠大于單萜烯和OVOC,這一特點在Chang 等[58]對杭州地區的研究和Li 等[55]對全國范圍的研究中也有體現,說明異戊二烯更易受到環境因子的影響,對太陽輻射具有較高的敏感性.

2.5 O3 和SOA 生成潛勢及其對空氣質量的影響

北京森林源BVOCs 對OFP 的總貢獻量為240.51×109g.其中異戊二烯的OFP貢獻量為222.86×109g,占總量的92.66%,是主要的貢獻者;單萜烯的OFP 貢獻量為17.65×109g,占總量的7.34% (表6).森林源 BVOCs 對 SOA 生成潛勢的總貢獻量為1.73×109g,其中異戊二烯和單萜烯的貢獻量分別為0.42×109g 和1.31×109g,分別占SOA 生成總量的24.26%和75.73%,單萜烯是主要的貢獻者.各主要森林類型的O3和SOA 生成潛勢具體見表6, O3生成潛勢最大的樹種為楊樹, SOA 生成潛勢較大的樹種為楊樹、櫟樹和油松.

表6 北京市森林源BVOCs 的O3 和SOA 生成潛勢Table 6 Ozone and SOA formation potential of BVOCs from forest in Beijing

運用“箱式”模型的估算結果顯示,森林源BVOCs 對夏季O3污染具有顯著貢獻(圖6(a)).BVOCs 排放會導致全市范圍內大氣年均O3濃度上升9.01μg/m3,由于BVOCs 排放具有明顯的季節動態,其引起的O3濃度增加也表現出明顯的季節變化.從圖6(a)中可以看出,在排放量最高的7 月份,BVOCs 排放導致O3濃度上升了31.30μg/m3,占環境中O3濃度的17.7%.在其他的研究中,也發現夏季BVOCs對O3形成具有相對較高的貢獻率,如Li等[59]對長江三角洲的研究和Ren 等[21]對北京城區的研究發現,在夏季,BVOCs 排放貢獻了大氣O3濃度的19%～30%.此外,雖然相關的研究表明北京地區AVOCs 的排放量是BVOCs 的6～16 倍,但是AVOCs的OFP 只有BVOCs 的3～8 倍[22],這說明BVOCs 具有更高的反應活性,因此其對O3污染的貢獻不容忽視[60].除此以外, BVOCs 排放還會造成全市大氣年均SOA 濃度增加0.94μg/m3,占2015 年北京市年均PM2.5濃度的1.2%.圖6(b)顯示,BVOCs 對大氣PM2.5污染的貢獻相對較低,這與前人的研究結果一致[9,21].

圖6 BVOCs 排放對環境O3 和PM2.5 污染的貢獻Fig.6 The contribution of BVOCs emissions to ambient ozone and PM2.5 pollution

近年來北京市大力開展綠色生態空間建設,森林覆蓋率逐年上升,同時采取了一系列措施嚴格控制人為源VOCs 的排放[61-62],因而植物源VOCs 對大氣質量產生的影響日漸突出.由于每種BVOCs 的反應速率和途徑不同,其對O3和PM2.5形成的影響也存在差異.一般來說,異戊二烯排放量高的樹種對O3污染的貢獻較大,而單萜烯排放量高的樹種對PM2.5污染的促進作用較強[63].在北京市主要的綠化樹種中,異戊二烯排放潛力最高的是楊樹,其次是櫟樹和國槐等,強單萜烯排放樹種主要是油松,調整優化樹種組成也許是控制BVOCs 排放最可行的途徑之一.Ren 等[21]的研究結果顯示,通過采取主動管理措施(如改變樹種組成等),到2050 年可以減少61%的BVOCs 排放和50%的由BVOCs 引起的健康損害.因此,未來進行景觀規劃及城市綠化時,應充分考慮樹種的VOCs 排放能力,這將有助于改善環境空氣質量并減少對人體健康的損害.

2.6 不確定性分析

在本文估算過程中,森林BVOCs 排放量及其對區域空氣質量影響的不確定性主要來自以下幾個方面.

環境校正因子.在計算溫度影響因子時,本研究使用環境空氣溫度代替葉表面溫度.由于葉片溫度通常高于氣溫,尤其在闊葉樹中這一現象更為明顯,因此森林源BVOCs 排放量可能會被低估[64].此外,光合有效輻射在森林冠層內部的傳輸過程中會出現衰減,因此位于冠層下方葉片的光合有效輻射影響因子較上方葉片會偏低[23].

葉生物量的計算.由于目前對植被葉生物量的直接調查較少,本文中葉生物量由樹種蓄積量、樹干密度、樹葉和樹干與喬木層生物量比重等植物學參數推算所得,在用資料反推的過程中必然存在一定誤差,這些參數的準確度也會影響葉生物量的估算結果.

BVOCs 對區域空氣質量的影響.在估算BVOCs對空氣質量的影響時,使用了一個簡化的“箱式”模型,這個模型假設BVOCs 的排放只對箱內的空氣質量產生影響.但事實上, O3及其前體物等會隨著空氣的流動被輸送到下風向地區[65],導致估算結果存在誤差.

3 結論

3.1 基于樹種蓄積量的北京市森林源BVOCs 排放量平均值為27.97×109g C/a,變化域值范圍為9.46～76.45×109g C/a,其中異戊二烯、單萜烯和OVOCs平均排放量分別為21.00, 4.37和2.60×109g C/a,分別占排放總量的75.09%、15.62%和9.29%.不同樹種間BVOCs 排放量差異較大,楊樹和櫟樹是主要的異戊二烯排放源,貢獻率分別為63.16%和25.92%,油松是主要的單萜烯排放源,貢獻率為40.90%.將多項研究進行對比發現,基于森林蓄積量估算的BVOCs 排放量與其它研究方法的結果可以較好地匹配.

3.2 北京森林源BVOCs 排放呈現明顯的季節變化,排放量夏季(77.48%) ＞ 春季(12.55%) ＞秋季(9.76%)＞ 冬季(0.21%),僅7 月份排放量就占全年的28.83%.異戊二烯、單萜烯和OVOC 呈現不同的日變化規律:異戊二烯受溫度和太陽輻射的影響,僅白天排放,日變化幅度大,峰值出現在13:00～14:00 時;單萜烯和OVOC 排放主要受溫度影響,晝夜均有排放,日變化相對平緩,峰值出現在15:00 時前后.

3.3 森林源BVOCs 對O3生成潛勢的貢獻量為240.51×109g,其中異戊二烯的貢獻量為222.86×109g,占總量的92.66%,是主要的貢獻者; BVOCs 對SOA生成潛勢的總貢獻量為1.73×109g,異戊二烯和單萜烯的貢獻量分別為0.42×109g和1.31×109g,分別占總量的24.26%和75.73%. BVOCs 排放會導致全市范圍內大氣年均O3濃度上升9.01μg/m3, SOA 濃度增加0.94μg/m3. BVOCs 對夏季O3污染具有顯著貢獻,未來應通過優化樹種組成等管理措施減少BVOCs排放,以改善環境空氣質量.