王俊鴻 郭福濤 吳鵬飛 周麗麗 蘇漳文 馬祥慶

摘要 在收集杉木分布區(qū)中帶不同發(fā)育階段杉木人工林生物量資料基礎(chǔ)上，選擇10種不同的生物量估測(cè)模型，對(duì)杉木幼齡林、中齡林、近熟林和成熟林各器官（地上部分、葉、枝、皮、干、根）生物量與主要測(cè)樹因子進(jìn)行擬合，篩選不同發(fā)育階段杉木不同器官生物量估測(cè)模型。結(jié)果表明：擬合效果最好的為冪函數(shù)模型，其次為指數(shù)函數(shù)模型，再次為多項(xiàng)式模型；共篩選出估算杉木幼齡林、中齡林、近熟林和成熟林各器官和總生物量的最優(yōu)模型42個(gè)（包括36個(gè)不同器官生物量模型、6個(gè)全株生物量模型）；從杉木各器官生物量的擬合效果看，擬合度最高的模型均是以胸徑和樹高為自變量的模型。這些模型為分布區(qū)中帶不同發(fā)育階段杉木人工林生物量的確定和碳儲(chǔ)量評(píng)價(jià)提供科學(xué)依據(jù)。

關(guān)鍵詞 杉木；分布區(qū)中帶；生物量；估測(cè)模型；發(fā)育階段

中圖分類號(hào) S791.27 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2014）16-05104-05

生物量能很好地反映人工林與周圍環(huán)境在能量流動(dòng)和物質(zhì)循環(huán)上的相關(guān)關(guān)系，在評(píng)價(jià)人工林對(duì)全球碳平衡和溫室氣體影響方面有重要作用[1-3]。傳統(tǒng)的森林生物量估算方法主要有氣體交換法和直接收獲法，這些方法存在儀器昂貴、費(fèi)時(shí)費(fèi)力、對(duì)研究對(duì)象破壞較大及觀測(cè)不連續(xù)等問題[3-4]。應(yīng)用模型方法，通過若干易測(cè)因子如胸徑、樹高的測(cè)定來估算各樹種生物量成為森林生物量測(cè)定的重要手段[2，5]。

杉木（Cunninghamia lanceolata）是我國(guó)南方重要的用材樹種，在我國(guó)南方林區(qū)林業(yè)生產(chǎn)中占重要地位。長(zhǎng)期以來有許多學(xué)者進(jìn)行了杉木生物量的研究，并取得了大量研究成果[6-9]。由于杉木廣泛分布于我國(guó)南方16個(gè)?。▍^(qū)），不同產(chǎn)區(qū)的杉木生長(zhǎng)差異較大，以不同產(chǎn)區(qū)杉木資料建立的杉木生物量模型存在較大差異，目前缺乏適合不同分布區(qū)應(yīng)用的生物量預(yù)測(cè)模型，并且很多杉木生物量模型對(duì)不同發(fā)育階段杉木林的生物量考慮較少[6，10]，多采用D2H為自變量進(jìn)行杉木林生物量估算[11-12]。這在很大程度上限制了杉木生物量預(yù)測(cè)模型的應(yīng)用范圍，因此建立適合不同分布區(qū)應(yīng)用的杉木生物量模型成為當(dāng)前林業(yè)科學(xué)研究中需要解決的重大課題。鑒于此，筆者在收集杉木分布區(qū)中帶不同發(fā)育階段杉木人工林生物量資料的基礎(chǔ)上，選擇10種不同生物量預(yù)測(cè)模型，對(duì)杉木幼齡林、中齡林、近熟林和成熟林各器官生物量與主要測(cè)樹因子進(jìn)行擬合，篩選不同發(fā)育階段杉木不同器官生物量估測(cè)模型，為杉木分布區(qū)中帶不同發(fā)育階段杉木林生物量確定和碳儲(chǔ)量評(píng)價(jià)提供科學(xué)依據(jù)。

1 研究地區(qū)概況

