側(cè)柏人工幼齡林單木生長(zhǎng)的競(jìng)爭(zhēng)效應(yīng)模型1)

馬士友 楊華 張恒 尚旭陽(yáng) 張?jiān)?張永福

(省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室(北京林業(yè)大學(xué))，北京，100083) (北京市十三陵林場(chǎng))

?

側(cè)柏人工幼齡林單木生長(zhǎng)的競(jìng)爭(zhēng)效應(yīng)模型1)

馬士友 楊華 張恒 尚旭陽(yáng) 張?jiān)?張永福

(省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室(北京林業(yè)大學(xué))，北京，100083) (北京市十三陵林場(chǎng))

為量化采伐對(duì)林木生長(zhǎng)的影響，以北京郊區(qū)側(cè)柏人工幼齡林為研究對(duì)象，以競(jìng)爭(zhēng)指數(shù)、采伐指數(shù)作為自變量，伐后5 a胸高斷面積生長(zhǎng)量為因變量，構(gòu)建側(cè)柏人工幼齡林競(jìng)爭(zhēng)效應(yīng)模型，采用判斷系數(shù)(R2)、剩余標(biāo)準(zhǔn)差(RMSE)、赤池信息準(zhǔn)則(AIC)對(duì)模型進(jìn)行評(píng)價(jià)。結(jié)果表明：斷面積比競(jìng)爭(zhēng)指數(shù)(CI2)、采伐指數(shù)(RI2)優(yōu)于其余4種指數(shù)；得到了以5種競(jìng)爭(zhēng)指數(shù)和采伐指數(shù)為自變量的競(jìng)爭(zhēng)效應(yīng)模型，根據(jù)模型R2、RMSE、AIC值，推斷出側(cè)柏人工幼齡林的最佳競(jìng)爭(zhēng)范圍的半徑為4 m；以CI2r4、RI2r4為自變量的模型比只有CI2r4的模型R2提高了33.8%，說(shuō)明采伐對(duì)側(cè)柏單木生長(zhǎng)的影響顯著。

側(cè)柏人工幼齡林；競(jìng)爭(zhēng)范圍；采伐；競(jìng)爭(zhēng)指數(shù)；采伐指數(shù)

Platycladusorientalisyoung plantation were selected to quantify the cutting influence on tree growth in Beijing suburb. We used the competition index, cutting index and 5-year basal area growth after cutting as dependent variables to constructP.orientalisyoung plantation competition effect model. The model was evaluated by the judgment coefficient (R2), the residual standard deviation (RMSE) and the Akaike Information Criterion (AIC). The basal area ratio competition indexCI2, cutting indexRI2were superior to the other four indexes. Then the competition effect model with five kinds of competition index and cutting index were got as independent variables. According to theR2,RMSE,AIC, the best competition range of semi diameter ofP.orientalisyoung plantation was 4 m. TheR2of the model, which usedCI2r4andRI2r4as the independent variables, was increased by 33.8%, indicating that cutting has a significant influence on individual tree growth.

競(jìng)爭(zhēng)是生物間相互作用的一個(gè)重要方面，是指2個(gè)或多個(gè)生物體在對(duì)同一環(huán)境資源和能量的爭(zhēng)奪中所發(fā)生的相互作用。由于林分內(nèi)存在大小不同的樹(shù)木，它們所占據(jù)的生長(zhǎng)空間不同，各自承受著不同的競(jìng)爭(zhēng)壓力，這種壓力必然導(dǎo)致林木個(gè)體生長(zhǎng)的差異性。為描述這種單木的生長(zhǎng)動(dòng)態(tài)，引入了單木競(jìng)爭(zhēng)指數(shù)[1]，為了更好的描述不同林分的競(jìng)爭(zhēng)狀態(tài)，建立了各種各樣的競(jìng)爭(zhēng)指數(shù)模型[2-12]。

