中国生态农业学报(中英文)  2021, Vol. 29 Issue (1): 154-162  DOI: 10.13930/j.cnki.cjea.200707
0

引用本文 

刘静, 朱鑫宇, 李顺江, 康凌云, 马茂亭, 杜连凤. 地下水位波动对不同施氮量农田土壤硝态氮运移影响[J]. 中国生态农业学报(中英文), 2021, 29(1): 154-162. DOI: 10.13930/j.cnki.cjea.200707
LIU J, ZHU X Y, LI S J, KANG L Y, MA M T, DU L F. Effects of groundwater fluctuation on nitrate nitrogen transport after nitrogen application in cropland soil[J]. Chinese Journal of Eco-Agriculture, 2021, 29(1): 154-162. DOI: 10.13930/j.cnki.cjea.200707

基金项目

国家“十三五”重点研发计划项目(2016YFD0800104)、国家自然科学基金项目(41907078)和北京市自然科学基金项目(8202020)资助

通信作者

杜连凤, 研究方向为农田面源污染防控。E-mail:dulianfengyzs@163.com

作者简介

刘静, 研究方向为农田面源污染防控。E-mail:liujing9712@163.com

文章历史

收稿日期:2020-08-31
接受日期:2020-09-29
地下水位波动对不同施氮量农田土壤硝态氮运移影响*
刘静, 朱鑫宇, 李顺江, 康凌云, 马茂亭, 杜连凤     
北京市农林科学院植物营养与资源研究所 北京 100097
摘要:明确地下水位波动对农田土壤剖面和地下水NO3--N运移的影响,可为减少土壤氮素淋失、降低地下水硝酸盐污染风险提供依据。本研究采用大型土柱温室种植甘蓝,研究2种水位波动(水位不变、水位每隔10 d波动20 cm)和3种施氮量[0 kg(N)·hm-2、225 kg(N)·hm-2、450 kg(N)·hm-2]对土壤含水量、土壤溶液NO3--N浓度、地下水NO3--N浓度和作物产量的影响。结果表明,水位波动和施氮肥对NO3--N运移的影响与土壤剖面深度有关。0~20 cm包气带土壤NO3--N含量受施氮量影响,过量施氮肥[450 kg(N)·hm-2]导致该剖面NO3--N累积。20~60 cm水位波动带土壤NO3--N含量受施氮量和水位波动的共同作用:施氮量增加提高NO3--N含量;水位波动降低剖面土壤NO3--N含量,水位上升和下降均促进土壤NO3--N随着水流运动向下层迁移;剖面土壤硝态氮含量高,增加NO3--N进入地下水的风险。60~80 cm淹水区剖面土壤NO3--N含量较低。作物产量受水位波动影响不显著。在地下水位埋深较浅的农业区进行氮素污染防控时,不可忽视水位波动对NO3--N运移的影响。
关键词水位波动    施氮量    土壤溶液    硝态氮    运移    地下水    产量    
Effects of groundwater fluctuation on nitrate nitrogen transport after nitrogen application in cropland soil*
LIU Jing, ZHU Xinyu, LI Shunjiang, KANG Lingyun, MA Maoting, DU Lianfeng     
Institute of Plant Nutrition and Resources, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China
Abstract: Understanding the effects of groundwater fluctuations on nitrate nitrogen (NO3--N) transport in the soil vertical profile is important for reducing nitrogen (N) leaching and nitrate pollution in croplands with shallow groundwater. This study investigated effects of groundwater fluctuations and nitrogen application rate in a greenhouse cabbage with large soil column experimental devices. Two groundwater fluctuation levels (W0, stable; W1, 20 cm fluctuation per 10 days) and three levels of nitrogen application (N0, 0 kg(N)·hm-2; N1, 225 kg(N)·hm-2; N2, 450 kg(N)·hm-2) were tested to determine water contents and NO3--N concentration changes at different soil depths, NO3--N concentration in groundwater, and crop yield. The results showed that fluctuating groundwater affected NO3--N transport depending on the soil depth and nitrogen application rate. At 0–20 cm depth (part of the vadose zone), excessive N application led to NO3--N accumulation in soil, but there was no correlation with groundwater fluctuation. At 20–60 cm depth (groundwater fluctuation zone), increased N fertilizer increased the soil NO3--N, and groundwater level changes promoted NO3--N migration into deeper soil, increasing the NO3--N groundwater pollution risk (especially in soil with already high NO3--N concentrations). At 60–80 cm depth (the flooding area), there was less NO3--N in soil, mainly due to denitrification. Crop yield did not significantly correlate with groundwater level changes. Groundwater fluctuations affected = NO3--N transport and should be considered in agricultural areas with shallow groundwater levels to prevent and control nitrogen pollution.
Keywords: Groundwater fluctuation    Nitrogen application rate    Soil solution    Nitrate nitrogen    Transport    Groundwater    Yield    

