阿力曼，张 杰，曾 辉，丁天宇，罗志英，2，李丽丽，胡克林，刘 刚**

(1.中国农业大学 北京 100193； 2.中山火炬高技术产业开发区农业服务中心 中山 528436；3.清华大学 北京 100084)

由于土壤水分蒸发而导致的土壤干缩开裂，普遍存在于自然界中。失水导致的土壤干缩开裂，对土壤学、环境科学、生态科学等相关研究具有重要的学术价值和实际应用价值。首先，裂缝的存在会导致垂向宏观大裂隙的存在，并继而诱发优先流[1]，造成地表污染物随水分快速地进入地下水[2]，并继而污染地下含水层。此外，大量存在的裂隙还会加快包气带水分和溶质向深层迁移，这会降低农业灌溉用水以及肥料的利用效率[3-5]。

我国农田存在氮肥利用率低，氮肥损失严重等问题，大量的氮肥随着土壤水分向下运移。而通过优先流途径的土壤水分和溶质更容易运移到土体深层，增加了污染物向下迁移的几率，从而造成土体和地下水污染。由土壤开裂造成的裂隙流是优先流的一种[6]，开展有关土壤干缩开裂过程的研究，有助于了解土壤开裂过程的主要影响因素，并促进对开裂过程的预测与控制[7-9]。目前，针对土壤裂缝引发的优先流，及其对水分和溶质乃至农田氮磷等养分向深层的淋洗和迁移规律被绝大多数研究者忽略。因此研究优先流存在对于农田土壤的溶质运移具有重要的环境和经济效应。为了定量化研究裂缝及其优先流对土壤水分和溶质及养分迁移的影响，需要提取裂缝三维结构形态。目前，对于野外裂缝形态的提取，有沙子填充、乳胶灌注、CT 扫描等方法。Dasog 等[10]利用沙子填充法，测量开裂裂缝的体积，并用钢丝和麻绳测量宽度和深度等特征，该方法测量精度较低，并且繁琐； Abou 等[11]利用乳胶灌注开裂的土体，乳胶凝固后将其挖出，获得裂隙的三维结构。该方法由于乳胶容易变形，不能准确反映裂缝三维形态，且乳胶易黏附土壤颗粒等杂质，造成失真。Zhang 等[12]利用CT 扫描方法测定裂隙的三维结构，但该方法不适合提取大尺寸样品(通常样品尺寸不足10 cm 直径)且测试分析昂贵。此外，Peng等[2]利用图像处理技术研究土壤裂缝收缩。由于石蜡具有易凝固成型，可以填充土壤细小裂缝、附着力小于乳胶、具有一定的硬度不易变形的特点，本研究选用石蜡，通过将其熔化后灌入土壤裂缝中来提取裂缝空间结构。

在获取土壤裂隙的形态结构参数后，借助软件模拟[13]，即可定量化研究土壤裂隙诱发的优先流对农田水氮运移的影响。目前模拟氮淋失的模型主要有 LEACHM (Leaching Estimation and Chemistry Model)[14]、SWAT (Soil Water Assessment Tool)[15]和RZWQM (Root Zone Water Quality Model)模型[16]等。LEACHM 由康奈尔大学开发，主要描述土壤中水分、氮素及农药等污染物的迁移转化； SWAT 模型是利用自由水中硝酸盐的浓度差来获取地表径流和侧向流中损失的硝酸盐，主要用SCS 径流曲线数法和Green-Ampt 下渗法； RZWQM 模型是由美国农业部农业系统研究所(USDA-ARS，Great Plain System Research Unit)于1992年推出的农业系统作物和环境管理模型，综合考虑了作物根区所有对作物生长的物理、生物和化学过程的影响。以上模型虽然在国内外获得良好的应用，但它们仍然没有考虑到土壤中裂隙对氮淋失的影响。现有裂隙优先流模型大多数是基于理想化裂缝假设，把裂缝近似为楔状体。例如，Hoogmoed 等[17]模拟含有裂隙的土壤的入渗情况，它能较好地预测表层土壤的垂直入渗，但对于裂缝对水分入渗的贡献仍不明朗。Janssen 等[18]提出PADDY-FLUX 模型来评估稻田中裂隙变化对于水分渗漏的影响，但未能量化土体中裂缝存在对于溶质运移的影响。WHCNS (Soil Water Heat Carbon and Nitrogen Simulation)是梁浩等[19]在已有土壤-作物系统模拟的基础上，构建的适用于分析农田土壤水氮运移、作物生长的模型，模拟人工排水农田的地下排水和氮损失，表明地下排水和氮损失总量的模拟效果较好。Liang 等[20]应用该模型对中国西北绿洲地区不同水肥管理模式的灌溉制度进行了优化；Wang 等[21]利用该模型成功模拟了北京地区冬小麦(Triticum aestivum)和夏玉米(Zea mays)轮作系统的水氮运移过程； Li 等[22]应用WHCNS 模型对华北平原小麦-玉米典型轮作区的不同水肥与栽培管理组合模式进行了校验和评估。此外，WHCNS 模型在设施菜地上也有较好的应用[23]。因此，本研究通过WHCNS 模型结合野外裂缝三维形态结构参数，来量化由于土壤裂隙存在而引起的农田土壤中水分和硝态氮的变化。

