李非凡， 孙 冰， 裴男才， 闫玮明， 罗鑫华

（1. 中国林业科学研究院 热带林业研究所， 广东 广州510520； 2. 南京林业大学 林学院， 江苏 南京210037； 3. 广东乐昌杨东山十二度水省级自然保护区 管理处， 广东 乐昌512231）

生态化学计量学己成为连接分子、 细胞、 种群、 群落和生态系统等不同尺度生物学的新工具［1］。 生态化学计量学也在消费者驱动的养分循环、 限制性养分元素判别以及全球碳、 氮、 磷循环等方面发挥重要作用［2］。 氮、 磷是陆地生态系统中限制植物生长最普遍的2 种元素［3-4］， 对植物各种功能起着非常重要的作用［5］， 而碳氮比、 碳磷比、 氮磷比化学计量特征可以为生态系统中营养结构变化［6］、 生物多样性和生物地球化学循环研究提供基本依据。 森林凋落叶是土壤动物、 微生物的重要物质和能量来源， 也是土壤养分循环的重要参与者［7］。 有90%以上植被的氮、 磷养分通过凋落叶向土壤归还［8］。 根系不仅是植物固定和机械支持的器官， 也在水分和养分吸收、 同化物分配方面至关重要［9］。 中国人工林面积居世界之首， 因林分结构单一、 培育措施不健全导致生物多样性较低和林分生产力下降等问题日益严重［10］。 目前， 国内外已经开展了大量有关人工林植物和土壤的生态化学计量特征研究， 主要是探究不同森林类型、 区域、 演替阶段的植物叶片、 凋落叶、 土壤单独或者两两甚至是三者耦合的生态化学计量特征［11-14］。但是， 根系作为连接地上部分和地下部分重要的植物器官， 往往被忽略重要性， 大部分是凋落叶-土壤的研究［15-16］。 目前， 凋落叶-根系-土壤三者的耦合在国内已在杉木Cunninghamia lanceolata中被研究［17］。研究不同林分凋落叶-根-土生态化学计量特征， 有助于全面、 系统地了解森林生态系统养分循环规律，对当地的土壤养分和林分健康状况具有指导意义。 红锥Castanopsis hystrix是中国南亚热带重要的乡土树种， 速生树种， 也是优良的家具用材。 关于红锥的生态化学计量已有报道， 但多集中在枝叶［18］、 土壤［19］等。 乐昌含笑Michelia chapensis是木兰科Magnoliaceae 含笑属Michelia植物， 具有很高的经济价值， 是南方重要的阔叶树种， 种子含有丰富的油脂， 且人体中必须的亚油酸含量高达64.53%［20］。 已有学者在乐昌含笑不同植物器官的养分含量［21］、 土壤养分特征［22］等方面做了相关研究。 有关南亚热带次生林和人工林凋落叶-根系-土壤的碳、 氮、 磷含量和生态化学计量特征的研究较少。 本研究以粤北乐昌市龙山林场为研究区， 对红锥林、 乐昌含笑林和次生林为研究对象， 通过对不同林分养分质量分数及化学计量的定量比较， 分析人工林和次生林在养分分布方面的异同点及其诱因， 掌握养分现状和变化规律及其在不同部位的存在情况， 为该地区的林分经营提供理论依据。

1 研究区域与研究方法

1.1 研究区概况

研究区位于广东省乐昌市（24°57′～25°31′N， 112°51′～113°34′E）龙山林场， 该区属中亚热带季风气候， 主要为流水地貌和岩溶地貌。 年平均气温为19.9 ℃， 无霜期约300.0 d， 雨水丰沛， 年平均降水量为1 460.0 mm。 土壤类型为红壤土， 偏酸性。 在2018 年3 月选取红锥林、 乐昌含笑林和次生林3 种林分， 每种林分3 个重复， 海拔高度和坡面基本一致， 样地面积为20 m × 30 m。

红锥林林下植被灌木层主要有红锥、 杉木、 深山含笑Michelia maudiae； 草本层主要有石楠藤Piper wallichii、 阔片短肠蕨Allantodia matthewii、 扇叶铁线蕨Adiantum flabellulatum。 乐昌含笑林林下植被灌木层主要有石楠藤、 粗叶榕Ficus hirta、 杜茎山Maesa japonica、 常山Dichroa febrifuga； 草本层主要有华山姜Alpinia chinensis、 火炭母Polygonum chinense、 柳叶箬Isachne globosa。 次生林乔木层主要有润楠Machilus pingii、 木荷Schima superba、 酸枣Ziziphus jujubevar.spinosa、 新木姜子Neolitsea aurata、鹅掌柴Schefflera octophylla； 林下植被灌木层主要有石楠藤、 爬藤榕Ficus sarmentosavar.impressa、 草珊瑚Sarcandra glabra、 小叶女贞Ligustrum quihoui； 草本层主要有石楠藤、 华山姜、 草珊瑚、 小叶女贞。

