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

引用本文 

牛新胜, 张翀, 巨晓棠. 华北潮土冬小麦-夏玉米轮作包气带氮素淋溶机制[J]. 中国生态农业学报(中英文), 2021, 29(1): 53-65. DOI: 10.13930/j.cnki.cjea.200644
NIU X S, ZHANG C, JU X T. Mechanism of nitrogen leaching in fluvo-aquic soil and deep vadose zone
in the North China Plain[J]. Chinese Journal of Eco-Agriculture, 2021, 29(1): 53-65. DOI: 10.13930/j.cnki.cjea.200644

基金项目

国家重点研发计划项目(2016YFD0800102, 2017YFD0200105)和国家自然科学基金项目(41830751, 31861133018)资助

通信作者

巨晓棠, 研究方向为农田生态系统氮素循环与温室气体。E-mail: juxt@cau.edu.cn

文章历史

收稿日期:2020-08-04
接受日期:2020-09-01
华北潮土冬小麦-夏玉米轮作包气带氮素淋溶机制*
牛新胜1, 张翀2, 巨晓棠2,3     
1. 中国农业大学曲周实验站 曲周 057250;
2. 海南大学热带作物学院 海口 570228;
3. 中国农业大学资源与环境学院 北京 100193
摘要:合理水氮管理可以实现作物目标产量和品质、维持土壤肥力和降低环境污染。然而, 自20世纪90年代以来, 我国农田过量施氮和大水漫灌等问题突出, 引起农业面源污染日趋加重, 地下水硝酸盐污染成为一个普遍现象。本文以华北潮土区冬小麦-夏玉米体系为研究对象, 采用数据整合和文献分析的方法, 阐明了典型农田硝态氮淋溶的时空特征及影响因素, 研究了地表裂隙和土壤大孔隙对硝态氮淋溶的影响, 定量了氮素在地表-根层-深层包气带-地下水的垂直迁移通量及过程。结果表明, 农户常规管理的冬小麦-夏玉米轮作体系氮素盈余较高(299~358 kg·hm-2·a-1), 导致土壤根区和深层包气带累积了大量的硝态氮。冬小麦季硝态氮的迁移主要受灌溉影响, 以非饱和流为主, 且迁移距离较短; 春季单次灌溉量低于60 mm, 可以有效控制水和硝态氮淋溶出根区。冬小麦耕作和灌溉引起的地表裂隙对水氮运移的贡献不大。雨热同期的夏玉米季, 土壤水分经常处于饱和状态, 再降雨就可以导致硝态氮淋溶出根层进入深层包气带。夏玉米季极易发生硝态氮淋溶事件(占全年总淋溶事件的81%左右), 硝态氮淋溶量占全年总淋溶量的80%左右, 且单次淋溶事件的淋溶量较高。大孔隙优先流对夏玉米季根区硝态氮淋溶的贡献率在71%左右, 这些硝态氮脱离了作物根系吸收范围, 反硝化作用对硝态氮去除具有一定作用。在华北气候-土壤条件下, 特别应注意冬小麦收获后土壤不应残留过多硝态氮, 以避免夏玉米季降雨发生大量淋溶; 夏玉米季需要注意施氮与作物需氮的匹配。由于夏玉米追肥困难, 生产上提倡一次性施肥措施, 控释肥应该能够发挥更大作用。未来气候变化, 导致夏季极端高强度降雨事件的频率增加, 将会加剧包气带累积硝态氮通过饱和流或优先流向地下水的迁移。合理的水氮管理是从源头上减少硝态氮向深层包气带和地下水迁移的主要措施。
关键词潮土    冬小麦-夏玉米体系    氮素淋溶    硝态氮    包气带    裂隙    大孔隙    优先流    
Mechanism of nitrogen leaching in fluvo-aquic soil and deep vadose zone
in the North China Plain*
NIU Xinsheng1, ZHANG Chong2, JU Xiaotang2,3     
1. Quzhou Experimental Station, China Agricultural University, Quzhou 057250, China;
2. College of Tropical Crops, Hainan University, Haikou 570228, China;
3. College of Resources and Environmental Sciences, China Agricultural University, Beijing 100193, China
Abstract: Rationally managing nitrogen (N) and water results in high crop yield and quality, maintains (or improves) soil fertility, and reduces environmental pollution. However, since the 1990s, excessive use of N fertilizer and flood irrigation has created problems in Chinese croplands, causing agricultural nonpoint source pollution and groundwater nitrate contamination. Data integration and literature review of winter-wheat summer-maize farmlands in fluvo-aquic soil in the North China Plain was used to investigate the temporal and spatial variation of N leaching, the contribution of cracks and macropores to N leaching, and N movement through soil (along the surface to groundwater continuum). The results showed that the N surplus was very high (299-358 kg·hm-2·a-1) when conventional management was used, resulting in high nitrate accumulation in the root and deep vadose zones. Nitrate movement in the winter wheat season was primarily caused by unsaturated flow and affected by irrigation; the nitrate movement distance was short. Water and nitrate loss from the root zone was negligible if the irrigation amount was lower than 60 mm. In the winter wheat season, tillage- and irrigation-induced cracks contributed minimally to nitrate and water movement out of the root zone. In the wet and hot summer maize season, the soil was frequently water-saturated, and small precipitation amounts lead to nitrate leaching, accounting for 81% of the annual leaching events and 80% of the annual nitrate leaching. In the summer maize season, the leached nitrate amounts were much higher than in the winter wheat season, and 71% of the total nitrate leaching was preferential flow caused by macropores. Nitrate from the root zone could be partially removed by denitrification in the deep vadose zone. In the North China Plain, avoiding high nitrate accumulation after winter wheat harvest was effective at decreasing nitrate leaching in the summer maize season. Matching the N fertilizer supply with crop demand and controlled release fertilizer in summer maize season (to avoid costly topdressing N fertilizer) may also play important roles in leaching reduction. Frequent and heavy rainfall accelerates the movement of nitrate via saturated and preferential flow to groundwater. Therefore, rational water and N management is key to reducing nitrate movement to the deep vadose zone and groundwater.
Keywords: Fluvo-aquic soil    Winter wheat-summer maize rotation    Nitrogen leaching    Nitrate nitrogen    Vadose zone    Crack    Macropore    Preferential flow    

合理的水氮管理可以实现作物目标产量和品质、维持土壤肥力和降低环境污染[1]。然而, 自20世纪90年代以来, 我国农田过量施用氮肥和大水漫灌等问题突出, 水肥投入量往往超过作物需求, 施用技术相对粗放, 由此引起的农业面源污染日益加重, 在北方集约化农区主要表现为氮素淋溶及地下水硝酸盐污染[2-5]。2016年水利部《地下水动态月报》报道, 北京、河北、陕西、东北等18个行政区2103眼水井监测结果表明, 80%以上的监测位点地下水污染严重, 其中“三氮”(氨氮、亚硝态氮和硝态氮)污染最为严重。目前, 我国地下水硝酸盐污染已经成为普遍问题, 大多数浅层地下水硝态氮含量接近或超过世界卫生组织(WHO)或美国国家环境保护局(US EPA)的限量标准(10 mg·L-1)[6]。我国潮土、褐土和黑土是氮素淋溶易发区, 典型厚包气带农区超过50%的浅层地下水中硝酸盐浓度超过上述标准, 深层地下水污染程度也在逐年加剧[7]