杉木分布范圍大體上與我國(guó)氣候、植被區(qū)劃中的3個(gè)亞熱帶劃分范圍一致。根據(jù)杉木栽培區(qū)域自然環(huán)境、生長(zhǎng)狀況和生產(chǎn)潛力等因素的不同，我國(guó)杉木產(chǎn)區(qū)可分為北、中、南3個(gè)區(qū)劃帶[13]。杉木分布區(qū)中帶大致相當(dāng)于《中國(guó)植被》區(qū)劃的中亞熱帶地區(qū)，是我國(guó)杉木分布和栽培區(qū)域最廣的地帶。該區(qū)域北部界限由大巴山南坡以北、巫山北坡、大別山南坡，經(jīng)過黃山、天目山到杭州灣一線；西以我國(guó)西南高山林區(qū)為界線；南界從滇東南國(guó)境線起，沿桂滇邊界到右江、紅水河、沿南嶺山脈南麓、戴云山南麓直到海濱；東至東海海域。中帶分布區(qū)范圍包括浙江省、江西省、湖南省、福建省、貴州省、四川省、重慶市、云南省等。該區(qū)域氣候溫暖濕潤(rùn)，年平均溫度16 ℃，1月平均氣溫3.5～5.5 ℃，年均降水量1 200～1 700 mm，土壤為紅壤和黃棕壤[13]。

2 研究方法

2.1 數(shù)據(jù)收集 在國(guó)家林業(yè)局杉木工程研究中心長(zhǎng)期積累的杉木生物量調(diào)查數(shù)據(jù)基礎(chǔ)上，收集國(guó)內(nèi)外發(fā)表的杉木分布區(qū)中帶106篇文獻(xiàn)中有關(guān)杉木生物量數(shù)據(jù)579組。首先，對(duì)收集的生物量數(shù)據(jù)進(jìn)行預(yù)處理，刪除缺失和不全的數(shù)據(jù)。根據(jù)文獻(xiàn)記錄，確定杉木分布地點(diǎn)及林齡。其次，根據(jù)林分年齡劃分幼齡林（1～10年）、中齡林（11～20年）、近熟林（21～25年）、成熟林（26～35年）[13-14]。同時(shí)根據(jù)杉木自身生長(zhǎng)特性，將幼齡林分為幼齡林前期（1～5年）和幼齡林后期（6～10年），將中齡林分為中齡林前期（11～15年）和中齡林后期（16～20年）[15-17]。

2.2 數(shù)據(jù)統(tǒng)計(jì)與分析 利用SPSS18.0軟件和不同預(yù)測(cè)模型，對(duì)杉木各器官及總生物量與樹高、胸徑的關(guān)系進(jìn)行擬合，并對(duì)擬合效果進(jìn)行檢驗(yàn)。筆者選取線性回歸模型、多項(xiàng)式模型、冪函數(shù)模型等10種不同模型（表1）對(duì)不同發(fā)育階段杉木林生物量進(jìn)行估算[17-19]。模型中自變量包括胸徑（D）、樹高（H）等[20]。

采用判斷系數(shù)R2：1-（殘差平方和/平方矯正的總和）對(duì)生物量方程進(jìn)行擬合度評(píng)價(jià)，R2值越接近1，說明模型的擬合預(yù)測(cè)效果越好，并依此對(duì)各器官組分?jǐn)M合出的模型進(jìn)行篩選[9，21]。

3 結(jié)果與分析

3.1 分布區(qū)中帶杉木幼齡林生物量估測(cè)模型

3.1.1 幼齡林前期階段。表2為杉木幼齡林前期各器官及總生物量的擬合模型參數(shù)。從表2可看出，總生物量、地上部分生物量、葉生物量、枝生物量、皮生物量和根生物量的最優(yōu)模型均為Y=aDbHc，而干的最優(yōu)生物量模型為Y=a+bD2H+c（D2H）2。其中杉木皮的模型擬合效果最好，擬合度達(dá)0.970，干模型的擬合度為0.765、地上部分生物量模型的擬合度為0.744，而根模型的擬合效果最差。

3.1.2 幼齡林后期階段。 對(duì)幼齡林后期階段各器官及總生物量進(jìn)行擬合，結(jié)果見表3。從表3看出，不同器官生物量和總生物量的最優(yōu)模型存在差異，其中總生物量、枝生物量和根生物量的最優(yōu)模型為Y=a+bDH+c（DH）2，葉和皮的生物量最優(yōu)模型為Y=a+bD+cD2，地上部分和干的生物量最優(yōu)模型為Y=aDbHc。其中地上部分生物量模型的擬合效果最好，皮的模型擬合效果最差。

