川西亚高山天然次生林不同演替阶段土壤团聚体组成及有机碳分布特征

蔡蕾, 候笔锋, 冯秋红, 刘前程, 李旭华, 潘红丽, 刘兴良*

1.四川省林业科学研究院，森林和湿地生态恢复与保育四川重点实验室，四川 成都 610081；

2.四川卧龙森林生态系统国家定位观测研究站，四川 成都 610081；

3.阿坝州松潘国有林保护局，四川 阿坝 623300；

4.四川农业大学，四川 成都 611130

土壤是全球陆地生态系统中最大的碳库[1]，气候变化和生态系统中植被类型的不同会引起土壤呼吸过程的变化，因此演替进程是影响森林生态系统碳库动态的一个极为重要的影响因素[2]。随着演替阶段的变化，森林生态系统的群落结构、物种组成以及物种丰富度都会产生巨大的变化，生物量的积累以及碳分配也会发生变化，进而影响森林生态系统土壤碳库变化和碳吸存潜力[3-4]。

土壤团聚体作为土壤结构的基本单元，是反映土壤肥力的一项重要指标[5]。不同粒级的团聚体在营养元素的保持、供应及转化能力等方面发挥着不同的作用。对不同演替阶段次生林土壤团聚体的稳定性进行研究可以反映出土壤生态特征以及土壤有机碳的储存能力，揭示土壤特征与演替进程的关系。有机碳在团聚体中的分布形式决定了土壤储存和保留有机碳的能力[6]，有机碳水平的提高能促进土壤结构的形成并增强其稳定性。有研究表明，不同粒级团聚体的有机碳含量存在明显差异，一部分认为团聚体有机碳含量随粒级的增加而降低[7-8]，另一部分研究则提出了与之相悖的观点，指出大部分土壤有机碳储存在大粒径团聚体中，＜0.25 mm 粒级团聚体有机碳含量相对较低[9-10]。

川西亚高山森林是我国西南亚高山林区水源涵养林的重要组成部分，以冷杉为主的原始暗针叶林在经历大规模采伐利用后，天然次生林已成为该区域的主要森林类型[11-12]，因此，研究不同演替阶段的次生林类型的土壤结构及有机碳分布特征，对于准确评估本区森林生态系统的碳收支及大尺度的碳循环模型构建具有非常重要的意义。

1 材料与方法

1.1 研究区概况

研究区位于四川省理县米亚罗林区，地理坐标31º24′—31º55′N，102º35′—103º4′E。该区位于青藏高原东缘褶皱带最外缘部分，具有典型的高山峡谷地貌。气候受太平洋、印度洋及青藏高原三大气团影响，年均温 6 ℃～12 ℃，极端最高气温 32 ℃，极端最低气温−16 ℃。夏季温凉多雨，冬季寒冷干燥，最冷月（1 月）平均气温为−8 ℃，最热月（7 月）均温为 12.6 ℃，≥10 ℃年积温为 1 200 ℃～1 400 ℃，年蒸发量为 1 000～1 900 mm，常年无霜期 200 d，年均降水量为 600～1 100 mm，具有降水次数多、强度小的特点。米亚罗林区植被垂直成带明显，其类型和生境随海拔及坡向而分异。原生森林分布于海拔2 400～4 200m 之间，以亚高山暗针叶林为主，主要优势树种为岷江冷杉（Abies faxoniana），桦木（Betulasp.）、云杉（Picea asperata）、野樱（Pruns tatsienensis）等乔木有分布；灌木主要有臭樱（Maddenia hypoleuca）、红刺悬钩子（Rubus aurantiacus）、红脉忍冬（Lonicera nervosa）、绣线菊（Spiraea salicifolia）等；草本植物主要有野草莓（Fragaria vesca）、婆婆纳（Veronica didyma）、柳叶菜（Epelobium hirsutum）、六叶葎（Galium asperuloides）、橐吾（Ligularia sibirica）等。

该区成土母岩主要为千枚岩、板岩、白云岩等的残坡积风化物，极易风化，主要土壤类型为山地棕色森林土，土壤养分情况见表1。

表1 不同演替阶段次生林0～20 cm 土壤养分情况Tab.1 0−20 cm soil nutrient content of secondary forests in different successional stages

1.2 土壤样品采集

在研究区内依据林分组成，选取桦木纯林、桦木林+冷杉混交林、冷杉原始林作为次生林不同演替阶段的研究对象，分别代表初级演替、中级演替、高级演替三个演替阶段。每个演替阶段设置3～5 个20 m×20 m 的样方进行植物群落调查并在其中设1 m×1 m 的小样方3～5 个，采用5 点取样法对0～20 cm的表土进行取样，土壤样品用硬质塑料盒盛装，带回实验室，风干、剔除石块和根系后用于分析土壤团聚体结果，并测定土壤碳含量。

