不同灌溉方式下施用生物炭对 土壤水盐运移规律及玉米水分利用效率的影响

贾咏霖，屈忠义*，丁艳宏，杨 威，马贵仁，李争争

（1.内蒙古农业大学 水利与土木建筑工程学院，呼和浩特 010018； 2.西北农林科技大学 水利与建筑工程学院，陕西 杨凌 712100）

0 引 言1

【研究意义】水资源是维持生态系统的基本要素，同时也是支撑人类社会发展的重要战略资源[1]。我国西北干旱区面积占全国面积24.5%，但水资源仅占全国5%左右，农业一直面临严峻的水资源短缺问题。内蒙古河套灌区是我中国的粮食主产区，为保障我国的粮食安全做出了重要贡献。然而该区土壤蒸发强烈，土壤盐碱化问题严重，且灌溉水利用效率较低[2]。灌区目前主要依靠黄河水进行漫灌和畦灌，这种灌溉方式尽管能有效减轻土壤盐碱化问题[3]，但灌水量较大，农田水分利用效率较低，水资源浪费严重[4]。滴灌是河套灌区目前广泛推广应用的节水灌溉方式，其不仅能提高作物水分利用效率和产量，同时也能改变土壤水盐运移规律[5-6]。减轻土壤盐碱化程度、提高农田水分利用效率。【研究进展】生物炭是由植物生物质在高温缺氧条件下裂解碳化而成的一种土壤改良剂[7]。已有研究表明，生物炭不仅能改变土壤理化性质[8]，同时还具有较好的固碳节肥能力[9-10]。生物炭具有的多孔结构及较大表面积使得其具有较强的阳离子交换能力[11]，应用生物炭对盐渍土进行改良，对有效缓解土壤含盐量积累，减少盐胁迫对作物生长的不利影响[12-13]。国内外学者也对施加生物炭条件下土壤盐分、水分运动进行了大量研究。朱成立等[14]在滨海垦区研究畦灌条件下咸淡交替灌溉对土壤盐分分布时发现，生物炭可以促进玉米光合作用、减轻水分胁迫、避免离子毒害，最终提高产量和收获指数。岳燕等[15]在通过土柱模拟畦灌条件下施加生物炭对土壤盐分淋洗作用时，发现施加生物炭可提高土壤洗盐效率，将电导率淋洗至5 mS/cm 的时间缩短了41～100 d。也有研究表明，施加生物炭可提高土壤蓄水持水能力[16]，显著提高沙土土壤有效含水率，减少黏土有效含水率[17]。安艳等[18]、李昌见等[19]研究表明，在地下水滴灌条件下施加生物炭可有效提高土壤团聚体的吸湿系数，增加了土壤总孔隙度和田间持水率。岑瑞等[20]研究发现，黏土中施加生物炭30 t/hm2时土壤含水率增加9%、累积入渗量增加45.5%。魏永霞等[21]通过野外大田试验发现施加生物炭可提高土壤持水保水能力，土壤含水率提高7.1%～14.5%，随着施碳年限越长，其效果越弱。综上所述，施加生物炭不仅可以减低土壤含盐量，加快盐分淋洗，同时也可以有效提高土壤持水保水、土壤入渗能力。【切入点】目前河套灌区针对土壤水盐运移规律的研究主要集中在不同覆膜方式、节水灌溉等方面，对比不同灌溉方式下水盐运移规律的研究较少，针对不同灌溉条件下单次施加生物炭对土壤水盐运移特征及作物水分利用效率的中长期综合影响效应尚不清楚。【拟解决的关键问题】本文旨在研究不同灌溉方式下，秸秆生物炭添加对土壤水盐运动规律的中长期影响效应，提出生物炭田间应用参数，为河套灌区水资源高效利用及盐碱化耕地生态治理提供理论依据和技术支持。

1 材料与方法

1.1 试验区概况

