典型中亞熱帶天然闊葉林各林層直徑分布及其變化規律*

莊崇洋 黃清麟 馬志波，2 鄭群瑞 王 宏

(1.中國林業科學研究院資源信息研究所 國家林業局林業遙感與信息技術重點實驗室 北京 100091； 2.中國林業科學研究院森林生態環境與保護研究所 北京 100091； 3. 建甌萬木林省級自然保護區管理處 南平353105)

?

典型中亞熱帶天然闊葉林各林層直徑分布及其變化規律*

莊崇洋1黃清麟1馬志波1，2鄭群瑞3王 宏1

(1.中國林業科學研究院資源信息研究所 國家林業局林業遙感與信息技術重點實驗室 北京 100091； 2.中國林業科學研究院森林生態環境與保護研究所 北京 100091； 3. 建甌萬木林省級自然保護區管理處 南平353105)

【目的】 研究典型中亞熱帶天然闊葉林各林層直徑分布及其變化規律，為揭示天然林結構、輔助天然林相關測樹因子的調查和指導天然林經營規劃奠定理論基礎。【方法】 采用最大受光面法劃分5塊中亞熱帶天然闊葉林林層，利用Shapiro-Wilk檢驗(S-W檢驗)對各林層(包括全林分和各亞層)直徑分布進行正態性檢驗，采用偏度和峰度描述各林層直徑分布圖形特征，利用Meyer負指數函數和Weibull分布函數對各林層直徑分布進行擬合，根據卡方檢驗結果選擇擬合效果較好的函數擬合各林層直徑分布，分析各林層直徑分布及其變化規律。【結果】 S-W檢驗表明，各標準地全林分、第Ⅲ亞層和第Ⅱ亞層直徑分布均不服從正態分布； 1-3號標準地第Ⅰ亞層直徑分布服從正態分布，4和5號標準地第Ⅰ亞層直徑分布不服從正態分布； 各標準地內，各亞層的S-W值隨亞層高度的升高而增大。從總體上看，除5號標準地第Ⅱ、Ⅲ亞層的偏度和峰度及2號標準地第Ⅰ、Ⅱ亞層的峰度外，各亞層偏度和峰度的絕對值均隨林層高度的升高而減小。利用Meyer負指數函數擬合全林分(所有標準地均未通過卡方檢驗)和第Ⅱ亞層(2和4號標準地通過卡方檢驗，1、3和5號未通過)時適應性較差； 而在擬合第Ⅰ亞層(所有標準地均通過卡方檢驗)和第Ⅲ亞層(3、4和5號標準地通過卡方檢驗，1和2號未通過)時有較好的適應性。Weibull分布函數對各林層的直徑分布具有較好的擬合效果，除了1號標準地全林分直徑分布外，1號標準地剩余亞層和其他標準地各林層均通過Weibull分布的卡方檢驗。利用Weibull分布函數擬合典型林分各標準地直徑分布后發現，隨著各亞層高度的升高，各亞層直徑分布由倒“J”形曲線向右偏山狀曲線和正態分布曲線過渡，頂峰向右移動，同時峰值隨之減小。【結論】 劃分林層后，Weibull分布函數在擬合中亞熱帶天然闊葉林各林層(包括全林分和各亞層)直徑分布時比Meyer負指數函數具有更好的擬合效果； 各亞層S-W檢驗中P值、偏度和峰度以及Weibull函數中各亞層直徑分布圖形波峰的變化情況說明各亞層直徑分布有隨亞層平均胸徑和平均高度增大而逐漸向正態分布過渡的趨勢； 各林層(包括全林分和各亞層)直徑分布間的差異說明分林層研究直徑分布的必要性。

中亞熱帶； 典型天然闊葉林； 分層； 直徑分布； 變化規律

直徑分布是指林分內各種大小直徑林木按徑階的分配狀態(孟憲宇， 2006)。林分直徑分布規律是林分結構的基本規律之一，是林分環境長期相互作用的結果(周國模等， 1994)，研究林分直徑分布規律能為造林設計(Meyer， 1952； Leak， 1964)、確定林分發展階段(Goelzetal.， 2002)、建立林分生長和收獲模型(Clutteretal.， 1984； Bordersetal.， 1987； 1990)提供依據； 同時，林分直徑分布在林分結構、森林調查和研制林分材種出材量表中具有重要作用，林分直徑分布規律可作為營林技術評價指標(孟憲宇， 1991)。研究直徑分布以同齡林直徑分布為主，天然林因受林分自身的演替過程、樹種組成、樹種特性、立地條件、更新過程、自然災害等因素的影響而變得復雜多樣(孟憲宇， 2006)。