氮素对促进作物生长, 保证作物高产、稳产具有重要作用[1]。中国是世界上最大的化肥生产和消费国, 每年氮肥施用量为3360万t, 占全球的33%, 导致农田氮盈余量达175 kg·hm-2[2]。旱地土壤中强烈的硝化作用使氮素主要以NO3--N形式存在[3], NO3--N带负电荷, 不易被土壤吸附, 可以通过质流方式随土壤水分迁移进入地下水, 对人类健康产生危害[4-5]。NO3--N在农田土壤中的迁移及对地下水污染问题一直是研究热点。许多文献已表明土壤含水量、pH、温度、土壤质地以及施肥、气候、地理条件、耕作方式、作物管理等因素影响NO3--N在农田土壤剖面的迁移、再分布过程[6-7], 水分是NO3--N移动的载体、淋洗的驱动因素[8]。然而, 在地下水埋深较浅的潮土地区或背河洼地, 大气降水与农田灌溉用水的垂向补给或河流的侧向补给会造成地下水水位上升或下降。地下水的波动导致溶解氧、含水量等土-水环境指标改变[9-10], 这将对土壤及地下水中NO3--N的迁移转化产生影响。因此, 为防控该区域农田面源污染, 亟需认识地下水位波动下农田土壤NO3--N的迁移转化规律。

近几年, 水位波动对土壤环境或地下水中NO3--N的影响已有少许报道[11-13]。李翔等[12]利用填装中砂的柱子研究得知水位波动有利于氮素垂向迁移, NO3--N浓度随着水位上升或下降过程发生降低或升高。颜晓元等[13]研究发现长三角地区农田不存在明显的NO3--N累积, 主要原因是1~4 m地下水中发生的反硝化作用有效去除了NO3--N。华北地区夏季降雨量占全年的60%~80%, 是该地区露地农田土壤氮淋溶的防控关键期[14]。夏季集中降雨导致该地区近岸农田或地下水位埋深较浅农田出现地下水波动现象, 目前对该区域农田由于地下水位波动导致NO3--N在土壤和地下水运移及淋失风险的研究尚少见报道。此外, 以往的研究大多采用填装沙子的土柱或者野外采集原状土, 不能反映不同施氮量农田种植状况, 亟待对现有试验装置进行改进。为此, 本研究为同时开展种植农作物、调节地下水位、获取剖面土壤淋溶液及地下水样品等提供试验装置, 对不同施氮量和水位波动下农田土壤体积含水率、土壤溶液NO3--N浓度及地下水NO3--N浓度变化和作物产量进行研究, 以便明确水位波动和不同施氮量下农田土壤NO3--N运移规律及淋失风险, 对控制土壤残留氮素淋失, 降低地下水污染风险、制定农业面源污染防控措施具有实际意义, 也为该地区农田高效优质和可持续发展提供理论基础。

1 材料与方法 1.1 试验装置

为研究地下水位波动对土壤氮素运移的影响, 开发研制了一套模拟水位波动的土柱装置(图 1)。装置由外筒、内筒、底盘、进水阀、排水阀等组成; 外筒和内筒之间为可控升降水槽, 根据“U”型管原理, 通过进水阀和排水阀控制水位的上升和下降, 模拟地下水消长。内筒直径50 cm, 高120 cm, 填装农田土壤, 填充容量0.2 m3。土柱由下至上依次填充10 cm粗陶粒、5 cm细陶粒、100 cm农田土壤。在土柱0~20 cm、20~40 cm、40~60 cm和60~80 cm均匀布设土壤溶液提取器和水分传感器(型号MCS203, 数据采集器型号MC502L)。利用土壤溶液提取器采集土壤淋溶水, 土柱出水阀开启流出水作为地下水样。

图 1 水位波动模拟试验土柱装置图 Fig. 1 Schematic diagram of soil column device for simulation test of groundwater level fluctuation
1.2 试验设计

2020年4—6月于北京市农林科学院温室开展试验, 试验期间温度为20~35 ℃。供试作物为甘蓝(Brassica pekinensis), 品种为‘日本铁头甘蓝’, 种植密度为35 000棵·hm-2, 4月5日定植, 6月10日收获。试验设置2个因素:水位波动和施氮量。水位波动2个处理: W0, 水位始终保持在80 cm; W1, 水位每隔10 d变动20 cm, 试验期内W1处理水位依次为80 cm、100 cm、80 cm、60 cm、40 cm和60 cm, 水位高低通过装置的进水阀和出水阀控制, 以模拟地下水下降-上升-下降的过程。施氮量设3个: 0 kg(N)·hm-2(N0)、225 kg(N)·hm-2(N1)和450 kg(N)·hm-2(N2)。共6个处理, 每个处理重复2次, 随机区组排列。各处理所施氮肥为尿素; 磷钾养分总量一致, 磷肥90 kg(P2O5)·hm-2、钾肥180 kg(K2O)·hm-2。N1和N2处理组施用鸡粪有机肥(N 45 kg·hm-2、P2O5 92.33 kg·hm-2、K2O 37.92 kg·hm-2)。有机肥、磷肥和钾肥全部基施, 氮肥40%做基肥, 60%做追肥。追肥2次, 追施比例均为30%。试验监测共60 d, 每10 d作为一个监测期。在每个监测期的第1 d完成水位波动, 稳定9 d后, 第10 d测定剖面水分含量、收集剖面淋溶水样和地下水样等。第21 d和第41 d追施氮肥。灌溉、病虫害防治及日常管理按照当地农民习惯进行。