本研究拟从田间土壤开裂的动态监测出发，提取华北农田裂缝三维几何结构形态特征参数； 通过WHCNS 模型结合野外裂缝三维形态结构参数，量化由于土壤裂隙存在而引起的农田土壤中水分和硝态氮的变化，分析裂缝、施肥方式、灌溉方式以及强降雨对农田氮素运移的影响，为农田水肥的高效管理提供理论依据，并进一步减少肥料的浪费以及地下水的污染。

1 材料与方法

1.1 试验地点和土壤样品

室外试验地点选在中国农业大学上庄实验站。该实验地点位于北京市海淀区，属于温带大陆性气候，年降水量650 mm，土壤以潮土为主。土壤矿物经X 射线衍射仪测定，发现主要成分为石英与长石(表1)。经中国科学院南京土壤研究所液塑限分析仪(GYS-2 型)测试得到： 土壤的液限为0.75 g·g−1，塑限为0.44 g·g−1(本文中使用的含水量均为质量含水量)。利用激光粒径仪(马尔文2000 激光粒度仪)测试土壤样品粒径大小分布(PSD)，其砂粒、粉粒和黏粒(美国制)含量分别为78.49%、19.15%和2.35%。

表1 试验地土壤X 射线衍射仪测定的矿物组成Table 1 Soil mineral composition of the experimental field determined by X-ray diffractometer %

本研究于2019年8—10月在上庄实验站进行自然条件下的土壤开裂试验。试验地由4 个5 m×1.2 m的地块构成，每个地块划分成4 个1 m×1 m 的小区域用于拍摄。其中一个5 m×1.2 m 地块用于裂缝发育动力学特性研究。其他3 块地用于原位提取裂缝三维结构。用相机记录裂缝瞬态发育过程，为了消除后期图像畸变，在拍摄前对每个拍摄区域的4 个角用白色塑料棒标记。为了消除树叶、废塑料、作物秸秆、石子等杂质对裂缝发育的影响，将整个区域30 cm 厚的土壤过10 mm 的筛以去除杂质。

试验田去除杂质并平整之后注水，注水时间为3 h，保证注水结束后地表残留15 cm 左右深的积水，以保证地表以下50 cm 土层完全饱和。然后启动单反相机(600D，佳能)拍摄，拍摄高度为140 cm，并设置拍摄间隔为1 h。

1.2 图像处理

由于视角的原因，相机在拍摄过程中容易导致图像畸变失真以及摄入不相关的背景杂散元素。因此，对拍摄好的照片进行后处理，消除几何变形。如图1A 所示，本文统一截取1 m×1 m 的区域作为研究对象，借助专业图像处理软件Adobe Photoshop CS6中的透视裁剪功能对开裂图片区域进行裁剪、变形消除，最终图像处理效果如图1B 所示。

