徐 虎，蔡岸冬，周怀平，Colinet Gilles，张文菊，徐明岗*

（1 中国农业科学院农业资源与农业区划研究所/耕地培育技术国家工程实验室，北京 100081；2 Gembloux Agro-Bio Tech，University of Liege，Passage des deportes 2，Gembloux 5030，Belgium；3 山西农业大学资源环境学院，太原 030031；4 中国农业科学院农业环境与可持续发展研究所，北京 100081）

土壤有机碳 (SOC) 作为农田肥力的重要指标，是农作物高产稳产的基础[1]。土壤有机碳库也是陆地生态系统最大的碳库，小幅度的变化就会对全球气候产生巨大影响[2]。SOC不仅对土壤肥力和作物生产力有直接影响，在缓解气候变化方面也有极其重要作用[1-2]。我国农作物秸秆资源十分丰富，2015年农作物秸秆总量约为9.31 × 108t[3]。作物秸秆含丰富有机碳和养分，秸秆还田既可减少就地焚烧产生的温室气体[4]，也能增加农田碳输入和作物归还，有效维持和提升土壤有机碳储量[5–7]。因此，秸秆还田在土壤有机碳固定方面有巨大潜力[8]。Yang等[9]研究表明，5年秸秆还田0—30 cm土层SOC提高4.45%～9.00%。刘玲等[5]发现，秸秆粉碎还田和过腹还田下黑土表层SOC分别增加28.70%和12.40%。目前，秸秆还田对有机碳的影响研究多集中在表层土壤[5–7,9]，然而底层SOC储量是表层的2～3倍，所以探究底层SOC对秸秆还田的响应对深入理解农田碳循环和缓解气候变化也非常重要[10]。SOC提升是一个长期缓慢的过程，因此，只有在长时间尺度的研究才能获得更可靠结果[11]。马力等[12]研究表明，17年秸秆还田增加红壤水稻土0—45 cm土层SOC。解丽娟等[13]研究发现，20年秸秆还田配施化肥显著增加黑土20—40 cm 土层 SOC (41%)。然而，徐虎等[14]在旱地红壤和盖霞普等[15]在潮土研究表明，长期秸秆还田可显著增加40 cm以下土层SOC。褐土作为我国主要土壤类型之一，总面积约 2516 × 104hm2，该地区是我国小麦和玉米主产区之一，作物秸秆资源十分丰富，探究农田剖面土壤有机碳对秸秆还田响应的特征，对该地区农作物秸秆的高效利用、解决粮食安全及农业环境可持续性问题等具有指导意义[11]。此外，褐土多具团粒结构且土质疏松，其底层SOC对秸秆还田的响应特征还不明确，尤其是不同方式秸秆还田[11,16]。因此，本研究依托山西省寿阳县长期(21年) 秸秆还田试验平台，通过对不同秸秆还田方式下褐土0—100 cm剖面SOC和土壤养分及作物产量开展研究，以期明确不同秸秆还田方式下褐土剖面有机碳的变化特征及影响因素，为该地区秸秆还田措施优化、农田培肥和固碳减排等提供科学依据。

1 材料与方法

1.1 试验地概况

长期秸秆还田试验位于我国山西省寿阳县 (东经113°06'、北纬37°58')，地处中纬度暖温带半湿润偏旱区，年均温为7.4℃，年均降水量为501 mm (降雨集中在每年6～9月份) 且年际间波动较大。供试土壤为褐土，成土母质为马兰黄土，试验初始 (1992年)表层 (0—20 cm) 土壤理化指标：有机碳 15.71 g/kg、全氮 1.07 g/kg、全磷 0.76 g/kg、全钾 23.7 g/kg、有效磷 4.84 mg/kg、速效钾 95.0 mg/kg和pH 8.30[11]。

1.2 试验设计与田间管理

试验开始于1992年，供试作物为春玉米，种植制度为一年一熟，玉米生育期内 (4月15—25日播种，9月20日—10月10日收获)，主要田间管理措施是除草和防治病虫害。试验采用裂区设计，主区为化学氮磷肥春季施用和秋季施用，副区为4个秸秆还田方式处理：1) 秸秆不还田 (CK)；2) 秸秆覆盖还田(SM)；3) 秸秆粉碎后直接还田 (SC)；4) 秸秆过腹还田 (CM，腐熟鲜牛粪)。小区面积为54 m2，无重复。供试化肥为尿素 (含 N 46%)和过磷酸钙 (含 P2O514%)。玉米秸秆 (风干) 有机碳、全氮、全磷 (P2O5)和全钾(K2O) 含量分别为443.0、7.39、0.44和27.50 g/kg。CM 处理的牛粪 (风干) 有机碳、全氮、全磷 (P2O5)和全钾 (K2O) 含量分别为52.50、3.93、1.37和14.10 g/kg[11]。