1.3 样品采集与测定

土壤淋溶水和地下水的NO3--N浓度用连续流动分析仪(A3)分析, 土壤剖面水分含量采用水分传感器实时监测。本研究以剖面土壤溶液NO3--N含量变化表征不同施氮量和水位波动下NO3--N沿农田土壤剖面运移状况; 以地下水NO3--N浓度作为NO3--N沿土壤剖面进入地下水的状况。甘蓝成熟后收获测产, 根据种植密度计算各处理产量(kg·hm-2)。

1.4 数据分析

运用Microsoft Excel 2010进行数据处理及绘图。运用R统计软件(3.0.3)进行相关分析, 明确水位波动和施氮量对剖面土壤溶液NO3--N浓度、地下水NO3--N浓度和作物产量的影响。

2 结果与分析 2.1 水位波动和施氮量对剖面土壤体积含水率的影响

图 2可知, 随着土壤深度增加, 剖面土壤体积含水率显著递增。60 d与10 d相比, W0处理组0~20 cm、20~40 cm、40~60 cm和60~80 cm剖面平均体积含水率变化幅度分别为-3.60%、0.01%、0.55%和0.75%;除0~20 cm, 其余剖面土壤体积含水率变化幅度较小。W1处理组0~20 cm、20~40 cm、40~60 cm和60~80 cm剖面平均体积含水率第60 d比第10 d增加-1.71%、7.83%、4.01%和1.73%, W1变化幅度大于W0处理组; 20~40 cm和40~60 cm剖面平均体积含水率变化最大, 第20 d与第10 d相比, 即水位从80 cm下降至100 cm, 两剖面平均土壤体积含水率分别下降5.16%和3.46%;当水位从100 cm (第20 d)逐步上升至40 cm (第50 d)时, 20~40 cm和40~60 cm剖面土壤体积含水率持续上升17.14%、7.29%;水位从40 cm (第50 d)下降至60 cm (第60 d)时, 20~40 cm剖面土壤体积含水率下降4.11%, 而40~60 cm剖面略微增加; 20~40 cm和40~60 cm剖面处于水位波动带, 其土壤体积含水率的变化较大。0~20 cm剖面位于包气带层, W0和W1处理下该剖面土壤体积含水率出现间隔性的上升和下降, 第10 d、30 d和50 d时土壤体积含水率升高, 这与施肥后灌溉有关。60~80 cm剖面基本处于饱和水层, 水位由80 cm (第10 d)下降至100 cm (第20 d)时, 土壤体积含水率降低; 水位从100 cm逐步上升时, 土壤体积含水率缓慢增加。总体上, 除了0~20 cm剖面以外, 其余剖面土壤体积含水率随着水位的远离而逐渐下降; 随着水位逐步靠近, 土壤体积含水率逐渐上升。

图 2 水位波动下不同施氮量农田的土壤剖面体积含水率变化 Fig. 2 Volumetric water contents at different depths in soil vertical profile under different treatments of nitrogen application rate and groundwater level W0:水位稳定在80 cm; W1:水位每10 d变化20 cm, 10 d、20 d、30 d、40 d、50 d、60 d的水位分别为80 cm、100 cm、80 cm、60 cm、40 cm、60 cm; N0:施氮量为0 kg(N)·hm-2; N1:施氮量为225 kg(N)·hm-2; N2:施氮量为450 kg(N)·hm-2。N1和N2处理第1 d施底肥, 第21 d和41 d追施氮肥。W0: groundwater level is 80 cm; W1: groundwater level is 80 cm, 100 cm, 80 cm, 60 cm, 40 cm and 60 cm on the 10th, 20th, 30th, 40th, 50th and 60th days, respectively. N0: 0 kg(N)·hm-2 nitrogen fertilizer; N1: 225 kg(N)·hm-2 nitrogen fertilizer; N2: 450 kg(N)·hm-2 nitrogen fertilizer. For N1 and N2, base fertilizer was applied on the 1st day, and topdressing nitrogen fertilizer was applied on the 21th and 41th days, respectively.
2.2 水位波动和施氮量对剖面土壤溶液NO3--N含量的影响

图 3可知, 随着施氮量增加各剖面NO3--N含量升高。W0处理组中N1和N2土壤剖面溶液平均NO3--N含量是N0的1.32倍和1.82倍; W1处理组中N1和N2土壤剖面溶液平均NO3--N含量是N0的1.32倍和1.86倍。N0处理下W1土壤剖面溶液平均NO3--N含量是W0的0.90倍。

