玉米是我国产量最高的粮食作物,其种植面积不断扩大
[1]。目前,随着农业经济的发展,各种类型化肥的投入对农作物产量有着一定的促进作用,但与此同时,高强度的化肥投入也给农田土壤和环境带来一系列变化,其中氮素平衡问题尤为突出。氮素平衡在玉米生产中起着至关重要的作用,因此,研究如何有效调控玉米地的氮素平衡对玉米可持续生产尤为关键
[2]。大量研究已表明,土壤氮素以硝态氮为主要形式发生淋溶和累积
[3-4]。硝态氮淋失导致的氮素失衡与土壤水分、作物产量密切相关,这些因素的相互作用最终会显著影响农业生产。
在农业生产中,氮素失衡是影响作物生长和土壤环境的一个关键问题。地表覆盖技术作为一种有效的调控手段,能够在一定程度上缓解这一问题。例如,秸秆覆盖可增加土壤氮的固定,减少氮素淋溶,提高氮肥利用率
[5]。地膜覆盖可显著降低氮肥向深层土壤的淋溶、减少氨挥发等气态损失,延长氮肥在土壤中的有效性,从而提高氮肥利用率
[6]。目前,地膜覆盖与秸秆覆盖是干旱与半干旱地区农业生产中2个主要的地表覆盖技术,二者在玉米生育期均具有良好的蓄水保墒作用,还能在不同程度上提高玉米产量,如,与无覆盖相比,地膜覆盖和秸秆覆盖2年后春玉米的产量分别增加了24.9%和11.8%
[1,7]。覆膜在旱地农业中能够明显促进作物生长,提高氮素吸收利用效率,并对土壤中硝态氮的分布和累积产生影响
[8-9]。有学者认为,覆膜条件下,高温高湿环境会加速土壤有机氮的矿化速率,并在旱地条件下导致以硝态氮(NO
3--N)为主的矿质氮累积增加
[10-12]。根据张悦等
[13]的研究,采用全膜双垄沟栽培模式显著提升了春玉米产量和春玉米植株吸氮量及氮素回收利用率,同时也有效减少土壤剖面硝态氮的累积和淋洗。秸秆覆盖除能有效抑制土壤水分蒸发并调节地温外,还可腐解并为土壤提供额外的养分来源,增强土壤肥力,促进地上部氮素吸收以及作物产量形成
[14]。但是,在生产实践中,地膜和秸秆覆盖分别存在地力下降和苗期地温较低的不利现象
[15-17]。因此,尽管这2种技术各自具有明显优点,它们在长期应用中的综合效益仍需深入研究。
地膜和秸秆覆盖技术各有优缺点,在不同试验周期中展现出不同的效果。例如,1年的短期试验表明,将秸秆覆盖时间推迟到作物三叶期时,土壤中氮素含量显著增加,除40~60 cm土层外,其它土层的土壤硝态氮含量均随着秸秆覆盖量的增加而上升
[18]。3年的中期试验表明,在不同种植方式下,覆膜种植相比裸地平作更能有效减少施肥导致的硝态氮残留,尤其是在高肥处理下,全膜双垄沟播表现出显著的优势
[19]。若能继续延长试验周期,有助于揭示不同措施的累积效应,从而为农业生产的可持续性提供更加科学的指导,例如,通过对氮素动态和土壤健康的长期持续观察,更加准确地评估覆膜与秸秆覆盖的长期效益,以实现作物可持续生产。Wang等
[20]统计2000-2022年的不同生态区域的无土覆盖、地膜覆盖和秸秆覆盖处理的田间试验报道后,发现地膜覆盖在单作玉米系统中产量增幅更高,为30.3%;秸秆覆盖则在双作种植系统中表现较好,使得双作玉米的产量增加了21.7%;该研究认为地膜覆盖和秸秆覆盖处理联合施用对作物可持续生产具有潜在的促进作用。根据陈桂平等
[21]的研究,发现秸秆地膜二元覆盖对于促进玉米生长发育而稳产高产具有明显的增效作用,但试验周期较短,还不清楚2种覆盖方式相结合的长期影响。因此,本研究综合了秸秆覆盖在增加有机物料输入方面的优势,以及地膜覆盖在蓄水保墒方面的优势,深入探讨这2种覆盖方式及其结合应用的效果。通过在山西省旱塬地上种植春玉米并长期进行不同覆盖方式影响的定点观测,研究种植第15年不同覆盖方式对玉米籽粒产量、氮肥利用效率的影响,解析不同覆盖方式下土壤硝态氮的分布和累积特征;以及分析种植第0、7、10、15年时不同覆盖方式下各处理硝态氮在土壤剖面中的年际移动变化,旨在探寻适宜当地玉米生产和生态环境的覆盖类型,进而为该区玉米生产中提高氮肥利用效率、有效降低硝态氮在土壤中的淋失提供理论依据。
1 材料和方法
1.1 试验地概况
长期定位试验于2008-2023年在山西省寿阳县宗艾村北坪旱塬区(37°45′21″N,113°54′2″E)进行,海拔为1130 m,年均气温7.4 ℃,年均降水量501.1 mm,70%左右的降水集中在6-9月,属暖温带半湿润偏旱区。试验地表面平坦,供试土壤为褐土性土,表层土壤的理化性质见
表1[22-23],土层深厚,作物为春玉米,种植制度为一年一季。
1.2 试验设计
长期定位试验设无覆盖(CK)、秸秆覆盖(J)、普通地膜覆盖(P)、宽幅地膜覆盖(S)、宽幅地膜与行间秸秆覆盖(SJ)和普通地膜与行间秸秆覆盖(PJ)6个处理(
表2)
[24],2008-2023年种植春玉米,每年种植并收获1次,每个处理设3次重复,随机区组排列,小区面积为6 m×8 m,以不施肥为空白。地膜覆盖时间为春玉米全生育期(4-9月),每年覆盖1次。所有处理均施用晨雨调控肥(N∶P
2O
5∶K
2O=19∶11∶10)作基肥,施用量为900 kg/hm
2。全生育期不进行灌溉,以自然降水为水分来源。
1.3 样品采集与指标测定
1.3.1 植物与土壤样品采集
于种植第15年,即2023年10月,在每个小区内选取10株生长高度相对一致、籽粒成熟度相似、未发生明显病虫害的玉米植株,收集植株地上部分;收获玉米后用五点采样法采集土壤样品,采样深度为0~500 cm,每20 cm为1层,每小区收集1个土样,重复3次。
1.3.2 玉米产量与全氮含量测定
于种植第15年,即2023年10月,将植株地上部分划分为茎、叶、穗轴和籽粒,经过65 ℃、24 h干燥后进行称重,最后计算其生物产量。粉碎植株样品后,采用H
2SO
4-H
2O
2法进行消煮处理,并通过凯氏定氮仪测定植株全氮含量
[24]。各指标重复3次。产量公式如下:
玉米产量(kg/hm2)=处理小区成穗数(穗/48 m2)×平均单穗风干籽粒重(kg/穗)×10000 m2/hm2
1.3.3 土壤硝态氮含量的测定
于种植第0、7、10、15年收集土壤样品,样品风干后使用1 mol/L的KCl溶液,以水土比10∶1对土壤样品进行浸提,恒温震荡1小时后进行过滤,通过流动分析仪测定滤液中的硝态氮浓度,以计算种植第0、7、10、15年的土壤硝态氮含量及种植第15年的氮素利用相关指标
[5,10],各指标重复3次。常见的氮素利用相关指标参照解文艳
[24]和张秀芝等
[25]的方法进行计算,公式如下:
植株吸氮量=植株干重×植株全氮含量
氮素回收利用率=(施氮肥处理地上部吸氮量-不施氮肥处理地上部吸氮量)/氮投入量×100%
土壤硝态氮累积量(kg/hm2)=土层厚度×土壤容重×土壤硝态氮浓度/10
氮肥贡献率=(施氮肥处理籽粒产量-不施氮肥处理籽粒产量)/施氮肥处理籽粒产量×100%
土壤氮素依存率=不施氮肥处理作物吸氮量/施氮肥处理作物吸氮量×100%
氮肥偏生产力=施氮肥处理籽粒产量/施氮量
1.4 数据处理
采用Excel 2019软件计算数据的平均值和标准差;采用SPSS 2019软件进行差异显著性分析,使用单因素方差分析(ANOVA)中的最小显著性差异(Duncan)法;采用Origin 2019软件进行作图。
2 结果与分析
2.1 不同覆盖方式对玉米产量及氮素吸收利用的影响
由
表3可见,与无覆盖(CK)处理相比,不同覆盖方式均能提高玉米产量,但不同覆盖方式下玉米产量存在明显差异,排序为SJ>S>P>PJ>J>CK。与CK处理相比,5种覆盖方式的产量增加了5.98%~19.86%,SJ处理增幅最大,显著高于CK处理。不同覆盖方式处理显著提高了玉米地上部的总吸氮量,各处理的总吸氮量从大到小的排序与玉米产量一致;与CK处理相比,5种覆盖方式的吸氮量增加了5.53%~31.17%,SJ处理增幅最大,显著高于CK处理。在同等施氮量的条件下,氮素盈亏量排序为CK>J>P>PJ>S>SJ。与CK处理相比,5种覆盖方式的氮素盈亏量减少了12.37%~69.80%,SJ处理减幅最大,显著低于CK处理。各处理氮素回收利用率从大到小的排序与产量和吸氮量一致,SJ处理的氮素回收利用率显著高于其它处理。氮素盈亏的变化趋势与氮素回收利用率相反,相较于CK处理,5种覆盖方式的氮素回收利用率提高了3.82%~21.55%,SJ处理增幅最大,显著高于CK。
2.2 不同覆盖方式对氮肥贡献率、土壤氮素依存率和氮肥偏生产力的影响
氮肥贡献率可反映当年氮素吸收利用对植株干重影响的贡献程度。从
表4中可见,SJ处理氮肥贡献率最高,达到51.96%,显著高于其它覆盖处理以及CK;CK氮肥贡献率最低,为42.42%,与CK相比,5种覆盖方式的氮肥贡献率提高了3.25%~9.54%,各处理氮肥贡献率的顺序为SJ>S>PJ>P>J>CK。土壤氮素依存率可反映春玉米对土壤中氮素资源的依赖程度,体现土壤供氮的能力。CK处理的土壤氮素依存率达到了29.59%,显著高于其它覆盖处理;与CK相比,5种覆盖方式的土壤氮素依存率降低了1.55%~7.03%,排序为CK>J>P>PJ>SJ>S。从氮肥偏生产力来看,SJ处理最高,达到39.98%,显著高于其它覆盖处理以及CK;CK处理氮肥偏生产力最低,仅为33.36%,与CK相比,5种覆盖方式的氮肥偏生产力增加了5.98%~19.86%,排序为SJ>S>PJ>P>J>CK。通过氮肥贡献率、土壤氮素依存率和氮肥偏生产力的分析结果,发现SJ处理效果最好,能增强肥效,提高氮素的利用率与生产力。
2.3 不同覆盖方式下土壤剖面硝态氮分布特征
长期定位试验条件下,土壤硝态氮在土体中向下迁移(
图1),且在不同覆盖方式下,0~500 cm土层中硝态氮的分布存在显著差异(
P<0.05),各处理硝态氮含量呈先下降后上升再下降的变化趋势。各处理在0~100 cm土层的变化趋势是先上升后下降,在100~300 cm和300~500 cm土层的变化趋势均是先上升后下降。各处理在200 cm以下土层含水量总体维持在较高水平,且不同覆盖处理的含水量显著高于CK。不同覆盖方式下硝态氮在0~100 cm土层变化趋势相似,均先下降后上升;但在100~500 cm土层,各处理硝态氮含量分布的变化幅度较大,处理间差异显著(
P<0.05)。在长期定位试验中,各处理土壤中硝态氮出现了明显的积累现象。CK处理的硝态氮含量出现2个峰值,分别在100~120 cm和220~240 cm土层。相比之下,不同覆盖方式处理的硝态氮峰值较CK处理有所降低,且,向250~400 cm土层迁移。P、S、SJ和PJ处理在100~500 cm土层的硝态氮含量分布出现多个峰值,在300~340 cm土层硝态氮含量最高。除CK处理外,各处理0~200 cm土层的硝态氮含量分布依次为SJ>S>PJ>J>P,而在200 cm以下土层,硝态氮含量分布依次为J>S>P>SJ>PJ。
2.4 不同覆盖方式对土壤硝态氮累积量的影响
经过长期施肥后,不同覆盖方式处理的土壤硝态氮累积量如
表5所示,CK处理0~500 cm土层硝态氮累积最多,高达2 955.28 kg/hm
2,表现出较高的淋溶风险。各覆盖处理硝态氮积累主要出现在200~500 cm土层,其次为80~200 cm土层,在0~80 cm土层各处理硝态氮的累积量均极低。CK、P、S、SJ和PJ处理0~200 cm土层硝态氮累积量在0~500 cm土层中占比不超过30%;CK、J、P、PJ、S、SJ处理200~500 cm土层硝态氮累积量分别占0~500 cm土层总累积量的70%、90%、93%、78%、82%、77%。长期的定位试验结果表明,6个处理的硝态氮累积量呈现明显增加趋势,种植15年后各处理总累积量有77%~93%集中在200 cm以下的土层,各处理在200 cm以下的比例大都超过70%,可视为氮素的淋失。SJ在200 cm下的累积量占比小于其它覆盖处理,且总量小于无覆盖处理。因此,SJ处理能够减少深层土壤中硝态氮的积累,从而降低硝态氮的淋失风险。
2.5 不同覆盖方式下土壤剖面硝态氮的阶段性迁移
图2以种植第0、7、10、15年500 cm土层中的硝态氮含量变化为例,展示了6个处理下土壤中硝态氮含量的年际变化,可用于分析土壤硝态氮的阶段性迁移。CK的硝态氮含量在种植第7年,峰值出现在40~60 cm和120~140 cm土层;种植第10年,峰值出现在120~140 cm土层;种植第15年,峰值出现在100~120 cm和220~240 cm土层,即种植第7年硝态氮淋洗到140 cm以下的土层,至种植第15年时已淋洗到240 cm以下的土层。与CK相似,5个覆盖方式下硝态氮含量在种植第7年的峰值为40~180 cm土层,但在种植第15年时已下移至120~360 cm土层。种植第15年与第7年相比,各处理的硝态氮峰值均明显下移,8年间下移范围为80~180 cm,且呈现向500 cm土层以下淋洗的趋势。水分是硝态氮迁移的主要载体,硝态氮会随水分在土壤中迁移,尤其是在降雨或灌溉后,硝态氮可能被淋溶到深层土壤中。可见,与CK相比,SJ处理氮素盈余量减幅最大,显著低于CK处理;在年际变化中,该处理总迁移量显著小于CK以及其它覆盖处理(
P<0.05)。
此外,不同覆盖方式下土壤剖面硝态氮的阶段性迁移和降雨量密切相关。例如,沈新磊等
[26]研究表明,试验区降水少时,硝态氮主要分布在0~80 cm土层,降水多时,硝态氮进入下层土壤比例增加,这与本研究结果相似。在15年的试验期间,通过计算降雨量得出,正常年有9年(2008、2009、2010、2011、2012、2015、2017、2021、2023年),偏旱年有2年(2014、2016年),大旱年有1年(2019年);丰水年有4年(2013、2018、2020、2022年)。不同覆盖方式处理土壤剖面硝态氮峰值随种植时间延长明显下移,且均呈现向500 cm土层以下淋洗的趋势,可能也与种植期间正常与丰水年份较多有关。与CK相比,SJ处理氮素盈余量降幅最大,显著低于CK处理;在年际变化中,其总迁移量显著小于无覆盖处理以及其它覆盖处理(
P<0.05),氮素迁移量小有利于春玉米产量增加。
3 讨论
在干旱半干旱地区,覆膜和施肥处理可明显提高玉米产量,还能提高氮肥利用率和作物对土壤氮素的吸收量
[27],本研究表明宽幅地膜与秸秆覆盖处理和无覆盖处理以及其它覆盖处理相比,玉米产量增加了35.86~1 055.76 kg/hm
2;玉米植株吸氮量增加了1.90~40.42 kg/hm
2;氮素回收利用率提高了0.84%~17.73%。这可能是由于宽幅地膜相比普通尺寸地膜具有更好的蓄水保墒作用,同时秸秆覆盖增加了有机物料的输入,改善了土壤肥力。但经过长期的耕作生产,我国北方干旱和半干旱地区的土壤中有大量的硝态氮累积,这种累积不仅会影响农业生产的可持续性,还导致潜在的环境风险增高。氮素平衡作为一个重要指标,不仅关系到农业生产发展的可持续性,还直接影响到环境效益
[28]。长期试验可通过揭示不同覆盖方式的累积效应及其对氮素动态和土壤健康的长期影响,为评估农业管理措施的可持续效益提供科学依据,从而实现作物可持续生产。本研究结果表明,种植第15年6个处理的春玉米产量和氮素回收利用率的排序相同,氮素盈亏量和硝态氮含量分布、累积量的排序相反,这也与多项研究结果相似
[29-31],即硝态氮的分布和累积与春玉米产量和氮素盈亏量大小密切相关。
不同深度土层硝态氮的含量反映了硝态氮在土体中的积累与迁移趋势,可用于评价土壤硝态氮的淋失风险
[25]。郭星宇等
[32]研究发现随着年份的增加,土壤硝态氮的累积量逐年增加且向更深的土层迁移,并在土层加深的过程中呈现下降趋势,这与本试验结果相似,各处理硝态氮含量在种植第7年的峰值为40~180 cm土层,但在种植第15年时已下移至120~360 cm土层。种植第15年与第7年相比,各处理的硝态氮峰值均明显下移,8年间下移范围为80~180 cm,且呈现向500 cm土层以下淋洗的趋势。这种变化可能导致土壤氮素失衡,进而影响作物对其它养分的吸收,最终对作物产量产生负面影响。覆盖可抑制氮淋失,但不同覆盖方式对氮素淋失的长期影响尚少有报道。本研究表明,施用化肥用量一致的条件下,种植第15年不同覆盖方式处理0~500 cm土体土壤硝态氮累积量低于CK处理,各处理的总累积量有77%~93%集中在200~500 cm土层,且各处理在200 cm以下的累积量占比大都超过70%,可视为氮素的淋失。而SJ处理在200~500 cm的累积量占比小于其它覆盖处理,且总量小于无覆盖处理,意味着SJ处理能够减少深层土壤中硝态氮的积累,从而降低硝态氮的淋失风险。根据胡锦昇等
[33]研究,发现在旱塬区采用免耕基础上的秸秆结合地膜覆盖能有效地调节土壤的水分运动,可以减少硝态氮的淋溶累积,从而改善农田的水肥状况,进而增加作物产量。本研究将秸秆与地膜覆盖后秸秆结合地膜覆盖能够有效提高土壤保水保肥能力,减少和氮素流失,同时改善土壤结构,促进作物生长和增产,也与此结论相似。
不仅是长期施氮可导致硝态氮向下层土壤迁移,水分变化也对土壤氮素的移动和分布有着明显影响。土壤中的硝态氮分布特征的年际变化与水分变化密切相关,降雨会增加土壤水分,而土壤水分越多,硝态氮淋失越多
[25,34]。不同覆盖处理的0~100 cm土层土壤含水量比CK处理的含水量高,即相同降雨量条件下,CK处理的土壤水分下渗量大于不同覆盖方式处理。但不同覆盖方式处理下,种植第15年与第10年相比,硝态氮峰值向下移动的幅度较大,这可能与种植期间的降水情况有关。在旱作农业区,降雨对氮素淋溶起着决定性的作用,而且降雨年际间的变化会造成土壤中硝态氮淋失的变化
[35]。在降雨量或者灌溉量较高时,土壤中的硝态氮易于向下淋洗;而在干旱条件下,由于深层土壤水分向上移动,较深土层中的硝态氮也随之上移
[36]。本研究中,与种植第7年、第10年的峰值位置相比,种植第15年不同覆盖方式处理硝态氮峰值向下移动的幅度较大。这可能是由于种植第10年至种植第15年间有多个丰水年,造成土壤水分的增加,从而导致硝态氮的淋失增加,向下层土壤迁移。同时,由于地表覆盖具有很好的保水效果,可以改善深层土壤水分状况,减少不同土层水分变化,使得土壤水分含量较高且相对稳定,未被利用的氮素以硝态氮的形式随水分下渗淋移到土壤深层并累积
[37]。此外,地下水位下降或区域性缺水也会导致土壤的保水能力下降,进而导致硝态氮不断向下淋溶
[33]。未来可继续深入研究长期覆盖对氮素淋失的影响机理,为可持续农业发展提供科学支持。
4 结论
本研究依托长期定位试验,于种植第15年对无覆盖(CK)、秸秆覆盖(J)、普通地膜覆盖(P)、宽幅地膜覆盖(S)、宽幅地膜与行间秸秆覆盖(SJ)和普通地膜与行间秸秆覆盖(PJ)6个处理下的氮素回收利用率、氮肥贡献率和土壤氮素依存率及氮肥偏生产力这4个指标进行综合分析,发现与CK相比,不同覆盖方式处理均显著提高了玉米产量和氮素吸收效率,SJ处理在玉米产量、氮素回收利用率、氮肥贡献率和土壤氮素依存率以及氮肥偏生产力方面均表现最佳。6个处理下有77%~93%的硝态氮在200~500 cm土层累积,SJ处理在200 cm以下土层的硝态氮累积量占比低于其它覆盖处理,且总量低于CK。在年际变化方面,不同覆盖方式处理土壤剖面硝态氮峰值随种植时间延长明显下移,且均呈现向500 cm土层以下淋洗的趋势。此外,土壤剖面硝态氮的阶段性迁移变化(年际变化)结合降水年型分析可发现,经过多个丰水年份的累积,与种植第7年、第10年的硝态氮峰值位置相比,种植第15年不同覆盖方式处理硝态氮峰值向下移动的幅度较大。综合分析得出,施用化肥用量一致的条件下,与其它处理相比,宽幅地膜与行间秸秆覆盖相结合,显著提高了氮素回收利用效率和氮肥贡献率,同时有效减少了土壤氮素依存率,并降低了土壤硝态氮的淋失风险,从而促进了玉米产量的提升。
国家重点研发计划项目(2023YFD1900503-02)
科技基础资源调查专项(2021FY100501)