吳建國,徐天瑩 (1.中國環(huán)境科學研究院,北京 100012；2.甘肅農(nóng)業(yè)大學林學院,甘肅 蘭州 734000)

大氣揮發(fā)性有機物(VOCs)中很多成分具有高度化學活性,極易與大氣中NO和NO2、氧化劑和 OH自由基等發(fā)生反應,在對流層臭氧光化學過程中具有重要作用[1-3].另外, VOCs為臭氧形成的重要前體物之一,對氣候變暖有間接的影響[3-5].準確確定 VOCs排放特征及影響因素, 對科學認識對流層臭氧光化學過程及大氣環(huán)境變化和氣候變暖有重要的理論意義.大氣 VOCs中總量 90%以上的非甲烷碳氫化合物(NMHC)來自陸地生物圈排放,單萜烯類是植物排放最主要的 VOCs種類之一,是對流層中最主要的非甲烷烴種類之一,約占植物VOCs釋放總量的11%[5-9].植物排放單萜烯速率隨葉片溫度增加而增加[5,10].全球氣候正經(jīng)歷以變暖為主要特征的顯著變化[11].隨著全球氣候變化,植物葉溫將隨之改變,進而將使植物單萜烯排放速率改變[9,12].這些變化將對大氣化學成分和大氣環(huán)境產(chǎn)生一定的影響[5,13].因此,系統(tǒng)分析氣候變化對植物單萜烯排放速率的影響,對準確確定氣候變化下植物VOCs排放特征及大氣環(huán)境變化具有重要的理論和現(xiàn)實意義.

國際上對植物 VOCs排放研究廣泛展開,包括發(fā)展植物單萜烯排放測定方法[14],分析植物單萜烯排放[15-20],或測定不同植物單萜烯排放并計算清單[21],或比較不同植物下土壤單萜烯排放速率[22],或分析植物單萜烯排放季節(jié)或年際變化[23-27].另外,一些研究分析影響植物單萜烯排放的因素,如分析溫度升高[28]或有效水變化[29]或光照改變[30]及葉片發(fā)育[31]和植物氣孔[32]及植物地下作用對 VOCs排放影響[33],或草食動物和氣候因素等[34-35]及霉病[36]等對植物 VOCs排放影響等.近年來,氣候變化對植物單萜烯排放影響頗受關(guān)注,包括分析氣候變化對區(qū)域植被 VOCs排放影響[37],特別是氣候變暖對極地植物和亞極地植物 VOCs排放的影響[38-41],以及水分變化對植物 VOCs排放的影響[42-47]等.另外,也觀測了氣候變暖與 CO2濃度或 O3濃度增加對植物解剖結(jié)構(gòu)、生理過程及VOCs排放速率的共同影響[48-57].此外,也模擬分析了氣候變化引起植被改變[58]或種類組成變化對單萜烯排放總量的影響[59-60],以及景觀和氣候變化對植物單萜烯排放影響[61]等總體上,這些研究極大提高了對氣候變化對植物VOCs排放影響的認識,但研究還集中極地或亞極地或北美和歐洲部分地區(qū)氣候變化對個別植物單萜烯排放影響,或分析氣候變化引起的植被組成變化對單萜烯排放影響方面,對其它地區(qū)氣候變化對植物單萜烯排放的影響研究還不多.我國對植物 VOCs排放研究也有許多報道,包括測定亞熱帶竹林[62]、長白山森林[63]、西雙版納三葉橡膠[64]、沈陽典型樹種[65]、黃連木[66]、北京植物[67]單萜烯排放速率及南京紫金山植物群落VOCs組成[68],以及亞熱帶城鄉(xiāng)復合系統(tǒng) BVOC排放[69]或中國天然源 VOCs排放清單[70-71]等.這些研究對認識我國植被單萜烯排放規(guī)律有一定意義,但氣候變化對植物單萜烯排放影響方面還沒有研究報道.

油松(Pinus tabuliformis)為松科針葉常綠喬木,是我國特有樹種,分布廣泛.油松是典型的單萜烯排放植物[72-75],目前已對這種植物的單萜烯排放特征進行了一些觀測分析[72-75],但氣候變化對這種植物單萜烯排放影響還沒有研究報道.為此,本文以全球氣候模式 NorESM1-M 產(chǎn)生的RCP2.6,RCP4.5,RCP6.0和RCP8.5氣候變化情景數(shù)據(jù)和植物 VOCs排放計算模型,模擬分析了未來氣候變化對油松單萜烯排放速率的影響,希望為科學認識氣候變化對植物 VOCs排放影響提供一定的依據(jù),為相關(guān)的研究提供一定的參考.

1 材料與方法

1.1 研究區(qū)

為了集中分析氣候變化對油松單萜烯排放速率的影響,本文選擇了油松分布集中區(qū)域—山西省沁源縣東部的太岳山林區(qū)中部馬泉林場(36.42～36.83N,112.25～112.55E)為例,該區(qū)氣候類型屬暖溫帶半干旱大陸性季風氣候,四季分明,近幾十年氣候變化趨勢明顯[77-78].

1.2 油松單萜烯排放計算

基于文獻確定油松在標準狀態(tài)下單萜烯排放速率為 4.3μg/(g·h)C[76].遵循植物排放單萜烯的一般計算模型,油松單萜烯排放通量按照公式(1)進行計算.

式中：E為植物單萜烯和其它排放通量,μg/(g·h);Es為標準條件下(T=303K) 植物單萜烯排放量,μg/(g·h);δ為考慮長期排放活性的控制因子,由于油松為常綠樹種,δ取值為1;B為葉片生物量密度,由于本文是模擬計算單位重量油松葉生物量單萜烯排放速率,生物量密度取單位重量值為 1g;γ為活性排放因子,它為葉溫的函數(shù),具體計算利用公式(2).

