賈 超

(山西省桑干河楊樹豐產(chǎn)林實(shí)驗(yàn)局，山西 大同 037006)

樟子松(Pinussylvestrisvar.mongolicaLitv.)為歐洲赤松的一個(gè)變種，具有耐寒、抗旱、適應(yīng)性強(qiáng)的特點(diǎn)。20世紀(jì)60年代，在雁北地區(qū)營(yíng)造了大面積小葉楊人工林，該林分多年后進(jìn)入成熟期，林木生長(zhǎng)幾乎停止，部分萎蔫近死亡。從70年代開始,山西省楊樹豐產(chǎn)林實(shí)驗(yàn)局利用樟子松、油松等樹種進(jìn)行小葉楊退化林分修復(fù)，取得了良好的效果。筆者以用于雁北地區(qū)退化林分修復(fù)的樟子松人工林為研究對(duì)象，收集樣木數(shù)據(jù)，通過備選模型進(jìn)行擬合優(yōu)選，研究林木胸徑、樹高及樹冠因子之間關(guān)系，確定雁北地區(qū)樟子松人工林立木材積方程，以期為該地區(qū)退化林分修復(fù)后的森林資源調(diào)查、林木蓄積量估算及科學(xué)經(jīng)營(yíng)提供一定的科學(xué)依據(jù)。

1 研究區(qū)概況

山西省雁北地區(qū)地處我國(guó)黃土高原東側(cè)，雁門關(guān)外內(nèi)外長(zhǎng)城之間，跨大同、朔州兩市，屬半干旱森林草原過渡地帶。氣候類型為半干旱溫帶氣候，海拔1 000 m～1 500 m.年均氣溫3 ℃～7 ℃，極端最低氣溫-40.4 ℃.無霜期100 d～300 d，八級(jí)以上大風(fēng)平均天數(shù)20 d～25 d.土壤以栗鈣土類和栗褐土類為主，還有大片鹽化潮土和潮土,一些地方分布有小片風(fēng)沙土、草甸鹽土和堿化鹽土。

2 材料與方法

2.1 樣本采集

為使樟子松立木材積模型在雁北地區(qū)范圍內(nèi)具有通用性，綜合分析樟子松在此區(qū)域的分布范圍，充分考慮不同氣候和立地條件對(duì)其的影響，選擇大同市云州區(qū)，朔州市朔城區(qū)、懷仁縣、應(yīng)縣、右玉縣5個(gè)(區(qū))縣作為樣木調(diào)查區(qū)，采取典型抽樣的方法設(shè)置調(diào)查樣地，見表1.

表1 調(diào)查樣地情況

本研究共調(diào)查樟子松365株，胸徑1.60 cm～27.00 cm，樹高1.50 m～13.70 m，冠長(zhǎng)2.20 m～10.00 m，冠幅0.95 m～7.75 m.樟子松樣木數(shù)量按照3 cm，6 cm，9 cm，12 cm，14 cm，16 cm，18 cm，20 cm，22 cm及27 cm共10個(gè)徑階均勻配置。每株樣木實(shí)測(cè)胸徑、樹高、冠幅和冠長(zhǎng)，并采用伐倒木區(qū)分求積法得到帶皮的樹干材積。

2.2 建立模型

利用采集的樣木數(shù)據(jù)，按胸徑大小排序，通過回歸估計(jì)方法，分析樟子松胸徑-樹高、胸徑-冠幅和冠長(zhǎng)-樹高之間的關(guān)系。擬建模型如表2，3，4所示。

表2 胸徑-樹高備選模型

2.3 建立二元立木材積模型

利用origin 9.1數(shù)據(jù)統(tǒng)計(jì)分析軟件進(jìn)行模型回歸和優(yōu)選，采用加權(quán)最小二乘法進(jìn)行模型擬合。擬建的二元立木材積模型如表5所示。

表3 胸徑-冠幅備選模型

表4 冠長(zhǎng)-樹高備選模型

表5 二元立木備選材積模型

表5中，模型1～7來源于我國(guó)林業(yè)行業(yè)標(biāo)準(zhǔn)《二元立木材積表編制技術(shù)規(guī)程》(LY/T 2012—2013)推薦的二元立木材積模型，8～16用于了解樟子松樹木各項(xiàng)因子與立木材積之間的相關(guān)程度。

3 結(jié)果與分析

3.1 模型擬合結(jié)果與評(píng)價(jià)指標(biāo)分析

雁北地區(qū)樟子松胸徑-樹高、胸徑-冠幅和冠長(zhǎng)-樹高之間的模型擬合結(jié)果如表6，7，8所示。

表6 胸徑-樹高模型擬合效果評(píng)價(jià)

表7 胸徑-冠幅模型擬合效果評(píng)價(jià)