图 3 水位波动下不同施氮量农田土壤剖面溶液NO3--N含量变化 Fig. 3 NO3--N contents of soil solution at different depths in soil vertical profile under different treatments of nitrogen application rate and groundwater level W0:水位稳定在80 cm; W1:水位每10 d变化20 cm, 10 d、20 d、30 d、40 d、50 d、60 d的水位分别为80 cm、100 cm、80 cm、60 cm、40 cm、60 cm; N0:施氮量为0 kg(N)·hm-2; N1:施氮量为225 kg(N)·hm-2; N2:施氮量为450 kg(N)·hm-2。N1和N2处理第1 d施底肥, 第21 d和41 d追施氮肥。W0: groundwater level is 80 cm; W1: groundwater level is 80 cm, 100 cm, 80 cm, 60 cm, 40 cm and 60 cm on the 10th, 20th, 30th, 40th, 50th and 60th days, respectively. N0: 0 kg(N)·hm-2 nitrogen fertilizer; N1: 225 kg(N)·hm-2 nitrogen fertilizer; N2: 450 kg(N)·hm-2 nitrogen fertilizer. For N1 and N2, base fertilizer was applied on the 1st day, and topdressing nitrogen fertilizer was applied on the 21th and 41th days, respectively.

随着剖面深度的增加, 土壤溶液NO3--N含量呈下降趋势。W0处理, 0~20 cm剖面土壤溶液平均NO3--N含量分别是20~40 cm、40~60 cm和60~80 cm的1.62倍、3.15倍和7.30倍; W1处理, 0~20 cm剖面土壤溶液平均NO3--N含量分别是20~40 cm、40~60 cm和60~80 cm的1.96倍、3.61倍和7.65倍。W1处理0~20 cm平均NO3--N含量是W0处理的0.98倍。

各处理土壤溶液NO3--N含量随着取样时间在剖面分布差异性较大。0~20 cm剖面土壤溶液NO3--N含量在W0N0、W0N1以及W1N0、W1N1处理下随着取样时间总体呈下降趋势, 第60 d与第10 d相比, 分别降低33.59%、14.98%、33.26%和13.33%。W0N2和W1N2处理的NO3--N含量随着取样时间总体呈上升趋势, 第60 d与第10 d相比, NO3--N含量分别增加105.55%和100.79%;从第20 d到第50 d, W0N2和W1N2的NO3--N含量持续增加, 其中第21 d和41 d为氮肥追施期。20~40 cm、40~60 cm和60~80 cm剖面6个处理土壤溶液NO3--N含量随着取样时间均呈下降趋势。第60 d与第10 d相比, W0N0、W0N1、W0N2、W1N0、W1N1和W1N2处理下20~40 cm剖面土壤溶液NO3--N含量分别下降40.66%、39.49%、26.14%、49.40%、44.05%和45.67%, 40~60 cm剖面分别下降74.91%、57.03%、47.61%、88.17%、87.44%和85.60%, 60~80 cm剖面分别下降37.17%、57.81%、56.01%、36.79%、71.46%和69.07%。与W0相比, W1水位波动显著相加了NO3--N含量下降幅度。

W1处理水位上升和下降过程中, 各剖面土壤溶液NO3--N含量变化差异较大。在0~20 cm剖面, 水位从80 cm (10 d)下降到100 cm (20 d)、80 cm (30 d)上升到60 cm (40 d)以及水位从40 cm (50 d)下降到60 cm (60 d), W1N0、W1N1和W1N2处理土壤溶液NO3--N含量均下降; 水位从100 cm (20 d)依次上升到80 cm (30 d)、60 cm (40 d)和40 cm (50 d), W1N2的NO3--N含量持续增加。在20~40 cm剖面, 水位从80 cm (10 d)下降到100 cm (20 d)以及水位从40 cm (50 d)下降到60 cm (60 d), W1N0、W1N1和W1N2处理土壤溶液NO3--N含量均大幅度下降。在40~60 cm剖面, 水位从80 cm (10 d)下降到100 cm (20 d)以及水位从60 cm (40 d)上升到40 cm (50 d)后又下降到60 cm (60 d), W1N0、W1N1和W1N2处理土壤溶液NO3--N含量均大幅度下降。在60~80 cm剖面, 3个施氮处理土壤溶液NO3--N含量下降幅度与水位波动的相关性不明显。

2.3 水位波动和施氮量对地下水NO3--N浓度的影响