研究林分直徑分布規律主要從概率分布函數和理論模型入手，經歷了從正態分布、負指數分布、Sb分布(Baileyetal.， 1973； Hafleyetal.， 1977)、β分布(Clutteretal.， 1965)、Г分布(Nelson， 1964)、Weibull分布(Baileyetal.， 1973； 孟憲宇， 1985; Cao， 2004)等為主的靜態模擬方法，過渡到轉移矩陣、參數預測(Bailey， 1980)和參數回收(Burketal.， 1984)等為主的動態擬合過程。從模型參數看，直徑分布模型研究又可分為參數法(Bordersetal.， 1987)和非參數法(Haraetal.， 1997)，具體包括最相似回歸法、核密度估計法(Dressieretal.， 1989； 徐健君， 1999； 郭芳等， 2008)、相對直徑法(孟憲宇等， 1996)、理論方程法(Ishikawa， 1998； 段愛國等， 2003)、聯立方程組法(Bordersetal.， 1987； Maltamoetal., 2000； 劉福香， 2013)、概率密度函數法(Baileyetal., 1973)、有限混合分布方法(Liuetal.， 2002； 劉福香， 2013)等。張建國等(2004)從林分直徑分布模型的選擇和參數求解2方面對直徑分布的研究進展做了詳細綜述。

從林層角度研究直徑分布的報道較少，劉健等(1996)分林層研究天然針闊混交林的直徑分布、許彥紅等(2004)探討西雙版納熱帶雨林植被亞型各林層的直徑分布時，只是簡單介紹了各林層(包括全林分和各亞層)的直徑分布，沒有探討直徑分布隨林層的變化規律； Knoke等(2001)只是給出了不同林齡的直徑分布圖，并未對其進行進一步探討。亞熱帶是我國植物資源最豐富的地區之一(中國植被編輯委員會， 1980)。 我國亞熱帶常綠闊葉林是世界上較為罕見的植被類型，橫跨11個緯度，面積約2 280 000 km2， 其中有2/3(約1 540 000 km2)為中亞熱帶常綠闊葉林(黃清麟， 1998a； 黃清麟等， 1999a)，林分結構較為復雜。對于復雜的天然異齡林，如何更深入地了解林分結構一直是研究的熱點，而林分直徑結構是林分中最重要、最基本的結構，在復雜的天然林異齡林中分林層研究直徑分布對揭示天然林結構、輔助天然林相關測樹因子的調查和指導天然林經營規劃具有重要意義。本文以典型中亞熱帶天然闊葉林為研究對象，在劃分林層的基礎上，探討各林層直徑分布及其隨林層的變化規律。

1 研究區概況

按以下標準選擇典型中亞熱帶天然闊葉林為試驗林分： 1) 達到地帶性頂級群落； 2) 基本未經過人為干擾； 3) 基本符合中亞熱帶天然闊葉林理想結構的標準(黃清麟等， 2003)； 4) 群落類型多樣且相對集中。

建甌萬木林省級自然保護區是1957年林業部根據人大代表提議劃定的全國首批19個天然森林禁伐區之一，是我國最早的自然保護區，經過持續600多年的封禁保護(何友釗， 1989)，現已自然演替成最為典型的中亞熱帶森林生態系統(常綠闊葉林)。 保護區位于閩北建甌市境內，地理坐標為118°08′22″—118°09′23″E，27°02′28″—27°03′32″N，屬武夷山南坡低山丘陵，海拔230～556 m； 中生代燕山運動侵入的花崗巖為主要成土母巖，紅壤，立地類型以Ⅱ類地為主； 保護區屬中亞熱帶海洋性季風氣候，熱量豐富，降水充沛，季風顯著，四季分明； 年平均氣溫18.7 ℃，1月平均氣溫13.8 ℃，極端最低氣溫-5.9 ℃，7月平均氣溫28.3 ℃，極端最高氣溫40.7 ℃； 年平均降水量1 700 mm，6月最多，雨日年平均166天。保護區主要有殼斗科(Fagaceae)、樟科(Lauraceae)、山茶科(Theaceae)、木蘭科(Magnoliaceae)、金縷梅科(Hamamelidaceae)和杜英科(Elaeocarpaceae)等常綠闊葉樹種，并有珍稀樹種觀光木(Tsoongiodendronodorum)、樂東擬單性木蘭(Parakmerialotungensis)、沉水樟(Cinnamomummicranthum)、閩楠(Phoebebournei)、紅豆樹(Ormosiahosiei)、藍果樹(Nyssasinensis)等(高峻等， 1998)。

2 數據調查與研究方法

2.1 數據調查

經全面踏查，在建甌萬木林省級自然保護區內較平坦的地段中選擇典型中亞熱帶天然闊葉林作為試驗林分，于2015年設置5塊50 m×50 m的標準地，并進行1次數據調查。對標準地內胸徑≥5.0 cm的林木進行定位，調查林木的樹種、胸徑(cm)、樹高(m)、枝下高(m)、冠長(m)、冠幅長(m)、冠幅寬(m)和冠形等因子，采用最大受光面法(maximum light receiving plane，MLRP)判讀林木的林層歸屬(莊崇洋， 2016)。標準地數據統計如表1所示。

表1 各標準地調查因子和林層高度Tab.1 Measure factors and storey height of sample plots