對(duì)象木最佳競(jìng)爭(zhēng)范圍指對(duì)某一株對(duì)象木而言，周?chē)徑灸軌驅(qū)ζ洚a(chǎn)生影響的空間范圍。在這空間范圍以外的樹(shù)木對(duì)對(duì)象木生長(zhǎng)的影響很小，可以忽略不計(jì)[13]。根據(jù)研究對(duì)象的不同，對(duì)最佳競(jìng)爭(zhēng)范圍的研究可分為2種。1)以對(duì)象木表現(xiàn)為研究對(duì)象，以對(duì)象木的表現(xiàn)作為因變量、競(jìng)爭(zhēng)木的相關(guān)信息作為自變量，建立回歸方程，通過(guò)不斷地增加更遠(yuǎn)距離競(jìng)爭(zhēng)木的信息，用回歸方程取得最大R2值或最小剩余誤差的方法來(lái)確定與對(duì)象木距離最遠(yuǎn)的有效競(jìng)爭(zhēng)木，間接地估計(jì)林分中單木的最佳競(jìng)爭(zhēng)范圍[4，11，13-14]。2)以鄰近木的表現(xiàn)為研究對(duì)象，林木采伐后產(chǎn)生林窗，鄰近木會(huì)相應(yīng)的產(chǎn)生生長(zhǎng)釋放，楊光等[15]、夏冰等[16]通過(guò)判定伐樁周?chē)a(chǎn)生生長(zhǎng)釋放的最遠(yuǎn)競(jìng)爭(zhēng)木與伐樁所對(duì)應(yīng)的距離來(lái)判定最佳競(jìng)爭(zhēng)范圍。

側(cè)柏是北京的鄉(xiāng)土樹(shù)種，分布面積最廣，是北京山區(qū)低山針葉林的主要樹(shù)種，北京地區(qū)側(cè)柏林以中幼齡林為主，需要進(jìn)行撫育。研究采伐對(duì)側(cè)柏幼齡林生長(zhǎng)的影響，對(duì)探索北京地區(qū)側(cè)柏人工林合理的經(jīng)營(yíng)模式，促進(jìn)森林生態(tài)系統(tǒng)功能的完善和提高有著重要意義。

前人對(duì)于林木采伐的相關(guān)研究中，既有在林分水平上，研究擇伐后，林分生長(zhǎng)和空間結(jié)構(gòu)的變化[17-18]，也有在單木水平上，利用樹(shù)木年代學(xué)方法，分析采伐后，單木生長(zhǎng)的變化[15-16,19]。楊光等[15]、夏冰等[16]以伐樁為中心，研究鄰近木的生長(zhǎng)變化。從文獻(xiàn)研究來(lái)看，因?yàn)榱址植煞ツ甏眠h(yuǎn)，無(wú)法準(zhǔn)確獲取伐樁的大小和位置信息，采伐強(qiáng)度只是粗略地估計(jì)，缺少對(duì)象木周?chē)煞ツ厩闆r的量化指標(biāo)，本文將圍繞這一問(wèn)題，量化對(duì)象木周?chē)?，分析采伐?duì)對(duì)象木生長(zhǎng)產(chǎn)生的影響。

1 研究區(qū)概況

研究區(qū)位于北京市西北郊昌平區(qū)的十三陵林場(chǎng)，該地屬于燕山系低山丘陵區(qū)，地理坐標(biāo)為東經(jīng)115°50′17″～116°29′49″，北緯40°2′18″～40°23′13″，山地海拔68.0～954.2 m，平均海拔400 m，林場(chǎng)經(jīng)營(yíng)面積8 561.5 hm2，林地面積8 553.8 hm2，有林地面積6 926.5 hm2，側(cè)柏面積占有林地總面積的62.7%。研究區(qū)屬季風(fēng)氣候區(qū)，全年平均氣溫11.6 ℃，1月份平均氣溫-4.1 ℃，7月份平均氣溫25.8 ℃，全年≥0 ℃積溫為4 500 ℃，≥10 ℃積溫為4 200 ℃。年平均日照時(shí)間2 669 h，無(wú)霜期為202 d，平均生長(zhǎng)期為200 d。年平均降水量584 mm，6—8月份降水量占全年的75%以上。山地大部分巖石裸露，少量風(fēng)化土層的厚度一般在20～40 cm，含石礫量達(dá)40%以上。土壤pH值呈堿性或中性反應(yīng)的碳酸鹽褐色土壤，其水分條件差，肥力低，保墑能力弱。山地坡向以陽(yáng)坡、半陽(yáng)坡為主，坡度一般在30°～45°。

