姜 麗，張雄清，段愛國，張建國

（中國林業(yè)科學(xué)研究院林業(yè)研究所，國家林業(yè)和草原局林木培育重點實驗室，北京 100091）

林分密度指標(biāo)是評定一個時期林分生長密度的尺度指標(biāo)，而確定合適的指標(biāo)是開展林分密度研究的重要前提。Dean[1]認(rèn)為林分密度指標(biāo)則至少需要滿足以下3個條件：（1）與林分生長有關(guān)；（2）與立地質(zhì)量、林分年齡等無內(nèi)在聯(lián)系；（3）較為容易測量得到。林分密度指標(biāo)有多種，如：Reineke林分密度指數(shù)（SDI）[2]，Nilson密度指數(shù) （SD）[3-4]，優(yōu)勢高-營養(yǎng)面積比[5]，相對密度RD指數(shù)[6]，以及相對植距RS指數(shù)[7]等。這些密度指標(biāo)的提出均是處于一定的條件下，每個指標(biāo)都有其使用條件及適用范圍。選用不同的林分密度指標(biāo)會直接影響著所建模型的精度與預(yù)估結(jié)果，因此選取合適的密度指標(biāo)建立單木生長模型顯得尤為重要。

基于林木長期生長數(shù)據(jù)，預(yù)估單木年生長量時，之前常用的固定生長率法是假設(shè)在森林的整個生長期內(nèi)，林木的生長率是固定不變的這一假設(shè)進(jìn)行的研究。顯而易見，這種假設(shè)并不符合現(xiàn)實規(guī)律，隨著林分生長時間的推移，所用于預(yù)估模型中的特征因子（林分優(yōu)勢高、林齡、林分密度）和林分中單木的特征因素（樹高、胸徑）等在每一個生長階段都會有變化，必然也會導(dǎo)致單木生長量的變化。為克服固定生長率法的這個缺陷，不少方法被陸續(xù)提出，如內(nèi)插法和迭代法等。Cao等[8]基于迭代法提出的可變生長率法。該方法將林分生長期間林分與單木因子變化引起的單木年生長量變化納入分析，并且用時比迭代法更少。張雄清等[9]比較了可變生長率法與固定生長率法，發(fā)現(xiàn)可變生長率法表現(xiàn)更好。

杉木（Cunninghamia lanceolata （Lamb.）Hook.）是中國南部地區(qū)最重要的用材林樹種。本研究基于可變生長率法構(gòu)建包含不同林分密度指標(biāo)的杉木單木直徑年生長模型，選出適宜單木直徑年生長模型的最優(yōu)林分密度指標(biāo)，以期實現(xiàn)杉木人工林單木直徑生長的精準(zhǔn)預(yù)測，為實現(xiàn)杉木人工林質(zhì)量精準(zhǔn)提升提供密度管理支撐依據(jù)。

1 研究數(shù)據(jù)與方法

1.1 研究區(qū)域及數(shù)據(jù)

樣地設(shè)置在福建武夷山北部的邵武市，位于117°43′ E，27°05′ N。地貌特征主要是高山和低丘陵地區(qū)，海拔為250～700 m，坡度為25°～35°，為亞熱帶季風(fēng)氣候，年平均溫度17.7 ℃，1月平均溫度6.8 ℃，7月平均溫度28 ℃，最低極端氣溫為?7.9 ℃，年日照時間1 740.7 h，平均霜凍期為95 d。年降水量1 768 mm，年平均相對濕度為82%。林下植被有中華杜英（Elaeocarpus chinensis （Gardn.et Champ.） Hook.F.ex Benth.）、胡頹子（Elaeagnus pungens Thunb.）、狗脊（Woodwardia japonica （L.f.） Sm.）、芒萁（Dicranopteris pedata （Houtt.） Nakaike）、烏毛蕨（Blechnum orientale L.）和扇葉鐵線蕨（Adiantum flabellulatum L.）等。

試驗林使用1年生苗木于1982年造林，完全隨機區(qū)組設(shè)計，分5種造林密度：A：2 m×3 m（1 667 株·hm?2），B：2 m×1.5 m（3 333株·hm?2），C：2 m×1 m（5 000 株·hm?2），D：1 m×1.5 m（6 667 株·hm?2）和E：1 m×1 m（10 000 株·hm?2）。每個樣地大小為 20 m×30 m，每種造林密度均重復(fù)3次，總計15個樣地。標(biāo)記了一共4 800棵樹，在冬季測量樹高和胸徑，從1984年至1990年，每年進(jìn)行1次測量；從1992年至2010年，每隔1或者2年進(jìn)行1次測量。具體統(tǒng)計數(shù)據(jù)見表1。研究數(shù)據(jù)隨機抽取60%用于模型構(gòu)建，剩余40%用于模型驗證。

