高寒草地灌丛化对土壤团聚体稳定性及其铁铝氧化物分异的研究
陈红, 马文明*, 周青平, 杨智, 刘超文, 刘金秋, 杜中曼
西南民族大学青藏高原研究院,四川 成都 610041
*通信作者Corresponding author. E-mail: Mawmtf@swun.edu.cn

作者简介:陈红(1990-),女,四川泸州人,在读硕士。E-mail: 450551430@qq.com

摘要

团聚体对有机碳的物理保护是土壤碳固存的重要机制。本研究以青藏高原东缘灌丛化草地和未灌丛化草地为研究对象,采用Cambardella与Elliott湿筛法和沉降虹吸法分离和提取团聚体,测定了灌丛化草地和未灌丛化草地中土壤团聚体含量,分析了团聚体稳定性、不同形态铁铝氧化物的分异及铁铝氧化物与团聚体稳定性的关系。结果显示:灌丛化显著降低了<0.053 mm微团聚体含量( P<0.05),对团聚体稳定性没有显著影响( P>0.05);灌丛化显著增加了>0.25 mm大团聚体内游离态铝的含量(Ald的含量增加了14.09%)( P<0.05)和<0.25 mm微团聚体内游离态铁铝的含量(Fed、Ald的含量分别增加了9.75%和15.31%)( P<0.05),显著降低了<0.25 mm微团聚体内络合态铁铝的含量(Fep、Alp的含量分别降低了31.03%和9.41%)( P<0.05)。通过增强回归树分析,游离态铁铝氧化物对团聚体稳定性的作用最强,络合态铝氧化物对团聚体稳定性作用最弱。表明灌丛化有助于增强团聚体的稳定,促进土壤有机碳的固存。

关键词: 青藏高寒草地; 灌丛化; 团聚体稳定性; 铁铝氧化物; 增强回归树
Shrub encroachment effects on the stability of soil aggregates and the differentiation of Fe and Al oxides in Qinghai-Tibet alpine grassland
CHEN Hong, MA Wen-ming*, ZHOU Qing-ping, YANG Zhi, LIU Chao-wen, LIU Jin-qiu, DU Zhong-man
Qinghai-Tibet Plateau Research Institute, Southwest Minzu University, Chengdu 610041, China
Abstract

The physical protection of soil aggregates is an important mechanism of soil carbon sequestration. In this study, shrub-grassland and shrubless-grassland in East margin of Qinghai-Tibet Plateau were investigated. We used wet sieving proposed by Cambardella and Elliott and a settling siphon to separate soil aggregates, which were then measured. Data obtained included: Soil aggregate stability, differentiation of different forms of Fe and Al oxides, and the relationship between Fe and Al oxides and soil aggregate stability. It was found that: Shrub encroachment significantly reduced ( P<0.05) the proportion of microaggregates (<0.053 mm), but had no significant effect on soil aggregate stability ( P>0.05). Shrub encroachment significantly increased the content of dithionthite-citrate-bicarbonate extractable Al oxides in macroaggregates (>0.25 mm), which was increased by 14.09% ( P<0.05). Shrub encroachment also significantly increased the content of dithionthite-citrate-bicarbonate extractable Fe and Al oxides in microaggregates (<0.25 mm), which was increased by 9.75% and 15.31%, respectively ( P<0.05). However, shrub encroachment significantly reduced the content of complex Fe and Al oxides in microaggregates (<0.25 mm), the content of sodium pyrophosphate extractable Fe and Al oxides decreased by 31.03% and 9.41%, respectively ( P<0.05). Boosted regression tree analysis indicated that the free Fe and Al oxides had the strongest effect on soil aggregate stability, while the complex Al oxides had the weakest effect.The results indicate that shrub encroachment is helpful to enhance soil aggregate stability and promote the sequestration of soil organic carbon.

Keyword: Qinghai-Tibet alpine grassland; shrub encroachment; aggregate stability; Fe and Al oxides; boosted regression tree

全球草地面积约占陆地总面积的40%, 储存了全球约10%的有机碳[1], 相对于其他草地生态系统, 高寒草地生态系统对气候变化更为敏感[2]。青藏高原被称为“ 世界第三极” , 是气候变化的敏感区和启动区, 高寒草地土壤中储存着丰富的有机碳[3]。青藏高原草地土壤有机碳储量约为30~40 Gt, 占中国总土壤有机碳储量的20%以上, 占全球土壤碳储量的2%~3%[4]。如果草地有机碳库发生轻微的变化, 都将会对区域乃至全球碳循环产生重要影响[5, 6]

土壤团聚体是土壤结构的基本单元, 也是土壤碳的固存场所, 团聚体的稳定对有机碳的固存起着重要作用[7]。一般认为土壤有机碳的稳定机制主要包括3种:即团聚体的物理保护机制, 与微团聚体、黏土矿物和金属氧化物结合的化学稳定机制及抗微生物分解的生物稳定机制[8, 9]。表层土壤中约90%有机碳位于团聚体中, 而超过70%的有机碳位于< 0.053 mm微团聚体中, 被微团聚体包裹的有机碳稳定性更高[7]。铁铝氧化物是土壤中重要的无机胶体, 其表面活性强, 在有机无机复合过程中充当“ 桥” 的作用, 与有机碳紧密结合, 通过配位交换或形成离子键, 使有机碳难以被分解, 从而对团聚体的形成和稳定产生重要作用[10, 11]。同时, 有学者认为, 土壤团聚体的稳定关键在于土粒之间胶结物质的生物稳定; 有机胶结剂被微生物分解利用的速率是影响团聚体生物稳定性的重要因素[11]。但通常团聚体的稳定是有机和无机稳定剂共同作用的结果, 因此, 研究草地土壤团聚体稳定性及其影响因素的变化, 对土壤碳循环的研究具有重要意义。