式中：Ts為標準條件下葉溫(303K);T為葉溫(K);β一般取0.09/K.

以上公式計算獲得單萜烯和其它排放速率是每小時值,且葉溫也是小時值. 由于目前缺乏葉溫的觀測資料,用能量守衡方法計算葉溫又非常復雜,在考慮到葉溫和氣溫平均差別不是很大的情況下,本文參照文獻[76,79]中做法,直接用氣溫代替葉溫計算葉片每小時葉溫.計算溫度小時值方法具體見1.3.2溫度日變化計算部分.

1.3 氣候變化對油松單萜烯排放的影響分析

1.3.1 氣候變化情景數(shù)據(jù) 采用全球氣候模式NorESM1-M產(chǎn)生的RCP2.6,RCP4.5,RCP6.0和RCP8.5情景數(shù)據(jù)(RCP2.6,RCP4.5,RCP6.0和RCP8.5分別代表大氣中輻射強迫為 2.6,4.5,6.0,8.5W/m2).在 RCP2.6情景中,未來人類將采取大力減排措施,使輻射強迫在21世紀達到頂峰并下降,到 2100年大氣中 CO2濃度將達到 421×10-6,輻射強迫達2.6W/m2,全球升溫幅度達0.3～1.7℃;RCP8.5是最高的溫室氣體排放情景,到 2100年大氣中 CO2濃度將達 936×10-6,輻射強迫將持續(xù)上升,到2100年將達8.5W/m2,全球升溫幅度將達2.6～4.8℃; RCP4.5和RCP6.0情景介于RCP2.6和RCP8.5之間,在RCP4.5情景下2100年輻射強迫達 4.5W/m2,升溫幅度將達 1.1～2.6℃,在RCP6.0情景下2100年輻射強迫達6.0W/m2,升溫幅度將達1.4～3.1℃[81-84].這些情景數(shù)據(jù)是由中國科學院地理與資源研究所從 ISI-MIP(多領(lǐng)域間影響模型比較計劃)獲取的經(jīng)過空間降尺度,偏差校正的RCP情景下的全球氣候模式模擬數(shù)據(jù),由中國農(nóng)業(yè)科學院進行了訂正[85].

1.3.2 溫度日變化計算 氣候變化情景中溫度是日值,在計算油松單萜烯排放中按照小時值進行計算,需要把溫度日值轉(zhuǎn)化為溫度每小時值.參考文獻[86],白天每小時溫度值按公式(3)計算.

式中：Tmin為最低溫度,Tmax為最高溫度;s(t)為計算溫度變化函數(shù),按公式(4)計算.

式中：π為常數(shù);th為一天太陽時,按公式(5)計算;LSH為最大太陽高度時刻,一般取值 12h;DL為白晝長度,按照公式(8)計算;P為最高氣溫與最大太陽高度出現(xiàn)的時差,一般取值2.

式中：bj為北京時,J為地理經(jīng)度;η為時差,參考文獻[87],按公式(6)計算.

式中：Q按公式(7)計算.

式中：DAY為1年中日序, 每年1月1日DAY=1,1月2日DAY=2,依次類推,12月31日DAY=365.

式中：L為地理緯度;δ中太陽赤緯,參考文獻[87],按公式(9)計算.

白天時刻按照公式(10)判斷.

夜間時刻按照公式(11)判斷.

夜間每小時溫度值按公式(12)計算.

式中： tmin為日最低氣溫;th為一天太陽時,按公式(5)計算;ts為日落時間,取值;η為夜間長度,取值h=24-DL;τ為時間系數(shù),取值4h.

1.3.3 氣候變化影響下油松單萜烯排放速率計算 把基準情景和未來氣候變化情景下每年逐日逐時溫度值,代入公式(1)～(2),計算油松單萜烯每年逐日逐時的排放速率,在此基礎(chǔ)上,再進行24h的累加,計算油松單萜烯每年逐日排放速率;基于油松單萜烯的每年逐日排放速率,經(jīng)過全年逐日累加,再計算油松單萜烯每年排放速率值.通過比較基準情景和未來氣候變化情景下不同時段油松單萜烯日排放速率及按年計算的排放速率的差異,分析氣候變化對油松單萜烯排放速率影響.所有計算都是在Visual Fortran6.6平臺中經(jīng)過編程進行.

2 結(jié)果與分析

2.1 氣候變化特征

以1951～2000年為基準,圖1表明,2001～2050年太岳山林區(qū)中部年均氣溫呈現(xiàn)明顯上升趨勢,在4種情景下(RCP2.6,RCP4.5,RCP6.0和RCP8.5)變化幅度不同.到2050年,相比基準年,太岳山林區(qū)中部年均氣溫在 4種情景下分別增加 1.23,1.28,1.23,1.52℃.相比1951～2000年,在2001～ 2025年時段,太岳山林區(qū)中部年均氣溫在4種情景下分別增加 0.90,0.71,0.84,0.99℃;到 2026～2050 年時段,太岳山林區(qū)中部年均氣溫在 4種情景下分別增加1.56,1.85,1.61,2.04℃.

圖1 2001～2050年相比基準情景太岳山林區(qū)中部平均氣溫變化趨勢Fig.1 Changes in annual mean air temperature in 2001～2050 compared with baseline periods (1951～2000) in middle of Taiyue Mountains

以 1951～2000年為基準,圖 2表明,2001～2050年太岳山林區(qū)中部年降水量呈較大波動,在RCP2.6,RCP4.5,RCP6.0和RCP8.5四種情景下變化不同.到 2050年,相比基準年,太岳山林區(qū)中部年降水量在 4種情景下分別減少11.61mm,增加1.38,15.84,11.71mm.相比基準年,到 2001～2025年時段,太岳山林區(qū)中部的年降水量在4種情景下分別減少23.77,2.40mm,增加23.54,16.56mm; 到 2026～2050年,太岳山林區(qū)中部年降水量在4種情景下分別增加0.55,5.16,8.13,6.86mm.