由于拍照时不可避免地出现土壤表面反光及明暗阴影的干扰，图片需经过滤波、平滑、二值化处理才能用于进一步的分析。本文采用开源图像处理软件Image J 对几何校正后的图片进行二值化处理、平滑去噪等处理，得到了清晰且信噪比高的二值化图像(图1C)。利用Python 程序中自带的morphology子模块中的Skeletonize 函数对土壤裂缝的骨架进行提取，骨架提取后将图像中的像素个数经尺度换算即得到裂缝长度。图像处理过程如图1 所示。

1.3 裂缝发育形态特征的提取

为保证三维裂缝完整、不变形，选用石蜡开展野外试验。将石蜡熔化后，灌入注水后72 h 的干缩开裂土壤的裂缝中，冷却定型后，选择原始裂缝边长为1 m 区域，采取破坏性取样方式将石蜡取出，冲洗表面杂质，得到土壤裂缝形态。本研究利用三维激光扫描仪(天远三维扫描仪OKIO-5M)对取出的石蜡进行三面扫描、处理、合并后，得到土壤裂缝的三维形态。

1.4 WHCNS 模型模拟农田裂缝对水氮运移的影响

使用WHCNS 模型，以天为时间步长，对水分分布和入渗过程分别采用 Green-Ampt 模型以及Richard’s 方程量化。WHCNS 模型中的氮素运移模拟，选用对流-弥散方程，并考虑了矿化、挥发、硝化、反硝化、吸收等源汇项。有关WHCNS 模型详细的参数校验、模型参数输入等细节参考 Liang等[20]。本研究利用校验后的WHCNS 模型模拟了2014年10月—2016年10月期间玉米-小麦轮作的农田表层有无裂隙情况下土壤水分和溶质运移的情况。试验地土壤物理性质和水力学参数如表2 所示，其中的土壤大孔隙度，来自于三维激光扫描仪提取的裂隙等价孔隙度。表层0～10 cm 土层存在裂隙时，其饱和导水率与Cook[24]在2003年提出的介于裂隙之间的饱和导水率的计算方法一致。

表2 WHCNS 模型中土壤水力学性质参数Table 2 Parameters of soil hydraulics in WHCNS model

2 结果与讨论

2.1 土壤干缩裂缝三维形态特征

野外提取土壤三维裂缝形态如图2 所示。裂缝表面没有呈现明显的等水势线及羽状发育的轮廓线，裂缝发育尖端的形态趋于复杂化与无序化。利用软件Geomagic DesignX 测量裂缝的发育的深度、表面宽度及横截面积，并利用Origin 处理数据。对所有提取的三维裂缝几何结构参数的分析显示： 地表处裂缝宽度的最大值为7.19 mm，最小值为3.77 mm，平均值为5.72 mm； 裂缝深度是所得到裂缝在地面以下所能发育的最大值，最大值为101.50 mm，最小深度为76.20 mm，平均值为90.60 mm，其中深度为8～10 cm 的裂缝占绝大部分。通过统计照片中裂缝开裂长度数据，得到每平方米内的平均裂隙长度为4.58 m。

2.2 WHCNS 模型模拟土壤裂缝、施肥方式、灌溉方式以及降雨强度对硝态氮淋洗量的影响

利用WHCNS 模型模拟了2014—2016年的农田生长数据，农田种植制度为冬小麦和夏玉米轮作。分别设置了有裂缝与无裂缝、传统施肥与优化施肥、喷灌与漫灌以及强降雨与正常降雨的4 个对比情景进行模拟，具体模拟设置如表3 所示。其中裂缝对比设置分为传统施肥条件下的模拟1 和2 以及优化施肥条件下的模拟3 和4； 施肥模式对比设置了裂缝存在条件下的模拟1 和3 以及无裂缝条件下的模拟2 和4，传统施肥和优化施肥量分别为560 kg·hm−2和380 kg·hm−2； 同时模拟1 和5 为灌溉模式对比，模拟中设置漫灌次数为4 次，每次漫灌水量为135 mm，喷灌次数为8 次，喷灌水量为45 mm； 模拟1 和6为降雨对比，其中正常年降水量为635 mm，强降雨分别在6、7、8月份设置了3 次强降雨，每次降水量为250 mm。

