楊盛揚(yáng)，曾思齊，龍時(shí)勝，王帥玲

（中南林業(yè)科技大學(xué) 林學(xué)院，湖南 長(zhǎng)沙 410004）

在森林結(jié)構(gòu)中樹(shù)高和直徑是最基本的林分因子，高徑比是樹(shù)高和直徑的比值，它是說(shuō)明林木樹(shù)高和直徑相關(guān)關(guān)系的一個(gè)統(tǒng)計(jì)指標(biāo)[1]，也反映樹(shù)干圓滿度的好壞，直接影響到木材的出材率和經(jīng)濟(jì)價(jià)值[2]。高徑比還可用于確定立木的生長(zhǎng)率，在營(yíng)林工作中可根據(jù)不同的經(jīng)營(yíng)目標(biāo)控制樹(shù)高和直徑的生長(zhǎng)，從而達(dá)到最佳的經(jīng)營(yíng)效果。廖澤釗等[1]利用七坡林場(chǎng)的二調(diào)數(shù)據(jù)得出林木高徑比的增減變化程度與樹(shù)種的生物學(xué)特性、立地條件和林分密度有關(guān)，高徑比與直徑成反比，與年齡成正比；王彩云等[3]分別立地指數(shù)、林齡和疏密度對(duì)云南松天然林的高徑比進(jìn)行研究，提出不同立地指數(shù)林分的高徑比大致相等的結(jié)論；李旭等[4]用高徑比法確定落葉松人工林的生長(zhǎng)率，認(rèn)為此方法外業(yè)操作簡(jiǎn)單、測(cè)量精度較高，可以滿足生產(chǎn)的要求，也提出了各種立地條件下的高徑比之間差異不顯著；張更新[5]探討了油松林高徑比的變化規(guī)律，并得出了高徑比的變動(dòng)范圍。很多學(xué)者研究高徑比時(shí)大都從宏觀上探討高徑比的總體變化趨勢(shì)，其結(jié)論不一，較少涉及到單木高徑比的具體變化。樹(shù)高和直徑相關(guān)關(guān)系的研究一直是林業(yè)界的研究熱點(diǎn)，常見(jiàn)的研究樹(shù)種有杉木、馬尾松、落葉松等[6-10]，國(guó)外學(xué)者也致力于使用清查和樣地?cái)?shù)據(jù)研究樹(shù)高-直徑關(guān)系模型等[11-12]，國(guó)內(nèi)的學(xué)者如李際平等基于人工神經(jīng)網(wǎng)絡(luò)，建立單木樹(shù)高-直徑生長(zhǎng)的BP 模型[13]，其模型使用方便，精度較高；李海奎等[14]使用5 種樹(shù)高-直徑曲線模型模擬全國(guó)主要樹(shù)種的樹(shù)高-直徑關(guān)系，其結(jié)果表明模型總體穩(wěn)定參數(shù)可靠。研究樹(shù)高-直徑關(guān)系模型應(yīng)從實(shí)際需要出發(fā)，操作簡(jiǎn)便且精度較高，太過(guò)于復(fù)雜的模型不利于推廣應(yīng)用。

馬尾松Pinus massoniana屬松科Pinaceae松屬Pinus植物，它適應(yīng)力強(qiáng)，分布廣，常被用于荒山造林的先鋒樹(shù)種[15]，是我國(guó)重要的用材樹(shù)種和工業(yè)原料林樹(shù)種[16]。目前國(guó)內(nèi)少見(jiàn)有關(guān)馬尾松人工林高徑比變化的研究，以馬尾松人工林為研究對(duì)象，探討樹(shù)高、直徑、年齡與高徑比的關(guān)系，以及研究總結(jié)在不同立地指數(shù)、不同林分密度下馬尾松人工林高徑比的變化規(guī)律。建立馬尾松人工林的樹(shù)高-直徑關(guān)系模型，為馬尾松人工林的經(jīng)營(yíng)培育提供科學(xué)的理論指導(dǎo)，也為外業(yè)精確、簡(jiǎn)便地測(cè)定樹(shù)高提供新的方法。

