坡耕地因其特殊的地形条件,在降雨和地表径流的作用下,成为土壤侵蚀、水土流失的重灾区,其对我国水土流失的贡献率高达28.3%
[1]。水土流失不仅使表层肥沃土壤的剥离与厚度降低,还会造成泥沙和洪涝灾害,导致土壤养分流失与水体污染。大量的营养物质如氮、磷、钾等关键土壤养分随泥沙和径流流失,导致该地区的土壤肥力下降,影响农作物产量和质量
[2]。
紫色土
[3]作为我国长江中上游的重要耕地资源,占土地总面积的18%
[4],因其多分布于丘陵坡地,土层较薄,多为初育土
[5],是我国水土流失、土壤侵蚀最为严重的地类之一
[6-7]。为减缓紫色土坡耕地土壤营养物质流失,国内学者进行了大量研究。研究发现,刘自飞等
[8]研究表明,降雨是造成地表径流最主要的驱动因子。孙理阳等
[9]研究表明,降雨是坡面产流产沙的关键驱动因素,且产流产沙量是影响氮磷流失的重要因素。史广宇等
[10]研究表明,施加生物炭缓解土壤氮素流失的机理主要为缓解铵态氮和硝态氮的流失、抑制N
2O的逸散、提高植物吸收氮素的能力。生物炭不仅拥有较强的铵态氮吸附能力,还具有提高土壤碳含量、pH值和土壤含水量等能力,从而降低因地表径流和水分入渗导致的铵态氮和硝态氮的流失量。刘懿枢等
[11]研究表明,施肥模式显著影响径流中氮磷流失特征,其中综合优化施肥模式对降低流失负荷的效果最优。与常规施肥相比,有机物料还田(秸秆)与矿质肥料(灰肥)协同施肥可有效减少稻田氮磷流失负荷。周晓航等
[12]研究表明,施加猪粪导致硝态氮流失量大于施加牛粪的流失,猪粪主要表现为硝态氮流失。目前关于紫色土坡耕地氮流失过程的研究主要集中在模拟降雨
[13-14]下土壤氮素随径流及泥沙的流失,虽然便于控制变量,但难以反映紫色土区自然降雨的复杂特征。亚热带湿润季风气候下的紫色土区自然降雨具有雨强变异性大、雨型多样、间隔时间不确定等特点,这些因素共同影响氮素流失过程。同时施肥处理单一,多以传统碳源为主,模拟降雨与自然降雨产生的径流亦有显著差异,实际估算氮素流失量方面有较大的误差。因此,在自然降雨条件下开展多碳源添加对氮素流失影响效应的研究,对于揭示紫色土坡耕地真实环境中的氮素流失阻控效应具有重要意义。
因此,本研究依托四川省农业科学院中江仓山综合试验示范基地,以紫色土坡耕地为研究对象,在自然降雨条件下,进行不同新型碳源(动物源、植物源)添加,对降雨产生径流泥沙以及氮素流失进行相关分析,探究不同处理下径流中不同形态氮素迁移特征,为提高紫色土坡地土壤肥力和减少土壤氮素流失提供理论依据。
1 材料与方法
1.1 试验区概况
研究地点位于四川省德阳市中江县仓山镇的生态试验基地,该地区经度为105°3′1.04″E,纬度为30°37′59.43″N。属于典型的紫色土坡耕地区域。气候类型为亚热带湿润季风气候,雨量充沛,年平均气温16.7 ℃,年降水量882 mm左右。试验区以丘陵和山地地形为主,主要种植油菜、玉米和小麦。土壤母质来源于侏罗纪遂宁组,发育形成紫色土,土层厚度约60 cm。由于长期耕作,土壤有机质含量较低,质地为砂壤,整体抗侵蚀能力较弱。
表1为该土壤的理化性质。
1.2 试验设计
用混凝土浇筑多个坡度为10°的土壤区域,土壤的表面面积为16 m
2(3.2 m×5 m),深度60 cm,在各个小区末尾设计一个汇流沟,使得降雨产生的地表径流汇集到一起,并在最后安放集流桶,为防止水流过多,在集水桶旁再添加一个分流桶,具体模型见
图1。采用夏季种玉米的种植方式,玉米使用“仲玉3号”,试验在玉米生长期进行,2023年5月17日播种,2023年9月20日收获。每个小区添加不同的碳源有机肥和土壤改良材料等N替代,一共设置6种处理,分别是化肥+秸秆生物炭(T
1)、化肥+猪粪有机肥料(T
2)、化肥+牛粪有机肥料(T
3)、化肥+高有机质含量生物有机肥(T
4)、化肥+复合微生物肥料(T
5)、施用常规化肥NPK处理作为对照组(T
6),每种处理进行3次重复,进行等N量施肥。玉米季氮肥分两次进行添加,第一次添加在玉米种植时作为基肥添加,添加量为总量的70%,第二次在7月上旬添加,在玉米种植前,把有机肥均匀地覆盖在土壤表面,并把土壤进行翻耕使有机肥进入0—20 cm土壤。每种处理的3个小区碳源添加量和管理都相同,每种处理的施肥量见
表2。
1.3 样品采集
在玉米播种前,收集基础土样。监测期间,玉米生育期无机氮质量分数变化比较大。把玉米各生育期期间的0—20 cm土层作为代表性样本,作为当月土壤无机氮含量背景值,玉米收获后再行取样。用梅花形多点法在0—20 cm土层进行采样,四分法混匀后装入自封袋,迅速带回实验室。把土样中的石块、根系、枯落物、土壤动物及秸秆等杂物清理干净,对土样进行自然风干,按分析需求研磨并通过20目和100目筛,用于测定土壤有机碳(SOC)含量和基本理化性质。每次降雨之后,对地表径流、壤中流、泥沙进行采样。采样时,观测各处理集水桶和分流桶的水位,计算各处理地表径流的产量。将各个处理集水桶里的水搅拌均匀,分别装入500 ml聚乙烯瓶,水样经过0.45 μm过滤,之后把过滤膜放到铝盒里,等至完全干燥,进行称重,计算流失泥沙量。
1.4 样品测定
SOC采用干烧法(multiN/C2100S)测定,土壤基本理化指标的测定方法见《土壤农化分析与环境监测》。水样测定指标包括硝态氮(NO3--N)、铵态氮(NH4+-N)和总氮(TN),总氮(TN)用碱性过硫酸钾消解-紫外分光光度法测定,铵态氮(NH4+-N)用纳氏试剂比色法测定,硝态氮(NO3--N)用紫外分光光度法测定。采用凯氏定氮法测试土壤的固态氮素含量(T)。土壤团聚体为水稳性团聚体,采用湿筛法测定,容重采用环刀法测定。
1.5 数据处理与分析
式中:F为氮湿沉降通量(mg/m2);c为降水事件中氮化合物的浓度(mg/L);p为降雨事件的降水量(mm)。
式中:Q为径流中氮元素迁移量(mg/m2);Ci 为径流中不同氮化合物的浓度(mg/L);q为小区单位面积产生的径流量(L/m2)。
本文采用Origin 2024进行图表绘制,采用SPSS 27版本进行影响氮素流失的相关性分析。用RStudio中的PLSPM软件包进行分析,采用偏最小二乘法路径模型(PLS-PM)来模拟降雨量、土壤理化性质、碳源添加、各种氮化合物的流失通量和土壤有机碳的关系,得出路径系数和R2。
2 结果与分析
2.1 自然降雨下氮湿沉降通量
如
图2所示,2023年总降雨量为983.1 mm,2024年总降雨量为835.2 mm。降雨主要集中在4—10月份,其中两年的7月份降雨量都达到顶峰,分别为386.3,317.1 mm。2023年10月为2023年降雨次数最多的月份,降雨14次,但总降雨量只有54.7 mm,只有2023年7月份降雨量的14.16%。2024年4月为2024年降雨次数最多月份,为17次,总降雨量为146.3 mm,只有2024年7月降雨量的46.13%。研究期间,小区共产生8次径流,分别为2023年的7月27日、8月7日、8月26日和9月11日以及2024年的7月8日、7月10日、7月14日和7月17日,具体降雨量见
图3A。
图3B是总氮湿沉降通量变化,2023年7月27日和9月11日产生的总氮湿沉降远高于其他6天,其中9月11日各处理的氮湿沉降达到峰值。T
1处理沉降通量达到最大值,为6.28 kg/hm
2。T
2,T
3,T
4,T
5与T
6处理沉降通量,分别为4.91,3.76,4.46,3.13,4.38 kg/hm
2。
2.2 降雨产流产沙的影响分析
图4是降雨产流产沙拟合图,可以看出降雨与产流与产沙均呈正相关,降雨量越大,产流与产沙越多。由
图4A可以看出,T
5处理的斜率(0.86)最低,表明T
5处理可以较好地延缓径流产生,但是T
5处理的截距较高,较低降雨量时延缓径流产生的效果并不如其他几种处理。T
2处理斜率(3.00)最大,对抑制径流产生的能力最弱。在较大降雨时,各处理对抑制径流产生的能力为T
5>T
6>T
1>T
4>T
3>T
2。
图4B是降雨量与泥沙流失量的关系,降雨量与泥沙流失量均呈正相关,降雨量越大,泥沙流失量越大。对各处理的斜率进行比较,各处理抑制产沙能力为T
5>T
6>T
3>T
1>T
2>T
4,T
5相对其他几种处理为最有效的减沙措施。这表明化肥+复合微生物肥料不仅可以减少径流产生,还可能通过增强土壤团聚体稳定性,降低径流对表土的冲刷力。在T
6和T
1,T
2,T
3,T
4处理比较中,常规化肥处理(T
6)的产沙量拟合线具有最高截距,表明其在中小降雨事件中具有最高的土壤侵蚀风险,这与其土壤结构恶化、表面物质松散的特征相符。
图5是产流产沙关系图,表示产流产沙之间的相互关系,各处理产流—产沙回归曲线均呈显著正相关(
p<0.001),产流增大,产沙量也增大。其中增长速率为T
1>T
4>T
5>T
3>T
6>T
2,其中T
4,T
5幂指数分别为1.157,0.979,产流产沙相对稳定,T
6处理的幂指数虽然不是最大,但是其对应的系数远大于其他5种处理,在本研究所设计的降雨量中,依旧是产沙量相对较多的处理。
2.3 不同形态氮在地表径流中的流失量
图6A是各种处理在产生地表径流时的铵态氮(NH
4+)的流失通量,在2023年7月27日除T
6处理外,其他5种处理铵态氮流失通量均达到峰值,2023年的铵态氮流失先减后增,分析原因是降雨量先减后增,泥沙流失量先减后增,铵态氮易吸附于土粒随之流出,2024年先增后减,4次试验数据降雨量相近,但7月14日的铵态氮流失显著变多,推测原因是经过多天的降雨,土壤水饱和,不易向下渗透,地表径流强度增大,泥沙流失也变多,导致铵态氮流失量变多,其中T
4处理流失总量最少,该处理可以更好地改良土壤环境,使土壤的抗侵蚀能力提高。
图6B是各种处理在产生地表径流时的硝态氮(NO
3-)的流失通量,通过与铵态氮对比,可以看出,氮素的流失大多是硝态氮的流失,2023年硝态氮流失先减后增,分析原因是硝态氮易随水流失,降雨量大,产生地表径流大,所流失的硝态氮越多,其中T
5处理流失硝态氮通量最少,T
2处理硝态氮流失量远大于其余5种处理。
图6C是各种处理在产生地表径流时的总氮(TN)的流失通量,从各处理的TN流失量总和来看,各处理占地表径流总流失量为T
5(11.97%)<T
4(12.97%)<T
1(15.12%)<T
3(15.32%)<T
6(17.54%)<T
2(27.08%)。
2.4 不同形态氮在壤中流中的流失量
图7A是壤中流中铵态氮(NH
4+-N)的流失通量,在2024年7月14日,T
4,T
5处理铵态氮流失高于其余4种处理,7月17日,T
4处理铵态氮流失远低于其余5种处理。
图7B是壤中流硝态氮(NO
3--N)的流失通量,与铵态氮流失相比,硝态氮是壤中流里氮素流失的主要形式,且流失通量远高于铵态氮。
图7C是壤中流总氮流失通量,2024年7月14日时,T
3处理流失量最小,T
6处理流失量最大,7月17日,T
5处理流失量最少,T
3与T
6处理流失通量最大。各处理壤中流氮素流失占比T
4(13.8%)<T
3(15.3%)<T
1(15.32%)<T
5(16.82%)<T
2(17.5%)<T
6(21.26%)。
2.5 碳源泥沙SOC与地表径流氮素流失的相关性
图8是添加碳源,泥沙SOC与地表径流不同形态氮素流失通量的相关性分析,6种处理中,除T
2处理外,其余5种处理中,硝态氮流失通量与总氮流失通量呈显著正相关(
p<0.01)。除T
6处理外,5种处理铵态氮流失通量与总氮流失通量呈显著正相关,其中T
2,T
4处理显著正相关(
p<0.05),其余3种处理显著正相关(
p<0.01)。T
4,T
5处理,碳源添加量与SOC显著正相关(
p<0.05),T
3,T
5处理,固态氮素含量(
T)与SOC显著正相关(
p<0.05),T
6处理碳源添加量与SOC显著正相关(
p<0.01)。
2.6 不同碳源添加对地表径流氮素流失的影响因素
PLS-PM分析进一步探究了不同碳源添加下降雨、土壤团聚体、容重、SOC、氮素流失之间的作用关系(
图9)。6种处理中,降雨量与氮流失通量均呈正向路径,其中T
3处理降雨量与氮流失通量路径最大(路径系数0.76,
p<0.05),这说明T
3处理(化肥+牛粪有机肥料)较其他处理对于保氮的能力较弱。碳源添加与SOC呈正向路径,其中T
2(路径系数0.577,
p<0.05)、T
4(路径系数0.570,
p<0.05)处理路径较大。降雨量与容重均呈负向路径且容重与土壤团聚体呈负向路径。T
4处理是化肥+高有机质含量生物有机肥,T
5处理是化肥+复合微生物肥料,这两种处理相比其他4种处理,降雨对土壤团聚体的路径系数(T
4路径系数0.93,
p<0.001;T
5路径系数0.98,
p<0.001)显著高于其他4种处理,但对大团聚体(>2 mm)形成起到了负面作用,降雨破坏了土壤大团聚体,进一步减少了土壤氮素随泥沙和径流的流失。
3 讨 论
3.1 自然降雨对土壤氮流失的驱动机制
氮流失与降雨冲刷、产流以及冲刷产生的泥沙量有关,氮湿沉降是土壤氮流失重要驱动力,是氮素地球化学循环的重要途径
[17-18],氮素通过沉降进入水体后,对生态环境系统的生产力及稳定性有着重要的影响。本研究期间,共有8次降雨产生地表径流,且都为强降雨(>30 mm)级别。通过对地表径流的各种氮元素流失通量进行分析,得出硝态氮是地表径流氮元素流失的主要部分,占地表径流总氮流失的52%,这与罗艺
[19]和任秀文
[20]等研究结果一致。硝态氮也是壤中流流失的主要部分,占壤中流总氮流失的80.67%。降雨量与产沙量呈正相关关系,在T
6和T
1,T
2,T
3,T
4处理比较中,常规化肥处理(T
6)的产沙量拟合线具有最高截距,表明其在中小降雨事件中即具有最高的土壤侵蚀风险,这与其土壤结构恶化、表面物质松散的特征相符。较多铵态氮随地表径流流失,在壤中流流失较少。这与熊子怡等
[21]研究结果一致。推测原因是,土壤颗粒会吸附铵态氮但几乎不吸附硝态氮,导致铵态氮大量滞留在土壤剖面的上中层,而硝态氮主要富集在下层,因此在本文中,各处理的地表径流铵态氮流失浓度均高于壤中流。
碳源添加下,各处理壤中流中的铵态氮流失通量差距不大,表明碳源对壤中流中铵态氮流失的阻控效果较弱,T
4处理壤中流氮素流失占所有处理总氮流失的13.8%。梁亚宇
[22]、刘平
[23]等研究表明降雨量是影响氮湿沉降通量的重要因素,呈正相关关系,但本研究2023年7月27日和2023年9月11日降雨量和氮湿沉降通量与上述研究结果有明显差异。7月27日降雨量205.2 mm,降雨历时达1 080 min,雨强11.4 mm/h;9月11日降雨量147.4 mm,降雨历时达1 800 min,雨强4.91 mm/h。9月11日降雨量总体小于7月27日,但是氮湿沉降通量总体大于27日。经过多日连续降雨,氮湿沉降通量有所下降,是由于雨水对空气中的氮素多次冲刷,导致地表径流中氮素的浓度远不如9月11日的氮素浓度。氮湿沉降对于土地氮循环是一个重要过程,我们需要进一步加强对大气氮沉降的监测和研究,优化田间施肥,并采取科学措施减少农田氮损失及人为氮排放,降低氮沉降带来的环境问题,为区域可持续发展提供科学依据。
3.2 土壤增碳对土壤固氮的影响机制
本研究通过添加多种碳源,减少地表径流中的氮素流失,碳源添加处理中,当降雨量大于50 mm(暴雨)时,T5处理化肥+复合微生物肥料对地表径流、泥沙量的产生具有较好的抑制效果。T1,T2,T3处理抑制产流产沙的能力弱于T6施加常规化肥处理,T4处理抑制产流能力弱于T6处理,抑制产沙能力高于T6处理,推测原因是,T6处理团聚体部分,大于2 mm部分与小于0.053 mm部分占比小于其他处理,导致地表径流量与泥沙产量较低,但6种处理中,T6处理产生的壤中流最多。除T2处理(化肥+猪粪有机肥料)外,均有效地减少地表径流氮素的流失,推测原因是猪粪分解较快,短期内不能较快地改善土壤团聚体结构,其次,猪粪的快速矿化使得氮素浓度较高,导致氮素流失量变大。本研究所产生降雨量大多数小于50 mm,碳源的添加在抑制产流产沙上均有成效,其中T5处理的地表径流总氮流失最少。T5处理氮元素流失量远少于只添加化肥的T6处理,说明添加复合微生物肥料,使得微生物产生多种可以促进土壤团粒结构形成的物质,使土壤变得疏松透气,保水保肥能力增强。且T5处理对产生径流和泥沙的抑制效果也高于其他5种处理。碳源添加下,各处理壤中流氮素流失占比T4(13.8%)<T1(15.29%)<T3(15.31%)<T5(16.8%)<T2(17.5%)<T6(21.3%)。5种处理均有效地抑制了壤中流中氮素的流失。
T
1处理在固氮上具有较好的稳定性,周期性较长
[24],提供了持续而非爆发性的碳源,避免了因碳源快速耗竭导致的固氮功能衰退。T
2与T
3处理在快速提升土壤整体肥力的同时,其对固氮功能的促进效率随时间波动的特征较为明显。T
4处理通过引用高质量碳源,与SOC正向路径较大,输入的稳定有机质能更高效地转化为持久性SOC。这一持续扩增的SOC库,为固氮微生物提供了长期稳定的碳源,确保了固氮功能的持续性,综合效率较高。T
5处理通过直接添加微生物,微生物主导了氮的矿化和固定,有效削减了氮素的流失
[25]。本研究通过PLS-PM模型和对影响氮元素流失因素的相关性进行分析,6种处理中,降雨量对土壤氮素流失均表现为稳定的正向直接影响,同时对土壤容重产生稳定的负向影响。这印证了降雨是氮素损失的主要环境驱动力。同时通过改变土壤团聚体结构和容重间接影响农田氮素流失。碳源添加通过改变土壤容重和土壤团聚体结构间接影响农田氮素流失。
4 结 论
降雨显著影响地表径流的产生,降雨与泥沙流失呈正相关关系,与产流产沙之间呈正相关关系。在降雨量大于50 mm时,T5(化肥+复合微生物肥料)处理可以较好地抑制地表径流和泥沙量的流失,T5,T4处理产生的氮素流失相对较少,占地表径流总流失的11.97%和12.97%,壤中流流失为13.8%和16.8%。
地表径流与壤中流中,硝态氮均是氮素的主要流失且铵态氮大多在地表径流中流失。各处理地表径流中氮素流失占比T5(11.97%)<T4(12.97%)<T1(15.12%)<T3(15.32%)<T6(17.54%)<T2(27.08%),各处理壤中流氮素流失占比为T4(13.8%)<T1(15.29%)<T3(15.31%)<T5(16.8%)<T2(17.5%)<T6(21.3%)。