草地灌丛化(shrub encroachment)是指草原生态系统中灌木或木本植物植株密度、盖度和生物量增加的现象[12], 是全球草地生态系统普遍面临的生态现象之一, 也是影响碳库储量的重要因素[13]。近年来, 在全球气候变化及人类活动的干扰下, 青藏高原高寒灌丛化现象日益严重[14, 15], 从1990-2009年已有39%的高山草甸被灌丛入侵[16]。目前关于草地灌丛化对团聚体稳定性和铁铝氧化物对团聚体稳定性的研究已多有报道[10, 16, 17, 18, 19, 20], 但关于青藏高原灌丛化条件下团聚体稳定性及铁铝氧化物分异的研究还鲜见报道。研究表明, 土壤中的铁铝氧化物在植物-土壤体系中起着重要作用[21, 22], 其不同形态之间的转化受土壤母质、气候、土壤有机碳含量等因素的影响[23, 24], 而铁铝氧化物在植物与土壤体系的转化往往与根系分泌物、微生物等因素密切相关。涂书新等[25]认为, 不同植物的根系分泌物明显不同, 不仅可以改善土壤物理结构, 促进矿物风化, 还可影响土壤矿物表面吸附性能和土壤生物学特性。王建林等[22]的研究表明植物根系和微生物等因素能使土壤中老化的氧化铁活化。黄永奎[26]指出铁铝氧化物对植物残体腐解产生的酚酸的吸附速率和容量远高于土壤矿物。贾龙[24]发现不同植被恢复下土壤微团聚体中铁铝氧化物的转化方向不同, 表明有机质是影响铁铝氧化物转化的重要因素。灌丛化是否会影响团聚体含量和稳定性及铁铝氧化物的分异, 进而影响土壤团聚体的稳定性尚未揭示。基于此, 本研究以青藏高原东缘灌丛化草地和未灌丛化草地为研究对象, 分析了团聚体含量及其稳定性和团聚体中不同形态铁铝氧化物含量的分异, 以及铁铝氧化物对团聚体稳定性的影响, 以期为高寒草地土壤碳循环研究提供理论依据。

1 材料与方法
1.1 样地概况

研究样地位于四川省阿坝州红原县西南民族大学青藏高原研究基地旁(N 32° 49'38″, E 102° 34'21″, 平均海拔3485 m), 该区域地处青藏高原东缘, 为川西北山地向高原过渡地带, 地貌类型主要由浅丘山地和丘间低地构成。该区域气候为大陆性高原气候, 常年寒冷, 四季无明显界限。日温差大、霜冻期长, 年均气温1.4 ℃, 最冷月为1月, 平均温度-36 ℃, 最暖月为7月, 平均温度26 ℃; 干湿季分明, 雨热同期, 年均降水量791 mm, 集中在5-10月, 年均蒸发量648.2 mm。采样区已灌丛化生长20年, 灌丛以高山绣线菊(Spiraea alpina)为建群种, 伴生种为窄叶鲜卑花(Sibiraea angustata), 灌木的分盖度约50%, 草本层植物种类较少, 分盖度约30%, 复合群落的总盖度约70%。未灌丛化草地主要物种为禾本科。样地土壤类型为高山灌丛草甸土, 成土母质以坡积物为主。

1.2 样品采集与处理

野外取样在2017年9月初进行。随机选取两块50 m× 50 m的无灌木生长草地群落, 在其周边(直线距离10 m)随机选取3块已灌丛化样地(50 m× 50 m)。每块灌丛化样地按对角线采样法设置3个1 m× 1 m的样方, 在每个样方的灌木丛边设置30 cm× 30 cm的样方, 灌丛化草地和未灌丛化草地的土壤理化性质见表1。用剪刀收集所有地上生物量, 用土壤铲挖取原状土(0~10 cm), 削去边缘受挤压的土壤后采集2 kg团聚体土样, 装于硬质塑料盒, 同时用环刀采集土壤容重样品。在2块随机选取的未灌丛化样地内按对角线采样法设置3个1 m× 1 m的样方, 收集样方内地上生物量, 团聚体土样采集方法同上。将采集的原状土样装于硬质塑料盒, 环刀取样后装于自封袋, 地上生物量样品装于纸质信封袋带回实验室。土样经自然风干, 期间沿纹理轻轻地将大块土掰成10~12 mm的小土块, 收集的地上生物量样品带回实验室立即归类。

表1 灌丛化草地和未灌丛化草地的土壤理化性质 Table 1 Soil physical and chemical properties of shrub-grassland and unshrub-grassland
1.3 土壤水稳性团聚体的分离提取

本研究中的团聚体指水稳性团聚体, 团聚体的分级采用Cambardella和Elliott的湿筛法[27]和沉降虹吸法。称取原状风干土样100 g, 将孔径分别为2、0.25和0.053 mm的套筛按从上到下的顺序组合好, 将称量好的土样均匀铺撒于最上层, 然后将套筛置于盛有蒸馏水的筛分桶内, 调整桶内水面高度, 缓缓湿润直至刚好淹没土样, 保持最顶层筛的上边缘始终高于水面, 在室温条件下浸润5 min后, 以30次· min-1的频率, 上下振幅为3 cm, 振荡2 min。筛分结束收集各筛层上的团聚体并分别转移至铝盒中, 50 ℃烘干称重, 获得> 2 mm、2~0.25 mm、0.25~0.053 mm三级土壤团聚体, 然后用沉降虹吸法分离沉降桶内的土壤悬液, 分离出0.053 mm和< 0.053 mm团聚体。并计算各级水稳性团聚体组成及其稳定性。团聚体稳定性采用平均重量直径(mean weight diameter, MWD)评价, 平均重量直径计算公式如下[28]:

MWD=ind̅×mi

式中:MWD为团聚体平均重量直径(mm), i为第i级团聚体; n为团聚体总级数, n=5; d̅为第i级团聚体颗粒的平均直径(mm), mi为第i级团聚体组成(%)。

1.4 土壤理化性质测定

土壤pH用酸度计测定, 水土比为2.5∶ 1, 土壤有机碳(soil organic carbon, SOC)和全氮(total nitrogen, TN)采用元素分析仪Elementar Variomax CNS Analyser (Elementar Corp., Germany)测定, 土壤容重采用环刀法测定[29]。土壤机械组成采用激光粒度仪分析(Mastersizer 2000, Malvern Corp., UK), 分为砂粒(2~0.05 mm)、粉粒(0.002~0.05 mm)、黏粒(< 0.002 mm)3个粒级[10, 29]

1.5 铁铝氧化物的提取测定

土壤团聚体中不同形态铁铝氧化物的提取和测定采用鲁如坤[29]的方法:游离态铁铝氧化物(Fed, Ald)采用DCB(连二亚硫酸钠-柠檬酸钠-重碳酸钠, dithionite-citrate-bicarbonate)法提取, 无定形态铁铝氧化物(Feo, Alo)采用草酸-草酸铵缓冲液提取, 络合态铁铝氧化物(Fep, Alp)采用焦磷酸钠溶液提取。提取、稀释后, 采用试铁灵铁铝联合比色法, 用分光光度计测定。

