赵友朋, 孟苗婧, 张金池, 马洁怡, 刘胜龙,2

(1.南方现代林业协同创新中心 江苏省水土保持与生态修复重点实验室 南京林业大学 林学院,江苏 南京 210037; 2.浙江省凤阳山—百山祖国家级自然保护区凤阳山管理处, 浙江 龙泉 323700)

长期以来，土壤退化问题日益严重，土壤退化的核心是土壤结构的退化[1]，其显著特征表现在土壤团聚体粒级组成变化及其稳定性下降[2]。土壤团聚体是组成土壤结构的基本单元[3]，也是土壤肥力的物质基础[4]，土壤团聚体的稳定性很大程度上决定了土壤结构的好坏[5]。土壤中的胶结物质是土壤中最活跃的部分，是土壤具有肥力的物质基础[6]，也是土壤团聚体形成的基础和核心[7]。作为胶结物质的重要组成部分[8]，铁铝氧化物表面常带不同类型和数量的电荷[9]，可以吸附在黏粒表面而形成氧化物胶膜，促进土粒更加紧密，形成稳定的团聚体[10]，从而改变土壤结构。铁铝氧化物对土壤中有机质也具有保护作用，由此对土壤团聚体产生间接的影响[11]。目前已有较多学者[12-13]对地带性土壤铁铝氧化物进行研，但铁铝氧化物对土壤团聚体的影响的相关研究较少，而且不同学者的研究结果差异较大，聂文婷等[14]的研究结果表明，团聚体的稳定性随着游离度增大而增强；王小红等[15]的研究表明，无定形态及络合态的铁铝氧化物比游离态铁铝氧化物更能促进大团聚体的形成与稳定。不同形态铁铝氧化物对土壤团聚体如何作用仍尚无定论，关于铁铝氧化物对土壤团聚体稳定性的影响机制仍缺乏足够的研究。由于森林管理的发展，森林植被类型的变化经常发生在全球范围内。例如，在亚热带地区，自然常绿森林转化为种植园林，如混合种的种植园或竹林，是林分类型变化的例子[16]。凤阳山区自然保护区，最初是由原生态的阔叶林所主导，20世纪70年代，保护区通过间伐、飞机播种等形式进行大规模的造林，形成了以针阔混交林、常绿阔叶林、针叶林及杉木林为主的植被类型[17]。上述林分类型变化在我国东部的浙江、福建、江西等亚热带地区广泛分布，具有一定的典型性和代表性[18]。本研究以亚热带地区凤阳山国家自然保护区4种不同类型林地作为研究对象，对3个土层深度的土壤团聚体粒级分布、不同形态的铁铝氧化物含量以及土壤铁铝氧化物与土壤团聚体及其稳定性的关系进行分析，旨在探讨各形态铁铝氧化物对土壤团聚体稳定性的影响机制，为揭示亚热带地区土壤团聚体的形成与稳定机制和恢复土壤资源提供依据，并进一步丰富团聚体形成及稳定理论。

1 材料与方法

1.1 研究区概况

本研究试验区位于浙江龙泉市凤阳山国家自然保护区(119°06′—119°15′E，27°46′—27°58′N)，凤阳山保护区由福建武夷山脉向东伸展而成，属于亚热带温暖湿润气候区，雨量充沛，年降水量在2 000 mm以上，年平均气温11.8 ℃。土壤类型以黄棕壤为主，分布在800 m以上的高海拔山坡地，土层发育完整，人为干扰较少。植物资源异常丰富，植被类型较为复杂，森林覆盖率为90.8%。天然分布的木本植物主要有：木荷(Schimasuperba)、杉木(CunninghamiaLanceolata)、短柄枹(Quercusglandulifera)、马尾松(Pinusmassoniana)、黄山松(Pinustaiwanensis)、多脉青冈(Cyclobalanopsismultinervis)等。

1.2 样地调查和样品采集

