邵殿坤，欒志慧，吳慧書，王 剛，王延紅

(1.吉林省林業(yè)調(diào)查規(guī)劃院，吉林 長(zhǎng)春 130022;2.通化師范學(xué)院，吉林 通化 134002)

長(zhǎng)白落葉松(Larix olgensis)是一種具生長(zhǎng)速度快、耐寒、耐旱、耐水濕特點(diǎn)的優(yōu)良大徑級(jí)用材樹種，為長(zhǎng)白山林區(qū)主要用材樹種之一。隨著碳循環(huán)研究深入開展，陸地森林生態(tài)系統(tǒng)碳儲(chǔ)量受到研究者普遍關(guān)注。森林碳庫(kù)可分為地上和地下兩部分［1～3］。地上生物量作為生態(tài)系統(tǒng)生產(chǎn)力重要組成部分，是整個(gè)生態(tài)系統(tǒng)運(yùn)行的能量基礎(chǔ)和營(yíng)養(yǎng)物質(zhì)來源［4～5］。準(zhǔn)確測(cè)定森林生物量一直是困擾生態(tài)學(xué)家和林學(xué)家的難點(diǎn)問題［6］。國(guó)內(nèi)學(xué)者對(duì)長(zhǎng)白山林區(qū)森林生物量進(jìn)行大量研究，如李文華等［7］分析長(zhǎng)白山主要生態(tài)系統(tǒng)生物量，郭忠玲等［8］研究長(zhǎng)白山各植被帶主要樹種凋落物分解速率，并建立模型，代力民等［9］分析紅松林倒木貯量的變化規(guī)律，楊麗韞等［10］研究暗針葉林中倒木貯量和分解速率。作者對(duì)吉林省東部山地林區(qū)立地條件下長(zhǎng)白落葉松生物量異速生長(zhǎng)模型進(jìn)行研究，并對(duì)長(zhǎng)白落葉松林碳匯儲(chǔ)量進(jìn)行估測(cè)。

1 研究區(qū)域自然概況

長(zhǎng)白落葉松自然分布在我國(guó)長(zhǎng)白山、張廣才嶺，地理位置:41°30'～45°20'N，126°30'～131°40'E。本研究選擇長(zhǎng)白山區(qū)延吉、敦化、和龍、龍井、蛟河五縣市12 個(gè)林場(chǎng)長(zhǎng)白落葉松林分為研究對(duì)象(42°01'～44°04'N;127°56'～131°08'E)，海拔405 ～1 131 m，山巒起伏、河流密布。該區(qū)氣候?qū)偌撅L(fēng)影響溫帶大陸性山地氣候，冬季漫長(zhǎng)寒冷，春季風(fēng)大干燥，夏季短暫、溫暖多雨。全區(qū)年平均氣溫2.7℃～4.9℃之間。年降水量一般為525.9 ～607.1 mm。研究區(qū)主要樹種為長(zhǎng)白落葉松，林分為經(jīng)過擇伐的天然純林，主要灌木和草本植物有蕨類植物(Selaginella spp.)、莎草(Cyperus tenuispica)、忍冬(Lonicera japonica)、衛(wèi)矛(Euonymus alatus)、榛子(Ostryopsis davidiana)等。

2 研究方法

2.1 數(shù)據(jù)收集

利用全站儀在試驗(yàn)區(qū)域內(nèi)選擇典型地段設(shè)定標(biāo)準(zhǔn)地，標(biāo)準(zhǔn)地大小為60 m×40 m，面積0.24 hm2;按10 m 間距分割成10 m×10 m 樣方，調(diào)查標(biāo)準(zhǔn)地內(nèi)所有林木特征值。利用x、y直角坐標(biāo)系，在標(biāo)準(zhǔn)地內(nèi)所編寫行列號(hào)即該調(diào)查單元單元號(hào)。以每個(gè)調(diào)查單元西南角為坐標(biāo)原點(diǎn)，用皮尺測(cè)量每株樹木在該調(diào)查單元內(nèi)坐標(biāo)(x，y)。測(cè)量每株樹木地徑(D0)、胸徑(D)、樹高(H)、枝下高(h)、冠幅(cr)等因子。在標(biāo)準(zhǔn)地內(nèi)按平均胸徑和平均高選擇1 株落葉松作為平均木，按徑階選取解析木。將解析木所有樹枝齊樹干鋸下，分別測(cè)定其樹干、樹冠、樹枝、樹葉生物量，并取樣做進(jìn)一步分析。選取不同徑級(jí)大小長(zhǎng)白落葉松樣木共計(jì)30 株。將各器官樣品帶回試驗(yàn)室，置于烘干箱，在105℃下烘干，稱取樣品干物質(zhì)質(zhì)量，通過計(jì)算各器官的含水率，得到相應(yīng)干重。長(zhǎng)白落葉松樣木生物量測(cè)定結(jié)果見表1。