化肥春季施用是结合春播进行化肥穴施或浅条施 (深度为4—7 cm)，化肥秋季施用是结合秋季深耕翻地进行化肥条施或全耕层深施 (深度10—30 cm)，生育期内不再追肥。SM处理是每年5月下旬，将前茬全部玉米秸秆均匀覆盖地表，当季玉米收获后，再将未腐解秸秆深翻还田；SC处理是当季玉米收获后，直接将秸秆铡碎 (长 15 cm) 结合深耕还田；CM处理是将前茬玉米秸秆 (与SM和SC处理秸秆等量) 喂牛，将产生的牛粪腐熟，在当季玉米收获后均匀撒施并深翻还田[11]。每个处理具体化肥施用量和秸秆还田量详见表1。

表 1 长期定位试验各处理年化肥施用量和秸秆还田量Table 1 Annual input of chemical fertilizers and different straw incorporation methods in long-term experiment

1.3 测定项目与方法

在2013年玉米收获后采集土壤剖面样品。采样方法为：在每个小区选3个代表性样点，用土钻分5 个层次 (每层 20 cm) 采集 0—100 cm 剖面样品。土壤样品风干后研磨过2 mm筛，测定有效磷和有效钾；磨细过0.25 mm筛，测定土壤有机碳、全氮、全磷和全钾。各指标测定方法如下：土壤有机碳和全氮含量分别采用重铬酸钾容量法和半微量凯氏法，土壤全磷和全钾含量采用氢氧化钠熔融法，土壤有效磷和速效钾含量分别采用Olsen法和乙酸铵浸提—火焰光度法，土壤容重采用环刀法。同时，在每年玉米收获后，测定各小区作物籽粒产量和秸秆生物量及其有机碳和养分含量。

1.4 数据处理及计算

农田碳输入主要包括植物根系碳输入 (Croot)、根茬碳输入 (Cstubble)和肥料碳输入 (Cstrawormanure)。根据蔡岸冬等[8]报道的方法，估算各种来源碳输入量。根据Jobbagy等[10]报道的农作物根系生物量在0—100 cm剖面分配比例，估算各层土壤根系碳输入量。

土壤有机碳储量计算公式为：

SOCs=(Ci× BDi× Di)/10

式中，SOCs为第i层土壤有机碳储量 (t/hm2)；Ci为第i层土壤有机碳含量 (g/kg)；BDi为第i层土壤容重 (g/cm3)；Di为第i层土壤厚度 (cm)。

所有数据运用SigmaPlot 12.5软件进行作图，Canoco 5 软件进行冗余 (RDA) 分析。运用 SPSS 22.0软件，采用最小显著差数法 (LSD) 检验处理间差异显著性 (P＜0.05)。

2 结果与分析

2.1 长期不同秸秆还田方式下褐土作物产量和累积碳输入量

在春季和秋季 (图1)，4个处理玉米产量随种植时间延长总体均呈现相似的上升趋势。CK处理的玉米产量最低且年际间波动较大 (产量范围为2.21～10.91 t/hm2)。与CK相比，玉米产量的中位值 (1992—2013年) SC处理分别增加7.73%～14.69%，SM和CM处理分别增加17.12%和25.81%。

土壤碳输入量受秸秆还田方式和作物产量双重影响 (图1)。SM和SC处理秸秆碳输入量是CM处理有机肥碳输入量的3.7倍。与CK相比，SM、SC和CM处理0—100 cm土体累积根茬和根系的碳输入量平均分别增加16.12%、16.51%和31.48%，而总累积碳输入量平均分别增加1.4、1.4和1.0倍。

图 1 春季和秋季施肥下各秸秆还田处理玉米产量和土壤累积碳输入量 (1993—2012)Fig.1 Maize yield and cumulative carbon input under different straw incorporation methods and chemical fertilization seasons (1993–2012)

2.2 长期不同秸秆还田方式下褐土各土层有机碳含量和储量

化肥施用季节和秸秆还田方式显著影响剖面SOC含量 (表2)。与春季施肥相比，秋季施肥CM处理0—40 cm土层SOC含量显著提高，而SM和SC处理40—80 cm显著降低。不论春季施肥还是秋季施肥，与CK相比，CM、SM和SC处理表层(0—20 cm) SOC含量显著提高，而SM和SC处理40—60 cm和80—100 cm 显著降低。

表 2 春季和秋季施肥下各秸秆还田处理剖面有机碳含量Table 2 Soil organic carbon content in soil profile under different straw incorporation methods and chemical fertilization seasons