样地概况和土壤理化性质见表1 和表2。

表1 不同林分乔木层基本情况Table 1 Basic conditions of arbor layers in different forest stands

表2 不同林分土壤（0～20 cm）的理化性质Table 2 Soil physical and chemical properties of different forest stands

1.2 样品采集与处理

在每个样地内按照 “梅花” 五点法设置凋落物收集框， 收集框规格由孔径1 mm 的尼龙网布制成，上端1 m × 1 m， 下端距离地面0.25 m， 用聚氯乙烯（PVC）管固定在泥土中的容器。 根系和土壤样品均采用6 cm 根钻及“梅花” 五点取样法采集： 根系用钻土芯法在样地内采集0～20 土层的所有根系， 土壤的采样也采用“梅花” 五点取样法用直径6 cm 的根钻采集0～20 cm 土层的土壤， 编号保存带回实验室。凋落叶放入烘箱中70 ℃烘至恒量， 研磨过40 目筛， 用于全碳、 全氮、 全磷养分分析。 带有根系的土壤用40 目的网筛在流水下冲洗， 收集所有根系， 未分级。 土壤一部分用于测量物理性质， 另一部分过60目筛， 用于测定全碳、 全氮、 全磷养分。

1.3 样品测定

依据《土壤农化分析》［23］对植物和土壤进行元素全碳、 全氮、 全磷质量分数的检测。 植物、 土壤有机碳采用重铬酸钾-外加热法测定； 植物全氮、 全磷采用浓硫酸-过氧化氢（H2SO4-H2O2）消煮法， 消解后的待测液供全氮、 全磷的测定， 其中全氮采用蒸馏滴定法测定， 全磷采用钒钼黄比色法测定； 土壤全氮、全磷分别采用开氏-蒸馏滴定法和氢氧化钠融熔-钼锑抗比色法。 数据为平均值±标准误。

1.4 数据处理

用Excel 2007 进行数据记录与整理， SPSS 19.0 进行数据分析。 对3 种林分优势树种进行正态检验，不符合正态分布的数据进行对数转换， 使之符合正态分布， 再检验方差齐性。 对样品全碳、 全氮、 全磷以及碳氮比、 碳磷比、 氮磷比采用单因素方差分析（one-way ANOVA）的最小显著（LSD）法进行差异性检验， 利用Pearson 相关性分析和回归分析对样品碳氮比、 碳磷比、 氮磷比进行相关性分析和线性回归， 用Origin 9.1 作图。

2 结果与分析

2.1 不同林分类型凋落叶-根系-土壤的全碳、 全氮、 全磷质量分数特征

由表3 可知： 红锥林全碳、 全氮质量分数在凋落叶、 根系和土壤中差异均显著（P＜0.05）， 而全磷质量分数无显著性差异（P＞0.05）； 乐昌含笑林全碳质量分数在凋落叶、 根系和土壤中均存在显著差异（P＜0.05）； 次生林全碳、 全氮、 全磷质量分数在凋落叶、 根系和土壤中均为差异显著（P＜0.05）。 凋落叶中的全碳和全磷质量分数在3 种林分间均无显著性差异（P＞0.05）； 3 种林分的凋落叶、 根系和土壤的全碳质量分数均为根系＞凋落叶＞土壤； 而全氮质量分数则均为凋落叶＞根系＞土壤； 红锥林和次生林的全磷质量分数均为凋落叶＞根系＞土壤， 而乐昌含笑的全磷质量分数则为根系＞凋落叶＞土壤。 对于全碳质量分数， 3 种凋落叶的全碳为乐昌含笑林＞次生林＞红锥林， 根系和土壤中的全碳则分别表现为次生林＞乐昌含笑＞红锥林和红锥林＞次生林＞乐昌含笑； 对于全氮质量分数， 乐昌含笑林的凋落叶和根系无显著差异（P＞0.05）， 而红锥林和次生林差异显著（P＜0.05）； 凋落叶和土壤中的全氮质量分数均为红锥林＞次生林＞乐昌含笑林， 根系中全氮质量分数大小为乐昌含笑林＞次生林＞红锥林人工林； 对于全磷质量分数， 凋落叶中的全磷为次生林＞红锥林＞乐昌含笑林， 根系中的全磷为乐昌含笑林＞次生林＞红锥林， 土壤中的全磷为红锥林＞乐昌含笑和次生林。

