阿爾泰山薩彥嶺4種優(yōu)勢樹種徑向生長對氣候因子的響應(yīng)

康 劍,蔣少偉,黃建國,*

1 中國科學(xué)院華南植物園,退化生態(tài)系統(tǒng)植被恢復(fù)與管理重點實驗室,廣州 510650 2 中國科學(xué)院大學(xué)資源與環(huán)境學(xué)院,北京 100049 3 中國科學(xué)院核心植物園,廣州 510650

近年來,全球變化趨勢明顯,人類活動導(dǎo)致的全球變暖比工業(yè)化前的水平高出0.8—1.2℃之間。如果全球變暖繼續(xù)以目前的速度增長,預(yù)計在2030年至2052年間,全球溫度變暖增幅可能會達(dá)到1.5℃[1]。森林作為重要的碳匯,極易受到全球變化的影響。樹木作為森林的組成元素,其生長除受其自身生理因素影響外,還與氣候因子緊密相關(guān)。因此,氣候變化對全球范圍內(nèi)的森林生態(tài)系統(tǒng)產(chǎn)生重要影響[2],特別是作為全球最大生物群落之一的北方森林。研究氣候變化背景下北方森林優(yōu)勢樹種間徑向生長-氣候變化關(guān)系的差異,有助于反映森林生態(tài)系統(tǒng)動態(tài)發(fā)展,為森林生態(tài)管理與可持續(xù)發(fā)展提供科學(xué)依據(jù)。由于具有年分辨率高以及可靠性強的特點,樹木年輪學(xué)成為研究氣候變化條件下樹木生長-氣候關(guān)系的重要手段[3]。

研究發(fā)現(xiàn)除了海拔、緯度等因子是影響樹木生長-氣候因子間關(guān)系的重要因素外,同一生境中不同樹種的差異也會導(dǎo)致響應(yīng)的分異[4- 9]。黃建國等通過對加拿大東部的北方森林中的4種主要樹種研究發(fā)現(xiàn),黑云杉(Piceamariana)和北美短葉松(Pinusbanksiana)受到生長季溫度的影響,而生長季前的降水對闊葉樹種山楊(Populustremuloides)、北美白樺(Betulapapyrifera)的生長有積極作用[10]。張衛(wèi)國等對玉龍雪山3個針葉樹種對氣候響應(yīng)的研究發(fā)現(xiàn),大果紅杉(Larixpotaninii)和長苞冷杉(Abiesgeorgei)主要受到生長季溫度、水分的限制,而麗江云杉(Picealikiangensis)則主要受限于生長季前的水熱條件[11]。王婷等通過對比油松和華山松年表與氣候變化的相關(guān)關(guān)系,發(fā)現(xiàn)2個樹種與溫度、降水的響應(yīng)發(fā)生在不同時期[12]。上述研究結(jié)果表明,徑向生長-氣候變化關(guān)系在優(yōu)勢樹種種間存在差異。

阿爾泰山位于北半球的中高緯度地區(qū),是中亞地區(qū)最大的山系之一,也是中亞氣候變化的敏感區(qū)域,在全球氣候變化研究中具有特殊地位[13-14]。其森林生態(tài)系統(tǒng)被普遍認(rèn)為是對氣候變化響應(yīng)最敏感的生態(tài)系統(tǒng)[15]。阿爾泰山西薩彥嶺地區(qū)森林覆蓋廣泛,為典型的北方森林,包括西伯利亞落葉松(L.sibirica)、西伯利亞紅松(P.sibirica)、西伯利亞冷杉(A.sibirica)以及西伯利亞云杉(P.obovata)等優(yōu)勢種,且研究區(qū)域人為干擾較少,是研究氣候變化-樹木生長關(guān)系的理想?yún)^(qū)域。作為全球重要的生物多樣性保護(hù)地區(qū)之一,該區(qū)域目前已經(jīng)開展了部分對樹木徑向生長-氣候變化關(guān)系的研究,且發(fā)現(xiàn)樹木徑向生長主要受到溫度、水分的影響[16-17]。但目前還沒有針對優(yōu)勢樹種間響應(yīng)的報道,而徑向生長-氣候關(guān)系在樹種間的差異往往影響了氣候變化下森林的響應(yīng)。因此,研究優(yōu)勢樹種間生長對氣候變化的敏感度差異,對正確認(rèn)識氣候變化下林分生長動態(tài)及分布格局非常重要。本文通過測定阿爾泰山西薩彥嶺的4種優(yōu)勢樹種的樹輪寬度數(shù)據(jù),建立標(biāo)準(zhǔn)年表,對76年間(1940—2016)4種優(yōu)勢樹種與6個氣候因子(平均溫度(T)、最高溫度(TMAX)、最低溫度(TMIN)、降水(PCP)、干旱指數(shù)(PDSI)、水汽壓(VAP))的響應(yīng)做了深入分析,探究不同樹種徑向生長對氣候的響應(yīng)異同。該結(jié)果有助于完善該區(qū)域優(yōu)勢樹種徑向生長-氣候變化關(guān)系,預(yù)測未來林分動態(tài)以及分布格局,旨在為我國“一帶一路”戰(zhàn)略的重要生態(tài)屏障保護(hù)提供科技支撐。