试验于2016—2018 年每年5—10 月在内蒙古河套灌区临河区双河镇进步村九庄农业合作社（107°18′E，40°41′N，海拔1 042 m）进行。该地属于中温带半干旱大陆性气候，云雾少、日照充足、降水量少、风大气候干燥。多年年均降水量140 mm 左右，年平均气温6.8 ℃，平均日照时间3 229.9 h，无霜期130 d 左右。经计算[22]，玉米生育期蒸散量为891.3 mm，非生育期蒸散量457.8 mm。生育期降雨量126.6 mm，非生育期降雨量23.4 mm。生育期地下水平均埋深为3.04 m。试验地地势东高西低，地面坡降1/6 000。供试土壤为砂壤土（0～60 cm），基本理化性质见表1。

1.2 试验材料

生物炭来源于辽宁金和福农业开发有限公司，采用玉米秸秆为原料，在450 ℃无氧条件烧制而成（粒径1.5～2.0 mm），其基本理化性质见表1。

1.3 小区试验设计

以灌溉方式、玉米秸秆生物炭用量为二因素进行完全随机区组小区试验设计，其中灌溉方式为地下水滴灌、地下水畦灌、黄河水畦灌，生物炭用量为0、30 t/hm2，共设置6 个处理：地下水滴灌处理（D）、黄河水畦灌处理（H）、地下水畦灌处理（J）、地下水滴灌+生物炭处理（DB）、黄河水畦灌+生物炭处理（HB）、地下水畦灌+生物炭处理（JB）。每个处理设置3 个重复，共18 个试验小区，每个小区面积90 m2（15 m×6 m）。2016 年4 月将生物炭均匀撒施于土壤表层，而后用旋耕机将生物炭与耕层20 cm 土壤均匀混合，2017 年和2018 年不再重复施加生物炭。

表1 土壤、生物炭的基本理化性质 Table 1 Principal physical and chemical properties of soil and biochar

表2 灌溉水的水质状况 Table 2 Irrigation water quality

供试玉米品种为西蒙6 号，种植的株距为30 cm，行距为60 cm，密度为56 667 株/hm2。播种前基施磷酸二铵（N 质量分数18%，P 质量分数46%）450 kg/hm2，氮磷钾复合肥（N 质量分数15%，P 质量分数15%，K 质量分数15%）337.5 kg/hm2。拔节期前、中、后期追肥3 次，抽雄期1 次，灌浆期1 次，每次追尿素（N质量分数46%）75 kg/hm2，共追施尿素375 kg/hm2。膜下滴灌的水源为黄河水，滴灌带为内镶片式滴灌带，管径16 mm，滴头流量1.38 L/h，滴头间距0.3 m。滴灌带铺设于地膜中间，每条滴灌带控制2 行玉米。通过埋于膜内滴头正下方20 cm 处的张力计控制灌溉，灌水下限-25 kPa，灌溉定额270 m3/hm2。分别于2018年5 月30 日（150 d）、6 月27 日（178 d）、8 月6 日（218 d）进行灌溉黄河水畦灌、黄河水滴灌，共灌溉3 次，灌溉量均参照当地农业灌溉实际情况，滴灌单次灌水22.5 mm，畦灌单次灌水100 mm。

1.4 样品采集与数据分析

2018 年玉米生长季每隔10～15 d，使用土钻分别在膜上（滴头正下方位置）、膜中（距滴头17.5 cm 处）、膜间（距滴头35 cm 处）、膜外（距滴头60 cm 处）4个位置进行土壤取样（图1），其中0～40 cm 土层每10 cm 采集1 次，40～100 cm 土层每20 cm 采集1 次，每个处理重复3 次，每次共取1 512 个土样（6 个处理×3 个重复×4 个取样点×7 层）。取回的土壤样品分成二部分，一部分通过烘干法测定土壤含水率，另一部经自然风干后研磨过筛（2 mm），使用便携式全自动电导率分析系统（DS-307，上海霖诺环保科技有限公司）测定土壤浸提液的电导率根据土壤含盐量与电导率的经验转换公式S=0.349×EC1:5，估算土壤全盐量[23]。并在试验过程中监测与分析土壤的水盐平衡状况，计算脱盐率和作物水分利用效率，各指标的计算式如下：

图1 土壤水盐取样位置示 Fig.1 Soil water-salt sampling points

1）盐分平衡分析

对1 m 土层进行盐分平衡分析是土壤盐分累积分析的重要方法，根据不同取样点进行加权平均计算[24]，公式如下：