表1 杉木人工林林分和單木變量統(tǒng)計Table 1 Summary statistics of stand and tree variables of Chinese fir plantation

1.2 研究方法

利用常用的7種林分密度指標(biāo)，分別是每公頃株數(shù)密度N、每公頃胸高斷面積Ba、林分密度SDI指數(shù)、SD指數(shù)、優(yōu)勢高—營養(yǎng)面積比Z指數(shù)、相對密度RD指數(shù)和相對植距RS指數(shù)構(gòu)建杉木單木直徑生長模型。

1.2.1 林分密度指標(biāo)

（1） SDI指數(shù)

SDI指數(shù)[2]是通過單位面積上林木的株數(shù)與林木的平均胸徑的關(guān)系確定的，是反映林分直徑分布與林分內(nèi)單位面積的株數(shù)之間關(guān)系的綜合指標(biāo)。因此，SDI密度指數(shù)在林木密度管理中使用較常見，其數(shù)學(xué)表達(dá)式為：

式中：Dq為林分平方平均胸徑cm；D0—林分標(biāo)準(zhǔn)直徑，杉木一般取值20 cm；β為Reinek自稀疏系數(shù)，一般取值?1.605。

（2） SD指數(shù)

根據(jù)Nilson[3-4]對林木之間的平均距離L和胸徑Dq的定義，林分相對密度SD指數(shù)的公式如下[10]：

式中：a，b為參數(shù)，k=?a，D0—取20 cm。

（3）優(yōu)勢高—營養(yǎng)面積比Z指數(shù)

劉金福等[5]采用單位面積內(nèi)株數(shù)與優(yōu)勢木的平均高來反映林分中林木的相對密度的大小。其表達(dá)式如下：

式中：N—每公頃株數(shù)；Hd—優(yōu)勢木平均高。

（4） 相對密度RD指數(shù)

Curtis[6]提出利用林分?jǐn)嗝娣e和平方平均胸徑計算林分密度指標(biāo)，該指標(biāo)簡單，容易計算[11]，其公式如下：

（5）相對植距RS指數(shù)

相對植距是跟立地、胸徑等因子無關(guān)的林分密度指標(biāo)，表達(dá)式如下[7]：

1.2.2 單木直徑年生長模型 構(gòu)建單木直徑生長模型，可以引入多個自變量，其中包括單木特征因子（樹高、胸徑等）和林分特征因子（林分密度、林齡、林分優(yōu)勢高等）。一般這些因子歸為3類：林齡、競爭因子和立地因子。而立地因子中，一般以林分優(yōu)勢高表示。因此，在單木的生長過程中，主要是受到競爭、立地和林齡的影響。在本研究中，選取了林齡（A）、7種林分密度指標(biāo)和林分優(yōu)勢高（Hd）構(gòu)建杉木單木直徑生長模型?；诳勺兩L率法，利用遞歸的方式推導(dǎo)過程如下[9]：

（t+1）年時:

（t+2）年時:

（t+q）年時:

式子中：Di,t為t年時第 i株林木的直徑（cm）；

At為t年時林分的平均年齡；

K為7種林分密度指標(biāo)：包括林分胸高斷面積（BA）、優(yōu)勢高營養(yǎng)面積比（Z指數(shù)）、株數(shù)（N）、相對植距（RS）、林分密度指數(shù)（SDI）、相對直徑（RD）和SD指數(shù)；

Hdt為t年時的林分優(yōu)勢高；

α1～α4為模型參數(shù)；

將林分變量因子，如林分優(yōu)勢高等每年的變化納入分析，通過建立模型預(yù)估對應(yīng)的林分變量，再將預(yù)估出來的林分優(yōu)勢高用來建立杉木單木直徑生長模型。林分優(yōu)勢高模型形勢如下：

其中β1、β2為待估參數(shù)。

本研究中，杉木單木直徑模型的參數(shù)估計均利用SAS中非線性回歸模塊來完成。利用可變生長率法估計模型參數(shù)時，應(yīng)用了循環(huán)運算，直到循環(huán)停止。

1.2.3 模型選擇評價 林分優(yōu)勢高模型和單木直徑生長模型通過統(tǒng)計量平均絕對偏差（MAD）、均方根誤差（RMSE）和決定系數(shù)（R2） 進(jìn)行評價。

2 結(jié)果與分析

杉木林分優(yōu)勢高模型的參數(shù)估計結(jié)果及模型的均方根誤差、決定系數(shù)見表2。從模型結(jié)果中參數(shù)的標(biāo)準(zhǔn)誤差可以得知，模型的參數(shù)估計值均為有效值，并且模型的決定系數(shù)很高，達(dá)到了0.953 9，均方根誤差RMSE較小，為1.343 1。因此林分優(yōu)勢高模型可用來預(yù)估杉木林分優(yōu)勢高。