1 研究區(qū)域

研究區(qū)域為俄羅斯境內(nèi)的西薩彥嶺,屬于阿爾泰山東支,位于蒙古高原的北沿、西伯利亞的南緣(42°N—53°N,93°E—94°E)(圖1)。該地區(qū)屬于寒冷帶常濕冷溫氣候,受到極地氣旋和西伯利亞寒流的影響,夏季溫暖相對較短,冬季嚴(yán)寒漫長[18-19]。該區(qū)域年均氣溫-5.38 ℃,其中7月份月均溫度最高,為13.92 ℃；1月份月均溫度最低,為-27.17 ℃,1—4月、10—12月平均氣溫均為零下。年均總降水量925.55 mm,且集中在5—8月,占全年總降水量的50%(圖2)。研究區(qū)域的植被具有高度的垂直地帶性,受地形影響,北部被大面積北方森林所覆蓋。海拔2200—2500 m達(dá)到樹線,樹線以上則主要是灌木與亞高山草甸等植被類型[20]。

圖1 研究區(qū)位置示意及樣點位置Fig.1 Study area location and sites locationRUABSI：西伯利亞冷杉樣點,Russia Abies sibirica；RUPIOB：西伯利亞云杉樣點,Russia Picea obovata；RUPISI：西伯利亞紅松樣點,Russia Pinus sibirica；RULASI：西伯利亞落葉松樣點,Russia Larix sibirica

圖2 研究區(qū)1940—2016年區(qū)間月平均氣溫、月平均最低氣溫、月平均最高氣溫和降水的變化趨勢Fig.2 The variation trend of monthly mean temperature,monthly mean minimum temperature,monthly mean maximum temperature and precipitation in the research area from 1940 to 2016

2 研究方法

2.1 樣本采集

西伯利亞冷杉喜濕耐寒,多生長于陰坡半陰坡區(qū)域,常與西伯利亞落葉松或西伯利亞云杉等樹種組成混交林[21]。西伯利亞云杉耐寒耐陰,多為純林,分布于海拔1200—2400 m范圍內(nèi),是北方森林的重要建群種之一[22]。西伯利亞紅松的生境與冷杉較為相似,同樣不耐干旱[23]。而西伯利亞落葉松為落葉針葉樹種,該樹種喜光,相對較為耐旱[24](圖3)。研究區(qū)域內(nèi),分布有針葉樹混交林,針闊葉混交林,此外還有部分以純林形式存在。2017年7月于研究區(qū)海拔1200—1300 m之間的區(qū)域內(nèi)選取了4個樣地,樣地分別為各樹種與疣枝樺(Betulapendula)組成的混交林,每個樣地分別選擇一個優(yōu)勢樹種,并在每個樣地中均設(shè)置 1 個 20 m×20 m 的無人為干擾、未經(jīng)過火災(zāi)、蟲害等的林分作為樣方,樣方內(nèi)有可采樣的樹木15株以上(表1),樣方內(nèi)平均樹木高度22.5 m,平均胸徑35.7 cm,立木密度約為每公頃1875株。用生長錐在每株樹木胸徑處(樹高1.3 m處)從平行于所在坡面的兩個不同方向分別鉆取一個樣芯,封裝于塑料管中帶回實驗室。