2 研究方法

2.1 數(shù)據(jù)獲取

在研究區(qū)內(nèi)設(shè)置密度為1 000～2 000株/hm2的人工側(cè)柏幼齡林樣地共5塊，樣地在2010年進(jìn)行了不同強(qiáng)度的撫育間伐，采伐強(qiáng)度在20%～30%，側(cè)柏林年齡在30～40 a，樣地面積0.03～0.072 hm2。對(duì)樣地內(nèi)的所有喬木(胸徑>5 cm)進(jìn)行編號(hào)，記錄其坐標(biāo)值，測(cè)定每木的胸徑、樹(shù)高、第一活枝高、東西向和南北向2個(gè)方向的冠幅、伐樁位置及根徑。在樣地邊界緩沖區(qū)5 m以?xún)?nèi)，隨機(jī)選取側(cè)柏對(duì)象木20株進(jìn)行解析，以1 m區(qū)分段截取圓盤(pán)，帶回實(shí)驗(yàn)室進(jìn)行內(nèi)業(yè)測(cè)定。用砂紙對(duì)圓盤(pán)進(jìn)行打磨，直到年輪界限清晰為止，對(duì)不清楚的年輪利用顯微鏡進(jìn)行校正。利用LinTab6系統(tǒng)測(cè)量1.3 m處圓盤(pán)東西南北4個(gè)方向的年輪寬度，精確到0.01 mm，取4個(gè)方向年輪寬度的平均值作為當(dāng)年徑向生長(zhǎng)量，用國(guó)際年輪庫(kù)的COFECHA[20]交叉定年質(zhì)量控制程序進(jìn)行交叉定年檢驗(yàn)，直到滿(mǎn)足要求為止。以采伐當(dāng)年2010年到2014年作為5 a間隔期，計(jì)算5 a間隔期胸高斷面積生長(zhǎng)量。

木材干縮率計(jì)算：林木在采伐后失水，會(huì)有一定程度的收縮。待解析木圓盤(pán)完全干燥后測(cè)量年輪寬度，用伐前胸徑值除以風(fēng)干后胸徑值得到側(cè)柏的平均干縮率為1.047。將測(cè)得的圓盤(pán)年輪寬度乘以1.047得到伐前年輪寬度。

2.2 相關(guān)競(jìng)爭(zhēng)指數(shù)的選定

2.2.1 確定側(cè)柏最佳競(jìng)爭(zhēng)范圍

以對(duì)象木為圓心，以1 m為步長(zhǎng)，競(jìng)爭(zhēng)半徑從1～5 m逐步擴(kuò)大，共定義5個(gè)范圍，分別用r1、r2、r3、r4、r5表示，計(jì)算各競(jìng)爭(zhēng)范圍內(nèi)的競(jìng)爭(zhēng)指數(shù)、采伐指數(shù)(圖1)。

圖1 對(duì)象木所受影響范圍示意圖

2.2.2 競(jìng)爭(zhēng)指數(shù)的選擇

本文采取與胸徑相關(guān)的5個(gè)競(jìng)爭(zhēng)指數(shù)，計(jì)算對(duì)象木不同競(jìng)爭(zhēng)范圍下的競(jìng)爭(zhēng)強(qiáng)度。5個(gè)競(jìng)爭(zhēng)指數(shù)計(jì)算公式如下：

(1)

(2)

(3)

(4)

(5)

式中:CI為競(jìng)爭(zhēng)指數(shù)；Di表示競(jìng)爭(zhēng)范圍內(nèi)第i株競(jìng)爭(zhēng)木胸徑；D為對(duì)象木胸徑；Li為第i株競(jìng)爭(zhēng)木與對(duì)象木的距離。

CI3為Hegyi(1974)年提出的競(jìng)爭(zhēng)指數(shù)，該指數(shù)易于計(jì)算應(yīng)用廣泛[2，7，11，14]，CI1、CI2、CI4是以Hegyi競(jìng)爭(zhēng)指數(shù)為基礎(chǔ)，進(jìn)行相應(yīng)的變形，而CI5為張躍西[21]提出的競(jìng)爭(zhēng)指數(shù)模型，該模型根據(jù)生態(tài)學(xué)原理及自疏規(guī)律而提出[4]。