淋溶是农田氮素损失的重要途径, 主要受气候、土壤和管理因素的影响。在湿润的灌溉作物体系中, 大量的水和硝态氮容易移出作物根区向土壤深层渗漏。华北夏玉米(Zea mays)季降水集中, 极易发生硝态氮淋溶事件[8-9]。土壤性质是影响硝态氮淋溶的重要因素, 如土壤剖面中的黏土层很容易形成饱和层, 有利于减缓硝态氮向深层土壤和地下水的迁移, 也有利于硝态氮通过反硝化去除, 从而减少向地下水的迁移量[10]。土地利用方式和包气带厚度也显著影响土壤硝态氮及地下水硝酸盐含量[11-12]。农田管理措施是影响硝态氮淋溶的重要因素, 高硝态氮淋溶多发生在施氮量较高且根系较浅的作物体系中[13-14]

土壤深层包气带反硝化作用可以减少硝态氮向地下水淋溶。但反硝化脱氮量受土壤理化性质和微生物活性等(如有机碳供应)的限制[15]。厚包气带普遍存在氧气, 也限制了反硝化的发生[12]。因此, 确定硝态氮淋溶速率与深层包气带的反硝化脱氮量, 对评估农作区地下水污染的风险具有重要意义。在深层包气带, 硝态氮储量较大, 硝态氮从移出根区到进入地下水的时间较长, 以上过程在全球氮循环中还未做准确定量[16-17]。由农田耕层向包气带到地下水的全过程氮素淋溶通量(饱和流、非饱和流、优先流)及机制, 是国际农田氮素淋溶和污染控制研究的重点[18]。厚包气带农田在我国农业主产区分布广泛, 长期大水大肥造成土壤硝态氮在厚包气带的大量累积[2, 12, 19]。国内研究主要集中在农田氮素淋溶及其影响因素, 研究方法主要采用渗滤计(Lysimeter)对根层氮素淋溶进行观测和定量, 也分析了气候、土壤类型、种植制度、水肥管理、耕作方式、秸秆还田等对农田土壤氮素迁移及淋溶过程的影响[8, 20-21], 但对氮素在地表-根层-深层包气带-地下水的垂直迁移过程及消减机制, 还没有形成清晰的认识。

华北平原是我国高产农作区[22], 为了满足高产需求, 地下水被过度开采用于灌溉, 同时过量施氮现象比较普遍[2, 23]。浅层地下水的过度开采导致地下水位以0.5~1.0 m·a-1的速度下降[24], 同时地下水硝态氮浓度呈增加趋势[5, 25]。过去的研究主要集中于0~4 m的上层包气带土壤硝态氮累积、转化及淋溶[12, 26], 很少关注深层包气带硝态氮运移及去除过程[15, 27]。在深层包气带, 土壤质地和结构可能存在较大异质性, 影响硝态氮的淋溶、反硝化和地下水补给。例如, 在排水良好的土壤, 土壤硝态氮可能很快地移至地下水, 而被作物再吸收、转化为有机态氮和反硝化的几率很小。而在排水不良, 蒸散量高且地下水补给量有限的土壤, 硝态氮向下运移较慢; 并且在含有较高有机碳或完全饱和条件下发生反硝化脱氮[28-29]

已有研究表明, 硝态氮在包气带迁移的平均速率为38~41 cm·a-1, 且与施氮量关系不密切[24]。硝态氮向下迁移的速率显著低于地下水位下降的速度(100 cm·a-1左右), 这意味着不会有大量的硝态氮通过基质流进入地下水。硝酸盐同位素结果表明, 有机肥和化粪池废弃物是当前地下水硝态氮的主要来源[24]。另外, 华北地下水硝酸盐浓度增加可能是由于砂质土层优先流或者上游地区地下水横向移动造成的。反硝化(尤其是0~4 m土壤)可在一定程度上减少硝态氮向地下水的迁移量, 但深层包气带大量累积的硝态氮以及较弱的反硝化能力, 对未来地下水硝酸盐污染造成了严重的威胁。

那么, 在华北典型的土壤-作物体系下, 水氮在根区、深层包气带是怎么累积和运移的?运移的时空变异是什么?是怎样迁移到地下水的?影响这些过程和机制的主要土壤、气候和农田管理因素是什么?本文结合文献资料和笔者的研究结果, 用大量试验数据定量论述氮素淋溶的时空特征, 土壤裂隙发育对氮素淋溶的影响, 大孔隙特征及优先流对氮素淋溶的贡献, 硝态氮在地表-根层-深层包气带-地下水的垂直迁移过程及消减机制。为依据土壤-气候等不可控因素的时空变异特征, 提出相应的农艺措施, 同时为实现生产目标、环境目标和维持土壤肥力提供科学依据。

1 华北潮土冬小麦-夏玉米轮作氮素淋溶的时空特征

华北平原属典型的温带大陆性季风气候区, 年平均气温14~15 ℃, 年平均降雨量500~900 mm, 其中约80%的降雨集中于夏季的6—10月份[30-31]。近年来极端强降雨事件频发, 如2016年7月发生了单次降雨253 mm的事件, 明显高于1981—2016年间单次降雨的最大量(32~156 mm)[32]。该区域主要种植体系为冬小麦(Triticum aestivum)-夏玉米轮作, 是我国重要的小麦和玉米生产区, 贡献了全国2/3的小麦[33]和1/3的玉米产量[34]。夏玉米季雨量充足, 仅在玉米苗期需要补充灌溉(在土壤含水量较高情况下不需灌水)。由于冬小麦季作物需水量(约450 mm)远高于降雨量(平均280 mm), 需要在冬小麦不同生育期分3~4次灌水200 mm以上, 才能满足冬小麦高产的需求[35-37]。农户常规灌溉方式为漫灌, 灌水量根据经验确定, 且取决于公用灌溉设施(机井、水泵和电力等)的可利用时长, 一般情况下每次灌水量高于100 mm[38], 导致过量灌溉引起硝态氮淋溶出根区, 需要投入更多氮肥以满足作物高产需求[38]。冬小麦-夏玉米轮作全年施氮量为550~660 kg·hm-2·a-1[2], 在一些地区施氮量可高达600~800 kg·hm-2·a-1[39], 远高于冬小麦或夏玉米合理施氮量范围(150~250 kg·hm-2)[40]

除化肥氮外, 氮沉降也是该轮作体系重要氮素输入来源。尤其是近年来, 氮沉降随着工农业向大气中排放氮的增加而增加[41]。华北平原农田的氮沉降量自21世纪初的28 kg·hm-2·a-1[42]增加至21世纪20年代的63 kg·hm-2·a-1[43]。其他氮素输入, 包括有机肥氮、生物固氮、种子氮和灌溉水氮之和也占总氮素投入的10%左右(表 1), 高氮输入导致轮作体系极高的氮素盈余(299 kg·hm-2·a-1, 表 1)。由于华北潮土具有较强的硝化能力, 作物收获后盈余的氮素主要以硝态氮为主, 造成农田根层及以下土壤累积了大量的硝态氮。Meta分析表明, 1980—2015年间华北冬小麦和夏玉米收获后0~4 m土体硝态氮累积量分别高达453 kg·hm-2和749 kg·hm-2[12, 26]。这些硝态氮很难再被作物利用, 且潮土较弱的反硝化强度导致硝态氮在根区反硝化去除的量较小, 大部分随夏季降雨向深层土壤迁移, 直至迁移至地下水[26, 45], 未来气候变化所导致的极端降雨事件也将加速硝态氮的向下迁移过程[12]

