栗生枝
(遼寧省林業調查規劃監測院,遼寧 沈陽 110122)
林分結構是指在未遭受嚴重干擾、經過長期的自然生長、枯損與演替的情況下,林分內許多特征因子,如直徑、樹高、形數、材積、樹冠等都具有一定的分布狀態,而且表現出較為穩定的林分結構規律[1]。林分結構內含著這些反映林分特征因子的變化規律,探討這些規律對森林經營技術、編制經營數表及林分調查都具有重要意義。
隨著林業科技進步及統計分析科學的發展,林分結構模型朝著復雜化、多樣化方向發展,從而從整體上提升了林分結構模擬與預測系統的性能及準確度,更能為科學營林和準確預估材積提供翔實的數字依據。
在林分內各種大小直徑林木按徑階的分配狀態,稱作林分直徑結構(Stand diameter structure),亦稱林分直徑分布(Stand diameter distribution)。在林分中,不同樹高的林木按樹高組的分配狀態,稱作林分樹高結構(Stand height structure),亦稱林分樹高分布(Stand height distribution)[2]。林分直徑結構和林分樹高結構都是最基本、最重要的林分結構,將直接影響樹木因子的變化。因此研究同齡純林的直徑及樹高結構規律在營林技術方面具有重要的指導意義。
西豐縣郜家店鎮位于鐵嶺東北部,距縣城22 km。郜家店鎮屬中溫帶亞濕潤區大陸性季風氣候,四季分明,氣候溫和,雨量充沛。年平均氣溫5.1 ℃。1月平均氣溫-17 ℃,最低氣溫-41.0 ℃;7月平均氣溫23.2 ℃,最高氣溫35.2 ℃。年平均降水量738 mm, 無霜期135 d,土層厚度50 cm左右,腐殖質層厚3 cm左右。
試驗地、對照地全部為人工落葉松純林,林相整齊,易于作業,便于管理。
在西豐縣郜家店鎮林區選擇坡度平緩、有代表性的長白落葉松(Larixolgensis)人工林,分不同林齡(18年生、22年生、38年生)設置樣地,樣地面積25 m×25 m。共設樣地132塊。其中18年生林分設置樣地36塊,22年生林分設置樣地48塊,38年生林分設置樣地48塊,樣地概況見表1。

表1 樣地概況
在標準地內對所有喬木進行每木檢尺,詳細記錄每株樹的胸徑、樹高、冠幅、枝下高等林木因子以及坡度、坡向、海拔等立地因子。并按胸徑2 cm為一個徑階,樹高為1 m為樹高距分組記錄數據。
2.3.1 正態分布法
(1)正態分布概率密度:標準正態分布的概率密度為:
(1)

m——林木實測直徑/樹高值;


(2)
在正態分布中,以數學期望及方差為分布參數,只要林分算術平均直徑(樹高)及標準差計算出來之后,然后借助已有的正態分布函數表計算出林分直徑正態分布中各徑階株數理論值。最終與實測株數對比結果進行卡方檢驗,研究擬合效果。
2.3.2 卡方檢驗 林分株數按徑階和樹高分布的現實結構需進行卡方檢驗才能確定是否真正符合擬合曲線。
設x= (x1,x2,…xn) 為長白落葉松林分各徑階(樹高)實際株數所構成的數列,y= (y1,y2…yn) 為利用擬合曲線分布得到長白落葉松林分各徑階(樹高)理論株數所構成的數列。則χ2統計量為:
根據徑階中值及其對應株數值繪制直徑分布曲線,如圖1。

由圖1可以看出,3個林齡的林分直徑分布曲線的具體形狀雖有些差異,但都形成一條以林分平均直徑為峰點、中等大小的林木株數占多數、向其兩端徑階的林木株數逐漸減少的單峰山狀曲線,類似服從正態分布曲線。因此,可以用正態分布對直徑分布曲線進行擬合。
利用林分的算術平均直徑及標準差(表2)就可以擬合出林分直徑的正態分布函數,根據該函數,即可求出林分直徑正態分布中各徑階株數理論值(表3)。

表2 各林分直徑的期望和標準差值

表3 長白落葉松人工林直徑分布擬合結果 株
對擬合效果進行分析及卡方檢驗,結果見表4。