對(duì)杉木幼齡林全期生物量各器官級(jí)總生物量進(jìn)行擬合，結(jié)果見表4。從表4可以看出，不同器官生物量和總生物量的最優(yōu)模型存在差異，其中總生物量、地上、枝、皮生物量的最優(yōu)模型均為多項(xiàng)式模型，葉、干、根生物量最優(yōu)模型為Y=aDbHc。其中，地上部分生物量模型擬合度最高，為0.425，各器官及總生物量的擬合效果均不理想，因此對(duì)幼齡林進(jìn)行階段劃分以獲得更好的擬合效果。

杉木幼齡林前期和后期的生物量預(yù)測(cè)模型擬合結(jié)果（表2、3）表明：幼齡林前期和后期各器官最優(yōu)模型的差異較大，除地上部分生物量模型均為冪函數(shù)Y=aDbHc外，其他生物量模型均不同。研究顯示，前期地上部分生物量最優(yōu)模型擬合度高于后期模型擬合度。皮在前期最優(yōu)模型為冪函數(shù)Y=aDbHc，其模型擬合度明顯高于后期的多項(xiàng)式模型。此外，枝、總生物量最優(yōu)模型在前期為冪函數(shù)Y=aDbHc，不同于后期的多項(xiàng)式模型Y=a+bDH+c（DH）2 ，且前期模型擬合度高于后期。干在前期為多項(xiàng)式Y(jié)=a+bD2H+c（D2H）2，擬合度為0.765，優(yōu)于后期Y=aDbHc模型。此外，葉、根生物量后期最優(yōu)模型均為多項(xiàng)式，并且在模型擬合度方面要高于前期。

3.2 分布區(qū)中帶杉木中齡林生物量估測(cè)模型

3.2.1 中齡林前期階段。杉木中齡林前期各器官及總生物量的擬合結(jié)果見表5。從表5可看出，不同器官和總生物量的最優(yōu)模型存在顯著差異。其中，枝、根生物量的最優(yōu)模型為Y=a+bD2H+c（D2H）2，地上部分生物量的最優(yōu)模型為Y=a+bDH+c（DH）2，皮的生物量最優(yōu)模型為Y=a+bD+cD2，總生物量、葉生物量的最優(yōu)模型均為Y=aDbHc，干的生物量最優(yōu)模型為Y=a（D2H）b。皮的最優(yōu)模型擬合效果最好，葉的擬合效果最差。

3.2.2 中齡林后期階段。表6為杉木中齡林后期各器官及總生物量的的擬合結(jié)果。從表6看出，總生物量、地上部分生物量、皮生物量和干生物量的最優(yōu)模型均為Y=a+bDH+c（DH）2，其中干生物量的擬合度最高，枝生物量擬合度最低。葉、根、枝的生物量最優(yōu)模型不同，分別為Y=a+bD+cD2、Y=aDbHc和Y=a（DH）b。

從表7看出，總生物量、地上、枝、皮、干、根生物量的最優(yōu)擬合模型均為Y=aDbHc，葉生物量的最優(yōu)擬合模型為Y=aDb。其中皮、干生物量模型擬合度最高，僅為0.384，全期各器官及總生物量模型擬合效果均較差，因此，對(duì)中齡林進(jìn)行階段劃分以獲得更好的擬合效果。

杉木中齡林前期和后期階段的生物量預(yù)測(cè)模型擬合結(jié)果（表5、6）表明：前期和后期各器官最優(yōu)模型的差異較大，除地上部分、皮生物量模型均為多項(xiàng)式函數(shù)外，其他生物量模型均不同。研究顯示，前期階段地上部分、皮生物量最優(yōu)模型擬合度為0.518、0.749，明顯高于后期模型擬合度。葉生物量在后期階段為Y=a+bD+cD2，優(yōu)于前期的冪函數(shù)Y=aDbHc。枝生物量在前期的最優(yōu)模型為Y=a+bD2H+c（D2H）2，略高于后期的冪函數(shù)模型。根生物量在后期最優(yōu)模型為冪函數(shù)Y=aDbHc，不同于前期多項(xiàng)式函數(shù)Y=a+bD2H+c（D2H）2，且模型擬合度方面要高于后期。此外，干、總生物量前期的生物量最優(yōu)模型均為冪函數(shù)模型，并且在模型擬合度方面要高于后期。