1.6 数据处理

数据分析采用SPSS 22.0软件, 差异显著性分析采用最小显著差异法(LSD)和邓肯法(Duncan), 增强回归树分析采用R 3.5.0软件, 采用Excel 2016作图。

增强回归树(boosted regression tree, BRT)是一种以分类回归树算法为基础发展而来的自学习方法[30]。在运算过程中, 对一定量数据进行多次随机抽取, 分析自变量对因变量的影响程度, 剩余数据用于拟合结果检验, 最后生成多重回归树, 取平均值输出结果。BRT通过产生多重回归树, 可提高模型稳定性和预测精度, 得出自变量与因变量的相互关系及贡献率, 相关计算方法目前已被广泛应用于生态研究中[31]

2 结果与分析
2.1 不同草地中土壤团聚体含量变化

灌丛化草地土壤团聚体含量表现为以2~0.25 mm粒径最高(P< 0.05), 其次是0.25~0.053 mm粒径, < 0.002 mm粒径最低(P< 0.05), 所占比例分别为37.31%、22.31%和7.22%(图1)。未灌丛化草地土壤团聚体含量表现为与灌丛化草地相似的趋势, 也表现为2~0.25 mm粒径最高(P< 0.05), < 0.002 mm粒径最低(P< 0.05), 所占比例分别为38.31%和10.14%, 但0.053~0.002 mm含量高于> 2 mm含量。进一步分析得出, 灌丛化草地中> 0.25 mm大团聚体含量显著高于< 0.25 mm微团聚体含量(P< 0.05), 而未灌丛化草地中二者差异不显著(P> 0.05)。

图1 灌丛化草地和未灌丛化草地土壤团聚体含量分布
不同小写字母表示不同粒径团聚体间差异显著(P< 0.05), 不同大写字母表示不同草地类型间差异显著(P< 0.05), 下同。
Fig.1 Distribution of soil aggregates contents in shrub-grassland and unshrub-grassland
Different lowercase letters represent significant differences among aggregates with different particle sizes (P< 0.05), different capital letters represent significant differences between different grassland types (P< 0.05), the same below.

灌丛化草地与未灌丛化草地土壤团聚体含量差异表现为:< 0.053 mm含量显著低于未灌丛化草地(P< 0.05), 且在< 0.002 mm粒径中, 灌丛化草地土壤团聚体含量比未灌丛化草地高1.4倍(P< 0.05)。

2.2 不同草地中土壤团聚体稳定性的变化

灌丛化草地> 0.25 mm比例比未灌丛化草地高1.09倍, 但两者间差异不显著(P> 0.05)(图1)。灌丛化后团聚体平均重量直径(MWD)增加了15.07%, 但灌丛化和未灌丛化草地MWD间差异也不显著(P> 0.05)(图2)。

图2 灌丛化草地和未灌丛化草地土壤团聚体平均重量直径(MWD)Fig.2 The mean weight diameter of soil aggregates in shrub-grassland and unshrub-grassland

2.3 不同草地土壤中团聚体内不同形态铁铝氧化物含量变化

灌丛化草地和未灌丛化草地土壤各粒径团聚体中铁铝氧化物的分布为:Fed > Ald> Feo> Alo> Alp> Fep(表2)。在灌丛化草地中, 铁氧化物的含量为0.05~20.15 g· kg-1, 铝氧化的含量为0.22~2.99 g· kg-1; 在未灌丛化草地中, 铁氧化物的含量为0.09~18.54 g· kg-1, 铝氧化物的含量为0.24~2.65 g· kg-1

团聚体内不同形态铁铝氧化物含量在灌丛化草地和未灌丛化草地的差异总体表现为:灌丛化显著增加了> 0.25 mm大团聚体中游离态铝及< 0.25 mm微团聚体中游离态铁铝的含量, 显著降低了< 0.25 mm微团聚体中络合态铁铝的含量。> 0.25 mm大团聚体内游离态铝的含量在灌丛化草地比未灌丛化草地增加了14.09%(P< 0.05), < 0.25 mm微团聚体内游离态铁铝的含量在灌丛化草地比未灌丛化草地分别增加了9.75%和15.31%(P< 0.05), 且0.25~0.053 mm和< 0.053 mm微团聚体内游离态铁铝的含量均表现为灌丛化草地显著高于未灌丛化草地(P< 0.05); < 0.25 mm微团聚体内络合态铁铝的含量在灌丛化草地比未灌丛化草地比分别降低了31.03%和9.41%(P< 0.05)。灌丛化草地中> 0.25 mm大团聚体内游离态铁的含量也增加了3.68%, 但差异不显著(P> 0.05); 灌丛化草地土壤团聚中络合态铁铝的含量总体呈下降趋势, 无定形态铁铝的含量没有显著变化(图3和图4)。

图3 灌丛化草地和未灌丛化草地土壤团聚体内不同形态铁氧化物分布Fig.3 Distribution of different Fe oxides contents in soil aggregates under shrub-grassland and unshrub-grassland

图4 灌丛化草地和未灌丛化草地土壤团聚体内不同形态铝氧化物分布Fig.4 Distribution of different Al oxides contents in soil aggregates under shrub-grassland and unshrub-grassland

表2 灌从化草地和未灌丛化草地土壤团聚体内铁铝氧化物含量 Table 2 The contents of Fe and Al oxides in soil aggregates in shrub-grassland and unshrub-grassland (g· kg-1)
2.4 团聚体中不同形态铁铝氧化物与团聚体稳定性的关系

运用增强回归树(BRT)模型分析不同形态铁铝氧化物对土壤团聚体稳定性主要指标MWD的影响。采用以下调整设置进行生态建模:树复杂性(tree complexity)=5, 学习速率(learning rate)=0.005, 袋分数(bag fraction)=0.5, 通过BRT分析, 最终模型经过10000次分形达到最优, 训练数据的相关性为0.996, 证数据相关性为0.992, 两者接受者操作特征曲线(receiver operating characteristic curve, ROC)值均大于0.8, 说明BRT分析结果可靠[32]