表3 3 种林分凋落叶-根系-土壤的全碳、 全氮、 全磷质量分数Table 3 C, N and P contents of litter-root-soil of three forest stands

2.2 不同林分类型凋落叶-根系-土壤的全碳、 全氮、 全磷化学计量比特征

表4 3 种林分凋落叶-根系-土壤的碳氮比、 氮磷比、 碳磷比Table 4 C/N, N/P, C/P stoichiometry of three forest stands of litter-root-soil

由表4 可知： 红锥林的碳氮比和氮磷比在凋落叶、 根系和土壤中均差异显著（P＜0.05）， 碳磷比在红锥林凋落叶与根系无显著差异（P＞0.05）， 而在土壤与根系和凋落物之间显著性差异（P＜0.05）； 乐昌含笑林的碳氮比和碳磷比在根系和凋落叶中差异不显著（P＞0.05）， 而在土壤与根系和凋落物中差异显著（P＜0.05）， 氮磷比在凋落叶、 根系和土壤中均差异显著（P＜0.05）； 次生林的碳氮比、 碳磷比和氮磷比在凋落叶、 根系和土壤中均差异显著（P＜0.05）。 碳氮比在3 种林分中均为根系＞凋落叶＞土壤； 碳磷比在红锥林和次生林中均为根系＞凋落叶＞土壤， 而在乐昌含笑林中为凋落叶＞根系＞土壤； 氮磷比在红锥林、 乐昌含笑林和次生林中均为凋落叶＞根系＞土壤。 土壤碳氮比在红锥林、 乐昌含笑林和次生林中均值分别为24.48、 24.60 和23.41， 且在3 种林分间均为差异不显著（P＞0.05）； 凋落叶、 土壤和根系的碳磷比在3 种林分间均为差异不显著（P＞0.05）； 土壤氮磷比在红锥林、 乐昌含笑林和次生林均值分别为4.19、 4.95 和5.39， 且在3 种林分间差异不显著（P＞0.05）， 而凋落叶和根系氮磷比在3 种林分间均为差异显著（P＜0.05）。

图1 3 种林分凋落叶-根系-土壤全氮和全磷质量分数的相关关系Figure 1 Correlations between N and P concentrations in litter-root-soil of three forest stands

图2 3 种林分凋落叶-根系-土壤氮磷比和碳氮比（对数）的相关关系Figure 2 Correlations between log-transformed N/P and C/N in litter, roots and soil of three forest stands

图3 3 种林分凋落叶-根系-土壤氮磷比与碳磷比（对数）的相关关系Figure 3 Correlations between log-transformed N/P and C/P in litters, roots and soil of three forest stands

2.3 凋落叶-根系-土壤全碳、 全氮、 全磷与化学计量比的相关性分析

由图1A 可知： 研究区域红锥林根系全氮和全磷质量分数呈显著正相关关系（R2=0.507 9，P＜0.05），凋落叶与土壤的全氮和全磷相关不显著（R2=0.135 4，P＞0.05；R2=0.274 7，P＞0.05）。 乐昌含笑林凋落叶和根系的全氮和全磷均极显著正相关（R2=0.716 9，P＜0.01；R2=0.765 1，P＜0.01）（图1B）。 次生林凋落叶和根系中的全氮和全磷也呈极显著正相关（R2=0.821 6，P＜0.01；R2=0.771 9，P＜0.01）（图1C）。 3种林分的土壤全氮和全磷质量分数相关性均不显著（P＞0.05）。

经过对数转换， 红锥林土壤碳氮比与氮磷比呈极显著正相关（R2=0.721 2，P＜0.01）， 凋落叶和根系的碳氮比与氮磷比相关关系不显著（图2A）； 凋落叶、 根系和土壤三者的碳磷比与氮磷比均呈极显著正相关（R2=0.747 9，P＜0.01；R2=0.893 6，P＜0.01；R2=0.663 5，P＜0.01）（图3A）。 乐昌含笑林和次生林凋落叶、 根系和土壤的碳氮比与氮磷比相关关系均不明显（图2B 和图2C）， 而碳磷比与氮磷比均呈显著正相关（R2=0.602 5，P＜0.05；R2=0.527 6，P＜0.05；R2=0.911 6，P＜0.01）（图3B）。 次生林凋落叶和土壤碳磷比与氮磷比呈显著正相关（R2=0.676 4，P＜0.05；R2=0.824 3，P＜0.01）（图3C）。