表 1 华北平原冬小麦-夏玉米体系氮素收支平衡与盈余 Table 1 Nitrogen budget and surplus in the winter wheat-summer maize rotation system in the North China Plain 

华北农田周年降雨量和气温变化整体表现为“冬季寒冷干燥, 夏季高温多雨”的特点(图 1)[46]。冬小麦和夏玉米土壤水分和硝态氮含量的时空分布受灌水、降雨和作物吸收的共同调控。图 1中的等值线显示了周年根层土壤水分和硝态氮含量的时空分布特征[47]。自冬小麦苗期至分蘖期(1—12月), 耕层土壤含水量为13%~15%;随着土壤深度增加, 土壤含水量呈下降趋势。尽管在冬小麦播种期施用氮肥导致土壤硝态氮含量升高, 但土壤水分较低不易发生优先流和饱和流, 未发生明显的硝态氮淋溶事件。冬小麦越冬期至返青期(1—3月), 根层土壤含水量随蒸发而减少, 硝态氮不易淋溶。冬小麦拔节期(4月初)追氮方式为撒施后灌水, 土壤含水量提高; 由于冬小麦拔节期对氮素需求旺盛, 追肥并未显著提高根层土壤硝态氮含量, 但发生了硝态氮淋溶事件。冬小麦拔节期追肥灌水、抽穗期和灌浆期灌水导致的硝态氮淋溶量分别为1.8 kg·hm-2、2.6 kg·hm-2和2.8 kg·hm-2。夏玉米季频繁的强降雨事件, 导致土壤含水量得到较大增加, 尤其是8月中旬土壤含水量(20%~25%)接近于田间持水量(华北轻壤土的田间持水量范围为20%~24%[48]), 从而使冬小麦季土壤累积的硝态氮, 以及夏玉米季施入的氮肥极易随土壤水分以对流(convective transport)方式迁移出根层土壤[49]。笔者的监测结果表明, 从根层向深层包气带淋溶主要发生在夏玉米季, 北京上庄试验站全年所观察到的21次硝态氮淋溶事件中, 有17次发生在夏玉米季; 夏玉米季单次淋溶事件的硝态氮淋溶量普遍高于冬小麦季, 夏玉米季硝态氮淋溶占全年总硝态氮淋溶量的80%左右[8]

图 1 华北潮土冬小麦-夏玉米轮作周年气候、田间管理、土壤质量含水量、硝态氮累积量及硝酸盐淋溶量 Fig. 1 Annual climate, field management, soil gravimetric water content, nitrate accumulation and leaching in the winter wheat-summer maize rotation in the North China Plain 红色箭头代表发生硝态氮向1 m以下土体迁移, 粗细和数字代表硝态氮淋溶量(kg·hm-2)。降雨量和大气温度来自Gao等[46]。土壤含水量和硝态氮累积量分布图引自李洪亮[47], 硝态氮淋溶数据引自Huang等[8]。The red arrows denote the moving down of nitrate from root zone (0-1 m), the width of arrows and numbers denote the amount of nitrate leaching. Precipitation and air temperature were cited from Gao et al. [46]. Soil water content and soil nitrate accumulation were cited from Li[47], nitrate leaching was cited from Huang et al.[8].
2 华北潮土裂隙发育对氮素淋溶的影响

土壤裂隙发育是土壤内在属性在外界条件变化综合作用的结果, 与土壤物理、化学及生物性质密切相关。裂隙主要是干湿裂隙、冻融裂隙[50]。其中干湿裂隙影响因素主要是黏粒、有机质、干湿循环和人为(如耕作)干扰[50-51]。田间土壤裂隙的存在影响土壤水入渗、能够作为优先流的路径[51]。对华北潮土地表裂隙的发育、特征及对水氮运移的影响研究较少。笔者在中国农业大学曲周实验站, 采用数码相机(Nikon D720)对土壤裂隙拍摄, Photoshop CC 2017对图像中裂隙进行提取(包括几何校正、二值化和降噪处理), 最后用Image J (National Institutes of Health, USA)软件量化拍摄区域的裂隙面积, 计算裂隙率[37, 52]。对于翻耕或旋耕播种的冬小麦田块, 在无冬水灌溉及雨量较小(2016年10月20日冬小麦播种后至2017年4月1日, 单日降雨量在0.1~9.3 mm·d-1)情况下, 没有观测到明显的地表裂隙产生。冬小麦返青期灌水75 mm, 土壤水分含量变化与裂隙呈显著负相关, 与土壤温度变化没有显著关系; 灌水50 h左右土壤裂隙率基本趋于稳定(图 2)。耕作后土壤大水漫灌, 当地表积水落干后, 地表迅速形成裂隙, 裂隙宽度不断增加。华北潮土地表裂隙的发育主要由耕作扰动和灌水(或降雨)后的干湿交替产生。耕层土壤经过耕作、播种等扰动后, 造成土壤紧密结构疏松化[51, 53], 单位体积内的孔隙体积增加; 当有一定量的水输入, 土壤团粒在水力的作用下再互相吸附和凝聚, 扰动疏松增加的空间体积就会形成裂隙, 裂隙发育与土壤水分损失密切相关[51, 53-56]

图 2 模拟灌水后土壤水和温度对裂隙发育的影响[52] Fig. 2 Effects of soil water and temperature on soil crack development after the simulated irrigation[52] 裂隙率=裂隙面积/灌溉小区的总面积×100%。Crack area density=crack area/plot area×100%.

笔者进一步定量了冬小麦季裂隙的宽度、深度和容积。裂隙的宽度采用数码相机(Nikon D720)拍照, 用Photoshop CC 2017按照上述方法对图像进行处理, 随后用OriginPro 9.0软件提取图像数据获取裂隙宽度[57]。裂隙深度采用乳胶灌注, 将环氧树脂和聚硫醇按1:1(体积比)均匀制成乳胶, 随后将其沿着地表裂隙缓慢连续倒入, 直至裂隙灌满为止。填充24 h后, 挖开土壤并取出成型的乳胶, 将其平展地铺放在带有毫米刻度尺子的一张白色背景上, 按照上述定量裂隙宽度的方法, 拍照、图像处理及提取裂隙宽度[57]。裂隙容积定量采用混凝土灌注法, 将过1 mm筛子的水泥和细砂按1:1.2配制混凝土, 并立刻在所选定的具有典型代表性裂隙的1 m×1 m范围内灌注所有地表裂隙, 记录所用灌注混凝土体积。结果表明, 小麦冬季和春季裂隙的平均宽度为(13.5±7.1) mm(数据范围2.3~34.7 mm)(n=85), 平均深度(数据范围75.0±34.8) mm(数据范围28.1~154 mm)(n=85, 旋耕土壤), 平均容积为(971±589) cm3·m-2(数据范围448~2304 cm3·m-2)(n=9)。受乳胶或水泥灌注流动性的限制, 得到的裂隙深度应略低于实际深度。田间观察发现, 农田土壤裂隙深度与土壤扰动深度相当, 笔者推测深耕时土壤裂隙深度一般为25 cm左右, 旋耕时土壤裂隙深度在10 cm左右。

