森林是陆地生态系统的主体，具有涵养水源、保育土壤、固碳制氧、保护生物多样性等多种功能^[1]，森林生态系统的稳定性程度直接关系到其生态服务功能的发挥。森林生态系统稳定性评价是森林生态系统研究的重要领域。稳定性是植物群落结构与功能的综合特征^[2]。自群落稳定性概念首次被提出后，至今已有70多种描述生态系统不同性质的稳定性概念，其内涵可归结为系统对外界干扰的抵抗力(resistance)和受扰动后返回初始状态的恢复力(resilience)等2个方面^[3-4]。稳定性维持机制同样引起了国内外学者的深入探讨，已有理论中较完善的有多样性或复杂性理论、食物网理论、冗余理论等^[5]。由于森林群落结构和功能的复杂性以及稳定性与研究区域的气候、种源和干扰等外界因素有关^[2]，所以至今尚没有统一的稳定性测度方法。学者们根据研究对象所选取的测度方法和指标体系不尽相同。如郑元润^[6]在对大青沟森林植物群落的研究中，提出了群落演替与比较分析相结合的方法，并对Godron稳定性测定方法作出改进，近年的稳定性研究中也广泛运用到改进后的Godron法^[7-8]。除此之外，许多研究利用野外调查和实验数据，包括群落组成、更新潜力、物种多样性、土壤理化性质、人为干扰等，建立稳定性评价指标体系；运用主成分分析法^[9]、模糊数学中隶属函数值法^[10-11]、综合指数法^[12]等对森林群落稳定性进行评价。生态系统结构决定其功能。有研究证明：林分空间结构指标也可以较好地反映森林群落的稳定性^[13]。综上所述，稳定性是一个比较复杂的问题^[2]，在研究方法上应综合考虑森林生态系统的群落特征、空间结构、生物量以及外部环境等因素，构建全面且客观的稳定性评价指标体系，并选用合适的方法，以对森林群落稳定性作出定量评价^[13]。辽宁仙人洞国家级自然保护区是以森林生态系统为保护对象的综合性自然保护区，地处华北植物区系与长白植物区系的交错地带，植物地理成分复杂。区内保存有目前亚洲面积最大的原生型赤松Pinus densiflora-栎Quercus林顶级群落，同时还有银杏Ginkgo biloba、东北红豆杉Taxus cuspidata、杜仲Eucomnia ulmoides、野大豆Glycine soja和三桠乌药Lindera obtusiloba等国家级保护植物以及东北罕见的亚热带植物^[14]。该地区森林生态系统在涵养水源、保持水土，调节小峪河、英那河水量水质，促进生态旅游开发和周边第三产业发展等方面具有重大作用。以往的研究多关注保护区内植物区系^[15]、林下土壤^[16]等，对其森林生态系统稳定性尚未有相关研究。因此，本研究以辽宁仙人洞国家级自然保护区内原生的典型森林群落为研究对象，集合种群更新潜力、群落基本特征、物种多样性、种群生态位重叠和凋落物特性等群落自身特征，以及森林土壤肥力等指标，运用模糊综合评判的隶属函数值法对各森林群落的稳定性进行综合评价，以期为该保护区森林生态系统的有效保护及管理提供科学依据。

1.   研究区概况

辽宁仙人洞国家级自然保护区位于辽宁省大连市庄河境内的北部山区(39°54′00″~40°03′00″N，122°53′24″~123°03′30″E)，总面积3 574.7 hm²。地貌为侵蚀、剥蚀低山和丘陵。地带性土壤为石灰岩棕壤性土棕壤(60%)和石灰岩典型棕壤(30%)。年平均气温8.9 ℃，年平均降水量799.0 mm，属于暖温带湿润季风气候，四季温和，雨热同季，并具有一定的海洋性气候特点^[14]。研究区内植物种类丰富，森林植被主要以赤松为主，此外还广泛分布有蒙古栎Quercus mongolica、栓皮栎Q. variabilis、麻栎Q. acutissima、槲树Q. dentata等；灌木层主要有青花椒Zanthoxylum schinifolium、榛Corylus heterophylla、牛叠肚Rubus crataegifolius、白檀Symplocos paniculata、照山白Rhododendron micranthum、胡枝子Lespedeza bicolor等；草本层主要有矮丛薹草Carex callitrichos var. nana、内折香茶菜Isodon inflexus、变色白前Cynanchum versicolor、辽宁堇菜Viola rossii、京芒草Achnatherum pekinense、菴闾Artemisia keiskeana等。

2.   研究方法

2.1   样地设置及调查