以1951～2000年為基準,圖3表明2001～2050年太岳山林區(qū)中部輻射強度呈較大的波動,并呈下降趨勢,在RCP2.6,RCP4.5,RCP6.0和RCP8.5 4種情景下變化程度不同.相比基準,到 2001～2050年在4種情景下太岳山林區(qū)中部的輻射強度分別平均減少0.10,0.18,0.35,0.25MJ/m2(圖3).相比基準年,在2001～2025年太岳山林區(qū)中部的輻射強度在4種情景下平均減少幅度分別為 0.28,0.33,0.32,0.37MJ/m2;相比基準年,到 2026～2050年太岳山林區(qū)中部輻射強度在4種情景下平均減少幅度分別為0.08,0.03,0.37,0.14MJ/m2.

圖2 2001～2050年相比基準情景太岳山林區(qū)中部年降水量變化趨勢Fig.2 Changes in annual precipitation amount in middle of Taiyue Mountains in 2001～2050 compared with baseline periods (1951～2000)

圖3 2001～2050年相比基準情景太岳山林區(qū)中部太陽總輻射變化趨勢Fig.3 Changes in total solar radiation intensity in the middle of Taiyue Mountains in 2001～2050 compared with baseline periods (1951～2000)

2.2 氣候變化影響下油松單萜烯日排放速率變化

圖4 為1951～2000年基準情景和2001～2050年RCP2.6,RCP4.5,RCP6.0和RCP8.5 4種情景下50a平均的油松單萜烯日排放速率值.圖4顯示,在基準情景和 4種未來情景下,油松單萜烯全年日排放速率變化較大.在1～210d油松單萜烯日排放速率呈現(xiàn)上升趨勢, 在210～365d呈下降趨勢.按照對應的月份,油松單萜烯日排放速率在 5～9月較高,在1～4月和9～12月較小.在基準情景和4種未來情景下,油松單萜烯全年的日排放速率平均 值 分 別 為 21.22,23.59,23.47,23.26,23.85μg/(g·d).顯然,油松單萜烯全年日排放速率平均值在未來氣候變化情景下比基準情景下要高約 2μg/(g·d),并且在RCP8.5情景下最高.

圖4 基準和不同氣候變化情景下50a平均的油松單萜烯日排放速率Fig.4 Daily Monoterpenes emission rate from leaves of Pinus tabuliformis under baseline and different climate change scenarios

圖5 為2001～2050年RCP2.6,RCP4.5, RCP6.0和RCP8.5情景下與1951～2000年基準情景下50a平均的油松單萜烯日排放速率差值.如圖 5所示,油松單萜烯日排放速率在4種未來情景下與基準情景下差值變化很大.在 1～95d,油松單萜烯日排放速率在未來氣候變化情景與基準情景下的差異較小;在 95～295d,油松單萜烯日排放速率在未來氣候變化情景與基準情景下差值呈現(xiàn)較大波動趨勢;在 296～365d,油松單萜烯日排放速率在不同氣候變化情景與基準情景下差異不大.總體上,油松單萜烯日排放速率在未來氣候變化情景與基準情景下的差異平均在 2μg/(g·d)以上.另外,油松單萜烯日排放速率在 RCP2.6,RCP4.5,RCP6.0情景與基準情景下的差異相對較小(分別平均為2.37,2.24,2.03μg/(g·d)),在RCP8.5情景與基準情景下差異較大(平均 2.63μg/(g·d)).

圖 6為 2001～2050年 4種未來情景相比1951～2000年基準情景下50a的油松單萜烯日排放速率相對變化百分率.相比基準情景,在 4種未來情景下,在 1～190d,油松單萜烯日排放速率相對變化百分率較高(相比基準情景平均增加12%～14%以上);在 191～315d,油松單萜烯日排放速率相對變化百分率較小(相比基準情景平均增加9%～13%以上);在316～365d,油松單萜烯日排放速率相對變化百分率在12%～18%以上.另外,與基準情景相比,在 4種未來情景下,油松單萜烯日排放速率變化幅度平均分別為 11.68%、12.37%、11.45%和13.99%, 在RCP8.5情景下較大.

圖5 不同氣候變化情景與基準情景下50a平均的油松單萜烯日排放速率差異Fig.5 The difference of daily monoterpenes emission rate from leaves of Pinus tabuliformis under baseline and different climate change scenarios

圖6 相比基準情景不同氣候變化情景50a平均油松單萜烯日排放速率變化Fig.6 The change in daily monoterpenes emission rate from leaves of Pinus tabuliformis under different climate change scenarios compared to the baseline scenario

2.3 氣候變化影響下油松單萜烯年排放速率變化

圖 7顯示,與基準情景下相比,在 RCP2.6,RCP4.5,RCP6.0和 RCP8.5 4 種情景下,2001～2050年油松單萜烯年排放速率呈現(xiàn)增加趨勢.在基準情景下,油松單萜烯年排放速率平均為7746μg/(g·a);在 4 種未來情景下,50a油松單萜烯年排放速率分別平均為 8610,8565,8488,8705μg/(g·a).總體上,未來氣候變化情景下,50年油松單萜烯年排放速率平均比基準情景下高出約 1000μg/(g·a)以上,在 RCP8.5 情景下油松單萜烯年排放速率最大.相比基準情景,2001～2050年油松單萜烯年排放速率在4種未來情景下平均增加幅度為 9%～13%,50a在 4種未來情景下分別增加 11.15%、10.57%、9.58%和 12.37%.顯然,在 RCP8.5情景下變化幅度最大.

圖7 不同氣候變化情景與基準情景下油松單萜烯年排放速率變化Fig.7 Changes in annual emission rate of monoterpenes from leaves of Pinus tabuliformis under different climate change scenarios compared to the baseline scenario