进一步采用亮蓝染色法指示冬小麦季裂隙的深度[52, 58]。设置45 mm、90 mm、135 mm的灌溉量, 灌水48 h后获得染色面积。通过地表灌溉, 亮蓝能够移动到的最深处为30~50 cm, 其中90%以上的染色面积主要集中在0~20 cm(图 3a, 3b)。该试验选择的是一块长期深耕的田块, 耕作深度20~25 cm, 这与前述土壤裂隙深度接近于耕作深度的推测一致。染色深度随着灌溉量的增加而增加, 单次灌溉量45 mm、90 mm、135 mm的染色深度分别达25 cm、30 cm和50 cm。地表裂隙定量和亮蓝示踪结果表明, 地表裂隙对优先流的影响与灌溉量具有交互作用, 但裂隙对优先流的作用不会超过根区, 裂隙对水分和硝态氮移出根区的贡献可能较小。

图 3 模拟灌溉对优先流的影响 Fig. 3 Effect of irrigation on preferential flow under simulated rainfall 图a是亮蓝在土壤垂直剖面的染色情况, 图b是亮蓝在不同层次水平方向的染色面积比例。IS45、IS90和IS135分别指模拟灌水45 mm、90 mm和135 mm, -1、-2和-3分别指3个重复。figure a is the vertical distribution of brilliant blue in the soil profile, figure b is the percentage of dyeing area in the horizontal direction of each soil layer. IS45, IS90 and IS135 denote irrigaiton of 45 mm, 90 mm and 135 mm, respectively.-1, -2 and-3 denote three replications of the same irrigaiton treatment.

笔者在北京上庄实验站对已经产生裂隙的冬小麦田灌水时发现, 灌溉水很快会充满地表裂隙并形成瞬间积水, 也说明了裂隙深度较浅, 没有形成优先流而消耗大量的灌溉水。对开裂土壤进行水分入渗试验表明[59], 开裂土壤水分吸渗过程分为3个阶段:初始裂隙流阶段、稳渗基质流阶段和吸渗末期阶段, 裂隙仅对吸渗第一阶段产生影响, 且开裂土壤的累积入渗量与未开裂土壤无显著差异。因此, 农田土壤扰动和干湿交替发育形成的裂隙, 对水氮向深层土壤运移的影响不大。但是, 冬灌水所造成的地表裂隙, 导致的蒸发面增加对水分散失的影响较大, 尤其是在华北平原多风和干旱的冬春季。

3 华北潮土大孔隙特征及其对氮素淋溶的贡献

土壤溶液的运移是通过溶质扩散和优先流进行的, 饱和基质流是土壤溶液在土壤基质中扩散运动的结果, 而优先流是通过土壤中形成的大孔隙产生的。土壤大孔隙是优先流的通道, 土壤中70%~85%的水分运动源于优先流的贡献[60-61]。土壤中的大孔隙包括像植物根痕那样的管状孔、小动物的洞穴、干燥收缩产生的大裂隙、化学风化溶解产生的空洞、耕作形成的暗洞以及产生于破碎心土层的孔隙, 也包括土壤团粒间较大的孔隙, 其直径大小为0.03~3 mm[62], 也有人认为大于1 mm[63], 孔径大小没有一个确定的标准。

在典型潮土区的河北省曲周县境内, 通过开挖土壤剖面, 在不同层次横截面拍照, 用图像软件Image J提取不同级别的大孔隙面积和数量的方法[64-65], 定量了冬小麦-夏玉米轮作农田根区大孔隙的特征(图 4)。在0~100 cm根区的各层, 细孔(0.008~0.020 cm2)数量和总面积在各等级孔隙中比例最高, 约占80%。细孔隙数量占优势, 而中等或更大的孔隙少。孔隙数量在整个100 cm土体上下分布无显著差异。细孔隙数量或面积均表现为上少下多, 可能与上层土壤频繁扰动有关。

图 4 土壤垂直剖面不同级别大孔隙单位面积数量及比例 Fig. 4 Percentage and numbers of macropores of different sizes in soil profile Ind.表示个数。Ind. denotes individual.

进一步在曲周县第四疃镇选择8个长期冬小麦-夏玉米轮作的农户田块, 采用亮蓝示踪法研究优先流在土壤剖面的迁移及分布。于2017年6月20日, 在马疃村西、马疃村东、第四疃、焦庄、刘庄、王庄、杏园和东魏村8个农户夏玉米田块, 起垄做成3个1 m×1 m的微区, 在微区内均匀撒施亮蓝(每个样方25 g)。待玉米收获后, 以每5 cm为一层, 小心铲平土壤水平面, 用数码相机(Nikon D720)对水平面拍照, 采用上述方法获得亮蓝染色面积和染色百分比。在农户常规管理条件下, 大多数田块优先流主要发生在根区以内, 与前述灌溉模拟亮蓝研究结果一致。亮蓝主要分布在0~50 cm的土层, 染色率在10%以上也集中在0~50 cm深度。但有些田块(王庄和刘庄)优先流可达95~155 cm土层, 到达深层包气带(图 5), 表明优先流发生的空间变异性。

图 5 农户常规管理农田土壤剖面亮蓝染色面积占比 Fig. 5 Preferential flow in the cropland under farmer's conventional management 图例代表不同采样点。Villages in the legend denote sampling sites.

为定量大孔隙优先流对氮素淋溶的贡献, 在中国农业大学曲周实验站用田间土柱渗滤计开展了研究。试验设置4个处理:无大孔隙模拟土柱(CK, 筛分均匀土壤)、少孔隙模拟土柱(AC1, 筛分土壤放1根麻绳)、多孔隙模拟土柱(AC2, 筛分土壤放3根麻绳)和原状土柱(UC, 代表土壤原有孔隙结构, 包括大孔隙)。于2019年6月29日向每个土柱施硝酸钾220 g (相当于280 kg·hm-2的施氮量), 7月6日收集到第1次淋溶液, 随后每天下午17:00收集1次淋溶液, 直至9月18日试验结束。试验过程中观察到, 由于土壤对麻绳的挤压作用, 放置麻绳并没有起到创设大孔隙的作用, 该处理的淋溶液体积和硝态氮淋溶量与CK相当。由于土柱CK的孔隙结构被破坏, 其水氮淋溶不包含优先流途径, 只有基质流。UC是原状土柱, 其水氮淋溶包含优先流和基质流。因此我们推算优先流对硝态氮淋溶贡献率的计算方法为: (UC淋溶量-CK淋溶量)/UC淋溶量×100。