辽宁仙人洞国家级自然保护区的典型地带性植被为赤松林、赤松栎林和多种栎组成的落叶阔叶林，集中分布于核心区内，受林窗干扰较小^[14]。于2018年6—7月，在辽宁仙人洞国家级自然保护区核心区内选取9个典型森林群落(赤松林、赤松栎林和落叶阔叶林3种森林植被类型的代表性群落各3个)为研究对象。在每个森林群落内设置20 m×30 m的调查样地1块(表1)，其经纬度、海拔高度、坡度、坡向等信息利用易测宝地理信息系统(GIS)和手持罗盘测定。对样地内所有胸径≥3 cm的树木进行每木检尺，记录其树种名称、树高、胸径、冠幅、树龄等；在每块样地的中心和四角共设置5 m×5 m的灌木样方5个，调查记录其种类组成、基径、丛高、冠幅和总覆盖度等；在每个灌木样方中心各设置1个1 m×1 m的草本样方，调查记录草本层植物的总覆盖度、种类组成及各植物种的盖度、株丛数和平均高度。在每块20 m×30 m样地内沿对角线等分选取3个土样、凋落物采集点，采集0~20 cm层扰动土样和表层原状土样(容积100 cm³的环刀)，其中扰动土样用于测定土壤比重、有机质质量分数和主要养分质量分数，原状土样用于测定土壤容重；同时采用1 m×1 m样方收获法区分未分解层和半分解层，采集凋落物样品，记录其鲜质量和风干质量，以风干质量测算其蓄积量，并采用浸泡实验法测定了其最大持水率和最大持水量。

表 1  森林群落调查样地基本情况

Table 1  Basic situation of the sample plots in the forest communities

样地 编号	群落名称	群落类型	北纬 (N)	东经 (E)	海拔/ m	坡度/ (°)	坡向/ (°)	郁闭度
1	赤松-青花椒-矮丛薹草群落	针叶林	39°58′55.5″	122°57′48.5″	176.6	14	198	0.70
2	蒙古栎+赤松-牛叠肚+照山白-矮丛薹草群落	针阔混交林	39°59′24.4″	122°56′19.1″	231.8	26	230	0.85
3	蒙古栎-牛叠肚-矮丛薹草群落	落叶阔叶林	39°59′22.5″	122°56′18.1″	229.8	26	278	0.89
4	赤松-牛叠肚+青花椒-矮丛薹草群落	针叶林	39°59′19.4″	122°56′15.6″	243.4	21	308	0.68
5	赤松-三桠乌药+照山白-矮丛薹草群落	针叶林	39°59′42.8″	122°56′05.3″	239.3	19	259	0.65
6	赤松+麻栎-地锦Parthenocissus tricuspidata +青花椒- 两型豆Amphicarpaea edgeworthii +白花碎米荠 Cardamine leucantha群落	针阔混交林	39°59′09.1″	122°57′33.7″	161.9	30	113	0.75
7	槲树+麻栎-花木蓝Indigofera kirilowii +榛-菴闾+ 矮丛薹草群落	落叶阔叶林	39°59′01.4″	122°57′46.0″	175.0	25	218	0.75
8	麻栎+栓皮栎-青花椒+榛-矮丛薹草群落	落叶阔叶林	39°59′02.5″	122°57′58.2″	211.5	27	219	0.65
9	槲树+赤松-榛-矮丛薹草	针阔混交林	39°58′58.6″	122°57′58.6″	201.8	31	220	0.70

下载: 导出CSV 

| 显示表格

2.2   评价指标体系

2.2.1   种群更新潜力

选取乔木树种的幼苗、幼树占乔木总植株数的比例作为评价因子，规定胸径D_BH≤1 cm的个体为幼苗，1＜D_BH＜3 cm为幼树，D_BH≥3 cm为成树^[10]。

2.2.2   群落基本特征

选取郁闭度、平均胸径、平均树高、平均年龄和物种种数作为评价因子，其中林下物种种数包括林下灌木层和草本层的所有植物种类。

2.2.3   物种多样性

选取多样性指数(Shannon-Wiener指数和Simpson指数)、均匀度指数(Pielou指数和Alatalo指数)和丰富度指数(Margalef指数)作为评价因子。具体的计算公式如下^[17-18]：Shannon-Wiener指数 $H = - \displaystyle \sum\limits_{i = 1}^s {{P_{{i}}}} \ln {P_i}$。Simpson指数$D = 1 - \displaystyle \sum\limits_{i = 1}^s {P_i^2} $。Pielou指数 $J = H/\ln S$。Alatalo指数 $E_{\rm{a}} = $$ \left[ {{{\left( {\displaystyle \sum\limits_{i = 1}^s {P_i^2} } \right)}^{ - 1}} - 1} \right]/\left[ {\exp {H} - 1} \right]$。Margalef指数$M_{\rm{a}} = (S - 1)/\ln N$。其中：S表示物种数；P_i表示第i个物种的相对重要值，即物种i占样地内所有物种重要值之和的比例；N表示样地内全部个体总数。

2.2.4   种群生态位重叠

