-
森林是陆地生态系统的主体,具有涵养水源、保育土壤、固碳制氧、保护生物多样性等多种功能[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]等,对其森林生态系统稳定性尚未有相关研究。因此,本研究以辽宁仙人洞国家级自然保护区内原生的典型森林群落为研究对象,集合种群更新潜力、群落基本特征、物种多样性、种群生态位重叠和凋落物特性等群落自身特征,以及森林土壤肥力等指标,运用模糊综合评判的隶属函数值法对各森林群落的稳定性进行综合评价,以期为该保护区森林生态系统的有效保护及管理提供科学依据。
-
辽宁仙人洞国家级自然保护区位于辽宁省大连市庄河境内的北部山区(39°54′00″~40°03′00″N,122°53′24″~123°03′30″E),总面积3 574.7 hm2。地貌为侵蚀、剥蚀低山和丘陵。地带性土壤为石灰岩棕壤性土棕壤(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等。
-
辽宁仙人洞国家级自然保护区的典型地带性植被为赤松林、赤松栎林和多种栎组成的落叶阔叶林,集中分布于核心区内,受林窗干扰较小[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 cm3的环刀),其中扰动土样用于测定土壤比重、有机质质量分数和主要养分质量分数,原状土样用于测定土壤容重;同时采用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 -
选取乔木树种的幼苗、幼树占乔木总植株数的比例作为评价因子,规定胸径DBH≤1 cm的个体为幼苗,1<DBH<3 cm为幼树,DBH≥3 cm为成树[10]。
-
选取郁闭度、平均胸径、平均树高、平均年龄和物种种数作为评价因子,其中林下物种种数包括林下灌木层和草本层的所有植物种类。
-
选取多样性指数(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表示物种数;Pi表示第i个物种的相对重要值,即物种i占样地内所有物种重要值之和的比例;N表示样地内全部个体总数。 -
选取Schoener指数计算的灌木层和草本层种群生态位重叠作为评价因子,其计算公式如下[19]:Schoener指数
${O_{ik}} = 1 - \dfrac{1}{2}\displaystyle \sum\limits_{j = 1}^r {\left| {{P_{ij}} - {P_{kj}}} \right|} $ 。其中:Oik为种i与种k的生态位重叠,0≤Oik≤1,Pij和Pkj分别为种i与种k在第j个资源水平上的重要值占该种在所有资源水平上重要值之和的比例;r代表资源水平数,本研究采用样地数。 -
选取未分解层和半分解层凋落物的蓄积量、最大持水量、最大持水率,以及凋落物总蓄积量和总最大持水量作为评价因子。以风干质量估算不同林分凋落物的蓄积量,采用室内浸泡法测定凋落物的最大持水量及其吸水速率。
-
选取土壤容重、土壤比重和有机质、有效钾、碱解氮、有效磷质量分数作为评价因子。土壤容重和比重分别采用环刀法和比重计法测定,土壤有机质质量分数采用重铬酸钾氧化-外加热法测定,土壤有效钾、有效氮和有效磷质量分数采用X荧光光谱法测定。
-
采用主成分分析法比较6项评价指标对研究区9个典型森林群落稳定性的影响程度[20]。基于以上6项指标应用模糊数学中隶属函数值法进行稳定性综合评价。参与评价因子原始数据的标准化处理方法参照文献[10]。隶属函数值的计算公式如下:
$(U){X_{ij}} = \left( {{X_{ij}} - {X_{j\min }}} \right)/\left( {{X_{j\max }} - {X_{j\min }}} \right)$ 。其中:(U)Xij为第i个群落第j项评价指标的隶属函数值,且(U)Xij∈[0,1];Xij为第i个群落第j项评价指标所有因子标准化数据的平均值;Xjmin和Xjmax分别为所有参与评价群落第j项评价指标所有因子标准化数据平均值的最小值和最大值[10]。通过计算每个森林群落各项指标隶属函数值的平均值,得出各森林群落稳定性的高低排名。采用单因素方差分析(one-way ANOVA)和最小显著差法(LSD)对不同森林群落中物种多样性、凋落物特性和土壤肥力指标进行显著性检验。显著性水平设定α=0.05。 -
由表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 -
由表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 -
森林群落具有明显的垂直分层结构。由表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), 相同字母表示差异不显著 -
由表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表示群落内种群生态位重叠值在此范围的种对 占全部种对的比例 -
由表6可知:9个森林群落的凋落物蓄积总量为2.87~16.33 t·hm−2,差异极显著(P<0.01)。蒙古栎-牛叠肚-矮丛薹草群落的半分解层凋落物蓄积量最高,槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落未分解层凋落物蓄积量和总蓄积量最高,赤松-三桠乌药+照山白-矮丛薹草群落和赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落的凋落物总蓄积量也较高。凋落物蓄积量是影响凋落物层持水性能的特征之一[21],最大持水量决定最大持水率。9个森林群落凋落物总最大持水量为70.90~589.04 t·hm−2,差异极显著(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) -
由表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) -
综合比较各森林群落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 -
对种群更新潜力等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 表 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 -
森林群落的稳定性由群落内外多种因素共同决定,群落内的种群差异可以看作稳定性差异的表现。乔木层作为森林群落的优势层,对群落具有支配作用[22]。不仅林冠郁闭度会影响林下植被的种类与数量,同时乔木幼苗和幼树是林分天然更新的基础,因此群落的种群更新潜力高,森林群落持续存在且健康发展的能力好,在一定程度上提高了群落的稳定性。本研究中,赤松和麻栎群落的种群更新潜力相对最高,更新层中栎属Quercus乔木如栓皮栎和槲栎的幼苗、幼树较多,其根系较发达且耐荫性强,能够较好地适应林下的荫蔽环境,逐渐生长并进入优势林层,因此就种群更新潜力而言,该群落稳定性最高。这与张梦弢等[23]的研究结果一致,即云杉Picea jezoensis-冷杉Abies nephrolepis针阔混交林更新潜力明显好于其他2种森林群落,稳定性也最好。
-
物种多样性是描述群落结构和功能的数量性指标,多样性和稳定性的关系也一直备受关注。研究者们通过生物多样性实验、数学模型和理论综述等方法对其展开研究,但结论并不完全一致[24]。大量实证研究表明:物种多样性通过功能多样性和功能冗余的增加实现对群落稳定性的正效应[25]。物种多样性高的森林群落,不仅包含的种群数量多,而且各种群对生态因子的需求存在差异,因此面临干扰时,群落抵抗波动的能力强,受到扰动后的恢复力也强[26],有利于群落的稳定。但部分学者结合数学模型得出结论,多样性和稳定性呈负相关关系或者两者之间并不存在明显的联系[27-28]。本研究中,各群落的稳定性高低与物种多样性大小不一致,如赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落、槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落和麻栎+栓皮栎-青花椒+榛-矮丛薹草群落的物种多样性递减,其稳定性也逐渐降低;而蒙古栎-牛叠肚-矮丛薹草群落物种多样性最低,稳定性却高于赤松-牛叠肚+青花椒-矮丛薹草群落和槲树+赤松-榛-矮丛薹草群落。该研究结果与于晓文等[2]对冀北辽河源河岸林群落的研究结果相似。这表明物种多样性是群落稳定的必要条件,但物种多样性与群落稳定性并非一定是正相关关系。对这一情况的解释主要有以下2种:①群落各功能层对稳定性的影响程度不同,导致不同层次物种多样性与群落总体稳定性的关系表现并不一致。例如张明霞等[9]对秦岭松栎混交林的研究结果显示:松栎混交林的稳定性与乔木层的物种多样性呈负相关,与林下植被层的物种多样性呈正相关;②可能存在着1个多样性阈值,在该值以下,物种多样性表现为调节控制和种间协调功能,因此物种多样性与稳定性呈正相关关系;高于该值后,物种多样性变化对群落结构和功能的影响不大,群落稳定性差异不明显[2]。另外,由于稳定性存在多种定义,同一系统可能表现出因不同稳定性定义而产生不同的多样性—稳定性关系。
-
土壤物理因素和化学因素等共同决定了土壤肥力大小。本研究中,槲树+赤松-榛-矮丛薹草群落土壤容重较大,土壤有机质质量分数仅2.85%,同时群落土壤碱解氮、有效磷、有效钾质量分数低,土壤肥力较差,制约着森林植被长期持续发展,导致群落稳定性较低。土壤作为森林群落必不可少的组成部分,为植被的存在和生长提供水、肥、气、热等必需条件。土壤容重与土壤质地、结构和有机质质量分数等有关,影响着土壤中水分、空气和养分的运输;土壤有机质作为土壤肥力的重要指标之一,是植物营养元素的给源,也是制约土壤理化性质的关键因素;土壤中氮、磷、钾是植物生长最主要的三大营养物质[29],它们的含量与有机质含量显著正相关。凋落物层作为森林群落中联结土壤与植物的纽带,是森林土壤有机碳的主要来源,影响土壤微生物群落活性和功能发挥[30],不同分解和积累程度凋落物下,土壤的养分组成和腐殖质性质都具有明显的不同。因此,积累较厚且分解程度高的凋落物层能有效减弱土壤淋溶作用,利于氮磷富集和有机质积累,土壤肥力得到改善[31],有利于植被生长和群落的稳定;反之,若林地内凋落物层稀薄,则导致土壤肥力较低,不利于森林群落的持续发展。
-
本研究综合考虑森林群落自身特征和土壤肥力,应用模糊综合评判的隶属函数值法,对辽宁仙人洞国家级自然保护区9个森林群落的稳定性进行了综合评价,取得以下结论:①辽宁仙人洞国家级自然保护区的9个典型森林群落稳定性差异显著,由高到低依次为赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落、槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落、麻栎+栓皮栎-青花椒+榛-矮丛薹草群落、赤松-三桠乌药+照山白-矮丛薹草群落、赤松-青花椒-矮丛薹草群落、蒙古栎+赤松-牛叠肚+照山白-矮丛薹草群落、蒙古栎-牛叠肚-矮丛薹草群落、赤松-牛叠肚+青花椒-矮丛薹草群落、槲树+赤松-榛-矮丛薹草群落。②种群更新潜力影响群落的恢复力,物种多样性提高了群落的抵抗力,土壤肥力为群落发展提供良好基础。三者是影响森林群落稳定性的主要因素。
Evaluation of forest community stability in Xianrendong National Nature Reserve, Liaoning
-
摘要:
目的 研究辽宁仙人洞国家级自然保护区9个典型森林群落稳定性,探讨其主要影响因素。 方法 以保护区9个典型森林群落的样地调查数据为数据源,比较各森林群落种群更新潜力、群落基本特征、物种多样性、种群生态位重叠、凋落物特性和土壤肥力差异,以此6项指标共44个因子构建群落稳定性评价体系,应用模糊数学中隶属函数值法对各森林群落稳定性状况作出评价。 结果 ①赤松Pinus densiflora +麻栎Quercus acutissima-地锦Parthenocissus tricuspidata +青花椒Zanthoxylum schinifolium-两型豆Amphicarpaea edgeworthii +白花碎米荠Cardamine leucantha群落种群更新潜力、物种多样性和土壤肥力最好,赤松-牛叠肚Rubus crataegifolius +青花椒-矮丛薹草Carex callitrichos var. nana群落生长发育基本特征平均值最高。9个群落种群生态位重叠情况相似,槲树Q. dentata +麻栎-花木蓝Indigofera kirilowii +榛Corylus heterophylla-菴闾Artemisia keiskeana +矮丛薹草群落凋落物层蓄积量和持水性最高。②9个森林群落的稳定性排名由高到低依次为赤松+麻栎-地锦+青花椒-两型豆+白花碎米荠群落、槲树+麻栎-花木蓝+榛-菴闾+矮丛薹草群落、麻栎+栓皮栎Q. variabilis-青花椒+榛-矮丛薹草群落、赤松-三桠乌药Lindera obtusiloba +照山白Rhododendron micranthum-矮丛薹草群落、赤松-青花椒-矮丛薹草群落、蒙古栎Q. mongolica +赤松-牛叠肚+照山白-矮丛薹草群落、蒙古栎-牛叠肚-矮丛薹草群落、赤松-牛叠肚+青花椒-矮丛薹草群落、槲树+赤松-榛-矮丛薹草群落。③6项森林群落稳定性评价指标中,第1主成分中种群更新潜力、物种多样性和土壤肥力的荷载系数最大。 结论 辽宁仙人洞国家级自然保护区9个典型森林群落稳定性存在显著差异,其主要影响因素为森林群落种群更新潜力、物种多样性以及土壤肥力。表10参31 Abstract:Objective The objective of this study is to evaluate the stability of 9 typical forest communities in Xianrendong National Nature Reserve, Liaoning Province, and to analyze the main influencing factors. Method Taking the survey data of 9 typical forest communities in the reserve as the data source, the population regenerative potential, basic community characteristics, species diversity, population niche overlap, litter characteristics and soil fertility of 9 forest communities were compared, and a stability evaluation system was constructed based on the above 6 indicators and 44 factors. The subordinate function value method in fuzzy mathematic was used to evaluate the stability of each forest community. Result (1) Pinus densiflora+Quercus acutissima-Parthenocissus tricuspidata+Zanthoxylum schinifolium-Amphicarpaea edgeworthii+Cardamine leucantha had the best regeneration potential, species diversity and soil fertility. The average value of basic characteristics of growth and development of P. densiflora-Rubus crataegifolius+Z. schinifolium-Carex callitrichos var. nana community was the highest. The niche overlap of the 9 communities was similar, and the litter accumulation and water holding capacity of Q. dentata+Q. acutissima-Indigofera kirilowii+Corylus heterophylla-Artemisia keiskeana+C. callitrichos var. nana community were the best. (2) The stability ranking of the 9 forest communities from high to low was P. densiflora+Q. acutissima-P. tricuspidata+Z. schinifolium-A. edgeworthii+C. leucantha, Q. dentata+Q. acutissima-I. kirilowii+C. heterophylla-A. keiskeana+C. callitrichos var. nana, Q. acutissima+Q. variabilis-Z. schinifolium+C. heterophylla-C. callitrichos var. nana, P. densiflora-Lindera obtusiloba+Rhododendron micranthum-C. callitrichos var. nana, P. densiflora-Z. schinifolium-C. callitrichos var. nana, Q. mongolica+P. densiflora-R. crataegifolius+R. micranthum-C. callitrichos var. nana, Q. mongolica-R. crataegifolius-C. callitrichos var. nana, P. densiflora-R. crataegifolius+Z. schinifolium-C. callitrichos var. nana and Q. dentata+P. densiflora-C. heterophylla-C. callitrichos var. nana. (3) Among the 6 forest community stability indicators, the load coefficient of population regenerative potential, species diversity and soil fertility was the largest in the first principal component. Conclusion There are significant differences in the stability of the 9 typical forest communities in Xianrendong National Nature Reserve, and the main influencing factors are regenerative potential, species diversity and soil fertility. [Ch, 10 tab. 31 ref.] -
毛竹Phyllostachys edulis是中国面积最广、经济价值最高的竹种[1],但其主要分布区土壤有效磷含量偏低[2],严重影响毛竹的生长和产量。研究表明:低磷胁迫降低了毛竹叶片叶绿素和金属含量[3],导致磷组分[4]和内源激素含量[5]等也发生变化,进而降低了整体生物量。常用施磷肥来增加土壤磷有效性,但过度施磷会引起水体富营养化、磷矿资源枯竭等问题[6]。因此,为应对缺磷带来的生产危机,研究植物对低磷胁迫的响应机制显得尤为重要[7]。低磷胁迫下,植物根系会分泌化学物质,影响根际土壤养分有效性[8];同时毛竹等植物会调整根冠比[5,9]等生物量分配格局来适应低磷环境。低磷胁迫还会促使植物改变自身对养分的需求,协调植物体内养分平衡[10]。但不同生长时期植物对养分的需求和获取存在差异[11]。快速生长时期,植物对养分需求较大,土壤供应可能不足,会促使植物发挥养分重吸收作用[12]。LU等[13]研究表明:低磷胁迫使紫花苜蓿Medicago sativa的磷素利用效率增强,但其效率随育苗时间延长而降低。张文元[14]研究表明:不同生长阶段毛竹对养分的需求及所处根际土壤化学性质存在差异。低磷胁迫会打破土壤原有的养分平衡,使植物体内养分调控机制发生变化。但施磷肥后不同生长时期土壤养分动态变化及毛竹实生苗适应性响应的研究尚无报道。综上,本研究模拟不同土壤磷环境,旨在探究低磷胁迫对不同生长时期根际土壤养分环境、毛竹幼苗生长和养分生理的影响及其持续效应,为植物对低磷环境的适应机制提供理论基础。
1. 材料与方法
1.1 材料
试验地设在浙江农林大学果木园大棚。毛竹种子采自广西桂林灵川县,用质量浓度为0.3%高锰酸钾消毒和无菌水冲洗并浸种24 h后催芽。供试土壤挖取于浙江省杭州市临安区青山湖仙人洞毛竹林地1 m以下土层,自然风干过20目筛后,以V(土壤)∶V(珍珠岩)=3∶1混匀成育苗基质。土壤理化性质:有机质5.02 g·kg−1,总氮0.46 g·kg−1,全磷0.24 g·kg−1,全钾13.01 g·kg−1,碱解氮93.63 mg·kg−1,有效磷2.58 mg·kg−1,速效钾54.99 mg·kg−1,pH 4.83。育苗容器采用高14.0 cm、底径11.8 cm、上口径14.0 cm的BN170型PP5塑料育苗软盆,盆底有16个透水孔,且各配套1个塑料托盘。
1.2 方法
以毛朝明等[2]对松阳县毛竹林土壤养分分级为依据,采用单因素随机区组试验设计,共设4个初始土壤有效磷水平,分别模拟极低磷(2.5 mg·kg−1,P1)、低磷(5.0 mg·kg−1,P2)、中磷(10.0 mg·kg−1,P3)和适磷(20.0 mg·kg−1,P4)土壤环境。每盆装基质1.4 kg,磷肥采用磷酸二氢钾(分析纯),溶解于水后施入基质并均匀搅拌,每盆浇施100 mL,其质量浓度分别为0、153.5、460.6和1 074.8 mg·L−1,折合每盆分别添加磷为0、3.5、10.5和24.5 mg。本底氮肥已足够支撑幼苗生长所需,不再额外施加氮肥。对于低磷土壤,施加分析纯硫酸钾以补齐钾元素,保证单一磷变量。2019年6月14日,每盆播8粒饱满健康的毛竹种子,长出4片真叶后间苗,留2株健壮且长势一致的幼苗,每处理设4个重复,每重复150株苗木,共2 400株。进行常规栽培管理,为减少边缘效应,隔半月随机移动苗木1次,并在培育过程中摒弃异常植株。
1.3 样品采集和测定
根据毛竹幼苗苗高增长变化及普遍生长节律[15],于2019年12月中旬(当年生长季末,T1)和2020年6月中旬(翌年快速生长期,T2)进行破坏性取样,分为根、茎、叶3个部分在105 ℃烘箱中杀青30 min后,70 ℃烘至恒量,测定生物量后,粉碎并过100目筛,用于测定养分。测试样品采用H2SO4-H2O2消煮[16],分别采用凯氏定氮法、钼锑抗比色法和火焰光度法测定全氮、全磷、全钾质量分数,并计算毛竹幼苗根冠比、养分积累量、养分分配比例和磷素利用效率[17]。每个指标重复测定4次,每次重复分别随机取6株均匀混合样。
同时将附着在根系上的土壤作为根际土壤。采用pH计电位法测定根际土壤pH,高温外热重铬酸钾-容量法测定有机质质量分数,H2SO4消煮-凯氏定氮法测定全氮质量分数,乙酸铵浸提法测定速效钾质量分数,盐酸-氟化铵-钼锑抗比色法测定有效磷质量分数[16],并计算土壤有效磷变化幅度。有效磷变化幅度=(某时期有效磷质量分数−初始有效磷质量分数)/初始有效磷质量分数×100%。
1.4 数据处理与分析
采用Excel 2010和SPSS 19.0统计和分析数据,SigmaPlot 12.5作图。采用单因素方差分析(one-way ANOVA),Duncan法在α=0.05水平上进行多重比较。
2. 结果与分析
2.1 低磷胁迫对根际土壤养分环境的影响
表1所示:T1时期,与初始有效磷水平相比,P3和P4处理下的根际土壤有效磷质量分数降低,但低磷处理组(P1和P2)处理下反而增加;此时低磷处理组的根际土壤pH分别较P4显著降低了10.60%和8.55% (P<0.05);但根际土壤全氮质量分数随初始有效磷水平降低而增加,其中P1较P4显著提高了19.23% (P<0.05);此时根际土壤有机质和速效钾质量分数在不同土壤有效磷水平间均无显著差异(P>0.05)。T2时期,根际土壤有效磷质量分数与初始设置水平相比均降低,但低磷处理组降幅较P4显著减小(P<0.05);此时根际土壤pH仅在P1下较P4显著降低3.53% (P<0.05);而低磷处理组的根际土壤全氮质量分数分别较P4显著增加了30.76%和11.82% (P<0.05);此时根际土壤速效钾质量分数在不同有效磷水平间仍无显著差异(P>0.05),但低磷处理组的根际土壤有机质质量分数分别较P4均显著增加了10.70% (P<0.05)。
表 1 初始土壤有效磷水平对根际土壤养分环境的影响Table 1 Effects of initial soil available phosphorus level on soil nutrient environment in rhizosphere生长时期 初始土壤有效磷水平 有效磷/(mg·kg−1) 有效磷变化幅度/% pH 有机质/(g·kg−1) 全氮/(mg·kg−1) 速效钾/(mg·kg−1) T1 P1 3.05±0.11 d 22.13±4.56 a 4.81±0.07 b 5.53±0.33 a 902.26±37.64 a 65.47±0.89 a P2 6.57±0.04 c 31.47±0.75 a 4.92±0.03 b 5.90±0.31 a 805.16±38.17 ab 65.94±1.63 a P3 8.47±0.21 b −15.30±2.09 b 5.40±0.03 a 5.83±0.42 a 868.47±4.46 ab 64.94±1.13 a P4 14.98±0.87 a −25.10±4.33 b 5.38±0.04 a 5.17±0.07 a 756.77±42.93 b 63.58±0.64 a T2 P1 2.32±0.04 D −7.33±1.76 A 5.19±0.03 B 7.45±0.21 A 1013.74±29.97 A 93.92±3.87 A P2 4.34±0.12 C −13.27±2.48 A 5.42±0.04 A 7.45±0.19 A 866.91±24.09 B 82.20±1.46 A P3 7.53±0.30 B −24.73±2.99 B 5.43±0.01 A 7.21±0.03 AB 836.34±17.01 BC 75.43±2.63 A P4 12.16±0.28 A −39.20±1.39 C 5.38±0.03 A 6.73±0.23 B 775.27±6.81 C 78.91±10.97 A 说明:数据为平均值±标准误。大小写字母分别表示T2和T1时期不同处理间差异显著(P<0.05);有效磷变化幅度中的−表示土壤有效 磷质量分数与初始水平相比降低,否则表示增加 2.2 低磷胁迫对毛竹幼苗生物量及共分配的影响
表2所示:T1时,低磷处理组(P1和P2)的毛竹幼苗单株生物量分别较P4显著降低了27.09%和15.93% (P<0.05);根冠比较P4分别显著降低22.22%和14.29% (P<0.05);但各土壤有效磷水平下的毛竹幼苗生物量分配率从大到小均依次为根、叶、茎。T2时,低磷处理组的单株生物量较P4分别显著降低了48.96%和43.40% (P<0.05);但根冠比较P4分别显著增加了44.30%和37.97% (P<0.05);此时毛竹幼苗生物量分配率在各土壤有效磷水平下均从大到小依次为根、茎、叶。
表 2 初始土壤有效磷水平对毛竹幼苗生物量及分配的影响Table 2 Effects of initial soil available phosphorus level on biomass and allocation of Ph. edulis seedlings生长时期 初始土壤
有效磷水平生物量/(g·株−1) 生物量分配/% 根冠比 根 茎 叶 T1 P1 0.197±0.005 c 49.40±1.19 d 23.35±0.74 a 27.25±0.98 a 0.98±0.05 c P2 0.227±0.003 b 52.00±0.54 c 20.14±0.04 b 27.86±0.58 a 1.08±0.02 c P3 0.264±0.005 a 60.03±0.43 a 17.00±0.42 c 22.97±0.04 b 1.50±0.03 a P4 0.270±0.009 a 55.82±0.65 b 19.50±0.46 b 24.68±0.24 b 1.26±0.03 b T2 P1 0.243±0.020 B 54.56±1.39 A 29.17±2.40 AB 16.27±1.05 B 1.14±0.10 A P2 0.270±0.025 B 52.09±1.33 A 32.06±0.67 A 15.85±1.12 B 1.09±0.06 A P3 0.307±0.012 B 53.34±1.05 A 25.02±0.83 B 21.64±1.31 A 1.14±0.05 A P4 0.477±0.050 A 43.97±1.35 B 32.37±1.43 A 23.66±0.61 A 0.79±0.04 B 说明:数据为平均值±标准误。大小写字母分别表示T2和T1时期不同处理间差异显著(P<0.05) 2.3 低磷胁迫对毛竹幼苗养分生理的影响
2.3.1 毛竹幼苗养分质量分数
由图1可见:T1时,低磷处理组(P1和P2)较P4均显著降低了毛竹各组织氮质量分数(P<0.05),也显著降低了茎和叶的磷质量分数(P<0.05),但根的磷质量分数仅在P1处理下较P4显著降低了17.01% (P<0.05);此时P1和P4处理的毛竹各组织钾质量分数均无显著差异(P>0.05),但P2处理根和叶的钾质量分数较P4显著降低了6.18%和8.05% (P<0.05)。T2时,毛竹各组织氮质量分数也均随土壤有效磷水平降低而减小,其中低磷处理显著降低了根和茎的氮质量分数(P<0.05),但叶的氮质量分数在不同土壤有效磷水平间无显著差异(P>0.05);此时低磷处理显著降低了根磷质量分数(P<0.05),叶中则显著增加(P<0.05),但茎上无显著差异(P>0.05);而低磷处理组较P4均显著降低了该时期各组织钾质量分数(P<0.05)。
2.3.2 毛竹幼苗养分积累和分配
由表3显示:无论T1还是T2时期,低磷处理(P1和P2)均显著降低了毛竹幼苗各组织和整株氮、磷、钾积累量(P<0.05)。其中T1时,P1和P2处理的氮积累量较P4分别显著降低了50.28%和38.74% (P<0.05),磷积累量则分别显著降低了42.59%和28.70% (P<0.05),钾积累量则分别显著减小了27.10%和21.55% (P<0.05);T2时,P1和P2处理的氮、磷、钾养分积累量分别依次较P4显著降低了66.86%和57.03%、57.11%和43.16%、60.94%和60.42% (P<0.05)。
表 3 初始土壤有效磷水平对毛竹幼苗养分积累量的影响Table 3 Effects of initial soil available phosphorus level on nutrient accumulation of Ph. edulis seedlings生长
时期初始土壤有效磷水平 根养分积累量/(mg·株−1) 茎养分积累量/(mg·株−1) 氮 磷 钾 氮 磷 钾 T1 P1 0.74±0.08 b 0.025±0.001 b 2.53±0.07 b 0.31±0.02 b 0.013±0.001 b 0.89±0.03 b P2 0.91±0.07 b 0.033±0.001 b 2.82±0.06 b 0.33±0.03 b 0.015±0.000 b 0.87±0.01 bc P3 1.50±0.08 a 0.046±0.004 a 3.50±0.09 a 0.34±0.03 b 0.015±0.001 b 0.78±0.01 c P4 1.60±0.18 a 0.047±0.003 a 3.85±0.19 a 0.63±0.01 a 0.021±0.001 a 1.01±0.05 a T2 P1 1.57±0.18 B 0.073±0.009 B 2.62±0.26 B 0.89±0.02 C 0.048±0.002 B 1.49±0.03 B P2 1.81±0.06 B 0.093±0.004 B 2.61±0.18 B 1.50±0.19 B 0.065±0.008 B 1.68±0.17 B P3 2.01±0.23 B 0.133±0.003 A 2.79±0.06 B 1.66±0.05 B 0.066±0.003 B 1.48±0.06 B P4 3.36±0.42 A 0.160±0.016 A 5.08±0.62 A 3.59±0.11 A 0.109±0.015 A 4.03±0.32 A 生长
时期初始土壤有效磷水平 叶养分积累量/(mg·株−1) 整株养分积累量/(mg·株−1) 氮 磷 钾 氮 磷 钾 T1 P1 0.77±0.02 c 0.024±0.001 c 1.18±0.02 b 1.81±0.10 c 0.062±0.001 d 4.60±0.09 c P2 0.99±0.11 bc 0.029±0.001 b 1.26±0.04 b 2.23±0.19 c 0.077±0.001 c 4.95±0.08 c P3 1.16±0.07 b 0.031±0.001 b 1.29±0.04 b 2.99±0.08 b 0.092±0.004 b 5.57±0.13 b P4 1.41±0.02 a 0.039±0.001 a 1.45±0.06 a 3.64±0.20 a 0.108±0.003 a 6.31±0.28 a T2 P1 1.06±0.08 C 0.042±0.006 B 1.14±0.16 C 3.51±0.24 C 0.163±0.016 C 5.25±0.39 B P2 1.24±0.26 BC 0.059±0.010 B 1.03±0.17 C 4.55±0.46 BC 0.216±0.022 BC 5.32±0.51 B P3 2.02±0.31 B 0.058±0.005 B 2.26±0.24 B 5.69±0.57 B 0.257±0.003 B 6.54±0.33 B P4 3.65±0.38 A 0.111±0.014 A 4.33±0.54 A 10.59±0.86 A 0.380±0.045 A 13.44±1.43 A 说明:数据为平均值±标准误。大小写字母分别表示T2和T1时期不同处理间差异显著(P<0.05) 由图2显示:T1时期,氮养分主要分配在根和叶上,分别占40%,且低磷处理组(P1和P2)较P3显著降低了氮在根的分配比例(P<0.05);此时磷养分在各组织分配比例从大到小均依次为根、叶、茎,且低磷处理组也较P3显著降低了磷分配在根的比例(P<0.05);钾养分与磷相似,分配比例从大到小均依次为根、叶、茎,此时低磷处理组较P3和P4均显著降低了其在根的分配比例(P<0.05)。T2时期,氮分配在各组织比例相当,且低磷处理组的氮分配在根的比例较P4均分别显著增加了41.05%和28.85% (P<0.05);钾养分分配比例因土壤有效磷水平不同而异,其中低磷处理组的钾分配比例从大到小依次为根、茎、叶,但中适磷则为根、叶、茎,此时低磷处理组的钾养分分配在根的比例较P4分别显著增加了31.79%和30.86% (P<0.05)。
2.3.3 毛竹幼苗磷素利用效率
由表4可见:T1时期,低磷处理(P1和P2)显著增加了毛竹茎、叶及整株磷素利用效率,但根的磷素利用效率仅在P1下较P4显著增加了20.31%(P<0.05)。T2时期,低磷处理显著增加了根的磷素利用效率(P<0.05),叶中则显著降低(P<0.05),茎上无显著差异(P>0.05);此时仅P1的整株磷素利用效率较P4显著增加了19.05%,P2与P4间无显著差异(P>0.05)。
表 4 初始土壤有效磷水平对毛竹幼苗磷素利用效率的影响Table 4 Effects of initial soil available phosphorus level on phosphorus utilization efficiency of Ph. edulis seedlings生长
时期初始土
壤有效
磷水平磷素利用效率/(g·mg−1) 根 茎 叶 整株 T1 P1 3.85±0.09 a 3.57±0.07 a 2.22±0.03 a 3.16±0.02 a P2 3.58±0.15 ab 3.07±0.11 b 2.15±0.08 ab 2.93±0.07 b P3 3.44±0.21 ab 2.97±0.08 b 1.98±0.06 b 2.86±0.07 b P4 3.20±0.09 b 2.48±0.15 c 1.71±0.03 c 2.50±0.02 c T2 P1 1.83±0.03 a 1.46±0.06 a 0.95±0.01 c 1.50±0.03 a P2 1.51±0.05 b 1.35±0.11 a 0.74±0.01 d 1.25±0.04 b P3 1.23±0.04 c 1.17±0.07 a 1.15±0.02 a 1.19±0.04 b P4 1.31±0.02 c 1.43±0.09 a 1.02±0.01 b 1.26±0.02 b 说明:磷素利用效率=生物量/磷素积累量;数据为平均值±标准误;大小写字母分别表示T2和T1时期不同处理间差异显著(P<0.05) 3. 讨论
3.1 低磷胁迫对毛竹幼苗根际土壤的影响机制
根际土壤环境直接影响土壤养分的转化、迁移及被植物吸收的能力[18],根际土壤有效磷作为表征土壤供磷能力的重要指标,是植物吸收磷素的直接形态[19]。本研究中,低磷胁迫显著降低了当年生长季末根际土壤pH,但其有效磷质量分数较初始设置水平增加,这可能是磷素缺乏时,毛竹根系会分泌苹果酸、草酸等有机酸物质进入土壤[20],增强根际酸化能力,并与土壤中磷酸铁、磷酸铝等难溶性磷酸盐作用,释放出磷酸盐,可相对提高磷的有效性;此时根际土壤全氮质量分数随有效磷水平降低而增加,但有机质和速效钾质量分数无显著差异,这表明当年生长季末时,极低磷胁迫会相对更阻碍毛竹幼苗吸收根际土壤氮素,抑制幼苗生长发育。随育苗时间推移,低磷胁迫显著增加了翌年快速生长期的根际土壤有机质质量分数,但有效磷降幅较小,这可能是由于低磷胁迫下的毛竹根系会随着生长分泌更多有机物质,缓解难溶性磷酸盐在土壤中的沉积,维持土壤磷肥力。此时,低磷环境均会显著限制毛竹幼苗吸收土壤氮素。
3.2 低磷胁迫对毛竹幼苗生长和养分生理的影响机制
低磷胁迫显著抑制了当年生长季末毛竹幼苗生物量的积累,这可能是由于磷素不足会导致植物体内核酸、核蛋白合成受阻,抑制细胞形成和分裂[21]。植物养分平衡作为生物量积累的基础,其养分浓度和积累量受土壤养分有效性影响[22]。本研究中,养分浓度与积累量随土壤磷水平降低变化规律并不完全一致,主要原因是各组织养分积累量和生物量均随土壤磷水平降低而减小,但两者降幅不完全同步,从而对某些养分浓度造成相对稀释或浓缩[23],这是植物养分协调供应的表现形式之一。养分利用效率反映了植物对养分的利用和适应情况,本研究表明:低磷胁迫显著提高了当年生长季末毛竹整体磷素利用效率,与目前主流观点一致:缺磷素环境下的植物养分利用效率高于富态营养环境[23]。低磷胁迫显著降低了当年生长季末毛竹幼苗根冠比,也相对减少了氮、磷、钾养分在根系的分配比例,这表明生长季末时,低磷胁迫会提高毛竹叶片储存养分的比例,增强光合作用,并尽可能多的制造有机物质,在来年春季新生长时,作为源将养分输送到其他组织促进生长发育。
3.3 低磷胁迫对毛竹幼苗生长和养分影响的持续效应
低磷胁迫显著降低毛竹生物量和养分积累量的作用持续到翌年快速生长期,且降幅较当年生长季末更高,这表明了低磷胁迫对毛竹幼苗吸收根际土壤养分的抑制作用随育苗时间延长有所增强,并进一步抑制了生物量的积累。在翌年快速生长期,低磷胁迫显著降低了根的磷质量分数,但叶相反,这可能是毛竹幼苗为应对低磷胁迫采用的“被动”策略,即通过活化细胞壁中所储存的磷素[24],促进磷素内循环,从而提高叶片光合效率[25],以维持正常生理生化活动所需。毛竹整体磷素利用效率至翌年快速生长期时仅在极低磷胁迫下显著增加,这表明毛竹幼苗对低磷胁迫的适应性随育苗时间延长而增强。此外,至翌年快速生长期,低磷胁迫显著提高了毛竹幼苗根冠比和氮、钾养分分配到根系的比例,这可能是由于低磷土壤迫使快速生长期的毛竹幼苗分配更多养分用于支持根系生长,提高根系活力;也可能是当年生长季末储存在叶片中的养分被更多地转运到根系,以维持苗木正常生长发育。
4. 结论
低磷胁迫降低了当年生长季末毛竹幼苗根际土壤pH,但维持了根际土壤高氮质量分数,这种作用持续到了翌年快速生长期,此时低磷胁迫还显著增加了根际土壤有机质质量分数。此外,低磷胁迫显著抑制了毛竹生物量和养分的积累,但提高了整株磷素利用效率;低磷抑制作用随育苗时间推移增强,但毛竹幼苗会通过提高根冠比、根系养分分配比例来提高对低磷胁迫的适应性。
-
表 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 表 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 表 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 表 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), 相同字母表示差异不显著 表 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表示群落内种群生态位重叠值在此范围的种对 占全部种对的比例 表 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) 表 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) 表 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 表 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 表 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 -
[1] 刘胜涛, 高鹏, 刘潘伟, 等. 泰山森林生态系统服务功能及其价值评估[J]. 生态学报, 2017, 37(10): 3302 − 3310. LIU Shengtao, GAO Peng, LIU Panwei, et al. An ecosystem services assessment of Tai Mountain [J]. Acta Ecol Sin, 2017, 37(10): 3302 − 3310. [2] 于晓文, 宋小帅, 康峰峰, 等. 冀北辽河源典型森林群落稳定性评价[J]. 干旱区资源与环境, 2015, 29(5): 93 − 98. YU Xiaowen, SONG Xiaoshuai, KANG Fengfeng, et al. Evaluation on the stability of typical forest communities in sourse region of Liaohe River in North Hebei [J]. J Arid Land Resour Environ, 2015, 29(5): 93 − 98. [3] 高东, 何霞红. 生物多样性与生态系统稳定性研究进展[J]. 生态学杂志, 2010, 29(12): 2507 − 2513. GAO Dong, HE Xiahong. Research advances on biodiversity and ecosystem stability [J]. Chin J Ecol, 2010, 29(12): 2507 − 2513. [4] GRIMM V, SCHMIDT E, WISSEL C. On the application of stability concepts in ecology [J]. Ecol Modell, 1992, 63(1/4): 143 − 161. [5] 杨瑞红, 赵成义, 王新军, 等. 固沙植物群落稳定性研究[J]. 中南林业科技大学学报, 2015, 35(11): 128 − 135. YANG Ruihong, ZHAO Chengyi, WANG Xinjun, et al. Research progress of stability of sand-binding plant communities [J]. J Cent South Univ For Technol, 2015, 35(11): 128 − 135. [6] 郑元润. 森林群落稳定性研究方法初探[J]. 林业科学, 2000, 36(5): 28 − 32. ZHENG Yuanrun. Comparison of methods for studying stability of forest community [J]. Sci Silv Sin, 2000, 36(5): 28 − 32. [7] 陆龙龙, 郭忠玲, 范春楠, 等. 吉林磨盘山次生落叶阔叶林群落特征和稳定性分析[J]. 应用生态学报, 2018, 29(7): 2079 − 2087. LU Longlong, GUO Zhongling, FAN Chunnan, et al. Community characteristics and stability analysis of secondary deciduous broad-leaved forest in Mopan Mountains, Jilin Province, China [J]. Chin J Appl Ecol, 2018, 29(7): 2079 − 2087. [8] 李伟成, 盛海燕, 金孝锋, 等. 浙江铜铃山国家森林公园次生常绿阔叶林和针阔混交林乔木群落特征研究[J]. 植物资源与环境学报, 2019, 28(3): 66 − 77. LI Weicheng, SHENG Haiyan, JIN Xiaofeng, et al. Study on arbor community characteristics of secondary evergreen broad-leaved forest and coniferous and broad-leaved mixed forest in Tonglingshan National Forest Park of Zhejiang Province [J]. J Plant Resour Environ, 2019, 28(3): 66 − 77. [9] 张明霞, 王得祥, 彭舜磊, 等. 秦岭松栎混交林群落的稳定性[J]. 生态学报, 2015, 35(8): 2564 − 2573. ZHANG Mingxia, WANG Dexiang, PENG Shunlei, et al. Community stability analysis for the oak-pine mixed forest in Qinling Mountains [J]. Acta Ecol Sin, 2015, 35(8): 2564 − 2573. [10] 郭其强, 张文辉, 曹旭平. 基于模糊综合评判的森林群落稳定性评价体系模型构建——以黄龙山主要森林群落为例[J]. 林业科学, 2009, 45(10): 19 − 24. GUO Qiqiang, ZHANG Wenhui, CAO Xuping. Establishment of an evaluation model of the forest community stability based on fuzzy synthetic evaluation: a case study of main forest communities in Huanglong Mountains [J]. Sci Silv Sin, 2009, 45(10): 19 − 24. [11] 康希睿, 李晓刚, 张涵丹, 等. 不同混交措施下杉木人工林群落稳定性特征[J]. 生态学杂志, 2020, 39(9): 2912 − 2920. KANG Xirui, LI Xiaogang, ZHANG Handan, et al. Community stability characteristics of Cunninghamia lanceolata plantations with different mixing measures [J]. Chin J Ecol, 2020, 39(9): 2912 − 2920. [12] 任平, 洪步庭, 程武学, 等. 长江上游森林生态系统稳定性评价与空间分异特征[J]. 地理研究, 2013, 32(6): 1017 − 1024. REN Ping, HONG Buting, CHENG Wuxue, et al. Stability evaluation of forest ecosystem and study of spatial differential features in the upper Yangtze River [J]. Geogr Res, 2013, 32(6): 1017 − 1024. [13] 周超凡, 王博恒, 王蔚炜, 等. 基于空间结构指数的黄龙山林区不同森林群落稳定性评价[J]. 中南林业科技大学学报, 2018, 38(7): 76 − 82. ZHOU Chaofan, WANG Boheng, WANG Weiwei, et al. Stability evaluation of different forest communities in Huanglongshan mountain based on the spatial structure index [J]. J Cent South Univ For Technol, 2018, 38(7): 76 − 82. [14] 王鹏, 邱英杰. 辽宁仙人洞国家级自然保护区科学考察集[M]. 北京: 中国林业出版社, 2002: 1 − 5. WANG Peng, QIU Yingjie. Scientific Investigation in Xianrendong National Nature Reserve, Liaohing [M]. Beijing: China Forestry Publishing House, 2002: 1 − 5. [15] 张梅. 辽东半岛仙人洞自然保护区植物区系多样性分析[J]. 西北植物学报, 2014, 34(8): 1693 − 1698. ZHANG Mei. Plant floristic diversity in Xianrendong Natural Reserve from Liaodong Peninsula [J]. Acta Bot Boreal-Occident Sin, 2014, 34(8): 1693 − 1698. [16] 张乐, 孙向阳, 李素艳, 等. 辽宁仙人洞典型林分森林土壤碳氮分布特征[J]. 吉林农业大学学报, 2017, 39(2): 183 − 188. ZHANG Le, SUN Xiangyang, LI Suyan, et al. Characteristics of carbon and nitrogen distribution across typical forest types in Xianrendong Nature Reserve of Liaoning Province [J]. J Jilin Agric Univ, 2017, 39(2): 183 − 188. [17] 张维伟, 薛文艳, 杨斌, 等. 桥山栎林群落结构特征与物种多样性相关关系分析[J]. 生态学报, 2019, 39(11): 3991 − 4001. ZHANG Weiwei, XUE Wenyan, YANG Bin, et al. Differences in community characteristics, species diversity, and their coupling associations among three Quercus acutissima forests in Qiaoshan [J]. Acta Ecol Sin, 2019, 39(11): 3991 − 4001. [18] 李楚均, 陈小梅, 温小浩, 等. 城市化背景下珠江三角洲常绿阔叶林群落结构及植物多样性[J]. 生态学杂志, 2019, 38(11): 3298 − 3305. LI Chujun, CHEN Xiaomei, WEN Xiaohao, et al. Effects of urbanization on community structure and plant diversity in evergreen broad- leaved forests of the Pearl River Delta [J]. Chin J Ecol, 2019, 38(11): 3298 − 3305. [19] 伏捷, 孙才志, 张华, 等. 辽东山地老秃顶子古冰缘地貌植物种群生态位特征[J]. 地理研究, 2018, 37(4): 731 − 741. FU Jie, SUN Caizhi, ZHANG Hua, et al. Niche characteristics of plant populations in the paleo-periglacial areas of Mt. Laotudingzi, eastern Liaoning mountainous region [J]. Geogr Res, 2018, 37(4): 731 − 741. [20] 王志泰, 马瑞, 马彦军, 等. 利用隶属函数法分析胡枝子抗旱性[J]. 干旱区资源与环境, 2013, 27(9): 119 − 123. WANG Zhitai, MA Rui, MA Yanjun, et al. The drought resistance of 5 species of Lespedeza [J]. J Arid Land Resour Environ, 2013, 27(9): 119 − 123. [21] 孙欧文, 蔡建国, 吴家森, 等. 浙江省典型森林类型枯落物及林下土壤水文特性[J]. 水土保持研究, 2019, 26(1): 118 − 123. SUN Ouwen, CAI Jianguo, WU Jiasen, et al. Hydrological characteristics of litter and foret soil of typical forest types in Zhejiang Province [J]. Res Soil Water Conserv, 2019, 26(1): 118 − 123. [22] 谢立红, 曹宏杰, 黄庆阳, 等. 五大连池火山森林群落多样性与稳定性[J]. 浙江农林大学学报, 2021, 38(2): 235 − 245. XIE Lihong, CAO Hongjie, HUANG Qingyang, et al. On the diversity and stability of forest communities in Wudalianchi Volcanoes [J]. J Zhejiang A&F Univ, 2021, 38(2): 235 − 245. [23] 张梦弢, 张青, 亢新刚, 等. 长白山云冷杉林不同演替阶段群落稳定性[J]. 应用生态学报, 2015, 26(6): 1609 − 1616. ZHANG Mengtao, ZHANG Qing, KANG Xingang, et al. Community stability for spruce-fir forest at different succession stages in Changbai Mountains Northeast China [J]. Chin J Appl Ecol, 2015, 26(6): 1609 − 1616. [24] 张景慧, 黄永梅. 生物多样性与稳定性机制研究进展[J]. 生态学报, 2016, 36(13): 3859 − 3870. ZHANG Jinghui, HUANG Yongmei. Biodiversity and stability mechanisms: understanding and future research [J]. Acta Ecol Sin, 2016, 36(13): 3859 − 3870. [25] 王恒方, 吕光辉, 周耀治, 等. 不同水盐梯度下功能多样性和功能冗余对荒漠植物群落稳定性的影响[J]. 生态学报, 2017, 37(23): 7928 − 7937. WANG Hengfang, LÜ Guanghui, ZHOU Yaozhi, et al. Effects of functional diversity and functional redundancy on the stability of desert plant communities under different water and salt gradients [J]. Acta Ecol Sin, 2017, 37(23): 7928 − 7937. [26] 宋启亮, 董希斌. 大兴安岭不同类型低质林群落稳定性的综合评价[J]. 林业科学, 2014, 50(6): 10 − 17. SONG Qiliang, DONG Xibin. Comprehensive evaluation of forest community stability of different types of low-quality forest stands in the Greater Higgnan Mountains [J]. Sci Silv Sin, 2014, 50(6): 10 − 17. [27] MAY R M. Will a large complex system be stable [J]. Nature, 1972, 238(5364): 413 − 414. [28] McNAUGHTON S J. Diversity and stability of ecological communities: a comment on the role of empiricism in ecology [J]. Am Nat, 1977, 111(979): 515 − 525. [29] 朱鹤健, 陈健飞, 陈松林, 等. 土壤地理学[M]. 2版. 北京: 高等教育出版社, 2010: 44 − 45. ZHU Hejian, CHEN Jianfei, CHEN Songlin, et al. Soil Greography [M]. 2nd ed. Beijing: Higher Education Press, 2010: 44 − 45. [30] 王利彦, 周国娜, 朱新玉, 等. 凋落物对土壤有机碳与微生物功能多样性的影响[J]. 生态学报, 2021, 41(7): 1 − 10. WANG Liyan, ZHOU Guona, ZHU Xinyu, et al. Effects of litter on organic carbon and microbial diversity [J]. Acta Ecol Sin, 2021, 41(7): 1 − 10. [31] 邵国栋, 艾娟娟, 孙启武, 等. 昆嵛山不同林分类型土壤质量状况及评价[J]. 林业科学研究, 2018, 31(6): 175 − 184. SHAO Guodong, AI Juanjuan, SUN Qiwu, et al. Soil quality status and assessment of different forest types in Kunyushan Mountains [J]. For Res, 2018, 31(6): 175 − 184. 期刊类型引用(5)
1. 刘娅惠,徐瑾,雷蕾,万一,曾立雄,肖文发. 不同磷质量分数下马尾松幼苗根的生理生化特征. 浙江农林大学学报. 2023(01): 126-134 . 本站查看
2. 诸炜荣,黄云峰,柴庆辉,吴佳,靳晓东,余礼华,毕毓芳,吴雪玲,王璐,唐辉. 毛竹林下不同树种的生长、光合作用及生理特性. 东北林业大学学报. 2023(05): 26-29 . 百度学术
3. 赵海民. 缺磷胁迫对旱柳生长及光合作用的影响. 河北林业科技. 2023(02): 16-20 . 百度学术
4. 高宁,邢意警,熊瑞,史文辉. 丛枝菌根真菌和溶磷细菌协调植物获取磷素的机制. 浙江农林大学学报. 2023(06): 1167-1180 . 本站查看
5. 孟超敏,耿翡翡,卿桂霞,周佳敏,张富厚,刘逢举. 陆地棉磷高效基因GhMGD3的克隆与表达分析. 浙江农林大学学报. 2022(06): 1203-1211 . 本站查看
其他类型引用(4)
-
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20210330