1 材料與方法

1.1 樣地設(shè)置和數(shù)據(jù)處理

試驗(yàn)材料取自江西、廣西、重慶、河南、福建5 個(gè)?。▍^(qū)、市）的馬尾松人工林樣地，設(shè)置臨時(shí)樣地255 塊并取平均木解析木255 株，解析木概況如表1。

設(shè)置10 至18 共5 個(gè)指數(shù)級(jí)，每個(gè)指數(shù)級(jí)分別在255 塊樣地中選擇疏密度相近、年齡相同的5株解析木數(shù)據(jù)，并取其各齡階樹(shù)高、直徑、高徑比的平均數(shù)，分別記為B1～B5，用以研究不同立地指數(shù)對(duì)高徑比的影響；設(shè)置三種株數(shù)密度級(jí)，每種株數(shù)密度分別選擇立地指數(shù)相同、年齡相同的5 株解析木數(shù)據(jù)，并且同一密度級(jí)的5 株解析木所在樣地的株數(shù)密度相差不超過(guò)每公頃200 株，并取其平均數(shù)，分別記為C1、C2、C3，用以研究不同林分密度對(duì)高徑比的影響，樣地概況表見(jiàn)表2。

表2 馬尾松人工林樣地概況Table 2 General situation of Pinus massoniana sample plots

1.2 林分相關(guān)因子的確定方法

對(duì)于相同立地指數(shù)和相同年齡的林分，為了準(zhǔn)確研究密度對(duì)高徑比的影響，以單位面積上的株數(shù)作為密度指標(biāo)。在研究不同立地指數(shù)對(duì)高徑比的影響時(shí)，由于立地指數(shù)對(duì)直徑有一定影響，使同株數(shù)密度同年齡林分的單木直徑有差異，勢(shì)必對(duì)高徑比也造成影響，因此采用的林分密度為疏密度，疏密度是現(xiàn)實(shí)林分每公頃直徑（或蓄積）與相同立地下標(biāo)準(zhǔn)林分每公頃直徑斷面積（或蓄積）之比，引用1958年由原林業(yè)部編制的《森林調(diào)查員手冊(cè)》中疏密度1.0 的馬尾松林分?jǐn)嗝娣e蓄積量表計(jì)算出各樣地疏密度，立地指數(shù)是根據(jù)林分中優(yōu)勢(shì)木標(biāo)準(zhǔn)年齡時(shí)的平均高來(lái)確定。

1.3 模型精度評(píng)價(jià)

用決定系數(shù)R2、預(yù)估精度P值兩項(xiàng)評(píng)價(jià)指標(biāo)對(duì)模型擬合結(jié)果進(jìn)行評(píng)價(jià)，R2和P值越接近1 說(shuō)明模型擬合精度越高。表達(dá)式如下：

式中：yi為實(shí)際值；為理論值；為平均值；n為樣本數(shù)，k為模型中的參數(shù)個(gè)數(shù)。

1.4 分析方法

繪制解析木的高徑比與年齡、直徑、樹(shù)高關(guān)系散點(diǎn)圖，分析這三種關(guān)系的變化趨勢(shì)。根據(jù)高徑比與樹(shù)高、直徑、年齡關(guān)系的帶平滑線的散點(diǎn)圖，研究三個(gè)因子的變化規(guī)律。

繪制高徑比與年齡、樹(shù)高，年齡與樹(shù)高、直徑的關(guān)系圖；分別對(duì)不同立地指數(shù)的高徑比、樹(shù)高、直徑進(jìn)行單因素方差分析和LSD 檢驗(yàn)，由于對(duì)不同立地指數(shù)下的高徑比進(jìn)行方差分析時(shí)，其方差同質(zhì)性檢驗(yàn)的顯著性為0，因此采用方差分析中的Kruskal-Wallis 檢驗(yàn)。