式中：S 为盐分总量（kg/hm2）；Sxy为土壤取样点含盐量（g/kg）；Ly为取样控制土层宽度（m）；hx为取样点土层厚度（m）；γx为相应点土壤体积质量（g/cm3），i为计算深度内取样层数；n 为控制宽度1.2 m。

2）脱盐率

式中：RSDR 为脱盐率（%），RSDR 为正值时表示土层盐分增加，脱盐率为负值时表示土层盐分减少，S1为2016年播种前土层盐分初始值（kg/hm2），S2为2018年对应土层秋收后含盐量（kg/hm2）。

3）土壤贮水量

式中：H 为土壤贮水量；Δθi为土层i 的土壤质量含水率（%）；Ci为土层i 的土层体积质量（g/cm3）；Zi为土层i 厚度（mm）。

4）作物生育期耗水量

式中：P 为作物全生育期的降雨量（mm）；I 为作物生育期的灌溉定额（mm）；K 为地下水补给量（mm）；RO 为地表径流（mm），因为试验用地平整，故不考虑地表径流；D 为深层渗漏量（mm），其中地表滴灌，每次灌水量较少，顾不考虑深层渗漏量；畦灌深层渗漏量经计算为11.23 mm，ΔH 为0～100 cm 土壤剖面贮水量变化量（mm）。

5）玉米水分利用效率及灌溉水分利用效率

式中：WUE 为水分利用效率（kg/m3）；Y 为玉米籽粒产量（kg/hm2）；ETc为玉米全生育期的耗水量（m3/hm2）。

式中：IWUE 为灌溉水水分利用效率（kg/m3）；Y 为玉米籽粒产量（kg/hm2）；I 为玉米全生育期的灌溉水量（m3/hm2）。

6）收获指数

式中，HI 为收获指数；Y 为玉米籽粒产量（kg/hm2）；DB 为玉米全生育期的干物质累积量（kg/hm2）。

1.5 数据处理方法

采用Microsoft Excel 2016 软件对原始观测数据进行整理与分析，使用Origin 9.1 和Sufer 12.0 软件绘图，利用SPSS 20.0 对整理的数据进行二因素方差分析（P<0.05）。

2 结果与分析

2.1 土壤水分时空动态变化规律

不同观测日期0～100 cm 土层土壤含水率动态变化规律如图2 和图3（注：箱线图中箱体上黑线表示均值，黑点表示特异质，箱体代表75%～25%，误差线代表最大值和最小值）所示。从图2 可以看出，不同灌溉方式下土壤含水率变化规律不同，在不施用生物炭的条件下，生育前期表层地下水滴灌处理（D）的土壤含水率明显高于地下水畦灌（J）和黄河水畦灌处理（H），在生育中后期滴灌处理与畦灌处理含水率变化趋势相近，整体趋于平稳状态。相同灌溉条件下，生物炭处理的土壤含水率整体上高于未施用生物炭的空白对照小区。继续分析发现（图3），不同灌溉方式下，生物炭处理施用后0～20、60～80 cm 的土壤含水率波动范围明显高于未施加生物炭处理，20～60 cm 土层施加生物炭处理的土壤含水率变化范围与未施加生物炭处理的基本相同，而80～100 cm 土层各处理的土壤含水率波动范围较小。

图2 不同处理下0～100 cm 土壤含水率动 Fig.2 Dynamic of soil moisture content of 0～100 cm under different treatment

图3 2018年不同处理下全生育期0～100 cm土壤含水率箱线图 Fig.3 Results of 0～100 cm soil water content in the whole growth period under different treatment

2.2 土壤盐分运移规律

玉米收获期不同灌溉条件下土壤盐分分布情况如图4 所示。由图4 可知，在相同灌溉方式下，施加生物炭处理的小区较未施加生物炭的小区更能减少0～100 cm 的盐分累积。地下水滴灌方式下，DB 处理积盐主要集中在膜外位置，D 处理盐分主要集中在膜上距滴头60 cm处和膜外表层土壤。地下水畦灌处理，JB 处理的含盐量变化较小，J 处理盐分主要集中膜间位置50 cm 处。黄河水畦灌处理，HB 处理剖面盐分小于90.7 g/kg，H 处理的盐分在滴头正下方60 cm 处达到最大值。畦灌条件下土壤盐分分布较为均匀，而地下水滴灌处理盐分主要集中于膜外及膜间位置，这主要是因为地下水滴灌处理单次灌水量较小，单次灌水入渗区域在0～30 cm 土层，无法均匀将盐分淋洗至深层土壤。