雨季时, 淋溶事件总伴随着降雨事件的发生, 甚至在日降雨量仅为1.4 mm·d-1时, 所有土柱均能收集到淋溶液(图 6)。日降雨量小于10 mm·d-1时, 淋溶一般在降雨的第2天发生且在当天淋溶完毕, 即淋溶具有不连续的特点。当降雨量大于10 mm·d-1, 淋溶可多日连续发生(图 6)。在40次淋溶事件中, UC土柱淋溶量都显著高于CK土柱。在降雨量较低的雨季初期(8月2日以前), UC土柱的淋溶量低于CK土柱; 随着8月2日后降雨量增大, UC土柱淋溶量高于CK土柱(图 6)。不同土柱淋溶液体积与降雨量呈显著正相关, 原状土柱的水和硝态氮淋溶量对降雨量的响应最敏感(拟合直线斜率最大)(图 7)。整个试验期间的降雨量为274 mm, 每个土柱的氮素投入量为1978 mg, UC、AC2、AC1和CK淋溶液体积分别为159.8 mm、39.4 mm、33.2 mm和45.9 mm; 硝态氮的淋溶量分别为89 mg、37 mg、45 mg和24 mg。因此, 水分淋溶率分别为61%、15%、13%和17%, 硝态氮淋溶率分别为4.5%、1.9%、2.3%和1.2%, 优先流对硝态氮淋洗的平均贡献率为71%左右。优先流对硝态氮淋溶的贡献率与整个土体起始水分含量相关, 如果在降雨发生前, 整个土体含水量较低, 降雨首先要润湿土壤, 优先流未达到根区末端, 淋溶事件已经结束, 贡献率较低; 相反, 若土体已达到饱和状态, 则少量的降雨就可发生优先流淋溶。

图 6 优先流对土壤水(a)和硝态氮淋溶(b)的贡献 Fig. 6 Contribution of preferential flow to water (a) and nitrate leaching (b) CK、AC1、AC2和UC分别代表无大孔隙模拟土柱、少孔隙模拟土柱、多孔隙模拟土柱和原状土柱, 折线图代表发生连续的淋溶事件, 点状图代表不连续淋溶事件。CK, AC1, AC2 and UC denote the artificial soil column with no, little and much macropores, and undisturbed soil column, respectively. Line plot denote the continuous leaching event, scatter plot denote the non-continuous leaching event.
图 7 降雨量与土壤淋溶液体积和硝态氮淋溶量的关系 Fig. 7 Relationship among the losses of water and nitrate nitrogen and rainfall CK、AC1、AC2和UC分别代表无大孔隙模拟土柱、少孔隙模拟土柱、多孔隙模拟土柱和原状土柱。K, AC1, AC2 and UC denote the artificial soil column with no, little and much macropores, and undisturbed soil column, respectively.
4 氮素在地表-根层-深层包气带-地下水的垂直迁移过程及消减机制

为了减少硝态氮向地下水的迁移, 需要解析氮素在地表-根层-深层包气带-地下水的垂直迁移过程及通量。包气带是连接地表-根层-地下水的生态屏障, 其上层是根区或者混合区, 水分输入农田后, 与溶质在这一区域内混合, 水分要么以蒸发方式消耗或者被作物根系吸收, 要么以重力形式向下运移[66]。在包气带中, 降雨、土壤水和地下水相互作用频繁, 是土壤水和溶质(如农田化肥氮)运移至地下水的通道[67]

在华北气候-土壤条件下, 冬小麦-夏玉米农户常规施氮量, 通过淋溶移出1 m根区的硝态氮每年可高达209 kg·hm-2(图 8)[24, 68], 主要是因为施氮量超过作物需氮量, 导致大量的氮素盈余[12, 69]。根据上文估算, 土壤大孔隙产生的优先流对硝态氮淋溶贡献率可达71%左右。移出根区进入深层包气带的硝态氮中约有51% (106 kg·hm-2)被反硝化脱掉氮[24], 其余103 kg·hm-2的硝态氮继续向深层土壤(>10 m)和地下水迁移。由于10 m以下包气带土壤的反硝化较弱[24], 移出10 m深度土壤的硝态氮对地下水污染的威胁加大。优化施氮极大减少了氮素移出根区(100 kg·hm-2), 进入深层包气带的硝态氮约有74 kg·hm-2通过反硝化脱氮, 其余仅有26 kg·hm-2的硝态氮向深层土壤或地下水迁移, 该淋溶量应该在地下水不受硝态氮污染可承受的范围[70]

图 8 氮素在地表-根层-深层包气带-地下水的垂直迁移过程及消减机制 Fig. 8 Mechanism of nitrogen movement and removal in the surface soil-root zone-deep vadose zone-groundwater continuum 所有数字单位均为kg·hm-2·a-1, 氮素投入、作物吸收氮引自Wang等[24]和Ju等[2]的结果, 氨挥发依据Ju等[2]的氨挥发损失率估算。根层和深层包气带反硝化根据Wang等[24]的结果估算。根区氮淋溶为氮投入减去作物吸收、氨挥发和根层反硝化。深层包气带淋溶为根区氮淋溶减去深层包气带反硝化。Units for all numbers in the figure are kg·hm-2·a-1. Input, uptake, NH3, Deni.(root), Leaching (root), Deni. (deep), Leaching (root) denote total nitrogen input, crop nitrogen uptake, ammonia volatilization, denitrification and nitrate leaching in the root zone, denitrification and nitrate leaching in the deep vadose zone (> 1 m). Nitrogen input and crop nitrogen uptake was cited from Wang et al. [24] and Ju et al. [2], ammonia volatilization was calculated based on the ammonia volatilization rate of Ju et al. [2], denitrification in the root zone (0-1 m) and deep vadose zone (>1 m) were estimated from Wang et al. [24]. Nitrate leaching from root zone was calculated by deducting crop nitrogen uptake, nitrogen losses via ammonia volatilization and denitrification in root zone from nitrogen input. Denitrification in deep vadose zone was calculated by deducting denitrification in deep vadose zone from nitrate leaching from root zone.

氮素在包气带的迁移通量受农田管理措施的影响。在给定的生产条件下, 作物产量随施氮量增加而达到最高值, 随后不再增加或降低, 但硝态氮淋溶量随施氮量呈现指数增长[1, 71-72]。英国洛桑实验站经典的Broadbalk小麦试验也得出了类似结果[73]。在合理施氮量范围内, 分次施氮可以实现氮素供应与作物需氮在时间和空间上的匹配, 从而提高作物对氮素的吸收[74], 避免土壤在短期内累积大量的硝态氮, 减少硝态氮淋溶的风险[49]。在作物生长的不同时期, 必须保证足够的土壤含水量以满足作物对水分的需求。由于降雨等气候因素不可控, 灌溉(灌水量、频次以及方式)是协调作物需水和硝态氮移出根区之间矛盾的重要手段。灌水量一般需保持土壤含水量在45%~80%的作物可利用水分之间(根据田间持水量和萎蔫点确定)[35, 75]。将传统的漫灌转变为喷灌、滴灌或水肥一体化, 也是减少硝态氮淋溶的重要措施。

5 结论与建议