選擇同立地指數(shù)同年齡的不同株數(shù)密度的解析木數(shù)據(jù)各5 組，相同密度級(jí)的各組數(shù)據(jù)株數(shù)密度相差不超過(guò)每公頃200 株，取其平均數(shù)。繪制高徑比與年齡、樹(shù)高關(guān)系，年齡與樹(shù)高、直徑的關(guān)系圖；分別對(duì)不同立地指數(shù)的高徑比、樹(shù)高、直徑進(jìn)行方差分析。

2 結(jié)果與分析

2.1 樹(shù)高、直徑、年齡與高徑比的關(guān)系

從圖1可知，從總體上高徑比隨樹(shù)高、直徑、年齡的增大而呈“L”型；在樹(shù)高約6 m、直徑約4 cm 之前，高徑比值分布范圍約在0.5～3.5，不同解析木在相同樹(shù)高或相同直徑時(shí)，高徑比的最大差值達(dá)到2.78；高徑比值主要集中在0.5～1.5之間，高徑比最小值隨樹(shù)高、直徑、年齡的增大而減小，最后有穩(wěn)定的趨勢(shì)。

圖1表明，高徑比與樹(shù)高、直徑、年齡的關(guān)系都可分為兩個(gè)階段，第一階段主要發(fā)生在幼齡林時(shí)期，此時(shí)林分尚未郁閉或剛郁閉不久，高徑比急劇減小至最小值；第二階段主要發(fā)生在中齡林，此時(shí)林分已完全郁閉，高徑比規(guī)律地緩慢增大。

圖1 高徑比與樹(shù)高/直徑/年齡關(guān)系Fig.1 Relationship between HD and tree height/DBH/age

2.2 不同立地指數(shù)對(duì)高徑比的影響

從圖2看出，相同疏密度相同年齡的林分第8年以前，各指數(shù)級(jí)的高徑比在總體上與年齡成反比，其曲線的斜率也不盡相同且無(wú)規(guī)律，高徑比大都在第8年達(dá)到最小值。在第8年后，采用單因素方差分析中Kruskal-Wallis 檢驗(yàn)得出，10 與14、16、18，12 與18 指數(shù)級(jí)的高徑比達(dá)到極顯著差異（表3），指數(shù)級(jí)差距越大高徑比差異越大，并且這種差異隨年齡的增大而進(jìn)一步加大。

圖2表明，第8年至20年立地指數(shù)越大，樹(shù)高生長(zhǎng)曲線的斜率越大，不同指數(shù)級(jí)的樹(shù)高差異就越顯著（表4），18 和10、12 指數(shù)級(jí)的樹(shù)高顯著性P值小于0.05；在第4年時(shí)各指數(shù)級(jí)的樹(shù)高差異較小，其差異程度隨年齡的增大而增大，10 指數(shù)級(jí)和18 指數(shù)級(jí)的樹(shù)高在第4年時(shí)的差值為0.8 m，到第20年時(shí)差值達(dá)到了6.6 m。在圖3中，立地指數(shù)與直徑的關(guān)系近似成正比，各指數(shù)級(jí)直徑之間的差異較小，10 指數(shù)級(jí)和18 指數(shù)級(jí)的直徑在第4年時(shí)的差值為1 cm，到第20年時(shí)的差值僅為1.2 cm，表3顯示各指數(shù)級(jí)的直徑無(wú)顯著性差異。因此可認(rèn)為造成不同立地指數(shù)的高徑比出現(xiàn)差異的原因是樹(shù)高的生長(zhǎng)的差異，與直徑生長(zhǎng)關(guān)系較小。

圖2 高徑比/樹(shù)高與年齡關(guān)系Fig.2 Relationship between HD/ tree height and age

圖3 直徑與年齡的關(guān)系Fig.3 Relationship between DBH and age

表3 高徑比非參數(shù)檢驗(yàn)分析Table 3 Non-parametric test analysis of height-diameter ratio