在未灌丛化草地中, 对土壤团聚体平均重量直径(MWD)影响最大的是Ald和Fed, 其贡献率分别为25.0%和18.5%; 其余形态的铁铝氧化物对MWD的贡献率变化较小, 为12.8%~15.4%。Ald的含量在> 2.5 g· kg-1后与MWD呈正相关关系, 并逐渐趋于稳定, 表明Ald的含量在> 2.5 g· kg-1后, 其含量越大, 团聚体稳定性越强, 之后不再变化(图5)。Fed的含量在17.0~18.2 g· kg-1时, 与MWD呈显著正相关关系, 此时, Fed含量越多, 团聚体越稳定。

图5 未灌丛化草地中不同形态铁铝氧化物对团聚体平均重量直径(MWD)的增强回归树分析Fig.5 The boosted regression tree analysis of different forms of Fe and Al oxides on MWD in unshrub-grassland

在灌丛化草地中, 对MWD影响最大的也是Fed和Ald, 其贡献率分别为24.3%和23.8%; 其次是Fep和Feo, 贡献率分别为18.7%和16.6%; Alo和Alp的贡献率较小, 分别为9.7%和6.9%。Fed的含量在18.5~19.5 g· kg-1时, 与MWD呈正相关关系; Ald含量在2.5~2.7 g· kg-1和3.0~3.2 g· kg-1时, 与MWD呈正相关关系(图6)。

图6 灌丛化草地中不同形态铁铝氧化物对团聚体平均重量直径(MWD)的增强回归树分析Fig.6 The boosted regression tree analysis of different forms of Fe and Al oxides on MWD in shrub-grassland

3 讨论
3.1 灌丛化对团聚体含量及其稳定性的影响

团聚体的组成和含量对土壤结构起着重要作用, 不同粒径团聚体由不同胶结物质粘合形成。Tisdall等[33]认为根系和菌丝有利于> 0.25 mm大团聚体的形成, 持久有机质和金属氧化物有利于< 0.25 mm微团聚体的形成, 而关连珠等[34]认为小粒径微团聚体主要由紧结合态有机质和黏粒胶结而成。一般认为, > 0.25 mm大团聚体比例越高, 土壤结构越好, 团聚体稳定性越强[35]。本研究区内灌丛化草地和未灌丛化草地团聚体含量均以> 0.25 mm大团聚体为主, 含量均大于50%, 其中2~0.25 mm占优势, 这与Guan等[5]在青藏高原草地土壤团聚体的研究结果一致, 表明研究区土壤结构较好。本研究中灌从化草地中> 0.25 mm团聚体含量增加, 而< 0.053 mm团聚体含量显著降低, 表明灌丛化后更多微团聚体被胶结成了大团聚体, 且团聚体稳定性增加, 大团聚体分解释放的微团聚体数量减少[36]。这可能与灌丛化后团聚体的有机、无机胶结物质含量有关。一方面, 灌丛化草地中< 0.053 mm微团聚体中游离态铁铝的含量显著增加(图2和图3), 而黏粒含量减少(表1), 这说明游离态铁铝胶结了更多黏粒形成大团聚体。同时, 灌丛化草地中根系生物量、凋落物输入量和根际活动增加, 根系和真菌菌丝进一步促进了大团聚体的形成[36]。另一方面, 灌丛化草地中> 0.25 mm大团聚体中游离态铁铝含量显著增加(图2和图3), 能胶结更多有机质, 形成稳定的大团聚体, 降低微生物分解速率, 在大团聚体中分散的微团聚体也相应减少。

MWD是评价团聚体稳定性的重要指标, MWD值越大, 团聚体稳定性越强[37]。本研究中灌丛化草地与未灌丛化草地中MWD均表现为灌丛化草地大于未灌丛化草地, 这可能与有机、无机胶结剂的数量和类型有关。苏静等[37]认为, > 0.25 mm团聚体含量与MWD成正比, 本研究也得到了一致的结果。这可能是由于灌丛化草地中凋落物数量、根系生物量、根系活动及游离态铁铝氧化物含量的增加, 提供了更多有机、无机胶结剂, 促进了土壤团聚体的形成和稳定。Bird等[38]的研究表明, 灌丛化草地土壤团聚体稳定性显著高于未灌丛化草地(P< 0.05), 但本研究中, 两种草地间团聚体稳定性差异并不显著, 这可能与有机质数量和质量及团聚体稳定机制有关。不同灌木种类对有机质的积累和微生物活性的影响存在较大差异, 灌木入侵的时间、灌木密度和人类活动等因素都也会影响有机质含量[13], 进而影响团聚体的稳定性。本研究区灌丛化草地中, 有机质和黏粒含量均低于未灌丛化草地, 但团聚体稳定性高于未灌丛化草地土壤, 表明团聚体的稳定不仅与有机质和黏粒的含量有关, 还与黏土的活性及矿物与有机质的相互作用有关[39]

3.2 灌丛化对团聚体内铁铝氧化物含量的影响

铁铝氧化物是土壤中最活跃的组分, 也是团聚体重要的无机胶结物质, 对微团聚体的形成起着更重要的作用[40]。铁铝氧化物的形态与生物、气候、土壤pH等因素有关[41, 42, 43]。本研究中, 灌丛化草地和未灌丛化草地各粒径土壤团聚体中铁铝氧化物的分布整体上均表现为:Fed> Ald > Feo> Alo> Alp> Fep, 这与赵友朋等[19]对黄棕壤和王小红等[10]对中亚热带山区红壤的研究结果相似, 且游离态铁的含量也远高于其他形态铁铝氧化物。研究发现, 游离态铁的含量可以表征土壤风化程度, 与生物和气候条件密切相关[41]。游离态氧化铁往往在雨热充沛的热带、亚热带气候下较为多见[41], 但研究区气候为大陆性高原气候, 其含量高可能主要受青藏高原冻融交替作用的影响。有学者认为, 冻融对土壤结构影响剧烈, 使微生物活性增强, 促进了土壤中硅酸盐矿物的风化, 从而增加游离态铁的含量[44]。pH也是影响铁氧化物形态的重要因素之一, 研究区土壤pH为弱酸性, 有利于结晶良好的赤铁矿形成[45]。本研究中< 0.25 mm微团聚体内各形态铁铝氧化物的含量显著高于> 0.25 mm大团聚体中, 表明铁铝氧化物主要富集在微团聚体中, 这与已有的一些研究结果一致。衡利沙等[46]对水稻土黄棕壤团聚体中铁铝氧化物含量的研究表明, 铁铝氧化物的含量随粒级降低而增加, 且主要集中在微团聚体中。贾龙[24]对退化红壤区植被恢复土壤的研究表明, 不同恢复模式下微团聚体中不同形态铁铝氧化物主要向细粒级富集。这可能与微团聚体中阳离子交换量更大和表面负电荷的数量更多有关[47]。此外, 微团聚体颗粒组表面上更多以多层吸附为主, 吸附黏粒的能力更强[47], 从而也更有利于铁铝氧化物的富集。