圖3 主要樹種Fig.3 Four dominant tree species

2.2 年表建立

采集的樹輪樣本在實驗室經(jīng)過干燥、固定、打磨等處理后,運用樹輪測量專業(yè)儀器 LinTab樹輪寬度測量儀及配套TSAP軟件(LinTab 6.0,Frank Rinntech,Heidelberg,Germany)測量年輪寬度(精度為0.001 mm),采用目視初步定年、骨架圖定年等方式交叉定年,確定每一輪的精確年份,并使用 COFFECHA 程序檢測交叉定年的準(zhǔn)確度以及缺失輪判斷的可靠性[25]。運用R軟件中的“dplR”包[26],使用響應(yīng)截止頻率為50%,級數(shù)長度67%的樣條函數(shù)對年表進(jìn)行去趨勢處理[27],以消除與樹齡有關(guān)的生長趨勢及非氣候因子所引起的樹木生長波動,最終建立各樹種的標(biāo)準(zhǔn)年表(STD,Standard chronology)(圖4)。同時,為了更好的評估年表的質(zhì)量,計算了部分參數(shù),包括平均敏感度(MS,Mean sensitivity)、一階自回歸系數(shù)(AR1,First-order serial autocorrelation)、樣本相關(guān)系數(shù)(rbar,Correlation coefficient for all series)、第一主分量方差(PC1,Variance in first eigenvector)、信噪比(SNR,Signal-to-noise ratio)以及樣本總體解釋量(EPS,Expressed population signal)。

表1 采樣基本信息表

2.3 氣象資料