图4 不同处理下生育末期0～100 cm 土壤全盐含量/(g/kg) Fig.4 Total salt content in 0～100 cm soil at the end of growing stage under different treatment

对根系主要集中的0～40 cm 土层及0～100 cm 盐分较2016 年改变量进行进一步统计分析（表3），发现不同处理0～40 cm 土层盐分的改变量相差较大（31.5～345.2 kg/hm2），脱盐效果存在显著差异（3.1%～33.6%）。不同灌溉方式下，施加生物炭处理0～40 cm 和0～100 cm 的脱盐率明显小于未施加生物炭处理，其中地下水滴灌+生物炭处理降盐效果最优，其次为地下水畦灌+生物炭处理，黄河水畦灌+生物炭处理的脱盐效果最差。进一步分析发现，地下水滴灌条件下土壤施用生物炭的相对脱盐率（0～100 cm）比地下水畦灌、黄河水畦灌下施用生物炭处理高出1.24倍和3.30 倍。

图5 为0～40 cm 及0～100 cm 各处理膜内、膜中、膜间、膜外盐分较2016 年改变量。0～40 cm 土层不同灌溉方式下生物炭处理对不同取样位置土壤盐分的影响规律存在差异，D 处理在膜上、膜中、膜间、膜外位置的积盐量分别比DB 处理少158.3、1 263.4、776.2、-499.7 kg/hm2，JB 处理在膜上、膜中、膜间、膜外位置的积盐量分别比J 处理少481.9、523.9、719.6、-692.1 kg/hm2，其中在膜上、膜中、膜间位置差异显著。HB 处理在膜上、膜中、膜间、膜外位置的积盐量分别比H处理少1 036.7、709.9、-87.3、-903.9 kg/hm2，其中在膜上、膜中、膜外位置较H 处理差异显著。不同灌溉处理中，滴灌处理盐分累积主要集中在膜外位置，盐分从膜外向膜上位置逐渐减少。地下水畦灌处理盐分累积量较小，盐分分布较稳平均，而黄河水畦灌处理整体上呈脱盐状态。

表3 不同处理下0～40 和0～100 cm 土层盐分的改变量及脱盐率计算 Table 3 Calculation of salinity change and desalination rate in 0～40 and 0～100 cm soil layers under different treatment

图5 不同处理下0～40 cm 和0～100 cm 土层剖面脱盐量 Fig.5 Desalination of 0～40 cm and 0～100 cm soil profiles under different treatments

2.3 玉米产量性状与水分利用效率

由灌水量、蒸散量、产量及产量构成要素推算水分利用效率、灌溉水分利用效率及收获指数如表4 所示。不同灌溉方式下施加生物炭与未施加生物炭处理的作物蒸散量没有显著差异。地下水滴灌、黄河水畦灌条件下施加生物炭与未施加生物炭处理水分利用效率、灌溉水利用效率差异显著。作物产量在施加生物炭与未施加生物炭处理间差异不显著，但是地下水滴灌条件下产量显著高于未施加生物炭处理。而产量构成由百粒质量及穗粒数共同决定，不同处理间百粒质量差异不显著，而滴灌处理穗粒数显著高于畦灌处理。干物质量除黄河水畦灌条件下施加生物炭处理与未施加生物炭处理无显著差异外，地下水滴灌、地下水畦灌条件下均表现出显著差异。不同灌溉方式下，滴灌处理较畦灌处理可减少蒸散量的同时提高水分利用效率及灌溉水利用效率，减少水资源的浪费。同时由于滴灌处理少量高频的灌溉方式，使得其产量显著高于畦灌处理。

3 讨 论

3.1 不同灌溉方式下施加生物炭对作物水分利用效率的影响