表8 冠長(zhǎng)-樹高模型擬合效果評(píng)價(jià)

根據(jù)各個(gè)評(píng)價(jià)指標(biāo)綜合評(píng)判，可以看出在10個(gè)胸徑-樹高預(yù)選模型中，模型7擬合效果最好，決定系數(shù)為0.745 71，均方根誤差0.579 81.10個(gè)胸徑-冠幅預(yù)選模型中，模型1最佳，決定系數(shù)為0.587 72，均方根誤差0.625 86；8個(gè)冠長(zhǎng)-樹高預(yù)選模型中，模型3最佳，決定系數(shù)為0.651 93，均方根誤差1.106 55.樟子松3種模型最優(yōu)擬合結(jié)果的參數(shù)見表9.3種模型擬合效果中胸徑-樹高擬合效果最優(yōu)，可以認(rèn)為樟子松胸徑、樹高間存在最大相關(guān)性，胸徑-樹高最佳模型為H=(3.113 44×41 746.22+10.925 27D4.658 1)/(41 746.22+D4.658 1)；而胸徑-冠幅、冠長(zhǎng)-樹高模型擬合效果較差。

表9 樟子松3種模型最優(yōu)擬合結(jié)果

3.2 樟子松立木材積模型擬合結(jié)果與評(píng)價(jià)分析

樟子松二元立木材積模型擬合效果評(píng)價(jià)及參數(shù)結(jié)果見表10和第16頁表11.

表10 二元立木材積模型擬合效果評(píng)價(jià)

表11 樟子松立木材積模型擬合參數(shù)結(jié)果

由表10可知，16個(gè)樟子松立木材積模型中，前7個(gè)模型均具有較好的模擬效果。其中，以模型3擬合效果最好，決定系數(shù)為0.992 46，均方根誤差7.402 40×10-6，方程為V=1.978 78×10-16+4.830 19×10-5D2+1.610 06×10-5D2H-4.485 45×10-18H2+2.725 08×10-19DH2. 模型8，9，10，11為一元立木材積模型，擬合效果最好的為模型8，可以認(rèn)為立木材積與樹木胸徑關(guān)系最密切，其次為樹高。由模型12，13，14的擬合結(jié)果可以看出，樟子松樹木冠幅和樹高組合與立木材積擬合效果最好，冠幅與樹高相互配合能夠更好地反映樹木材積。比較模型13與16認(rèn)為，在樹高、冠幅基礎(chǔ)上增加冠長(zhǎng)對(duì)材積擬合影響不大。

3.3 模型適用性檢驗(yàn)

通過origin 9.1統(tǒng)計(jì)分析軟件，利用參與建模的樣本對(duì)立木材積模型進(jìn)行適用性檢驗(yàn)。以林木實(shí)際材積為因變量y，以模型計(jì)算出的理論材積為自變量x，建立一元線性回歸方程y=a+bx.經(jīng)計(jì)算，一元線性回歸方程為y=0.000 986+0.998 990x，F(xiàn)統(tǒng)計(jì)量1.266，小于F(0.05)=4.085.檢驗(yàn)樣本和實(shí)測(cè)樣本間沒有顯著差異，此模型適用于樟子松立木材積的計(jì)算。

4 結(jié)論與討論

本研究在雁北地區(qū)樟子松人工林分布區(qū)域調(diào)查的基礎(chǔ)上，研究了樹木胸徑、樹高、樹冠因子之間的關(guān)系。得到了3個(gè)最佳模型，其中，胸徑-樹高最佳模型為H=(3.113 44×41 746.22+10.925 27D4.658 1)/(41 746.22+D4.658 1)，能夠較好地反映胸徑樹高間關(guān)系。而樹高、冠幅與樹冠因子模型模擬效果較差，但冠長(zhǎng)與樹高存在較為緊密的聯(lián)系，擬合模型可以較好地反映此地區(qū)樟子松樹冠狀況。通過對(duì)16個(gè)立木材積模型的擬合，得到最優(yōu)立木材積模型為V=1.978 78×10-16+4.830 19×10-5D2+1.610 06×10-5D2H-4.485 45×10-18H2+2.725 08×10-19DH2，經(jīng)適用性檢驗(yàn)，可以為雁北地區(qū)樟子松人工林蓄積量調(diào)查提供科學(xué)依據(jù)。

截至2018年底，三北地區(qū)林分退化面積達(dá)到323.56×104hm2，約占工程保存林分面積的11.3%.退化林種中，以楊樹退化最為突出。在退化楊樹殘次林中，選用樟子松、油松等針葉樹種構(gòu)建針闊葉混交林，能夠形成穩(wěn)定的防護(hù)林，發(fā)揮防護(hù)功能。