2016年8月，在凤阳山保护区海拔1 300～1 400 m处，选择4种主要林分类型：阔叶林、针阔混交林、杉木林及竹林。阔叶林为人工林，主要树种是木荷和短柄枹；针阔混交林多为残存的黄山松、杉木等，以及石楠(Photiniaserrulata)、柃木(Euryajaponica)、青冈等阔叶树种；杉木林和竹林分别以杉木和毛竹为主，其他树种较少，竹林在春、冬季节可供笋用。每种林分设置3个重复样地(20 m×20 m)，每个样地内按S形设置5个样点。去除土壤表层枯落物，挖取土壤剖面，根据研究区土壤结构分0—10 cm，10—20 cm，20—30 cm取土样，并剔除树根、石块等杂物。同一样方内5个样点采集的土样按土层混合均匀后作为该样地的土壤样品，共计36份。所取土样带回实验室，经风干过筛后用于土壤理化性质测定；另用环刀和铝盒在各层取原状土样，装入硬质塑料盒中带回实验室。当土样含水量为20%左右(可塑性上限)时，将采集原状土沿自然结构掰成直径1 cm左右的小土块，在室温下风干备用。样地基本情况详见表1。

表1 研究区样地基本特征

注：BMF代表阔叶混交林; CBMF代表针阔混交林; FF代表杉木林;BF代表竹林。下同。

1.3 试验方法

(1) 土壤团聚体的测定。水稳定性团聚体的测定采用Six等[19]的方法，并稍作修改。即称取50 g风干土样在团聚体分析仪上进行湿筛分析，先将孔径为2，1，0.5，0.25，0.106，0.053 mm的套筛放入水桶中，加水至没过筛组上缘，再将称量的土样小心的放入套筛，震动30 min；振荡结束后，将留在筛子上的各级团聚体用水分别洗入铝盒中，待澄清后倒去上面的清液，于低温电热板上烘干，再在空气中平衡24 h后称量(精确至0.01 g)。 (2) 铁铝氧化物的测定：游离态铁铝氧化物(Fed，Ald)采用DCB法(连二亚硫酸钠—柠檬酸钠—重碳酸钠法)提取，无定形铁铝氧化物(Feo，Alo)用草酸胺缓冲溶液提取，络合态铁铝氧化物(Fep，Alp)采用焦磷酸钠溶液提取，提取液经稀释5倍后，用分光光度计测定[20]。

1.4 数据处理

土壤团聚体平均重量直径(MWD)计算公式[21]为：

采用Excel，SPSS 20，R语言等软件进行数据统计分析，Origin 9.0作图。

2 结果与分析

2.1 不同林分类型土壤团聚体及其稳定性分析

2.1.1 土壤团聚体粒级分布分析 水稳定性团聚体根据粒级可分为水稳定性大团聚体(>0.25 mm)和水稳定性微团聚体(<0.25 mm)[22]。由表2可以看出，阔叶混交林、针阔混交林和杉木林水稳性团聚体比例先以>2 mm粒径含量最多，范围在33.32%～45.23%；竹林水稳性团聚体比例在1～2 mm粒径含量达到最大值，为45.07%～47.49%。同一林分类型不同土层之间的各粒级水稳定性团聚体均差别不显著。水稳定性团聚体根据粒级可分为水稳定性大团聚体(>0.25 mm)和水稳定性微团聚体(<0.25 mm)[26]。4种林分类型中，>0.25 mm粒级土壤水稳性大团聚体比例均大于90%，范围在90.56%～96.34%之间，远大于<0.25粒级土壤水稳定性微团聚体。>0.25 mm的各粒级中，林分类型对>2 mm和1～2 mm粒径的水稳性团聚体有较大影响，>2 mm粒径团聚体比例在20—30 cm土层差异显著(p<0.05)，表现为杉木林>阔叶混交林>竹林，>2 mm粒径在其他土层和1～2 mm粒径在各土层均表现为阔叶混交林、针阔混交林及杉木林显著小于竹林(p<0.05)；林分类型对0.25～0.5和0.5～1 mm粒级水稳定性团聚体影响较小。<0.25 mm的各粒级水稳定性团聚体中，0.053～0.106 mm粒级土壤团聚体比例在0—10和20—30 cm土层表现为竹林显著大于杉木林(p<0.05)，其中0—10 cm土层大小依次为:竹林>针阔混交林>阔叶混交林>杉木林，20—30 cm土层大小顺序则为：竹林>阔叶混交林>针阔混交林>杉木林。

