高 靈，王瑞輝，王 睿，謝耀堅(jiān)，吳志華＊

(1.中南林業(yè)科技大學(xué) 森林培育國家級(jí)重點(diǎn)實(shí)驗(yàn)室，湖南 長沙410000；2.國家林業(yè)局桉樹研究開發(fā)中心，廣東 湛江 524022)

幼年赤桉胸徑與冠幅、樹高、材積的相關(guān)性分析

高 靈1,2，王瑞輝1，王 睿1,2，謝耀堅(jiān)2，吳志華2＊

(1.中南林業(yè)科技大學(xué) 森林培育國家級(jí)重點(diǎn)實(shí)驗(yàn)室，湖南 長沙410000；2.國家林業(yè)局桉樹研究開發(fā)中心，廣東 湛江 524022)

對(duì)試驗(yàn)地的1 470株赤桉進(jìn)行樹高、胸徑、冠幅的測量，并計(jì)算出樣木的單株材積。把胸徑分別和冠幅、樹高及株材積進(jìn)行相關(guān)性分析并且建立數(shù)學(xué)模型，用SPSS軟件對(duì)所選模型進(jìn)行曲線估計(jì)。結(jié)果表明：其中冪方程的 R2最大，F(xiàn)值亦為最大，說明赤桉胸徑與樹高的冪關(guān)系顯著，可確定赤桉胸徑與樹高的最優(yōu)回歸方程為H=1.804 D0.673。胸徑—冠幅，胸徑—材積的最優(yōu)模型分別為CW=0.674 D0.561，V=0.000 161 4 D2.341。分別對(duì)3組最優(yōu)模型進(jìn)行適應(yīng)性檢驗(yàn)，結(jié)果表明：材積的3個(gè)最優(yōu)回歸模型預(yù)測誤差均在±3%以內(nèi)，方程預(yù)測精度較高，可用于估算立木樹高、冠幅、材積。

赤桉；胸徑；冠幅；樹高；立木材積；回歸模型

桉樹(Eucalyptus)是世界人工造林的第二大樹種，因其具有速生豐產(chǎn)、適應(yīng)性強(qiáng)、用途廣、輪伐期短等優(yōu)良特性，在我國南亞熱帶作為短周期工業(yè)用材林的高產(chǎn)商品材樹種已迅速發(fā)展起來，且有顯著成效[1]。赤桉(E. camaldulensis)是澳大利亞大陸分布最廣泛的樹種，具有優(yōu)良的適應(yīng)性、豐富的遺傳多樣性。盡管其生長速度、干形的優(yōu)良程度不及熱帶、亞熱帶的主要造林樹種尾葉桉(E. urophylla)、巨桉(E. grandis)，但在抗風(fēng)、耐旱、耐瘠薄方面具有獨(dú)特的優(yōu)勢，是桉樹遺傳育種的重要材料，也是世界上栽培較廣的桉樹種之一[2]。如何對(duì)赤桉的各項(xiàng)生長指標(biāo)進(jìn)行有效快捷的測量并用于生產(chǎn)研究評(píng)估成了一個(gè)值得研究的內(nèi)容。胸徑和樹高是2個(gè)重要的生長性狀，也是單株材積計(jì)算、林分生長和收獲評(píng)估中不可或缺的數(shù)據(jù)[3]。研究樹高與胸徑的關(guān)系意義重大[4-7]。立木材積通常都是利用樹高和胸徑計(jì)算而得[8]。樹冠是樹木進(jìn)行光合作用和呼吸作用的重要場所，冠幅和冠長率直接影響樹木的生活力和生產(chǎn)力[9]，因此冠幅是林木生長中重要的性狀。胸徑的測量比較簡單，數(shù)據(jù)可靠，而樹高、冠幅和單株材積則相對(duì)較難準(zhǔn)確獲得。在以往的研究中，樹高和胸徑的相關(guān)性已得到多次驗(yàn)證，樊艷文等[10]對(duì)北京栓皮櫟(Quercus variabilis)林胸徑、樹高相關(guān)生長關(guān)系的分析中，證明了各種生境因子和林分因子(地形、土壤、光照、林齡等)對(duì)胸徑、樹高生長關(guān)系的冪指數(shù)均無顯著影響；魏京[11]對(duì)胸徑、樹高的非線性模型進(jìn)行了擬合性研究；劉敏等[8]研究了遼寧樟子松(Pinus sylvestris)人工林樹木胸徑與樹高、材積的相關(guān)關(guān)系，得出胸徑與材積最優(yōu)回歸方程為V=0.000 165 D2.27886。本研究用胸徑分別與冠幅、樹高、單株材積進(jìn)行相關(guān)性分析并建立數(shù)學(xué)模型，為赤桉林地通過測量胸徑而推算出樹高、冠幅和立木材積的估計(jì)值提供依據(jù)，同時(shí)為赤桉的遺傳育種研究和生產(chǎn)提供一定的技術(shù)支持。