1.3 土壤团聚体分类及有机碳的测定

有机碳采用重铬酸钾氧化-外加热法测定。

土壤团聚体分类采用干筛法分离出＞5 mm、5～2 mm、2～1 mm、1～0.5 mm、0.5～0.25 mm、0.25～0.053 mm、＜0.053 mm 共7 级土壤团聚体，并计算各粒径土壤团聚体对土壤碳的贡献率（%）[9]。计算公式如下：

1.4 数据分析

所用测定结果用采Microsfot Excel 2012 进行整理和初步分析，用SPSS 16.0 进行方差分析，多重比较分析采用LSD 检验。

2 结果与分析

2.1 不同演替阶段土壤团聚体的分布特征

3 个演替阶段次生林粒径含量与粒径大小基本上呈正相关趋势，含量随粒径减小而降低，且均表现为＞5 mm 粒径团聚体分布均占主导地位。其中 ＞5 mm和＜0.053 mm 粒径的土壤团聚体在中级演替阶段（桦木+冷杉混交林）分布含量最高，2～1 mm 粒径的团聚体在高级演替阶段（冷杉原始林）分布大于其他两个演替阶段，而1～0.5 mm、0.5～0.25 mm、0.25～0.053 mm 3 个粒径的团聚体含量均是初级演替（桦木林）＞高级演替（冷杉原始林）＞中级演替（桦木+冷杉混交林）（见表2）。单因素方差分析发现，同一演替阶段内各径级团聚体含量分布存在显著差异；而同一粒径不同演替阶段间除了5～2 mm 粒径外，其余6 个粒径均存在显著差异（见表3）。

表2 不同演替阶段次生林土壤团聚体分布及组成Tab.2 Distribution and composition of soil aggregates in secondary forests at different successional stages

表3 同一粒径不同演替阶段间土壤团聚体含量单因素方差分析Tab.3 Single factor variance analysis of soil aggregates content at different successional stages with the same particle size

对同一演替阶段不同粒径团聚体含量进行多重比较发现初级演替阶段（桦木林）的团聚体在＞1 mm 粒 径 含 量 间 不 存 显 著 差 异，1～0.5 mm 与0.5～0.053 mm、0.25～0.053 mm 与＜0.053 mm 粒径存在显著差异；中级演替阶段（桦木+冷杉混交林）、高级演替阶段（冷杉原始林）＞0.25 mm 粒径大团聚体与 ＜0.25 mm 粒径微团聚体间均存在显著差异，其中冷杉林的大团聚体（＞0.25 mm）间各粒径均有显著差异，混交林在＞5 mm 与5～1 mm、2～1 mm 与1～0.5 mm、1～0.5 mm 与0.5～0.25 mm 间存在显著差异，而微团聚体（＜0.25 mm）在高级演替阶段（冷杉原始林）和中级演替阶段（桦木+冷杉混交林）均无显著差异。

虽然单因素方差分析结果显示除了5～2 mm 粒径外，其余6 个粒径均存在显著差异，但对同一粒径不同演替阶段团聚体含量进行多重比较发现，同一粒径内并不是3 个演替阶段均存在显著差异，如＞5 mm 粒径是桦木+冷杉混交林与桦木林、冷杉原始林存在显著差异，而桦木林与冷杉原始林间无显著差异。

由图1 可见，各演替阶段次生林均以土壤大团聚体（＞0.25 mm）为主，大团聚体比例皆达90%左右，微团聚体（＜0.25 mm）分布较少。各演替阶段次生林土壤大团聚体含量（＞0.25 mm）表现为冷杉+桦木混交林最高，微团聚体含量（＜0.25 mm）表现为冷杉+桦木混交林＞冷杉原始林＞桦木林，且桦木林与混交林、冷杉原始林间存在显著差异，这说明土壤团聚性呈“∩型”随着演替阶段上升在中级演替阶段达到最高，随着演替进程的发展到高级阶段又下降。

图1 不同演替阶段次生林团聚体分布情况Fig.1 Soil aggregates distribution of secondary

2.2 土壤团聚体的有机碳（含量）分布

3 个演替阶段土壤有机碳总含量介于159.24～561.28 g·kg−1之间，对不同演替阶段有机碳总含量进行比较，冷杉原始林＞桦木+冷杉混交林＞桦木林，但差异不显著。不同演替阶段各级团聚体的土壤有机碳含量均是冷杉原始林最高，桦木+冷杉混交林次之，桦木林最低，且差异性显著。土壤大团聚体（＞0.25 mm）有机碳的含量百分比分别为70.32%、70.21%、68.72%，微团聚体（＜0.25 mm）有机碳的含量百分比分别为29.68%、29.79%、31.28%，其中0.053—0.25 mm 粒径有机碳含量为14.61%、15.15%、15.05%，＜0.053 mm 有机碳的含量百分比为15.37%、15.58%、14.60%，微团聚体的两个径级间有机碳含量比例相当，差异不显著（见图2）。