华北潮土氮素在地表-根层-深层包气带-地下水的垂直迁移过程受气候、土壤和农田管理措施的综合影响, 且具有强烈的时空变异。超过高产作物需求的氮素投入导致土壤-作物体系的氮素盈余, 主要表现为硝态氮在包气带的累积和迁移。冬小麦季硝态氮的迁移主要受灌溉影响, 以非饱和流为主, 且迁移距离较短; 春季单次灌溉量低于60~90 mm, 可以有效控制水和硝态氮淋溶出根区。冬小麦灌溉引起的地表裂隙对水氮运移的贡献不大, 但对冬春季土壤水分的散失值得进一步研究。雨热同期的夏玉米季, 造成土壤水分经常处于饱和状态, 再降雨就可以导致硝态氮淋溶出根层进入深层包气带。大孔隙优先流对夏玉米季根区硝态氮淋溶的贡献率在71%左右。这些硝态氮脱离了作物根系吸收范围, 反硝化作用对硝态氮去除具有一定作用, 但去除量还需要进一步研究。未来气候变化, 导致夏季极端高强度降雨事件的频率增加, 将会加剧包气带累积硝态氮通过饱和流或优先流向地下水的迁移过程。

生产上应同时考虑作物每个生育期的水氮管理, 使土壤水氮供应既能满足高产优质的需求, 又不至于形成过剩而溢出作物根区, 成为对作物无效的水氮而向下继续迁移。在华北的气候-土壤条件下, 特别应注意冬小麦收获后, 土壤根区不应该残留很多的硝态氮, 避免在夏玉米的雨季发生大量淋溶; 在夏玉米季, 需要注意施氮与作物需氮的匹配。由于夏玉米追肥的困难, 生产上提倡一次性施肥措施, 控释肥应该能够发挥更大作用[21]。合理的水氮管理是从源头上减少硝态氮向深层包气带和地下水迁移的主要措施。