选取Schoener指数计算的灌木层和草本层种群生态位重叠作为评价因子，其计算公式如下^[19]：Schoener指数${O_{ik}} = 1 - \dfrac{1}{2}\displaystyle \sum\limits_{j = 1}^r {\left| {{P_{ij}} - {P_{kj}}} \right|} $。其中：O_ik为种i与种k的生态位重叠，0≤O_ik≤1，P_ij和P_kj分别为种i与种k在第j个资源水平上的重要值占该种在所有资源水平上重要值之和的比例；r代表资源水平数，本研究采用样地数。

2.2.5   凋落物特性

选取未分解层和半分解层凋落物的蓄积量、最大持水量、最大持水率，以及凋落物总蓄积量和总最大持水量作为评价因子。以风干质量估算不同林分凋落物的蓄积量，采用室内浸泡法测定凋落物的最大持水量及其吸水速率。

2.2.6   土壤肥力

选取土壤容重、土壤比重和有机质、有效钾、碱解氮、有效磷质量分数作为评价因子。土壤容重和比重分别采用环刀法和比重计法测定，土壤有机质质量分数采用重铬酸钾氧化-外加热法测定，土壤有效钾、有效氮和有效磷质量分数采用X荧光光谱法测定。

2.3   研究方法

采用主成分分析法比较6项评价指标对研究区9个典型森林群落稳定性的影响程度^[20]。基于以上6项指标应用模糊数学中隶属函数值法进行稳定性综合评价。参与评价因子原始数据的标准化处理方法参照文献[10]。隶属函数值的计算公式如下：$(U){X_{ij}} = \left( {{X_{ij}} - {X_{j\min }}} \right)/\left( {{X_{j\max }} - {X_{j\min }}} \right)$。其中：(U)X_ij为第i个群落第j项评价指标的隶属函数值，且(U)X_ij∈[0,1]；X_ij为第i个群落第j项评价指标所有因子标准化数据的平均值；X_jmin和X_jmax分别为所有参与评价群落第j项评价指标所有因子标准化数据平均值的最小值和最大值^[10]。通过计算每个森林群落各项指标隶属函数值的平均值，得出各森林群落稳定性的高低排名。采用单因素方差分析(one-way ANOVA)和最小显著差法(LSD)对不同森林群落中物种多样性、凋落物特性和土壤肥力指标进行显著性检验。显著性水平设定α=0.05。

3.   结果与分析

3.1   种群更新潜力

由表2可知：9个森林群落的幼苗、幼树比例存在明显差异。赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落幼树比例最高、幼苗比例较高，树种丰富，以花曲柳Fraxinus rhynchophylla、栓皮栎和槲栎为主，种群更新潜力最好。麻栎+栓皮栎-青花椒+榛-矮丛薹草群落幼苗比例最高(55.97%)，但幼树很少，种群更新潜力次之。槲树+赤松-榛-矮丛薹草群落和赤松-青花椒-矮丛薹草群落幼苗、幼树比例中等，但组成存在差异；槲树+赤松-榛-矮丛薹草群落的幼树以槲树为主，赤松-青花椒-矮丛薹草群落的幼树种类相对丰富，如赤松、花曲柳和核桃楸Juglans mandshurica等。两者种群更新潜力相近，其他群落仅有少量幼苗或幼树，且种类简单，种群更新潜力较低。

表 2  森林群落种群更新潜力

Table 2  Regenerative potential of the forest population

样地 编号	幼苗比 例/%	幼树比 例/%	标准化数据 平均值		样地 编号	幼苗比 例/%	幼树比 例/%	标准化数据 平均值		样地 编号	幼苗比 例/%	幼树比 例/%	标准化数据 平均值
1	12.62	17.54	460.26		4	14.49	4.08	210.34		7	17.07	0.00	152.52
2	9.72	1.34	113.37		5	16.13	1.06	165.13		8	55.97	7.83	655.09
3	7.29	2.08	106.29		6	26.14	25.24	733.48		9	16.79	16.69	480.77

下载: 导出CSV 

| 显示表格

3.2   群落基本特征

由表3可知：蒙古栎+赤松-牛叠肚+照山白-矮丛薹草群落、蒙古栎-牛叠肚-矮丛薹草群落和赤松-牛叠肚+青花椒-矮丛薹草群落平均树高均大于10 m，说明群落内乔木的生存空间处于较上层，但由于蒙古栎+赤松-牛叠肚+照山白-矮丛薹草群落和蒙古栎-牛叠肚-矮丛薹草群落均以蒙古栎为优势林层且株数多，蒙古栎分枝多叶片大，树冠覆盖范围广，故两者林冠郁闭度均大于0.80，而赤松-牛叠肚+青花椒-矮丛薹草群落的优势树种赤松株数较少，故林冠郁闭度仅0.68；赤松-青花椒-矮丛薹草群落和赤松-三桠乌药+照山白-矮丛薹草群落郁闭度、平均树高等与其他群落差异不大，但草本层植物种类较丰富，如尾叶香茶菜Isodon excisus、甘菊Chrysanthemum lavandulifolium和长萼鸡眼草Kummerowia stipulacea等仅发现于赤松-青花椒-矮丛薹草林群落，毛缘宽叶薹草Carex siderosticta var. pilosa、中华蹄盖蕨Thyrium sinense和紫草Lithospermum erythrorhizon等仅发现于赤松-三桠乌药+照山白-矮丛薹草群落，说明林下植被的种类除受林冠层影响，还可能与土壤、坡度等环境因素有关。