樹種蓄積根據福建省主要樹種二元材積表計算得到，林層下限值由最大受光面法得到。Shannon-Winener指數按彭少麟等(1983)提出的公式計算：

(1)

式中： SW表示Shannon-Wiener多樣性指數； 平均高為斷面積加權平均高；S表示物種數量；ni表示某一物種的個體數；N表示全部物種的個數；Pi表示第i個種個體數量占所有個體數量的百分比。

整理各林層的直徑數據，刪除枯死木、斷梢木等，并以4 cm為徑階對各亞層的林木進行歸組，計算各林層的平均胸徑、標準差、最大徑階和最小徑階。由于4徑階內只有5～6 cm的林木，為不完整徑階，可能會對直徑分布規律的分析造成干擾，因此在討論直徑分布時將其舍棄。具體數據見表2。

2.2 研究方法

采用最大受光面法(MLRP)劃分林分林層，并在此基礎上利用Shapiro-Wilk檢驗(S-W檢驗)對各林層(包括全林分和各亞層)直徑分布進行正態性檢驗，采用偏度(Skewness，SK)和峰度(Kurtosis，KT)描述各林層直徑分布圖形特征，利用Meyer負指數函數和Weibull分布函數擬合各林層直徑分布，并根據卡方檢驗結果選擇擬合效果較好的函數擬合各林層直徑分布，分析各林層直徑分布及其變化規律。

2.2.1 最大受光面法 最大受光面法(黃清麟， 1998b； 黃清麟等， 1995； 1999b； 2003)指依據林木樹冠能否接受到光照和是否明顯突出而進行的典型中亞熱帶天然闊葉林林層劃分方法。首先依據典型林分的林木樹冠(林隙內的林木樹冠除外)能否直接接受到垂直光照將典型林分的所有林木樹冠劃分為受光層(所有能接受到光照的林木樹冠組成的一個層次，包括第Ⅰ和第Ⅱ亞層2個亞層)和非受光層(所有不能接受到光照的林木樹冠組成的一個層次，指第Ⅲ亞層)，其次在受光層中依據其林木樹冠是否明顯突出劃分為林木樹冠明顯突出的受光層(第Ⅰ亞層)和林木樹冠不明顯突出的受光層(第Ⅱ亞層)。由于受光層與非受光層之間的交界面是其高度以上所有能接受到光照的林木樹冠垂直投影面積(受光面積)最大的水平截面(簡稱最大受光面)，所以將這種方法稱為最大受光面法。

表2 各林層直徑數據Tab.2 Diameter data of storeys cm

2.2.2 Shapiro-Wilk檢驗 利用Shapiro-Wilk檢驗(S-W檢驗)對各林層直徑分布進行正態性檢驗。原假設H0： 總體服從正態分布； 備擇假設H1： 總體不服從正態分布。利用S-W檢驗原假設H0的過程如下。

先將檢驗樣本值升序排列：

x1≤x2≤……≤xn。

再計算S-W檢驗的統計量W：

(2)

在95%的置信區間內，如果W≥0.05，則表示接受原假設，服從正態分布； 否則，拒絕原假設，不服從正態分布。

2.2.3 偏度和峰度 采用偏度(SK)和峰度(KT)描述各林層直徑分布圖形特征：

SK和KT是描述分布圖形對稱性和陡緩程度與正態分布圖形差異程度的指標，其絕對值越小，數據分布圖形與正態分布圖形越接近。

2.2.4 Meyer負指數函數 Meyer負指數函數在描述異齡林直徑分布時有較好的適應性，其公式為：

(5)

式中：Y表示每個徑階的林木株數；x表示徑階； e表示自然對數的底；a，K表示直徑分布特征常數。

2.2.5 Weibull分布函數 Weibull分布函數的密度函數為：

(6)

式中：a為位置參數；b為尺度參數；c為形狀參數。

對利用Meyer負指數函數和Weibull分布函數求得的理論株數進行χ2檢驗，公式如下：

(7)

選擇最優模型擬合各林層直徑分布，分析各林層直徑分布及其變化規律。

3 結果與分析

3.1 S-W檢驗結果分析

由表3可知，全林分直徑分布S-W檢驗的顯著度P均小于0.05，表明典型試驗林分全林分直徑分布不服從正態分布，這是異齡林直徑分布的特征之一(孟憲宇， 2006)。從分層角度看，對于第Ⅰ亞層，1、2和3號標準地直徑分布S-W檢驗的顯著度P分別為0.166、0340和0.590，均大于0.05，表明這3個標準地第Ⅰ亞層直徑分布服從正態分布假設； 而4和5號標準地直徑分布S-W檢驗的顯著度P分別為0.014和0.013，都小于0.05，表明這2個標準地第Ⅰ亞層直徑分布不服從正態分布。所有試驗林分第Ⅱ亞層和第Ⅲ亞層直徑分布S-W檢驗的顯著度P均小于0.05，表明典型試驗林分的第Ⅱ亞層和第Ⅲ亞層直徑分布均拒絕正態分布假設。從典型試驗林分各林層的顯著度結果可以看出，PⅠ(第Ⅰ亞層顯著度)>PⅡ(第Ⅱ亞層顯著度)>PⅢ(第Ⅲ亞層顯著度)。