2.2.3 采伐指數(shù)

為衡量對(duì)象木周?chē)煞?duì)對(duì)象木生長(zhǎng)的影響，本文在前人單木競(jìng)爭(zhēng)指數(shù)模型的基礎(chǔ)上，提出采伐指數(shù)，根據(jù)公式(1)～(5)進(jìn)行計(jì)算，其中Di表示2010年采伐后伐樁的根徑，Li表示第i株伐樁與對(duì)象木的距離，得到的采伐指數(shù)分別用RI1、RI2、RI3、RI4、RI5表示。

2.3 競(jìng)爭(zhēng)效應(yīng)模型

考慮到競(jìng)爭(zhēng)對(duì)生長(zhǎng)的影響，在無(wú)競(jìng)爭(zhēng)狀態(tài)下的樹(shù)木生長(zhǎng)量為R0，當(dāng)競(jìng)爭(zhēng)強(qiáng)度加大時(shí)，由競(jìng)爭(zhēng)指數(shù)進(jìn)行修正。本文假設(shè)由于采伐造成對(duì)象木競(jìng)爭(zhēng)強(qiáng)度迅速降低，會(huì)引起對(duì)象木生長(zhǎng)超過(guò)本身競(jìng)爭(zhēng)強(qiáng)度對(duì)應(yīng)的生長(zhǎng)量，因此，引入采伐指數(shù)對(duì)模型進(jìn)行修正。林木間的競(jìng)爭(zhēng)影響具有可加性[14]，本研究選取Weiner[22]提出的原模型的基礎(chǔ)上引入采伐指數(shù)RI，作為側(cè)柏生長(zhǎng)的競(jìng)爭(zhēng)效應(yīng)模型(6)。

(6)

為了便于回歸分析，將原模型取倒數(shù)轉(zhuǎn)化為線性方程，用ln(BAI+1)代替生長(zhǎng)量G，得到公式(7)：

(7)

式中：G為生長(zhǎng)量，BAI為對(duì)象木伐后5年間隔期胸高斷面積生長(zhǎng)量，R0為待估參數(shù)，CI為競(jìng)爭(zhēng)范圍內(nèi)競(jìng)爭(zhēng)指數(shù)，RI為競(jìng)爭(zhēng)范圍內(nèi)采伐指數(shù)，CIri為不同競(jìng)爭(zhēng)范圍內(nèi)的競(jìng)爭(zhēng)指數(shù)，RIri為不同競(jìng)爭(zhēng)范圍內(nèi)的采伐指數(shù)。

2.4 模型的評(píng)價(jià)與檢驗(yàn)

對(duì)建模樣本進(jìn)行檢驗(yàn)，計(jì)算檢驗(yàn)指標(biāo)，評(píng)價(jià)模型優(yōu)度。用于模型評(píng)價(jià)與檢驗(yàn)的指標(biāo)主要有R2(判斷系數(shù))、RMSE(剩余標(biāo)準(zhǔn)差)、AIC(赤池信息準(zhǔn)則)。R2越接近1越好；RMSE越接近0越好；AIC越小越好。

(8)

(9)

(10)

3 結(jié)果與分析

3.1 設(shè)定側(cè)柏對(duì)象木競(jìng)爭(zhēng)范圍

對(duì)樣地調(diào)查數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，得到樹(shù)冠的平均半徑為2.7 m，2株樹(shù)木之間的平均距離為2.3 m。有研究表明[23]，競(jìng)爭(zhēng)范圍一般取對(duì)象木樹(shù)高的一半，本文側(cè)柏對(duì)象木平均樹(shù)高為8.0 m，所以將競(jìng)爭(zhēng)范圍限定在5 m之內(nèi)，分別以1、2、3、4、5 m作為對(duì)象木影響圈半徑。

3.2 側(cè)柏最佳競(jìng)爭(zhēng)范圍