施用生物炭可明显提高土壤含水率，有利于增加土壤储水量。本研究运用水量平衡方程对作物耗水量进行分析，发现施加生物炭有效减少地下水滴灌条件下土壤蒸散量，这是因为施加生物炭增加了土壤孔隙度，使得土壤水分更多地被生物炭吸附留存在土体中，从而减少了土壤无效蒸散量，提高了土壤储水量，可被作物利用的有效水增加，这与魏永霞等[25]研究结果一致。畦灌条件单次灌水量大，灌水间隔时长，表层土壤含水率变化较小，使得畦灌条件下施加生物炭处理与未施加生物炭处理土壤蒸散量差异较小。各灌溉方式下灌水量相同，产量直接影响灌溉水利用效率。滴灌处理显著增加了作物产量，这主要由于穗粒数增加显著，少量高频的灌溉模式可以使土壤保持相对稳定的含水率，在玉米需水量较大的抽雄期很好的保障了土壤水分供应，从而进一步促进作物了生长，充足的花粉及伸展的花丝最终导致了玉米秃尖长的减少及穗粒数的增加[26]。地下水滴灌条件下施加生物炭可有效减少土壤蒸散量，增加作物产量，从而提高水分利用效率及灌溉水利用效率。

3.2 不同灌溉方式下施用生物炭对土壤盐分运移的中长期影响

研究表明，盐分胁迫是影响作物正常生长的主要因素之一[27]，生物炭独特多孔隙结构、较大的氧化表面积使得其具有较强的吸附能力[28]、通过促进土壤矿物复合体形成以提高土壤团聚体[21]，生物炭施用多年后对土壤盐分的减少主要有以下2 个原因：①生物炭本身含有大量的可溶性阳离子，大量的K+、Ca2+、Mg2+等离子交换了盐渍化土壤胶体吸附的Na+离子[29]，同时大量的可溶性阳离子作为作物矿物来源的一部分被植物吸收利用[30]。②由于生物炭自身含有大量有机质，施入土壤后有效改善了土壤物理化学性质，增加有机质量，进而直接参与土壤形成、改变土壤吸水吸肥过程，促进土壤团聚体的形成[31]，降低土壤体积质量，改变土壤持水性，在生物炭处理中土壤体积质量小、孔隙度大，促进土壤水分运动，所以盐分更容易被淋洗到土壤下层。本试验研究发现，地下水滴灌条件下施加生物炭处理3 a 后土壤盐分显著低于未施加生物炭处理，滴灌处理少量高频的灌水方式使得深层渗漏量减少，地下水无法得到及时补充，单次灌水也未能充分将盐分淋洗到深层土壤，使得在滴灌条件下盐分主要被截留在40～80 cm 土层，而在施加生物炭处理中由于生物炭具有多孔结构及较大表面积等物理特性，使得土壤孔隙度增加，更有益于盐分运移[32]。畦灌处理单次灌水量大、灌水次数少、灌水间隔长，充足的水分使得盐分被淋洗到土壤深层，土壤体积质量增加、孔隙度降低、土壤透气性变差。黄河水畦灌降盐量明显地下水畦灌，这是因为地下水盐分较高，使得盐分不能被充分淋洗到深层土壤。同时畦灌条件下施加生物炭处理较未施加生物炭处理有显著差异，地下水畦灌处理中，施加生物炭处理较表现为脱盐效果，而未施加生物炭处理则表现为积盐效果；黄河水畦灌处理中施加生物炭处理与未施加生物炭处理均表现出脱盐效果，但未施加生物炭处理脱盐量显著高于施加生物炭处理。

表4 不同处理下玉米产量构成要素和水分利用效率 Table4 Yield components and water use efficiency of maize under different treatments

4 结 论

1）滴灌方式下施用生物炭的节水增产效果最为显著，其蒸散量分别比滴灌处理和畦灌处理降低47.3、41.5 mm，水分利用效率分别提高0.2 和0.3 kg/m3，灌溉水利用效率分别提高0.5 和1.1 kg/m3。

2）地下水滴灌+生物炭处理的脱盐效果最优（11.43%），黄河水畦灌+生物炭处理的脱盐效果最差（7.91%），地下水畦灌+生物炭处理的脱盐效果居中（-1.55%）。其中地下水滴灌处理盐分主要集中在膜外40 cm 处，而畦灌处理灌水较为均匀，膜内膜外脱盐效果较为接近。

3）生物炭+地下水滴灌是河套灌区玉米农田优先选择的田间管理措施。