3.2 偏度、峰度結果分析

分析典型試驗林分各亞層的偏度(SK)和峰度(KT)可以看出，除了2號標準地第Ⅰ、Ⅱ亞層的KT(第Ⅰ亞層為-0.40，第Ⅱ亞層為-0.12)及5號標準地第Ⅱ、Ⅲ亞層的SK和KT(第Ⅱ亞層SK和KT分別為1.65、2.82； 第Ⅲ亞層SK和KT分別為1.34、1.03)外，其他各標準地內各亞層SK和KT的絕對值均隨亞層高度的升高而減小，即SKⅢ(第Ⅲ亞層的SK)>SKⅡ(第Ⅱ亞層的SK)>SKⅠ(第Ⅰ亞層的SK)，KTⅢ(第Ⅲ亞層的KT)>KTⅡ(第Ⅱ亞層的KT)>KTⅠ(第Ⅰ亞層的KT)(表3)。

3.3 Meyer負指數函數結果分析

由卡方檢驗結果(表3)可知，對于典型試驗林分全林分直徑分布，只有2號標準地服從Meyer負指數函數，其余4塊標準地均拒絕。從亞層的角度看，各標準地第Ⅰ亞層直徑分布均服從Meyer負指數函數，但其參數值均較小，說明各徑階內株數變化情況不大；對于各標準地的第Ⅱ亞層，2和4號標準地服從Meyer負指數函數，其余3塊標準地均拒絕； 對于各標準地的第Ⅲ亞層，3、4和5號3塊標準地服從Meyer負指數函數，其余2塊標準地拒絕。

3.4 Weibull分布函數結果分析

由卡方檢驗結果(表3)可知，除了1號標準地全林分直徑分布不服從χ2檢驗外，1號標準地亞層和其他標準地各林層均通過卡方檢驗，說明無論是全林分還是分亞層，Weibull分布均可很好地擬合典型試驗林分各林層的直徑分布。

表3 各研究方法結果數據①Tab.3 Value of research methods

① *表示服從假設分布。*means obeying the hypothesis distribution.

3.5 各林層直徑分布圖

相較于Meyer負指數函數，Weibull分布函數可以較好地描述標準地內各林層(包括全林分和各亞層)的直徑分布，因此選擇Weibull函數擬合各標準地內各林層的直徑分布情況，結果如圖1所示。

圖1中4徑階內只有5～6 cm的林木，為不完整徑階，在本次直徑分布的擬合時已將其舍棄。從分布圖的形狀上看，各標準地表現基本一致，隨著各亞層高度的升高，各亞層直徑分布由倒“J”形曲線向右偏山狀曲線和正態分布曲線過渡，頂峰向右移動，同時峰值隨之減小。各亞層直徑分布與全林分直徑分布有一定的重疊區域，第Ⅲ亞層小于12徑階和第Ⅰ亞層大于40徑階區域與全林分直徑分布基本一致。全林分和各亞層直徑分布差異最大的是12～40徑階，在這個區間內除了5號標準地外，其他4個標準地3個亞層的直徑分布都有交集，說明在該段區間內，林分結構(如樹高關系)較為復雜，可能會對森林經營、數表的編制和林分調查造成干擾，且各亞層的直徑分布也差別較大，說明從亞層角度研究直徑分布是有必要的。

圖1 各標準地直徑分布和Weibull分布函數擬合Fig.1 Line charts of diameter class and fitting curves of Weibull distribution function of sample plots

4 討論

S-W檢驗的顯著度P是判斷函數是否服從正態分布的數量化指標，當P≥0.05時，函數為正態分布，當P<0.05時，P越接近0.05越趨向于正態分布。本研究在典型林分中各亞層的顯著度P隨亞層高度的升高而增大，說明越高的亞層，層內林木的直徑分布越趨向于正態分布，各亞層SK和KT的絕對值隨林層高度的升高而減小的變化規律同樣說明這個趨勢。

從生態學特性看，林齡較大、個體較大的林木個體一般占據上層空間，林齡較小、個體較小的林木個體一般處于下層空間(徐化成， 2004)。本研究在分林層后，各亞層平均樹高HⅠ>HⅡ>HⅢ，平均胸徑DⅠ>DⅡ>DⅢ，因此，可判斷各亞層平均林齡為AⅠ>AⅡ>AⅢ。在同齡林中，直徑分布的偏度和峰度隨著林齡和林分算術平均直徑的增大而減少，其直徑分布由非正態分布過渡到正態分布(孟憲宇， 2006)。而在典型試驗林分(異齡復層林)中，劃分林層后各標準地S-W檢驗的顯著度PⅠ最大，SKⅠ和KTⅠ絕對值最小，表明各標準地內各亞層直徑分布可能表現出與同齡林相似的結果，即隨著平均樹高、平均直徑和林齡增大，其直徑分布趨向于正態分布。