表 3  森林群落基本特征

Table 3  General characteristics of the forest communities

样地 编号	郁闭度	平均树 高/m	平均胸 径/cm	平均树 龄/a	物种种 数/种	标准化数据 平均值
1	0.70	9.27	14.22	25.71	40	794.38
2	0.85	10.33	12.71	24.64	31	783.55
3	0.89	11.63	13.16	22.26	35	824.52
4	0.68	11.26	18.03	36.19	35	902.02
5	0.65	8.64	13.35	37.20	45	848.27
6	0.75	9.55	17.36	33.07	36	868.13
7	0.75	9.57	13.93	33.13	28	795.22
8	0.65	8.80	12.96	25.16	31	717.87
9	0.70	7.41	10.75	24.82	25	652.28

下载: 导出CSV 

| 显示表格

3.3   物种多样性

森林群落具有明显的垂直分层结构。由表4可以看出：研究区各森林群落不同功能层的物种多样性表现不一致。各群落乔木层物种多样性差异明显，槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落、赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落和槲树+赤松-榛-矮丛薹草群落的Shannon-Wiener指数和Simpson指数较大，最小的为蒙古栎-牛叠肚-矮丛薹草群落；Pielou指数和Alatalo指数的平均值最大的为槲树+赤松-榛-矮丛薹草群落，蒙古栎-牛叠肚-矮丛薹草群落和赤松-牛叠肚+青花椒-矮丛薹草群落的均匀度指数较小；赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落的Margalef指数最大，蒙古栎+赤松-牛叠肚+照山白-矮丛薹草群落和蒙古栎-牛叠肚-矮丛薹草群落乔木层仅有蒙古栎和赤松2个种群，Margalef指数较小。不同森林群落间灌木层Shannon-Wiener指数差异极显著(P＜0.01)，Margalef指数差异不显著(P＞0.05)；草本层仅Pielou指数差异显著(P＜0.05)。将森林群落各层次综合来看，9个森林群落物种多样性指标差异从大到小依次表现为赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落、赤松-三桠乌药+照山白-矮丛薹草群落、槲树+赤松-榛-矮丛薹草群落、槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落、麻栎+栓皮栎-青花椒+榛-矮丛薹草群落、赤松-青花椒-矮丛薹草林群落、蒙古栎+赤松-牛叠肚+照山白-矮丛薹草群落、赤松-牛叠肚+青花椒-矮丛薹草群落、蒙古栎-牛叠肚-矮丛薹草群落。这与各功能层物种多样性变化在一定程度上具有一致性。

表 4  森林群落物种多样性指数

Table 4  Species diversity indexes of the forest communities

林层	指数	不同样地的多样性指数									F	P
		1	2	3	4	5	6	7	8	9
乔木层	H	0.91	0.62	0.11	0.34	0.78	1.28	1.39	0.93	1.09
	D	0.48	0.43	0.05	0.16	0.40	0.65	0.72	0.57	0.66
	J	0.66	0.89	0.16	0.25	0.44	0.66	0.86	0.85	0.99
	Ea	0.62	0.87	0.41	0.47	0.57	0.71	0.86	0.88	0.98
	Ma	0.70	0.24	0.22	0.74	1.12	1.61	0.95	0.51	0.45
灌木层	H	1.54±0.19 ab	1.51±0.12 b	1.49±0.42 b	1.38±0.25 b	1.94±0.10 a	1.91±0.31 ab	0.83±0.47 c	1.55±0.53 ab	1.28±0.34 b	4.915	＜0.001
	D	0.75±0.04 ab	0.74±0.03 ab	0.66±0.12 b	0.65±0.10 b	0.83±0.02 a	0.79±0.08 ab	0.46±0.26 c	0.74±0.14 ab	0.64±0.12 b	4.002	0.002
	J	0.90±0.08 ab	0.86±0.07 ab	0.74±0.11 b	0.77±0.10 b	0.90±0.03 ab	0.82±0.07 ab	0.84±0.11 ab	0.92±0.04 a	0.80±0.08 b	2.847	0.015
	Ea	0.83±0.11 ab	0.81±0.10 ab	0.61±0.08 b	0.66±0.12 b	0.82±0.04 ab	0.71±0.14 b	0.80±0.16 ab	0.87±0.08 a	0.74±0.11 ab	3.068	0.010
	Ma	1.71±0.70 ab	1.49±0.30 ab	1.74±0.80 ab	1.56±0.58 ab	2.26±0.25 a	2.27±0.80 a	0.92±0.66 b	1.71±0.92 ab	1.28±0.49 b	2.203	0.051