1 試驗(yàn)地與試驗(yàn)林概況

1.1 試驗(yàn)地概況

試驗(yàn)地在南方國家級(jí)林木種苗示范基地內(nèi)，位于廣東省遂溪縣嶺北鎮(zhèn)附近，地理位置為21°30′N，111°38′E，處于濕潤大區(qū)雷瓊熱帶北緣,為熱帶北緣季風(fēng)氣候；年均溫23.6℃，絕對(duì)最低氣溫2.8℃，絕對(duì)最高氣溫37.0℃，月平均溫度18.0℃，夏季降水,年均降水量1 567 mm，5—9月降水占全年的85%，臺(tái)風(fēng)頻繁，平均每年遭受10級(jí)以上臺(tái)風(fēng)襲擊2 ~ 3次，最大風(fēng)力大于35 m·s-1[12]。

1.2 試驗(yàn)林概況

2012年5月造林，實(shí)際造林為115個(gè)家系，每個(gè)家系4個(gè)重復(fù)，共1 470株?試驗(yàn)設(shè)計(jì)為完全隨機(jī)區(qū)組，4個(gè)區(qū)組，株行距2 m×3 m?面積約14.3 hm2?本試驗(yàn)林地土壤由花崗巖發(fā)育而成，為肥沃的磚紅壤，土層厚度大于1 m?試驗(yàn)林地較為平緩，采用機(jī)械全墾整地(墾深30 cm)，鏵犁開溝(行距3 m開溝，溝深40 cm)?

2 研究方法

2014年4月，對(duì)赤桉試驗(yàn)林進(jìn)行樹高、胸徑、冠幅的林地調(diào)查。分別對(duì)樣木的胸徑和冠幅進(jìn)行東西、南北兩個(gè)方向的測量，取兩個(gè)方向的平均值得到胸徑和冠幅的值。

2.1 測定方法

以Vertex IV測高測距儀對(duì)樹高(H)進(jìn)行測定，分別以卷尺和胸徑尺對(duì)冠幅(CW)和胸徑(D)進(jìn)行測定。

單株材積(V)采用廣西林業(yè)勘查設(shè)計(jì)院總結(jié)的廣西速生桉單株材積公式，具體如下：

式(1)中的V為單株材積，D為胸徑，H為樹高，C0=0.000 109 154 150，C1=1.878 923 70，C2=0.005 691 855 03，C3=0.652 598 05，C4=0.007 847 535 07。

2.2 數(shù)據(jù)整理

通過立木材積公式計(jì)算出赤桉單株材積，分別以統(tǒng)計(jì)軟件SPSS 17(Statistical Product and Service Solutions)作胸徑—樹高、胸徑—冠幅、胸徑—材積的散點(diǎn)圖，去除個(gè)別異常的數(shù)據(jù)，共得到1 470株赤桉樣木用來建立回歸模型和模型檢驗(yàn)。樣木各個(gè)變量因子的取值范圍分別如下：胸徑0.5 ~ 10.0 cm，樹高0.9 ~ 9.6 m，冠幅0.20 ~ 4.15 m。為保證建模樣木和檢驗(yàn)樣木在選取時(shí)均勻合理，將所有樣木按胸徑大小排序，每4株家系為1組，抽取每組第1株作為檢驗(yàn)樣木。剩余樣木則選作建模樣木。最終總樣木為1 470，建模樣木和檢測樣木分別為1 102 (占3/4)和368(1/4)。

2.3 模型的選擇