3 讨论

3.1 凋落叶-根系-土壤全碳、 全氮、 全磷质量分数特征

凋落叶养分作为土壤外来养分的主要来源， 在植物生长过程中发挥着重要作用。 凋落叶通过物理作用淋溶阶段和化学作用地下微生物的分解将有机物分解成无极小分子供给植物吸收。 凋落叶养分与分解速率及初始养分含量有关， 甚至通过新鲜叶片养分含量的转移影响凋落叶的养分含量， 也可能与叶片形态和解剖结构有关。 比如叶片厚度、 气孔密度、 叶肉细胞栅栏组织和海绵组织的层数及厚度等都会影响养分的分布情况。 红锥林、 乐昌含笑林和次生林3 种林分凋落叶全氮、 全磷质量分数均低于在全球尺度上研究的10.93 和0.85 mg·g-1［24］， 但是全磷质量分数明显低于全球平均全磷。 这可能与不同树种的生理生态特征以及当地的小气候等有关。 凋落叶的全碳、 全磷在3 种林分间无显著差异， 但红锥林凋落叶全氮质量分数明显高于含笑林和次生林， 凋落叶作为植物和土壤养分之间的基本载体， 说明红锥林在养分归还能力方面较含笑林和次生林强。

根系是植物吸收氮、 磷等养分的重要器官［25］。 VOGT 等［26］发现： 通过细根周转归还到土壤中的全氮、 全磷超过了地上凋落叶。 红锥林、 含笑人工林和次生林的根系全氮和全磷质量分数平均值低于中国陆地植物根系平均水平（9.2 和1.0 mg·g-1）。 乐昌含笑林根系的全碳、 全氮、 全磷质量分数均高于红锥林， 这可能与根系的形态特征（表面积、 体积等）有很大的关系。 本研究中红锥林根系全碳质量分数为（465.10±2.81） mg·g-1， 而邵梅香［19］研究的红锥细根（直径＜0.5 mm）全碳质量分数平均值为（514.25±9.65）mg·g-1， 这与两者研究选择的根系直径范围不同以及不同的环境基底有关。 红锥林根系的全碳、 全氮、全磷质量分数低于其他2 种林分， 而土壤的全碳、 全氮、 全磷质量分数高于其他2 种林分， 与红锥林根系养分状况相比， 其他2 种林分呈相反关系。 这可能与红锥是菌根营养性树种有关［27］， 可以利用菌根菌获取养分， 相应的根系获得的养分较少。 郭大立等［28］对土壤氮的有效性与细根生产量和周转率研究表明： 随着土壤全氮有效性的增加， 细根生产量和周转率都提高或下降， 具体提高还是下降， 可以将其与养分含量相联系。 本研究可能会从侧面论证随着土壤全氮有效性的增加， 细根生产量和周转率都下降这一观点的可靠性， 因此可以初步确定乐昌含笑林的根系吸收并储存养分的能力比红锥林和次生林要好。

土壤作为植物养分直接来源库， 养分含量的多少关系到植物生长状况、 森林健康情况以及景观生态评估质量的好坏。 红锥林全碳、 全氮、 全磷质量分数均值分别为38.40、 1.56 和0.38 mg·g-1； 乐昌含笑林土壤全碳、 全氮、 全磷质量分数均值分别为29.74、 1.21 和0.27 mg·g-1； 次生林土壤全碳、 全氮、 全磷质量分数均值分别为33.89、 1.44 和0.27 mg·g-1， 与杉木林［29］（土壤碳为22.57 mg·g-1， 土壤氮为1.33 mg·g-1， 土壤磷为0.46 mg·g-1）相比， 呈高全碳低全磷的格局， 说明红锥林、 乐昌含笑林和次生林对土壤碳库的维持较强， 而对土壤磷的固持较弱。 土壤全碳、 全氮、 全磷质量分数均为红锥林最高， 说明红锥林的养分维持和养分归还能力较含笑林和次生林强。

3.2 凋落叶-根系-土壤生态化学计量比特征