3 討論

本研究說明,未來氣候變化影響下山西太岳山林區(qū)中部氣溫呈上升的趨勢,這與全球趨勢比較一致[11].另外,本研究也表明,未來氣候變化下山西太岳山林區(qū)中部年降水量變化波動明顯,總體上呈現(xiàn)增加或下降趨勢.同時,氣候變化下山西太岳山林區(qū)中部輻射強度波動也比較大,呈現(xiàn)下降趨勢,這些趨勢與全球氣候變化預測中一致[11].當然,這里主要反映了山西太岳山林區(qū)中部區(qū)域的變化趨勢,不同區(qū)域?qū)⒂胁煌淖兓卣?

一些研究表明未來氣候變化下將使植物單萜烯排放速率增加.如Lindwall等[12]發(fā)現(xiàn),夏天變暖4℃使極地凍原VOCs排放增加280%,積雪增加影響下VOCs排放變化卻不大;Kramsh?j等[41]發(fā)現(xiàn),氣候變暖使極地植物 VOCs排放增加,變暖使生態(tài)系統(tǒng) VOCs排放速率增加 260%,植物VOCs排放速率增加 90%,減少太陽輻射使系統(tǒng)VOCs排放速率下降69%,植物VOCs排放速率下降61%～65%; Faubert等[60]發(fā)現(xiàn)氣候變化使亞北極苔原 VOCs排放加倍,Rinnan等[40]發(fā)現(xiàn)氣候變化引起灌木VOCs排放增加,Vedel-Petersen等[39]發(fā)現(xiàn)氣候變化導致植物 VOCs排放增加.還有研究發(fā)現(xiàn)溫度升高使植物 VOCs排放增加[28]或區(qū)域氣候變化使VOCs排放增加[13].本研究也表明,氣候變化下油松單萜烯日和年排放速率都呈現(xiàn)增加趨勢,其中日增加幅度為 9%～18%,年增加幅度為 9%～13%.這些增加幅度沒有前述文獻[12,41,60]中高,這主要可能因為這些文獻中基于控制試驗結(jié)果,相比本研究模擬時間短,而本研究以 1950～2000年為基準情景,模擬未來 2001～2050年50a的變化趨勢,并且文獻中試驗增溫幅度比本研究模擬情景增溫幅度要大,試驗中植物與本研究選擇的油松生理代謝等也不同.另外,本研究結(jié)果也表明,氣候變化下油松單萜烯排放存在明顯日期和氣候變化情景的差異. 如油松單萜烯日排放速率 1～210d呈上升趨勢,210～365d呈下降趨勢, 在RCP8.5情景下較高;未來氣候變化情景與基準情景下油松單萜烯日排放速率1～95d差異較小,95～295d波動較大,296～365d 差異不大;油松單萜烯日排放速率相對變化百分率1～190d較高, 191～315d變化較小, 316～365d變化在12%～18%以上,在RCP8.5情景下變化較大.這可能與研究區(qū)氣溫變化和植物生理活動時間差異有關(guān).在本研究區(qū),在基準和未來氣候變化情景下,全年中 1～210d氣溫呈上升趨勢,在 211～365d氣溫呈現(xiàn)下降趨勢,植物單萜烯日排放速率也隨之改變,但在1～90d和296～365d氣溫都不是很高,并且也不是植物最活躍生理活動期,而在 91～295d氣溫變化范圍大,同時也是植物最活躍的生理活動期,并且相比基準情景,未來氣候變化情景下春夏季氣溫增加相對明顯,加之在春夏季植物生理活動活躍,所以使油松單萜烯日排放速率在基準和未來氣候變化情景下差異增大.另外,從氣候變化下溫度增加相對程度來看,在春季和冬季氣溫增加相對幅度更高,使油松單萜烯日排放速率相對變化百分率在1～190d和315～365d較高,在191～315d變化較小.這些結(jié)果意味著未來氣候變化強度和時間差異將導致植物單萜烯排放速度改變程度和時間差異.但也有研究發(fā)現(xiàn),氣候變化對植物 VOCs排放影響不確定.如 Schollert等[80]研究發(fā)現(xiàn),氣候變化改變了極地植物解剖結(jié)構(gòu),但對植物 VOCs排放卻沒有影響.說明植物發(fā)育過程對植物 VOCs排放影響較大,對此需要在以后研究中進行考慮.

一些研究認為氣候變化下植物單萜烯排放改變與植物生物量改變有關(guān),氣候變化下植物生物量增加導致單萜烯排放增加[37].本研究沒有考慮氣候變化對油松生物量的影響而導致對單萜烯排放的改變.如果考慮氣候變化對油松生物量的影響,可以推斷氣候變化下油松單萜烯排放增加程度將更明顯.這些影響需要在以后研究中進一步考慮.同時,氣候變化下植物單萜烯排放速率改變與植物生理活動改變有關(guān).如研究發(fā)現(xiàn)在高溫和CO2濃度升高下植物VOCs排放與光合作用和葉綠素結(jié)構(gòu)關(guān)系密切[54]或植物單萜烯排放與光合作用,氣孔通導性關(guān)系密切[32].另外,氣候變化下植物單萜烯排放速率改變與植物解剖結(jié)構(gòu)變化也關(guān)系密切.如研究發(fā)現(xiàn)溫度升高和O3濃度變化下植物 VOCs排放與葉片結(jié)構(gòu)關(guān)系密切[55]或植物 VOCs排放與植物地下結(jié)構(gòu)變化關(guān)系密切[33],氣候變化下植物葉片發(fā)育不同使植物VOCs排放季節(jié)變化[31].但也有研究發(fā)現(xiàn)氣候變化雖然引起植物解剖結(jié)構(gòu)改變,但對 VOCs排放沒有影響[88];還有研究發(fā)現(xiàn)氣候變暖和積雪增加使植物光合作用速率下降,但對植物 VOCs排放沒有影響,主要是因為葉片結(jié)構(gòu)改變[50].本研究中,假設氣候變化下油松葉片生理活動與解剖結(jié)構(gòu)沒有變化,這些假設是否正確有待檢驗.