以建模樣木胸徑為自變量，樹高(H)、冠幅(CW)、單株材積(V)分別為因變量，作D-H、D-CW、D-V的散點(diǎn)圖。通過對(duì)散點(diǎn)圖的分析，選取線性函數(shù) y=b0+b1x、對(duì)數(shù)函數(shù) y=b0+b1lnx、倒數(shù)函數(shù)y=b0+b1/x、二次函數(shù) y=b0+b1x+b2x2、三次函數(shù)y=b0+b1x+b2x2+b3x3、復(fù)合函數(shù) y=b0b1x、冪函數(shù)y=b0xb1、S函數(shù)y=e(b0+b1/x)、生長函數(shù)y =e(b0+b1/x)、指數(shù)函數(shù)y=b0eb1x8個(gè)曲線模型作為擬合模型。

2.4 適應(yīng)性檢驗(yàn)方法

將各檢驗(yàn)樣木實(shí)測的胸徑值導(dǎo)入所建立的胸徑—樹高、胸徑—冠幅、胸徑—材積最優(yōu)模型中，通過計(jì)算，可以得出檢驗(yàn)樣木的樹高、冠幅和材積的估計(jì)值。通過對(duì)比檢驗(yàn)樣木中樹高、冠幅和材積的實(shí)際值，分別計(jì)算出各個(gè)最優(yōu)模型的總相對(duì)誤差、平均相對(duì)誤差和預(yù)估精度，用來確定所建立的最優(yōu)模型是否符合標(biāo)準(zhǔn)。各指標(biāo)估計(jì)值檢驗(yàn)計(jì)算按孟憲宇[13]的方法。

3 結(jié)果與分析

3.1 D-H回歸模型

以建模樣木胸徑(D)為自變量、樹高(H)為因變量，用SPSS 17.0分析軟件作D-H散點(diǎn)圖(圖1)。在SPSS中以所選模型進(jìn)行曲線估計(jì)，當(dāng)R2最大、F值亦最大時(shí)，擬合的模型為最優(yōu)模型，可確定為赤桉D-H的最優(yōu)回歸方程。通過所選模型對(duì)D-H關(guān)系在軟件中進(jìn)行擬合，各個(gè)擬合曲線的參數(shù)見表1。結(jié)果顯示，所有方程除倒數(shù)方程R2外均大于0.700，P值均小于0.05，F(xiàn)值遠(yuǎn)大于F(1,1 100)=3.84，表明所有的擬合模型均有意義。其中冪方程的R2最大，為0.840，F(xiàn)值為5 758.798，比其他擬合模型的F值都高，說明赤桉D-H的冪關(guān)系顯著。二次、三次方程的R2均為0.819，但是F值卻偏低。線性方程的R2和F值均較高，是回歸性比較好的方程?？紤]到方程的簡捷，取D-H的最優(yōu)模型為冪方程H=1.804 D0.673。

圖1 赤桉D-H散點(diǎn)圖以及各擬合線

表1 D-H各類模型曲線估計(jì)結(jié)果

3.2 D-CW回歸模型

把建模樣木的胸徑(D)、冠幅(CW)分別選作自變量、因變量，用SPSS軟件作D-CW散點(diǎn)圖(圖2)。在SPSS中以所選模型進(jìn)行曲線估計(jì)，當(dāng)R2較大、F值亦大時(shí)，擬合的模型的可靠性主要是看方差檢驗(yàn)結(jié)果(P值)，模型的優(yōu)選性預(yù)測誤差小的模型。最后確定赤桉D-CW的最優(yōu)回歸方程。將所選模型對(duì)D-CW關(guān)系在軟件中進(jìn)行擬合，各個(gè)擬合曲線的參數(shù)詳見表2。結(jié)果顯示，所有方程除倒數(shù)方程R2外均大于0.500，P值均小于0.05，F(xiàn)值遠(yuǎn)大于F(1,1 100)=3.84，表明擬合出的所有模型均可靠。其中冪方程的R2最大，為0.663，F(xiàn)值為2 161.848，比其他擬合模型的F值都高，說明赤桉D-CW呈顯著冪函數(shù)關(guān)系。僅次于冪關(guān)系的是S方程，二者的R2均大于0.6。取D-CW的最優(yōu)模型為冪方程CW=0.674 D0.561。

圖2 赤桉D-CW散點(diǎn)圖以及各擬合線

表2 D-CW各類模型曲線估計(jì)結(jié)果

3.3 D-V回歸模型