表4 長白落葉松人工林直徑分布擬合效果分析
注:χ2值中帶* 的為符合分布的標準地
由表4可知,3個林分的直徑分布曲線均可以用正態分布很好的擬合。與相對直徑法所得結果相同,18年、22年、38年生長白落葉松林分直徑的偏度系數均為正值,其絕對值均分布在[0.2,0.4] 區間,說明該標準地林木直徑分布為左偏。從峰度分布來看,18年、22年生林分的峰度為正值, 說明它們比正態分布曲線略為陡峭,38年生林分的峰度為負值,說明該曲線比正態分布曲線略為平坦。其中,18年生林分的直徑分布曲線偏度、峰度最大,38年生林分的偏度、峰度最小,其曲線更近似于正態分布。說明隨著林分林齡的增長,林分平均直徑逐漸增大,直徑分布曲線的峰度由大變小,形成寬而平的分布曲線,更接近于正態分布曲線。出現這種情況的原因可能是:幼齡階段,林分密度低,光照、營養、空間充足,林木生長必須的條件都可以滿足,競爭壓力小,所以林木生長速率差別不大,其胸徑均集中在林分平均直徑左右,中等大小的林木株數占絕大多數。但隨著年齡的增加,林木生長所需的空間和養分增多,而林地所能提供的空間、養分有限,林木競爭逐漸加劇,此時,林木生長速率出現差異,林木生長量大小開始分化,隨著時間的推移,這種分化越來越明顯,因而出現林木胸徑分布分散化,即林分直徑分布曲線峰度減小,形成寬而平的分布曲線。
根據樹高值及其對應株數值繪制樹高分布曲線,如圖2。

由圖2分析表明,3塊長白落葉松人工林的樹高分布具有明顯的結構及變化規律,呈現出近似于以該林分平均高的株數最多的非對稱性的山狀曲線。

表5 各林分樹高的期望和標準差值

表6 長白落葉松林樹高分布擬合結果(株)
近似服從正態分布曲線。因此,可以用正態分布對長白落葉松林樹高分布曲線進行擬合。
利用林分的算術平均樹高及標準差(表5)就可以擬合出林分樹高的正態分布函數,根據該函數,即可求出林分樹高正態分布中各徑階株數理論值(表6)。
對擬合結果進行分析和卡方檢驗,結果如表7。

表7 長白落葉松林樹高分布擬合效果檢驗
注:χ2統計量值中帶* 的為符合分布的標準地
依表7所示,正態分布對3個林分的樹高分布曲線擬合效果非常好。3個林齡的長白落葉松林分樹高的偏度系數均為正值,其絕對值分布在[0.1,0.5] 區間,說明該標準地林木樹高分布為左偏,說明林分內樹高中等和偏低的林木較多,而對林分直徑結構的分析表明林分內中小徑階的林木占多數,胸徑、樹高剛好可以對應,說明3個林齡的林分內中等及偏小的樹占多數。但相對樹高法得出的結論是38年生林分的樹高分布為右偏,可能是由于相對樹高法擬合效果不理想,導致誤差產生所致,應以正態分布所得結果為準。從峰度分布來看,峰度為負值, 說明它們比正態分布曲線略為平坦。與林分直徑分布曲線不同的是,樹高分布曲線峰度隨著林分年齡的變化基本保持不變,并不隨著林分年齡的變化而有明顯增大或減小的趨勢。這也說明,與胸徑生長不同,樹高生長受林分密度影響較小,其很大程度上取決于立地條件的優劣。
根據以上分析和檢驗結果,得出正態分布對林分直徑和樹高結構的擬合效果好,靈活性大、適應性強,且隨林分林齡的增大,擬合效果越來越好。根據正態分布的擬合結果得出,西豐縣郜家店鎮地區不同林齡的長白落葉松人工林林分直徑和樹高分布均為單峰左偏正態分布,即各個林齡的林分內中等及偏小的樹占多數。直徑和樹高結構分析是林分蓄積量和出材量表編制的基礎,即各徑階單株木材積乘以相應徑階株數得徑階蓄積合計,各徑階蓄積合計得林分蓄積合計;各徑階蓄積乘以相應徑階不同材種的出材率得徑階材種出材量,徑階材種出材量合計得林分材種出材量合計。因此,本研究為建立西豐縣郜家店地區長白落葉松人工林收獲預估體系提供一定的理論依據。
本文僅以正態分布法對樣地直徑樹高結構進行分析,分析方法較為單一,應進一步運用相對直徑法和相對樹高法進行比較分析,得出更能反映西豐縣郜家店地區長白落葉松人工林林分結構的方法。