3.3 分布區(qū)中帶杉木近熟林生物量估測(cè)模型 杉木近熟林各器官及總生物量的最優(yōu)模型擬合結(jié)果見表8。從表8可看出，杉木地上部分生物量、葉生物量、干生物量和根生物量的最優(yōu)模型均為Y=a+bD2H+c（D2H）2。總生物量、枝生物量的最優(yōu)模型為Y=a+bDH+c（DH）2，皮生物量最優(yōu)模型為Y=aDbHc。其中皮生物量的擬合度最高，干生物量模型擬合度為0.723，葉生物量模型的擬合效果較差。

3.4 分布區(qū)中帶杉木成熟林生物量估測(cè)模型 杉木成熟林各器官及總生物量的擬合結(jié)果見表9。從表9看出，杉木地上部分生物量、葉生物量和皮生物量的最優(yōu)模型均為Y=aDbHc?？偵锪亢透锪康淖顑?yōu)模型為Y=aebD2，干生物量的最優(yōu)模型為Y=aebDH，枝生物量的最優(yōu)模型為Y=a+bD2H+c（D2H）2。其中干模型的擬合效果最好，皮、地上部分、總生物量、根生物量的模型擬合度分別為0.838、0.822、0.818、0.811，葉的模型擬合度相對(duì)較低。

杉木近熟林和成熟林的生物量預(yù)測(cè)模型擬合效果（表8、9）表明：近熟林和成熟林各器官最優(yōu)模型的差異較大，除枝、皮生物量模型均為多項(xiàng)式、冪函數(shù)外，其他部分生物量模型均不同。研究顯示，成熟林階段枝、皮生物量最優(yōu)模型擬合度高于近熟林時(shí)期模型擬合度?？偵锪吭诔墒炝謺r(shí)期最優(yōu)模型為指數(shù)函數(shù)Y=aebD2，不同于近熟林時(shí)期的多項(xiàng)式函數(shù)，且模型擬合度高于近熟林時(shí)期。地上部分、葉生物量最優(yōu)模型在成熟林時(shí)期為均為冪函數(shù)Y=aDbHc，其模型擬合度明顯高于近熟林時(shí)期的多項(xiàng)式模型。干生物量在成熟林時(shí)期的最優(yōu)模型為指數(shù)函數(shù)Y=aebDH，不同于近熟林時(shí)期的多項(xiàng)式函數(shù)，且模型擬合度優(yōu)于近熟林時(shí)期。根生物量成熟林時(shí)期最優(yōu)模型為Y=aebD2，其模型擬合度明顯優(yōu)于近熟林時(shí)期。

4 結(jié)論

（1）在不同生物量估測(cè)模型中，冪函數(shù)模型在杉木各部位生物量的擬合效果最好，其次為指數(shù)函數(shù)模型，再次為多項(xiàng)式模型。

（2）幼齡林前期冪函數(shù)對(duì)杉木皮生物量的擬合效果最好，表達(dá)式為Y= 0.002D10.378H-8.244；多項(xiàng)式函數(shù)、冪函數(shù)對(duì)干、地上部分也有較好的擬合效果。幼齡林后期冪函數(shù)對(duì)總生物量模型的擬合效果最好，表達(dá)式為Y=2.086D0.286H0.912。杉木幼齡林處于速生階段，樹高和胸徑生長(zhǎng)旺盛，樹形是影響杉木各器官及總生物量的主要因素。

（3）中齡林前期多項(xiàng)式對(duì)杉木皮的擬合效果最好，表達(dá)式為Y=25.292-3.999D+0.198D2；中齡林后期多項(xiàng)式對(duì)干的模型擬合效果最好，表達(dá)式為Y=60.922-0.29DH+0.001（DH）2。中齡林胸徑和樹形對(duì)各器官及總生物量有較大影響。

（4）近熟林階段冪函數(shù)模型對(duì)杉木皮的生物量模型擬合效果最好，表達(dá)式為Y=0.359D0.13H1.185。多項(xiàng)式函數(shù)對(duì)干、地上部分的生物量也有較好的擬合效果。此階段樹高、胸徑生長(zhǎng)緩慢，樹形對(duì)杉木各器官及總生物量影響較大。

（5）成熟林階段指數(shù)函數(shù)對(duì)杉木干生物量的模型擬合度最好，表達(dá)式為Y=53.821e0.002D2。冪函數(shù)對(duì)皮、地上部分也有較好的擬合度，指數(shù)函數(shù)對(duì)總生物量、根生物量的擬合效果也較好。