2.3 不同林分密度對(duì)高徑比的影響

從圖4可知，同立地指數(shù)同年齡不同林分密度的林分，其高徑比變化的趨勢(shì)較為一致，林分密度與高徑比成正比，且林分株數(shù)密度差距越大高徑比差異越顯著。

如表5所示，C1 和C2 樣地高徑比的顯著性P值為0.010，達(dá)到極顯著差異，C1 和C3 樣地、C2 和C3 樣地的高徑比的P值都小于0.05，也達(dá)到了極顯著水平。高徑比由樹(shù)高和直徑共同決定，并且由圖4可知不同林分密度的直徑差異隨年齡的增大而緩慢增大，方差分析顯示C1 和C2 樣地、C2 和C3 兩種林分密度下直徑的P值分別為0.341和0.387，但C1 和C3 樣地2 種林分密度下直徑的P值達(dá)到了0.076，接近了極顯著水平，這表明林分株數(shù)密度差距越大，直徑差異越顯著。

圖4 高徑比/直徑與年齡的關(guān)系Fig.4 Relationship between HD/DBH and age

表4 不同立地指數(shù)下樹(shù)高、直徑多重比較分析Table 4 Multiple comparative analysis of tree height and DBH under different site indexes

表5 不同林分密度下樹(shù)高、直徑和高徑比多重比較分析Table 5 Multiple comparative analysis of tree height,DBH and HD under different stand densities

從圖5可知，不同林分密度的樹(shù)高差異較小，C1、C2 和C3 樣地的樹(shù)高曲線斜率接近相等，且C1 和C2 樣地樹(shù)高曲線幾乎重合，由表4可知，這3 種林分密度下的樹(shù)高無(wú)顯著性，因此，直徑的差異是造成不同林分密度下高徑比出現(xiàn)差異的主要原因。

圖5 樹(shù)高與年齡的關(guān)系Fig.5 Relationship between tree height and age

2.4 樹(shù)高-直徑的相關(guān)性研究

在圖4中，樹(shù)高與直徑的關(guān)系為較規(guī)律的拋物線，其斜率由大變小，主要是直徑受林分密度效應(yīng)影響較大，在5 a 以前林分未郁閉，樹(shù)高直徑自由生長(zhǎng)，直徑隨年齡變化曲線的斜率較大；5 a以后林分開(kāi)始郁閉，直徑的橫向擴(kuò)張受到抑制，生長(zhǎng)速度減緩，斜率降低。

圖6表明直徑的平方與樹(shù)高、單木斷面積與樹(shù)高近似呈直線的關(guān)系，分別擬合C1、C2 和C3樣地在6 a 郁閉后的平均木解析木的樹(shù)高-直徑模型，如表6所示，所有模型的R2都大于0.99，預(yù)估精度P值也在0.99 左右，這說(shuō)明兩種模型的擬合精度極高。

圖6 直徑的平方/單木斷面積與樹(shù)高關(guān)系Fig.6 Relationship between the square of DBH/ individual-tree basal area and tree height

表6 模型參數(shù)估計(jì)結(jié)果Table 6 Parameter estimation results

3 討 論

3.1 樹(shù)高、直徑、年齡與高徑比的關(guān)系

馬尾松幼齡時(shí)期10 a 以前，首先是樹(shù)高開(kāi)始生長(zhǎng)，然后是直徑生長(zhǎng)[17]，而高徑比是樹(shù)高與直徑的比值，因而在幼林時(shí)期馬尾松的高徑比較大。隨著樹(shù)高生長(zhǎng)速度的減緩和直徑的生長(zhǎng)速度加快，高徑比急劇減小，當(dāng)林分開(kāi)始郁閉一段時(shí)間后，高徑比降至最小值；進(jìn)入中齡林時(shí)期11～20 a，由于林分已經(jīng)郁閉，林分密度對(duì)樹(shù)高生長(zhǎng)有明顯的促進(jìn)作用[18]，此時(shí)期是樹(shù)高生長(zhǎng)的旺盛期，樹(shù)高生長(zhǎng)率稍大于直徑生長(zhǎng)率，因而高徑比小幅度增大。