一些研究發(fā)現(xiàn)水分改變對植物 VOCs排放影響明顯[29],氣候變化引起的干旱對植物 VOCs排放影響較大[38,43,45],干旱使植物 VOCs排放增加[42,44],或洪水[36]或干旱再加水[46]或降水增加使植物 VOCs排放增加[23].同時,一些研究發(fā)現(xiàn)大氣CO2濃度變化[49,51]及溫度和CO2濃度變化[48,52-54]或 O3和 CO2濃度變化[56-57]對植物單萜烯排放都有一定影響.另外,土地利用、CO2濃度、植被組成和氣候變化對單萜烯排放有不同影響[5,61].此外,植物霉變[36]或昆蟲取食[34,57]或氣候和昆蟲取食共同對植物單萜烯排放產(chǎn)生一定的影響[35].本研究主要考慮了溫度變化對油松單萜烯排放的影響, 沒有考慮其它因素變化對油松單萜烯排放的影響.事實上,氣候變化下隨著降水和蒸發(fā)的改變將引起水分的變化,同時大氣中CO2和O3濃度也將改變,并且將間接引起霉變或昆蟲取食改變,這些變化將不同程度直接或者間接影響植物 VOCs排放.在本研究中,如果考慮水分或其它因素變化對油松單萜烯排放的影響,可以推斷油松單萜烯排放變化將可能更加復雜,這需要在以后研究中進一步詳細分析.

氣象因素的誤差可能對植物單萜烯排放模擬結(jié)果造成誤差[89-91].在目前模擬分析植物單萜烯排放中,主要考慮了溫度因子.雖然植物葉片溫度與氣溫有關(guān),但存在一定差異.在本文中,參考相關(guān)文獻[76-77],忽略了葉片溫度和氣溫差異,這可能對研究結(jié)果造成一些誤差.但從分析氣候變化對植物單萜烯排放影響宏觀趨勢來看,這些差異也可以忽略.

需要指出的是,本研究只采用了單一氣候模式模擬的氣候變化情景結(jié)果,不同氣候模式下結(jié)果可能有一定差異,需要在以后研究中進一步分析.另外,一些研究發(fā)現(xiàn)氣候變化引起植被變化導致 VOCs排放變化[58]或氣候變化改變植物組成使 VOCs排放改變[60].而本研究只考慮了在山西省太岳山分布的油松,沒有考慮其他區(qū)域油松單萜烯排放特征. 同時,本研究也沒有考慮氣候變化可能導致油松分布范圍改變而對單萜烯排放影響.這可能使本研究結(jié)果存在一些不足,需要在以后研究中加強這些分析.盡管如此,本研究利用了目前新的氣候變化情景數(shù)據(jù),模擬分析了未來氣候變化對油松單萜烯排放速率影響,對科學認識氣候變化對植物 VOCs排放及大氣環(huán)境的影響還是有一定的參考意義.

4 結(jié)論

4.1 在氣候變化影響下,山西太岳山中部年均氣溫呈明顯上升趨勢,年降水量和日輻射強度都呈下降或者增加趨勢,波動較大.

4.2 油松單萜烯日排放速率在未來氣候變化情景與基準情景下差值在1～95d和296～365d較小,在 95～295d波動較大.油松單萜烯日排放速率相對變化百分率在 1～190d增加 12%～14%以上,在191～315d 增加 9%～13%以上,在 316～365d 增加12%～18%以上,在RCP8.5情景下較大.

4.3 相比基準情景,在未來氣候變化下油松葉片單萜烯年排放速率和變化程度呈增加趨勢,在氣候變化情景下油松單萜烯年排放速率平均增幅在9%～13%,在RCP8.5情景下最大.

[1]Collins W J, Derwent R C, Johnson C E,etal. The oxidation of organic compounds in the troposphere and their global warming potentials [J]. Climatic Change, 2002,52：453-479.

[2]Day M C, Pandis S N. Predicted changes in summertime organic aerosol concentrations due to increased temperatures [J].Atmospheric Environment, 2011,45(36)：6546-6556.

[3]唐孝炎,張遠航,邵 敏,等.大氣環(huán)境化學 [M]. 第 2 版.北京：高等教育出版社, 2010：74-77.

[4]Fehsenfeld F, Callvert J, Fall R, et al. Emissions of volatile of organic compounds from vegeration and the implications for atmospheric chemistry [J]. Global Biogeochemical Cycles, 1992,6(4)：389-430.

[5]Laothawornkitkul J, Taylor J E, Paul N D, et al. Biogenic volatile organic compounds in the Earth system [J]. New Phytologist,2009,183(1)：27-51.

[6]Guenther A, Hewitt C N, Erickson D, et al. A global model of natural volatile organic compounds emissions [J]. Journal of Geophysical Research, 1995,100：8873-8892.

[7]Monson R K, Holland E A. Biospheric trace gas fluxes and their control over tropospheric chemistry [J]. Annual Review of Ecology and Systematics, 2001,32：547-576.

[8]Levis S, Wiedinmyer C, Bonan G B, et al. Simulating biogenic volatile organic compound emissions in the Community Climate System Model [J]. Journal of Geophysical Research-Atmospheres, 2003,108(D21)：ACH 1-9.

[9]Lathière J, Hauglustaine D A, Noblet-DucoudréDe N, et al. Past and future changes in biogenic volatile organic compound emissions simulated with a global dynamic vegetation model [J].Geophysical Research Letters, 2005,32(20)：L20818 1-4.