图 4可知, 各处理地下水NO3--N浓度随着取样时间均呈上升趋势, 第60 d与10 d相比, W0N0、W0N1、W0N2、W1N0、W1N1和W1N2处理的NO3--N浓度分别增加121.95%、86.72%、69.01%、142.42%、79.82%和45.60%, W1N1和W1N2地下水NO3--N浓度增幅小于W0N1和W0N2。各处理地下水NO3--N浓度随着施肥量增加而增大, W0处理组中, N1和N2各阶段地下水平均NO3--N浓度分别是N0的1.49倍和3.00倍; W1处理组中, N1和N2各阶段地下水平均NO3--N浓度是N0的1.34倍和2.79倍; N0处理下, W1各阶段地下水平均NO3--N浓度是W0的0.90倍。

图 4 水位波动下不同施氮量农田地下水NO3--N浓度变化 Fig. 4 NO3--N concentrations of farmland groundwater under different treatments of nitrogen application rate and groundwater level W0:水位稳定在80 cm; W1:水位每10 d变化20 cm, 10 d、20 d、30 d、40 d、50 d、60 d的水位分别为80 cm、100 cm、80 cm、60 cm、40 cm、60 cm; N0:施氮量为0 kg(N)·hm-2; N1:施氮量为225 kg(N)·hm-2; N2:施氮量为450 kg(N)·hm-2。N1和N2处理第1 d施底肥, 第21 d和41 d追施氮肥。W0: groundwater level is 80 cm; W1: groundwater level is 80 cm, 100 cm, 80 cm, 60 cm, 40 cm and 60 cm on the 10th, 20th, 30th, 40th, 50th and 60th days, respectively. N0: 0 kg(N)·hm-2 nitrogen fertilizer; N1: 225 kg(N)·hm-2 nitrogen fertilizer; N2: 450 kg(N)·hm-2 nitrogen fertilizer. For N1 and N2, base fertilizer was applied on the 1st day, and topdressing nitrogen fertilizer was applied on the 21th and 41th days, respectively.

水位从100 cm (20 d)上升到60 cm (40 d), W1处理的地下水NO3--N浓度持续下降, 但水位从60 cm (40 d)上升到40 cm (50 d)过程中, W1处理的地下水NO3--N浓度上升。W0处理地下水NO3--N浓度变化也呈现一定规律, 在第20 d、30 d和50 d NO3--N浓度增加, 这与21 d和41 d的追施氮素有关。

2.4 水位波动和施氮量对作物产量的影响

图 5可知, 甘蓝产量的变化范围在46 000~67 000 kg·hm-2。施氮量对甘蓝产量的影响较大, 随着施氮量增加甘蓝产量呈现递增趋势。N2和N1处理组分别比N0处理组增产37.4%和27.9%。但是, 相同的施氮量下不同水位波动处理的甘蓝产量差异不大。相关性分析结果得知, 甘蓝产量与施氮量呈显著正相关(P < 0.05), 与水位波动不相关。

图 5 水位波动下不同施氮量农田的作物产量 Fig. 5 Crop yields under different treatments of nitrogen application rate and groundwater level W0:水位稳定在80 cm; W1:水位每10 d变化20 cm, 10 d、20 d、30 d、40 d、50 d、60 d的水位分别为80 cm、100 cm、80 cm、60 cm、40 cm、60 cm; N0:施氮量为0 kg(N)·hm-2; N1:施氮量为225 kg(N)·hm-2; N2:施氮量为450 kg(N)·hm-2。N1和N2处理第1 d施底肥, 第21 d和41 d追施氮肥。W0: groundwater level is 80 cm; W1: groundwater level is 80 cm, 100 cm, 80 cm, 60 cm, 40 cm and 60 cm on the 10th, 20th, 30th, 40th, 50th and 60th days, respectively. N0: 0 kg(N)·hm-2 nitrogen fertilizer; N1: 225 kg(N)·hm-2 nitrogen fertilizer; N2: 450 kg(N)·hm-2 nitrogen fertilizer. For N1 and N2, base fertilizer was applied on the 1st day, and topdressing nitrogen fertilizer was applied on the 21th and 41th days, respectively.
2.5 水位波动和施氮量与土壤溶液和地下水NO3--N浓度的相关性

相关性分析得知, 土壤溶液NO3--N浓度与施氮量呈显著正相关(P < 0.05), 相关系数0.299;土壤溶液NO3--N浓度与土壤体积含水率呈显著负相关(P < 0.05), 相关系数为-0.717。地下水NO3--N浓度与施氮量呈显著正相关(P < 0.05), 相关系数为0.866;地下水NO3--N浓度与水位呈负相关, 但相关性不显著。

3 讨论

