浙江五岭坑原生常绿阔叶林群落的动态特征
2022-09-07周荣飞吴义松蔡焕满吴敏吴素美刘金亮吴友贵
周荣飞,吴义松,蔡焕满,吴敏,吴素美,刘金亮,吴友贵
(1.钱江源-百山祖国家公园 庆元保护中心,浙江 庆元 323800;2.温州大学 生命与环境科学学院,浙江 温州 325035)
亚热带常绿阔叶林是世界上主要森林植被类型之一,主要分布在中国,以中亚热带的常绿阔叶林最为典型[1-4]。亚热带常绿阔叶林区域约占中国国土面积的1/4,其结构复杂,生物多样性丰富,生态系统服务功能强,对我国经济社会的可持续发展及生态安全有着重要作用。由于长期的人类干扰,中国亚热带常绿阔叶林绝大部分退化为次生林或被改造为人工林,老龄林或原生林几乎丧失殆尽[1,5-7]。百山祖国家公园内的五岭坑区块海拔为560~1 300.7 m,海拔跨度较大,该山地常绿阔叶林长久未受较明显的人为干扰,是目前我国罕见的保存完好的中亚热带原生常绿阔叶林,具有很高的研究价值。鉴于对原生常绿阔叶林的研究较少,研究原生常绿阔叶林的动态,有助于我们掌握原生常绿阔叶林的群落结构、动态特征、物种多样性水平及其维持机制,对于常绿阔叶林的保护和管理具有重要的意义,目前已有文献对五岭坑原生常绿阔叶林进行了长期监测研究[8-11]。
本文通过对五岭坑监测样地的数据分析,以期了解群落的木本植物死亡和补员总体情况;主要优势种群的年死亡率、年补员率、种群大小变化率、物种相对适合度;群落的胸高断面积及平均胸径的变化,为制定国家公园重要自然生态系统的原真性和完整性保护措施提供科学依据,也能为其他单位的森林管理提供参考,同时,还为森林资源二类调查常绿阔叶林小班测树因子判读以及森林资源年度监测数据更新提供参考。
1 材料与方法
1.1 研究区自然概况
研究区位于百山祖国家公园内的五岭坑区块,地理位置为119°3′10′′~119°4′38′′ E,27°32′26′′~27°34′1′′ N,面积为439.3 hm2。地处我国东南沿海的闽浙丘陵中山区,以侵蚀地貌为主,山地坡度大多在30°以上,基岩为侏罗纪火成岩,海拔800 m 以下为红壤,800 m 以上以黄壤为主。该地区主要气候类型为亚热带湿润季风气候,属“南岭闽瓯中亚热带”气候区,但由于区内相对较高的海拔且植被较为茂密,再加上其独特的地形和水文地理环境,气候具有特殊性,主要特点为温暖湿润、雨量充沛、温度低、雾日多、相对湿度大[8]。
1.2 研究方法
1.2.1 样地设置 2012 年8 月,参照CTFS(Center for Tropical Forest Science)样地建设方案,在五岭坑原生常绿阔叶林的主要分布区内,建立3 个1 hm2(100 m×100 m)的样地,分别命名为1、2、3 号样地。将3 号样地划分为400 个5 m×5 m 小样方,对所有胸径≥1 cm 的木本植物进行挂牌并记录树种名称,测量其在所在小样方(5 m×5 m)内的坐标,在距离地表1.3 m 处用围尺测量胸径,并目测高度、冠幅、生长状况等。若植株距离地表1.3 m 以下有分枝,则同样测量分枝的胸径、高度、冠幅等[8]。调查时发现,1 hm2内对角(西南和东北)的两个30 m×30 m 样方(划分为36 个5 m×5 m 小样方)的物种数占总物种数的90%,几乎可以反映整个1 hm2内的物种组成,因此,在保证取样充分的前提下,为了节约调查成本,对于1、2 号样地,仅选择在每个样地的西南角和东北角分别调查30 m×30 m 的样方,西南角样方命名为01,东北角样方命名为02(表1)。
表1 五岭坑常绿阔叶林样地自然概况Table 1 Site condition of sample plots at evergreen broad-leaved forest in Wulingkeng Mountain
2017 年7 月,按2012 年的测量方法重新测量所有的样方个体胸径、树高等指标,并调查记录新增胸径≥1 cm 的个体(补员个体),记录死亡的个体。用罗盘仪测量坡向、坡度,用鱼眼相机法测定平均郁闭度,用GPS测量样地中心海拔高度[8]。
1.2.2 数据处理与统计
(1)数据来源 木本植物采用调查获得的树种数量和胸径数据,分析群落主要树种的补员和死亡动态及种群大小变化趋势。其中,补员个体是指在第二次调查中新出现的胸径(DBH)≥1 cm 的个体;死亡个体是指在复查时找到后发现已经死亡或无法找到的第一次调查时记录到的胸径≥1 cm 的植株。
(2)相关关系回归分析 利用WPS 表格进行死亡个体胸径(2012 年初次调查时的胸径)与径级相关关系分析,计算公式和相关系数,并用F值和相关系数表进行显著性检验。
(3)功能型划分 根据树种的生长型,将其分为乔木、小乔木、灌木三类;并根据树种的生长习性,将其分为常绿和落叶两类[12]。
(4)年死亡率、年补员率、种群大小变化率、种群适合度计算 将2012 年时样地中株数大于50 株的树种(即主要树种,共42 个树种)作为研究对象,用标准化法计算年死亡率、年补员率、种群数量的变化和相对适合度。通过计算各个物种的年死亡率(M)和年补员率(R)参数来比较物种间的差异。
式中,N0和Nt分别为某一树种在第一次和第二次调查时的个体数,St是该树种在第一次和第二次调查时都存活的个体数,T为两次调查的时间间隔。种群大小变化率采用年补员率和年死亡率的差值来表示,种群适合度用年补员率和年死亡率的比值来表示[13-14]。
2 结果与分析
2.1 死亡和补员的总体情况
2012—2017 年,在5 个样方中,胸径≥1 cm 的木本植物死亡株数为1 277 株,占2012 年个体总数(10 459 株)的12.21%,隶属于93 个物种,占物种总数(142 种)的65.49%。死亡木本植物的最大胸径为44.8 cm,最小胸径为1.0 cm。样地内消失的胸径≥1 cm 的物种有5 种(垂株花Styrax dasyanthus、毛果槭Acer nikoense、苦槠Castanopsis sclerophylla、宜昌荚蒾Viburnum erosum和老鸦糊Callicarpa giraldii),占物种总数的2.08%,这5 种物种在2012年调查时皆为1 株。
将5个样方内所有死亡树木的胸径按2 cm径阶进行整化,并绘制胸径-株数分布图(图1)。从图1 中可看出,胸径为1.0~3.0 cm 的死亡树木株数占绝大多数,占死亡树木总株数的64.18%;胸径为>3.0~5.0 cm 的死亡树木株数占死亡树木总株数的20.14%;胸径为>5.0~7.0 cm 的死亡树木株数占死亡树木总株数的6.13%。幼树(胸径1.0~7.0 cm)的死亡树木总数超过90%。死亡树木的径级呈倒“J”形分布。将同径级的死亡树木株数相加,将死亡株数(N)和径级(D:cm)做回归分析,结果表明,死亡株数的自然对数与径级呈负线性相关(即死亡株数与径级呈指数负相关关系),经F检验和相关系数检验,线性回归关系均为极显著相关,其关系式可表示为:lnN=5.061-0.157D,(R2=0.769,F=43.90,P<0.01)。
图1 死亡木本植物径级分布Figure 1 Diameter class distribution of dead trees
2012—2017 年,5 个样方合计补员412 株,隶属于51 个物种,分别占2012 年个体总数和物种总数的3.94%和35.92%。补员木本植物的最大胸径为5.4 cm,最小胸径为1.0 cm。补员木本植物胸径1.0~2.0 cm、>2.0~3.0 cm、>3.0~4.0 cm、>4.0~5.0 cm、>5.0~6.0 cm 的株数分别为345、53、10、2、1 株。
2.2 不同功能型树种的死亡和补员分析
2012 年,将5 个样方中树木胸径≥1 cm 的个体按功能型分,常绿树种和落叶树种分别为10 316 株和143 株;乔木、小乔木和灌木分别为3 847、3 776和2 836 株。
2017 年不同功能型树种的死亡和补员情况见图2。由图2 可知,常绿树种和落叶树种的死亡比例分别为12.22%和11.19%,两者差别不大;灌木、小乔木和乔木的死亡比例分别为15.27%、10.75%和11.39%,其中灌木的死亡比例明显高于乔木和小乔木。
图2 死亡和补员木本植物的生活型分布Figure 2 Different life form distribution of dead and new trees
常绿树种和落叶树种的补员比例分别占3.96%和2.10%,常绿树种的补员比例高于落叶树种;灌木、小乔木和乔木的补员比例分别为6.14%、3.73%和2.52%,灌木的补员比例明显高于乔木和小乔木,小乔木的补员比例高于乔木。
2.3 主要树种的变化
胸径≥1 cm 个体数量大于50 株的有42 种(树种),其中乔木树种19 种、小乔木树种14 种、灌木树种9种。它们的年死亡率、年补员率、种群数量变化率和相对适合度见表2。年死亡率最大的为黑山山矾(9.27%),此外,年死亡率大于5%的树种还有矩叶鼠刺、乌药和短尾越桔3 种灌木。年死亡率最小的为木荷(0.52%),年死亡率低于1%的还有厚皮香、树参和石木姜子,这4 个树种均为乔木树种。主要树种的平均年死亡率按功能型分,灌木(3.46%)>小乔木(2.82%)>乔木(1.98%)。
表2 常绿阔叶林42 个主要树种的年死亡率、年补员率、年变化率和适合度Table 2 The annual mortality,new trees,population change rate and suitability of 42 main tree species in the sample plots
表2 (续)
补员率最高的为三花冬青(2.45%),此外,补员率大于1%的还有3 个小乔木树种(赤楠、凤凰润楠、细枝柃)和5 个灌木树种(具毛常绿荚蒾、乌药、窄基红褐柃、石斑木、短尾越桔)。没有补员的物种有7 种乔木树种(木荷、虎皮楠、少叶黄杞、米槠、尾叶冬青、硬壳柯、山矾)。主要树种的平均年补员率按功能型分,灌木(1.27%)>小乔木(0.70%)>乔木(0.48%)。
除了细枝柃和石木姜子的种群数量没有发生变化外,其余40 种的种群数量均在减少,其中海桐山矾种群的变化率最大,为-8.77%,种群变化率低于-4%的还有4 种,分别为杨梅叶蚊母树、少叶黄杞、米槠、矩叶鼠刺。主要树种的平均种群变化率按功能型分,灌木(-2.19%)>小乔木(-2.11%)>乔木(-1.50%)。
种群适合度最大的为细枝柃和石木姜子(1.00),种群适合度大于0.9 的还有5 种,分别为赤楠、厚皮香、三花冬青、具毛常绿荚蒾和窄基红褐柃。没有补员的7 种树种的种群适合度为0。主要树种的平均适合度按功能型分,灌木(0.42%)>小乔木(0.36%)>乔木(0.35%)。
从功能型上看,主要树种的年死亡率、年补员率、种群变化率、适合度的平均值都表现出:灌木>小乔木>乔木。
2.4 胸高断面积和平均胸径的变化
将5 个样方中胸径≥1 cm 树木的胸高断面积和平均胸径变化情况进行统计,结果见表3。由表3 可知,胸径≥1 cm 树木的胸高断面积增长比例为4.59%,其中,树木胸高断面积增长比例最大的为101 样方(6.82%)。102、201 和202 号样方的树木胸高断面积增长比例都较小。102 号样方是由于树木的枯死比例较高(19.78%),在一定程度上影响了树木的胸高断面积增长比例。201 号样方则是因为死亡树木胸径偏大,死亡个体胸径≥10 cm的有12 株。202 号样方虽然树木的死亡率较低,但样地内死亡树木有径阶26 cm、28 cm、30 cm 各1 株,仅这3 株大树的死亡,就使面积900 m2样地的胸高断面积减少了5.70%。
表3 样方活立木的胸高断面积和平均胸径变化情况Table 3 Changes of basal area and average DBH of living wood in the sample plots
5 个样方中胸径≥1 cm 树木的平均胸径增长范围为0.44~1.21 cm,平均增长0.55 cm,增长比例为6.78%。其中,树木平均胸径增长最大的为102 号样方,增长了1.21 cm,增长比例为11.53%。虽然102 号样方树木有近20%的死亡比例,但枯死木主要集中在较小的径级,小径级树木死亡多反而提高了活立木平均胸径增长比例。另外,102 号样方中的树木少,在较大死亡比例的情况下,树木之间在空间、营养等方面的竞争相对不激烈,保留木得以较快生长,这也与样方中树木平均胸径增长较快存在一定的关系。201、202 和3 号样方树木胸径的平均增长较小(0.44~0.52 cm),其中202 号样方中树木平均胸径增长仅为0.44 cm,这还是因为多株大树死亡的原因,较大地减少了“密而小”小样方树木的平均胸径。3 号样方中的死亡个体胸径≥10 cm 的树木有29 株,其中胸径最大树木达到44.8 cm(甜槠),从而较大地影响了样方树木的平均胸径。
为了给森林资源二类调查常绿阔叶林小班测树因子判读提供参考,计算5 个样方中胸径≥5 cm 的活立木胸高断面积、平均胸径等变化情况,结果见表3。由表3 可知,5 个样方中胸径≥5 cm 的活立木胸高断面积增长比例为1.04%,其中,活立木胸高断面积增长比例最大的为101 号样方(8.93%),增长比例最小的为3 号样方(0.23%),201 号样方和202 号样方活立木胸高断面积的增长比例较小(均为0.86%)。5 个样方中胸径≥5 cm的活立木的平均胸径增长范围为0.25~1.31 cm,平均增长0.47 cm,增长比例为3.12%。相较于胸径≥1 cm 的活立木,胸径≥5 cm 活立木的平均胸径增长范围更大,增长比例更小。其中,活立木平均胸径增长最大的为102号样方,增长了1.31 cm;增长最小的为3 号样方,增长了0.25 cm。
3 讨论与结论
3.1 讨论
3.1.1 原生常绿阔叶林的死亡率和补员率 群落内个体的死亡原因较为复杂,除了种群自身的生长规律外,还有种间竞争、环境因素、病虫害和人为干扰等诸多原因[15-17]。自然稀疏是指在个体密度高的植物群落中,随着群落的生长,个体密度逐渐降低,直至达到平衡的现象,自然稀疏是群落的自我调节机能之一,通过密度调控把个体适应、种群数量动态、生物群落演替与生态系统稳定性联系起来[18-20],只有通过死亡才能降低群落的密度,为存活个体提供足够的生长空间,维持生物量的增长,为了使群落获得长期的稳定性,死亡比例必须逐渐升高,并最终超过出生比例,正是受到“密度制约效应”和“非对称竞争”的影响,群落大小才会具有稳定的平衡状态[21]。群落具有较高的枯死比例,表明该群落具有较强的自我调控能力,具有一定的稳定性[18]。林冠层乔木的死亡对于群落结构的形成与维持起着重要作用,它能为某些动植物提供独特的生境,并影响群落的径级、年龄结构、林隙特征和不同层次植被的发育等[22]。
本研究中,群落具有较高的死亡率和较低的补员率,这与百山祖中山常绿阔叶林群落具有基本相同的特征[13],对于具体原因,本文没有深入研究。但这可能是暂时的,因为较高的死亡率,尤其是大树的死亡,会在原来较密的森林里形成较大的“林窗”,为周围的种子萌发、幼苗生长创造条件,必将迎来新一轮的较高补员率,只不过“林窗效应”在短期内未能显现,需要较长的时间,可能出现在第二次或第三次复查中。相反,较高的补员率,随着树木的生长,因营养、空间等可利用资源有限,“密度制约效应”显现,非对称竞争加剧,也必将有新一轮的较高死亡率。因而,可以辩证地理解,死亡率和补员率是相对的,同时也是互为条件的。补员是森林更新的基础,只有新的个体补员,才能实现森林更新。同时,死亡是森林更新的条件,只有老的个体不断死亡,才能为新的个体生长提供后续生长条件。树木的死亡和更新是自然现象,作为森林生态系统类型的自然保护区,一般不需要采取特殊的措施进行干预,但要减少人为干扰,除了科研监测人员外,严禁其他人员进入。
3.1.2 不同径级个体的死亡率 该常绿阔叶林群落死亡树木的径级呈倒“J”形分布,与活立木存在基本相同的分布[8]。幼树(胸径1.0~7.0 cm)的死亡总数超过90%,这与百山祖常绿阔叶林、武夷山米槠林、海南岛霸王岭热带季雨林在中小径级时的死亡率较高的结论相吻合[15,23-24]。这可能是因为幼树的大量存在会造成激烈的种内种间竞争,随着径级的增大,个体数量快速减少,即“活的多死的也多”。还可能是由于上层树木的优势不断的积累和放大,处于劣势的下层林木个体密度较大,光照不足,种内种间竞争强度较大,所以导致较高的死亡数量,这也表明群落内有较强的自疏和它疏作用[18,20]。然而,常晓敏等研究认为对于不同径级的林木,各影响因子对树木死亡的影响权重不同,胸径<10 cm 时,林木相对断面积对树木死亡的影响最大[16]。灌木树种补员率较高,而乔木树种补员率较低,这也说明灌木树种耐荫性强,在林下能够更好地更新,这是森林下层恢复力高的表现[25]。
3.1.3 对森林蓄积量调查的参考 胸高断面积是计量群落生物量的主要依据之一[12],而平均胸径是衡量森林群落的一个重要指标,同时也是计算林分蓄积量、林木生物量的依据之一。在森林资源二类调查中,经常采用估算树种组成、平均胸径、林分平均树高、单位面积胸径≥5 cm 的树木株数的方法来测算林分蓄积量。以上测树因子估算尤其是平均胸径估算往往具有较大的难度,对于没有丰富经验的工作人员来说,只有通过在各种类型的森林中建立若干标准地,通过数据分析,来对标准地所在森林的测树因子加以确定,在此基础上对类似森林的测树因子进行估算。森林中除了活立木胸径在增大外,还有树木死亡和树木进径(胸径<5 cm 者长成胸径≥5 cm)影响着林分的平均胸径。死亡树木分布如果大径级比例较高,则会使平均胸径有减少的趋势,如本研究的201、202 和3 号样方;相反,死亡树木分布如果小径级比例较高,则使平均胸径有增加的趋势,如本研究的102 号样方。
一般除了人为干扰(如采伐)、自然灾害(如台风、病虫害发生)导致大树被清除(或者死亡),天然异龄林不会有较大比例的进径,因而进径对平均胸径的影响不会像死亡树木的影响那么大。标准地是判读林分测树因子的参考,选择标准地时不仅其树高、胸径要有代表性,树木死亡情况也要有代表性。对于空间尺度不大的样地来说,需要综合考虑树木死亡对林分平均胸径的影响,做标准地时需要对树木的死亡情况加以分析,这样才能更加准确地判读整个林分的测树因子,这在具体实践中要引起注意。当一个生态系统的活生物量与总生物量开始围绕一个平均值波动时,即达到稳态[12,26]。一个成熟稳定的常绿阔叶林,其生长相对较缓慢,在本研究中,5 年间,无论是对胸径≥1 cm 还是胸径≥5 cm 树木的统计,其平均胸径增长仅在(0.50±0.05)cm 左右,这也为森林资源年度监测数据更新提供了参考。
3.2 结论
2012—2017 年,在浙江五岭坑原生常绿阔叶林群落5 个样方中胸径≥1 cm 的树木死亡1 277 株,占2012年个体总数(10 459 株)的12.21%,隶属于93 个物种,占物种总数(142 种)的65.49%。死亡植株个体数量与径阶之间存在极显著的指数负相关关系。从功能型上看,灌木的死亡比例和补员比例明显高于乔木和小乔木,小乔木的补员比例高于乔木。主要树种的年死亡率、年补员率、种群变化率、适合度的平均值都表现出:灌木>小乔木>乔木。胸径≥1 cm 的活立木胸高断面积增长比例为4.59%,平均胸径增长0.55 cm;而胸径≥5 cm的活立木胸高断面积增长比例为1.04%,平均胸径增长0.47 cm。以上研究结果说明,该亚热带原生常绿阔叶林群落个体死亡率较高,补员率较低,林下更新困难,林木生长缓慢,种群动态存在较大的种间差异。