相对于CK处理，秸秆还田处理碳储量变化量在还田方式间存在显著差异，但在化肥施用季节间没有显著差异。在春季和秋季 (图2)，与CK相比，SM、SC和CM处理0—20 cm土层SOC储量平均分别增加2.32、5.42和12.60 t/hm2，CM处理的增加量显著高于SM和SC处理；而对于亚表层 (20—40 cm)，与CK相比，CM处理SOC储量略有增加 (0.56～0.99 t/hm2)，SM处理无明显变化，SC处理在春季和秋季施肥下分别降低2.98和4.62 t/hm2；与CK相比，SM、SC和CM 处理底层 (40—100 cm) SOC 储量平均分别降低3.98、6.99和3.76 t/hm2，3个处理间无显著差异。在0—100 cm土层，与CK相比，CM处理SOC储量平均增加9.62 t/hm2，增加量小于表层土壤 (12.60 t/hm2)；而SM和SC处理平均分别降低 1.81和5.36 t/hm2。

2.3 长期不同方式秸秆还田下褐土剖面养分含量

不同方式秸秆还田也显著影响褐土剖面养分含量 (表 3)。与 CK相比，春季和秋季施肥SM、SC和CM处理0—20 cm土层全氮、有效磷和速效钾含量均有不同程度增加，且CM处理显著高于SM和SC处理。对于亚表层，春季施肥CM处理全磷和有效磷含量显著高于SM处理，而秋季施肥CM处理有效磷含量显著高于其他处理；而对于底层，春季施肥CM处理全氮和全钾含量显著高于SC处理，秋季施肥SM处理全钾和有效磷含量显著高于CK处理。

表 3 春季和秋季施肥下各秸秆还田处理的土壤养分含量Table 3 Soil nutrient content under different straw incorporation methods and chemical fertilizer application seasons

2.4 长期不同秸秆还田方式下褐土有机碳储量变化的影响因素

以土壤累积碳输入量以及养分元素含量 (全氮、全磷、全钾、有效磷、速效钾) 作为影响SOC储量变化的因素进行冗余分析 (图3)。结果表明，碳输入和土壤养分含量对表层和下层SOC储量变化量的总体解释率分别为90.10%和31.80%。其中，影响表层SOC储量变化量的主要因子为有效磷、全氮含量和累积碳输入量，其解释率分别为80.10%、4.40%和2.80%；而影响下层SOC储量变化量的主要因子为全氮含量、碳输入量和全磷含量，其解释率分别为25.28%、3.67%和4.20%。

图 3 累积碳输入量和土壤养分含量对0—20 cm土层和20—100 cm土层土壤有机碳储量变化的解释率Fig.3 The explanation rate of cumulative organic carbon input and soil nutrient contents to the change of soil organic carbon stock in 0–20 cm and 20–100 cm soil

3 讨论

3.1 长期不同方式秸秆还田对褐土养分和作物产量的影响

本试验表明，与CK相比，秸秆还田均显著增加表层全氮、有效磷和速效钾含量，且均以秸秆过腹还田最优。这主要因为鲜牛粪和秸秆施用在表层土壤，其被微生物分解利用后养分元素逐渐释放，且鲜牛粪比秸秆含有更多碳和丰富的营养元素[17]。此外，所有处理玉米产量随种植时间延长总体均呈现升高趋势，这也间接表明单施化肥和秸秆还田配施化肥均能提升褐土地力，并改善作物生长环境[16]。此外，CK处理玉米产量在大多数年份最低且年际间波动较大，这可能是因为单施化肥降低土壤系统稳定性。同时，CM处理作物产量高于SM和SC处理，这表明CM处理的培肥效果比SM和SC处理好，本研究CM处理表层有机碳和养分含量更高也证明了这一点。然而，本研究还发现，SM处理增产效果比SC处理好，这是因为本研究地属半湿润半干旱偏旱地区，秸秆覆盖能调节土温并减少地表径流和水分的蒸发，能更好保墒和保肥[18]。

3.2 长期不同秸秆还田方式下褐土表层和底层有机碳储量变化差异及其影响因素分析

