近年来，随着农业集约化程度不断提高，长期高强度利用和不合理的土壤耕作方式导致土壤结构破坏、肥力下降等耕地土壤质量问题凸显。为恢复和提高土壤地力，实现“藏粮于地”战略目标，2016年农业部出台了《探索实行耕地轮作休耕制度试点方案》^[1]。轮作是一定地块在一定耕作周期内以种类多样性和转换性种植农作物的耕作方式，也是促使土地利用自身的“内源技术”来增加土壤肥力的生态方式^[2]。轮作通过不同作物之间的时序和空间配置，更好地利用不同作物对环境、水分和养分等生态因素需求的差异性，以改善土壤结构、均衡利用土壤养分、提高土壤肥力，从而提高生态效益和经济效益^[3]。土壤团聚体是土壤结构的基本单元^[4]，其形成是土壤微生物、植物根系及复杂的物理化学过程的共同作用结果，它的数量和分布反映了土壤结构的稳定性和抗蚀性^[5]，主要受到土壤有机碳含量、种植制度和轮作模式等的影响^[6]。土壤有机碳是土壤养分转化的重要参与者，也是土壤环境质量演变的核心。研究表明，疏松熟化层表土中近90% 的土壤有机碳存在于团聚体内^[7]。土壤团聚体和有机碳相互作用^[8]，土壤团聚体为有机碳提供物理保护^[9]，有机碳促进团聚体的形成和稳定^[10]。

贵州是全国唯一没有平原支撑的典型喀斯特省份，人均耕地面积少且土壤质量低^[11]，合理的轮作技术在保持土壤养分的前提下提高有限耕地的生产力对贵州粮食安全和农业可持续发展意义重大。研究表明，不同轮作模式会对不同粒级团聚体之间的转化和再分布造成不同程度的影响^[12]，进而影响有机碳分布、土壤结构稳定性和抗侵蚀能力，且不同作物轮作下水稳性大团聚体含量、平均重量直径、各粒级团聚体有机碳含量高于单一作物种植。关于土壤有机碳在团聚体中的分布研究结果不尽一致，一部分研究认为土壤有机碳主要分布在 > 0.25 mm的大团聚体中^[13-14]，另一部分研究则认为土壤有机碳主要分布在 < 0.25 mm的微团聚体中^[15-16]。黄壤是贵州分布最广且面积最大的耕地土壤，玉米是贵州旱地种植面积最大的农作物，玉米单作或者与油菜、小麦轮作是农业生产中最主要的种植模式。目前，关于施肥对黄壤养分影响的研究较多^[17-20]，而关于不同轮作模式对黄壤养分，尤其是对黄壤团聚体影响的研究鲜见报道。因此，本研究依托国家黄壤肥力与肥效长期定位试验平台，分析比较不同轮作模式对土壤团聚体稳定性及有机碳含量的影响，旨在探寻有效改善土壤结构、维持土壤质量的轮作模式，为贵州黄壤旱地的可持续利用提供理论依据。

1 材料与方法 1.1 试验区概况

试验地位于贵州省贵阳市花溪区贵州省农业科学院内(106°39′52″ E、26°29′49″ N)，地处黔中黄壤丘陵区，海拔1 071 m，年均气温15.3℃，年均降雨量1 100 ~ 1 200 mm，年均日照时数1 354 h，相对湿度75.5%，全年无霜期270 d左右。试验地土壤为黄壤土类黄泥土土种，成土母质为三叠系灰岩与砂页岩残积物，耕层(0 ~ 20 cm)基础土样有机质36.2 g/kg、全氮1.96 g/kg、全磷0.96 g/kg、全钾10.5 g/kg、pH 6.87。

1.2 试验设计

黄壤肥力与肥效长期定位试验始于1995年，共设12个处理，本研究选取其中3个处理作为研究对象，分别是：玉米单作(maize monoculture system，MM)；小麦||绿肥–玉米轮作(wheat (intercropping green manure)-maize rotation，WMR)；油菜–玉米轮作(rape-maize rotation，RMR)。采用大区对比试验，每处理小区面积340 m²，未设重复，无灌溉设施。试验用氮肥为尿素(含N 460 g/kg)，磷肥为过磷酸钙(含P₂O₅ 160 g/kg)，钾肥为氯化钾(含K₂O 600 g/kg)。玉米季每年施N 165 kg/hm²、P₂O₅ 82.5 kg/hm²、K₂O₅ 82.5 kg/hm²，小麦季作物每年施N 82.5 kg/hm²、P₂O₅ 41.3 kg/hm²、K₂O₅ 41.3 kg/hm²。小麦和油菜每年10月播种，4月初收获。绿肥品种为光叶苕子，盛花期干物质量为1 200 kg/hm²，每年10月播种，小麦收获后翻压还田。玉米每年4月播种，9月底收获。玉米、小麦、油菜秸秆地上部分均移除，根茬还田。