高徑比在幼齡林時(shí)期與年齡成反比，在中齡林時(shí)期與年齡成正比，與溫佐吾等[2]的結(jié)論一致，而張更新等[5]認(rèn)為高徑比總體上與年齡成反比，忽略了高徑比在中齡林時(shí)期的較大變化，所以高徑比與年齡的關(guān)系應(yīng)分為兩個(gè)階段。直徑和樹(shù)高的異速生長(zhǎng)關(guān)系是造成高徑比變化的根本原因，當(dāng)樹(shù)高和直徑的生長(zhǎng)速度差異隨著年齡的增大而逐漸減小，高徑比的變化幅度也逐漸減小，直至穩(wěn)定。

3.2 不同立地指數(shù)對(duì)高徑比的影響

本研究表明，在相同林分密度、相同年齡且林分郁閉的情況下，立地指數(shù)與高徑成正比，各指數(shù)級(jí)的高徑比之間的差異隨年齡的增大而進(jìn)一步加大。平均生長(zhǎng)量可用來(lái)比較同一樹(shù)種在不同立地條件下生長(zhǎng)的快慢，黃春指出立地指數(shù)越大樹(shù)高生長(zhǎng)量越大，樹(shù)高生長(zhǎng)越快[19]，由于相同林分密度相同年齡不同立地指數(shù)的直徑差異很小，因此立地指數(shù)越大，高徑比越大，這與廖澤釗等[1]的研究結(jié)果一致。王彩云等[3]認(rèn)為高徑比只要在林齡相同的情況下，無(wú)論地位質(zhì)量差異如何，林木的高徑比無(wú)太大差異，李旭等[4]認(rèn)為立地條件對(duì)高徑比的影響很小，張更新等[5]指出除了幼齡，相同的齡組內(nèi)，不論立地條件如何，高徑比無(wú)太大變化。本研究通過(guò)數(shù)據(jù)分析得出相鄰指數(shù)級(jí)間的高徑比差異較小，出現(xiàn)高徑比顯著差異的指數(shù)級(jí)多為12 m 以下的低指數(shù)級(jí)與16 m 以上的高指數(shù)級(jí)之間，指數(shù)級(jí)間相差約4 m 以上的高徑比可到極顯著差異。不同立地指數(shù)下的高徑比出現(xiàn)差異，根本原因是立地指數(shù)對(duì)直徑無(wú)顯著影響，是通過(guò)對(duì)樹(shù)高的顯著影響使高徑比發(fā)生變化。

3.3 不同林分密度對(duì)高徑比的影響

在相同立地指數(shù)相同年齡的情況下，林分密度與高徑比成正比，結(jié)論與王彩云、張更新和廖澤釗等一致[1-5]。林業(yè)界的學(xué)者們?cè)诹址置芏葘?duì)樹(shù)高影響方面做了大量研究，由于林分密度對(duì)樹(shù)高的影響較為復(fù)雜，所得結(jié)論也不盡相同，經(jīng)方差分析，本研究中3種林分密度下的樹(shù)高無(wú)顯著差異，這與諶紅輝、丁貴杰等[20-21]研究表明在一個(gè)相當(dāng)大的中等密度范圍內(nèi)林分密度對(duì)樹(shù)高無(wú)影響的結(jié)論一致。林分密度對(duì)直徑有顯著影響[21-22]，密度與直徑的生長(zhǎng)關(guān)系主要是通過(guò)密度對(duì)冠幅的影響發(fā)生的，密度越大冠幅越小，直徑與冠幅呈正相關(guān)[23]，密度與直徑呈負(fù)相關(guān)。本研究顯示，不同林分密度的樣地間的直徑未達(dá)到極顯著性，但是直徑的顯著性差異隨著林分密度的增大而增大。由于密度效應(yīng)對(duì)直徑有較大影響，而對(duì)樹(shù)高無(wú)顯著影響，當(dāng)樹(shù)高差異很小，林分密度越高，高徑比越大，不同密度間的差異達(dá)到了顯著或極顯著水平[2,24]，因此不同林分密度對(duì)高徑比的影響，其本質(zhì)是直徑的差異造成高徑比的差異。