2.1.2 土壤团聚体稳定性分析 土壤团聚体稳定性是指团聚体在外部环境改变或外力作用下保持其原有形态的能力，是土壤结构和土壤性状评价的一个重要指标。平均重量直径(MWD)是评价土壤团聚体稳定性的重要指标[23]，其通过对团聚体直径取对数加权求和，弱化团聚体直径的影响份额，提高团聚体含量在评价指标体系中的影响份额，更能准确地反映土壤团聚体稳定性特征[24]，在整体研究土壤团聚体时通常运用它来综合评定土壤结构的优劣[25]。由图1可以看出，不同林分类型土壤团聚体平均质量直径MWD呈现的总体趋势为：杉木林>针阔混交林>阔叶混交林>竹林，同一林分类型土壤团聚体MWD在土层之间无显著性差异。随土层深度增加，杉木林的土壤团聚体MWD逐渐增加，阔叶混交林、针阔混交林和竹林土壤团聚体MWD均先减少后增加。竹林土壤团聚体MWD范围在1.12～1.15 mm，与其他3个林分类型的差异比较显著(p<0.05)，具体表现为：在0—10和10—20 cm土层，阔叶混交林、针阔混交林及杉木林土壤团聚体MWD显著大于竹林(p<0.05)；在20—30 cm土层，杉木林、阔叶混交林和竹林土壤团聚体MWD彼此之间差异显著(p<0.05)。

表2 研究区不同林分类型土壤水稳定性团聚体数量及粒级分布

注：同土层内同列数据后不同小写字母表示林分类型间差异显著(p<0.05)。下同。

注：误差线表示标准误差，小写字母表示同土层内林分类型差异显著(p<0.05)。下同。

图1研究区不同林分类型土壤团聚体平均质量直径(MWD)

2.2 不同林分类型土壤铁铝氧化物分析

根据铁铝氧化物存在形态的不同，可将其分为游离态(Fed，Ald)、无定形态(Feo，Alo)和络合态(Fep，Alp)3种。不同形态的铁铝氧化物对团聚体的影响作用也有所不同。同种形态的铁铝氧化物，其含量越多，作用越大。4种林分类型土壤铁铝氧化物含量中，不同形态氧化铁含量显示(图2)，Fed含量远高于Feo和Fep，Fep含量最少；不同形态氧化铝含量显示(图3)，不同形态氧化铝中，Ald含量最高，Alo次之，Alp最低。整体上表现为：Fed>Ald>Feo>Alo>Fep>Alp。不同林分类型铁铝氧化物含量在土层间的差异均不显著。随土层加深，4种林分类型各形态的氧化铁含量均增加，而各形态的氧化铝含量则未表现出一定的规律性。在0—10 cm土层，不同林分类型Fed含量呈现显著性差异(p<0.05)，表现为杉木林显著大于其他林分类型，阔叶混交林和竹林次之，针阔混交林最小，其他指标均无显著性差异。10—20 cm土层中，不同林分类型Alp含量呈现显著性差异(p<0.05)，表现为针阔混交林和杉木林远大于竹林，大小依次为：针阔混交林>杉木林>阔叶混交林>竹林，其他指标差异不显著。在0—10 cm和10—20 cm土层，Fed，Feo和Fep含量均是杉木林最大，竹林的Feo，Fep和Alp含量均是最小。20—30 cm土层中，杉木林的Fed，Feo含量最大，竹林仅在Ald含量略大于杉木林外，其他指标均是最小。

注：Fed：游离态铁铝氧化物；Feo：无定形铁铝氧化物；Fep：络合态铁铝氧化物。不同小写字母表示同形态氧化物在不同林分类型间差异显著(p<0.05)。下同。

图2研究区不同林分类型各土层土壤铁氧化物的含量