使用公式(1)～(5)，分別計(jì)算樣地內(nèi)隨機(jī)抽取的對(duì)象木競(jìng)爭(zhēng)范圍內(nèi)，競(jìng)爭(zhēng)指數(shù)、采伐指數(shù)；將ln(BAI+1)作為因變量，競(jìng)爭(zhēng)指數(shù)、采伐指數(shù)作為自變量，通過(guò)多元線性回歸分析建立人工側(cè)柏幼齡林的個(gè)體生長(zhǎng)競(jìng)爭(zhēng)效應(yīng)模型?；貧w結(jié)果見(jiàn)表1。

通過(guò)表1可以看出，隨著影響半徑的增大，5種模型的判斷系數(shù)R2逐漸增大，剩余標(biāo)準(zhǔn)差RMSE和AIC逐漸減小。采用CI3RI3、CI5RI5作為自變量時(shí)，影響半徑為5 m時(shí)，判斷系數(shù)R2最大，RMSE、AIC最小，而與影響半徑為4 m時(shí)相比差別不大；采用CI1RI1、CI2RI2、CI4RI4時(shí)，在各表中影響半徑為4 m時(shí)，判斷系數(shù)R2最大，RMSE、AIC最小。

綜合上述考慮，認(rèn)為側(cè)柏人工幼齡林的最佳競(jìng)爭(zhēng)范圍在4 m左右。影響半徑為4 m時(shí)，采用斷面積比競(jìng)爭(zhēng)指數(shù)CI2與采伐指數(shù)RI2作為自變量，模型的檢驗(yàn)指標(biāo)R2(0.637)最大，RMSE、AIC值最小，擬合方程效果最好。

表1 不同影響半徑下5種競(jìng)爭(zhēng)效應(yīng)模型擬合優(yōu)度對(duì)比

影響半徑/m擬合優(yōu)度(CI2RI2)R2RMSEPAIC10.1260.0710.318-44.0220.1340.0710.296-44.1930.4550.0560.006-53.4640.6370.0460-61.6150.4980.0540.003-55.12

影響半徑/m擬合優(yōu)度(CI3RI3)R2RMSEPAIC10.1390.0710.281-44.3120.0910.0730.444-43.2430.2480.0660.088-47.0340.2890.0640.055-48.1450.3010.0640.048-48.47

影響半徑/m擬合優(yōu)度(CI4RI4)R2RMSEPAIC10.1200.0720.337-43.8820.1590.0700.229-44.8030.3750.0600.018-50.7440.4550.0560.006-53.4550.4120.0590.011-51.94

影響半徑/m擬合優(yōu)度(CI5RI5)R2RMSEPAIC10.0810.0730.488-43.0120.0840.0730.476-43.0730.1670.0700.211-44.9840.2000.0680.150-45.7950.2070.0680.140-45.96

3.3 競(jìng)爭(zhēng)強(qiáng)度、采伐強(qiáng)度與生長(zhǎng)量相關(guān)性分析

由圖2可知，以1/ln(BAI+1)為因變量，僅采用競(jìng)爭(zhēng)指數(shù)CI2r4作為自變量，個(gè)別點(diǎn)的實(shí)際值與預(yù)估值差別較大，如位于圖2最上方的點(diǎn)競(jìng)爭(zhēng)指數(shù)CI2r4為15.56，與預(yù)估值偏離較大，而采伐指數(shù)RI2r4僅為1.23，可能由于采伐指數(shù)較小造成其生長(zhǎng)量較低。通過(guò)圖4可以看出，在影響半徑4 m以?xún)?nèi)，伐樁根徑斷面積之和與對(duì)象木胸高斷面積生長(zhǎng)量呈正相關(guān)。

圖2 競(jìng)爭(zhēng)指數(shù)(CI2r4)、采伐指數(shù)(RI2r4)與生長(zhǎng)量的相關(guān)性

圖3 影響半徑4 m內(nèi)競(jìng)爭(zhēng)木、伐樁斷面積之和與生長(zhǎng)量的相關(guān)性

3.4 競(jìng)爭(zhēng)效應(yīng)模型的建立

