不同演替程度下马尾松人工林生物多样性对生物量的影响

韩泽民，李 源，王 熊，菅永峰，周靖靖,，佃袁勇,,3，黄光体

（1. 华中农业大学 园艺林学学院，湖北 武汉 430070；2. 华中农业大学 湖北林业信息工程技术研究中心，湖北武汉 430070；3. 华中农业大学 农业部华中都市农业重点实验室，湖北 武汉 430070；4. 湖北省林业调查规划院，湖北 武汉 430079）

生物量作为反映森林生态系统结构优劣和功能高度最直接的指标之一，其大小是植物群落与周围环境共同作用的结果[1]，分析森林生物量变化的驱动因素，可为林业可持续经营提供理论依据[2]。目前，以物种丰富度为代表的物种多样性是最简单的多样性衡量标准，通常被应用于生物量和生产力的预测指标[3-4]。有研究认为：生物多样性对生态系统功能的影响，不仅归因于物种数量，还依赖于物种所具有的功能特征[5]。具有较高功能多样性的生态系统常被认为拥有较高的生产力、恢复力和抗干扰能力，因此，通过比较群落间物种多样性与功能多样性的高低，可分析生态系统的好坏[6-7]。除此之外，林分密度也是控制林木生长与生物量积累的关键因素[8]，而且作为最常见的林分经营措施，调整林分密度，可以有效促进纯林向天然林过渡的进程[9]。马尾松Pinus massoniana具有速生、丰产等优良特性，一直被作为中国南方重要的用材树种和荒山绿化造林的先锋树种。然而，由于长期采取单一的纯林经营模式，立地质量降低，生产力下降，严重威胁了林地的可持续经营。目前，通过疏伐补植的措施将针叶纯林改造成针阔混交的近自然林，正逐渐成为替代大面积人工针叶纯林最有希望的途径之一[10]。周东雄[11]以采用疏林结构培育马尾松大径材兼营苦竹Pleioblastus amarus的复合林分为研究对象得出：混交经营有利于马尾松大径材的培育；邹文魁[12]调查分析了40年生马尾松-木荷Schima superba混交林的生物量结构，发现混交乔木层的总生物量高于纯林。综上所述，虽然混交林对马尾松林生物量的积累有一定促进作用，但是存在非自然、非随机的缺点，而且大部分观测对象的周期较短。与此同时，在研究影响生物量变化机制的问题上，一些学者通过研究林下植被生物多样性间接说明生物量的变化[13-14]，但忽视了乔木层生物多样性对其生物量变化的影响。作为生态系统中森林地上部分生物量的主体[15]，乔木层的生物量在群落中占有很大比例[16-17]，研究其生物量与生物多样性的关系具有重要意义[18]。针对混交林与纯林生物多样性以及林分密度之间的变化规律，如何探究出因地制宜的改造措施依然存在很多难点。鉴于此，本研究通过对比不同演替阶段的马尾松人工林生物多样性以及林分密度变化规律，为不同生长阶段的马尾松纯林改造提供有效的数据支持。

1 研究地区与研究方法

1.1 研究区概况

太子山林场位于湖北省荆门市京山县西南部 (30°48′35″～31°02′40″N，112°49′05″～113°03′40″E)，北倚大洪山，南接江汉平原。其地貌为低山丘陵区，属中亚热带季风气候，夏秋多雨，冬春干旱，年平均降水量1 094.8 mm，年平均气温16.4 ℃。太子山土壤呈酸性，境内土壤主要为黄褐土、黄棕壤、山地黄棕壤、黄褐色石灰土，其中分布面积最大的是黄棕壤和黄褐色石灰土，矿质养分、有机质状况总体良好。太子山林场内资源丰富，分布着地球上仅能在此成片生长的国家级保护植物对节白蜡Fraxinus hupehensis在内的常绿或落叶针叶林、阔叶林、林下植被等138科、204属近400种。既有秃杉Taiwania flousiana、鹅掌楸Liriodendron chinense、杜仲Eucommia ulmoides、厚朴Magnolia officinalis、楠木Phoebe zhennan、刺楸Kalopanax septemlobus等国家级保护植物，还有第三纪孑遗植物银杏Ginkgo biloba、水杉Metasequoia glyptostroboides等。其中，典型林分有栓皮栎Quercus variabili、麻栎Quercus acutissima、杉木Cunninghamia lanceolata、马尾松Pinus massoniana、柏木Cupressus funebris、湿地松Pinus elliottii、枫香Liquidambar formosana、樟树Cinnamomum camphora等。