本研究中, W0 (水位不变)和W1 (水位波动)两种水位处理下0~20 cm、20~40 cm、40~60 cm和60~80 cm剖面土壤NO3--N含量均随着施氮量增加而升高, 随着土壤剖面深度的增加而降低。当施氮量远超过土壤固持能力和作物吸收量后, 土壤中会累积大量的无机态氮, 其累积量随着施氮量提高呈线性增大[15]。长期施肥会造成NO3--N在表层土壤的大量累积[16-17], 研究发现尿素的迁移转化主要发生在0~5 cm土层内, 5 cm以下土层无机氮含量与对照相比变化并不明显[18]; 也有研究表明土壤氮素养分累积及矿化、硝化等微生物过程主要发生在0~20 cm表层土壤[19]。逐渐累积的无机氮如果不能被作物吸收或转化为有机氮, 则对环境存在潜在威胁[20]。本研究中, 0~20 cm剖面土壤溶液NO3--N含量在N2处理下随着取样时间呈上升趋势, 发生NO3--N累积; N0和N1处理的NO3--N含量随着取样时间变化不明显, 但2次追肥后(第50 d)土壤溶液NO3--N含量有上升现象, 施肥量是影响该剖面土壤硝态氮累积的关键因素。除了追肥期, 土壤溶液高含量NO3--N出现在第10 d (第1次取样期), 这是由于作物生长初期对营养的需求少, 大量基肥的施用, 导致土壤中NO3--N的累积, 同时有机肥带入的大量有机物质还会促进表层土壤硝化过程[21], 导致土壤溶液中NO3--N较高。灌溉后水分入渗及再分布过程均对上层累积的硝态氮具有淋洗作用; 增大施氮量会增加硝态氮向下层淋溶的距离; 在运移过程中一部分硝态氮吸附在土壤中, 导致向下迁移的硝态氮含量降低[22]。这也是本研究中土壤溶液NO3--N含量随着剖面深度的增加呈下降趋势, 深层剖面土壤溶液NO3--N含量随着施氮量增加呈上升趋势的主要原因。

本研究中, W0和W1处理组20~80 cm剖面土壤溶液NO3--N含量均呈持续下降的趋势, 这可能和短期内反复抽取土壤溶液导致土壤溶液提取器(陶土头)周边NO3--N含量持续降低有关。但是在同一施氮量的2种水位处理下, 剖面土壤溶液NO3--N含量下降幅度差异较大, 水位波动(W1)处理显著降低土壤溶液NO3--N含量, 尤其是在20~40 cm和40~60 cm剖面。水位上升(从100 cm上升至40 cm)过程中, 各剖面土壤溶液NO3--N含量与W0相比呈下降趋势, 这与李翔等[12]对NO3--N浓度随着水位波动而变化的研究结果一致, 其认为水位上升导致溶解氧含量降低, 剖面由好氧阶段转变为缺氧阶段, 反硝化作用导致NO3--N浓度下降。然而, 水位下降(从80 cm下降至100 cm以及水位从40 cm下降至60 cm)过程中, 剖面土壤NO3--N随水流向下迁移, 土壤溶液NO3--N含量与W0相比也呈下降趋势。李翔等[12]研究结果与之不同, 其认为水位下降导致剖面土壤溶解氧浓度增大, 硝化作用增强, 土壤NO3--N浓度升高。此外, 本研究中水位下降过程中剖面土壤溶液NO3--N含量下降幅度均大于水位上升阶段, 这可能是由于水位下降速度较快, 剖面土壤NO3--N主要随水流运移至下层, 导致土壤硝态氮含量大幅降低。NO3--N在土壤中随水流运移时, 土壤含水量、水流运动状态、溶解氧含量等土壤理化性质以及微生物群落会影响其迁移转化过程, 包括吸附、反硝化等反应[13]。本研究中, 20~60 cm剖面土壤含水量、水流运动状态、溶解氧含量随着水位波动变化剧烈, 这也是该剖面土壤溶液NO3--N含量变化显著的主要原因。60~80 cm剖面处于淹水区域, 各处理的土壤溶液NO3--N含量均较低。淹水区域土壤NO3--N反硝化作用较强, NO3--N被运移至地下水之前已经被反硝化消耗掉, NO3--N淋溶损失非常少[13], 这可能是该剖面以及地下水NO3--N含量较低的原因。然而, 本研究中水位波动对0~20 cm剖面土壤溶液NO3--N含量影响较小, 可能是因为该剖面处于包气带, 土壤孔隙中溶解氧含量较高, 硝化反应较强, NO3--N含量高, 主要受施肥量影响, 对地下水位波动呼应较弱。

本研究中, W0处理(水位不变)地下水NO3--N浓度呈上升趋势, W1处理地下水NO3--N浓度波动较大, 这与进行水位波动试验时调节进水和出水有关, 未来应加强进入地下水NO3--N总量的研究, 以进一步明确水位波动引起NO3--N对地下水污染风险。水位从60 cm上升至40 cm过程中地下水NO3--N浓度显著高于水位从100 cm上升至60 cm过程。由此得知, 土壤剖面NO3--N浓度高低直接影响了进入地下水的NO3--N的量。此外, 地下水NO3--N浓度与施氮量呈显著正相关, 因此开展源头防控是保护地下水环境的主要措施。文献得知, 水位波动初期, 水分运动相对剧烈, NO3--N随着水分运动迅速在剖面发生分配; 随后, 土壤水分进入再分布过程, 水分运动相对缓慢, NO3--N主要以扩散作用发生迁移, 剖面NO3--N浓度均向下移动并进入地下水[21-22], 因此未来应加强同一水位波动处理下, 不同时间内土壤剖面及地下水NO3--N变化的研究, 明确水位波动下氮淋溶损失及对地下水污染的关键时间节点。