草本层	H	1.71±0.52 ab	1.50±0.16 b	1.55±0.25 ab	1.45±0.26 b	1.95±0.35 a	1.34±0.28 b	1.65±0.38 ab	1.44±0.38 b	1.53±0.26 ab	1.494	0.194
	D	0.72±0.18 a	0.71±0.06 a	0.72±0.09 a	0.73±0.09 a	0.82±0.07 a	0.71±0.08 a	0.75±0.10 a	0.65±0.13 b	0.74±0.07 a	0.953	0.487
	J	0.79±0.15 b	0.84±0.06 ab	0.82±0.08 b	0.81±0.06 b	0.89±0.05 ab	0.94±0.07 a	0.82±0.05 b	0.72±0.10 b	0.90±0.07 ab	3.140	0.008
	Ea	0.66±0.15 b	0.74±0.09 ab	0.73±0.11 ab	0.93±0.36 a	0.81±0.09 ab	0.90±0.12 a	0.74±0.04 ab	0.63±0.11 b	0.83±0.11 ab	2.040	0.069
	Ma	2.16±0.72 ab	1.71±0.26 ab	1.72±0.36 ab	1.38±0.44 b	2.32±0.73 a	1.61±0.43 b	1.79±0.69 ab	1.66±0.63 ab	1.65±0.50 ab	1.374	0.241
标准化数 据平均值		767.53	738.20	576.61	627.84	827.03	862.28	774.82	769.48	815.92
说明：H. Shannon-Wiener指数；D. Simpson指数；J. Pielou指数；Ea. Alatalo指数；Ma. Margalef指数。同一行中不同字母表示差异显著(P＜0.05)，　　　　　相同字母表示差异不显著

下载: 导出CSV 

| 显示表格

3.4   种群生态位重叠

由表5可知：多数群落灌木层种群生态位重叠值为0的种对占全部种对的比例低于草本层，说明灌木层种对重叠情况多于草本层。草本层存在一些生态位完全重叠的种对，说明草本层部分植物对生态因子的需求相似性较高。各群落种群生态位重叠值主要集中在0~0.6，说明群落内种间生态位分化程度低，种群生态位重叠现象普遍存在但程度较弱，不同物种共同利用或占有某一资源的竞争较缓和。从总体看，9个森林群落种群生态位重叠分布格局相似，且重叠程度不明显。

表 5  森林群落种群生态位重叠指数分布格局

Table 5  Distribution pattern of niche overlap index of the forest communities

样地 编号	灌木层						草本层					标准化数 据平均值
	Oik=0	0＜Oik≤0.3	0.3＜Oik≤0.6	0.6＜Oik≤0.9	Oik＞0.9		Oik=0	0＜Oik≤0.3	0.3＜Oik≤0.6	0.6＜Oik≤0.9	Oik＞0.9
1	35.90	23.08	29.49	6.41	5.13		55.27	10.54	15.10	10.26	8.83	620.16
2	37.18	20.51	32.05	6.41	3.85		58.82	16.34	15.69	3.27	5.88	556.04
3	20.00	30.48	30.48	16.19	2.86		46.32	18.42	24.21	2.63	8.42	617.82
4	35.16	31.87	21.98	6.59	4.40		61.43	12.86	11.90	3.33	10.48	568.91
5	28.68	26.47	30.88	13.24	0.74		56.08	15.87	15.34	4.50	8.20	573.29
6	48.94	18.78	19.05	7.67	5.56		68.63	5.88	6.54	3.92	15.03	581.55
7	47.22	13.89	27.78	2.78	8.33		35.09	24.56	21.64	11.11	7.60	684.30
8	19.70	22.73	42.42	15.15	0.00		45.03	16.37	26.32	7.02	5.26	602.97
9	21.82	45.45	20.00	7.27	5.45		42.86	20.88	29.67	2.20	4.40	598.57
说明：Oik=0表示群落内种群生态位重叠值为0的种对占全部种对的比例；0＜Oik≤0.3表示群落内种群生态位重叠值在此范围的种对　　　　　占全部种对的比例

下载: 导出CSV 

| 显示表格

3.5   凋落物特性

由表6可知：9个森林群落的凋落物蓄积总量为2.87~16.33 t·hm⁻²，差异极显著(P＜0.01)。蒙古栎-牛叠肚-矮丛薹草群落的半分解层凋落物蓄积量最高，槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落未分解层凋落物蓄积量和总蓄积量最高，赤松-三桠乌药+照山白-矮丛薹草群落和赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落的凋落物总蓄积量也较高。凋落物蓄积量是影响凋落物层持水性能的特征之一^[21]，最大持水量决定最大持水率。9个森林群落凋落物总最大持水量为70.90~589.04 t·hm⁻²，差异极显著(P＜0.01)，槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落凋落物各层最大持水量均为最高，赤松-青花椒-矮丛薹草群落各层最大持水率和最大持水量均较低，而各群落内部未分解层和半分解层的最大持水量和持水率相差不大。综合比较，槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落凋落物层表现最好，在拦蓄降水、改良土壤肥力方面作用显著。