（6）由于杉木生長(zhǎng)受氣候、立地條件、經(jīng)營(yíng)措施等多因素的影響，不同地區(qū)的杉木生物量存在一定差異[9]，這對(duì)最優(yōu)預(yù)測(cè)模型選擇和參數(shù)估算會(huì)產(chǎn)生一定影響，因此在對(duì)杉木生物量估算時(shí)，應(yīng)盡量保證研究對(duì)象的生境條件一致，使用同一分布區(qū)的生物量預(yù)測(cè)模型，以減少估測(cè)誤差。

參考文獻(xiàn)

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

[2] FANG J，CHEN A，PENG C，et al.Changes in forest biomass carbon storage in China between 1949 and 1998[J].Science，2001，292（5525）：2320-2322.

[3] 董利虎，李鳳日，賈煒瑋.黑龍江省紅松人工林立木生物量估算模型的研建[J].北京林業(yè)大學(xué)學(xué)報(bào)，2012（6）：16-22.

[4] 洪元程，徐偉強(qiáng)，葉功富，等.東南沿海木麻黃人工林生物量估測(cè)模型研究[J].浙江林業(yè)科技，2010，30（4）：66-69.

[5] PAN Y，LUO T，BIRDSEY R，et al.New estimates of carbon storage and sequestration in China's forests：effects of age.class and method on inventory.based carbon estimation[J].Climatic Change，2004，67（2/3）：211-236.

[6] 田大倫，盤宏華.第二代杉木人工林生物量的研究[J].中南林學(xué)院學(xué)報(bào)，1998，18（3）：11-16.

[7] 林開敏，鄭郁善，黃祖清，等.杉木和馬尾松幼林生物產(chǎn)量模型研究[J].福建林學(xué)院學(xué)報(bào)，1993（4）：351-356.

[8] 李燕，張建國(guó)，段愛國(guó)，等.杉木人工林生物量估算模型的選擇[J].應(yīng)用生態(tài)學(xué)報(bào)，2010，21（12）：3036-3046.

[9] 唐坤銀，唐代生.杉木生物量?jī)?yōu)化模型研究[J].林業(yè)調(diào)查規(guī)劃，2010（1）：47-52.

[10] 趙坤，田大倫.會(huì)同杉木人工林成熟階段生物量的研究[J].中南林學(xué)院學(xué)報(bào)，2000，20（1）：7-13.

[11] 鄧秋香，趙瑛，伍祿軍，等.廣西大青山12年生杉木人工林的生物生產(chǎn)力[J].廣西林業(yè)科學(xué)，2009，37（4）：187-190.

[12] 康冰，劉世榮，張廣軍，等.廣西大青山南亞熱帶馬尾松，杉木混交林生態(tài)系統(tǒng)碳素積累和分配特征[J].生態(tài)學(xué)報(bào)，2006，26（5）：1320-1329.

[13] 俞新妥.杉木栽培學(xué)[M].福州：福建科學(xué)技術(shù)出版社，1966.

[14] 王賽專.會(huì)同不同年齡階段第2代杉木人工林碳貯量動(dòng)態(tài)特征研究[D].長(zhǎng)沙：中南林業(yè)科技大學(xué)，2013.

[15] 傅軍，蔣建屏，彭立平，等.皖南杉木標(biāo)準(zhǔn)林分蓄積生長(zhǎng)階段劃分探討[J].浙江林學(xué)院學(xué)報(bào)，1992，9（1）：24-28.

[16] 鄭勇平，朱浩，徐煥農(nóng)，等.杉木生長(zhǎng)過程特征和生長(zhǎng)期劃分[J].浙江林學(xué)院學(xué)報(bào)，1991，8（2）：219-226.

[17] 張卓文.柳杉生長(zhǎng)過程分析及生長(zhǎng)階段劃分[J].中南林學(xué)院學(xué)報(bào)，2003，23（2）：46-51.

[18] 曾子純，劉思齊.馬尾松水保林經(jīng)營(yíng)模式與經(jīng)營(yíng)數(shù)表研究[M].北京：中國(guó)林業(yè)出版社，2003.

[19] 曾慧卿，劉琪璟，馬澤清，等.千煙洲灌木生物量模型研究[J].浙江林業(yè)科技，2006，26（1）：13-17.

[20] 王維楓，雷淵才，王雪峰，等.森林生物量模型綜述[J].西北林學(xué)院學(xué)報(bào)，2008，23（2）：58-63.

[21] PARRESOL B R.Assessing tree and stand biomass：a review with examples and critical comparisons[J].Forest Science，1999，45（4）：573-593.