利用分布函數擬合各林層(包括全林分和各亞層)的直徑分布時，Weibull分布函數表現出比Meyer負指數函數更好的適應性。采用Weibull分布函數擬合各林層直徑分布后發現，隨著亞層高度和平均直徑的增大，各亞層直徑分布的頂峰向右移動，同時峰值隨之減小，各亞層的徑階跨度隨平均直徑的增加而增大，這與Gingrich(1967)利用正態分布研究山地硬闊葉林直徑分布時的結論類似。

本研究同樣嘗試采用優勢高法劃分分層并討論各林層直徑分布的變化規律。從劃分結果看，優勢高法將林分喬木層分為3個亞層，第Ⅰ亞層下限值與最大受光面法較為接近，最大差異不超過2.3 m，株數最大差異不超過9株，最大受光面法具有較高的亞層下限值和較少的林木株數。但2種方法得出的第Ⅱ亞層的下限值和株數差距較大，最大受光面法的第Ⅱ亞層下限值在16.0～17.0 m之間，而優勢高法的下限值在10.8～12.8 m之間，同一標準地差異最大為6.2 m，最小為3.7 m，其余標準地差異在5 m左右，株數差異最小為32株，最大為61株，最大受光面法具有較低的亞層下限值和較多的林木株數。第Ⅲ亞層的株數差異同樣較大，與第Ⅱ亞層一致，最小為32株，最大為61株，最大受光面法具有較少的林木株數。2種劃分方法各亞層的平均胸徑和平均樹高同樣具有一定差異，第Ⅰ亞層總體差距最小，第Ⅱ、Ⅲ亞層差距則較為明顯，總體來說，同一標準地內采用最大受光面法具有較大的亞層平均胸徑和平均樹高。

采用優勢高法分層后，各林層的直徑分布規律基本與最大受光面法分層后的直徑分布規律類似，但優勢高法分層劃分亞層后第Ⅲ亞層的蓄積量均低于30 m3·hm-2，不符合國家標準《森林資源規劃設計調查技術規程》(GB/T 26424—2010)中林層分層的蓄積要求，且這種方法是相對機械的劃分方法，是高度上的簡單劃分，較難從生物特性方面來解釋林分亞層的自然分異。而最大受光面法劃分的各林層指標均符合國家標準(GB/T 26424—2010)中的要求，從光角度劃分林層也具有一定的生物學意義，且最大受光面法在野外便于操作，觀察者只要觀察樹冠是否接受直射光照射即可劃分受光層和非受光層，而受光層中根據林木樹冠受光情況也可大致分層。

在研究各亞層直徑分布時，第Ⅰ、Ⅱ亞層林木較少，可能會對研究結果造成干擾，因此，本研究也嘗試將相對位置較近、樹種組成類似的2和3號標準地、4和5號標準地合成2個新的標準地，用以分析第Ⅰ、Ⅱ亞層的直徑分布，結果發現其規律與原標準地基本相似，即隨著平均直徑的增大，偏度和峰度減小，S-W檢驗的顯著度增大，直徑分布向正態分布發展。

5 結論

本研究在采用最大受光面法劃分林層的基礎上，對典型中亞熱帶天然闊葉林各林層(包括全林分和各亞層)的直徑分布進行探討。劃分林層后，Weibull分布函數擬合各林層直徑分布的適用性比Meyer負指數函數要好。從典型中亞熱帶天然闊葉林各林層直徑分布分析結果可以看出，全林分和各亞層在某些徑階區間上的差異和各林層直徑分布曲線上的較大差異，證明了分層研究直徑分布的必要性。

典型中亞熱帶天然闊葉林劃分林層后，各亞層S-W檢驗的顯著度均隨亞層平均胸徑和平均樹高的增加而增大，偏度和峰度的絕對值隨亞層平均胸徑和平均樹高的增加而減小，Weibull函數圖形的頂峰右移、峰值下降等情況，說明亞層直徑分布有向正態分布過渡的趨勢。

段愛國,張建國,童書振. 2003. 6種生長方程在杉木人工林林分直徑分布上的應用.林業科學研究, 16(4): 423-429.

(Duan A G, Zhang J G, Tong S Z. 2003. Application of six growth equations on stands diameter structure of Chinese fir plantations. Forest Research,16(4): 423-429. [in Chinese])

高 峻,鄭群瑞,楊斌生,等. 1998. 福建萬木林自然保護區植被研究.上海師范大學學報：自然科學版,27(1): 77-82.

(Gao J, Zheng Q R, Yang B S,etal. 1998. Study on the vegetation Wanmulin Reserve of Jian’ou of Fujian Province. Journal of Shanghai Teachers University:Natural Sciences,27(1): 77-82. [in Chinese])

