研究区分别位于湖北省恩施州鹤峰县木林子国家级自然保护区、宣恩县七姊妹山国家级自然保护区和利川市团堡镇金子山国有林场。整个区域属典型亚热带季风性气候，夏季温湿多雨，冬季严寒，日照充足（表1）。研究区域内植物资源丰富，原生性强，区系组成复杂，具有保存较好的原生植被和受到人为干扰并处于不同恢复阶段的次生林^[29]。其主要林分类型为亚热带常绿落叶阔叶混交林。

2014 年6 月至2016 年8 月，参照CTFS（The Center for Tropical Forest Sciences）样地建设技术规范，使用全站仪先后在木林子、七姊妹山和金子山建立15 、6 和6 hm²大样地，并记录坡度、坡向、坡位、土壤、样地中心点经纬度和海拔等生境因子。样地内所有DBH ≥ 1 cm 的木本植物挂牌并进行每木检尺，记录树种名、胸径、坐标、树高、冠幅、萌生等林木特征值。

依据3个大样地基础调查数据，每个大样地以东西方向为X轴，南北方向为Y轴，左下角为坐标原点，依次向右上角以20 m为间距划分10块不同尺度样地，最大尺度为200 m × 200 m。

重要值( IV )是确定群落中优势种反映其作用和地位的指标。重要值( IV ) = 相对密度( RA ) + 相对频度( RF ) + 相对显著度( RP )。

基于林分空间结构单元，以参照树及其最近4株相邻木构成的结构单元选取角尺度、大小比数、混交度分析林分空间结构特征。具体计算公式见表2。

对物种丰富度、多样性以及均匀度的测定选用Simpson指数、Shannon-Weiner指数和Pielou指数计算（表3）。

筛选出DBH ≥ 5 cm的乔木，利用林木特征值计算乔木层的空间结构参数及多样性指数。采用林分空间结构分析软件Winkelmass^[30-32]计算核心区林木空间结构参数W、U和M，R 4.1.1计算物种多样性指数和物种重要值；利用SPSS 24、Excel进行数据分析。不同尺度下林分空间结构与物种多样性的关系使用R 4.1.1代码进行Pearson显著相关性检验。使用R 4.1.1和Origin 2017作图。

重要值代表了各树种在林分空间中的地位和优势程度。3个研究区内物种丰富，其中重要值排位前4的乔木树种如表4所示。七姊妹山自然保护区以川陕鹅耳枥（Carpinus fargesiana H. Winkl.）、多脉青冈（Cyclobalanopsis multinervis W.C.Cheng & T.Hong）为主；木林子自然保护区以多脉青冈、川陕鹅耳枥占优势种地位，石灰花楸（Sorbus folgneri (Schneid.) Rehd.）和小叶青冈（Cyclobalanopsis myrsinifolia (Blume) Oersted）为亚优势种。该群落中优势种和亚优势种的相对密度、重要值较其余研究区差距都较小，表明群落组成中这些树种的个体数量都较多。金子山国有林场中锥栗（Castanea henryi (Skan) Rehd. et Wils.）的重要值（22.38%），明显高于其他乔木物种，其次为川陕鹅耳枥（8.24%）。金子山林分内有较多的锥栗大径级木，所以相对显著度高达44.88%。总的来看，3个研究区主要优势树种有一定差异，但是川陕鹅耳枥的重要值在3个研究区群落中都较高，这说明该种群在鄂西南常绿落叶阔叶混交林群落构建中起重要作用，群落中常见优势乔木树种为川陕鹅耳枥、多脉青冈和锥栗。