1.3 样品采集和分析

由于历史原因，长期定位试验小区面积较大未设置重复，本研究将试验地延长边三等分，设置3个调查取样重复小区，于2017年玉米收获后采用梅花形5点取样法采集0 ~ 20 cm土层土样，带回实验室去除肉眼可见的植物残体和石块，再将大土块沿其自然结构掰开至1 cm左右，自然风干备用。

机械稳定性团聚体测定：干筛法^[12-15]。取400 g混合土样置于套筛(孔径依次为5、2、1、0.5、0.25 mm)顶部，用振动筛分仪(GRINDER SS200)进行振荡筛分，振幅2.0 mm，筛分时长为10 min，测定各孔径筛分后的土壤质量。

水稳性团聚体测定：湿筛法^[12-15]。按照干筛后土壤各粒级质量比称取50 g混合土样，将其置于套筛(孔径依次为5、2、1、0.5、0.25 mm)顶部，放于恒温土壤团粒分析仪(Daiki DIK-2012)的配套桶内，沿边缘缓慢加入去离子水至刻度线，确保最顶层筛的上边缘低于水面，静置10 min后开启振荡开关，以30次/min的频率振荡10 min。收集各级筛层团聚体并分别转移至铝盒中，烘干称重。

有机碳测定：重铬酸钾容量法(外加热法)^[12-15]。

1.4 计算公式

利用各粒级团聚体数据，计算平均重量直径(mean weight diameter, MWD)、几何平均直径(geometric mean diameter, GMD)、> 0.25 mm团聚体质量分数(R_0.25)^{[5, 14]}。

$ {\rm{MWD}} = \sum {\left( {{X_i} \times \frac{{{m_i}}}{{{m_t}}}} \right)} $

(1)

$ {\rm{GMD}} = {{\rm{E}}_{{\rm{xp}}}}\left[ {\frac{{\sum {({m_i} \times \ln {X_i})} }}{{\sum {{m_i}} }}} \right] $

(2)

$ {R_{0.25}} = \frac{{{m_{0.25}}}}{{{m_{\rm{t}}}}} \times 100\% $

(3)

式中：X_i为第i级团聚体平均直径(mm)；m_i为第i级团聚体质量(g)，m_t为共试土壤总质量(g)，m_0.25为 > 0.25 mm团聚体质量。

土壤团聚体破坏率(PAD)^[21]：

$ {\rm{PAD}} = \frac{{({\rm{D}}{{\rm{R}}_{0.25}} - {\rm{W}}{{\rm{R}}_{0.25}})}}{{{\rm{D}}{{\rm{R}}_{0.25}}}} \times 100\% $

(4)

式中：DR_0.25为 > 0.25 mm机械稳定性团聚体质量分数(%)；WR_0.25为 > 0.25 mm水稳性团聚体质量分数(%)。

土壤团聚体不稳定团粒指数(E_LT)^[22]：

$ {E_{{\rm{LT}}}} = \frac{{{m_t} - {m_{0.25}}}}{{{m_t}}} \times 100\% $

(5)

式中：m_0.25、m_t意义同前。

团聚体有机碳对土壤有机碳的贡献率(W)^[22]：

$ W = \frac{{o{c_i} \times {M_i}}}{{soc \times {M_t}}} \times 100\% $

(6)

式中：W为团聚体中有机碳对土壤中有机碳的贡献率(%)；OC_i为第i级团聚体中有机碳含量(g/kg)；SOC为土壤中有机碳含量(g/kg)；m_i、m_t意义同前。

1.5 统计分析

采用Excel 2010和SPSS 18.0进行数据处理和作图，采用一般线性模型进行方差分析，最小显著极差法(Duncan法)进行多重比较，显著性水平P < 0.05，文中数据为平均值±标准差(n=3)。