郭 芳,趙新振,殷鳴放,等. 2008. 非參數核密度估計在闊葉紅松天然林直徑分布研究中的應用.河南農業大學學報, 42(3): 303-307.

(Guo F, Zhao X Z, Yin M F,etal. 2008. Application of nonparametric kernel density estimation to research into diameter distribution of broadleavedPinuskoraiensisnatural stand. Journal of Henan Agricultural University, 42(3): 303-307. [in Chinese])

何友釗. 1989. 建甌縣萬木林保護區史事考.林史文集, 1(1): 139-140.

(He Y Z. 1989. Historical issue of Wanmulin Reserve of Jian’ou County. Forest History Study, 1(1): 139-140. [in Chinese])

黃清麟. 1998a. 亞熱帶天然闊葉林經營中的五大誤區.世界林業研究, 11(4): 31-34.

(Huang Q L. 1998a. Misunderstanding of the management of natural broadleaved forest in the subtropics.World Forestry Research, 11(4): 31-34. [in Chinese])

黃清麟. 1998b. 中亞熱帶天然闊葉林可持續經營技術研究. 北京: 北京林業大學碩士學位論文.

(Huang Q L. 1998b. Study on sustainable management technology of natural broad-leaved forest of the mid-subtropical zone. MS thesis of Beijing: Beijing Forestry University. [in Chinese])

黃清麟,李元紅. 1999a. 中亞熱帶天然闊葉林研究綜述.福建林學院學報, 19(2): 189-192.

(Huang Q L, Li Y H. 1999a. A brief review on the research of natural broad-leaved forest in the mid-subtropical zone. Journal of Fujian College of Forestry, 19(2): 189-192. [in Chinese])

黃清麟,董乃鈞,李元紅. 1999b. 福建中亞熱帶天然闊葉林的主要類型與特征. 山地學報, 17(4): 368-374.

(Huang Q L, Dong N J, Li Y H. 1999b. Main types and characteristics of natural broad leaved forest of the mid-subtropical zone in Fujian Province. Journal of Mountain Science, 17(4): 368-374. [in Chinese])

黃清麟,李志明,鄭群瑞. 2003. 福建中亞熱帶天然闊葉林理想結構探討. 山地學報, 21(1): 116-120.

(Huang Q L, Li Z M, Zheng Q R. 2003. A discussion on the ideal structure of natural broad-leaved forest in the mid-subtropical zone of Fujian Province. Journal of Mountain Science, 21(1): 116-120. [in Chinese])

黃清麟,羅發潘,鄭群瑞.1995.閩北天然闊葉林林層特征研究.福建林學院學報, 15(1):17-21.

(Huang Q L, Luo F P, Zheng Q R. 1995. Vertical structure of natural broad-leaved forest in North Fujian. Journal of Fujian College of Forestry, 15(1): 17-21. [in Chinese])

梁小筠. 1997. 正態性檢驗. 北京:中國統計出版社.

(Liang X J. 1997. Normality test. Beijing: China Statistics Press. [in Chinese])

劉福香. 2013. 大興安嶺地區天然落葉松-白樺混交林分樹種直徑分布研究. 哈爾濱:東北林業大學博士學位論文.

(Liu F X. 2013. Diameter distributions of individual species components of natural mixed forest ofLarixgmelinii-Betulaplatyphyllain Daxing’an Mountains. Harbin: PhD thesis of Northeast Forestry University. [in Chinese])

劉 健,陳平留,陳昌雄. 1996. 天然針闊混交林林分結構規律的研究.華東森林經理, 10(4): 1-4.

(Liu J, Chen P L, Chen C X. 1996. Study on the stand structure law of natural mixed stand of conifer and broad-leaved trees. East China Forest Management, 10(4): 1-4. [in Chinese])

孟憲宇. 2006. 測樹學.3版.北京:中國林業出版社.

(Meng X Y. 2006. Forest mensutation.3rd Edition. Beijing: China Forestry Publishing House. [in Chinese])

孟憲宇,張 弘. 1996. 閩北杉木人工林單木模型.北京林業大學學報, 18(2): 1-9.

(Meng X Y. Zhang H. 1996. An individual model for Chinese fir plantations in Fujian Province. Journal of Beijing Forestry University, 18(2): 1-9. [in Chinese])

孟憲宇. 1985. 使用Weibull分布對人工油松林直徑分布的研究.北京林學院學報, 7(1): 30-39.

(Meng X Y. 1985.A study on diameter distribution of some artificialPinustabulaeformisstands by using the Weibull function. Journal of Beijing Forestry College, 7(1): 30-39. [in Chinese])

孟憲宇.1991.削度方程和林分直徑結構在編制材種表中的重要意義.北京林業大學學報，13(2)：14-20.

(Meng X Y.1991.The significance of taper equations and stand diameter structures in development merchantable volume tables.Journal of Beijing Forestry University,13(2):14-20. [in Chinese])