為了進(jìn)一步說(shuō)明采伐對(duì)單木生長(zhǎng)的影響，對(duì)加入采伐指數(shù)的模型與只考慮競(jìng)爭(zhēng)指數(shù)的模型進(jìn)行對(duì)比(見(jiàn)表2)。由表2可知，影響半徑為4 m時(shí)，在引入采伐指數(shù)RI4r4之后，模型2比模型1判斷系數(shù)R2提高了33.8%，RMSE、AIC也有一定程度降低。加入采伐指數(shù)后，模型的擬合效果更好，說(shuō)明采伐對(duì)單木生長(zhǎng)的影響不可忽略。從模型的參數(shù)估計(jì)值分析，模型的參數(shù)估計(jì)在0.01水平下，差異顯著；a1為正值，說(shuō)明競(jìng)爭(zhēng)對(duì)對(duì)象木生長(zhǎng)起到抑制作用，b1為負(fù)值，說(shuō)明采伐對(duì)對(duì)象木生長(zhǎng)起到促進(jìn)作用，模型參數(shù)符合生態(tài)學(xué)原理。

根據(jù)建立的人工側(cè)柏幼齡林競(jìng)爭(zhēng)效應(yīng)模型繪出模型2的殘差分布圖，如圖4所示，殘差的散點(diǎn)分布是隨機(jī)的，說(shuō)明模型的擬合效果較好。

圖4 側(cè)柏競(jìng)爭(zhēng)效應(yīng)模型殘差分布

模型模型類(lèi)型參數(shù)估計(jì)值a0a1b1擬合優(yōu)度R2RMSEPAIC11/ln(BAI+1)=a0+a1+CI2r40.2945**0.0102**0.4760.054<0.001-56.2721/ln(BAI+1)=a0+a1CI2r4+b1RI2r40.3154**0.0137**-0.0101**0.6370.046<0.001-61.61

注：*表示差異顯著(P<0.05)，** 表示差異極顯著(P<0.01)。

4 討論

為衡量采伐對(duì)對(duì)象木生長(zhǎng)的影響，本文計(jì)算了與胸徑相關(guān)的直徑比、斷面積比、距離直徑比、距離斷面積比4種競(jìng)爭(zhēng)指數(shù)，以及張躍西[21]提出以生態(tài)學(xué)原理及自疏規(guī)律為基礎(chǔ)的競(jìng)爭(zhēng)指數(shù)，總共5種競(jìng)爭(zhēng)指數(shù)，討論了該林分內(nèi)適用的競(jìng)爭(zhēng)指數(shù)，并將競(jìng)爭(zhēng)指數(shù)計(jì)算公式推廣應(yīng)用構(gòu)建采伐指數(shù)，隨著采伐指數(shù)的增大，對(duì)象木的斷面積生長(zhǎng)量加大。在討論斷面積生長(zhǎng)量競(jìng)爭(zhēng)效應(yīng)模型時(shí)，采用模型的判斷系數(shù)(R2)、剩余標(biāo)準(zhǔn)差(RMSE)、赤池信息準(zhǔn)則(AIC)對(duì)模型進(jìn)行評(píng)價(jià)，5種競(jìng)爭(zhēng)與采伐指數(shù)中，斷面積比競(jìng)爭(zhēng)指數(shù)與采伐指數(shù)模型優(yōu)于其余4種，能更好的描述伐后林木的競(jìng)爭(zhēng)關(guān)系，這與郭恩瑩等[6]對(duì)杉木人工林4種緩沖區(qū)類(lèi)競(jìng)爭(zhēng)指數(shù)的研究結(jié)果相一致。

采用逐漸添加較遠(yuǎn)競(jìng)爭(zhēng)木、采伐木的信息，計(jì)算不同影響范圍內(nèi)的競(jìng)爭(zhēng)指數(shù)和采伐指數(shù)，并與采伐后5年間隔期胸高斷面積生長(zhǎng)量進(jìn)行回歸分析，根據(jù)生長(zhǎng)模型的R2、RMSE、AIC值判斷樹(shù)木個(gè)體間發(fā)生競(jìng)爭(zhēng)影響的最遠(yuǎn)距離，比通過(guò)R2值大小判斷最佳競(jìng)爭(zhēng)范圍更具有說(shuō)服力。