土壤类型、气候带、植被覆盖等均会影响铁铝氧化物的组成和含量, 而不同形态的铁铝氧化物对团聚体的形成和稳定影响程度不同[10, 19, 48]。研究表明, pH、有机质、微生物、植被等都会影响土壤中铁氧化物的形态转化[43]。灌丛化对铁铝氧化物含量的改变可能与有机质含量、根系活动和微生物等因素有关。铁氧化物形态的转化可概括为活化和老化两个过程, 用Feo/Fed表示氧化铁的活化度, 老化是由离子态-非晶型态-晶型态的过程, 活化即相反过程[41]。本研究中灌丛化草地团聚体内游离态铁铝的含量均增加, 而团聚体内无定态铁铝氧化物的含量没有显著变化, 说明灌丛化后铁铝氧化物的活化度降低, 趋向于结晶老化。氧化铁在土壤-植物体系中老化和活化过程在一定条件下可以相互转化, 老化的氧化铁在有机质和其他因素的作用下可被活化[22]。有学者认为, 矿物的晶相变化可能与根系活动引起的土壤pH改变、微生物数量和种类增加及根系分泌的有机酸有重要关系[25]。本研究中, 灌丛化后铁铝氧化物趋于老化, 这可能由于灌丛化草地土壤孔隙度更低, 微生物活性和数量降低[49], 活性有机质减少[22]。其次, 灌丛和草本植物根系分泌物的不同也可能导致铁铝氧化物的晶相变化; 此外, 根系分泌的有机酸还可能引起根际pH的改变, 从而导致灌丛化草地和未灌丛化草地中铁铝氧化物的形态转化[21]。灌丛化显著降低了< 0.25 mm微团聚体中络合态铁铝氧化物的含量, 且> 0.25 mm团聚体中络合态铁铝氧化物含量也呈下降趋势, 这可能与灌丛化草地土壤中有机质含量降低有关。研究表明, 有机质含量与络合态铁铝氧化物含量存在显著正相关关系[50], 而络合态铁铝本身为有机质与铁铝氧化物的胶结[10], 本研究中灌丛化草地土壤中有机质含量比未灌丛化草地低(表1), 所以其络合态铁铝氧化物含量也降低。总体而言, 灌丛化后团聚体中铁铝氧化物趋于结晶老化, 这一定程度上促进了团聚体的稳定[51]

3.3 不同形态的铁铝氧化物与团聚体稳定性的关系

铁铝氧化物对团聚体稳定性的机制主要有3种:在土壤溶液中充当絮凝剂; 充当黏粒和有机分子的胶结剂; 作为凝胶在黏粒表面沉淀[52]。铁铝氧化物的形态、结晶度、不同植被覆盖均会影响铁铝氧化物对团聚体的稳定性[19]。王小红等[10]认为, 铁铝氧化物的结晶度、比表面积、表面电荷等特性的差异会导致其结合腐殖质和黏土矿物的能力不同, 进而影响团聚体的稳定性。但谭文峰等[52]指出, 铁铝氧化物对不同类型土壤、矿物的电荷性质和表面积影响不同。本研究中, 灌丛化和未灌丛化草地, 不同形态的铁铝氧化物对团聚体稳定性的贡献率均表现为Fed、Ald的贡献率最高, Alp贡献率最低, 表明游离态铁铝对团聚体稳定性的作用最强, 络合态铝对团聚体稳定性作用最弱。这可能与本研究区土壤中游离态铁铝氧化物对黏粒和有机质的胶结能力更强有关。胡国成等[53]认为, 游离态铁氧化物对土壤颗粒具有较强的胶结能力, 通过强烈吸附在黏粒表面形成氧化物胶膜, 促进团聚体的稳定。黏粒对团聚体的影响取决于其黏土矿物的组成, 研究区土壤pH为酸性, 氧化铁主要带正电荷, 易与带负电荷的黏土矿物结合, 形成稳定的团聚体[52]。张琪[54]的研究表明, 游离态铁铝氧化物与黏粒直接存在极显著正相关关系(P< 0.01, 相关系数分别为0.98和0.95), 且与MWD相关性最强, 而无定形态和络合态铁铝氧化物与黏粒和MWD的相关性较弱或无显著相关性。游离态铁铝对有机质也有重要作用, Li等[55]在青藏高原草甸土壤的研究结果显示, 游离态铁铝与有机碳的相关关系大于无定形态铁铝, 表明游离态铁铝氧化物对有机碳的固存起着更重要的作用。但在Li等[55]的研究中, 无定形态铁铝与MWD的相关性却比游离态铁铝更强, 这可能与游离态铁铝氧化物的含量较低有关[56]。有学者认为[10, 57, 58], 无定形态和络合态铁铝氧化物更有利于团聚体的稳定性, 本研究结果与之相反, 这可能与研究区土壤性质、黏土矿物、土壤溶液的离子组成等因素有关[59]。此外, 不同pH下, 金属离子与有机质的配位方式也不同, 有机质的性质和种类也会影响金属矿物对有机质的吸附[60]。络合态铁对MWD的贡献率大于络合态铝, 表明络合态铁氧化物对团聚体稳定性作用比络合态铝氧化物更强, 这与它们和有机质结合的程度有关。研究表明, 络合态铁铝氧化物与有机质呈正相关关系, 络合态铝氧化物与有机质结合的牢固程度比络合态铁氧化物低[42]。铁铝氧化物与团聚体稳定性的关系, 还有待于进一步研究。

4 结论

本研究分析了青藏高原东缘灌丛化和未灌丛草地土壤团聚体稳定性及其铁铝氧化物的分异, 该区域土壤团聚体以2~0.25 mm为主, 说明其土壤结构较好; 灌丛化显著降低了< 0.053 mm微团聚体的含量, 但并未影响团聚体稳定性。各形态铁铝氧化物主要分布在< 0.25 mm微团聚体中; 灌丛化显著增加了> 0.25 mm大团聚体中游离态铝和< 0.25 mm微团聚体中游离态铁铝的含量, 显著降低了< 0.25 mm微团聚体中络合态铁铝的含量, 对无定形态铁铝氧化物含量没有显著影响, 说明灌丛化草地土壤中铁铝氧化物趋向于结晶老化。通过增强回归树分析了铁铝氧化物与团聚体稳定性的关系, 表明游离态铁铝氧化物对团聚体的稳定性作用最强, 络合态铝氧化物对团聚体稳定性的作用最弱。灌丛化后, 有助于增强团聚体稳定性, 促进有机碳的固存。