[10]Guenther A B, Zimmerman P R, Harley P C, et al. Isoprene and monoterpene emission rate variability： model evaluations and sensitivity analyses [J]. Journal of Geophysical Research, 1993,98(D7)：12609-12617.

[11]Stocker T F, Qin D, Plattner G K. et al. Climate Change 2013： the Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [R]. Cambridge, UK and New York, USA：Cambridge University Press, 2013：4-23.

[12]Lindwall F, Schollert M, Michelsen A, et al. Fourfold higher tundra volatile emissions due to arctic summer warming [J].Journal of Geophysical Research： Biogeosciences, 2016,121(3)：895-902.

[13]Zhang Y, Hu X-M, Leung L R, et al. Impacts of regional climate change on biogenic emissions and air quality [J]. Journal of Geophysical Research： Atmospheres, 2008,113：D18310 1-24.

[14]Mochizuki T, Saito T, Hirai G, et al. Development of a reliable method to determine monoterpene emission rate of plants grown in an open-top chamber [J]. Journal of Agricultural Meteorology,2015,71(4)：271-275.

[15]Janson R W. Monoterpene emissions from Scots pine and Norwegian spruce.Journal of Geophysical Research-Atmospheres,1993,98(D2)：2839—2850.

[16]Noe S M, Pe?uelas J, Niinemets ü. Monoterpene emissions from ornamental trees in urban areas： a case study of Barcelona, Spain[J]. Plant Biology, 2008,10(1)：163—169.

[17]R?is?nen T, Ryypp? A, Kellom?ki S. Monoterpene emission of a boreal Scots pine (Pinussylvestris L.) forest [J]. Agricultural and Forest Meteorology, 2009,149(5)：808-819.

[18]Son Youn-Suk, Kim Ki-Joon, Jung In-Ha, et al. Seasonal variations and emission fluxes of monoterpene emitted from coniferous trees in East Asia： focused on Pinus rigida and Pinus koraiensis [J]. Journal of Atmospheric Chemistry, 2015,72：27—41.

[19]Maja M M, Kasurinen A, Holopainen T, et al. Volatile organic compounds emitted from silver birch of different provenances across a latitudinal gradient in Finland [J]. Tree Physiology,2015,35(9)：975-986.

[20]Morrison E C, Drewer J, Heal M R. A comparison of isoprene and monoterpene emission rates from the perennial bioenergy crops short-rotation coppice willow and Miscanthus and the annual arable crops wheat and oilseed rape [J]. GCB Bioenergy, 2016,8：211—225.

[21]Stewart H E, Hewitt C N, Bunce R G H, et al. A highly spatially and temporally resolved inventory for biogenic isoprene and monoterpene emissions： Model description and application to Great Britain [J]. Journal of Geophysical Research-Atmospheres,2003,108(D20)：ACH 7 1-18.

[22]Ludley K E, Jickells S M, Chamberlain P M, et al. Distribution of monoterpenes between organic resources in upper soil horizons under monocultures of Picea abies, Picea sitchensis and Pinus sylvestris [J]. Soil Biology and Biochemistry, 2009,41(6)：1050-1059.

[23]Seco R, Karl T, Turnipseed A, et al. Springtime ecosystem-scale monoterpene fluxes from Mediterranean pine forests across a precipitation gradient [J]. Agricultural and Forest Meteorology,2017,150：237-238.

[24]Asensio D, Pe?uelas J, Prieto P, et al. Interannual and seasonal changes in the soil exchange rates of monoterpenes and other VOCs in a Mediterranean shrubland [J]. European Journal of Soil Science, 2008,59(5)：878-891.

[25]Holzke C, Hoffmann T, Jaeger L, et al. Diurnal and seasonal variation of monoterpene and sesquiterpene emissions from Scots pine (Pinussylvestris L.) [J]. Atmospheric Environment, 2006,40(17)：3174-3185.

[26]Lim Jun-Ho, Kim Jo-Chun, Kim Ki-Joon, et al. Seasonal variations of monoterpene emissions from Pinus densiflora in East Asia [J]. Chemosphere, 2008,73(4)：470-478.

[27]Llusia J, Pe?uelas J, Guenther A, et al. Seasonal variations in terpene emission factors of dominant species in four ecosystems in NE Spain [J]. Atmospheric Environment, 2013,70：149-158.

[28]Filella I, Wilkinson M J, LlusiàJ, et al. Volatile organic compounds emissions in Norway spruce (Picea abies) in response to temperature changes [J]. Physiologia Plantarum, 2007,130(1)：58-66.

[29]Filella I, Pe?uelas J, Seco R. Short-chained oxygenated VOC emissions in Pinus halepensis in response to changes in water availability [J]. Acta Physiologiae Plantarum, 2009,31(2)：311-318.

[30]Laffineur Q, Aubinet M, Schoon N, et al. Impact of diffuse light on isoprene and monoterpene emissions from a mixed temperate forest [J]. Atmospheric Environment, 2013,74：385-392.

[31]Brilli F, Gioli B, Fares S, et al. Rapid leaf development drives the seasonal pattern of volatile organic compound (VOC) fluxes in a‘coppiced’ bioenergy poplar plantation [J]. Plant, Cell &Environment, 2016,39(3)：539-555.

[32]Llusia J, Roahtyn S, Yakir D, et al. Photosynthesis, stomatal conductance and terpene emission response to water availability in dry and mesic Mediterranean forests [J]. Trees, 2016,30(3)：749-759.

[33]Henke A, Kunert M, Veit D, et al. Analysis of volatiles from Picea abies triggered by below-ground interactions [J]. Environmental and Experimental Botany, 2015,110：56-61.

[34]Vuorinen T, Reddy G V P, Nerg A-M, et al. Monoterpene and herbivore-induced emissions from cabbage plants grown at elevated atmospheric CO2concentration [J]. Atmospheric Environment, 2004,38(5)：675-682.

