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

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

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

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

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

2 研究方法

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

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

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

?

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

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

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

Y 為各器官生物量;

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

利用樹干直徑和樹高建立各器官生物量兩變量回歸模型。方程形式為

式中:X1為樹干直徑;

X2為樹高;

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

方程擬合優(yōu)度利用決定系數(shù)(R2)和信息指數(shù)(AIC)進行評價。

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

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

B 為地上生物量;

Cc為含碳率。

碳密度計算公式為

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

C 為碳儲量(t);

S 為森林面積(hm2)。

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

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

3 結(jié)果與分析

3.1 長白落葉松地上生物量異速生長模型

長白落葉松地上生物量分配如圖1 所示。

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

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

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

從表2 中可以看出，所有模型F 檢驗均顯著(p ＜0.001)，基于胸徑的單變量模型可解釋量均超過90%，其中樹干生物量、樹皮生物量和地上總生物量回歸模型可解釋量達到99%。利用胸徑和樹高建立的樹干生物量模型，決定系數(shù)Radj2 最大，信息指數(shù)AIC 最小。活枝和針葉生物量回歸模型可解釋量超過93%，與其他3 個回歸模型相比信息指數(shù)AIC 值最小。因此，利用胸徑和樹高可獲得樹干生物量最優(yōu)模型，模型方程為Y=0.059 DBH1.32.171H0.420。胸徑對樹皮、活枝、針葉和地上總生物量預測較重要，模 型 方 程 分 別 為:樹 皮，地上總生物量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值達到0.9 以上。方精云［18］和孫玉軍［19］研究顯示林木生物量與蓄積量之間存在線性關(guān)系，與本文結(jié)果一致。生物量與蓄積量之間的線性關(guān)系可以作為估測森林地上生物量的基礎。

3.3 長白落葉松林碳儲量與碳密度

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

表3 長白落葉松林碳儲量和碳密度Tab.3 Carbon storage and carbon density of Larix olgensis forest

從表3 中可以看出，長白山地區(qū)長白落葉松幼齡林面積為3.06×105hm2，蓄積量為1.66×106m3;中齡林面積為4.71×105hm2，蓄積量為2.72×106m3;近成熟林面積為5.42×105hm2，蓄積量為2.34×106m3。長白落 葉 松 林 生 物 量 范 圍 為2.38×106～7.52×106t，幼齡林、中齡林、近成熟 林 碳 儲 量 分 別 為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。隨著林齡增加，長白落葉松林單位面積生物量、碳儲量和碳密度增大。

4 結(jié)論

在吉林省東部山地林區(qū)立地條件下，長白山落葉松地上總平均生物量為215.021 kg，樹干生物量占地上總生物量67.15%，活枝、針葉、樹皮生物量分別占17.09%、6.00%和9.76%。樹干生物量最優(yōu)模型方程為Y=0.059 DBH1.32.171H0.420。樹皮、活枝、針葉和地上總生物量預測模型分別為:樹皮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。長白落葉松不同林分生物量范圍2.38×106～ 7.52×106t，碳 儲 量 范 圍1.19×107～3.76×107t，碳密度范圍38.89 ～69.37 t·hm－2。

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

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

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

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

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

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

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

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

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

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

［10］楊麗韞，代力民.長白山北坡苔蘚紅松暗針葉林倒木分解及其養(yǎng)分含量的研究［J］.生態(tài)學報，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］王效科，馮宗煒，歐陽志云.中國森林生態(tài)系統(tǒng)的植物碳儲量和碳密度研究［J］.應用生態(tài)學報，2001，12(1):13－16.

［15］劉國華，傅伯杰，方精云.中國森林碳動態(tài)及其對全球碳平衡的貢獻［J］.生態(tài)學報，2000，20(5):733－740.

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

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

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

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