本研究中作物产量与水位波动没有相关性, 这是由于甘蓝属于浅根系植物, 主要吸收表层土壤养分, 在相同灌溉条件下, 施肥量是决定甘蓝产量的主要因素。根据本研究结果, 在地下水埋深较浅的农田, 地下水位波动不会对作物产量产生影响。但是, 土壤剖面NO3--N会随着水位波动进入饱和水层。减少农田施氮量, 降低土壤剖面NO3--N含量, 可降低水位波动下进入地下水NO3--N的量, 保护地下水环境。

本研究采用土壤填装的土柱, 可进行作物种植试验; 根据“U”型管原理, 通过进水阀和出水阀控制水位的升降, 可模拟地下水位波动; 安装的土壤溶液提取器和水分传感器可获取不同剖面土壤淋溶液及土壤水分含量数据。该土柱可真实反映田间状况, 方便科学有效地开展地下水位波动下施氮量对农田土壤NO3--N运移影响的研究。试验设定10 d为一个水位波动期, 水位稳定处理作为水位波动处理的对照, 通过比较同一剖面土壤溶液NO3--N含量变化判断水位波动下不同施氮量农田土壤硝酸盐淋失风险的影响, 研究方法科学、合理, 研究结果可靠。

4 结论

水位波动对不同施氮量农田土壤溶液NO3--N含量有显著影响。包气带土壤溶液NO3--N含量受施氮量影响, 在一定范围内, 随着施氮量增加, 土壤溶液NO3--N含量也显著增加, 过量施氮肥导致NO3--N累积; 水位波动带土壤溶液NO3--N含量受施氮量和水位波动共同影响, 随施氮量增加土壤NO3--N含量增高, 但水位波动能明显降低土壤溶液NO3--N含量。水位上升和下降促进土壤NO3--N随着水流运移到下层, 增加了地下水NO3--N污染风险。地下水NO3--N浓度受施肥和水位变动的共同影响。水位波动对作物产量没有显著影响, 施肥量与作物产量呈显著正相关。