1.2 样地设置与调查

在2018年7月和2019年6月，依据太子山林场马尾松人工林近自然经营方案，针对早期疏伐后的纯林，采用空间代替时间的方法，分别选择自然演替30年生、40年生以及50年生的马尾松林作为调查对象。其中30年生马尾松林共调查12块样地，这些样地还处于演替早期阶段，并且在原有树种基础上新增了1～2个树种；40年生马尾松林共调查18块样地，这一演替程度的马尾松林包含的演替树种已经达3～4个；50年生的马尾松林共调查18块样地，此阶段的马尾松林在自然演替下新增了5～7个树种。所有样地共计48块，每个样地设置20 m× 20 m样方，记录坡度、坡向、林龄，并对胸径≥5 cm的木本植物每木检尺。

1.3 生物量计算

生物量计算过程是在统计样地林分组成概况之后，对样地的所有树种进行划分，主要包括马尾松、麻栎、樟树、柏木、杉木、其他阔叶树种。每种树种的生物量计算公式见表1。

表1 生物量计算方法Table 1 Biomass calculation methods

1.4 林木生长方程

引用了亚热带气候带下马尾松与阔叶树种的林木材积生长率方程，并通过林分生物量扩展因子方程将每个样地在10 a的总材积生长量转换为生物量[23-25]。由于柏木与杉木都分布在40年生的样地中，不需要对其进行转换，所以本研究没有考虑其转换的材积生长率公式，只引入了马尾松和其他阔叶树种的生产率方程来计算每木材积生长率(表2)，最后通过生物量扩展因子方程将材积换算为生物量(表 3)。

表2 单木材积生长率方程Table 2 Wood growth rate equation of P. massoniana and broadleaved trees

表3 生物量扩展因子方程Table 3 Biomass expansion factor equation of P. massoniana and mixed coniferous forests

1.5 物种多样性计算

为了描述群落内物种多样性的分布特征，借助R语言中vegan包的diversity函数计算了Shannon-Wiener指数、Simpson指数以及Pielou指数，计算公式参考文献[26]。

1.6 功能多样性计算

功能多样性参数测量采用R的FD包计算每个样地的功能丰富度(FRic)、功能均匀度(FEve)、功能离散度(FDis)、功能多样性(RaoQ)，计算公式及含义参见文献[27]。

1.7 数据分析

为了排除生物量随着年龄增加而变化的影响，在分析生物量与各因素之间的关系之前，通过林木生长模型将各演替阶段的生物量统一换算为40年生。同时，为了避免立地条件对生物量的影响，选择在相同土壤条件下位于同一坡向以及相近坡度的林分进行研究。通过单因素方差分析探究了不同演替程度下各参数的差异(显著性水平为α=0.05)。使用普通最小二乘(OLS)多元回归模型来检验多样性、林分密度对森林生物量的影响。在拟合模型前所有数据都进行了标准化处理。各参数的计算以及相关性分析主要借助R 3.5.1来实现。

2 结果与分析

2.1 不同演替程度下生物量、生物多样性和林分密度的分布特征

生物量在30年生和40年生、50年生的马尾松林之间均具有显著差异(P＜0.05)，然而40年生与50年生马尾松林差异不显著(图1A)。在物种多样性指数中，Shannon-Wiener指数在3种不同演替程度下差异显著(P＜0.05)(图1B)；Simpson指数的变化与Shannon-Wiener指数一致，在不同演替程度下差异显著(P＜0.05)(图1C)；Pielou指数也随着演替程度的增加呈递增趋势，30年生与40年生、50年生的马尾松林之间均差异显著(P＜0.05)，40年生与50年生马尾松林差异不显著(P＞0.05)(图1D)。功能多样性指数中，功能丰富度在50年生马尾松林中最高，40年生马尾松林次之，30年生马尾松林最低，30年生与40年生、50年生的马尾松林之间均差异显著(P＜0.05)，40年生与50年生马尾松林差异不显著(P＞0.05)(图1E)；功能离散度在50年生马尾松林中明显高于其他，40年生与30年生、50年生的马尾松林之间均差异不显著(P＞0.05)，30年生与50年生马尾松林显著差异(P＜0.05)(图1F)；功能均匀度随着演替的进行不断提高，且在3种不同演替程度下具有显著差异(P＜0.05)(图1G)；功能多样性30年生与40年生、50年生马尾松林之间均差异显著(P＜0.05)，40年生与50年生马尾松林无显著差异(P＞0.05)(图1H)。林分密度随着演替的进行不断提高，40年生与30年生、50年生马尾松林之间均无显著差异(P＞0.05)，30年生与50年生马尾松林差异显著(P＜0.05)(图1I)。