側(cè)柏人工幼齡林最佳競(jìng)爭(zhēng)范圍在4 m左右，與對(duì)象木樹(shù)高的一半相等，說(shuō)明采用樹(shù)高的一半作為對(duì)象木競(jìng)爭(zhēng)范圍具有參考意義[23-24]。這一指標(biāo)對(duì)北京地區(qū)中幼齡林撫育間伐強(qiáng)度選擇具有指導(dǎo)意義。史宇等[11]通過(guò)逐漸增大競(jìng)爭(zhēng)范圍觀察對(duì)象木競(jìng)爭(zhēng)指數(shù)的斜率變化，認(rèn)為天然側(cè)柏林最佳競(jìng)爭(zhēng)范圍為5 m；夏冰等[16]通過(guò)判定伐樁周?chē)a(chǎn)生生長(zhǎng)釋放的最遠(yuǎn)鄰體與伐樁所對(duì)應(yīng)的距離，得出山楊與青扦出現(xiàn)個(gè)體競(jìng)爭(zhēng)的空間范圍在4 m以?xún)?nèi)。

將采伐指數(shù)引入到競(jìng)爭(zhēng)效應(yīng)模型中，模型效果有明顯提高，影響半徑為4 m時(shí)，以CI2r4、RI2r4為自變量的模型比只有CI2r4的模型R2提高了33.8%，RMSE、AIC均有所降低，模型的參數(shù)估計(jì)在0.01水平下均差異顯著，競(jìng)爭(zhēng)指數(shù)與單木生長(zhǎng)量呈負(fù)相關(guān)，競(jìng)爭(zhēng)指數(shù)越大，單木生長(zhǎng)量越小，采伐指數(shù)與單木生長(zhǎng)量呈正相關(guān)，采伐指數(shù)越大，越受促進(jìn)對(duì)象木的生長(zhǎng)。本研究沒(méi)有對(duì)未采伐林分進(jìn)行對(duì)比分析，造成伐后單木生長(zhǎng)釋放的機(jī)理需要進(jìn)一步研究；沒(méi)有驗(yàn)證數(shù)據(jù)對(duì)模型進(jìn)行獨(dú)立檢驗(yàn)，模型的穩(wěn)定性與適用性需要進(jìn)一步研究驗(yàn)證。

[1] 孟憲宇.測(cè)樹(shù)學(xué)[M].3版.北京:中國(guó)林業(yè)出版社,2006，119(1):118-121.

[2] 李先琨，蘇宗明，歐祖蘭，等.元寶山冷杉群落種內(nèi)與種間競(jìng)爭(zhēng)的數(shù)量關(guān)系[J].植物資源與環(huán)境學(xué)報(bào),2002，11(1)20-24.

[3] 段仁燕，王孝安.太白紅杉種群鄰體范圍與鄰體競(jìng)爭(zhēng)強(qiáng)度的研究[J].西北植物學(xué)報(bào),2004，24(12):2335-2340.

[4] 段仁燕，王孝安.太白紅杉種內(nèi)和種間競(jìng)爭(zhēng)研究[J].植物生態(tài)學(xué)報(bào),2005,29(2):242-250.

[5] 段仁燕，王孝安，涂云博，等.太白紅杉鄰體競(jìng)爭(zhēng)研究方法的改進(jìn)[J].蘭州大學(xué)學(xué)報(bào)(自然科學(xué)版),2007，43(4):37-41.

[6] 郭恩瑩，吳保國(guó)，張瀚，等.杉木人工林緩沖區(qū)類(lèi)競(jìng)爭(zhēng)指數(shù)[J].東北林業(yè)大學(xué)學(xué)報(bào),2012，40(10):34-38.

[7] 鄒春靜，徐文鐸.沙地云杉種內(nèi)、種間競(jìng)爭(zhēng)的研究[J].植物生態(tài)學(xué)報(bào),1998，22(3):78-83.

[8] 劉強(qiáng)，李鳳日，董利虎.基于樹(shù)冠競(jìng)爭(zhēng)因子的落葉松人工林單木生長(zhǎng)模型[J].植物研究,2014,34(4):547-553.

[9] 李際平，房曉娜，封堯，等.基于加權(quán)Voronoi圖的林木競(jìng)爭(zhēng)指數(shù)[J].北京林業(yè)大學(xué)學(xué)報(bào),2015，43(3):61-68.

[10] 孜來(lái)比·買(mǎi)木提名，楊華，趙廣亮，等.單木競(jìng)爭(zhēng)指標(biāo)的研究進(jìn)展[J].西北林學(xué)院學(xué)報(bào),2012，27(6):152-158.