2 结果与分析 2.1 不同轮作模式下土壤机械稳定性团聚体组成

各轮作模式机械稳定性团聚体均以 > 5 mm的团聚体为主，且表现出随粒级减小含量(以质量分数计，下同)呈减少趋势(表 1)。与MM处理相比，WMR和RMR处理 > 5 mm的团聚体含量显著降低了20.7% 和29.0%，而5 ~ 2 mm的团聚体含量显著增加了62.67% 和64.18%，2 ~ 1 mm的团聚体含量显著增加了35.74% 和41.28%。此外，WMR和RMR处理 < 0.25 mm的团聚体含量均比MM处理降低了54.17%。。

2.2 不同轮作模式下土壤水稳定性团聚体组成

不同处理各粒级水稳性团聚体分布特征不尽相同(表 2)，MM处理以 < 0.25 mm及1 ~ 0.5 mm的团聚体为主，> 5 mm的团聚体含量最少；WMR处理以 > 5 mm的团聚体为主，其他各粒级团聚体含量相当；RMR处理各粒级团聚体含量相当。与MM处理相比，除WMR处理 > 5 mm的团聚体含量显著增加468.98%，1 ~ 0.5 mm的团聚体含量显著降低40.5% 外，不同轮作模式间其他各粒级团聚体含量无显著差异。WMR和RMR处理 < 0.25 mm的团聚体含量均比MM处理降低了35.89% 和31.25%。

2.3 不同轮作模式下土壤团聚体稳定性

MWD和GMD可表征土壤团聚体稳定性，其值越大表示平均粒径团聚度越高，稳定性越强。不同轮作模式下土壤机械稳定性显示，WMR处理GMD和R_0.25含量均较MM处理大(表 3)。水稳性团聚体MWD值以WMR处理最高，比MM处理显著提高了50%，各处理之间GMD和R_0.25值差异不显著。不同轮作模式下PAD和E_LT规律基本一致(图 1)，与MM处理相比，WMR和RMR处理的团聚体PAD分别显著降低了31.32% 和25.97%，E_LT分别显著降低了35.90% 和30.65%。可见，轮作尤其是小麦||绿肥–玉米轮作能有效提高土壤团聚体稳定性，降低土壤团聚体破坏率和不稳定团粒指数。

2.4 不同轮作模式下水稳性团聚体中有机碳分布

本研究各处理耕层有机碳含量为WMR(26.21 g/kg) > RMR(21.29 g/kg) > MM(19.26 g/kg)，且水稳性团聚体有机碳回收率均在95% 以上。不同耕作模式下土壤有机碳在不同团聚体粒径中分布存在差异(图 2)，不同轮作能不同程度地提高团聚体中有机碳含量。不同粒级团聚体中均以WMR处理的有机碳含量最高。与MM处理相比，WMR处理 > 5、1 ~ 0.5、0.5 ~ 0.25 mm团聚体有机碳含量分别显著提高了17.60%、34.41%、45.67%；RMR处理0.5 ~ 0.25 mm团聚体有机碳含量显著增加了41.32%，其他各粒级团聚体中有机碳含量无显著变化。可见，两种轮作模式均能增加各粒径团聚体尤其是大团聚体中有机碳含量，尤以小麦||绿肥–玉米轮作增幅最大。

2.5 不同轮作模式下各粒级水稳性团聚体中有机碳贡献率

不同轮作模式下不同粒级团聚体的有机碳贡献率差异较大(图 3)，其中 > 0.25 mm大团聚体中有机碳对土壤有机碳贡献率均较高，高达78% 以上，这与团聚体含量分布特征相似。> 0.25 mm大团聚体中，MM处理2 ~ 1 mm和1 ~ 0.5 mm团聚体中有机碳贡献率较高，为23.35% 和21.70%，而 > 5 mm团聚体中有机碳贡献率仅为5.55%；WMR处理 > 5 mm团聚体中有机碳贡献率高达28.73%，其他各粒级团聚体中有机碳贡献率为11.45% ~ 16.38%；RMR处理各粒级团聚体中有机碳贡献率相当，为12.29% ~ 18.93%。< 0.25 mm微团聚体中，WWR和RMR处理有机碳贡献率比MM处理降低了8.29和6.98个百分点。可见，轮作有利于提高大团聚体中尤其是 > 5 mm团聚体中有机碳贡献率，其中小麦||绿肥–玉米轮作增幅最大。