彭少麟,王伯蓀.1983.鼎湖山森林群落分析Ⅰ.物種多樣性.生態科學, (1):11-17.

(Peng S L, Wang B S. 1983. Analysis on the forest communities of DinghushanⅠ. species diversity. Ecologic Science, (1): 11-17. [in Chinese])

徐化成. 2004. 森林生態與生態系統經營. 北京:化學工業出版社.

(Xu H C. 2004. Forest ecology and ecosystem management. Beijing: Chemical Industry Press. [in Chinese])

徐健君. 1999. 非參數核密度估計在岷江冷杉天然林直徑分布研宄中的應用.甘肅林業科技, 24(4): 22-26.

(Xu J J. 1999. An application of non-parameter nucle-densiety estimation on research of diameter distribution ofAbiesfaxoniananatural forest. Journal of Gansu Forestry Science and Technology, 24(4): 22-26. [in Chinese])

許彥紅,楊宇明,杜 凡,等. 2004. 西雙版納熱帶雨林林分直徑分布研究.西南林學院學報, 24(2):16-18.

(Xu Y H, Yang Y M, Du F,etal. 2004. A study on stand structure of the tropical rainforest in Xishuangbanna of Yunnan Province, Journal of Southwest Forestry College, 24(2): 16-18. [in Chinese])

張建國,段愛國,童書振. 2004. 林分直徑結構模擬與預測研究概述.林業科學研究, 17(6): 787-795.

(Zhang J G, Duan A G, Tong S Z. 2004. Review on the modeling and prediction of stand diameter structure. Forest Research, 17(6): 787-795. [in Chinese])

周國模,王瑞銓,俞雙群,等. 1994. 慶元縣人工杉木林直徑分布的研究.華東森林經理, 10(1): 18-32.

(Zhou G M, Wang R Q, Yu S Q,etal. 1994. Study on the diameter distribution of artificial chinese fir forest in Qingyuan County. East China Forest Management, 10(1): 18-32. [in Chinese])

中國植被編輯委員會. 1980.中國植被. 北京： 科學出版社.

(Chinese Vegetation Editorial Board. 1980. Chinese vegetation. Beijing: Science Press. [in Chinese])

莊崇洋. 2016. 中亞熱帶天然闊葉林林層特征研究.北京:中國林業科學研究院博士學位論文.

(Zhuang C Y. 2016. Study on the storeys characteristics of natural broad-leaved forests in middle-subtropical zone. Beijing: PhD thesis of Chinese Academy of Forestory. [in Chinese])

Bailey R L, Dell T R.1973.Qualifying diameter distributions with the Weibull function.Forest Science,19(2):97-104.

Bailey R L. 1980. Individual tree growth derived from diameter distribution models. Forest Science, 26(4): 626-632.

Borders B E, Patterson W D. 1990. Projecting stand tables: a comparison of the Weibull diameter distribution method, a percentile-based projection method and a basal area growth projection method. Forest Science, 36(2): 413-424.

Borders B E, Surfer R A, Bailey R L,etal. 1987.Percentile-based distributions characterize forest stand tables. Forest Science, 33(2): 570-576.

Burk T E, Newberry J D. 1984. A simple algorithm for moment-based recovery of Weibull distribution parameters. Forence Science, 30(2): 329-332.

Cao Q V. 2004. Predicting parameters of a Weibull function for modeling diameter distribution. Forest Science, 50(5): 682-685.

Clutter J L, Bennett F A. 1965. Diameter distributions in old-field slash pine plantations(Report 13). Macon:Report of Georgia Forest Research Council.

Clutter J L, Harms W R, Brister G H,etal. 1984. Stand structure and yields of site-prepared loblolly pine plantations in the lower coastal plain of the Carolinas, Georgia, and north Florida. USDA Forest Service, Southeastern Forest Experiment Station, Asheville, General Technical Report SE-27.

Dressier T D, Burk T E. 1989. A test of nonparametric smoothing diameter distribution stand. Scandinavian Journal of Forest Research, 4(1/4): 407-415.

Gingrich S F. 1967. Mesuring and evaluating stock and stand density in upland hardwood forest in the Central States. Forest Science, 13(1):38-53.

Goelz J C G, Leduc D J. 2002. A. model describing growth and development of longleaf pine plantations: consequences of observed stand structures of structure of the model (U.S. Department of Agriculture,Forest Service, Gen Tech Rep SRS-48),Asheville: Report of Southern Research Station.

Hafley A M, Schreuder H T. 1977. Statistical distributions for fitting diameter and height data in even-aged stands. Canadian Journal of Forest Research, 7(3): 481-487.

Hara A, Malt M, Tokola T. 1997. Thek-nearest neighbour method for estimating basa1 area diameter distribution. Scandinavian Journal of Forest Research, 12(2): 200-208.

Ishikawa Y. 1998. Analysis of the diameter distribution using the Richards distribution function,3: relationship between mean diameter or diameter variance and parameter mork of uniform and even-aged stands. Japanese Journal of Forest Planning, 31: 15-18.