表 6  森林群落凋落物特性

Table 6  Litter characteristics of the forest communities

样地 编号	凋落物蓄积量/(t·hm−2)				最大持水率/%			最大持水量/(t·hm−2)			标准化数 据平均值
	未分解层	半分解层	合计		未分解层	半分解层		未分解层	半分解层	合计
1	2.26±0.79 d	2.08±0.26 c	4.34±1.05 d		175.54±25.73 b	176.24±36.46 b		3.97±1.49 d	3.67±0.56 c	7.63±1.86 d	285.51
2	6.14±1.21 b	9.64±2.66 a	15.78±2.10 a		348.50±154.58 ab	282.72±70.07 ab		21.40±7.91 b	27.25±10.05 a	48.65±13.64 ab	854.47
3	4.71±0.22 c	5.31±0.42 b	10.03±0.57 bc		343.30±18.57 ab	306.29±49.50 ab		16.17±1.12 bc	16.26±3.39 b	32.43±4.02 bc	666.81
4	1.86±0.35 d	1.01±0.38 c	2.87±0.65 d		191.29±21.48 b	349.69±21.47 a		3.56±0.79 d	3.53±1.14 c	7.09±1.50 d	344.49
5	5.60±1.34 bc	8.66±1.74 ab	14.26±2.19 ab		251.49±42.77 b	290.09±8.51 ab		14.08±3.42 c	25.12±4.71 a	39.20±3.04 b	734.58
6	5.14±0.12 bc	7.13±1.26 ab	12.27±1.21 b		348.29±41.26 ab	301.46±15.81 ab		17.90±1.64 bc	21.50±3.28 ab	39.40±4.55 b	746.50
7	8.02±0.27 a	8.31±3.66 ab	16.33±3.90 a		369.74±46.76 a	352.00±210.80 a		29.65±3.05 a	29.25±5.17 a	58.90±4.29 a	970.28
8	4.56±0.57 c	3.24±0.70 bc	7.80±0.45 c		385.47±38.37 a	293.46±49.04 ab		17.58±1.70 bc	9.51±0.73 bc	27.09±2.35 c	609.33
9	4.69±0.59 c	5.75±1.60 b	10.45±2.08 bc		284.81±17.51 ab	197.72±17.66 b		13.36±1.16 c	11.37±4.22 bc	24.73±5.36 c	553.62
F	19.649	8.804	19.148		4.994	1.724		18.487	11.971	24.750
P	＜0.01	＜0.01	＜0.01		＜0.01	0.161		＜0.01	＜0.01	＜0.01
说明：不同小写字母表示不同样地间差异显著(P＜0.05)

下载: 导出CSV 

| 显示表格

3.6   土壤肥力

由表7可知：9个森林群落的土壤比重差异不显著(P＞0.05)，容重差异极显著(P＜0.01)，反映出各群落土壤固相组成物质的种类及含量相似但紧实程度不同。槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落容重较大，孔隙较少，其储水性能较好，透水性差。不同森林群落土壤有机质、有效钾、碱解氮和有效磷质量分数差异极显著(P＜0.01)，赤松-三桠乌药+照山白-矮丛薹草群落和赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落的有机质质量分数较高且显著高于其他群落；槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落的有效钾质量分数最高，达2.65%；赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落和赤松-三桠乌药+照山白-矮丛薹草群落的碱解氮和有效磷质量分数较高，但与其他群落相差不大。综合比较，赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落的土壤肥力指标排名最高。

表 7  森林群落土壤肥力

Table 7  Soil fertilities of the forest communities

样地编号	土壤容重/(g·cm−3)	土壤比重/(g·cm−3)	有机质/%	有效钾/%	碱解氮/%	有效磷/%	标准化数据平均值
1	0.93±0.03 c	2.67±0.12 a	3.74±0.72 c	2.40±0.09 ab	0.50±0.06 bc	0.04±0.01 c	707.28
2	1.11±0.07 a	2.65±0.05 a	2.97±1.04 c	2.37±0.22 ab	0.46±0.06 c	0.03±0.01 c	683.19
3	0.98±0.07 bc	2.75±0.08 a	4.21±1.05 bc	2.00±0.14 bc	0.56±0.03 bc	0.03±0.00 c	706.72
4	0.92±0.05 c	2.72±0.06 a	5.14±0.44 bc	2.08±0.05 b	0.61±0.03 ab	0.04±0.01 bc	754.71
5	0.68±0.06 d	2.63±0.07 b	7.81±1.07 a	1.98±0.01 bc	0.71±0.05 a	0.05±0.00 b	805.70
6	0.91±0.10 c	2.68±0.02 a	7.38±3.31 a	1.84±0.24 bc	0.69±0.16 ab	0.09±0.02 a	888.54
7	1.22±0.15 ab	2.65±0.08 a	3.51±0.40 c	2.65±0.08 a	0.47±0.03 c	0.04±0.00 c	746.22
8	1.10±0.08 ab	2.67±0.16 a	5.32±0.42 b	2.44±0.07 a	0.60±0.02 b	0.05±0.00 b	812.33
9	1.06±0.01 b	2.74±0.06 a	2.85±0.05 c	1.66±0.49 c	0.34±0.03 d	0.02±0.01 c	601.62
F	11.778	1.320	11.626	7.354	9.645	23.214
P	＜0.01	0.265	＜0.01	＜0.01	＜0.01	＜0.01
说明：不同小写字母表示不同样地间差异显著(P＜0.05)