表3 WHCNS 模型模拟设置以及硝态氮淋洗量Table 3 WHCNS model simulation setting and nitrate leaching volume

2.2.1 施肥方式以及裂缝的存在对硝态氮淋溶的影响

图3a 为模拟期间降水随时间的变化情况以及灌溉和施肥的对应时间刻度(因为优化和传统施肥量不同，故此处仅对应施肥时刻，无施肥量设置)。图3b 和3d 显示了存在裂缝和无裂缝情况下，传统施肥土壤剖面硝态氮含量随时间的变化情况，结果显示裂缝的存在加速了硝态氮的淋洗，导致土壤剖面的硝态氮含量整体下降，年平均淋洗量较无裂缝情况下增加116.91 kg·hm−2。图3c 和3e 显示了存在裂缝和无裂缝情况下，优化施肥处理土壤剖面硝态氮含量随时间的变化情况，裂缝的存在同样也加速了硝态氮的淋洗，年平均淋洗量较无裂缝存在情况下增加142.88 kg·hm−2。由此可见裂缝的存在可明显增加硝态氮的淋洗，对作物生长产生影响。

土壤硝态氮含量同时也受施肥、土壤质地和温度等因素的影响[25]。对比图3b 和3c 可见，在裂缝的存在下，传统施肥造成土壤剖面的硝态氮含量比优化施肥条件下的低，二者硝态氮的年均淋洗量有39.33 kg·hm−2的差异； 在无裂缝存在的条件下，传统施肥的硝态氮年均淋洗量仍然比优化施肥多64.31 kg·hm−2。由此可见，优化田间施肥量对减少氮素淋失具有重要意义。

2.2.2 灌溉方式对硝态氮淋溶的影响

喷灌的水分利用率一般可达90%以上，相比于漫灌省水高达30%～50%，在透水性强的砂质土可省水高达70%[26]。因此，在灌溉方式的模拟中，根据当地农民的灌溉习惯和经验数据，设置每年漫灌4次，漫灌量为135 mm； 喷灌次数设置为8 次，喷灌水量为45 mm。图4b 和4c 分别为漫灌和喷灌条件下土壤剖面硝态氮浓度随时间的变化情况。因为土壤裂缝的存在下，0～10 cm 土层硝态氮含量在施肥后快速增加，漫灌条件比喷灌更加明显。模拟显示，漫灌的硝态氮年均淋洗量比喷灌多4.46 kg·hm−2。

2.2.3 强降雨对硝态氮淋溶的影响

当降雨发生时，土壤水分快速入渗，导致硝态氮在水的驱动下向土壤深层迁移，继而发生淋溶损失，因而降雨是影响土体中硝态氮含量的重要因素。为了探讨强降雨对硝态氮淋失的影响，增加了强降水情景模拟。图5 给出了模拟的一个夏玉米生长季的硝态氮淋洗情况，在正常降水的基础上，分别在6、7、8月份设置了3 次强降雨(250 mm)。通过对比图5b 和5c，强降水发生时，两种情况下的土壤硝态氮淋洗量还是有比较大的差距，强降水条件下硝态氮的淋洗量相较于正常降水条件下增加了83.61%，同时，裂隙的存在也加大了硝态氮淋失风险，使两者的差距进一步扩大。

3 结论与讨论

3.1 讨论