3.4 樹(shù)高-直徑的相關(guān)性研究

在測(cè)樹(shù)外業(yè)中有很多形式的樹(shù)高-直徑模型可選擇，以提高預(yù)測(cè)樹(shù)高的精度與簡(jiǎn)便性。本研究提出的兩種樹(shù)高-直徑模型表明了在馬尾松人工林中存在著直徑的平方與樹(shù)高、單木斷面積與樹(shù)高呈線性的關(guān)系。雖然樹(shù)高和直徑為異速生長(zhǎng)，其關(guān)系曲線在圖像上呈拋物線形狀，但是單木斷面積與樹(shù)高可呈線性關(guān)系，由于林分總斷面積是由單木斷面積與單位面積上的株數(shù)的乘積得出，當(dāng)總斷面積增加或者減少，可認(rèn)為是單位面積上的林木株數(shù)增加或者減少造成的。本文提出的兩種模型可以較精確地通過(guò)直徑預(yù)測(cè)出樹(shù)高，滿足外業(yè)測(cè)樹(shù)高的需要，根據(jù)此模型編制經(jīng)營(yíng)數(shù)表等還有待進(jìn)一步研究。

4 結(jié) 論

樹(shù)木的樹(shù)高和直徑在不同生長(zhǎng)階段有著不同的生長(zhǎng)狀態(tài)，造成高徑比變化的根本原因是樹(shù)高和直徑的異速生長(zhǎng)關(guān)系，總的來(lái)說(shuō)，在林分郁閉前，高徑比與樹(shù)高、直徑、年齡成反比，林分郁閉后，高徑比與樹(shù)高、直徑、年齡成正比。

本研究中立地指數(shù)對(duì)直徑影響較小，指數(shù)級(jí)間差距越大，樹(shù)高差異越明顯；在林分郁閉之后，立地指數(shù)越大的林分，其林木的高徑比越大，相鄰指數(shù)級(jí)間的高徑比差異較小，指數(shù)級(jí)間相差4 m以上的高徑比達(dá)到極顯著差異；不同指數(shù)級(jí)的高徑比變化程度主要由樹(shù)高生長(zhǎng)差異決定。

林分密度對(duì)樹(shù)高影響較小，林分株數(shù)密度差距越大，其直徑的差異越顯著；林分株數(shù)密度與高徑比成正比，在不同林分株數(shù)密度下高徑比出現(xiàn)差異的本質(zhì)是直徑生長(zhǎng)的差異。

直徑的平方與樹(shù)高的關(guān)系實(shí)質(zhì)是單木斷面積與樹(shù)高的關(guān)系，本文所提出的樹(shù)高-直徑關(guān)系模型成立，在圖像上為一條直線，且在林分郁閉前后，其斜率無(wú)變化。

本研究明晰了馬尾松人工林樹(shù)高、直徑和年齡與高徑比的關(guān)系，并分別探討了不同立地指數(shù)、林分密度對(duì)高徑比的影響，以及樹(shù)高-直徑關(guān)系模型，對(duì)研究馬尾松人工林的樹(shù)高和直徑的相關(guān)系數(shù)具有現(xiàn)實(shí)意義，為馬尾松人工林培育各類徑材提供理論指導(dǎo)，從而達(dá)到最佳的經(jīng)營(yíng)效果。提出的樹(shù)高-直徑關(guān)系模型可以應(yīng)用于外業(yè)測(cè)樹(shù)高，可提高樹(shù)高的預(yù)估精度。