参考文献
[1]
巨晓棠, 张翀. 论合理施氮的原则和指标[J]. 土壤学报, 2020.
JU X T, ZHANG C. The principles and indicators of rational N fertilization[J]. Acta Pedologica Sinica, 2020.
[2]
JU X T, XING G X, CHEN X P, et al. Reducing environmental risk by improving N management in intensive Chinese agricultural systems[J]. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(9): 3041-3046. DOI:10.1073/pnas.0813417106
[3]
LIU J G, YANG W. Water sustainability for China and beyond[J]. Science, 2012, 337(6095): 649-650. DOI:10.1126/science.1219471
[4]
LU J, BAI Z H, VELTHOF G L, et al. Accumulation and leaching of nitrate in soils in wheat-maize production in China[J]. Agricultural Water Management, 2019, 212: 407-415. DOI:10.1016/j.agwat.2018.08.039
[5]
王仕琴, 郑文波, 孔晓乐. 华北农区浅层地下水硝酸盐分布特征及其空间差异性[J]. 中国生态农业学报, 2018, 26(10): 1476-1482.
WANG S Q, ZHENG W B, KONG X L. Spatial distribution characteristics of nitrate in shallow groundwater of the agricultural area of the North China Plain[J]. Chinese Journal of Eco-Agriculture, 2018, 26(10): 1476-1482.
[6]
HAN D M, CURRELL M J, CAO G L. Deep challenges for China's war on water pollution[J]. Environmental Pollution, 2016, 218: 1222-1233. DOI:10.1016/j.envpol.2016.08.078
[7]
马洪斌, 李晓欣, 胡春胜. 中国地下水硝态氮污染现状研究[J]. 土壤通报, 2012, 43(6): 1532-1536.
MA H B, LI X X, HU C S. Status of nitrate nitrogen contamination of groundwater in China[J]. Chinese Journal of Soil Science, 2012, 43(6): 1532-1536.
[8]
HUANG T, JU X T, YANG H. Nitrate leaching in a winter wheat-summer maize rotation on a calcareous soil as affected by nitrogen and straw management[J]. Scientific Reports, 2017, 7: 42247. DOI:10.1038/srep42247
[9]
GALLEJONES P, CASTELLÓN A, DEL PRADO A, et al. Nitrogen and sulphur fertilization effect on leaching losses, nutrient balance and plant quality in a wheat-rapeseed rotation under a humid Mediterranean climate[J]. Nutrient Cycling in Agroecosystems, 2012, 93(3): 337-355. DOI:10.1007/s10705-012-9520-2
[10]
PRATT P F, JONES W W, HUNSAKER V E. Nitrate in deep soil profiles in relation to fertilizer rates and leaching volume[J]. Journal of Environmental Quality, 1972, 1(1): 97-101.
[11]
SUCHY M, WASSENAAR L I, GRAHAM G, et al. High-frequency NO3- isotope (δ15N, δ18O) patterns in groundwater recharge reveal that short-term changes in land use and precipitation influence nitrate contamination trends[J]. Hydrology and Earth System Sciences, 2018, 22(8): 4267-4279. DOI:10.5194/hess-22-4267-2018
[12]
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
[13]
HALVORSON A D, BARTOLO M E, REULE C A, et al. Nitrogen effects on onion yield under drip and furrow irrigation[J]. Agronomy Journal, 2008, 100(4): 1062-1069.
[14]
FAN J, HAO M, MALHI S S. Accumulation of nitrate-N in the soil profile and its implications for the environment under dryland agriculture in northern China: A review[J]. Canadian Journal of Soil Science, 2010, 90(3): 429-440.
[15]
CHEN S M, WANG F H, ZHANG Y M, et al. Organic carbon availability limiting microbial denitrification in the deep vadose zone[J]. Environmental Microbiology, 2018, 20(3): 980-992. DOI:10.1111/1462-2920.14027
[16]
ASCOTT M J, GOODDY D C, WANG L, et al. Global patterns of nitrate storage in the vadose zone[J]. Nature Communications, 2017, 8: 1416. DOI:10.1038/s41467-017-01321-w
[17]
JIA X X, ZHU Y J, HUANG L M, et al. Mineral N stock and nitrate accumulation in the 50 to 200 m profile on the loess plateau[J]. Science of the Total Environment, 2018, 633: 999-1006. DOI:10.1016/j.scitotenv.2018.03.249
[18]
SEBILO M, MAYER B, NICOLARDOT B, et al. Long-term fate of nitrate fertilizer in agricultural soils[J]. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(45): 18185-18189. DOI:10.1073/pnas.1305372110
[19]
LI X X, HU C S, DELGADO J A, et al. Increased nitrogen use efficiencies as a key mitigation alternative to reduce nitrate leaching in North China Plain[J]. Agricultural Water Management, 2007, 89(1/2): 137-147.
[20]
刘宏斌. 农田面源污染监测方法与实践[M]. 北京: 科学出版社, 2015.
LIU H B. Monitoring Method and Practive of Agricultural Non-point Source Pollution[M]. Beijing: Science Press, 2015.
[21]
ZHENG W K, WAN Y, LI Y C, et al. Developing water and nitrogen budgets of a wheat-maize rotation system using auto-weighing lysimeters: Effects of blended application of controlled-release and un-coated urea[J]. Environmental Pollution, 2020, 263: 114383. DOI:10.1016/j.envpol.2020.114383
[22]
YUAN Z J, SHEN Y J. Estimation of agricultural water consumption from meteorological and yield data: A case study of Hebei, North China[J]. PLoS One, 2013, 8(3): e58685. DOI:10.1371/journal.pone.0058685
[23]
MUELLER N D, GERBER J S, JOHNSTON M, et al. Closing yield gaps through nutrient and water management[J]. Nature, 2012, 490(7419): 254-257. DOI:10.1038/nature11420
[24]
WANG S Q, WEI S C, LIANG H Y, et al. Nitrogen stock and leaching rates in a thick vadose zone below areas of long-term nitrogen fertilizer application in the North China Plain: A future groundwater quality threat[J]. Journal of Hydrology, 2019, 576: 28-40. DOI:10.1016/j.jhydrol.2019.06.012
[25]
ZHANG W L, TIAN Z X, ZHANG N, et al. Nitrate pollution of groundwater in northern China[J]. Agriculture, Ecosystems & Environment, 1996, 59(3): 223-231.
[26]
JU X T, ZHANG C. Nitrogen cycling and environmental impacts in upland agricultural soils in North China: A review[J]. Journal of Integrative Agriculture, 2017, 16(12): 2848-2862. DOI:10.1016/S2095-3119(17)61743-X
[27]
MIN L L, SHEN Y J, PEI H W, et al. Water movement and solute transport in deep vadose zone under four irrigated agricultural land-use types in the North China Plain[J]. Journal of Hydrology, 2018, 559: 510-522. DOI:10.1016/j.jhydrol.2018.02.037
[28]
LIAO L X, GREEN C T, BEKINS B A, et al. Factors controlling nitrate fluxes in groundwater in agricultural areas[J]. Water Resources Research, 2012, 48(6): W00L09.
[29]
WELCH H L, GREEN C T, COUPE R H. The fate and transport of nitrate in shallow groundwater in northwestern Mississippi, USA[J]. Hydrogeology Journal, 2011, 19(6): 1239-1252. DOI:10.1007/s10040-011-0748-8
[30]
ZHANG Y T, WANG H Y, LEI Q L, et al. Optimizing the nitrogen application rate for maize and wheat based on yield and environment on the Northern China Plain[J]. Science of the Total Environment, 2018, 618: 1173-1183. DOI:10.1016/j.scitotenv.2017.09.183
[31]
WU D R, YU Q, LU C H, et al. Quantifying production potentials of winter wheat in the North China Plain[J]. European Journal of Agronomy, 2006, 24(3): 226-235.
[32]
SONG X T, JU X T, TOPP C F E, et al. Oxygen regulates nitrous oxide production directly in agricultural soils[J]. Environmental Science & Technology, 2019, 53(21): 12539-12547.
[33]
LU C H, FAN L. Winter wheat yield potentials and yield gaps in the North China Plain[J]. Field Crops Research, 2013, 143: 98-105.
[34]
HAN D R, WIESMEIER M, CONANT R T, et al. Large soil organic carbon increase due to improved agronomic management in the North China Plain from 1980s to 2010s[J]. Global Change Biology, 2018, 24(3): 987-1000.
[35]
GAO B, JU X T, MENG Q F, et al. The impact of alternative cropping systems on global warming potential, grain yield and groundwater use[J]. Agriculture, Ecosystems & Environment, 2015, 203: 46-54.
[36]
XIAO D P, SHEN Y J, QI Y Q, et al. Impact of alternative cropping systems on groundwater use and grain yields in the North China Plain region[J]. Agricultural Systems, 2017, 153: 109-117.
[37]
李东润.碳氮投入对潮土生产力及土壤肥力的长期影响[D].北京: 中国农业大学, 2018
LI D R. The effects of long-term C and N inputs on soil productivity and soil fertility of calcareous fluvo-aquic soil[D]. Beijing: China Agricultural University, 2018
[38]
MICHALCZYK A, KERSEBAUM K C, ROELCKE M, et al. Model-based optimisation of nitrogen and water management for wheat-maize systems in the North China Plain[J]. Nutrient Cycling in Agroecosystems, 2014, 98(2): 203-222.
[39]
ZHENG X J, YU Z W, ZHANG Y L, et al. Nitrogen supply modulates nitrogen remobilization and nitrogen use of wheat under supplemental irrigation in the North China Plain[J]. Scientific Reports, 2020, 10(1): 3305.
[40]
巨晓棠, 谷保静. 我国农田氮肥施用现状、问题及趋势[J]. 植物营养与肥料学报, 2014, 20(4): 783-795.
JU X T, GU B J. Status-quo, problem and trend of nitrogen fertilization in China[J]. Journal of Plant Nutrition and Fertilizers, 2014, 20(4): 783-795.
[41]
ZHANG X M, WU Y Y, LIU X J, et al. Ammonia emissions may be substantially underestimated in China[J]. Environmental Science & Technology, 2017, 51(21): 12089-12096.
[42]
张颖, 刘学军, 张福锁, 等. 华北平原大气氮素沉降的时空变异[J]. 生态学报, 2006, 26(6): 1633-1639.
ZHANG Y, LIU X J, ZHANG F S, et al. Spatial and temporal variation of atmospheric nitrogen deposition in North China Plain[J]. Acta Ecologica Sinica, 2006, 26(6): 1633-1639.
[43]
XU W, LUO X S, PAN Y P, et al. Quantifying atmospheric nitrogen deposition through a nationwide monitoring network across China[J]. Atmospheric Chemistry and Physics, 2015, 15(21): 12345-12360.
[44]
赵荣芳, 陈新平, 张福锁. 华北地区冬小麦-夏玉米轮作体系的氮素循环与平衡[J]. 土壤学报, 2009, 46(4): 684-697.
ZHAO R F, CHEN X P, ZHANG F S. Nitrogen cycling and balance in winter-wheat-summer-maize rotation system in Northern China[J]. Acta Pedologica Sinica, 2009, 46(4): 684-697.
[45]
WANG R, PAN Z L, ZHENG X H, et al. Using field-measured soil N2O fluxes and laboratory scale parameterization of N2O/(N2O+N2) ratios to quantify field-scale soil N2 emissions[J]. Soil Biology and Biochemistry, 2020, 148: 107904.
[46]
GAO B, JU X T, SU F, et al. Nitrous oxide and methane emissions from optimized and alternative cereal cropping systems on the North China Plain: A two-year field study[J]. Science of the Total Environment, 2014, 472: 112-124.
[47]
李洪亮.华北典型潮土区小麦-玉米生产体系氮磷淋失控制模式研究[D].北京: 中国农业大学, 2020
LI H L. Study on production patterns of wheat-maize rotation system for controlling of N and P leaching in typical fluvo-aquic soil of North China Plain[D]. Beijing: China Agricultural University, 2020
[48]
陈晓燕, 叶建春, 陆桂华, 等. 全国土壤田间持水量分布探讨[J]. 水利水电技术, 2004, 35(9): 113-116.
CHEN X Y, YE J C, LU G H, et al. Study on field capacity distribution about soil of China[J]. Water Resources and Hydropower Engineering, 2004, 35(9): 113-116.
[49]
CAMERON K C, DI H J, MOIR J L. Nitrogen losses from the soil/plant system: A review[J]. Annals of Applied Biology, 2013, 162(2): 145-173.
[50]
周明涛, 杨森, 秦健坤, 等. 土壤裂隙研究的回顾与展望[J]. 土壤通报, 2017, 48(4): 988-995.
ZHOU M T, YANG S, QIN J K, et al. Review on the research of soil cracks[J]. Chinese Journal of Soil Science, 2017, 48(4): 988-995.
[51]
张中彬, 彭新华. 土壤裂隙及其优先流研究进展[J]. 土壤学报, 2015, 52(3): 477-488.
ZHANG Z B, PENG X H. A review of researches on soil cracks and their impacts on preferential flow[J]. Acta Pedologica Sinica, 2015, 52(3): 477-488.
[52]
索志勇.典型潮土不同种植模式硝态氮运移及裂隙发育机制研究[D].北京: 中国农业大学, 2018
SUO Z Y. Study on NO3- files/image011.png">transport and crack development mechanism for typical calcareous soil in different cropping patterns[D]. Beijing: China Agricultural University, 2018
[53]
GRAY C W, ALLBROOK R. Relationships between shrinkage indices and soil properties in some New Zealand soils[J]. Geoderma, 2002, 108(3/4): 287-299.
[54]
ZHANG Z B, PENG X, WANG L L, et al. Temporal changes in shrinkage behavior of two paddy soils under alternative flooding and drying cycles and its consequence on percolation[J]. Geoderma, 2013, 192: 12-20.
[55]
熊东红, 周红艺, 杜长江, 等. 土壤裂缝研究进展[J]. 土壤, 2006, 38(3): 249-255.
XIONG D H, ZHOU H Y, DU C J, et al. A review on the study of soil cracking[J]. Soils, 2006, 38(3): 249-255.
[56]
BANDYOPADHYAY K K, MOHANTY M, PAINULI D K, et al. Influence of tillage practices and nutrient management on crack parameters in a vertisol of central India[J]. Soil and Tillage Research, 2003, 71(2): 133-142.
[57]
NAJM M R A, JABRO J D, IVERSEN W M, et al. New method for the characterization of three-dimensional preferential flow paths in the field[J]. Water Resources Research, 2010, 46(2): W02503.
[58]
郭婧妤.华北平原农田氮素污染评价方法及指标研究[D].北京: 中国农业大学, 2017
GUO J Y. Study on evaluation method and index of nitrogen pollution in North China Plain[D]. Beijing: China Agricultural University, 2017
[59]
张展羽, 范世敏, 王策, 等. 干缩开裂黄褐土非饱和水平吸渗特征研究[J]. 农业机械学报, 2017, 48(7): 214-221.
ZHANG Z Y, FAN S M, WANG C, et al. Characteristics of unsaturated horizontal infiltration in desiccation cracks yellow cinnamon soil[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(7): 214-221.
[60]
WILSON G V, LUXMOORE R J. Infiltration, macroporosity, and mesoporosity distributions on two forested watersheds[J]. Soil Science Society of America Journal, 1988, 52(2): 329-335.
[61]
WATSON K W, LUXMOORE R J. Estimating macroporosity in a forest watershed by use of a tension infiltrometer[J]. Soil Science Society of America Journal, 1986, 50(3): 578-582.
[62]
秦耀东, 任理, 王济. 土壤中大孔隙流研究进展与现状[J]. 水科学进展, 2000, 11(2): 203-207.
QIN Y D, REN L, WANG J. Review on the study of macropore flow in soil[J]. Advances in Water Science, 2000, 11(2): 203-207.
[63]
牛健植, 余新晓. 优先流问题研究及其科学意义[J]. 中国水土保持科学, 2005, 3(3): 110-116.
NIU J Z, YU X X. Preferential flow and its scientific significance[J]. Science of Soil and Water Conservation, 2005, 3(3): 110-116.
[64]
毕利东, 张斌, 潘继花. 运用Image J软件分析土壤结构特征[J]. 土壤, 2009, 41(4): 654-658.
BI L D, ZHANG B, PAN J H. Analysis of soil structural properties by using Image-J software[J]. Soils, 2009, 41(4): 654-658.
[65]
PASSONI S, DA SILVA BORGES F, PIRES L F, et al. Software Image J to study soil pore distribution[J]. Ciência e Agrotecnologia, 2014, 38(2): 122-128.
[66]
庞忠和, 黄天明, 杨硕, 等. 包气带在干旱半干旱地区地下水补给研究中的应用[J]. 工程地质学报, 2018, 26(1): 51-61.
PANG Z H, HUANG T M, YANG S, et al. The potential of the unsaturated zone in groundwater recharge in arid and semiarid areas[J]. Journal of Engineering Geology, 2018, 26(1): 51-61.
[67]
敦宇, 张焓, 杨征伦, 等. 再生水灌溉区农田包气带及地下水中硝酸盐分布特征及来源研究[J]. 中国生态农业学报(中英文), 2019, 27(11): 1725-1731.
DUN Y, ZHANG H, YANG Z L, et al. Distribution characteristics and sources of nitrate in the unsaturated zone and groundwater of farmlands in an area irrigated with reclaimed water[J]. Chinese Journal of Eco-Agriculture, 2019, 27(11): 1725-1731.
[68]
李晓欣, 胡春胜, 张玉铭, 等. 华北地区小麦-玉米种植制度下硝态氮淋失量研究[J]. 干旱地区农业研究, 2006, 24(6): 7-10.
LI X X, HU C S, ZHANG Y M, et al. Losses of nitrate-nitrogen from a wheat-corn rotation in North China[J]. Agricultural Research in the Arid Areas, 2006, 24(6): 7-10.
[69]
ZHANG C, JU X T, POWLSON D, et al. Nitrogen surplus benchmarks for controlling N pollution in the main cropping systems of China[J]. Environmental Science & Technology, 2019, 53(12): 6678-6687.
[70]
巨晓棠, 谷保静. 氮素管理的指标[J]. 土壤学报, 2017, 54(2): 281-296.
JU X T, GU B J. Indexes of nitrogen management[J]. Acta Pedologica Sinica, 2017, 54(2): 281-296.
[71]
CUI Z L, ZHANG H Y, CHEN X P, et al. Pursuing sustainable productivity with millions of smallholder farmers[J]. Nature, 2018, 555(7696): 363-366.
[72]
CHEN X P, CUI Z L, FAN M S, et al. Producing more grain with lower environmental costs[J]. Nature, 2014, 514(7523): 486-489.
[73]
GOULDING K W T, POULTON P R, WEBSTER C P, et al. Nitrate leaching from the broadbalk wheat experiment, rothamsted, UK, as influenced by fertilizer and manure inputs and the weather[J]. Soil Use and Management, 2000, 16(4): 244-250.
[74]
WANG X B, ZHOU W, LIANG G Q, et al. The fate of 15N-labelled urea in an alkaline calcareous soil under different N application rates and N splits[J]. Nutrient Cycling in Agroecosystems, 2016, 106(3): 311-324.
[75]
MENG Q F, SUN Q P, CHEN X P, et al. Alternative cropping systems for sustainable water and nitrogen use in the North China Plain[J]. Agriculture, Ecosystems & Environment, 2012, 146(1): 93-102.