论文部分内容阅读
摘 要：該文以湖北木林子国家级自然保护区15 hm2大样地同属物种 [木姜子（Litsea pungens）和黄丹木姜子（L. elongata）]为研究对象，通过成对相关函数和不同零模型（完全空间随机模型、异质泊松模型和先决条件模型）的方法分析两个物种的空间分布格局、种内和种间空间关联性，从同属物种的空间分布格局角度探讨群落的构建机制。结果表明：（1）种群空间分布格局研究中，两个种的空间分布格局大致为聚集分布。基于CSR零模型，木姜子和不同径级（幼树、小树、中树和成年树）的聚集尺度分别为0～114 m、0～103 m、0～56 m、0～34 m和0～16 m，黄丹木姜子和不同径级的聚集尺度分别为0～150 m、0～150 m、0～59 m、0～56 m和0～14 m;基于HP零模型，木姜子整体的空间分布格局聚集尺度为0～23 m，其中幼树、小树和中树均在0～22 m呈现聚集分布，成年树仅在0～3 m尺度上表现为聚集分布，黄丹木姜子整体的空间分布格局聚集尺度为0～36 m、140～150 m，各径级的聚集尺度分别为0～35 m、0～35 m、0～12 m和0～6 m。（2）种内关联性研究中，木姜子和黄丹木姜子不同径级个体种内关联在小尺度上基本呈现正关联，随着尺度的增大变为微弱的无关联或负关联。（3）种间关联研究中，木姜子对黄丹木姜子和黄丹木姜子对木姜子的空间关系大致相同。总体上，两物种之间均体现为在30 m尺度内对另外一个物种呈现正关联，40～68 m尺度内呈现负关联;不同径级之间都表现为无关联，偶尔出现正关联或负关联。研究认为，种群空间分布格局基本上为聚集分布，其主要受到生境异质性和扩散限制的影响，木林子大样地中木姜子和黄丹木姜子空间分布格局存在相似性，而不同的生境偏好可能是导致两同属物种共存的原因。
关键词：同属物种，点格局，空间分布格局，空间关联，物种共存
中图分类号：Q948.15
文献标识码：A
文章编号：1000-3142（2021）08-1333-13
Abstract： In this paper， species of the same congeneric species （Litsea pungens and L. elongata ） in large sample plot of 15 hm2 in Mulinzi National Nature Reserve of Hubei Province were taken as the study object. The spatial distribution pattern， intraspecies and interspecific spatial correlations of the two species were analyzed by pairwise correlation function and different zero models（complete spatial randomness， heterogeneous Poisson and antecedent condition）， and the construction mechanism of the community was discussed from the point of view of the spatial distribution pattern of the same congeneric species. The results were as follows： （1） In the study of population spatial distribution pattern， the spatial distribution pattern of the two species was generally aggregation distribution. The aggregation scales on the basis of the CSR zero model， the L. pungens and the different diameter stages（sapling， small tree， medium tree， adult tree） were 0-114 m， 0-103 m， 0-56 m， 0-34 m， and 0-16 m， respectively， and the aggregation scales of the L. elongata and the different diameter stages were 0-150 m， 0-150 m and 0-59 m， 0-56 m and 0-14 m; On the basis of the HP-zero model， the spatial distribution pattern of the whole L. pungens was 0-23 m， where in all of the sapling， small and medium trees were 0-22 m， and the aggregation distribution of adult trees was only on the scale of 0-3 m. The aggregation scale of the whole spatial distribution pattern of L. elongata was 0-36 m and 140-150 m， and the aggregation scale of each diameter grade was 0-35 m， 0-35 m， 0-12 m and 0-6 m. （2） In the study of intraspecific correlation， the intraspecies correlation of different diameter individuals showed positive correlation on a small scale， and became weakly unrelated or negative with the increase of scale. （3） In the study of interspecific correlation， the spatial relationship between the two species was about the same. On the whole， the two species showed positive correlation with another species within 30 m scale and negative correlation within 40-68 m scale， and there was no correlation between different diameter grades， and there was positive or negative correlation between the two species. It is considered that the spatial distribution pattern of the population is basically aggregation distribution， which is mainly affected by habitat heterogeneity and diffusion. The spatial distribution patterns of L. pungens and L. elongata in large sample plot are similar， and different habitat preferences may be the cause of the coexistence of the two species. Key words： same congeneric species， point pattern， spatial distribution pattern， spatial correlation， species coexistence
一般来说，同属物种拥有相同祖先，在形态学和生理特征上都具有一定的相似性，所需的资源可能相似，从生态位分化来看，同属物种可能因为这些相似性而产生竞争，而不能稳定共存（Ackerly & Donoghue， 1998; Webb & Ackerly， 2002），同属物种共存现象多发现于热带和温带植物群落中（张健， 2009）。因此，生态学家一直对同属植物的共存进行激烈的探讨。在生态研究中对同属物种进行比较，可以最大限度地减少系统差异造成的混乱影响，从同属物种的空间结构和相关性的角度进行探讨（Swenson et al.， 2006; 姚兰， 2016）。
植物种群的空间分布格局是指某一种群在一定水平内所有个体的分布位置或配置（曲仲湘， 1984; 张金屯， 1995），到现在为止是生态学研究的热点之一（Wiegleb， 1989; Nathan， 2006）。点格局方法由Ripley提出，它能够很好地解决传统格局分析方法的缺点（Condit， 2000; 张金屯和孟东平， 2004）。近些年，点格局分析越来越广泛，主要侧重于单物种的空间分布格局（范繁荣等， 2008; 黄小波等， 2018）、优势树种的空间分布格局（赵中华等， 2011; 张玉双等， 2016）、两种或两种以上物种的空间相关性（姚良锦等， 2018）以及空间分布格局与环境之间关系（刘旻霞等， 2018; 闫海冰， 2018）等研究较多，对于种群同科（何增丽等， 2017）或同属（陈俊等， 2018）研究较少。
樟科（Lauraceae）木姜子属（Litsea Lam） 的植物分布非常广泛，在中国大约有72种（陈幼竹， 2004）。亚热带常绿落叶阔叶混交林保存最为完好的地区包括木林子国家级自然保护区（黄永涛， 2016），其中木姜子属植物仅有2种：黄丹木姜子和木姜子，其数量占样地内物种总数量的4.8%。木姜子的药用价值很高，黄丹木姜子具有较高的商业价值（陈幼竹， 2004）。木姜子和黄丹木姜子的研究仅仅只限于药用价值等方面的研究，黄丹木姜子还见于叶凋落量及养分特征研究（董舒， 2016），对于更多的生态学方面的研究未见有报道，黄丹木姜子和木姜子虽不是建群种和先锋种，这类小乔木在种群空间分布格局的研究同样具有重要意义。本研究以湖北木林子大样地中同属的木姜子和黄丹木姜子为研究对象，采用g（r）函数对木姜子属两个种进行空间格局及空间关联分析，探讨其成因，探讨两物种是否存在相似性，同属物种共存机制，并为研究亚热带地区物种共存机制提供参考。
1 研究区概况
研究区（109°59′30″—110°17′58″ E、29°55′59″—30°10′47″ N）位于湖北省鹤峰县，面积为20 838 hm2，属于云贵高原的东北延伸区。该地区最低海拔为1 100 m，最高海拔为2 095.6 m，大陆性季风气候，冬季无严寒，夏季无极端炎热，雨热同季，日照充足，山地性气候明显，年均温15.5 ℃，年均降水量1 700～1 900 mm，年均相对湿度82%，无霜期270～279 d。土壤为黄棕壤、棕壤和黄壤。
木林子15 hm2大样地位于保护区核心区内，海拔为1 588～1 780 m，地势相对缓慢，属于亚热带山地常绿落叶阔叶混交林，主要树种有翅柃（Eurya alata）、多脉青冈（Cyclobalanopsis multinervis）、川陕鹅耳枥（C. fargesiana）、小叶青冈（C. myrsinifolia）、石灰花楸（Sorbus folgneri）、四照花（Dendrobenthamia japonica var. chinensis）、锥栗（Castanea henryi）、白栎（Quercus fabri）、山矾（Symplocos sumuntia）和光叶山矾（S. lancifolia）等（姚兰， 2016）。
2 研究方法
2.1 样地建设及野外调查
2013年，按照CTFS大样地建设与检测技术规范，建立15 hm2固定监测样地，分成375块20 m × 20 m小樣地，用防腐处理的不锈钢管将4角作为永久标识，每个小样地再分为4个10 m × 10 m样方，每个样方再分为4个5 m × 5 m的小样方，对DBH≥1.0 cm的木本植物进行挂牌，并用红色油漆标记（树高1.3 m），对DBH≥1 cm的所有木本植物进行每木检尺，识别物种（科、属、种），调查所有木本植物在样地内的坐标值（以大样地单元西南角为原点，测定x轴和y轴的坐标值）。
2.2 种群径级划分
运用径级代替龄级的方法，结合湖北木林子样地内木姜子和黄丹木姜子的径级特征，将胸径DBH划分为4个径级：Ⅰ级（1 cm≤DBH<2.5 cm），Ⅱ级（2.5 cm≤DBH<5 cm），Ⅲ级（5 cm≤DBH<10 cm），Ⅳ级（DBH≥10 cm），分别为幼树、小树、中树和成年树。
2.3 点格局分析方法
本研究采用点格局分析中成对相关函数g（r）进行分析。g（r）函数是基于对所有定位的成对个体之间距离的格局分析，消除了Ripiey’s K（r）函数在大尺度上产生的累积效应（Ripley， 1976; Thorsten et al.， 2007）。分析时，选取距离尺度为较短边长的1/2（150 m），圆环宽度为1 m，通过199次Monte Carlo随机模拟，利用模拟的最大值和最小值绘制上下包迹线，计算99%的置信区间（Wiegand et al.， 2006， 2009）。 2.4 零模型选择
采用点格局分析，需要检测生境异质性在物种分布中是否存在影响。一般来说，如果在大于10 m尺度上呈聚集分布，那就可以解释为受到生境异质性的影响;采用异质泊松模型和先决条件模型目的在于消除大尺度上生境异质性对空间分布格局的影响，在排除生境异质性后，如果小尺度的格局还是聚集，则说明扩散限制存在作用。所以在进行点格局分析时，不同的分析需要选择不同的零假设模型（null model），选取正确合适的零模型十分重要（梁爽等， 2014）。在本研究中单变量g（r）函数采用完全空间随机模型（complete spatial randomness， CSR）和异质泊松模型（heterogeneous Poisson， HP），双变量g（r）函数选择先决条件模型（antecedent condition， AC）。
文中数据的处理与分析使用R 3.4.1软件中的“spatstat”包中完成点格局分析，使用Sigmaplot 13软件绘制空间分布点，Origin 2018软件绘制点格局分析图。
3 结果与分析
3.1 两个物种空间分布格局
整个样地中（表1），木姜子1 701株，黄丹木姜子2 306株。由图1可知，两个物种的空间分布总体差异较大。木姜子主要分布于低海拔、山脊地区;黄丹木姜子主要分布于高海拔和山谷两侧区域，小部分零星分布于样地内。两个物种的整体空间分布格局的研究（图2）表明，基于CSR零模型：木姜子和黄丹木姜子的聚集尺度分别为0～114 m和0～150 m。基于HP零模型：木姜子的聚集尺度为0～23 m，其他尺度上为微弱的随机或均匀分布;黄丹木姜子的聚集尺度为0～36 m、140～150 m，随机分布尺度为48～139 m。
根据两个物种不同径级的空间分布格局研究（图3，图4）表明，基于CSR零模型：木姜子幼树、小树和中树的聚集尺度分别为0～103 m、0～56 m和0～34 m，成年树在0～16 m尺度表现为聚集分布;黄丹木姜子幼树、小树、中树和成年树分别在0～150 m、0～59 m、0～56 m和0～14 m尺度上均表现为聚集分布。基于HP零模型：木姜子幼树、小树和中树的聚集尺度均为0～22 m，而成年树仅在0～3 m尺度上表现为聚集分布;黄丹木姜子各径级的聚集尺度分别为0～35 m、0～35 m、0～12 m和0～6 m。由此可见，随着径级的增长，木姜子与黄丹木姜子均呈现聚集尺度逐渐减小的现象，特别是成年树仅在小尺度上表现为聚集分布，而随机分布格局现象比较明显。
3.2 两个物种不同径级种内关联性
木姜子不同径级种内关联见图5。基于AC零模型，成年树对幼树在所有尺度上均表现为无关联，而其他不同径级之间在较小尺度（r≤20 m）上表现为正相关，其中小树对幼树、中树对幼树、中树对小树、成年树对小树、成年树对中树正关联尺度分别为0～20 m、0～18 m、0～20 m、0～14 m和0～17 m。
黄丹木姜子不同径级种内关联见图6。基于AC零模型，黄丹木姜子成年树对中树在所有尺度上表现为无相关，小树对幼树、中树对幼树、成年树对幼树、中树对小树、成年树对小树分别在0～36 m、0～25 m、8～15 m、0～23 m和0～6 m尺度上呈现正关联。
3.3 两个物种种间关联性
用先决条件模型排除30 m以上生境异质性，结果显示：整体上，木姜子对黄丹木姜子在0～30 m呈现正关联，黄丹木姜子对木姜子在7～29 m呈现正关联（图7）。不同径级之间，木姜子不同径级对黄丹木姜子不同径级空间关联结果显示（图8）：木姜子幼树对黄丹木姜子幼树在4～34 m尺度表现为正关联;木姜子小树对黄丹木姜子幼树在13～29 m尺度表现为正关联，42～67 m尺度上为负关联;木姜子中树对黄丹木姜子幼树在40～50 m尺度上为负关联;木姜子小树对黄丹木姜子小树在0～25 m尺度上为正关联，44～60 m尺度上为负关联;木姜子中树对黄丹木姜子中树在0～5 m尺度上为正关联，在47～54 m和119～128 m尺度上表现为负关联;其余径级之间均呈现无关联。黄丹木姜子不同径级对木姜子不同径级空间结果显示（图9）：黄丹木姜子幼树对木姜子幼树在13～32 m尺度呈现正关联，42～76 m尺度呈现负关联;黄丹木姜子小树对木姜子幼树在0～13 m尺度呈现正关联;黄丹木姜子小树对木姜子小树在0～15 m呈现正关联，46～59 m尺度呈现负关联;其余径级之间均呈现无关联。
由此可见，木姜子对黄丹木姜子空间关系和黄丹木姜子对木姜子空间关系大致相同，仅有部分细微差别，其中两物种成年树与另一物种所有径级均没有显著的生态学关系。
4 讨论
4.1 两个物种空间分布格局形成的原因
植物自身生物学特性和环境因素共同作用于种群的空间分布格局，同一树种在不同的发育阶段将呈現不同的空间分布模式（Frost & Rydin， 2000）。本研究中基于CSR和HP零模型两个物种在研究尺度上主要表现为聚集分布，随着尺度增加表现为随机分布或均匀分布，结果与大部分研究一致（祝燕等， 2011; Guo et al.， 2013）。在很大程度上，生境的异质性被认为是植被空间格局的主要原因。通常，如果它的聚集尺度大于10 m，则可以解释为受到栖息地异质性的影响（Takyu et al.， 2002; Getzin et al.， 2008）。基于CSR零模型，在整体上木姜子和黄丹木姜子分别在0～114 m和0～150 m尺度范围内呈聚集分布，并使用HP零模型排除了30 m以上的生境异质性，聚集尺度减小成为随机分布，可以推测，两个种群空间分布格局受到生境异质性影响。从幼树、小树、中树、成年树时间序列来看，分布格局从聚集分布向随机分布发展趋势，符合一般规律（Vargas et al.， 2010; 张玉双等， 2011）。其中木姜子和黄丹木姜子不同径级之间均体现为幼树聚集尺度>小树聚集尺度>中树聚集尺度>成年树聚集尺度，可见生境异质性对两个物种的幼树影响最强，可能是在此阶段，种群的个体较小，生存力较低，需要通过聚集分布形成相互庇护，能够更好地维持种群自身的稳定性（张金屯， 1998; Wang et al.， 2014）。生境异质性对两个物种的成年树影响较小，两个物种负密度制约超过了生境异质性对聚集分布的影响，随着径级的增大，对资源、对空间的需求也增加，物种内部的竞争日益激烈，物种内部的稀疏性越来越强，这将导致一些个体死亡从而形成均匀或随机分布的状态（舒兰等， 2019）。 YAO LJ， YAO L， YI YM， et al.， 2018. Spatial patterns of dominant species Carpinus fargesiana and Cyclobalanopsis multinervis in a subtropical evergreen and deciduous broad-leaved mixed forest [J]. Sci Silv Sin， 54（12）： 1-11. [姚良锦， 姚兰， 易咏梅， 等， 2018. 亚热带常绿落叶阔叶混交林优势种川陕鹅耳枥和多脉青冈的空间格局 [J]. 林业科学， 54（12）： 1-11.]
ZHANG JT， 1995.Quantitative plant ecology [M]. Beijing： Science and Technology of China Press. [张金屯， 1995. 植被数量生态学方法 [M]. 北京： 中国科学技术出版社.]
ZHANG JT， 1998. Analysis of spatial point pattern for plant species [J]. Chin J Plant Ecol， 22（4）： 344-349. [张金屯， 1998. 植物种群空间分布的点格局分析 [J]. 植物生态学报， 22（4）： 344-349.]
ZHANG JT， MENG DP， 2004. Spatial pattern analysis of individuals in different age-classes of Larix principis-rupprechtii in Luya mountain reserve， Shanxi， China [J]. Acta Ecol Sin，（1）： 35-40. [张金屯， 孟东平， 2004. 芦芽山华北落叶松林不同龄级立木的点格局分析 [J]. 生态学报， （1）： 35-40.]
ZHANG YS， DUAN WB， DU S， et al.， 2016. Analysis on spatial distribution pattern for main populations and gap makers in Korean pine broad-leaved forest in Xiaoxing’anling Mountains of Northeast China [J]. For Res， 29（3）： 395-401. [張玉双， 段文标， 杜珊， 等， 2016. 阔叶红松林主要种群及林隙形成木的空间格局分析 [J]. 林业科学研究， 29（3）： 395-401.]
ZHAO ZH， HUI GY， HU YB， et al.， 2011. Spatial distribution patterns and associations of dominate population in two types of broadleaved Korean pine mixed forest [J]. For Res， 24（5）： 554-562. [赵中华， 惠刚盈， 胡艳波， 等， 2011. 2种类型阔叶红松林优势种群空间分布格局及其关联性 [J]. 林业科学研究， 24（5）： 554-562.]
ZHU Y， BAI F， LIU HF， et al.， 2001. Population distribution patterns and interspecific spatial associations in warm temperate secondary forests， Beijing [M]. Biodivers Sci， 19（2）： 252-259. [祝燕， 白帆， 刘海丰， 等， 2011. 北京暖温带次生林种群分布格局与种间空间关联性 [J]. 生物多样性， 19（2）： 252-259.]
ZONG G， BAI XJ， ZHANG SY， et al.， 2018. Spatial pattern and interspecific spatial association of tree seedlings in a secondary forest in montane region of eastern Liaoning Province， China [J]. Chin J Appl Ecol， 29（1）： 18-24. [宗国， 白雪娇， 张淑媛， 等， 2018. 辽东山区次生林乔木幼苗分布格局与种间空间关联性 [J]. 应用生态学报， 29（1）： 18-24.]
（责任编辑 周翠鸣）
关键词：同属物种，点格局，空间分布格局，空间关联，物种共存
中图分类号：Q948.15
文献标识码：A
文章编号：1000-3142（2021）08-1333-13
Abstract： In this paper， species of the same congeneric species （Litsea pungens and L. elongata ） in large sample plot of 15 hm2 in Mulinzi National Nature Reserve of Hubei Province were taken as the study object. The spatial distribution pattern， intraspecies and interspecific spatial correlations of the two species were analyzed by pairwise correlation function and different zero models（complete spatial randomness， heterogeneous Poisson and antecedent condition）， and the construction mechanism of the community was discussed from the point of view of the spatial distribution pattern of the same congeneric species. The results were as follows： （1） In the study of population spatial distribution pattern， the spatial distribution pattern of the two species was generally aggregation distribution. The aggregation scales on the basis of the CSR zero model， the L. pungens and the different diameter stages（sapling， small tree， medium tree， adult tree） were 0-114 m， 0-103 m， 0-56 m， 0-34 m， and 0-16 m， respectively， and the aggregation scales of the L. elongata and the different diameter stages were 0-150 m， 0-150 m and 0-59 m， 0-56 m and 0-14 m; On the basis of the HP-zero model， the spatial distribution pattern of the whole L. pungens was 0-23 m， where in all of the sapling， small and medium trees were 0-22 m， and the aggregation distribution of adult trees was only on the scale of 0-3 m. The aggregation scale of the whole spatial distribution pattern of L. elongata was 0-36 m and 140-150 m， and the aggregation scale of each diameter grade was 0-35 m， 0-35 m， 0-12 m and 0-6 m. （2） In the study of intraspecific correlation， the intraspecies correlation of different diameter individuals showed positive correlation on a small scale， and became weakly unrelated or negative with the increase of scale. （3） In the study of interspecific correlation， the spatial relationship between the two species was about the same. On the whole， the two species showed positive correlation with another species within 30 m scale and negative correlation within 40-68 m scale， and there was no correlation between different diameter grades， and there was positive or negative correlation between the two species. It is considered that the spatial distribution pattern of the population is basically aggregation distribution， which is mainly affected by habitat heterogeneity and diffusion. The spatial distribution patterns of L. pungens and L. elongata in large sample plot are similar， and different habitat preferences may be the cause of the coexistence of the two species. Key words： same congeneric species， point pattern， spatial distribution pattern， spatial correlation， species coexistence
一般来说，同属物种拥有相同祖先，在形态学和生理特征上都具有一定的相似性，所需的资源可能相似，从生态位分化来看，同属物种可能因为这些相似性而产生竞争，而不能稳定共存（Ackerly & Donoghue， 1998; Webb & Ackerly， 2002），同属物种共存现象多发现于热带和温带植物群落中（张健， 2009）。因此，生态学家一直对同属植物的共存进行激烈的探讨。在生态研究中对同属物种进行比较，可以最大限度地减少系统差异造成的混乱影响，从同属物种的空间结构和相关性的角度进行探讨（Swenson et al.， 2006; 姚兰， 2016）。
植物种群的空间分布格局是指某一种群在一定水平内所有个体的分布位置或配置（曲仲湘， 1984; 张金屯， 1995），到现在为止是生态学研究的热点之一（Wiegleb， 1989; Nathan， 2006）。点格局方法由Ripley提出，它能够很好地解决传统格局分析方法的缺点（Condit， 2000; 张金屯和孟东平， 2004）。近些年，点格局分析越来越广泛，主要侧重于单物种的空间分布格局（范繁荣等， 2008; 黄小波等， 2018）、优势树种的空间分布格局（赵中华等， 2011; 张玉双等， 2016）、两种或两种以上物种的空间相关性（姚良锦等， 2018）以及空间分布格局与环境之间关系（刘旻霞等， 2018; 闫海冰， 2018）等研究较多，对于种群同科（何增丽等， 2017）或同属（陈俊等， 2018）研究较少。
樟科（Lauraceae）木姜子属（Litsea Lam） 的植物分布非常广泛，在中国大约有72种（陈幼竹， 2004）。亚热带常绿落叶阔叶混交林保存最为完好的地区包括木林子国家级自然保护区（黄永涛， 2016），其中木姜子属植物仅有2种：黄丹木姜子和木姜子，其数量占样地内物种总数量的4.8%。木姜子的药用价值很高，黄丹木姜子具有较高的商业价值（陈幼竹， 2004）。木姜子和黄丹木姜子的研究仅仅只限于药用价值等方面的研究，黄丹木姜子还见于叶凋落量及养分特征研究（董舒， 2016），对于更多的生态学方面的研究未见有报道，黄丹木姜子和木姜子虽不是建群种和先锋种，这类小乔木在种群空间分布格局的研究同样具有重要意义。本研究以湖北木林子大样地中同属的木姜子和黄丹木姜子为研究对象，采用g（r）函数对木姜子属两个种进行空间格局及空间关联分析，探讨其成因，探讨两物种是否存在相似性，同属物种共存机制，并为研究亚热带地区物种共存机制提供参考。
1 研究区概况
研究区（109°59′30″—110°17′58″ E、29°55′59″—30°10′47″ N）位于湖北省鹤峰县，面积为20 838 hm2，属于云贵高原的东北延伸区。该地区最低海拔为1 100 m，最高海拔为2 095.6 m，大陆性季风气候，冬季无严寒，夏季无极端炎热，雨热同季，日照充足，山地性气候明显，年均温15.5 ℃，年均降水量1 700～1 900 mm，年均相对湿度82%，无霜期270～279 d。土壤为黄棕壤、棕壤和黄壤。
木林子15 hm2大样地位于保护区核心区内，海拔为1 588～1 780 m，地势相对缓慢，属于亚热带山地常绿落叶阔叶混交林，主要树种有翅柃（Eurya alata）、多脉青冈（Cyclobalanopsis multinervis）、川陕鹅耳枥（C. fargesiana）、小叶青冈（C. myrsinifolia）、石灰花楸（Sorbus folgneri）、四照花（Dendrobenthamia japonica var. chinensis）、锥栗（Castanea henryi）、白栎（Quercus fabri）、山矾（Symplocos sumuntia）和光叶山矾（S. lancifolia）等（姚兰， 2016）。
2 研究方法
2.1 样地建设及野外调查
2013年，按照CTFS大样地建设与检测技术规范，建立15 hm2固定监测样地，分成375块20 m × 20 m小樣地，用防腐处理的不锈钢管将4角作为永久标识，每个小样地再分为4个10 m × 10 m样方，每个样方再分为4个5 m × 5 m的小样方，对DBH≥1.0 cm的木本植物进行挂牌，并用红色油漆标记（树高1.3 m），对DBH≥1 cm的所有木本植物进行每木检尺，识别物种（科、属、种），调查所有木本植物在样地内的坐标值（以大样地单元西南角为原点，测定x轴和y轴的坐标值）。
2.2 种群径级划分
运用径级代替龄级的方法，结合湖北木林子样地内木姜子和黄丹木姜子的径级特征，将胸径DBH划分为4个径级：Ⅰ级（1 cm≤DBH<2.5 cm），Ⅱ级（2.5 cm≤DBH<5 cm），Ⅲ级（5 cm≤DBH<10 cm），Ⅳ级（DBH≥10 cm），分别为幼树、小树、中树和成年树。
2.3 点格局分析方法
本研究采用点格局分析中成对相关函数g（r）进行分析。g（r）函数是基于对所有定位的成对个体之间距离的格局分析，消除了Ripiey’s K（r）函数在大尺度上产生的累积效应（Ripley， 1976; Thorsten et al.， 2007）。分析时，选取距离尺度为较短边长的1/2（150 m），圆环宽度为1 m，通过199次Monte Carlo随机模拟，利用模拟的最大值和最小值绘制上下包迹线，计算99%的置信区间（Wiegand et al.， 2006， 2009）。 2.4 零模型选择
采用点格局分析，需要检测生境异质性在物种分布中是否存在影响。一般来说，如果在大于10 m尺度上呈聚集分布，那就可以解释为受到生境异质性的影响;采用异质泊松模型和先决条件模型目的在于消除大尺度上生境异质性对空间分布格局的影响，在排除生境异质性后，如果小尺度的格局还是聚集，则说明扩散限制存在作用。所以在进行点格局分析时，不同的分析需要选择不同的零假设模型（null model），选取正确合适的零模型十分重要（梁爽等， 2014）。在本研究中单变量g（r）函数采用完全空间随机模型（complete spatial randomness， CSR）和异质泊松模型（heterogeneous Poisson， HP），双变量g（r）函数选择先决条件模型（antecedent condition， AC）。
文中数据的处理与分析使用R 3.4.1软件中的“spatstat”包中完成点格局分析，使用Sigmaplot 13软件绘制空间分布点，Origin 2018软件绘制点格局分析图。
3 结果与分析
3.1 两个物种空间分布格局
整个样地中（表1），木姜子1 701株，黄丹木姜子2 306株。由图1可知，两个物种的空间分布总体差异较大。木姜子主要分布于低海拔、山脊地区;黄丹木姜子主要分布于高海拔和山谷两侧区域，小部分零星分布于样地内。两个物种的整体空间分布格局的研究（图2）表明，基于CSR零模型：木姜子和黄丹木姜子的聚集尺度分别为0～114 m和0～150 m。基于HP零模型：木姜子的聚集尺度为0～23 m，其他尺度上为微弱的随机或均匀分布;黄丹木姜子的聚集尺度为0～36 m、140～150 m，随机分布尺度为48～139 m。
根据两个物种不同径级的空间分布格局研究（图3，图4）表明，基于CSR零模型：木姜子幼树、小树和中树的聚集尺度分别为0～103 m、0～56 m和0～34 m，成年树在0～16 m尺度表现为聚集分布;黄丹木姜子幼树、小树、中树和成年树分别在0～150 m、0～59 m、0～56 m和0～14 m尺度上均表现为聚集分布。基于HP零模型：木姜子幼树、小树和中树的聚集尺度均为0～22 m，而成年树仅在0～3 m尺度上表现为聚集分布;黄丹木姜子各径级的聚集尺度分别为0～35 m、0～35 m、0～12 m和0～6 m。由此可见，随着径级的增长，木姜子与黄丹木姜子均呈现聚集尺度逐渐减小的现象，特别是成年树仅在小尺度上表现为聚集分布，而随机分布格局现象比较明显。
3.2 两个物种不同径级种内关联性
木姜子不同径级种内关联见图5。基于AC零模型，成年树对幼树在所有尺度上均表现为无关联，而其他不同径级之间在较小尺度（r≤20 m）上表现为正相关，其中小树对幼树、中树对幼树、中树对小树、成年树对小树、成年树对中树正关联尺度分别为0～20 m、0～18 m、0～20 m、0～14 m和0～17 m。
黄丹木姜子不同径级种内关联见图6。基于AC零模型，黄丹木姜子成年树对中树在所有尺度上表现为无相关，小树对幼树、中树对幼树、成年树对幼树、中树对小树、成年树对小树分别在0～36 m、0～25 m、8～15 m、0～23 m和0～6 m尺度上呈现正关联。
3.3 两个物种种间关联性
用先决条件模型排除30 m以上生境异质性，结果显示：整体上，木姜子对黄丹木姜子在0～30 m呈现正关联，黄丹木姜子对木姜子在7～29 m呈现正关联（图7）。不同径级之间，木姜子不同径级对黄丹木姜子不同径级空间关联结果显示（图8）：木姜子幼树对黄丹木姜子幼树在4～34 m尺度表现为正关联;木姜子小树对黄丹木姜子幼树在13～29 m尺度表现为正关联，42～67 m尺度上为负关联;木姜子中树对黄丹木姜子幼树在40～50 m尺度上为负关联;木姜子小树对黄丹木姜子小树在0～25 m尺度上为正关联，44～60 m尺度上为负关联;木姜子中树对黄丹木姜子中树在0～5 m尺度上为正关联，在47～54 m和119～128 m尺度上表现为负关联;其余径级之间均呈现无关联。黄丹木姜子不同径级对木姜子不同径级空间结果显示（图9）：黄丹木姜子幼树对木姜子幼树在13～32 m尺度呈现正关联，42～76 m尺度呈现负关联;黄丹木姜子小树对木姜子幼树在0～13 m尺度呈现正关联;黄丹木姜子小树对木姜子小树在0～15 m呈现正关联，46～59 m尺度呈现负关联;其余径级之间均呈现无关联。
由此可见，木姜子对黄丹木姜子空间关系和黄丹木姜子对木姜子空间关系大致相同，仅有部分细微差别，其中两物种成年树与另一物种所有径级均没有显著的生态学关系。
4 讨论
4.1 两个物种空间分布格局形成的原因
植物自身生物学特性和环境因素共同作用于种群的空间分布格局，同一树种在不同的发育阶段将呈現不同的空间分布模式（Frost & Rydin， 2000）。本研究中基于CSR和HP零模型两个物种在研究尺度上主要表现为聚集分布，随着尺度增加表现为随机分布或均匀分布，结果与大部分研究一致（祝燕等， 2011; Guo et al.， 2013）。在很大程度上，生境的异质性被认为是植被空间格局的主要原因。通常，如果它的聚集尺度大于10 m，则可以解释为受到栖息地异质性的影响（Takyu et al.， 2002; Getzin et al.， 2008）。基于CSR零模型，在整体上木姜子和黄丹木姜子分别在0～114 m和0～150 m尺度范围内呈聚集分布，并使用HP零模型排除了30 m以上的生境异质性，聚集尺度减小成为随机分布，可以推测，两个种群空间分布格局受到生境异质性影响。从幼树、小树、中树、成年树时间序列来看，分布格局从聚集分布向随机分布发展趋势，符合一般规律（Vargas et al.， 2010; 张玉双等， 2011）。其中木姜子和黄丹木姜子不同径级之间均体现为幼树聚集尺度>小树聚集尺度>中树聚集尺度>成年树聚集尺度，可见生境异质性对两个物种的幼树影响最强，可能是在此阶段，种群的个体较小，生存力较低，需要通过聚集分布形成相互庇护，能够更好地维持种群自身的稳定性（张金屯， 1998; Wang et al.， 2014）。生境异质性对两个物种的成年树影响较小，两个物种负密度制约超过了生境异质性对聚集分布的影响，随着径级的增大，对资源、对空间的需求也增加，物种内部的竞争日益激烈，物种内部的稀疏性越来越强，这将导致一些个体死亡从而形成均匀或随机分布的状态（舒兰等， 2019）。 YAO LJ， YAO L， YI YM， et al.， 2018. Spatial patterns of dominant species Carpinus fargesiana and Cyclobalanopsis multinervis in a subtropical evergreen and deciduous broad-leaved mixed forest [J]. Sci Silv Sin， 54（12）： 1-11. [姚良锦， 姚兰， 易咏梅， 等， 2018. 亚热带常绿落叶阔叶混交林优势种川陕鹅耳枥和多脉青冈的空间格局 [J]. 林业科学， 54（12）： 1-11.]
ZHANG JT， 1995.Quantitative plant ecology [M]. Beijing： Science and Technology of China Press. [张金屯， 1995. 植被数量生态学方法 [M]. 北京： 中国科学技术出版社.]
ZHANG JT， 1998. Analysis of spatial point pattern for plant species [J]. Chin J Plant Ecol， 22（4）： 344-349. [张金屯， 1998. 植物种群空间分布的点格局分析 [J]. 植物生态学报， 22（4）： 344-349.]
ZHANG JT， MENG DP， 2004. Spatial pattern analysis of individuals in different age-classes of Larix principis-rupprechtii in Luya mountain reserve， Shanxi， China [J]. Acta Ecol Sin，（1）： 35-40. [张金屯， 孟东平， 2004. 芦芽山华北落叶松林不同龄级立木的点格局分析 [J]. 生态学报， （1）： 35-40.]
ZHANG YS， DUAN WB， DU S， et al.， 2016. Analysis on spatial distribution pattern for main populations and gap makers in Korean pine broad-leaved forest in Xiaoxing’anling Mountains of Northeast China [J]. For Res， 29（3）： 395-401. [張玉双， 段文标， 杜珊， 等， 2016. 阔叶红松林主要种群及林隙形成木的空间格局分析 [J]. 林业科学研究， 29（3）： 395-401.]
ZHAO ZH， HUI GY， HU YB， et al.， 2011. Spatial distribution patterns and associations of dominate population in two types of broadleaved Korean pine mixed forest [J]. For Res， 24（5）： 554-562. [赵中华， 惠刚盈， 胡艳波， 等， 2011. 2种类型阔叶红松林优势种群空间分布格局及其关联性 [J]. 林业科学研究， 24（5）： 554-562.]
ZHU Y， BAI F， LIU HF， et al.， 2001. Population distribution patterns and interspecific spatial associations in warm temperate secondary forests， Beijing [M]. Biodivers Sci， 19（2）： 252-259. [祝燕， 白帆， 刘海丰， 等， 2011. 北京暖温带次生林种群分布格局与种间空间关联性 [J]. 生物多样性， 19（2）： 252-259.]
ZONG G， BAI XJ， ZHANG SY， et al.， 2018. Spatial pattern and interspecific spatial association of tree seedlings in a secondary forest in montane region of eastern Liaoning Province， China [J]. Chin J Appl Ecol， 29（1）： 18-24. [宗国， 白雪娇， 张淑媛， 等， 2018. 辽东山区次生林乔木幼苗分布格局与种间空间关联性 [J]. 应用生态学报， 29（1）： 18-24.]
（责任编辑 周翠鸣）