下载: 导出CSV 

| 显示表格

3.7   群落稳定性评价结果

综合比较各森林群落6项指标隶属函数值平均值的大小可以评判森林群落稳定性状况。从评价结果可知(表8)：9个森林群落稳定性从大到小依次为赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落、槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落、麻栎+栓皮栎-青花椒+榛-矮丛薹草群落、赤松-三桠乌药+照山白-矮丛薹草群落、赤松-青花椒-矮丛薹草群落、蒙古栎+赤松-牛叠肚+照山白-矮丛薹草群落、蒙古栎-牛叠肚-矮丛薹草群落、赤松-牛叠肚+青花椒-矮丛薹草群落、槲树+赤松-榛-矮丛薹草群落。

表 8  森林群落各指标隶属函数值及平均值

Table 8  Subordinate function values of community stability

样地编号	种群更新潜力			群落基本特征			物种多样性			种群生态位重叠			凋落物特性			土壤特性			平均值
	数值	排序		数值	排序		数值	排序		数值	排序		数值	排序		数值	排序		数值	排序
1	0.56	4		0.57	6		0.67	6		0.50	2		0.00	9		0.37	6		0.44	5
2	0.01	8		0.53	7		0.57	7		0.00	9		0.83	2		0.28	8		0.37	6
3	0.00	9		0.69	4		0.00	9		0.48	3		0.56	5		0.37	7		0.35	7
4	0.17	5		1.00	1		0.18	8		0.10	8		0.09	8		0.53	4		0.34	8
5	0.09	6		0.78	3		0.88	2		0.13	7		0.66	4		0.71	3		0.54	4
6	1.00	1		0.86	2		1.00	1		0.20	6		0.67	3		1.00	1		0.79	1
7	0.07	7		0.57	5		0.69	4		1.00	1		1.00	1		0.50	5		0.64	2
8	0.88	2		0.26	8		0.68	5		0.37	4		0.47	6		0.73	2		0.56	3
9	0.60	3		0.00	9		0.74	3		0.33	5		0.39	7		0.00	9		0.34	9

下载: 导出CSV 

| 显示表格

3.8   森林群落稳定性影响因素主成分分析

对种群更新潜力等6项森林群落稳定性评价指标的主成分分析结果显示(表9)：前3个主成分的累积贡献率达到81.17%，其中第1主成分特征值为1.90，贡献率为31.66%。由表10可见：第1主成分中种群更新潜力、物种多样性和土壤肥力的荷载系数最大，分别为0.79、0.82和0.75，表明这3项指标对森林群落稳定性的影响程度较大。

表 9  主成分分析结果

Table 9  Results of principal components analysis

主成分	特征值	贡献率/%	累积贡献率/%
1	1.90	31.66	31.66
2	1.61	26.85	58.51
3	1.36	22.66	81.17

下载: 导出CSV 

| 显示表格

表 10  森林群落稳定性影响因素载荷矩阵

Table 10  Component matrix for forest community stability factors

影响因素	主成分荷载
	1	2	3
种群更新潜力	0.79	−0.31	−0.42
群落基本特征	0.05	0.95	0.12
物种多样性	0.82	−0.33	0.16
种群生态位重叠	−0.14	−0.42	0.55
凋落物特性	0.14	−0.07	0.90
土壤肥力	0.75	0.56	0.20

下载: 导出CSV 

| 显示表格

4.   讨论与结论

4.1   讨论

4.1.1   种群更新潜力对群落稳定性的影响

森林群落的稳定性由群落内外多种因素共同决定，群落内的种群差异可以看作稳定性差异的表现。乔木层作为森林群落的优势层，对群落具有支配作用^[22]。不仅林冠郁闭度会影响林下植被的种类与数量，同时乔木幼苗和幼树是林分天然更新的基础，因此群落的种群更新潜力高，森林群落持续存在且健康发展的能力好，在一定程度上提高了群落的稳定性。本研究中，赤松和麻栎群落的种群更新潜力相对最高，更新层中栎属Quercus乔木如栓皮栎和槲栎的幼苗、幼树较多，其根系较发达且耐荫性强，能够较好地适应林下的荫蔽环境，逐渐生长并进入优势林层，因此就种群更新潜力而言，该群落稳定性最高。这与张梦弢等^[23]的研究结果一致，即云杉Picea jezoensis-冷杉Abies nephrolepis针阔混交林更新潜力明显好于其他2种森林群落，稳定性也最好。