参考文献
[1] Dlamini P, Chivenge P, Chaplot V. Overgrazing decreases soil organic carbon stocks the most under dry climates and low soil pH: A meta-analysis shows. Agriculture, Ecosystems : Environment, 2016, 221: 258-269. [本文引用:1]
[2] Zhang Y, Dong S, Gao Q, et al. Climate change and human activities altered the diversity and composition of soil microbial community in alpine grasslands of the Qinghai-Tibetan Plateau. Science of the Total Environment, 2016, 562: 353-363. [本文引用:1]
[3] Wang L, Zeng H, Zhang Y J, et al. A review of research on soil carbon storage and its influencing factors in the Tibetan Plateau. Chinese Journal of Ecology, 2019, 38(11): 3506-3515.
王荔, 曾辉, 张扬建, . 青藏高原土壤碳储量及其影响因素研究进展. 生态学杂志, 2019, 38(11): 3506-3515. [本文引用:1]
[4] Shang W, Wu X D, Zhao L, et al. Seasonal variations in labile soil organic matter fractions in permafrost soils with different vegetation types in the central Qinghai-Tibet Plateau. Catena, 2016, 137: 670-678. [本文引用:1]
[5] Guan S, An N, Zong N, et al. Climate warming impacts on soil organic carbon fractions and aggregate stability in a Tibetan alpine meadow. Soil Biology and Biochemistry, 2018, 116: 224-236. [本文引用:2]
[6] Qin Y, Yi S H, Li N J, et al. Advance in studies of carbon cycling on alpine grassland s of the Qinghai-Tibetan Plateau. Acta Prataculturae Sinica, 2012, 21(6): 278-288.
秦彧, 宜树华, 李乃杰, . 青藏高原草地生态系统碳循环研究进展. 草业学报, 2012, 21(6): 278-288. [本文引用:1]
[7] Liu Z L, Yu W T. Review of researches on soil aggregate and soil organic carbon. Chinese Journal of Eco-Agriculture, 2011, 19(2): 447-455.
刘中良, 宇万太. 土壤团聚体中有机碳研究进展. 中国生态农业学报, 2011, 19(2): 447-455. [本文引用:2]
[8] Sarkhot D V, Comerford N B, Jokela E, et al. Aggregation and aggregate carbon in a forested southeastern coastal plain spodosol. Soil Science Society of America Journal, 2007, 71(6): 1779-1787. [本文引用:1]
[9] Six J, Conant R T, Paul E A, et al. Stabilization mechanisms of soil organic matter: Implications for C-saturation of soils. Plant and Soil, 2002, 241(2): 155-176. [本文引用:1]
[10] Wang X H, Yang Z J, Liu X F, et al. Effects of different forms of Fe and Al oxides on soil aggregate stability in mid-subtropical mountainous area of southern China. Acta Ecologica Sinica, 2016, 36(9): 2588-2596.
王小红, 杨智杰, 刘小飞, . 中亚热带山区土壤不同形态铁铝氧化物对团聚体稳定性的影响. 生态学报, 2016, 36(9): 2588-2596. [本文引用:8]
[11] Xu X R, Wang J K. A review on different stabilized mechanisms of soil aggregates and organic carbon. Chinese Journal of Soil Science, 2017, 48(6): 1523-1529.
徐香茹, 汪景宽. 土壤团聚体与有机碳稳定机制的研究进展. 土壤通报, 2017, 48(6): 1523-1529. [本文引用:2]
[12] Van Auken O W. Causes and consequences of woody plant encroachment into western North American grasslands. Journal of Environmental Management, 2009, 90(10): 2931-2942. [本文引用:1]
[13] Li H, Shen H H, Chen L Y, et al. Effects of shrub encroachment on soil organic carbon in global grasslands. Scientific Reports, 2016, 6(1): 28974. [本文引用:2]
[14] Gao Q Z, Guo Y Q, Xu H G, , et al. Climate change and its impacts on vegetation distribution and net primary productivity of the 16 alpine ecosystem in the Qinghai-Tibetan Plateau. Science of the Total Environment, 2016, 554/555: 34-41. [本文引用:1]
[15] Wang Y X, Chen X J, Lou S N, et al. Woody-plant encroachment in grassland s: A review of mechanisms and aftereffects. Acta Prataculturae Sinica, 2018, 27(5): 219-227.
王迎新, 陈先江, 娄珊宁, . 草原灌丛化入侵: 过程、机制和效应. 草业学报, 2018, 27(5): 219-227. [本文引用:1]
[16] Brand t J S, Haynes M A, Kuemmerle T, et al. Regime shift on the roof of the world: Alpine meadows converting to shrublands in the southern Himalayas. Biological Conservation, 2013, 158(2): 116-127. [本文引用:2]
[17] Erfanzadeh R, Bahrami B, Motamedi J, et al. Changes in soil organic matter driven by shifts in co-dominant plant species in a grassland. Geoderma, 2014, 213(1): 74-78. [本文引用:1]
[18] Hewins D B, Sinsabaugh R L, Archer S R, et al. Soil-litter mixing and microbial activity mediate decomposition and soil aggregate formation in a sand y shrub-invaded Chihuahuan Desert grassland. Plant Ecology, 2017, 218(4): 459-474. [本文引用:1]
[19] Zhao Y P, Meng M J, Zhang J C, et al. Relationship between soil aggregate stability and different forms of Fe and Al oxides in different forest types. Bulletin of Soil and Water Conservation, 2018, 38(4): 75-81, 86.
赵友朋, 孟苗婧, 张金池, . 不同林地类型土壤团聚体稳定性与铁铝氧化物的关系. 水土保持通报, 2018, 38(4): 75-81, 86. [本文引用:4]
[20] Zhu X G. Response of soil properties changes in typical steppe to woody plant encroachment. Beijing: Institute of Botany, the Chinese Academy of Science, 2007.
祝晓光. 典型草原土壤性状变化对灌丛化的响应. 北京: 中国科学院植物所, 2007. [本文引用:1]
[21] Hu S P. Molecular speciation of heavy metals transformation in soil. Hangzhou: Zhejiang University, 2009.
胡少平. 土壤重金属迁移转化的分子形态研究. 杭州: 浙江大学, 2009. [本文引用:2]
[22] Wang J L, Liu Z Y. Chemical behavior of heavy metals in rhizosphere Ⅲ. The rhizosphere effect of iron oxide transformation in soils and its ecological significance. Chinese Journal of Applied Ecology, 1992, 3(4): 339-345.
王建林, 刘芷宇. 重金属在根际中的化学行为Ⅲ. 土壤中铁形态. 应用生态学报, 1992, 3(4): 339-345. [本文引用:4]
[23] Guo X M, Wu H H, Luo M, et al. The morphological change of Fe/Al-oxide minerals in red soils in the process of acidification and its environmental significance. Acta Petrologica Et Mineralogica, 2007, 26(6): 515-521.
郭杏妹, 吴宏海, 罗媚, . 红壤酸化过程中铁铝氧化物矿物形态变化及其环境意义. 岩石矿物学杂志, 2007, 26(6): 515-521. [本文引用:1]
[24] Jia L. Effects of vegetation restoration on microaggregates and their cementing substances in degraded red soil. Nanchang: Nanchang Institute of Technology, 2018.
贾龙. 植被恢复对退化红壤微团聚体及其胶结物质的影响. 南昌: 南昌工程学院, 2018. [本文引用:3]
[25] Tu S X, Sun J H, Guo Z F, et al. On relationship between root exudates and plant nutrition in rhizosphere. Soil and Environmental Sciences, 2000, 9(1): 64-67.
涂书新, 孙锦荷, 郭智芬, . 植物根系分泌物与根际营养关系评述. 土壤与环境, 2000, 9(1): 64-67. [本文引用:2]
[26] Huang Y K. Chemical abstracts of humic acids. Beijing: Science Press, 1982.
黄永奎. 腐植酸化学文摘. 北京: 科学出版社, 1982. [本文引用:1]
[27] Cambardella C A, Elliott E T. Carbon and nitrogen distribution in aggregates from cultivated and native grassland soils. Soil Science Society of America Journal, 1993, 57(4): 1071-1076. [本文引用:1]
[28] Barral M T, Bujan E, Devesa R, et al. Comparison of the structural stability of pasture and cultivated soils. The Science of the Total Environment, 2007, 378(1/2): 174-178. [本文引用:1]
[29] Lu R K. The analysis method of soil agricultural chemistry. Beijing: China Agricultural Science and Technology Press, 2000.
鲁如坤. 土壤农业化学分析方法. 北京: 中国农业科技出版社, 2000. [本文引用:3]
[30] She Y P, Chen S Q. Forecast and analysis of algae based on boosted regression tree. Journal of Changchun University, 2015, 25(6): 20-23.
佘玉萍, 陈淑清. 基于增强回归树的水藻预测分析. 长春大学学报, 2015, 25(6): 20-23. [本文引用:1]
[31] Wu Y W, Zhang J C, Guo X P, et al. Identification of efficient strain applied to mining rehabilitation and its rock corrosion mechanism: Based on boosted regression tree analysis. Environmental Science, 2017, 38(1): 283-293.
吴雁雯, 张金池, 郭晓平, . 应用于矿山修复的高效菌株鉴定与溶岩机制: 基于增强回归树分析. 环境科学, 2017, 38(1): 283-293. [本文引用:1]
[32] Elith J, Leathwick J R, Hastie T. A working guide to boosted regression trees. Journal of Animal Ecology, 2008, 77(4): 802-813. [本文引用:1]
[33] Tisdall J M, Oades J M. Organic matter and water-stable aggregates in soils. Journal of Soil Science, 1982, 33(2): 141-163. [本文引用:1]
[34] Guan L Z, Zhang B Q, Yan L. Composition of microaggregate and cementing substances in black soils and brown forest soils with different fertility levels. Acta Pedologica Sinica, 1991, 28(3): 260-267.
关连珠, 张伯泉, 颜丽. 不同肥力黑土, 棕壤微团聚体组成及其胶结物质的研究. 土壤学报, 1991, 28(3): 260-267. [本文引用:1]
[35] Nie W T, Hu B. The research progress of soil aggregate stability and its influence factors//Proceedings of the cross-strait symposium on soil and water conservation, Chinese society for soil and water conservation. Wuhan: Chinese Society of Soil and Water Conservation, 2014: 48-52.
聂文婷, 胡波. 土壤团聚体稳定性及其影响因素研究进展//中国水土保持学会海峡两岸水土保持学术研讨会论文集. 武汉: 中国水土保持学会, 2014: 48-52. [本文引用:1]
[36] Liao J D, Boutton T W, Jastrow J D. Organic matter turnover in soil physical fractions following woody plant invasion of grassland : Evidence from natural 13C and 15N. Soil Biology : Biochemistry, 2006, 38(11): 3197-3210. [本文引用:2]
[37] Su J, Zhao S W. Comparison of the analysis methods for soil aggregate stability. Bulletin of Soil and Water Conservation, 2009, 29(5): 114-117.
苏静, 赵世伟. 土壤团聚体稳定性评价方法比较. 水土保持通报, 2009, 29(5): 114-117. [本文引用:2]
[38] Bird S B, Herrick J E, Wand er M M, et al. Spatial heterogeneity of aggregate stability and soil carbon in semi-arid rangeland . Environmental Pollution, 2002, 116(3): 445-455. [本文引用:1]
[39] Rasmussen C, White D A, Bernard A W, et al. Vegetation effects on soil organic carbon quality in an arid hyperthermic ecosystem. Geoderma, 2010, 175(9): 438-446. [本文引用:1]
[40] Chen S. Stability of soil aggregates under different land use patterns and its relationships with organic matter and iron-alumina oxides. Wuhan: Huazhong Agricultural University, 2012.
陈山. 不同利用方式土壤团聚体稳定性及其与有机质和铁铝氧化物的关系. 武汉: 华中农业大学, 2012. [本文引用:1]
[41] Chen J F. Oxides in soil colloids. Chinese Journal of Soil Science, 1981, 18(2): 44-49.
陈家坊. 土壤胶体中的氧化物. 土壤通报, 1981, 18(2): 44-49. [本文引用:4]
[42] Xiong Y. Soil colloid, Book one: Material basis of soil colloid. Beijing: Science Press, 1983.
熊毅. 土壤胶体, 第一册: 土壤胶体中的物质基础. 北京: 科学出版社, 1983. [本文引用:2]
[43] Yu X L. Study on transformation of different forms of carbon and the coupling relationships with iron in Momoge wetland soil of Songnen Plain. Wuhan: Huazhong Agricultural University, 2016.
于秀丽. 松嫩平原莫莫格湿地土壤不同形态碳转化过程及其与土壤铁耦合关系研究. 武汉: 华中农业大学, 2016. [本文引用:2]
[44] Wang Z, Zhang Y L, Zou H T, et al. The effect of the freezing/thawing cycles on iron oxide forms in the soil. Chinese Journal of Soil Science, 2014, 45(1): 1-5.
王展, 张玉龙, 邹洪涛, . 冻融交替对土壤中不同形态氧化铁的影响. 土壤通报, 2014, 45(1): 1-5. [本文引用:1]
[45] Cudennec Y, André Lecerf. The transformation of ferrihydrite into goethite or hematite, revisited. Journal of Solid State Chemistry, 2006, 179(3): 716-722. [本文引用:1]
[46] Heng L S, Wang D Z, Jiang X, et al. Relationship between Fe, Al oxides and stable organic carbon, nitrogen in the yellow-brown soils. Environmental Science, 2010, 31(11): 2748-2755.
衡利沙, 王代长, 蒋新, . 黄棕壤铁铝氧化物与土壤稳定性有机碳和氮的关系. 环境科学, 2010, 31(11): 2748-2755. [本文引用:1]
[47] Dong C X. Study on the surface properties and adsorption and desorption of copper ions in the micro-aggragates of paddy soil (huangnitu). Nanjing: Nanjing Agricultural University, 2007.
董长勋. 水稻土(黄泥土)微团聚体表面性质及对铜离子吸附与解吸特性研究. 南京: 南京农业大学, 2007. [本文引用:2]
[48] Gu Z Y, Kang L, Luo M J, et al. Factors affecting stability of soil aggregate in eastern hunan Province. Bulletin of Soil and Water Conservation, 2018, 38(5): 58-63, 69.
谷忠元, 康黎, 罗梦娟, . 湘东地区典型土壤团聚体稳定性的影响因素. 水土保持通报, 2018, 38(5): 58-63, 69. [本文引用:1]
[49] Cheng L, Qin J H, Zhang L C, et al. Application of Le bissonnais method to study soil aggregate stability in red soils under different types of vegetation. http: //kns. cnki. net/kcms/detail/32. 1119. P. 20190820. 1619. 008. html, 2019-8-20.
程谅, 秦嘉惠, 张利超, . 应用Le Bissonnais法研究不同植被类型下红壤团聚体稳定性. http://kns.cnki.net/kcms/detail/32.1119.P.20190820.1619.008.html, 2019-8-20. [本文引用:1]
[50] Wei C F, Chen S Z. The composition of Fe and Al oxides in the organomineral complexes of purple paddy soils in Sichuan Basin. Journal of Southwest University, 1991, 13(5): 518-526.
魏朝富, 陈世正. 四川盆地紫色水稻土有机无机复合体铁铝氧化物组成的研究. 西南农业大学学报, 1991, 13(5): 518-526. [本文引用:1]
[51] He Q, Chen J F, Xu Z Y. Influence of transformation of iron oxides on soil structure. Acta Pedologica Sinica, 1981, 18(4): 326-334.
何群, 陈家坊, 许祖诒. 土壤中氧化铁的转化及其对土壤结构的影响. 土壤学报, 1981, 18(4): 326-334. [本文引用:1]
[52] Tan W F, Zhou S Z, Liu F, et al. Advancement in the study on interactions between iron-aluminum (hydro-) oxides and clay minerals in Soil. Soils, 2007, 39(5): 726-730.
谭文峰, 周素珍, 刘凡, . 土壤中铁铝氧化物与黏土矿物交互作用的研究进展. 土壤, 2007, 39(5): 726-730. [本文引用:3]
[53] Hu G C, Zhang M K. Mineralogical evidence for strong cementation of soil particles by iron oxides. Chinese Journal of Soil Science, 2002, 33(1): 25-27.
胡国成, 章明奎. 氧化铁对土粒强胶结作用的矿物学证据. 土壤通报, 2002, 33(1): 25-27. [本文引用:1]
[54] Zhang Q. The interaction mechanism of aggregates characteristics and physic-chemical properties of red soils and their effect on erosion process. Wuhan: Huazhong Agricultural University, 2003.
张琪. 红壤团聚体特征与物理化学性质相互作用机理及其对侵蚀过程的影响. 武汉: 华中农业大学, 2003. [本文引用:1]
[55] Li C L, Cao Z Y, Chang J J, et al. Elevational gradient affect functional fractions of soil organic carbon and aggregates stability in a Tibetan alpine meadow. Catena, 2017, 156: 139-148. [本文引用:2]
[56] Barral M T, Arias M, Guérif J. Effects of iron and organic matter on the porosity and structural stability of soil aggregates. Soil and Tillage Research, 1998, 46(3/4): 261-272. [本文引用:1]
[57] Igwe C A, Zarei M, Stahr K. Colloidal stability in some tropical soils of southeastern Nigeria as affected by iron and aluminium oxides. Catena, 2009, 77(3): 232-237. [本文引用:1]
[58] Duiker S W, Rhoton F E, Torrent J, et al. Iron (hydr) oxide crystallinity effects on soil aggregation. Soil Science Society of America Journal, 2003, 67(2): 606-611. [本文引用:1]
[59] Liu D. Composition of clay minerals in aggregates and effect on aggregate stability of red soils under different erosion. Wuhan: Huazhong Agricultural University, 2007.
刘冬. 不同侵蚀度红壤团聚体中粘粒矿物分布特征及其对团聚体稳定性的影响. 武汉: 华中农业大学, 2007. [本文引用:1]
[60] Wei S Y, Tan W F, Liu F. Advances on the study of mineral-humus interactions in soils. Soils and Fertilizers Sciences in China, 2009, (1): 1-6.
魏世勇, 谭文峰, 刘凡. 土壤腐殖质-矿物质交互作用的机制及研究进展. 中国土壤与肥料, 2009, (1): 1-6. [本文引用:1]