從圖3和表3可知，除了倒數(shù)、對(duì)數(shù)和S曲線外，其他模型的R2均在0.900以上，說明赤桉的胸徑(D)與材積(V)之間的相關(guān)關(guān)系非常緊密。F檢驗(yàn)統(tǒng)計(jì)量的觀察值均大于臨界值F(1,109)=3.930，P值均小于0.001，表明所建的模型符合統(tǒng)計(jì)學(xué)要求。其中冪函數(shù)的R2為0.992，最接近1.000，F(xiàn)值為最大，等于128 989.971，顯示D-V的冪關(guān)系極其顯著，冪函數(shù)擬合度最好，明顯優(yōu)于其他方程。故選擇冪曲線為赤桉D-V最優(yōu)模型為V=0.000 161 4 D2.341。

圖3 赤桉D-V散點(diǎn)圖以及各擬合線

表3 D-V各類模型曲線估計(jì)結(jié)果

3.4 樹高、胸徑、冠幅與材積相關(guān)性分析

對(duì)樹高、胸徑、冠幅與材積進(jìn)行相關(guān)分析，結(jié)果表明(表4)：材積分別與胸徑、冠幅、樹高成極顯著的正相關(guān)關(guān)系。

表4 胸徑、冠幅、樹高與材積兩兩相關(guān)性分析

3.5 模型適應(yīng)性檢驗(yàn)結(jié)果

經(jīng)過驗(yàn)算，赤桉D-H、D-CW、D-V模型精度檢驗(yàn)統(tǒng)計(jì)指標(biāo)如表5。結(jié)果顯示，赤桉D-H模型適應(yīng)性檢驗(yàn)的總相對(duì)誤差和平均相對(duì)誤差都在±1%以內(nèi)，精度高達(dá)98.76%，驗(yàn)證了赤桉D-H具有顯著的冪函數(shù)關(guān)系。D-CW的適應(yīng)性檢驗(yàn)誤差小于±2%，精度為98.30%，樹高和胸徑也具有明顯的冪函數(shù)關(guān)系。D-V最優(yōu)模型的適應(yīng)性檢驗(yàn)總相對(duì)誤差為0.0 262%，平均相對(duì)誤差為0.89%，均在±1%以內(nèi)，精度為98.96%，是精度最高的模型，說明胸徑和材積之間的關(guān)系十分顯著，可以利用赤桉的胸徑直接對(duì)單株材積進(jìn)行計(jì)算。

表5 各模型適應(yīng)性檢驗(yàn)評(píng)價(jià)統(tǒng)計(jì)指標(biāo) %

4 結(jié)論

本研究結(jié)果表明，南方國家級(jí)林木種苗示范基地內(nèi)赤桉胸徑分別與樹高、冠幅、材積呈顯著的冪函數(shù)關(guān)系。D-H、D-CW、D-V冪函數(shù)關(guān)系分別為H=1.804 D0.673、W=0.674 D0.561、V=0.000 161 4 D2.341。高于林業(yè)部頒布的技術(shù)規(guī)定的誤差控制范圍[14]，而模型中的誤差均在±2%以內(nèi)，精度在98%以上，模型擬合度相當(dāng)好，適應(yīng)性檢驗(yàn)誤差符合林業(yè)部規(guī)定范圍，可直接運(yùn)用到當(dāng)?shù)貙?shí)際生產(chǎn)中。這些模型不但可為赤桉的遺傳研究以及林業(yè)生產(chǎn)中估測生長量、采伐量、盜伐量、誤伐量等提供依據(jù)，也為林業(yè)部門提供了一種有效計(jì)算赤桉樹高、冠幅和材積的方法，即可利用赤桉的胸徑—材積模型，編制便于當(dāng)?shù)厥褂玫某噼褚辉貜讲姆e表，此法較查閱林木二元材積表計(jì)算材積更方便快捷[15]。

[1] 劉友全,劉加林,潘天玲.赤桉在湖南的抗寒與生長適應(yīng)性[J].中南林學(xué)院學(xué)報(bào),2000,20(3):86-89.

[2] Florence R G.對(duì)桉樹作為外來樹種的感性認(rèn)識(shí)[C]//中國林學(xué)會(huì)桉樹專業(yè)委員會(huì).國際桉樹學(xué)術(shù)研討會(huì)論文集.北京:中國林業(yè)出版社,1992.