表2 杉木林分優(yōu)勢高模型的參數(shù)估計、標(biāo)準(zhǔn)誤差、決定系數(shù)及均方根誤差Table 2 Parameter estimation,standard error,R2 and RMSE of dominant height model of Chinese fir.

本研究分別用7種不同林分密度指標(biāo)，利用可變生長率法建立了7種杉木單木直徑生長模型，得到的參數(shù)估計和標(biāo)準(zhǔn)誤差及模型評價見表3。由表3可知，大多數(shù)模型參數(shù)估計值的標(biāo)準(zhǔn)誤差都在接受范圍之內(nèi)，因此模型參數(shù)估計值都是有效的。從決定系數(shù)R2可以得到模型的擬合效果都相當(dāng)高，均在0.965以上。7類將密度指標(biāo)納入分析的杉木單木直徑模型精確度，大部分都高于不含密度指標(biāo)模型。只有以相對直徑RD為密度指標(biāo)的模型，其決定系數(shù)、均方根誤差和平均絕對偏差與不含密度指標(biāo)的模型一致。其次，從模型參數(shù)估計α3也可以發(fā)現(xiàn)，其標(biāo)準(zhǔn)誤比較高，參數(shù)估計在0.05水平不顯著，相對直徑參數(shù)估計不穩(wěn)定。因此可以得到在杉木直徑模型中，相對直徑RD并不適合作為模型的密度指標(biāo)。

表3 不同林分密度指標(biāo)應(yīng)用于杉木單木直徑生長模型的參數(shù)估計及模型評價Table 3 Parameter estimation and model evaluation of annual tree diameter growth model of Chinese fir with different density indices

所有含有密度指標(biāo)的模型中，模型決定系數(shù)最高的為以林分密度指數(shù)（SDI）為密度指標(biāo)的模型，為0.971 8，同時也有最低的平均絕對偏差以及均方根誤差，分別為0.511 3、0.836 1。模型精度其次為以林分?jǐn)嗝娣eBA和Nilson密度指數(shù)為密度指標(biāo)的模型，R2分別為0.971 7、0.971 6。其余5種林分密度指標(biāo)表現(xiàn)介于表現(xiàn)最佳SDI指標(biāo)（0.971 8）和最差RD指標(biāo)（0.967 8）之間，所得到的模型精度從大到小依次為林分?jǐn)嗝娣eBA（0.971 7）、Nilson指數(shù)SD（0.971 6）、優(yōu)勢高營養(yǎng)面積比Z（0.970 6）、每公頃株數(shù)N（0.970 5）和相對植距RS（0.970 4）。

基于最優(yōu)的含SDI指數(shù)建立的單木直徑生長模型，繪制了模型殘差圖。從圖1中可以看到，模型在直徑較小的時候，如0＜D≤5 cm，傾向于低估杉木單木直徑生長；在5＜D≤20 cm的時候，模型傾向于高估單木直徑生長；在D＞20 cm的時候，模型則更傾向于低估杉木單木直徑生長。總的來說殘差均勻分布在0附近，服從正態(tài)分布，總體模型情況較好。

圖1 包含SDI指數(shù)的單木直徑生長模型的殘差分布Fig.1 Residual distributions of annual tree diameter growth model including SDI

從所有模型的參數(shù)估計方面來看，α2在所有模型中均為負(fù)值，且α2與林齡A相關(guān)，可以看出杉木單木直徑增長與年齡在任何密度指數(shù)下均呈現(xiàn)負(fù)相關(guān)關(guān)系，即隨著林分年齡的增大，單木直徑生長量減小。α3在所有模型中均為負(fù)值，且是林分密度指標(biāo)因子的系數(shù)，可以得到杉木單木直徑生長隨著密度指標(biāo)的增大而減小。α4在所有模型中得到立地估計均為負(fù)值，且與林分優(yōu)勢高Hd的倒數(shù)負(fù)相關(guān)，因此隨著林分優(yōu)勢高的增大，杉木單木直徑生長量增大。林分優(yōu)勢高通?？梢源砹⒌刭|(zhì)量，即林分優(yōu)勢高越大立地質(zhì)量越好，杉木直徑生長量也就越大。

3 討論