表1 長(zhǎng)白落葉松30 株樣木各組分生物量Tab.1 Each component biomass of 30 Larix olgensis sample trees

?

2.2 數(shù)據(jù)分析方法

首先利用樹干直徑建立生物量異速生長(zhǎng)方程。方程形式為:

式中:a 和b 是方程的系數(shù);

Y 為各器官生物量;

X 為 樹 干 直 徑 (DBH1.3和DBH0)［11～13］。

利用樹干直徑和樹高建立各器官生物量?jī)勺兞炕貧w模型。方程形式為

式中:X1為樹干直徑;

X2為樹高;

a、b、c 為方程系數(shù)。

方程擬合優(yōu)度利用決定系數(shù)(R2)和信息指數(shù)(AIC)進(jìn)行評(píng)價(jià)。

本文利用森林資源二類調(diào)查數(shù)據(jù)，計(jì)算長(zhǎng)白落葉松林碳儲(chǔ)量。利用長(zhǎng)白落葉松林生物量與蓄積量關(guān)系估算生物量，估算長(zhǎng)白落葉松林碳儲(chǔ)量與碳密度。碳儲(chǔ)量計(jì)算公式為:

式中:C 為碳儲(chǔ)量(t);

B 為地上生物量;

Cc為含碳率。

碳密度計(jì)算公式為

式中:Cp為碳密度(t·hm－2);

C 為碳儲(chǔ)量(t);

S 為森林面積(hm2)。

喬 木 樹 種 含 碳 率 通 常 為 0. 45 或0.5［14～16］。馬欽彥研究認(rèn)為，喬木樹種平均含碳率均大于0.45，其中闊葉樹種平均含碳率值多數(shù)低于0.5，針葉樹平均含碳率多數(shù)等于或高于0.5，因此使用0.5 作為轉(zhuǎn)換系數(shù)要優(yōu)于0.45［17］。本研究含碳率采用0.5 進(jìn)行計(jì)算。

圖1 長(zhǎng)白落葉松地上生物量分配Fig.1 The above－ground biomass distribution of Larix olgensis

3 結(jié)果與分析

3.1 長(zhǎng)白落葉松地上生物量異速生長(zhǎng)模型

長(zhǎng)白落葉松地上生物量分配如圖1 所示。

從圖1 中可以看出，地上生物量以樹干生物量為主，生物量各組分所占比重隨著年齡和胸徑增大而增大。長(zhǎng)白落葉松總平均生物量為215.021 kg，各組分生物量大小差異明顯，樹干生物量占地上總生物量67.15%，活枝、針葉、樹皮生物量分別占17.09%、6.00%和9.76%。

分別采用單變量(DBH)和雙變量(DBH、H)模型建立長(zhǎng)白落葉松不同器官生物量最優(yōu)模型，見表2。

表2 長(zhǎng)白落葉松各器官生物量異速生長(zhǎng)模型Tab.2 Allometric models for each organ biomass of Larix olgensis