在不考虑小尺度样地（20 m × 20 m）出现极个别特殊数值的情况下，不同尺度上，3个研究区内林分多样性指数与尺度之间呈现出明显的非线性相关（图1），整体上Shannon指数和Simpson指数随着尺度增大呈现先快速增加后趋于缓慢增加或稳定的趋势，Pielou指数则随着尺度增加逐步下降然后趋于稳定，但在各个研究区呈现出一定差异。随着空间尺度的增加，七姊妹山多样性指数增幅逐渐趋缓；金子山波动变化较大，但整体趋势与七姊妹山接近；木林子多样性增幅逐渐趋于零，且在大尺度上（≥ 160 m × 160 m），多样性指数值显著低于其余2个研究区。金子山多样性指数在尺度40 m × 40 m、120 m × 120 m时骤降，这可能是由于金子山部分样地受到一定程度的人为干扰，林分内物种受到了损害，从而导致多样性指数的明显下降。鄂西南3个研究区的多样性指标在不同尺度上的变异系数整体呈现先减小后增大的趋势，其中Shannon-weiner指数变异系数在不同尺度上波动最剧烈，除在40 m × 40 m尺度时的Shannon-weiner指数，其余尺度不同指数都低于8%。Simpson指数和Pielou指数从尺度20 m × 20 m上升至140 m × 140 m，变异系数逐渐降低，随后达到160 m × 160 m以上，变异系数又出现上升趋势，且随着空间尺度增大，变异系数也随之增大。

图  1  不同尺度多样性指数变化

Figure 1.  Variations of diversity indices at different scales

七姊妹山自然保护区、木林子自然保护区和金子山国有林场在不同尺度上林分空间结构参数W、U、M变化如图2所示。从局部看，木林子角尺度曲线整体低于其他研究区，角尺度均值表现为七姊妹山（0.557）＞ 金子山（0.551）＞ 木林子（0.533），且各研究区角尺度均值达到0.536（> 0.517）以上，水平空间分布格局多属于轻微聚集分布，而木林子林木分布格局聚集程度最低。林分大小比数均值在0.500～0.510之间，随着尺度增加，各研究区差异变化不明显，均呈现中庸状态。说明各研究区林木大小分化处于中庸状态，空间结构单元中相邻木与参照树胸径大小差异化不明显，不同优势度树种数量相近，林分内树种生长竞争较为激烈。林分空间隔离程度在3个研究区有较明显差异，各研究区林分混交度均值达0.732以上，其中木林子整体处于极强度混交状态（平均混交度0.757），这表明研究区多数参照木周围相邻木不为同种，单种聚集发生概率低，林木整体混交良好。鄂西南3个研究区的林分空间结构参数在不同尺度上的变异系数变化，在20 m × 20 m至100 m × 100 m尺度范围内，结构参数指标变异系数波动较大（1.39%～10.82%），随着尺度上升至120 m × 120 m尺度，结构参数指标变异系数分别为：W（2.08%）、M（2.93%）和U（0.40%），逐渐趋于稳定。因此鄂西南常绿落叶阔叶混交林林分空间结构研究最小样地尺度范围为120 m × 120 m，该尺度能较全面的反映鄂西南林分空间结构变化规律。

图  2  不同尺度空间结构参数变化

Figure 2.  Variations of spatial structure parameters at different scales

由图3可知，在鄂西南常绿落叶阔叶混交林中，不同尺度上角尺度频率分布取值集中在0.5的林木个体数量超过了55%，处于0.5右侧的取值频率高于0.5左侧的频率值，且W = 0（非常均匀）的频率值均低于1%，表明林分中处于随机分布状态的林木个体数量最多，总体聚集程度不高，因此鄂西南常绿落叶阔叶混交林林分的水平分布格局在各尺度上呈现出微聚集分布状态，且随着尺度增加，分布状态未呈现明显变化。结合林分不同树种的空间结构单元的平均角尺度分析发现，主要优势树种的平均角尺度与整体林分角尺度接近，这进一步说明优势种对于整个林分林木的分布中起着主导作用。

图  3  不同尺度林分角尺度(W)一元分布

Figure 3.  The univariate distribution of stand uniform angle at different scales

3个研究区的结构参数取值频率分布表现出与尺度的极大相关性，即随着尺度的增加（40 m × 40 m → 200 m × 200 m），各研究区5个取值的频率值差异逐渐减小并趋于相等。例如：W = 0.25时，小尺度（40 m × 40 m）中，3个研究区频率值变异系数为28.4%，而在在大尺度（200 m × 200 m）中，3个研究区频率值变异系数为6.6%。这种现象同样适用于大小比数、混交度不同取值频率分布。可见，尺度观测增加到一定的范围时，鄂西南各研究区林分空间结构基本趋于一致。