对同一演替阶段土壤有机碳含量与各团聚体粒径的相关性进行分析，发现其二者相关性不显著，但有机碳的含量随团聚体粒径减小而基本呈增加趋势（见表4），3 个演替阶段均是所有＜0.25 mm 粒径的有机碳含量均值高于其他粒径，所有＞2 mm 粒径团聚体有机碳含量均值最低，这与刘敏英等[18]的研究结果一致。

图2 不同演替阶段次生林土壤团聚体有机碳百分比含量Fig.2 Percentage of organic carbon in soil aggregates forests at different successional stages of secondary forests at different succession stages

表4 不同演替阶段次生林土壤团聚体有机碳含量Tab.4 Percentage of organic carbon in soil aggregates of secondary forests at different successional stages

2.3 土壤团聚体对土壤有机碳的贡献率

如表5 所示，3 个演替阶段次生林类型的土壤团聚体对有机碳的贡献率均是＞5 mm 粒径最高，＜0.25 mm 粒径最低。土壤团聚体对有机碳的贡献率主要集中在大团聚（＞0.25 mm），微团聚体（＜0.25 mm）占的比重较少。单因素方差分析发现，同一演替阶段内各级团聚体对有机碳贡献率存在显著差异；而同一粒径不同演替阶段间分析结果显示，＞5 mm、2～1 mm、1～0.5 mm、0.5～0.25 mm、0.25～0.053 mm这5 个粒径对有机碳的贡献率在不同演替阶段间存在显著差异，5～2 mm、＜0.053 mm 这2 个粒径对土壤的贡献率在不同演替阶段间无显著差异（表6）。其中 ＞5 mm 粒径的土壤团聚体在中级演替阶段（桦木+冷杉混交林）分布含量最高，2～1 mm 粒径的团聚体在高级演替阶段（冷杉原始林）分布大于其他两个演替阶段，而1～0.5 mm、0.5～0.25 mm、0.25～0.053 mm 3 个粒径的团聚体含量均是初级演替（桦木林）＞高级演替（冷杉原始林）＞中级演替（桦木+冷杉混交林）（见表6）。

对同一演替阶段不同粒径团聚体的有机碳贡献率进行多重比较发现，初级演替阶段（桦木林）在＞0.5 mm 粒径时，团聚体的大小与有机碳贡献率是无显著差异的，而在0.5～0.053 mm 与＞0.5 mm 粒径间、0.25～0.053 mm 与＜0.053 mml 粒径间团聚体的机碳贡献率存在显著差异；中级演替阶段（桦木+冷杉混交林）除了5～2 mm 与2～1 mm 粒径间有机碳贡献率无显著差异，其余径级均具有显著差异；高级演替阶段（冷杉原始林）各粒径间团聚体对有机碳的贡献率均有显著差异。这说明随着演替阶段的上升，粒径的大小与土壤有机碳贡献率的联系更加紧密，并基本表现为随粒径增加有机碳贡献率增加。

表5 不同演替阶段次生林土壤团聚体对有机碳的贡献率Tab.5 Contribution rate of soil aggregates to organic carbon in secondary forests at different successional stages

表6 同一粒径不同演替阶段间土壤团聚体对有机碳贡献率单因素方差分析Tab.6 Single factor variance analysis of contribution rate of soil aggregates to organic carbon at different successional stages with the same particle size

对同一粒径不同演替阶段次生林土壤团聚体对有机碳的贡献率进行多重比较发现，＞5 mm 粒径时3 个演替阶段的有机碳贡献率分别存在显著差异，且中级演替（桦木+冷杉混交林）＞高级演替（冷杉原始林）＞初级演替（桦木林）；而在2～1 mm、1～0.5 mm、0.5～0.25 mm、0.25～0.053 mm、＜0.053 mm这5 个粒径间仅表现为桦木林与混交林、桦木林与冷杉林之间具有显著差异。说明演替阶段不同对＞5 mm 粒径的团聚体有机碳贡献率具有显著影响。

比较不同团聚体有机碳含量与贡献率发现，虽然＜0.25 mm 粒径的有机碳含量高但是贡献率低，说明微团聚体对有机碳的保护力更强，且桦木+冷杉混交林＞桦木林＞冷杉原始林（见表4 和表5）。

3 讨论与结论

3.1 不同演替阶段土壤团聚体的分布特征