[35]Trowbridge A M, Daly R W, Helmig D, et al. Herbivory and climate interact serially to control monoterpene emissions from pinyon pine forests [J]. Ecology, 2014,95：1591—1603.

[36]Copolovici L, Vaartnou F, Estrada M P, et al. Oak powdery mildew (Erysiphe alphitoides)-induced volatile emissions scale with the degree of infection in Quercus robur [J]. Tree Physiology,2014,34(12)：1399-1410.

[37]Constable J V H, Guenther A B, Schimel D S, et al. Modelling changes in VOC emission in response to climate change in the continental United States [J]. Global Change Biology, 1999,5(7)：791—806.

[38]Faubert P, Tiiva P, Rinnan ?, et al. Doubled volatile organic compound emissions from subarctic tundra under simulated climate warming [J]. New Phytologist, 2010,187(1)：199—208.

[39]Vedel-Petersen I, Schollert M, Nymand J, et al. Volatile organic compound emission profiles of four common arctic plants [J].Atmospheric Environment, 2015,120：117-126.

[40]Rinnan R, Rinnan ?, Faubert P, et al. Few long-term effects of simulated climate change on volatile organic compound emissions and leaf chemistry of three subarctic dwarf shrubs [J].Environmental and Experimental Botany, 2011,72：377-386.

[41]Kramsh?j M, Vedel-Petersen I, Schollert M, et al. Large increases in Arctic biogenic volatile emissions are a direct effect of warming [J]. Nature Geoscience, 2016,9(5)：349-352.

[42]Seco R, Karl T, Guenther A, et al. Ecosystem-scale volatile organic compound fluxes during an extreme drought in a broadleaf temperate forest of the Missouri Ozarks (central USA)[J]. Global Change Biology, 2015,21：3657—3674.

[43]Feller U. Drought stress and carbon assimilation in a warming climate： Reversible and irreversible impacts [J]. Journal of Plant Physiology, 2016,203：84-94.

[44]Geron C, Daly R, Harley P, et al. Large drought-induced variations in oak leaf volatile organic compound emissions during PINOT NOIR [J]. Chemosphere, 2016,146：8-21.

[45]Huang L, McDonald-Buller E C, McGaughey G, et al. Annual variability in leaf area index and isoprene and monoterpene emissions during drought years in Texas [J]. Atmospheric Environment, 2014,92：240-249.

[46]Pe?uelas J, Filella I, Seco R, et al. Increase in isoprene and monoterpene emissions after re-watering of droughted Quercus ilex seedlings [J]. Biologia Plantarum, 2009,53(2)：351-354.

[47]Copolovici L, Niinemets ü. Flooding induced emissions of volatile signalling compounds in three tree species with differing water logging tolerance [J]. Plant, Cell & Environment, 2010,33(9)：1582-1594.

[48]Constable J V H, Litvak M E, Guenther A B, et al. Monoterpene emission from coniferous trees in response to elevated CO2concentration and climate warming [J]. Global Change Biology,1999,5(3)：252—267.

[49]Loreto F, Fischbach R J, Schnitzler J-P, et al. Monoterpene emission and monoterpene synthase activities in the Mediterranean evergreen oak Quercus ilex L. grown at elevated CO2concentrations [J].Global Change Biology, 2001,7(6)：709-717.

[50]Schollert M, Kivim?enp?? M, Michelsen A, et al. Leaf anatomy,BVOC emission and CO2exchange of arctic plants following snow addition and summer warming [J]. Annals of Botany, 2017,119(3)：433-445.

[51]Baraldi R, Rapparini F, Oechel W C, et al. Monoterpene emission responses to elevated CO2in a Mediterranean-type ecosystem [J].New Phytologist, 2004,161(56)：17-21.

[52]R?is?nen T, Ryypp? A, Kellom?ki S. Effects of elevated CO2and temperature on monoterpene emission of Scots pine (Pinus sylvestris L.) [J]. Atmospheric Environment, 2008,42(18)：4160—4171.

[53]Snow M D, Bard R R, Olszyk D M, et al. Monoterpene levels in needles of Douglas fir exposed to elevated CO2and temperature[J]. Physiologia Plantarum, 2003,117(3)：352-358.

[54]Velikova V, Tsonev T, Barta C, et al. BVOC emissions,photosynthetic characteristics and changes in chloroplast ultrastructure of Platanus orientalis L. exposed to elevated CO2and high temperature [J]. Environmental Pollution, 2009,157(10)：2629-2666.

[55]Hartikainen K, Nerg A-M, Kivimaenpaa M, et al. Emissions of volatile organic compounds and leaf structural characteristics of European aspen (Populus tremula) grown under elevated ozone and temperature [J]. Tree Physiology, 2009,29(9)：1163-1173.

[56]Hartikainen K, Riikonen J, Nerg A-M, et al.Impact of elevated temperature and ozone on the emission of volatile organic compounds and gas exchange of silver birch (Betula pendula Roth) [J].Environmental and Experimental Botany, 2012,84：33-43.

[57]Kivim?enp?? M, Ghimire R P, Sutinen S, et al. Increases in volatile organic compound emissions of Scots pine in response to elevated ozone and warming are modified by herbivory and soil nitrogen availability [J]. European Journal of Forest Research,2016,135(2)：343-360.

[58]Valolahti H, Kivim?enp?? M, Faubert P, et al. Climate change-induced vegetation change as a driver of increased subarctic biogenic volatile organic compound emissions [J].Global Change Biology, 2015,21(9)：3478-3488.

[59]Kellom?ki S, Rouvinen I, Peltola H, et al. Impact of global warming on the tree species composition of boreal forests in Finland and effects on emissions of isoprenoids [J]. Global Change Biology, 2001,7(5)：531-544.