本文借助WHCNS 模型模拟了土壤裂缝、施肥方式、灌溉方式以及降雨强度对土壤硝态氮淋失的影响。通过模拟得知不管是在优化施肥还是传统施肥条件下，裂缝的存在都会对硝态氮的淋失产生重大影响，致使其淋洗量增加了一倍多。王彬俨等[27]为了研究优先流对农田土壤水分和溶质运移的影响，通过室内土柱试验分别模拟了存在优先流和平衡入渗的水分下渗过程，结果显示优先流能够使硝态氮大量、快速下渗，优先流引起的硝态氮累计淋出量占总量的97.60%，这表明优先流的存在提高了硝态氮在土壤中的运移，这与本文的结论具有一致性。随着我国经济发展，为了追求作物的高产，农业生产中一味地增加肥料用量，导致肥料利用率低下和地下水污染。卢数昌等[28]研究了传统施肥和优化施肥对设施芹菜(Apium graveolens)土壤硝态氮运移等方面的影响，结果显示优化施肥不仅能够减少硝态氮向土壤中下部运移，减轻环境风险，同时提高了芹菜品质。本文模拟结果同样显示优化施肥条件减少了硝态氮的淋洗量，与前人的研究结果具有一致性。相对于其他情景，模拟结果显示更改灌溉方式对硝态氮淋失的影响并不明显，这可能与灌溉强度以及模型设置的淋失土层(100 cm)有关。根据赵付江等[29]在不同灌溉方式对土壤硝态氮淋溶影响的研究中指出，在0～30 cm 土层，滴灌下的土壤硝态氮含量明显高于常规沟灌； 在30～90 cm 土层，常规沟灌的硝态氮含量高于滴灌，其中在30～60 cm 土层的差异不显著，60～90 cm 土层差异显著。因为本研究综合探讨100 cm 土层深度的硝态氮淋失量，所以两种灌溉方式所带来的差异并不明显。喷灌和滴灌作为节水灌溉措施，对硝态氮的淋洗效果也类似。喷灌和漫灌作为不同的灌溉方式能改变硝态氮在土层的分布，喷灌因为灌水量小，水分入渗速度慢，使硝态氮更多地保留在根区，供植物充分利用，而大水漫灌的水量大而急，更容易将土壤表层的硝态氮带入更深的土层，使植物无法使用，造成浪费。因此在后期的田间水分管理中，要更多地采用喷灌等方式，并结合合理的施肥方式，使之发挥省水、省肥的优点。高海鹰等[14]利用LEACHM 模型和室内试验分析了不同降水强度对农田土壤氮素淋失的影响，结果显示降水强度越大，硝态氮淋失速率也越大，当径流出现时，水分下渗速率也随着降水的增加而不断增大。但在本研究中，由于WHCNS 模型在模拟时仅计算了淋洗到地下的硝态氮，忽略了强降水带来的氮素径流损失，从而降低了强降水所带来的影响。降水强度是影响氮素淋失的关键因素之一，在田间管理时应合理安排施肥时间和施肥量，特别是降雨集中的夏季应尽量避免在强降雨发生前进行施肥活动。

3.2 结论

在野外通过对田间土壤开裂形态的动态监测以及裂隙的灌制，成功获取了裂隙的三维形态特征。利用三维扫描仪获取的裂缝平均宽度为5.72 mm，平均深度为9.06 cm； 利用图像分析得到试验区土壤每平米裂缝的平均长度为4.58 m。借助WHCNS 模型进行模拟显示，土壤干缩裂缝的形成增加了硝态氮的淋失风险。与无土壤裂缝相比，土壤裂缝的存在，导致硝态氮淋洗量在传统施肥条件下增加116.90 kg·hm−2，在优化施肥条件增加142.88 kg·hm−2；相较于优化施肥模式，传统施肥模式更容易造成硝态氮的淋失风险，且在有或无土壤裂缝存在的情况下表现一致； 模拟喷灌和漫灌模式对硝态氮淋洗影响差异不明显，漫灌模式的硝态氮淋失量仅比喷灌模式多4.46 kg·hm−2，这与喷灌在节水灌溉方面的作用有关； 强降雨的设置同样增加了硝态氮的淋失风险，导致硝态氮的年均淋洗量增加67%。土壤裂缝的存在极大提高了硝态氮的向下运移，因此为了减少硝态氮的淋失，在田间裂缝产生时可通过中耕等措施破坏裂缝，在农田管理时选择优化施肥模式以及喷灌等节水灌溉措施，同时避免在强降雨发生前对作物进行施肥。