[3] 李?？?法蕾.基于分級(jí)的全國主要樹種樹高－胸徑曲線模型[J].林業(yè)科學(xué),2011,47(10):83-90.

[4] 鄭傳英.四川榿木林分各主要生長因子間相關(guān)性研究[J].安徽農(nóng)學(xué)通報(bào),2010,16(17):161-162.

[5] 王曉林,郭斌.柞樹樹高與胸徑相關(guān)關(guān)系的研究[J].森林工程,2012,28(6):18-21.

[6] 王利,豐震,白世紅,等.麻櫟樹高與胸徑相關(guān)關(guān)系的研究[J].山東林業(yè)科技.2005(4):33-34.

[7] 趙敏,丁慧勇.上海水杉樹高與胸徑關(guān)系模式分析[J].上海師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,38(5):531-535.

[8] 劉敏,張日升,徐貴軍,等.遼寧樟子松人工林樹木胸徑與樹高、材積的相關(guān)關(guān)系[J].防護(hù)林科技,2010(1):36-37,41.

[9] 符利勇,孫華,張會(huì)儒,等.不同郁閉度下胸高直徑對(duì)杉木冠幅特征因子的影響[J].生態(tài)學(xué)報(bào),2013,33(8):2434-2443.

[10] 樊艷文,王襄平,曾令兵,等.北京栓皮櫟林胸徑-樹高相關(guān)生長關(guān)系的分析[J].北京林業(yè)大學(xué)學(xué)報(bào),2011,33(6): 146-150.

[11] 魏京.林木胸徑與樹高的關(guān)系研究[J].湖北民族學(xué)院學(xué)報(bào)(自然科學(xué)版),2014,32(2):190-192．

[12] 張華林,李天會(huì),吳志華,等.不同林分密度對(duì)桉樹幼林木材材性的影響[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2010,30(6): 85-91.

[13] 孟憲宇.測樹學(xué)[M].北京:中國林業(yè)出版社,1997.

[14] 中華人民共和國林業(yè)部.林業(yè)專業(yè)調(diào)查主要技術(shù)規(guī)定[M].北京:中國林業(yè)出版社,1990.

[15] 伍靜,吳英,龍楚云,等.尾葉桉地徑與胸徑、樹高、材積相關(guān)性分析[J].南方農(nóng)業(yè)學(xué)報(bào),2013,44(6):979-983.

Correlation Analyses Between DBH, Crown W idth, Height and Volume of Eucalyptus camaldulensis

GAO Ling1,2, WANG Rui-hui1, WANG Rui1,2, XIE Yao-jian2, WU Zhi-hua2
(1. The Key Lab of Silviculture of State, Central South University of Forestry and Technology, Changsha 410000, Hunan, China; 2. China Eucalypt Research Centre, Zhanjiang 524022, Guangdong, China)

Data were collected from 1 470 Eucalyptus camaldulensis trees including tree height, diameter at breast height and crown w idth measurements. Individual tree volume was calculated from these parameters for each individual tree. Analyses of the correlations between diameter at breast height, crown w idth, tree height and individual tree volume were carried out using SPSS software. The results show that the largest of which R2, F value of the maximum power equation also shows diameter at breast height was significantly higher power relationship can be determined E. camaldulensis optimal diameter and tree height regression equation H = 1.804 D0.673. Similarly come DBH - crown w idth, DBH - volume optimal model was CW = 0.674 D0.561, V = 0.000 161 4 D2.341. To ensure the accuracy of the research and analyses of, respectively, three groups of the far right model for adaptive test results showed that three best regression model to predict volume are in error w ithin ± 3.0%, higher prediction accuracy equation can be used to estimate the standing tree height, crown w idth and volume.

Eucalyptus camaldulensis; DBH; crown w idth; tree height; standing timber volume; regression model

S758.5

A

2014-11-14

廣東省林業(yè)科技創(chuàng)新專項(xiàng)資金項(xiàng)目 (2011KJCX019);林業(yè)重大公益性行業(yè)科研專項(xiàng)(201104003)

高靈(1990— )，女，碩士生，主要從事森林培育研究.E-mail:812854076@qq.com

＊通訊作者：吳志華(1974— )，男，碩士，高級(jí)工程師，主要從事林木逆境生理研究.E-mail:wzhua2889@163.com