所有樣點均位于CRU規(guī)范的同一個0.5° × 0.5°的格點當(dāng)中,格點歷史氣象插值資料(主要為月資料)從KNMI Climate Explorer網(wǎng)站CRU TS 4.02下載(https://climexp.knmi.nl/start.cgi?id=someone@somewhere),并整理用于后續(xù)統(tǒng)計分析。數(shù)據(jù)長度為1940—2016年,主要指標(biāo)包括月平均溫度(T)、月平均最低溫度(TMIN)、月平均最高溫度(TMAX)、月降水量(PCP)、水汽壓(VAP)、帕爾默干旱指數(shù)(PDSI)[28](PDSI是土壤干旱程度的重要表征指數(shù))。除當(dāng)年數(shù)據(jù)外,還選擇了上一年8—12月的上述氣候指標(biāo)作為分析對象。

2.4 數(shù)據(jù)分析

利用 R 軟件,計算各年表的參數(shù),運用R包“treeclim”對不同樹種與氣候因子之間的Pearson系數(shù)進(jìn)行分析,該函數(shù)根據(jù)年輪年表與每月的氣候數(shù)據(jù)計算響應(yīng)和相關(guān)系數(shù),同時在Bootstrap重采樣分析中獲得顯著性檢驗結(jié)果,所有圖件使用R、ArcMap等軟件制作[29-30]。

3 結(jié)果與分析

3.1 4種不同樹種年表的特征分析

通過對研究區(qū)域4種優(yōu)勢樹種樹輪年表統(tǒng)計參數(shù)特征分析發(fā)現(xiàn)(表2),年表序列長度的平均值為215年。平均敏感度(MS)較高,為0.11—0.14,其中西伯利亞落葉松(LY)的MS最高為0.14。一階自回歸系數(shù)(AR1)均較高,在0.41—0.71間。第一主分量方差(PC1)均在34.8%以上,樣本間相關(guān)系數(shù)均大于0.30,信噪比(SNR)的數(shù)值除了西伯利亞落葉松外均保持在18.77—19.26,表明年表中用于分析的環(huán)境信息量較大。樣本總體解釋量(EPS)的數(shù)值均高于0.85的臨界值[31]。因此,本研究建立的年表具有一定的敏感性,且包含了由氣候產(chǎn)生的低頻變化和樹木生理滯后效應(yīng)[26],氣候因子對樹木徑向生長影響較大,年表中含有較多的氣候信息[32],采集的樣本量中的信號可以代表總體特征。

對4種樹種的標(biāo)準(zhǔn)年表(1940—2016)間的相關(guān)分析發(fā)現(xiàn),西伯利亞云杉與西伯利亞冷杉年表的相關(guān)性最高(r=0.68,P<0.05),其次為西伯利亞云杉的年表與西伯利亞紅松之間的相關(guān)性(r=0.65,P<0.05),此外,西伯利亞云杉與西伯利亞落葉松的相關(guān)性也較高(r=0.47,P<0.05),相關(guān)系數(shù)反映了采樣點樹種之間的生長一致性(圖5)。

表2 4種優(yōu)勢樹種標(biāo)準(zhǔn)年表主要特征參數(shù)

圖4 4種優(yōu)勢樹種的標(biāo)準(zhǔn)年表和樣本量Fig.4 Standard chronologies and sample size of four dominant species

圖5 4種優(yōu)勢樹種的標(biāo)準(zhǔn)年表的相關(guān)性(1940—2016)Fig.5 Correlation of standard chronologies of four dominant species in common period from 1940 to 2016RUABSI：西伯利亞冷杉樣點,Russia Abies sibirica；RUPIOB：西伯利亞云杉樣點,Russia Picea obovata；RUPISI：西伯利亞紅松樣點,Russia Pinus sibirica；RULASI：西伯利亞落葉松樣點,Russia Larix sibirica

3.2 4種優(yōu)勢樹種對氣候因子的響應(yīng)

通過對4種樹種年表分別與6個氣候因子做相關(guān)分析發(fā)現(xiàn),4種優(yōu)勢樹種徑向生長的限制因子不盡相同。其中,西伯利亞冷杉的徑向生長與2—4月(生長季前期)(r=0.20,r=0.23,r=0.33；P<0.05)的降水量(PCP)顯著正相關(guān),與干旱指數(shù)(PDSI)有著很強的相關(guān)性,表現(xiàn)為上一年的10、11月(r=0.25,r=0.30；P<0.05)、當(dāng)年1—9月(r=0.25—0.43,P<0.05)該樹種的徑向生長與干旱指數(shù)呈顯著正相關(guān)。當(dāng)年1月份的平均溫度(T)、最高溫度(TMAX)與徑向生長有顯著負(fù)相關(guān)關(guān)系(r=-0.16,r=-0.16；P<0.05)。此外,當(dāng)年4、6月的水汽壓(VAP)與徑向生長顯著正相關(guān)(r=0.27,r=0.26；P<0.05)(圖6)。

圖6 西伯利亞冷杉樹輪寬度年表(1940—2016年)與月氣候因子的相關(guān)性Fig.6 Correlation between ring width chronology of Abies sibirica (1940—2016) and monthly climate factorsp8—p12表示上一年8—12月份,* 表示達(dá)到 0.05 水平的顯著相關(guān);PCP：降水,Precipitation；PDSI：干旱指數(shù),Palmer Drought Severity Index；T：溫度,Temperature；TMAX,最高溫度,Maximum temperature；TMIN：最低溫度,Minimum temperature；VAP：水汽壓,Vapour pressure

西伯利亞落葉松的徑向生長與降水量的相關(guān)性未達(dá)到顯著性水平,但與上一年8、9月的干旱指數(shù)顯著負(fù)相關(guān)(r=-0.25,r=-0.15；P<0.05)。此外,落葉松生長對溫度的響應(yīng)較為明顯,其中上一年8月和當(dāng)年8月的平均溫(r=-0.11,r=-0.19；P<0.05)與徑向生長顯著負(fù)相關(guān),當(dāng)年8月的平均最高溫也表現(xiàn)相同的響應(yīng)(r=-0.19,P<0.05),上一年8月和當(dāng)年8月的平均最低溫(TMIN)與徑向生長顯著負(fù)相關(guān)(r=-0.14,r=-0.16；P<0.05),而在當(dāng)年6月則呈現(xiàn)顯著正相關(guān)(r=0.23,P<0.05)。水汽壓在當(dāng)年8月的值也顯示與落葉松徑向生長顯著負(fù)相關(guān)(r=-0.13,P<0.05)(圖7)。

圖7 西伯利亞落葉松樹輪寬度年表(1940—2016年)與月氣候因子的相關(guān)性Fig.7 Correlation between ring width chronology of Larix sibirica (1940—2016) and monthly climate factors

4種樹種當(dāng)中,西伯利亞紅松徑向生長與氣候因子的相關(guān)性較低,僅與當(dāng)年3月降水、當(dāng)年7月平均最低溫、上一年10月的水汽壓存在顯著正相關(guān)的關(guān)系(r=0.24,r=0.25,r=0.20；P<0.05)(圖8)。

圖8 西伯利亞紅松樹輪寬度年表(1940—2016年)與月氣候因子的相關(guān)性Fig.8 Correlation between ring width chronology of Pinus sibirica (1940—2016) and monthly climate factors

西伯利亞云杉的徑向生長與當(dāng)年3、4月的降水量顯著正相關(guān)(r=0.23,r=0.21；P<0.05),與上一年10—11月、當(dāng)年2—4、9月的干旱指數(shù)有顯著正相關(guān)關(guān)系(r=0.24—0.31；P<0.05)。在溫度響應(yīng)方面,與6月的平均溫、平均最高溫都有顯著正相關(guān)關(guān)系(r=0.32,r=0.33；P<0.05)。此外,該樹種徑向生長還與上一年10月、當(dāng)年6月的水汽壓顯著正相關(guān)(r=0.25,r=0.39；P<0.05)(圖9)。

圖9 西伯利亞云杉樹輪寬度年表(1940—2016年)與月氣候因子的相關(guān)性Fig.9 Correlation between ring width chronology of Picea obovata (1940—2016) and monthly climate factors

4 討論

研究結(jié)果充分表明,對于同一氣候條件下,西伯利亞冷杉與西伯利亞云杉均與干旱指數(shù)有正響應(yīng)關(guān)系,6月份溫度對于西伯利亞云杉和西伯利亞落葉松有著相似的正響應(yīng),西伯利亞云杉與西伯利亞紅松均受到生長季前的降水影響。這可能是由于在大環(huán)境下,樹種的徑向生長在特定方面產(chǎn)生的一致性信號,樹種年表間的相關(guān)性分析結(jié)果與上述的樹種間徑向生長-氣候關(guān)系的共性是一致的。但是,不同的樹種徑向生長-氣候因子的關(guān)系也存在差異,包括對同一氣候因子的響應(yīng)方式、響應(yīng)時間、響應(yīng)程度等方面都存在一定的差異,這種差異可能來自于不同樹種的生理特征以及遺傳差異。

西伯利亞冷杉的徑向生長與2—4月降水量顯著正相關(guān),與上一年10、11月和當(dāng)年1—9月干旱指數(shù)顯著正相關(guān),而與1月溫度顯著負(fù)相關(guān),該結(jié)果說明該區(qū)域西伯利亞冷杉徑向生長主要受到干旱限制,為干旱敏感樹種。生長季前期的降雪增加,有利于對根部的保溫[12],同時,這一期間的降雪融化是生長季前期細(xì)胞水分的主要來源,充足的水分可以促進(jìn)徑向生長[33]。植物生長需要不斷通過根部從土壤吸收充足的水分來維持生理活動,冷杉屬于喜濕樹種,其生長極易受到干旱的限制[34]。生長季前期,樹木正處于休眠期,過高的溫度會導(dǎo)致樹木消耗增加,不利于樹木生長期前期的物質(zhì)積累,影響到即將到來的生長旺盛期。研究區(qū)降水在4月份開始逐漸增加,而此時植物還處于水分脅迫的狀態(tài),此時水汽壓力的增大并不會抑制葉片的氣孔活動,反而有助于緩解水分脅迫帶來的壓力。

而西伯利亞落葉松則表現(xiàn)不同,其生長主要受溫度的影響,主要與8月的溫度、水氣壓顯著負(fù)相關(guān),而與6月平均最低溫正相關(guān)。該區(qū)域落葉松的生長季在6—8月之間,10月落葉,6月為樹木萌發(fā)新葉、快速發(fā)育的初期,由于此時融雪、降雨都在增加,水分充足,溫度的升高有利于植物進(jìn)行光合作用,從而提供更多的營養(yǎng)物質(zhì),促進(jìn)樹木生長,該結(jié)果與此前在阿爾泰山關(guān)于落葉松的研究一致[35-37]。生長季末期,由于降水量依舊較大,8月溫度的升高會導(dǎo)致空氣濕度增加,使得葉片氣孔關(guān)閉,光合作用減弱,進(jìn)一步影響樹木生長[38]。

樹種生物學(xué)特征的不同,對氣候變化的響應(yīng)特征以及敏感性強弱也有一定的差異[39]。西伯利亞紅松與氣候因子的響應(yīng)較弱,其生長受到3月降水、7月平均最低溫的顯著正影響。在樹木生長季前期,氣溫整體較低,而較多的降雪可以有效的保護(hù)土壤溫度,從而降低低溫對根系的損傷。而生長季旺盛時期,溫度的上升能有效促進(jìn)酶活,提高植物光合作用速率,從而有利于樹木徑向生長。而有關(guān)于阿爾泰山中部的西伯利亞紅松與氣候變化的響應(yīng)研究表明,6月份溫度對樹木徑向生長起直接作用[40-42]。本研究結(jié)果與該研究結(jié)果不完全一致,可能是由于局部區(qū)域的氣候差異所導(dǎo)致的。

西伯利亞云杉生長與當(dāng)年3—4月的降水量,上一年10—11月、當(dāng)年2—4、9月的干旱指數(shù)正相關(guān),與6月的平均溫度、平均最高溫度、水氣壓顯著正相關(guān)。該結(jié)果與新疆阿爾泰山的西伯利亞云杉生長與生長季溫度負(fù)相關(guān)的研究結(jié)果不同,可能是由于研究區(qū)的生長季溫度較新疆阿爾泰山更低,因此呈現(xiàn)出不同的響應(yīng)結(jié)果[43]。生長季前期,充足的降雪可以有效地維持土溫。同時,在融雪初期,大量的水分也有利于樹木生長。伴隨生長季到來,氣溫逐漸升高,土壤水分蒸發(fā)量增加,水汽壓增大,植物葉片可以從空氣中吸收水分,提升了光合作用的速率,有利于有機碳物質(zhì)的積累,從而加快生長。由于該樹種喜濕,生長季期間(6—8月),土壤水分充足,可以保證樹木的快速生長,而在生長季前、后兩個階段由于土壤水分的不足,樹木生長會受到干旱帶來的限制。

5 結(jié)論

本研究表明,樹木徑向生長-氣候因子的關(guān)系在樹種間存在一定的相似性,西伯利亞冷杉與西伯利亞云杉、西伯利亞云杉和西伯利亞落葉松、西伯利亞云杉與西伯利亞紅松對于特定氣候因子的響應(yīng)有一致性。但同時,也存在著顯著的差異,這種差異可能主要緣于樹種生理特征的差異以及種間遺傳因素的差異。西伯利亞冷杉的生長對干旱的響應(yīng)最為強烈,因此主要受干旱限制。而溫度則是西伯利亞落葉松的主要限制因子,其徑向生長主要受到生長季低溫及生長后期高溫、高濕的不利影響。西伯利亞紅松對于氣候因子的敏感性較弱,主要受生長季前期降水、生長季溫度的限制。西伯利亞云杉在生長季受溫度限制,生長季之外主要受水分限制。綜上所述,西伯利亞冷杉和西伯利亞云杉主要對區(qū)域水分變化敏感,而西伯利亞落葉松和西伯利亞紅松則對溫度更敏感。本研究結(jié)果說明,在氣候變化的背景下,優(yōu)勢樹種種間對氣候響應(yīng)的差異會對林分結(jié)構(gòu)、森林動態(tài)等產(chǎn)生重要影響,進(jìn)一步證明了多樹種徑向生長-氣候關(guān)系的研究有助于正確反映森林動態(tài)。該區(qū)域優(yōu)勢樹種的可能會隨著氣候變化表現(xiàn)出不同的生長狀態(tài),從而影響該地區(qū)森林生態(tài)系統(tǒng)的群落組成以及分布格局,該結(jié)果可以為該區(qū)域森林管理與生態(tài)保護(hù)工作提供依據(jù)。