[11] 史宇，余新曉，岳永杰，等.北京山區(qū)天然側(cè)柏林種內(nèi)競(jìng)爭(zhēng)研究[J].北京林業(yè)大學(xué)學(xué)報(bào),2008，30(S2):36-40.

[12] 黃新峰，亢新剛，楊華，等.5個(gè)林木競(jìng)爭(zhēng)指數(shù)模型的比較[J].西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,40(7):127-134.

[13] 國(guó)慶喜，楊光.紅松天然種群鄰體影響半徑[J].應(yīng)用生態(tài)學(xué)報(bào),2006，17(12):2302-2306.

[14] 吳鞏勝，王政權(quán).水曲柳落葉松人工混交林中樹(shù)木個(gè)體生長(zhǎng)的競(jìng)爭(zhēng)效應(yīng)模型[J].應(yīng)用生態(tài)學(xué)報(bào),2000,11(5):646-650.

[15] 楊光，孫建，國(guó)慶喜.利用生長(zhǎng)釋放判定鄰體影響半徑[J].東北林業(yè)大學(xué)學(xué)報(bào),2009，37(6):20-22.

[16] 夏冰，鄧飛，周康，等.樹(shù)輪分析用于森林干擾強(qiáng)度推測(cè)的實(shí)例研究[J].植物資源與環(huán)境學(xué)報(bào),2001，10(1):1-6.

[17] 胡云云，閔志強(qiáng)，高延，等.擇伐對(duì)天然云冷杉林林分生長(zhǎng)和結(jié)構(gòu)的影響[J].林業(yè)科學(xué),2011，47(2):15-24.

[18] 陳輝榮，周新年，蔡瑞添，等.天然林不同強(qiáng)度擇伐后林分空間結(jié)構(gòu)變化動(dòng)態(tài)[J].植物科學(xué)學(xué)報(bào),2012，30(3):230-237.

[19] 蔣子涵，金光澤.擇伐對(duì)闊葉紅松林主要樹(shù)種徑向與縱向生長(zhǎng)的影響[J].生態(tài)學(xué)報(bào),2010，30(21):5843-5852.

[20] GRISSINO-MAYER H D. Evaluating crossdating accuracy: a manual and tutorial for the computer program COFECHA[J]. Tree-Ring Research,2001,57(2):205-221.

[21] 張躍西.鄰體干擾模型的改進(jìn)及其在營(yíng)林中的應(yīng)用[J].植物生態(tài)學(xué)與地植物學(xué)學(xué)報(bào),1993,17(4):352-357.

[22] WEINER J. Neighbourhood interference amongst pinus rigida individuals[J]. Journal of Ecology,1984,72(1):183-195.

[23] AMMER C, ZIEGLER C, KNOKE T. Assessing intra-and interspecfic competition in thickets of broadleaved tree species[J]. Allgemeine Forst und Jagdzeitung,2005,176(5):85-94.

[24] VON OHEIMB G, LANG A C, BRUELHEIDE H, et al. Individual-tree radial growth in a subtropical broad-leaved forest: The role of local neighbourhood competition[J]. Forest Ecology and Management,2011,261(3):499-507.

Individual Tree Growth-competition Model inPlatycladusorientalisYoung Plantation//

Ma Shiyou, Yang Hua(The Key Laboratory for Silviculture and Conservation, Beijing Forestry University, Beijing 100083, P. R. China); Zhang Heng, Shang Xuyang, Zhang Yong, Zhang Yongfu

(Beijing Ming Tombs Forest Farm)//Journal of Northeast Forestry University，2016,44(9)：9-13.

P.orientalisyoung plantation； Competition limits； Cutting； Competition index； Cutting index

1)國(guó)家林業(yè)局“948”引進(jìn)項(xiàng)目(2013-4-66)。

馬士友，男，1990年10月生，省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室(北京林業(yè)大學(xué))，碩士研究生。E-mail：majie131@163.com。

楊華，省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室(北京林業(yè)大學(xué))，副教授。E-mail：huayang8747@163.com。

2016年2月16日。

S758.1

責(zé)任編輯：王廣建。