從表2 中可以看出，所有模型F 檢驗(yàn)均顯著(p ＜0.001)，基于胸徑的單變量模型可解釋量均超過90%，其中樹干生物量、樹皮生物量和地上總生物量回歸模型可解釋量達(dá)到99%。利用胸徑和樹高建立的樹干生物量模型，決定系數(shù)Radj2 最大，信息指數(shù)AIC 最小?；钪歪樔~生物量回歸模型可解釋量超過93%，與其他3 個(gè)回歸模型相比信息指數(shù)AIC 值最小。因此，利用胸徑和樹高可獲得樹干生物量最優(yōu)模型，模型方程為Y=0.059 DBH1.32.171H0.420。胸徑對(duì)樹皮、活枝、針葉和地上總生物量預(yù)測(cè)較重要，模 型 方 程 分 別 為:樹 皮，地上總生物量Y=

3.2 生物量與蓄積量線性關(guān)系

利用線性回歸方程建立樹干、樹皮、樹枝、樹葉生物量與樹干蓄積量相關(guān)方程，如圖2 所示。

圖2 各組分生物量與蓄積量之間關(guān)系Fig.2 Relationship between biomass and stand volume of each component

從圖2 中可以看出，擬合方程分別為:樹干y=448.68 x+4.433 1，樹枝y=105.21 x+3.944 9，針葉y=32.89 x+2.639 8，樹皮y=57.39 x+3.099 7。各組分生物量與蓄積量均顯著線性相關(guān)(p ＜0.001)，R2值達(dá)到0.9 以上。方精云［18］和孫玉軍［19］研究顯示林木生物量與蓄積量之間存在線性關(guān)系，與本文結(jié)果一致。生物量與蓄積量之間的線性關(guān)系可以作為估測(cè)森林地上生物量的基礎(chǔ)。

3.3 長(zhǎng)白落葉松林碳儲(chǔ)量與碳密度

利用長(zhǎng)白落葉松生物量與蓄積量關(guān)系，估算不同林齡長(zhǎng)白落葉松林的碳儲(chǔ)量與碳密度，見表3。

表3 長(zhǎng)白落葉松林碳儲(chǔ)量和碳密度Tab.3 Carbon storage and carbon density of Larix olgensis forest

從表3 中可以看出，長(zhǎng)白山地區(qū)長(zhǎng)白落葉松幼齡林面積為3.06×105hm2，蓄積量為1.66×106m3;中齡林面積為4.71×105hm2，蓄積量為2.72×106m3;近成熟林面積為5.42×105hm2，蓄積量為2.34×106m3。長(zhǎng)白落 葉 松 林 生 物 量 范 圍 為2.38×106～7.52×106t，幼齡林、中齡林、近成熟 林 碳 儲(chǔ) 量 分 別 為1.19 ×107t·hm－2、2.68 ×107t·hm－2和3.76 ×107t·hm－2;碳密度分別為38.89 t·hm－2、56.09 t·hm－2和69.37 t·hm－2。隨著林齡增加，長(zhǎng)白落葉松林單位面積生物量、碳儲(chǔ)量和碳密度增大。

4 結(jié)論

在吉林省東部山地林區(qū)立地條件下，長(zhǎng)白山落葉松地上總平均生物量為215.021 kg，樹干生物量占地上總生物量67.15%，活枝、針葉、樹皮生物量分別占17.09%、6.00%和9.76%。樹干生物量最優(yōu)模型方程為Y=0.059 DBH1.32.171H0.420。樹皮、活枝、針葉和地上總生物量預(yù)測(cè)模型分別為:樹皮Y=0.046;活 枝Y=0.033;針葉Y=0.054;地 上 總 生 物 量Y=0.197。擬 合 方 程:樹 干 y =448.68x+4.4331，活枝y=105.21 x+3.944 9，針葉y=32.89 x+2.639 8，樹皮y=57.39 x+3.099 7。長(zhǎng)白落葉松不同林分生物量范圍2.38×106～ 7.52×106t，碳 儲(chǔ) 量 范 圍1.19×107～3.76×107t，碳密度范圍38.89 ～69.37 t·hm－2。

長(zhǎng)白落葉松是長(zhǎng)白山區(qū)主要用材樹種之一，適應(yīng)性強(qiáng)，生長(zhǎng)速度快，碳儲(chǔ)量積累迅速，因此，應(yīng)在適宜立地條件大力培育長(zhǎng)白落葉松林。