在农田土壤中，SOC含量是作物 (根茬和根系)及外源有机物 (秸秆和有机肥) 碳输入下，SOC累积与有机质分解导致SOC损失之间的净平衡[8]。本研究发现，与CK相比，SM、SC和CM处理均显著增加表层SOC储量，这与前人[14,19]研究结果基本一致。这是因为3个秸秆还田处理表层碳输入量是CK处理的2～3倍；而SOC饱和前，通常随碳输入增加而增加[8,20]。此外，CM处理表层SOC储量的增加量显著高于SM和SC处理，主要原因有：与秸秆相比，CM处理鲜牛粪主要是处于半分解状态有机质，其还田带入更多有效态氮，使其有较高碳转化效率[21]；同时，鲜牛粪促进团聚体形成和提升团聚体稳定性效果更好[20]。此外，SC处理表层SOC储量增加量显著高于SM处理，这表明等量碳输入下SC处理更有利于农田表层碳累积。Fan等[22]也发现，与覆盖还田相比，粉碎还田可增加耕层SOC含量。SC处理下秸秆与土壤颗粒充分接触，更易腐解转化并产生多糖等代谢物，促进土壤团聚体形成，进而增加碳吸附和固定[23]。SM处理下秸秆整株覆盖在干湿交替频繁且微生物活性高的地表，大部分秸秆碳被微生物分解，以CO2形式释放到大气[24]。

本研究还发现，与CK相比，春季和秋季施肥SM、SC和CM处理底层SOC储量均显著降低，与徐虎等[14]和盖霞普等[15]研究结果不一致，这表明长期秸秆还田措施下褐土底层碳的输入量并不能补偿有机质矿化造成的SOC损失量[8]。我们认为褐土底层碳耗竭可能有以下原因：1) 本研究区域属偏旱地区，在水分不足的条件下，致使根系向下延伸，但底层氮素含量低且供应不足，使微生物通过分解有机质提供作物生长所需氮素，引起底层SOC损失[25]。此外，本研究氮素对底层SOC变化解释率最高(25.28%)，也间接证明氮素供应不足可能是底层碳耗竭的主要原因。2) 秸秆和鲜牛粪还田带入大量可溶性养分，通过淋溶和促进根系生长增加底层新鲜碳源 (如可溶性有机质及根系分泌物)，使底层有机质因激发效应而矿化损失[26]。3) 秸秆和鲜牛粪还田促进作物根系生长，且深耕打破犁底层，改善土壤通气性，有利于作物根系下扎，导致土壤扰动及大团聚体崩解增加，加速底层有机质分解[27]。此外本研究还发现，SM、SC和CM处理均增加表层SOC储量，但大幅度降低底层SOC储量；对于0—100 cm土体，SM和SC处理SOC储量明显降低，而CM处理显著增加SOC储量 (9.62 t/hm2)，但增加量小于表层土壤 (12.60 t/hm2)。可见，仅考虑表层有机碳变化，可能会高估秸秆还田措施下SOC提升效果，尤其是SM和SC处理。因此，秸秆还田措施下褐土底层SOC变化应该被关注。

冗余分析发现，土壤累积碳输入和养分含量对表层SOC储量变化的总体解释率为90.10%，这表明养分和有机碳输入是调控表层SOC储量变化的主要因子。这是因为秸秆和腐熟鲜牛粪含有大量碳、氮及磷[4]。一方面，秸秆覆盖 (SM处理) 及粉碎和过腹还田 (SC和CM处理) 等措施主要作用在表层土壤，大幅度增加表层有机碳输入量，进而直接提升表层有机碳；另一方面这些措施降低了秸秆和有机肥转化过程中的微生物养分限制，更有利于表层有机碳增加[14,28-29]。然而，本研究还发现，碳输入和养分对下层SOC储量变化总体解释率仅为31.80%，这可能因为深层有机碳投入量和养分含量低且变化小，使其影响程度降低。此外，一些研究发现，下层SOC受粘粒含量、Ca2+、阳离子交换量、土壤孔隙度、铁和铝氧化物等影响可能更突出[25,30-31]。综上，长期秸秆还田措施下表层与下层SOC来源、含量和稳定性等存在差异，表层有机碳储量变化主要受养分含量和碳输入量调控；而在下层受这些因素影响相对较小，可能受到自身属性影响更大[25,28]。因此，在以后研究中，应进一步加强秸秆还田措施下底层SOC变化的驱动因素研究，为提升褐土农田有机碳固定和有机碳模型预测等提供基础理论。

4 结论

不论春季施肥还是秋季施肥，秸秆还田均提高褐土表层有机碳、全氮、有效磷、速效钾含量和作物产量，且均以秸秆过腹还田最优。秸秆还田均增加表层有机碳储量，大幅度降低底层有机碳储量；在0—100 cm土层，秸秆覆盖还田和粉碎还田降低有机碳储量，而秸秆过腹还田增加有机碳储量，但增加量小于表层土壤。可见，在评价秸秆还田措施对褐土有机碳影响时，仅考虑表层有机碳变化，可能会高估有机碳提升效果。综上，秸秆过腹还田是褐土农田培肥和固碳减排的重要措施，但秸秆还田措施下褐土底层有机碳变化也应该被关注。