典型中亚热带天然阔叶林各林层直径分布及其变化规律*

庄崇洋 黄清麟 马志波，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

*黄清麟为通讯作者。