有研究认为[13]，微团聚体的形成是大团聚体形成的前提条件，空间尺度上土壤团聚体由微团聚体向大团聚体逐级连续层次性形成，而在时间尺度上胶结物质从多糖（暂时稳定）向菌丝根系（短时间稳定）及芳香类物质（持久稳定）层次性变化。李玮[14]等的研究表明，种植年限的增加有利于团聚体胶结物质的积累，从而有利于大团聚体的形成与稳定；Ma[15]等也认为，中龄林和成熟林的土壤团聚性要强于幼龄林。本研究显示，处于中级演替阶段的冷杉+桦木混交林大团聚体含量最高，土壤团聚性呈“∩型”随着演替阶段上升在中级演替阶段达到最高，随着演替进程的发展到高级阶段又下降，这与李玮[14]、Ma[15]、Tisdall[13]的研究结论不一致。这可能是由于大团聚体的稳定性在很大程度上取决于植物根系和菌丝，当根系和菌丝被分解而得不到补给时，稳定性大团聚体的数量随有机质含量的上升而减少，本研究桦木冷杉混交林比冷杉原始林的有机质含量低，可能就是导致混交林大团聚体含量高于冷杉原始林的重要原因。

3.2 土壤团聚体的有机碳（含量）分布

3 个演替阶段土壤有机碳总含量介于159.24～561.28 g·kg−1之间，冷杉原始林＞桦木+冷杉混交林＞桦木林。不同演替阶段各级团聚体的土壤有机碳含量均是冷杉原始林最高，桦木+冷杉混交林次之，桦木林最低，且土壤大团聚体（＞0.25 mm）有机碳的含量百分比远高于微团聚体（＜0.25 mm）有机碳的含量百分比。说明各演替阶段土壤有机碳的含量主要是存在于大团聚体中，这与前人的研究结论[16-18]相同。对照土壤养分中铵态氮的含量来看，氮含量和有机碳变化规律趋于一致，这与黄天颖[2]等人对黄浦江上游水源涵养林碳、氮分布特征研究结论相同，原因可能在于土壤碳、氮变化通常相辅相成，团聚体氮元素含量随有机碳含量的变化而变化。

冷杉原始林的大团聚体分布较其他两个演替阶段稍低，但冷杉原始林有机碳总含量最高，出现这一结果可能是随着演替阶段的前进，干扰因子减少，植物及微生物群落结构趋于稳定，大团聚的更新降低，使得大团体内微团聚体的形成数量增大，而这种从大团聚体中心形成的微团聚体含有很高的且不易分解的有机碳含量，因此出现大团聚体含量低但有机碳总含量高的情况[19-20]。

对同一演替阶段土壤各粒径团聚体中有机碳含量分布的研究表明，有机碳含量与各团聚体粒径的相关性不显著，但有机碳的含量随团聚体粒径减小而基本呈增加趋势。3 个演替阶段均是＜0.25 mm 粒径的有机碳含量高于其他粒径，＞2 mm 粒径团聚体有机碳含量最低。本研究，＜0.25 mm 微团聚体有机碳含量最高，究其原因可能是小团聚体的固碳方式决定的。小团聚体主要靠有机和无机胶体紧密结合固持碳，固持的碳不易为微生物分解释放[21]，而大团聚体主要靠有机物、多糖等结合，因此大团聚体中的有机碳比微团聚体中的有机碳年轻，有机碳更容易矿化，而微团聚体中的有机碳则大多是高度腐殖化的惰性组[22-23]。

3.3 土壤团聚体对土壤有机碳的贡献率

3 个演替阶段次生林类型的土壤团聚体对有机碳的贡献率均是＞5 mm 粒径最高，＜0.25 mm 粒径最低。土壤团聚体对有机碳的贡献率主要集中在大团聚（＞0.25 mm），微团聚体（＜0.25 mm）占的比重较少。演替阶段不同对＞5 mm 粒径的团聚体有机碳贡献率具有显著影响。随着演替阶段的上升，粒径的大小与土壤有机碳贡献率的联系更加紧密，并基本表现为随粒径增加有机碳贡献率增加。比较不同团聚体有机碳含量与贡献率发现，虽然＜0.25 mm 粒径的有机碳含量高但是贡献率低，说明微团聚体对有机碳的保护力更强，且桦木+冷杉混交林＞桦木林＞冷杉原始林。

综上所述，对于川西山地天然次生林而言，大团聚体是有机碳的主要来源，而微团聚体能更加持久的稳固碳元素，比较三个演替阶段的次生林发现处于中级演替阶段的桦木+冷杉混交林对该区域森林固碳具有重要作用，而处于高级演替阶段的冷杉林的有机碳含量最高也说明了正向演替对森林固碳是具体促进作用的。