一般情況在預(yù)估單木的年生長量時，主要是基于樹干解析中年輪寬度分析[11]。但是利用樹干解析分析單木年生長量時需要伐倒樹木，而且所選擇的解析木一般是選擇平均木，并不能完全代表整個林分的生長狀況，其次樹干解析需要耗費大量的人力物力[12]。近年來，基于已有的連續(xù)調(diào)查數(shù)據(jù)，學(xué)者研究和利用單木生長模型來預(yù)估單木年生長量。例如，為了獲得單木年生長量（如胸徑、樹高和材積），通常用的方法是固定生長率法，其前提是假設(shè)在整個生長期，各單木的年生長量是保持不變的。一般常用定期平均生長量來代替年生長量[7]。很明顯該估計方法比較簡單，并不符合林木的實際生長規(guī)律。因為隨著森林的演替過程，一些林分的特征（立地，年齡，林分競爭等）和單木本身的一些因子（胸徑，樹高，冠幅等）都會隨著林齡變化而變化，因此也必然會導(dǎo)致單木每年的生長量也發(fā)生變化。為解決固定生長率法的缺點，Mcdill和Amateis[13]提出利用內(nèi)插法來估計單木年生長。之后Cao等[8]也基于內(nèi)插法模擬單木直徑、樹高和樹冠比的變化，發(fā)現(xiàn)內(nèi)插法比固定生長率法的模擬效果要好的多。然而，內(nèi)插法的前提是在單木年存活率不變。為進(jìn)一步提高模型的預(yù)測精度，Cao[14]提出了迭代法。該方法比固定生長率法和內(nèi)插法均表現(xiàn)出明顯的優(yōu)越性，因為該方法考慮了單木在生長期間，林分因子（如競爭和林分優(yōu)勢高）及單木本身的變化而引起的單木直徑年生長量的變化。Cao[15]基于迭代法提出了可變生長率法，該方法也是考慮了林木在生長期間，林分變量如競爭和林分優(yōu)勢高及單木本身的變化而引起的單木胸徑年生長量的變化。而且可變生長率法在估計模型參數(shù)時，比迭代法簡單，計算耗時較少[16]。張雄清等[9]利用可變生長率法構(gòu)建了油松單木直徑和存活率的預(yù)估模型，模型擬合效果很好。之后張雄清等[12]基于可變生長率法構(gòu)建了單木油松年生生長模型，發(fā)現(xiàn)該方法不僅能夠很好地預(yù)測單木斷面積生長，而且很好地解決了調(diào)查間隔期不一致而導(dǎo)致的相容性情況。

本研究發(fā)現(xiàn)密度指標(biāo)的估計值為負(fù)數(shù)，說明單木直徑生長隨著密度的增加而減少。杉木單木的生長主要依靠良好的光照，充沛的水資源以及土壤營養(yǎng)物等，在林分密度增大的情況下，必然會導(dǎo)致樹木間對光照、水資源、營養(yǎng)物質(zhì)等的競爭，競爭的后果則會導(dǎo)致單木直徑生長量的減少。即隨著競爭強度的增大，樹木的生長量會減少，這與之前學(xué)者們的結(jié)論是一致的[17]。此外表明，用林分密度指數(shù)SDI建模預(yù)測杉木單木直徑生長的模型是最佳選擇，這與高啟東等[18]研究結(jié)果一致。許多學(xué)者認(rèn)為SDI指數(shù)的主要優(yōu)點就是它整合了林分的株樹密度和平均直徑[19-20]。Zeide[21]報道稱在純林中，SDI指數(shù)是最簡單而且描述林分密度最有效的一個指標(biāo)。正是由于SDI指數(shù)的優(yōu)點，多數(shù)學(xué)者把SDI指數(shù)引入生長模型中來表示林分密競爭指標(biāo)[22-24]。也有學(xué)者在生長模型[25]和過程模型[26]中引入了SDI指數(shù)并得到了很好的擬合結(jié)果。此外，也有學(xué)者認(rèn)為Nilson密度指數(shù)（SD）是最優(yōu)林分密度指標(biāo)，如車少輝等[27]認(rèn)為林分密度指數(shù)SDI可以與相對植距RS相互轉(zhuǎn)換。模型表現(xiàn)最差的為相對直徑RD模型和不含密度指標(biāo)模型。因此，在構(gòu)建杉木單木直徑生長模型，將相對直徑RD以外的林分密度指標(biāo)納入模型是非常必要的，可顯著提高模型估計的精確度。

4 結(jié)論

本研究基于可變生長率法，應(yīng)用7種不同的林分密度指標(biāo)以及不含有密度指標(biāo)，構(gòu)建得到了8種杉木單木直徑年生長模型，模型的擬合精度比較高，R2均在0.96以上。相比較含有密度指標(biāo)的杉木單木直徑模型精度，高于不含密度指標(biāo)模型。在所有含有密度指標(biāo)的模型中，決定系數(shù)最高的是以林分密度指數(shù)（SDI）為密度指標(biāo)的模型，其次為以林分?jǐn)嗝娣eBA和Nilson密度指數(shù)為密度指標(biāo)的模型。因此，在估計杉木直徑年生長量模型時，應(yīng)該考慮引入Reinek提出的 SDI密度指數(shù)，以實現(xiàn)杉木人工林生長的精準(zhǔn)預(yù)測。