Knoke T, Plusczyk N. 2001. On economic consequences of transformation of a spruce (Piceaabies(L.) Karst) dominated stand form regular into irregular age structure. Forest Ecology and Management, 151(1/3): 163-179.

Leak W B. 1964. An expression of diameter distribution for unbalanced, uneven-aged stands and forests. Forest Science, 10(1): 39-50.

Liu C, Zhang L, Davis C J,etal. 2002. A finite mixture model for characterizing the diameter distribution of mixed-species forest stands. Forest Science, 48(4): 653-661.

Maltamo M, Kangas A, Uuttera J,etal. 2000. Comparison of percentile based prediction methods and the Weibull distribution in describing the diameter distribution of heterogeneous scots pine stands. Forest Ecology and Management, 133(3): 263-274.

Meyer H A. 1952. Structure, growth, and drain in balanced unevenaged forests. Journal of Forestry, 50(2): 85-92.

Nelson D M. 1964. Diameter distribution and growth of loblolly pine. Forest Science, 10(1): 105-114.

(責任編輯 石紅青)

Diameter Distribution in Each Storey and Law of Typical Natural Broad-Leaved Forest in Mid-Subtropical Zone

Zhuang Chongyang1Huang Qinglin1Ma Zhibo1 2Zheng Qunrui3Wang Hong1

(1.KeyLaboratoryofForestryRemoteSensingandInformationTechnology,StateForestryAdministrationResearchInstituteofForestResourceInformationTechniques,CAFBeijing100091; 2.ResearchInstituteofForestEcology，EnvironmentandProtection,CAFBeijing100091; 3.WanmulinProvincialNatureReserveManagementOfficeofJian’ou，FujianProvinceNanping353105)

【Objective】 The diameter distribution and its changing law at storeys of typical natural broad-leaved forest in mid-subtropical zone was studied to provide the theoretical basis for revealing the structure, assisting the investigation of relevant stand description factors and directing the management planning of natural forests.【Method】Based on storey identification of maximum light receiving plane(MLRP) method, the Shapiro-Wilk(S-W test) test was used to measure the normality of the diameter distribution of storeys(including the whole stand and each storey); the Skewness(SK) and Kurtosis(KT) were used to describe the characters of the diameter distributions; the Meyer negative exponential function and the Weibull distribution function were served to fit the diameter distribution of each storey; the diameter distribution of each storey was charted by the model with better result of Chi-square test, and the changing law among storeys were analyzed.【Result】The result of S-W test showed that all the whole stand, the substorey Ⅲ and the substorey Ⅱ of the typical natural broad-leaved forest in mid-subtropical zone were not normal distribution. As to the substorey Ⅰ, the diameter distribution of sample plot 1, 2 and 3 were the normal distribution, but the substorey Ⅰ of sample plot 4 and 5 were oppo; within each sample plot, the value of S-W test increased with the increasing substorey height. In general, the absolute value of SK and KT decreased with the height of substorey increased, except the SK and KT of the substorey Ⅱ and Ⅲ in sample plot 5 and the KT of the substorey Ⅰ and Ⅱ in sample plot 2. The Meyer negative exponential function could not fit the diameter distribution of the whole stand(all the sample plots disobeyed the Chi-square test) and substorey Ⅱ(the sample plot 2 and 4 obeyed, and sample plot1, 3 and 5 disobey)well, but it could fit the substoreys Ⅰ(all sample plots obeyed) and Ⅲ(sample plot 3, 4 and 5 obeyed, and sample plot 1 and 2 disobey) well. The Weibull distribution function could fit the diameter distribution of storeys(including the whole stand storey and each substorey) better, except for the whole stand in sample plot 1, the surplus storeys in sample plot 1 and the storeys of the other sample plots all passed the Chi-square test. Used the Weibull distribution function fitted the diameter distribution of each storey in every sample plot, suggested that with the increase of storey height, the charts of storey changed from reverse “J” shape curve to right partial mountain curve and normal distribution curve, peak of curve moved to right, and the height of peak decreased.【Conclusion】 The Weibull distribution function could fit the diameter distribution of each storey better than the Meyer negative exponential function at storeys(including the whole satand and each substorey) of typical natural broad-leaved forest in mid-subtropical zone. The results of S-W test, SK, KT and the status of the peak of curve shape at the chat of Weibull function indicated that the diameter distribution of each substorey in each sample plot was trend to transit gradually to normal distribution with the mean DBH and mean height of each substorey increased. The differences among the storeys(including the whole stand storey and each substorey) illustrated that the necessity of studying diameter distribution within storey identification.

mid-subtropical zone; typical natural broad-leaved forest; storey identification; diameter distribution; change law

10.11707/j.1001-7488.20170403

2016-05-20；

2016-12-29。

國家自然科學基金項目“中亞熱帶天然闊葉林林層特征研究”(31370633)。

S757.2

A

1001-7488(2017)04-0018-10

*黃清麟為通訊作者。