[60]Faubert P, Tiiva P, Michelsen A, et al. The shift in plant species composition in a subarctic mountain birch forest floor due to climate change would modify the biogenic volatile organic compound emission profile [J]. Plant & Soil, 2012,352(1/2)：199—215.

[61]Wang Y-F, Owen S M, Li Q-J, et al. Monoterpene emissions from rubber trees (Heveabrasiliensis) in a changing landscape and climate： chemical speciation and environmental control [J].Global Change Biology, 2007,13：2270—2282.

[62]白建輝,Tiffany D,余樹全,等.亞熱帶竹林揮發(fā)性有機物排放的模擬 [J]. 生態(tài)環(huán)境學報, 2015,24(12)：1923-1937.

[63]白建輝,林鳳友,萬曉偉,等.長白山溫帶森林揮發(fā)性有機物的排放通量 [J]. 環(huán)境科學學報, 2012,32(3)：545-554.

[64]王志輝,白郁華,王雪松,等.西雙版納地區(qū)三葉橡膠樹異戊二烯和單萜烯排放機理初步研究 [J]. 北京大學學報, 2003,39(4)：512-516.

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

[66]梁寶生,Baraldi R,Rapparini F.夏季黃連木排放萜烯類化合物濃度日變化及排放速率研究 [J]. 北京大學學報, 2003,39(4)：517-521.

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

[68]梁珍海,劉海燕,陳 霞,等.南京紫金山不同植物群落 VOCs組成 [J]. 南京林業(yè)大學學報, 2011,25(1)：4-38.

[69]常 杰,任 遠,史 琰,等.亞熱帶城鄉(xiāng)復合系統(tǒng) BVOC 排放清單-以臺州地區(qū)為例 [J]. 生態(tài)學報, 2012,32(2)：641-649.

[70]閆 雁,王志輝,白郁華,等.中國植被VOC排放清單的建立 [J].中國環(huán)境科學, 2005,25(1)：110-114.

[71]池彥琪,謝紹東.基于蓄積量和產(chǎn)量的中國天然源 VOC排放清單及時空分布 [J]. 北京大學學報, 2012,48(3)：475-482.

[72]李金龍,白郁華,胡建信,等.油松排放萜烯類化合物濃度的日變化及排放速率的研究 [J]. 中國環(huán)境科學, 1994,14(3)：165-169.

[73]侯慧波,李新崗,馬養(yǎng)民,等.2種吸附劑對油松揮發(fā)物的吸附效果比較 [J]. 西北林學院學報, 2006,21(2)：134-137.

[74]李繼泉,金幼菊,沈應柏,等.環(huán)境因子對植物釋放揮發(fā)性化合物的影響 [J]. 植物學通報, 2001,18(6)：649-656.

[75]呂 迪,王得祥,謝小洋,等.油松釋放萜烯類揮發(fā)性成分研究[J]. 西北林學院學報, 2016,31(1)：231-237.

[76]周秀驥,唐孝炎,馮宗煒,等.長江三角洲低層大氣與生態(tài)系統(tǒng)相互作用研究 [M]. 北京：氣象出版社, 2004：96-102.

[77]高文華,李忠勤,張明軍,等.西晉南地區(qū)近 50a氣候變化特征,突變與周期分析 [J]. 干旱區(qū)資源與環(huán)境, 2011,25(7)：124-127.

[78]趙桂香,趙彩萍,李新生,等.47a來山西省氣候變化分析 [J]. 干旱區(qū)研究, 2006,23(3)：499-505.

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

[80]Moss R H, Edmonds J A, Hibbard K A, et al. The next generation of scenarios for climate change research and assessment [J].Nature, 2010,463(7282)：746-756.

[81]Taylor K E, Stouffer R J, Meehl G A. An overview of CMIP5 and experiment design [J]. Bulletin of the American Meterological Society, 2012,93(4)：485-498.

[82]Meinshausen M, Smith S J, Calvin K, et al. The RCP greenhouse gas concentrations and their extensions from 1765 to 2300 [J].Climatic Change, 2011,109(16-2)：213-241.

[83]Warszawaski L, Frieler K, Huber V, et al. the inter-sectoral impact model intercomparison project (isi-mip)：project framework [J]. Proceedings of the National Academy Of Sciences,2013,119(9)：3228-3232.

[84]Hempel S, Frieler K, Warszawski L, et al. A trend-preserving bias correction-the isi-mip approach [J]. Earth Syst Dynamo,2013,42920：219-236.

[85]吳建國.氣候變化對生物多樣性影響與風險研究 [M]. 北京：科學出版社, 2017：73-116.

[86]Ephrath J E, Goudriaan J, Marani A. Modelling diurnal patterns of air temperature, radiation wind speed and relative humidity by equations from daily characteristics [J]. Agricultural Systems,1996,51(4)：377-393.

[87]盛裴軒,毛節(jié)泰,李建國,等.大氣物理學 [M]. 北京：北京大學出版社, 2003：88-89.

[88]Schollert M, Kivim?enp?? M,Valolahti H M, et al. Climate change alters leaf anatomy, but has no effects on volatile emissions from arctic plants [J]. Plant, Cell and Environment,2015,38：2048—2060.

[89]司徒淑娉,王雪梅,Guenther A,等.氣象模擬誤差對單萜烯排放估算的影響 [J]. 環(huán)境科學學報, 2010,30(12)：2383-2391.

[90]毛紅梅,張凱山,第寶鋒.四川省天然源VOCs排放量的估算和時空分布 [J]. 中國環(huán)境科學, 2016,36(5)：1289-1296.

[91]張 蕾,姬亞芹,趙 杰,等.烏魯木齊市天然源 VOCs排放量估算與時空分布特征 [J]. 中國環(huán)境科學, 2017,37(10)：3692-3698.