2.3 土壤铁铝氧化物与土壤团聚体稳定性关系分析

BRT是以分类回归树算法为基础发展而来[26]，BRT可产生多重回归树，提高模型稳定性及预测精度[27]。即在运算过程中，对一定量数据进行多次随机抽取，分析各自变量对因变量的影响程度，剩余数据用于拟合结果检验，最终生成的多重回归树，取平均值并输出结果。该分析方法可量化自变量对因变量的影响，分析不同影响因子的贡献率，目前被广泛应用于生态研究中。

注：Ald：游离态铁铝氧化物；Alo：无定形铁铝氧化物；Alp：络合态铁铝氧化物。

图3研究区不同林分类型各土层土壤铝氧化物的含量

运用增强回归树(Boosted regression tree)BRT模型分析不同形态铁铝氧化物对土壤团聚体稳定性主要指标MWD的影响。采用以下调整设置进行生态建模[28]：树复杂性(tree complexity)=5，学习速率(learning rate)=0.001，袋分数(bag fraction)=0.7。通过分析软件R3.3.0进行BRT分析，分析结果如图4所示，对土壤团聚体MWD影响最大的是Alp，其产生了43.6%的贡献率；其次是Ald和Fed，贡献率分别为23.5%，19.5%；Feo，Alo及Fep的影响程度较小，贡献率分别是6.9%，4%及2.4%。游离态铁铝氧化物Ald和Fed分别在17.2～18.1，27.9～29.6 g/kg浓度范围内与土壤团聚体MWD呈现正相关关系，即Ald和Fed浓度越大，土壤团聚体MWD越大；无定形铁铝氧化物Feo和Alo与土壤团聚体MWD存在一定程度的正相关，但相关性不大；络合态铁铝氧化物中Alp在0.8～1.46 g/kg浓度范围内与土壤团聚体MWD呈现正相关关系，即Alp浓度越大，土壤团聚体MWD越大，但Fep与土壤团聚体MWD的关系不明显。

图4 影响土壤团聚体MWD的主要铁铝氧化物的增强回归树BRT分析

3 讨 论

土壤团聚体的水稳定性与土壤结构关系密切[29]，水稳定性团聚体含量越大，表明土壤结构稳定性越好[30]。>0.25 mm的团聚体被称为土壤团粒结构体[22]，其含量是判定土壤质量好坏的重要指标，其含量越高，土壤结构的稳定性越大，抗蚀能力越高[21,31]。研究发现>0.25 mm粒级土壤水稳性大团聚体比例均大于90%，表明4种林分类型土壤团聚体的团聚体高，稳定性好，土壤结构优良，这与研究区土壤地质条件良好，水热资源丰富相匹配。杉木林土壤团聚体含量随着粒径的减小而降低，这与刘晓利等[32]对林分土壤团聚体的研究相同，而竹林由于会受到挖笋等人为破坏，导致大粒级团聚体的破坏，破碎为较小的粒级，随着粒径的减小呈现先增加后降低的趋势。MWD值越大，表明土壤团聚体的团聚度高，土壤团聚体稳定性就越强，土壤结构越好[33]。通常团聚体组成中大粒径的比例越高，MWD就越大，土壤团聚体也越稳定[34]。研究发现，不同林分类型土壤团聚体MWD呈现的总体趋势与水稳定性大团聚体相同，均表现为：杉木林>针阔混交林>阔叶混交林>竹林，也证实了谢贤健等[33]的研究。这可能是由于凤阳山地处亚热带，土壤中风化和成土产物的分解和淋溶水平高，有着明显的富铝化效应，盐基严重淋失，碳酸盐存在量偏少，铁铝氧化物是其重要的胶结物质，影响了林地的土壤团聚体分布及其稳定性。