图1 不同演替程度下生物量、生物多样性和林分密度的变化Figure 1 Change of biomass, biodiversity and stand density of different remolding modes

2.2 线性逐步回归建模筛选解释生物量的最佳子集

以生物量作为响应变量，多样性指数、林分密度为解释变量构建线性逐步回归模型，筛选出对生产力解释率最高的最佳子集。从表4可知：虽然多元线性模型R2=0.468并不高，但拟合结果差异极显著(P＜0.001)。对生物量影响的重要性从高到低依次为功能离散度(42.71%，P＜0.001)、Simpson指数(21.65%，P＜0.001)、Pielou指数(19.58%，P＜0.01)、Shannon-Wiener指数 (16.05%，P＜0.01)。

表4 逐步回归模型筛选结果Table 4 Stepwise regression model fitting results

3 讨论

3.1 马尾松人工林生物量与生物多样性的分布特征

本研究发现：40年生马尾松林生物量相较于30年生马尾松林有明显增加，与50年生马尾松林相比并没有较大差距；但Shannon-Wiener指数、Simpson指数以及Pielou指数随着演替的进行在不断提高。表明演替前期拥有不同性状的树种，可以提升资源的利用效率，这将对生物量的增加产生积极效应[28]。而演替后期，生态系统趋于稳定，树种的增加对林分生物量的影响逐渐减弱。本研究表明：不同演替程度之间功能多样性变化也较为显著。其中功能丰富度递增，说明群落所占据的功能空间有所提高[29]。而且从功能离散度和功能多样性的分布特征也不难看出，树种多的群落中不仅其生态位的互补程度较高，而且具有不同功能特性的树种相对增加，这为群落间个体生产力的提升提供了可能，并增加生态系统的功能。功能均匀度可以体现群落内树种对有效资源全方位的利用效率[29]。功能均匀度随着演替的进行不断提高，这有利于群落对资源利用效率的提高。随着演替的进行，林分密度开始增加，但在演替50年生的样地里也存在林分密度较小的情况，总体分布范围较30年生的有所增大，树种的增加使得单位面积上的树种数量增加。

3.2 马尾松人工林物种多样性与生物量之间的关系

在逐步回归模型筛选的最佳子集中，功能离散度对于生物量的影响占有重要权重。相较于物种多样性，功能多样性指标可以更有效地解释生态系统功能变化[30]。本研究显示：Shannon-Wiener指数、Simpson指数、Pielou指数都对生物量有显著影响，这与温带针叶林、亚热带人工林和热带森林的研究结果相似[31]，与TILMAN等[32]的物种多样性对生态系统生物量存在正效应相符，说明在一定研究尺度下，马尾松人工林物种多样性的适度提高会对生物量的增加有促进作用。就功能多样性与生物量之间的关系而言，功能丰富度与碳储量显著相关[33]，也有人只证明了功能离散度与碳储量的相关性[34]，本研究得出了相似的结论。综合分析模型中，林分密度对生物量没有表现出影响，这可能是生物多样性指标掩饰了林分密度的效应。

4 结论

马尾松人工林在演替进程中，随着树种的增加，物种多样性与功能多样性也随之增加，这使得群落能够占据更多的生态位，从而为获得更多的资源创造条件；同时由于树种间生长特性的不同，增加了群落利用资源的多样性，为个体生物量的积累提供了条件。回归模型结果表明：功能离散度可以更有效地解释生物量的变化，林分中的功能特性的离散化有利于减小群落内的竞争效应，为个体生物量的积累创造有利条件。