凋落叶的氮磷比影响其分解和养分归还速率， 且氮磷比与分解速率呈负相关关系， 氮磷比大于25时分解主要受全磷限制［30］。 红锥林凋落叶的氮磷比（27.93）大于25， 表明红锥凋落叶在分解过程中主要受磷元素的限制， 而乐昌含笑林和次生林凋落叶的氮磷比小于25， 说明乐昌含笑林和次生林的凋落叶分解时主要受氮限制。 该地区的土壤碳、 氮、 磷范围分别为29.74～38.40、 1.21～1.56 和0.27～0.38 mg·g-1，与全国土壤碳（11.12 mg·g-1）、 氮（1.06 mg·g-1）、 磷（0.65 mg·g-1）相比［31］， 这3 种林分的土壤表现出高氮低磷的现状， 与热带和亚热带地区养分状况相似。 “生长速率假说” 认为： 植物在快速生长过程中需要大量的磷使核糖体能够快速合成大量蛋白质， 因此会表现出较高的氮磷比［32］。 本研究的红锥林、 乐昌含笑林和次生林的根系氮磷比分别16.37、 17.66、 20.17， 高于中国的平均水平14.27。 说明生长速率为次生林＞乐昌含笑林＞红锥林， 而前期是3 种林分的土壤养分大致相同， 造成根系氮磷比不同的主要原因是林分树种的不同， 树种不同的本质间接就是根系形态特征不同， 从而影响对土壤养分吸收。 3 种林分的碳氮比和碳磷比均为根系＞凋落叶＞土壤， 而氮磷比均为凋落叶＞根系＞土壤。 凋落叶中的碳氮比和碳磷比以乐昌含笑林最高， 这可能与乐昌含笑植物叶面积较大， 呈革质叶较厚有关。 土壤的氮磷比也为次生林＞乐昌含笑林＞红锥林， 说明根系和土壤的氮磷比表现相同的变化趋势。 土壤中的碳氮比、 碳磷比和氮磷比在3 种林分中均无显著性差异（P＞0.05）， CHEN 等［33］对同一林龄不同人工林的土壤研究表明： 土壤全碳含量没有明显不同。 而土壤中的营养元素主要来源于岩石风化， 对于同一区域的土壤养分变化不会表现出过大的差异。 碳磷比在3 种林分的凋落叶、 根系和土壤中也无显著性差异。

3.3 凋落叶-根系-土壤全碳、 全氮、 全磷质量分数与化学计量比的相关性

全氮和全磷是各种蛋白质和遗传物质的重要组成元素， 植物体内的全氮和全磷往往随着植物体的生长发生变化［34］。 一般认为植物在代谢过程中， 全碳的固定需要大量蛋白酶和核酸的参与， 在全氮和全磷共同参与植物生理代谢过程中， 往往表现出一致性［35］。 植物养分的主要来源是凋落物的归还， 而凋落叶又是凋落物的主要组成部分， 碳、 氮、 磷养分循环在“凋落叶-根系-土壤” 间相互转换， 其中养分元素含量和计量特征以及相关性都可以反映森林生态系统养分状况、 养分循环和利用情况［36］。 3 种林分的根系全氮和全磷质量分数均显著正相关， 凋落叶和根系的全氮质量分数和碳氮比均显著相关； 3 种林分凋落叶和根系全磷质量分数和碳氮比、 碳磷比均为显著负相关， 红锥林和乐昌含笑林的凋落叶和根系全磷质量分数也与氮磷比呈显著负相关。 红锥林和乐昌含笑林的凋落叶、 根系和土壤的碳磷比和氮磷比均为显著正相关； 次生林的凋落叶和土壤中的碳磷比和氮磷比显著正相关， 而根系中的碳磷比和氮磷比无相关性， 这可能是次生林根系取样为混合样品， 涉及的树种种类较多， 没有一定的规律性。 但是凋落叶和根系的全氮和全磷质量分数在一定程度上表现出相似性， 说明凋落叶养分和根系养分有一定的关联。 大多数陆地生态系统植物吸收的养分， 有90%以上是凋落叶向土壤的归还， 而根系又是植物吸收养分的重要器官， 因此凋落叶养分和根系养分表现出类似的规律可以得到合理的解释［7］。 相比较次生林， 本研究的2 种人工林在生态化学计量相关性方面比次生林更强烈直观， 最重要的原因是次生林涉及到的树种变量更多， 人工林树种单一， 研究对象明确， 更能定量准确地反映化学计量特征相关性。

4 结论

同一地区不同林分的全碳、 全氮、 全磷生态化学计量特征主要与所用的树种有关。 本研究土壤碳氮比、 碳磷比、 氮磷比在3 种林分间无显著差异， 说明同一区域的碳∶氮∶磷维持在一定的范围。 通过比较凋落叶、 根系和土壤的养分情况可知： 研究地全氮偏高、 全磷明显偏低， 营林过程中可以适当地施用磷肥改善土壤养分。 红锥林土壤养分和凋落叶养分回归状况优于乐昌含笑林和次生林， 研究地可以选择红锥林作为合适的人工林树种。