［1］楊麗韞，李文華.長(zhǎng)白山不同生態(tài)系統(tǒng)地下部分生物量及地下C 貯量的調(diào)查［J］.自然資源學(xué)報(bào)，2003，18(2):204－209.

［2］姜航，高菲，崔曉陽. 帽兒山次生林區(qū)土壤有機(jī)碳儲(chǔ)量及地形因子的影響［J］. 森林工程，2015，31(3):15－20.

［3］孫偉，任軍，林玉梅，等. 吉林西部森林植被碳匯功能的研究［J］.吉林林業(yè)科技，2013，42(2):14－16.

［4］馮宗煒，王效科，吳剛.中國(guó)森林生態(tài)系統(tǒng)的生物量和生產(chǎn)力［M］.北京:科學(xué)出版社，1999.

［5］Roy J，Saugier B，Mooney H A. Terrestrial global productivity［M］.San Diego，California:Academic Press，2001.

［6］閆文德，田大倫，何功秀.湖南會(huì)同第2 代杉木人工林喬木層生物量的分布格局［J］.林業(yè)資源管理，2003，12(2):5－7.

［7］李文華，鄧?yán)っ?，李飛.長(zhǎng)白山主要生態(tài)系統(tǒng)生物量生產(chǎn)量的研究［J］.森林生態(tài)系統(tǒng)研究，1981，2:34－50.

［8］郭忠玲，鄭金萍，馬元丹，等.長(zhǎng)白山各植被帶主要樹種凋落物分解速率及模型模擬的試驗(yàn)研究［J］.生態(tài)學(xué)報(bào)，2006，26(4):1037－1046.

［9］代力民，徐振邦，陳華.闊葉紅松林倒木貯量的變化規(guī)律［J］.生態(tài)學(xué)報(bào)，1999，20(3):412－416.

［10］楊麗韞，代力民.長(zhǎng)白山北坡苔蘚紅松暗針葉林倒木分解及其養(yǎng)分含量的研究［J］.生態(tài)學(xué)報(bào)，2002，22(2):185－189.

［11］Socha J，Wezyk P.Allometric equations for estimating the foliage biomass of Scots pine［J］. European journal of forest research，2007，(26)2:263－270

［12］Kajimoto，T，Hitsuma，G，Masaki，T，et al.Growth pattern analysis and stemwood production in an unmanaged old plantation of Hiba，Thujopsis dolabrata，in northern Japan［J］.Journal of Forest Research，2006，11:107－116.

［13］Kajimoto，T，Matsuura，Y，Osawa，A，et al. Size－mass allometry and biomass allocation of two larch species growing on the continuous permafrost region in Siberia［J］.Forest Ecology and Management，2006，222:314－315.

［14］王效科，馮宗煒，歐陽志云.中國(guó)森林生態(tài)系統(tǒng)的植物碳儲(chǔ)量和碳密度研究［J］.應(yīng)用生態(tài)學(xué)報(bào)，2001，12(1):13－16.

［15］劉國(guó)華，傅伯杰，方精云.中國(guó)森林碳動(dòng)態(tài)及其對(duì)全球碳平衡的貢獻(xiàn)［J］.生態(tài)學(xué)報(bào)，2000，20(5):733－740.

［16］周玉榮，于振良，趙士洞.我國(guó)主要森林生態(tài)系統(tǒng)碳貯量和碳平衡［J］.植物生態(tài)學(xué)報(bào)，2000，24(5):518－522.

［17］馬欽彥，陳遐林，王娟等.華北主要森林類型建群種的含碳率分析［J］.北京林業(yè)大學(xué)學(xué)報(bào)，2002，24(5/6):96－100.

［18］方精云，劉國(guó)華，徐嵩齡.我國(guó)森林植被的生物量和凈生產(chǎn)量［J］.生態(tài)學(xué)報(bào)，1996，16(5):498－508.

［19］孫玉軍，張俊，韓愛惠，等.興安落葉松幼中齡林的生物量與碳匯功能［J］.生態(tài)學(xué)報(bào)，2007，27(5):1756－1762.