铁铝氧化物包括晶质、非晶质的铁铝氧化物、氢氧化物和偏氢氧化物，是土壤中最主要、最活跃的组分之一[35]，其往往强烈吸附在黏粒表面，形成氧化物胶膜，促进土粒之间紧密结合，形成稳定的团聚体[10]。南方土壤的各形态氧化铁多以游离氧化铁的形式存在[9]，本研究结果与之相似，4种林分类型土壤均以游离氧化铁含量最高，整体上大小依次为：Fed>Ald>Feo>Alo>Fep>Alp。聂文婷等[14]的研究指出，铁铝氧化物游离度的大小与团聚体的稳定性相关，团聚体的稳定性随着游离度增大而增强，王小红等[15]的研究也表明游离态氧化铁常与黏土矿物结合，这种复合物的形成可提高土壤团聚体的稳定。本研究证实游离态铁铝氧化物Ald和Fed是促进土壤团聚体稳定的主要因素，游离态氧化铝Ald对团聚体稳定性的促进作用强于Fed，更有助于大团聚体的形成及稳定，这可能与氧化铝在土壤中表面电荷密度等有关，Al3+离子电荷/半径的比率大于Fe3+，吸附能力更强[36-37]。但对于无定形铁铝氧化物的研究结果却发现，相比无定形氧化铝Alo，无定形氧化铁Feo对土壤团聚体稳定性的促进作用更大，这是由于无定形氧化铁Feo是各形态氧化铁中羟基释放量最高的，释放的羟基可以与其他配位体进行交换进入土壤溶液中充当絮凝剂[38]，从而促进土壤颗粒之间的团聚。而总体上无定形铁铝氧化物(Alo，Feo)对土壤水稳定性大团聚体的作用不明显，这可能是无定形氧化物对大粒径土壤颗粒的团聚作用较小，更多地影响小粒径土壤颗粒或者微团聚体[39-40]。络合态铁铝氧化物中Alp对土壤团聚体稳定性的促进作用显著大于Fep，Alp尽管含量较低，但与>0.25 mm大团聚体数量及MWD呈现正相关关系，对土壤团聚体稳定性的促进作用显著，是各形态铁铝氧化物中对土壤团聚体MWD影响最大的因素。这可能与络合态铁铝氧化物的特性有关，有机络合物可以强烈地吸附在铁铝氧化物沉淀产物的表面[41]，与有机质胶结[8]或与多功能团有机络合物结合成稳定剂从而促进土壤颗粒的团聚，形成有机无机复合体提高土壤团聚体稳定性[42]。

本研究主要分析了不同形态铁铝氧化物对土壤水稳定性团聚体稳定性的影响，通过土壤团聚体稳定性指标整体研究土壤团聚体及综合评价土壤结构的优劣[8]，而对凤阳山其他土壤类型和各粒级土壤团聚体所含铁铝氧化物的含量没有涉及，因此进一步的研究中，可以结合粒级和土壤类型分析土壤团聚体所含铁铝氧化物对土壤团聚体的影响，从而更全面地探究土壤铁铝氧化物对土壤团聚体稳定性的影响机理。

4 结 论

(1) 4种林分类型土壤团聚体均以土壤水稳性大团聚体(>0.25 mm)为主，比例为90.56%～96.34%。总体上，水稳定性大团聚体含量和土壤团聚体MWD呈现的总体趋势均为：杉木林>针阔混交林>阔叶混交林>竹林，表明土壤团聚体稳定性杉木林最好，针阔混交林和阔叶混交林次之，竹林最差。

(2) 4种林分类型不同形态铁铝氧化物含量整体上大小为：Fed>Ald>Feo>Alo>Fep>Alp。土层加深，各形态氧化铁含量均增加，各形态氧化铝含量未表现出一定的规律性。在0—10 cm和10—20 cm土层，林分类型间Fed和Alp含量分别呈现显著性差异(p<0.05)。

(3) 各形态铁铝氧化物中，对土壤团聚体MWD影响程度大小依次为：Alp>Ald>Fed>Feo>Alo>Fep，游离态铁铝氧化物Ald和Fed分别在17.2～18.1 g/kg，27.9～29.6 g/kg的浓度范围内与土壤团聚体MWD呈现正相关关系，络合态氧化铝Alp在0.8～1.46 g/kg浓度范围内与土壤团聚体MWD呈现正相关关系。