3 讨论

土壤学将当量粒径10 ~ 0.25 mm的团聚体称为大团聚体，其含量越高，说明土壤团聚性越好；< 0.25 mm的团聚体称为微团聚体，其是机械稳定性较差的团聚体，这一级别团聚体所占比重越高，表明土壤愈分散^[23]。MWD和GMD常用来反映土壤团聚体大小分布状况，GMD是对土壤团聚体主要粒级分布的描述，而MWD为各级团聚体的综合指标^[24]。一般认为，大团聚体含量、MWD和GMD值越大表示团聚体的平均粒径团聚度越高，团聚体稳定性越强^[25]。本研究结果表明，各轮作模式机械稳定性团聚体和水稳性团聚体均以大团聚体为优势团聚体，占比分别高达93.04% 及74.59% 以上；而水稳性微团聚体中 < 0.25 mm的团聚体含量比机械稳定性中 < 0.25 mm的团聚体含量高13.10 ~ 18.45个百分点，这与张钦等^[23]在黄壤上的研究结果一致。油菜–玉米轮作和小麦||绿肥–玉米轮作较玉米单作能显著提高5 ~ 2 mm和2 ~ 1 mm机械稳定性团聚体含量及 > 5 mm水稳定团聚体含量，并显著降低PAD和E_LT，尤以小麦||绿肥–玉米轮作效果更佳，说明轮作能改善团聚体组成，提高土壤团聚体稳定性，这与其他研究^{[23, 26-30]}结果一致。可见，黄壤上不同作物与绿肥间套轮作可在常规轮作模式的基础上进一步改善土壤团聚体组成并提高团聚体稳定性。李阳兵等^[31]研究表明，团聚体稳定性下降和水稳性团聚体减少的主要原因是有机质含量下降，而本研究中土壤的有机碳含量表现为小麦||绿肥–玉米轮作 > 油菜–玉米轮作 > 玉米单作，故土壤团聚体的稳定性也呈现出一致的规律。

土壤有机质与团聚体关系密切^[32-33]，其含量的提高有利于土壤结构形成，增强土壤结构稳定性，而团聚体形成反过来影响土壤有机碳分解^[34-35]。本研究结果表明，不同处理 > 0.25 mm的水稳性团聚体中有机碳对土壤有机碳贡献率高达78% 以上，说明土壤有机碳主要集中在 > 0.25 mm的水稳性团聚体中，微团聚体中含量较少，这与张凤华等^[36]研究结果一致。小麦||绿肥–玉米轮作和油菜–玉米轮作中，各粒级团聚体有机碳含量基本都高于玉米单作，尤以小麦||绿肥–玉米轮作增幅较大，说明玉米与其他作物轮作尤其是小麦季作物间套作绿肥有利于土壤团聚体的养分积累。这可能是因为不同轮作系统作物的根系分布、根系分泌物种类和数量、根系周转特征等存在差异，当不同的植物残茬增加时，可促进团聚体中颗粒有机质的形成，在有机质胶结作用下，微团聚体、颗粒有机质等相结合形成大团聚体，从而促进土壤大团聚体数量的增加并提高有机质含量^[37-40]。本研究中，轮作的冬季作物(小麦和油菜)根茬还田，增加了外源有机物的投入，从而增加了土壤有机质含量。小麦||绿肥–玉米轮作中，由于绿肥翻压还田进一步增加了外源有机物的投入，因而有机碳含量最高。另外，本研究中轮作显著提高了 > 5 mm的水稳性团聚体含量及其对土壤有机碳的贡献率，降低了土壤团聚体破坏率并提高了团聚体稳定性，稳定的团聚体能够对储存于其中的有机碳形成有效保护，且通常大团聚体能够储存更多的有机碳，故轮作可提高土壤团聚体中有机碳含量。

4 结论

轮作在改善土壤团聚体组成、提高土壤团聚体稳定性和团聚体中有机碳含量方面具有良好效果，小麦||绿肥–玉米轮作效果尤为显著。贵州冬绿肥资源丰富，在农业生产实践中，应选择合适的绿肥品种进行间套作轮作，充分发挥绿肥的培肥改土作用，提高土壤质量。小麦||绿肥–玉米轮作可作为贵州黄壤旱地较佳的轮作模式。