参考文献
[1]
ERISMAN J W, SUTTON M A, GALLOWAY J, et al. How a century of ammonia synthesis changed the world[J]. Nature Geoscience, 2008, 1(10): 636-639. DOI:10.1038/ngeo325
[2]
张卫峰, 马林, 黄高强, 等. 中国氮肥发展、贡献和挑战[J]. 中国农业科学, 2013, 46(15): 3161-3171.
ZHANG W F, MA L, HUANG G Q, et al. The development and contribution of nitrogenous fertilizer in China and challenges faced by the country[J]. Scientia Agricultura Sinica, 2013, 46(15): 3161-3171. DOI:10.3864/j.issn.0578-1752.2013.15.010
[3]
ZHOU J Y, GU B J, SCHLESINGER W H, et al. Significant accumulation of nitrate in Chinese semi-humid croplands[J]. Scientific Reports, 2016, 6: 25088. DOI:10.1038/srep25088
[4]
JU X T, KOU C L, ZHANG F S, et al. Nitrogen balance and groundwater nitrate contamination:Comparison among three intensive cropping systems on the North China Plain[J]. Environmental Pollution, 2006, 143(1): 117-125. DOI:10.1016/j.envpol.2005.11.005
[5]
赵同科, 刘宝存. 泛环渤海地区地下水硝酸盐格局与控制[M]. 北京: 科学出版社, 2014: 64-67.
ZHAO T K, LIU B C. Situation of nitrate in groundwater in regions of the Pan-Bohai Sea and its control[M]. Beijing: Science Press, 2014: 64-67.
[6]
SUGITA F, NAKANE K. Combined effects of rainfall patterns and porous media properties on nitrate leaching[J]. Vadose Zone Journal, 2007, 6(3): 548-553. DOI:10.2136/vzj2006.0095
[7]
KAOWN D, KOH D C, MAYER B, et al. Identification of nitrate and sulfate sources in groundwater using dual stable isotope approaches for an agricultural area with different land use (Chuncheon, mid-eastern Korea)[J]. Agriculture Ecosystems & Environment, 2009, 132(1/3): 223-231.
[8]
ZHANG L D, GAO L H, ZHANG L X, et al. Alternate furrow irrigation and nitrogen level effects on migration of water and nitrate-nitrogen in soil and root growth of cucumber in solar-greenhouse[J]. Scientia Horticulturae, 2012, 138: 43-49. DOI:10.1016/j.scienta.2012.02.003
[9]
ASHWORTH D J, SHAW G. Effect of moisture content and redox potential on in situ Kd values for radioiodine in soil[J]. Science of the Total Environment, 2006, 359(1/3): 244-254.
[10]
ASHWORTH D J, MOORE J, SHAW G. Effects of soil type, moisture content, redox potential and methyl bromide fumigation on Kd values of radio-selenium in soil[J]. Journal of Environmental Radioactivity, 2008, 99(7): 1136-1142. DOI:10.1016/j.jenvrad.2008.01.009
[11]
ZHANG B, GAO X P, LI L, et al. Groundwater depth affects phosphorus but not carbon and nitrogen concentrations of a desert phreatophyte in Northwest China[J]. Frontiers in Plant Science, 2018, 9: 338. DOI:10.3389/fpls.2018.00338
[12]
李翔, 杨天学, 白顺果, 等. 地下水位波动对包气带中氮素运移影响规律的研究[J]. 农业环境科学学报, 2013, 32(12): 2443-2450.
LI X, YANG T X, BAI S G, et al. The effects of groundwater table fluctuation on nitrogen migration in Aeration Zone[J]. Journal of Agro-Environment science, 2013, 32(12): 2443-2450. DOI:10.11654/jaes.2013.12.018
[13]
颜晓元, 周伟. 长江三角洲农田地下水反硝化对硝酸盐的去除作用[J]. 土壤学报, 2019, 56(2): 350-363.
YAN X Y, ZHOU W. Groundwater nitrate removal through denitrification under farmland in Yangtze River Delta[J]. Acta Pedologica Sinica, 2019, 56(2): 350-363.
[14]
LIU J, BI X Q, MA M T, et al. Precipitation and irrigation dominate soil water leaching in cropland in Northern China[J]. Agricultural Water Management, 2019, 211: 165-171. DOI:10.1016/j.agwat.2018.09.053
[15]
DAI J, WANG Z H, LI M H, et al. Winter wheat grain yield and summer nitrate leaching:Long-term effects of nitrogen and phosphorus rates on the Loess Plateau of China[J]. Field Crops Research, 2016, 196: 180-190. DOI:10.1016/j.fcr.2016.06.020
[16]
BASSO B, RITCHIE J T. Impact of compost, manure and inorganic fertilizer on nitrate leaching and yield for a 6-year maize-alfalfa rotation in Michigan[J]. Agriculture Ecosystem & Environment, 2005, 108(4): 329-341.
[17]
张慧霞, 周怀平, 杨振兴, 等. 长期施肥对旱地土壤剖面硝态氮分布和累积的影响[J]. 山西农业科学, 2014, 42(5): 465-469.
ZHANG H X, ZHOU H P, YANG Z X, et al. Effect of long-term fertilization to distribution and accumulation of NO3--N in dryland soil[J]. Journal of Shanxi Agricultural Sciences, 2014, 42(5): 465-469. DOI:10.3969/j.issn.1002-2481.2014.05.14
[18]
张朝, 车玉萍, 李忠佩. 水稻土模拟土柱中肥料氮素的迁移转化特征[J]. 应用生态学报, 2011, 22(12): 3236-3242.
ZHANG C, CHE Y P, LI Z P. Translocation and transformation characteristics of fertilizer nitrogen in paddy soil:A study with simulated soil column[J]. Chinese Journal of Applied Ecology, 2011, 22(12): 3236-3242.
[19]
盖霞普, 刘宏斌, 翟丽梅, 等. 长期增施有机肥/秸秆还田对土壤氮素淋失风险的影响[J]. 中国农业科学, 2018, 51(12): 2336-2347.
GAI X P, LIU H B, ZHAI L M, et al. Effects of long-term additional application of organic manure or straw incorporation on soil nitrogen leaching risk[J]. Scientia Agricultura Sinica, 2018, 51(12): 2336-2347. DOI:10.3864/j.issn.0578-1752.2018.12.010
[20]
胡春胜, 董文旭, 张玉铭, 等. 华北山前平原农田生态系统氮通量与调控[J]. 中国生态农业学报, 2011, 9(19): 997-1003.
HU C S, DONG W X, ZHANG Y M, et al. Nitrogen flux and its manipulation in the cropland ecosystem of the North China Plain[J]. Chinese Journal of Eco-Agriculture, 2011, 9(19): 997-1003.
[21]
刘亥扬.有机物料对稻田土壤硝化活性及N2O排放的影响[D].杭州: 浙江大学, 2019: 47-59
LIU H Y. The effects of organic materials on nitrification activity and N2O emissions in paddy soil[D]. Hangzhou: Zhe-jiang University, 2019: 47-59
[22]
刘秋丽.灌施条件下尿素在土壤中的迁移转化规律研究[D].太原: 太原理工大学, 2012: 1-7
LIU Q L. Migration and transformation of urea in the soil under the conditions of irrigation and fertilization[D]. Taiyuan: Taiyuan University of Technology, 2012: 1-7