4.1.2   物种多样性对群落稳定性的影响

物种多样性是描述群落结构和功能的数量性指标，多样性和稳定性的关系也一直备受关注。研究者们通过生物多样性实验、数学模型和理论综述等方法对其展开研究，但结论并不完全一致^[24]。大量实证研究表明：物种多样性通过功能多样性和功能冗余的增加实现对群落稳定性的正效应^[25]。物种多样性高的森林群落，不仅包含的种群数量多，而且各种群对生态因子的需求存在差异，因此面临干扰时，群落抵抗波动的能力强，受到扰动后的恢复力也强^[26]，有利于群落的稳定。但部分学者结合数学模型得出结论，多样性和稳定性呈负相关关系或者两者之间并不存在明显的联系^[27-28]。本研究中，各群落的稳定性高低与物种多样性大小不一致，如赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落、槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落和麻栎+栓皮栎-青花椒+榛-矮丛薹草群落的物种多样性递减，其稳定性也逐渐降低；而蒙古栎-牛叠肚-矮丛薹草群落物种多样性最低，稳定性却高于赤松-牛叠肚+青花椒-矮丛薹草群落和槲树+赤松-榛-矮丛薹草群落。该研究结果与于晓文等^[2]对冀北辽河源河岸林群落的研究结果相似。这表明物种多样性是群落稳定的必要条件，但物种多样性与群落稳定性并非一定是正相关关系。对这一情况的解释主要有以下2种：①群落各功能层对稳定性的影响程度不同，导致不同层次物种多样性与群落总体稳定性的关系表现并不一致。例如张明霞等^[9]对秦岭松栎混交林的研究结果显示：松栎混交林的稳定性与乔木层的物种多样性呈负相关，与林下植被层的物种多样性呈正相关；②可能存在着1个多样性阈值，在该值以下，物种多样性表现为调节控制和种间协调功能，因此物种多样性与稳定性呈正相关关系；高于该值后，物种多样性变化对群落结构和功能的影响不大，群落稳定性差异不明显^[2]。另外，由于稳定性存在多种定义，同一系统可能表现出因不同稳定性定义而产生不同的多样性—稳定性关系。

4.1.3   土壤肥力对群落稳定性的影响

土壤物理因素和化学因素等共同决定了土壤肥力大小。本研究中，槲树+赤松-榛-矮丛薹草群落土壤容重较大，土壤有机质质量分数仅2.85%，同时群落土壤碱解氮、有效磷、有效钾质量分数低，土壤肥力较差，制约着森林植被长期持续发展，导致群落稳定性较低。土壤作为森林群落必不可少的组成部分，为植被的存在和生长提供水、肥、气、热等必需条件。土壤容重与土壤质地、结构和有机质质量分数等有关，影响着土壤中水分、空气和养分的运输；土壤有机质作为土壤肥力的重要指标之一，是植物营养元素的给源，也是制约土壤理化性质的关键因素；土壤中氮、磷、钾是植物生长最主要的三大营养物质^[29]，它们的含量与有机质含量显著正相关。凋落物层作为森林群落中联结土壤与植物的纽带，是森林土壤有机碳的主要来源，影响土壤微生物群落活性和功能发挥^[30]，不同分解和积累程度凋落物下，土壤的养分组成和腐殖质性质都具有明显的不同。因此，积累较厚且分解程度高的凋落物层能有效减弱土壤淋溶作用，利于氮磷富集和有机质积累，土壤肥力得到改善^[31]，有利于植被生长和群落的稳定；反之，若林地内凋落物层稀薄，则导致土壤肥力较低，不利于森林群落的持续发展。

4.2   结论

本研究综合考虑森林群落自身特征和土壤肥力，应用模糊综合评判的隶属函数值法，对辽宁仙人洞国家级自然保护区9个森林群落的稳定性进行了综合评价，取得以下结论：①辽宁仙人洞国家级自然保护区的9个典型森林群落稳定性差异显著，由高到低依次为赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落、槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落、麻栎+栓皮栎-青花椒+榛-矮丛薹草群落、赤松-三桠乌药+照山白-矮丛薹草群落、赤松-青花椒-矮丛薹草群落、蒙古栎+赤松-牛叠肚+照山白-矮丛薹草群落、蒙古栎-牛叠肚-矮丛薹草群落、赤松-牛叠肚+青花椒-矮丛薹草群落、槲树+赤松-榛-矮丛薹草群落。②种群更新潜力影响群落的恢复力，物种多样性提高了群落的抵抗力，土壤肥力为群落发展提供良好基础。三者是影响森林群落稳定性的主要因素。