参照树与4株相邻木构成的最佳空间结构单元中，描述参照树与相邻木的相对优势程度，各优势程度林木占全部林木的比例，即大小比数的一元分布。图4可见，不同研究区不同尺度的大小比数频率值都在20%左右，且随着尺度增加到120 m × 120 m时，各大小比数频率趋于稳定，表明鄂西南常绿落叶阔叶混交林中不同优势度的树种数量大致相等，群落内竞争处于中庸状态，且随着尺度上升更加稳定。

图  4  不同尺度林分大小比数(U)一元分布

Figure 4.  The univariate distribution of stand dominance at different scales

混交度说明了结构单元中的树种空间隔离程度，以是否为同种树种计算出5个取值等级，每个取值等级在林分中所占比例，即为混交度一元分布。如图5所示，鄂西南3个研究区不同尺度上林木混交度频率分布主要集中在0.75（强度混交）和1（极强度混交）状态，频率和均在64%以上，零度混交和弱度混交所占比例较少，表明3个研究区林分呈现较强的混交状态，树种组成丰富。

图  5  不同尺度林分混交度(M)一元分布

Figure 5.  The univariate distribution of Mingling at different scales

以七姊妹山、木林子、金子山3个大样地调查数据为基础，将不同尺度林分空间结构参数与乔木物种多样性指标进行相关性分析（图6）。亚热带常绿落叶阔叶混交林在不同尺度上，混交度和大小比数与物种多样性不存在显著相关性（P > 0.05），角尺度与乔木Pielou指数存在中度负相关性（P < 0.01），角尺度与Shannon-Weiner指数存在中度正相关性（P < 0.01），角尺度与Simpson指数存在低度相关性（P < 0.05）。综合来看，不同尺度变化中角尺度与乔木物种多样性存在显著相关性，即角尺度是影响鄂西南亚热带常绿落叶阔叶混交林林分乔木多样性的一个重要因子，随着林分角尺度的增大，林分内物种多样性程度越高，群落也越复杂，林木分布越趋向于不均匀。

图  6  林分空间结构指标与乔木物种多样性指数之间的相关性

Figure 6.  Correlation between stand spatial structure indices and tree species diversity index

本研究基于鄂西南3个研究区大样地清查数据中DBH ≥ 5 cm的乔木林木属性和位置信息，分析林分空间结构指标和乔木物种多样性指数随着尺度的变化趋势与相关性。研究结果一定程度上弥补了以往对大尺度林分空间结构研究及其与多样性相互关系的空缺^[33-34]，与近几年相关研究对比发现：

通过对鄂西南3个研究区常绿落叶阔叶混交林乔木层物种多样性指标进行分析，发现各研究区物种多样性较高，且不同区域优势树种的重要值存在一定的相似性，因此鄂西南常绿落叶阔叶混交林常见优势乔木树种为川陕鹅耳枥、多脉青冈和锥栗等，这与相关学者研究结果基本一致^[28,35]。同时研究尺度的大小影响着物种空间分布格局，合适的样地尺度对于生物多样性和林分空间结构研究至关重要^[36]。鄂西南3个研究区物种多样性随着研究尺度的上升，多样性与尺度之间呈现出明显的非线性相关。金子山Shannon指数和Simpson指数在40 m × 40 m、120 m × 120 m尺度上骤降，主要是由于部分样地受到了较强的人为干扰。随着尺度增加到120 m × 120 m以上时，3个研究区的多样性指数波动幅度也逐渐减小，这与He F等人^[37]研究结果一致，表明小尺度范围内样地空间异质性表达较为明显，尺度增加带来更加多样的生境，但新生境出现速度也会减慢，异质性逐渐平均化，因此大尺度上很少出现剧烈波动。鄂西南作为喀斯特地貌发育最强烈的地区，创造了丰富多样的生境，使鄂西南地区森林生态系统维持着较高水平的物种多样性、异质性和复杂性^[21]，因此鄂西南七姊妹山、金子山、木林子有着鲜明的多样性差异，同时受限于目前的研究尺度，关于物种多样性最小面积还需进一步增加尺度以进行确定。

鄂西南地区不同尺度林分空间特征表现出微聚集分布、大小分化中庸且混交良好的状态，与薛卫星^[38-39]等的研究结果相一致。林分中主要优势树种的结构参数均值与整体林分均值接近，表明主要优势树种的林木分布格局很大程度上决定了林分空间特征^[40] 。同时3个空间结构参数存在区域性差异，木林子自然保护区处于微聚集、极强度混交状态，其余研究区林木聚集程度高于木林子，且呈强度混交状态。一般认为，天然林在演替初期会呈明显的聚集分布，这是由种源扩散和林中开放空间的形式决定的，随着演替的推进，林木个体对资源的需求导致了生存竞争，这种自然稀疏的过程促使林分朝着随机分布格局演替^[41]。同时聚集分布格局相较于随机分布格局，总体混交度更低^{[18, 42]}。因此木林子虽然聚集程度相对较低，但由于各树种种内排斥、同树种种间分布较均匀，因此林分整体混交度较高，演替阶段也高于七姊妹山和金子山。随着样地尺度的变化，其对应林木的空间分布格局也不同，这表明，取样面积应有一定的大小才能真实地反映林分特征^[42]。研究显示能够稳定反映鄂西南亚热带常绿落叶阔叶混交林林分空间结构变化规律的最小样地尺度为120 m × 120 m（1.44 hm²），而胡中洋等^[19]对新疆天山国家级自然保护区确立过80 m × 80 m（0.64 hm²）最小样地尺度，与本研究有一定差异。这可能是由于两个区域气候、地形、生境等有明显空间异质性，导致不同植被类型的样地尺度有所差别，并且鄂西南地区生物多样性非常丰富^{[19, 43]}，林分结构在尺度上表现出更复杂的特征，因此120 m × 120 m的尺度更适合鄂西南地区的空间结构样地调查。

林分空间结构与林分的物种多样性存在一定相关性^[4]。同一林分类型不同演替阶段中，多种空间结构指标将综合影响物种多样性，其中混交度、角尺度为关键因子，但各个演替阶段的影响又各不相同。鄂西南3个区域不同尺度上林分空间结构参数与乔木物种多样性指标存在一定的相关性，且仅表现为角尺度与林分多样性存在中度相关性，其他结构参数与林分多样性则不存在相关性，而黎芳等^[4]研究的飞播马尾松人工林，张亚昊等^[44]研究的不同演替阶段马尾松林则表现出空间结构指标如角尺度、混交度^[45]等对林分物种多样性有显著影响。导致研究结果差异首先是因为鄂西南常绿落叶阔叶混交林物种组成丰富，且处于基本无人为干扰状态，因此空间结构主要通过自然竞争、稀疏等形成，整体呈现结构稳定，多样性丰富维持的状态。此外，由于本研究的区域尺度较大，乔木物种分布格局影响因子如：地形、温度、光照、水分等发生异质性变化，物种多样性因此也呈现出相应变化，表现出一定的相关性。一般认为小尺度上，环境差异性较小，植物群落分布主要受种子扩散限制性影响^[46]；大尺度上植物群落分布更多是受到生境异质性的作用，尤其是喀斯特地区复杂的地形地貌，故在不同的尺度上反映出不同的生态学特征^[47]，这也成为影响物种多样性差异的关键因素之一^[47]。结合本研究的结果，大尺度上的乔木物种多样性变化不仅受到生境异质性的作用，还受到空间结构参数角尺度的影响。

鄂西南地区亚热带常绿落叶阔叶混交林中林分空间结构特征表现出微聚集分布、大小分化中庸且混交良好的状态，该区进行林分空间结构调查的最适尺度为120 m × 120 m。不同尺度上林分空间结构参数角尺度与林分物种多样性存在中度相关性，其他结构参数与林分物种多样性不存在相关性，可见，在天然林中不同研究尺度下角尺度是林分空间结构影响物种多样性的一个关键因子。在今后鄂西南亚热带常绿落叶阔叶混交林可持续发展过程中，建议考虑适当调整林分空间结构，优化角尺度，利于该区乔木物种多样性的提高，以进一步促进该区森林生态系统的稳定发展。