渭北旱塬免耕/深松轮耕麦田产量和土壤水分对施肥的响应模拟
张玉娇1, 李军1,*, 郭正2, 岳志芳2
1西北农林科技大学农学院, 陕西杨凌 712100
2西北农林科技大学林学院, 陕西杨凌 712100
* 通讯作者(Corresponding author): 李军, E-mail:junli@nwsuaf.edu.cn

第一作者联系方式: E-mail:zhangyujiao@nwsuaf.edu.cn, Tel: 15029904879

摘要

为探索不同肥力水平对渭北旱塬连作冬小麦田在长周期免耕/深松轮耕措施下土壤蓄水保墒和作物增产效应的影响, 在模拟精度验证基础上, 应用WinEPIC模型长周期定量模拟研究了1980-2009年渭北旱塬免耕/深松轮耕连作麦田5个不同施肥水平下(T1, N 75 kg hm-2+P2O5 60 kg hm-2; T2, N 120 kg hm-2+P2O5 90 kg hm-2; T3, N 150 kg hm-2+P2O5 120 kg hm-2; T4, N 180 kg hm-2+P2O5 150 kg hm-2; T5, N 255 kg hm-2+P2O5 90 kg hm-2)冬小麦产量和土壤水分效应。在30年模拟期间, 各处理的冬小麦产量、年度耗水量和水分利用效率均呈波动下降趋势, 下降幅度表现为T5>T4>T3>T2>T1。0~5 m土层土壤有效含水量呈季节性波动降低趋势, 且随施肥水平的升高而降低, 5个处理的麦田平均干燥化速率依次为每年13.5、17.1、17.4、20.1和23.9 mm。0~1.5 m土层土壤湿度随季节降水波动; 各处理在不同深度形成稳定的土壤干层, 其中T1在1.5~2.0 m, T2和T3在1.5~3.0 m, T4和T5在1.5~4.0 m。上述结果表明, 随着肥力水平的增加, 旱作冬小麦产量和耗水量也增加, 土壤干层加厚。综合考虑认为, 在渭北旱塬免耕/深松轮耕长期连作小麦田适宜的施肥量为纯氮150 kg hm-2+P2O5 120 kg hm-2

关键词: 渭北旱塬; 冬小麦产量; 施肥水平; 保护性耕作; WinEPIC模型模拟; 土壤水分
Simulating Wheat Yield and Soil Moisture under Alternative No-tillage and Subsoil Tillage in Response to Fertilization Levels in Weibei Highlands
ZHANG Yu-Jiao1, LI Jun1,*, GUO Zheng2, YUE Zhi-Fang2
1College of Agronomy, Northwest A&F University, Yangling 712100, China
2College of Forestry, Northwest A&F University, Yangling 712100, China
Abstract

This study aimed at understanding the responses of grain yield and soil moisture to fertilization level in dryland winter wheat under long-term alternative no-tillage and subsoil tillage. The WinEPIC model was employed on the basis of precision verification with five fertilization levels (T1, N 75 kg ha-1+P2O5 60 kg ha-1; T2, N 120 kg ha-1+P2O5 90 kg ha-1; T3, N 150 kg ha-1+P2O5 120 kg ha-1; T4, N 180 kg ha-1+P2O5 150 kg ha-1; and T5, N 255 kg ha-1+P2O5 90 kg ha-1) and the simulation period was from 1980 to 2009. During the 30-year period, water consumption in the growing season and water use efficiency of winter wheat tended to decrease in a fluctuating manner at different fertilization levels with the ranking sequence of T5>T4>T3>T2>T1. In the 0-5 m soil depth, the monthly available soil moisture tended to decrease in a seasonally fluctuating manner and reduced with the increase of fertilization amount. The soil desiccation rates from T1 to T5 were 13.5, 17.1, 17.4, 20.1, and 23.9 mm per year, respectively. During the simulation period, the soil humidity in 0-1.5 m soil layer fluctuated with the seasonal rainfall. A stable dry soil layer was found under all fertilization levels which was 1.5-2.0 m under T1, 1.5-3.0 m under T2 and T3, and 1.5-4.0 m under T4 and T5. These results indicate that grain yield and water consumption of winter wheat may increase with more fertilizer input, however, the dried soil layer is thickened. In a comprehensive consideration, we suggest N 150 kg ha-1+P2O5 120 kg ha-1 to be the optimal fertilization rates in winter wheat under long-term alternative no-tillage and subsoil tillage in Weibei Highlands.

Keyword: Weibei Highlands; Yield of winter wheat; Fertilization; Conservation tillage; WinEPIC simulation; Soil moisture

渭北旱塬是我国北方典型雨养旱作农区[1, 2], 干旱缺水是当地粮食生产主要限制因素[3]。冬小麦是该地区主要的粮食作物, 产量低而不稳[4, 5]。研究表明, 干旱胁迫不仅限制冬小麦的生长, 影响植株的株高、地上部干重、根干重和总生物量, 还导致冬小麦穗数、千粒重和穗粒数等产量构成因素的显著降低[6]。施肥能够有效增加土壤生产力和有机碳含量, 提高小麦产量[7]。施用氮肥对小麦产量和品质影响显著, 适宜的氮磷钾配比及用量可促进冬小麦对养分的吸收, 形成较大的群体数量和干物质, 提高产量和肥料利用率, 产量随施肥量增加而提高[8, 9]。目前, 关于施肥增产的相关研究较多, 产量随施肥量增加而提高, 但超过一定限度后, 增产不显著, 甚至减产[9, 10, 11], 且施肥量增加会导致产麦田深层土壤水分消耗, 发生持续性土壤干燥化[12]。保护性轮耕将翻耕、深松、免耕等保护性耕作措施合理组合与轮换, 能够克服各项单一土壤耕作措施弊端, 改善耕层土壤结构与通气状况, 有利于充分和持续均衡发挥保持水土和蓄水保墒效应, 促进作物生长发育和提高作物产量[13, 14, 15, 16, 17]。柏炜霞等[18]研究表明, 免耕/深松和深松/翻耕轮耕处理比传统连续翻耕增产9.1%和7.6%, 水分利用效率提高9.6%和11.0%。但由于轮耕试验周期比较长、实施难度较大, 相关土壤轮耕和施肥效应相结合的长期田间试验研究还不多见。

环境政策综合气候模型EPIC (Environmental Policy Integrated Climate)是美国农业部农业研究局研制的农田生产和水土资源管理综合评价动力学模型, 能够长周期连续定量模拟和评价“ 气候— 土壤— 作物— 管理” 综合连续系统响应, 可用来定量评价农田作物生产力和水土资源管理策略的效果[19, 20, 21, 22, 23]。本研究采用的WinEPIC3060模型, 对土壤水分动力学过程描述比较细致, 可以输出逐日分层土壤水分数据, 适用于作物生产系统综合性模拟分析和应用研究, 特别适合于旱地不同田间管理措施下土壤水分生态环境效应的模拟和分析研究。WinEPIC模型能够精确有效地模拟土壤中氮磷运转和不同耕作方式[24, 25], 可较为精确地模拟不同施肥条件下作物轮作的产量和土壤水分状况[26]。在模型中免耕、深松处理主要是通过各土壤耕作处理参数差异体现的, 免耕时不同的土壤耕作深度、秸秆和土壤混合效率、土壤紧实度与深松、翻耕处理等能影响夏闲期土壤耕层结构、土壤水分入渗和水分蒸发。本研究旨在定量分析和评价不同降水年型下不同施肥水平下免耕/深松轮耕麦田的蓄水保墒效果和增产效应, 筛选该旱作麦田最佳施肥量, 为渭北旱塬麦田蓄水保墒和增产增收提供科学依据。

1 WinEPIC模型模拟方法
1.1 参数设置

选择陕西省合阳县为试验点。合阳县位于渭北旱塬东部, 为暖温带半湿润易旱区, 海拔850 m, 无霜期192 d, 年平均气温11.8℃, 年平均降水量571.9 mm, 主要集中在7月至9月, 年际间变化较大。运行WinEPIC模型的必要参数包括当地逐日气象要素、土壤剖面理化特性、作物生长参数和肥料参数等模型数据集。由于缺乏合阳气象站多年逐日气象数据, 本研究借用临近的白水气象站数据。白水县与合阳县直线距离约52 km, 同属于渭北旱塬暖温带半湿润易旱气候区, 农业气象和农业生产条件基本相同。白水站逐日实时气象要素包括1980— 2009年逐日太阳辐射值、最高气温、最低气温、降水量、相对湿度、风速、风向等。根据渭北旱塬冬小麦田间试验测定数据和文献资料, 确定了渭北旱塬冬小麦生理生态参数[18, 22, 23]。合阳县农田代表性土壤为黑垆土, 土层疏松深厚, 土壤容重1.31 g cm-3, 田间持水量25.2%, 饱和含水量平均为33.9%。根据合阳大田试验土壤数据, 土壤剖面理化性状参数包括土壤质地、田间持水量、土壤含水量、容重和氮磷钾等, 参考《陕西土壤》[27]中黑垆土土壤普查数据来确定。本研究根据渭北旱塬农业生产的施肥习惯, 在EPIC模型中肥料数据库添加了尿素和磷酸二铵两种肥料, 并对这两种肥料的参数进行修订, 包括矿质氮和矿质磷等, 建立肥料参数数据库。

1.2 试验设计

试验区为冬小麦连作田, 前茬收获后秸秆全部粉碎还田, 播前采用免耕与深松交替的耕作模式, 其中深松为每隔60 cm宽深松土壤30~35 cm。根据渭北旱塬麦田长期定位试验施肥处理和大田生产施肥水平, 设置5种不同施肥水平, 即T1 (N 75 kg hm-2, P2O5 60 kg hm-2)、T2 (N 120 kg hm-2, P2O590 kg hm-2)、T3 (N 150 kg hm-2, P2O5 120 kg hm-2)、T4 (N 180 kg hm-2, P2O5 150 kg hm-2)、T5 (N 255 kg hm-2, P2O5 90 kg hm-2), 其中氮磷肥分别为尿素和磷酸二铵, 于每年9月25日播种施肥。供试冬小麦品种为晋麦47, 播量为150 kg hm-2

在实时气象条件下, 逐日定量模拟了T1、T2、T3、T4和T5施肥处理下旱作麦田作物生长和产量响应、土壤水分动态变化过程。为了充分反映麦田土壤水分消耗和降水补充平衡状况, 在模型输出的麦田0~5 m土层土壤湿度逐日数据序列中, 选择每月15日麦田0~5 m土层土壤有效含水量代表麦田该月土壤有效含水量特征, 以便比较不同年份和不同季节麦田土壤水分变化动态; 选择每年8月15日麦田土壤湿度剖面分布数据来代表该年麦田土壤湿度剖面分布特征, 比较分析不同耕作处理麦田土壤湿度剖面分布逐年变化特征。

1.3 模型精度验证

WinEPIC模型在长武旱塬麦玉轮作试验中, 不同施肥条件下产量和土壤水分模拟值和观测值相关系数均为0.818~0.886 (P< 0.01) [12, 26]; 在其他地区农田、草地和果园等不同土地利用类型条件下的模拟验证相关系数均在0.78以上(P< 0.05) [28, 29, 30, 31]。WinEPIC模型能够较为准确地模拟大田作物产量和土壤水分变化状况。

利用2008— 2013年合阳试点保护性耕作田间定位试验数据为观测值, 选取每年的产量和每隔3个月的土壤水分数据验证T1、T3和T5施肥处理下免耕/深松轮耕麦田冬小麦产量和0~3 m土层土壤有效含水量模拟值(表1)。比较模拟值和观测值发现, 对小麦产量和土壤有效含水量的模拟效果均令人满意, 模拟值与观测值相关系数均超过0.80, 达到极显著水平(表2)。冬小麦产量模拟值平均值低于其观测值, 相对误差介于-3.12%和-12.03%之间(表2), 模型的模拟精度是误差产生的重要原因; 对0~3 m的土壤有效含水量的吻合度较高(图1), 相对误差介于0.50%和2.60%之间(表2)。

表1 不同施肥水平下免耕/深松麦田0~3 m土层有效含水量模拟值与观测值变化 Table 1 Variation of simulated and observed available soil water amounts in 0-3 m soil layer of no-tillage/subsoil tillage system under different fertilization levels (mm)

图1 不同施肥水平下免耕/深松麦田产量和0~3 m土壤有效含水量模拟值与观测值比较Fig. 1 Comparison of simulated and observed wheat yields and available soil water amounts in 0-3 m soil layer under different fertilization levels in no-tillage/subsoil tillage system

表2 不同施肥水平下免耕/深松麦田冬小麦产量与0~3 m土层有效含水量模拟值与观测值 Table 2 Simulated and observed wheat yields and available soil water amount in 0-3 m soil layer under different fertilization levels in no-tillage/subsoil tillage system
2 结果与分析
2.1 不同施肥处理下免耕/深松轮耕麦田耗水量

30年模拟期内年度降水量平均值为569.8 mm, 变化范围为369.7~836.3 mm, 变异系数为21.9%, 且呈波动性下降趋势(图2-A)。前10年(1980— 1989)平均降水量为628.0 mm, 后10年(2000— 2009)平均降水量为573.3 mm, 后期较前期减少54.7 mm。

图2 模拟期间年度降水量(A)及各处理0~5 m土层耗水量模拟值(B)
年度指小麦生长季, 如1981表示1980年7月至1981年6月期间。土壤耗水量数据来自合阳试验点; 降水量为邻近的白水县数据, 由白水气象站提供。
Fig. 2 Annual precipitation during the simulation period (A) and simulated water consumption in 0-5 m soil layer (B)
Annual refers to wheat growing season, for example, the period from July 1980 to June 1981 is abbreviated with 1981. The water consumption data were from Heyang site and the precipitations were from neighboring Baishui, provided by Baishui weather station. T1: N 75 kg hm-2+P2O5 60 kg hm-2; T2: N 120 kg hm-2+P2O5 90 kg hm-2; T3: N 150 kg hm-2+P2O5 120 kg hm-2; T4: N 180 kg hm-2+P2O5 150 kg hm-2; T5: N 255 kg hm-2+P2O5 90 kg hm-2.

根据冬小麦生产年度降水量(Pi)及多年降水量的平均值( )和均方差(d), 将其分为丰水年( )、平水年( )和干旱年( ) 3种年型[32]。模拟期内3年为平水年, 年度降水量为528.6~611.0 mm; 2012年为丰水年, 年度降水量超过611.0 mm; 2014年为干旱年, 年度降水量低于528.6 mm。

冬小麦生产年度麦田耗水量与上年降水量有密切关系。与平水年相比, 干旱年型麦田耗水量急剧降低, 丰水年型耗水量明显升高(表3)。由于年降水量趋势性降低, 麦田耗水量也随着降水量年度变化呈现波动性降低趋势。干旱年型T5耗水量高于其他施肥水平, 差异显著; 平水年型, 各施肥水平的耗水量差异不显著; 丰水年型降水量明显增加, T4耗水量也明显增加, 与其他施肥水平差异显著。

表3 不同降水年型免耕/深松麦田平均耗水量模拟值比较 Table 3 Comparison of average simulated water consumptions in different rainfall years under no-tillage/subsoil tillage (mm)

随着降水量的减少, 冬小麦生产年度耗水量呈现波动性降低趋势(图2-B)。5个施肥水平间差异显著, T5的平均年耗水量最高, T4次之, T1最低(表4)。在30年模拟期内, 后10年的年度耗水量平均值低于前10年平均值, T1~T5分别降低85.0、84.3、87.3、72.9和76.0 mm。

表4 不同施肥水平0~5 m土层年度耗水量模拟值的变化 Table 4 Variation of simulated water consumption in 0-5 m soil layer during winter wheat growing season under different fertilization levels (mm)
2.2 不同施肥水平下冬小麦产量和水分利用效率

模拟期内5个施肥水平的旱作麦田冬小麦产量均随降水量年际变化呈现剧烈同步波动性变化趋势, 且施肥量越高, 冬小麦产量相应也越高(图3表5)。T2、T3、T4和T5较T1分别增产39.9%、77.7%、88.7%和102.4%, 各施肥水平间产量有显著差异。但是各施肥水平的小麦产量随时间均呈现明显的减产趋势, T1~T5前10年的平均产量分别为2.62、3.55、4.45、4.69和5.01 t hm-2, 而后10年的平均产量分别为1.81、2.56、3.16、3.37和3.60 t hm-2, 减产率依次为31.1%、27.9%、29.1%、28.0%和28.1%。

图3 不同施肥水平下冬小麦产量和水分利用效率模拟值动态Fig. 3 Dynamics of simulated average yield and water use efficiency in winter wheat under different fertilization levels
T1: N 75 kg hm-2+P2O5 60 kg hm-2; T2: N 120 kg hm-2+P2O5 90 kg hm-2; T3: N 150 kg hm-2+P2O5 120 kg hm-2; T4: N 180 kg hm-2+P2O5 150 kg hm-2; T5: N 255 kg hm-2+P2O5 90 kg hm-2.

表5 不同施肥处理下冬小麦产量和水分利用效率模拟值变化 Table 5 Variation of simulated yield and water use efficiency in winter wheat under different fertilization levels

水分利用效率(WUE)是指作物消耗单位水分所形成的经济产量。在模拟期内, T2、T3、T4和T5的WUE分别较T1提高14.2%、27.2%、38.1%和46.8%, 各施肥水平间差异显著, 且施肥量越高WUE增长率也越高(表5)。30年间小麦WUE均呈现波动性降低趋势, T1~T5前10年的WUE平均值分别为8.44、9.61、10.69、11.46和12.24 kg hm-2 mm-1 (图3), 而后10年降为5.25、5.99、6.69、7.37和7.83 kg hm-2 mm-1, 减少37.8%、37.6%、37.5%、35.7%和36.0%。

与平水年型相比, 干旱年型的小麦产量和WUE急剧降低, 而丰水年型明显升高(表6), 反映出小麦产量和WUE与年度降水量密切相关。无论哪种降水年型, T2、T3、T4和T5的小麦产量和WUE均明显高于T1。在30年模拟期间, 小麦产量和WUE均随着施肥量的增加而升高, 5个施肥水平间差异显著, 以T5的产量和WUE最高, T4次之, T1最低。

表6 不同降水年型冬小麦产量和水分利用效率模拟值比较 Table 6 Comparison of simulated wheat yield and water use efficiency (WUE) in different precipitation years
2.3 不同施肥处理下冬小麦水分胁迫和氮素胁迫

干旱胁迫是指由于土壤供水量低于作物需水量对作物生长产生的胁迫, WinEPIC模型通过计算逐日土壤有效含水量和逐日作物需水量判断作物生长是否遭受干旱胁迫。在模拟期内, 不同施肥水平间干旱胁迫日数差异显著, 并随施肥量增加而增加(表7)。不同施肥水平下干旱胁迫日数变化趋势类似(图4), 且与冬小麦生产年度降水量变化呈相反趋势, 即降水量呈逐年趋势性减少, 而干旱胁迫日数呈趋势性增加, T1~T5的平均干旱胁迫日数从前10年的5.5、8.7、10.5、15.6和21.4 d增加到后10年的36.2、41.5、44.8、48.4和56.7 d。

表7 不同施肥水平下冬小麦水分胁迫日数和氮素胁迫日数模拟值变化 Table 7 Variations in simulated water stress days and nitrogen stress days of winter wheat field under different fertilization levels (d)

图4 不同施肥水平下免耕/深松麦田水分胁迫和氮素胁迫天数模拟值动态Fig. 4 Dynamics of simulated water stress days and nitrogen stress days in different fertilization levels under no-tillage/subsoil tillage system
T1: N 75 kg hm-2+P2O5 60 kg hm-2; T2: N 120 kg hm-2+P2O5 90 kg hm-2; T3: N 150 kg hm-2+P2O5 120 kg hm-2; T4: N 180 kg hm-2+P2O5 150 kg hm-2; T5: N 255 kg hm-2+P2O5 90 kg hm-2.

随着作物持续高强度耗水作用, 土壤水库贮水量逐渐降低, 在模拟中后期阶段土壤供水量难以满足作物需水量, 干旱胁迫逐渐显现并逐年加剧。施肥量越高, 氮素胁迫(由于土壤氮素供应不足对作物生长产生的胁迫)日数越大。WinEPIC模型中通过计算逐日作物氮素需求量和逐日土壤氮素供给量判断作物生长是否受到氮素胁迫。在30年模拟期内, 不同施肥水平间冬小麦生育期氮素胁迫日数差异显著(表7), 且变化趋势类似(图4), 均呈现波动性上升趋势, 其主要原因是随着土壤氮素逐年消耗, 在模拟中后期采用全部化肥播种期基施的施肥方式, 至小麦生长旺盛期土壤供氮量无法满足作物需氮量, 氮素胁迫逐渐显现并呈现波动性加剧趋势。随着施肥量增加, 氮素胁迫程度逐渐降低, 增加施氮量能明显延迟氮素胁迫出现的时间, 降低氮素胁迫程度。

2.4 不同施肥处理免耕/深松麦田逐月土壤有效含水量

土壤有效含水量是土壤总含水量减去凋萎湿度含水量后可供作物吸收利用的有效水分含量。不同施肥水平下0~5 m土层逐月有效含水量呈现季节性波动变化趋势(图5)。在30年模拟时段内, T2、T3、T4和T5的土壤有效含水量明显低于T1 (表8), T2和T3的土壤有效含水量变化曲线接近重合(图5)。

表8 不同施肥水平下0~5 m土层逐月土壤有效含水量模拟值变化 Table 8 Variation of simulated monthly available soil water amount in 0-5 m soil layer under different fertilization levels (mm)

图5 不同施肥水平的麦田0~5 m土层逐月土壤有效含水量模拟值动态Fig. 5 Dynamics of simulated monthly available soil water amount in 0-5 m soil layer under different fertilization levels

渭北旱塬降水集中在7月至9月, 正值冬小麦的休闲期, 此时0~5 m土壤有效含水量多出现高峰。在模拟初期, 当地年降水量较高, 初始土壤水分和养分含量均较高, 因而模拟期第1~第2年(1980— 1981)的0~5 m逐月土壤有效含水量随着作物生长而剧烈波动性降低, 施肥水平对土壤有效含水量影响不大; 随着模拟年限延长, 土壤养分和水分亏缺现象逐渐显现, 且年降水量趋势性降低, 施肥对麦田土壤水分影响逐渐增大, 模拟中前期5个施肥水平间差异逐渐显现, 土壤有效含水量随着施肥量的增加而降低, 0~5 m土层土壤有效含水量表现为T1> T2> T3> T4> T5, 不同施肥水平变化趋势相似, 其中T2和T3施肥处理逐月有效含水量变化曲线十分接近趋向于重合(图5)。在模拟第1年, 5个施肥水平的0~5 m土壤有效含水量均为936.4 mm, 但模拟结束时(2009年6月), T1~T5的土壤含水量分别为543.9、441.1、433.1、353.1和244.3 mm, 比模拟开始时平均每年减少13.5、17.1、17.4、20.1和23.9 mm。T1、T2、T3和T4的土壤有效含水量明显高于T5, 且T5的土壤干燥化速度明显快于其他施肥水平, T1的土壤干燥化速度最慢, 麦田蓄水保墒效果较好。

5个施肥水平下, 以T1施肥处理夏闲期蓄墒效果好, T2和T3次之, T5最差(表8)。不同降水年型夏闲期土壤有效贮水量有明显变化, 在干旱年型各施肥处理夏闲期土壤有效贮水量表现为T1> T2和T3> T4> T5, 说明低肥有利于夏闲期麦田蓄水保墒。

2.5 不同施肥处理麦田土壤湿度剖面分布年度变化

在模拟初期(1981— 1987), 0~5 m土层土壤湿度逐年降低, 由模拟初期、中期(1991— 1997)到末期(2001— 2007), 0~5 m土壤湿度也呈现递减趋势, 各施肥水平下表现一致(图6)。0~5 m土壤湿度均随着施肥水平提高和模拟年限延长, 土壤干层逐年加深和加厚, 之后形成稳定的土壤干层。模拟初期麦田0~5 m土壤湿度随着年度降水而变化, T1、T2、T3、T4和T5施肥水平下, 土壤水分利用最大深度分别达到2 m (第4~第5年)、3 m (第5~第6年)、3 m (第5~第6年)、3 m (第7~第8年)和5 m (第8~第9年); 到模拟中后期, 5个施肥水平下, 土壤湿度长期保持在稳定低湿状态(0.085~0.089 m m-1左右)的土层深度分别是1.5~2.0、1.5~3.0、1.5~3.0、1.5~4.0和1.5~ 4.0 m。可见, 随模拟时间的延长深层土壤湿度也保持稳定状态, 已不随降水发生年度变化, 只有0~1.5 m土壤湿度随降水发生年度变化, 尤以距地表1 m内的土壤湿度变化剧烈, 在降水量较高的雨季或严重干旱年份, 1.0~1.5 m土层的土壤湿度也会发生明显变化。在30年模拟期内, 随着肥力水平的提高和降水量的趋势性降低, 麦田土壤干燥化程度逐渐加剧, 土壤干层逐渐加厚。

图6 模拟初、中、后期不同施肥水平下0~5 m土壤剖面分布特征变化Fig. 6 Dynamics of soil moisture distribution in 0-5 m soil profile under different fertilization levels during initial, metaphase, and telophase simulations

3 讨论

氮素作为最主要的作物大量营养元素, 对提高作物产量和改善品质具有重要作用[33]。轮耕有利于改善耕层土壤结构与通气状况, 促进作物生长发育[17], 起到增产的作用。本研究在1980— 2009年模拟研究期间, T1~T5肥力水平下冬小麦平均产量分别为1.99、2.78、3.53、3.75和4.03 t hm-2, WUE平均值分别为6.82、7.79、8.68、9.42和10.01 kg hm-2 mm-1, 随着施肥量增加, 免耕/深松麦田产量和水分利用效率随之升高, 这与王学春等[23]研究结论相似。免耕/深松能够有效地提高麦田夏闲期土壤蓄水能力, 故其产量较单一翻耕方式高。肥力越高, 作物耗水量越大, 土壤有效含水量越低[30]。在本研究30年模拟研究期间, 冬小麦生产年度耗水量因年降水量变化而有明显波动性降低趋势, T1~T5的0~5 m土层生长季耗水量平均值依次为577.0、582.7、583.0、599.2和608.8 mm。耗水量随着施肥量的增加而升高, 这一结果与在晋中半干旱区的研究结果[30]不同, 这可能是免耕/深松轮耕模式的蓄水效果及与施肥量的互作效应带来的影响。由于30年模拟期间冬小麦生产年度降水量和耗水量均呈现波动性减少趋势, 导致冬小麦产量和WUE模拟值均呈现波动性降低趋势。

在秸秆还田或覆盖条件下, 免耕、翻耕和深松等不同保护性耕作措施各有优缺点, 免耕和深松有利于减缓土壤侵蚀和蓄水保墒, 可使麦田降水蓄存率由翻耕的25%~35%提高到50%~65%, 土壤贮水量增加60~120 mm [34, 35, 36, 37, 38, 39, 40]。研究表明, 2年免耕1年深松能保持较高的入渗率, 增加入渗量和入渗速度, 增强土壤的水分入渗能力, 使土壤的蓄水、保墒能力增强[41, 42]。在本研究中, 土壤有效含水量随着肥力水平的增加而降低, 0~5 m土层逐月土壤有效含水量表现为T1> T2> T3> T4> T5, 进一步印证随着施肥量的增加, 作物耗水量增多, 土壤有效贮水量降低[38]。无论在整个冬小麦生产年度还是夏闲期, 麦田土壤贮水量均随着施肥量的增加而减少。王学春等[12]研究表明, 低肥麦田的土壤有效含水量比高肥麦田高100~200 mm, 土壤干燥化速率慢(每年46.7 mm)。我们通过研究30年0~5 m土层土壤湿度剖面分布变化动态, 发现0~1.0 m土层土壤湿度年度变化剧烈; 在降水量较高的雨季或严重干旱年份, 1.0~1.5 m土层土壤湿度也会发生明显变化, 5个施肥处理均在模拟过程中发生了土壤湿度降低和土壤干燥化程度加剧现象, 形成一定程度的稳定土壤干层。T1为1.5~2.0 m土层, T2和T3为1.5~3.0 m土层, T4和T5为1.5~4.0 m土层。随着施肥水平提高, 土壤干层逐年加深加厚[30]。降雨量减少和作物耗水量增加是土壤干层形成的主要原因, 合理控制农田施肥量是控制农田土壤干层形成的有效途径[37, 38]

关于麦田土壤贮水量、产量和WUE模拟值趋势性降低效应。在初期阶段, 麦田初始土壤含水量较高, 且同期年降水量也较高, 所以麦田土壤有效含水量较高; 但随着麦田肥力水平的提高, 麦田的耗水量增加以及作物多年生产耗水, 导致麦田土壤干燥化逐渐显现和加剧, 深层土壤供水能力逐渐削弱, 土壤干层加厚。小麦生长则主要依靠当季降水供给, 并且伴随自然降水趋势性减少导致麦田耗水量降低, 由于干旱胁迫程度逐渐加剧, 最终导致麦田产量和WUE剧烈波动性降低。本研究中, 各施肥水平下0~5 m土层逐月土壤有效含水量呈现季节波动性减少趋势, 30年间T1~T5的土壤有效含水量分别减少392.5、495.3、503.3、583.3和629.1 mm, 平均每年减少13.5、17.1、17.4、20.1和23.9 mm, 相当于麦田每年减少一场中等强度降水量。随着麦田土壤有效贮水量下降和土壤供水年际调剂能力削弱, 麦田耗水量趋势性降低导致冬小麦产量和WUE同步趋势性下降, 且产量和WUE年际波动性加剧, 这也是旱作高产麦田土壤水分过耗导致的不良土壤水分环境效应。

本研究采用WinEPIC模型长周期定量模拟不同施肥水平下免耕/深松麦田蓄水保墒和作物产量变化动态, 弥补了传统田间试验中试验周期长、试验难度较大等方面的不足, 能够较为快速、准确地获得不同施肥条件下渭北旱塬旱作麦田较长时间序列下免耕/深松轮耕麦田的小麦生产力和水肥利用规律。在本研究模型精度验证所采用的大田试验施肥处理中包含了钾肥, 但WinEPIC模型中并不包含钾素模块, 模拟结果对钾肥施用量并没有响应, 因而模拟试验施肥处理中没有涉及钾肥相关内容。由于黄土高原土壤富钾, 钾素尚不成为作物生长发育的限制因素。今后需要进一步修订EPIC模型, 补充和完善钾素模块, 以便详细描述氮磷钾肥互作效应及其对作物生长和土壤水分利用的影响。

4 结论

通过30年的模型模拟研究, 发现在免耕/深松轮耕模式下, 施肥量越高小麦产量和WUE越高, 作物耗水量越大, 导致土壤有效含水量降低, 土壤干燥化程度严重, 土壤干层加厚, 最终作物产量和WUE增幅缩小, 施肥节水和增产报酬递减。在5个施肥处理中, T3的土壤有效含水量和土壤干层与T2接近, 深层土壤水分消耗较为适度, 但小麦产量增幅明显, 节水和增产效果最好。从施肥对作物增产影响和土壤水分可持续利用角度统筹考虑, 渭北旱塬免耕/深松轮耕麦田较为适宜的施肥水平为纯氮150 kg hm-2、P2O5 120 kg hm-2

The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.

参考文献
[1] 杨文治, 邵明安. 黄土高原土壤水分研究. 北京: 科学出版社, 2000. pp 10-45
Yang W Z, Shao M A. Research of Soil Water on Loess Plateau. Beijing: Science Press, 2000. pp 10-45(in Chinese) [本文引用:1]
[2] 中华人民共和国国家统计局. 中国统计年鉴2012http: //www. stats. gov. cn/
National Bureau of Statistics of China. China Statistical Yearbook 2012 (http://www.stats.gov.cn(in Chinese) [本文引用:1]
[3] 李军, 邵明安, 张兴昌. 黄土高原旱塬地冬小麦水分生产潜力与土壤水分动态的模拟研究. 自然资源学报, 2004, 19: 738-746
Li J, Shao M A, Zhang X C. Simulation of water potential productivity of winter wheat and soil water dynamics on rainfed highland of the Loess Plateau. J Nat Res, 2004, 19: 738-746 (in Chinese with English abstract) [本文引用:1]
[4] 高焕文, 李问盈, 李洪文. 中国特色保护性耕作技术. 农业工程学报, 2003, 19(3): 1-4
Gao H W, Li W Y, Li H W. Conservation tillage technology with Chinese characteristics. Trans CSAE, 2003, 19(3): 1-4 (in Chinese with English abstract) [本文引用:1]
[5] 李军, 王立祥, 邵明安, 樊廷录. 黄土高原地区小麦生产潜力模拟研究. 自然资源学报, 2001, 16: 161-165
Li J, Wang L X, Shao M A, Fan T L. Simulation of wheat potential productivity on Loess Plateau region of China. J Nat Resour, 2001, 16: 161-165 (in Chinese with English abstract) [本文引用:1]
[6] 曹阳, 杨婕, 熊伟, 武永锋, 冯灵芝, 杨晓光. 1962-2010年潜在干旱对中国冬小麦产量影响的模拟分析. 农业工程学报, 2014, 30(7): 128-139
Cao Y, Yang J, Xiong W, Wu Y F, Feng L Z, Yang X G. Simulation of winter wheat yield influenced by potential drought in China during 1962-2010 Trans CSAE, 2014, 30(7): 128-139 (in Chinese with English abstract) [本文引用:1]
[7] 高会议, 郭胜利, 刘文兆, 车升国. 黄土旱塬区冬小麦不同施肥处理的土壤呼吸及土壤碳动态. 生态学报, 2009, 29: 2551-2559
Gao H Y, Guo S L, Liu W Z, Che S G. Soil respiration and carbon fractions in winter wheat cropping system under fertilization practices in arid-highland of the Loess Plateau. Acta Ecol Sin, 2009, 29: 2551-2559 (in Chinese with English abstract) [本文引用:1]
[8] 赵鹏, 陈阜. 秸秆还田配施化学氮肥对冬小麦氮效率和产量的影响. 作物学报, 2008, 34: 1014-1018
Zhao P, Chen F. Effects of straw mulching plus nitrogen fertilizer on nitrogen efficiency and grain yield in winter wheat. Acta Agron Sin, 2008, 34: 1014-1018 (in Chinese with English abstract) [本文引用:1]
[9] 陈毓君, 陈黎岭, 同延安, 杨宪龙, 王少杰, 马海洋, 梁连友. 不同施肥水平对冬小麦群体动态和产量形成的影响. 西北农业学报, 2012, 21(4): 47-53
Chen Y J, Chen L L, Tong Y A, Yang X L, Wang S J, Ma H Y, Liang L Y. Effect of different fertilizer levels on population dynamics and grain yield formation of winter wheat. Acta Agric Boreali-Occident Sin, 2012, 21(4): 47-53 (in Chinese with English abstract) [本文引用:2]
[10] 郝明德, 来璐, 王改玲. 黄土高原塬区旱地长期施肥对小麦产量的影响. 应用生态学报, 2003, 14: 1893-1896
Hao M D, Lai L, Wang G L. Effect of yield wheat with long-term fertilization in table dry-land of the Loess Plateau. Chin J Appl Ecol, 2003, 14: 1893-1896 (in Chinese with English abstract) [本文引用:1]
[11] 张少民, 郝明德, 陈磊. 黄土高原长期施肥对小麦产量及土壤肥力的影响. 干旱地区农业研究, 2006, 11(6): 85-89
Zhang S M, Hao M D, Chen L. Effects of long-term fertilization on yield of wheat and soil fertility in dry-land of Loess Plateau. Agric Res Arid Areas, 2006, 11(6): 85-89 (in Chinese with English abstract) [本文引用:1]
[12] 王学春, 李军, 郝明德. 施肥水平对冬小麦田土壤水分影响的模拟. 农业工程学报, 2008, 24(9): 38-43
Wang X C, Li J, Hao M D. Effects of fertilization levels on soil water in winter wheat field. Trans CSAE, 2008, 24(9): 38-43 (in Chinese with English abstract) [本文引用:3]
[13] 江晓东, 李增嘉, 侯连涛, 王芸, 王雪, 颜红. 少免耕对灌溉农田冬小麦/夏玉米作物水、肥利用的影响. 农业工程学报, 2005, 21(7): 20-24
Jiang X D, Li Z J, Hou L T, Wang Y, Wang X, Yan H. Impacts of minimum tillage and no-tillage systems on soil NO3--N content and water use efficiency of winter wheat/summer corn cultivation. Trans CSAE, 2005, 21(7): 20-24 (in Chinese with English abstract) [本文引用:1]
[14] Blevins R L, Frye W W, Baldwin P L, Robertson S D. Tillage effects on sediment and soluble nutrient losses from a maury silt loam soil. J Environ Qual, 1990, 19: 683-686 [本文引用:1]
[15] Edwards W M, Shipitalo M J, Traina S J, Edwards C A, Owens L B. Role of Lumbricus terrestris (L. ) burrows on quality of infiltrating water. Soil Biol Biochem, 1992, 24: 1555-1561 [本文引用:1]
[16] 李娟, 李军, 尚金霞, 贾志宽. 轮耕对渭北旱塬春玉米田土壤理化性状和产量的影响. 中国生态农业学报, 2012, 20: 867-873
Li J, Li J, Shang J X, Jia Z K. Effects of rotational tillage on soil physicochemical properties and spring maize yield in Weibei Highland s. Chin J Eco-Agric, 2012, 20: 867-873 (in Chinese with English abstract) [本文引用:1]
[17] 孙国锋, 陈阜, 肖小平, 伍芬琳, 张海林. 轮耕对土壤物理性状及水稻产量影响的初步研究. 农业工程学报, 2007, 23(12): 109-113
Sun G F, Chen F, Xiao X P, Wu F L, Zhang H L. Preliminary study on effects of rotational tillage on soil physical properties and rice yield. Trans CSAE, 2007, 23(12): 109-113 (in Chinese with English abstract) [本文引用:2]
[18] 柏炜霞, 李军, 王玉玲, 王丽. 渭北旱塬小麦玉米轮作区不同耕作方式对土壤水分和作物产量的影响. 中国农业科学, 2014, 47: 880-894
Bai W X, Li J, Wang Y L, Wang L. Effects of different tillage methods on soil water and crop yield of winter wheat -spring maize rotation region in Weibei Highland . Sci Agric Sin, 2014, 47: 880-894 (in Chinese with English abstract) [本文引用:2]
[19] 李保国, 龚元石, 左强. 农田土壤水分的动态模型及应用. 北京: 科学出版社, 2000. pp 25-55
Li B G, Gong Y S, Zuo Q. A dynamic Model of Soil Moisture and Applications in the Farmland . Beijing: Science Press, 2000. pp 25-55(in Chinese) [本文引用:1]
[20] 李军, 邵明安, 张兴昌. EPIC模型中农田水分运移与利用的数学模拟. 干旱地区农业研究, 2004, 22(2): 72-75
Li J, Shao M A, Zhang X C. Simulation equations for soil water transfer and use in the EPIC model. Agric Res Arid Areas, 2004, 22(2): 72-75 (in Chinese with English abstract) [本文引用:1]
[21] 李军, 邵明安, 张兴昌, 李世清. 黄土高原旱塬区高产玉米田土壤干燥化与产量波动趋势模拟研究. 中国生态农业学报, 2007, 15(2): 54-58
Li J, Shao M A, Zhang X C, Li S Q. Simulation of soil desiccation and yield fluctuation of high yield maize field on rain-fed highland of the Loess Plateau. Chin J Eco-Agric, 2007, 15(2): 54-58 (in Chinese with English abstract) [本文引用:1]
[22] 李军, 邵明安, 张兴昌, Williams J R. EPIC模型中作物生长与产量形成的数学模拟. 西北农林科技大学学报(自然科学版), 2004, 32(增刊): 25-30
Li J, Shao M A, Zhang X C, Williams J R. Simulation equations for crop growth and yield formation in the EPIC model. J Northwest Sci-Tech Univ Agric For (Nat Sci Edn), 2004, 32(suppl): 25-30 (in Chinese with English abstract) [本文引用:2]
[23] 王学春, 李军, 郝明德. 施肥水平对长武旱塬地冬小麦产量影响的模拟. 农业工程学报, 2008, 24(8): 45-50
Wang X C, Li J, Hao M D. Simulation of fertilization effect on winter wheat yield in Changwu dry highland . Trans CSAE, 2008, 24(8): 45-50 (in Chinese with English abstract) [本文引用:3]
[24] 李军, 邵明安, 张兴昌. EPIC模型中土壤氮磷运转和作物营养的数学模拟. 植物营养与肥料学报, 2005, 11: 166-173
Li J, Shao M A, Zhang X C. Simulation equations for soil nitrogen and phosphorus transfer and crop nutrition in the EPIC model. Plant Nutr Fert Sci, 2005, 11: 166-173 (in Chinese with English abstract) [本文引用:1]
[25] 张玉娇, 李军, 郭正, 岳志芳. 渭北旱塬小麦产量和土壤水分对保护性耕作的响应模拟. 应用生态学报, 2015, 26: 800-808
Zhang Y J, Li J, Guo X, Yue Z F. Simulated responses of winter wheat yield and soil moisture to different conservation tillage practices in Weibei Highland s, Northwest China. Chin J Appl Ecol, 2015, 26: 800-808 (in Chinese with English abstract) [本文引用:1]
[26] 王学春, 李军, 樊廷录. 黄土旱塬不同施肥水平下小麦玉米轮作的产量和土壤水分效应模拟研究. 植物营养与肥料学报, 2008, 14: 242-251
Wang X C, Li J, Fan T L. Modeling the effects of winter wheat and spring maize rotation under different fertilization treatments on yield and soil water in rain-fed highland of Loess Plateau. Plant Nutr Fert Sci, 2008, 14: 242-251 (in Chinese with English abstract) [本文引用:2]
[27] 陕西省土壤普查办公室. 陕西土壤. 北京: 科学出版社, 1992. pp 20-50
Shaanxi Soil Census Office. Shaanxi Soil. Beijing: Science Press, 1992. pp 20-50(in Chinese) [本文引用:1]
[28] 王亚莉, 李军, 王学春, 张社红. 陇东旱塬旱作苹果园水分生产力与土壤干燥化效应模拟. 西北农林科技大学学报(自然科学版), 2011, 39(8): 131-140
Wang Y L, Li J, Wang X C, Zhang S H. Simulation of water productivity and soil desiccation effects of dryland apple orchard on the plateau of eastern Gansu. J Northwest Sci-Tech Univ Agric For (Nat Sci Edn), 2011, 39(8): 131-140 (in Chinese with English abstract) [本文引用:1]
[29] 张社红, 李军, 王学春, 王亚莉. 渭北旱塬苹果园地产量和深层土壤水分效应模拟. 生态学报, 2011, 31: 3767-3777
Zhang S H, Li J, Wang X C, Wang Y L. Modeling the changes of yield and deep soil water in apple orchards in Weibei rainfed highland . Acta Ecol Sin, 2011, 31: 3767-3777 (in Chinese with English abstract) [本文引用:1]
[30] 王学春, 李军, 胡伟, 蒋斌. 黄土高原晋中半干旱区不同肥力水平下连作春玉米产量与干燥化效应模拟. 干旱地区农业研究, 2008, 26(1): 1-11
Wang X C, Li J, Hu W, Jiang B. Simulation of and yield soil desiccation effects of continuous spring maize under different level fertilization treatments in Jinzhong semiarid areas of areas of the Loess Plateau. Agric Res Arid Areas, 2008, 26(1): 1-11 (in Chinese with English abstract) [本文引用:4]
[31] 王学春, 李军, 郝明德. 长武旱塬草粮轮作田土壤水分可持续利用模式模拟. 农业工程学报, 2011, 27(增刊1): 257-266
Wang X C, Li J, Hao M D. Simulation of sustainable use of soil water in dry land for alfalfa-grain rotation system at Changwu arid-plateau of China. Trans CSAE, 2011, 27(suppl-1): 257-266 (in Chinese with English abstract) [本文引用:1]
[32] 陶林威, 马洪, 葛芬莉. 陕西省降水特性分析. 陕西气象, 2000, (5): 6-9
Tao L W, Ma H, Ge F L. Analysis on the characteristics of precipitation in Shaanxi province. Shaanxi Meteorol, 2000, (5): 6-9 [本文引用:1]
[33] 顾炽明, 郑险峰, 黄婷苗, 侯仰毅, 王朝辉. 秸秆还田配施氮肥对冬小麦产量及氮素调控的影响. 干旱地区农业研究, 2013, 31(5): 48-53
Gu C M, Zheng X F, Huang T M, Hou Y Y, Wang C H. Effects of straw returning combined with nitrogen fertilizer application on yield of winter wheat and nitrogen regulation. Agric Res Arid Areas, 2013, 31(5): 48-53 (in Chinese with English abstract) [本文引用:1]
[34] 李玲玲, 黄高宝, 张仁陟, 晋小军, Li G, Chan K Y. 不同保护性耕作措施对旱作农田土壤水分的影响. 生态学报, 2005, 25: 2326-2332
Li L L, Huang G B, Zhang R Z, Jin X J, Li G, Chan K Y. Effects of conservation tillage on soil water regimes in rainfed areas. Acta Ecol Sin, 2005, 25: 2326-2332 (in Chinese with English abstract) [本文引用:1]
[35] 付国占, 李潮海, 王俊忠, 王振林, 曹鸿鸣, 焦念元, 陈明灿. 残茬覆盖与耕作方式对土壤性状及夏玉米水分利用效率的影响. 农业工程学报, 2005, 21(1): 52-56
Fu G Z, Li C H, Wang J Z, Wang Z L, Cao H M, Jiao N Y, Chen M C. Effects of stubble mulch and tillage managements on soil physical properties and water use efficiency of summer maize. Trans CSAE, 2005, 21(1): 52-56 (in Chinese with English abstract) [本文引用:1]
[36] 方日尧, 同延安, 赵二龙, 梁东丽. 渭北旱塬不同保护性耕作方式水肥增产效应研究. 干旱地区农业研究, 2003, 21(1): 54-57
Fang R Y, Tong Y A, Zhao E L, Liang D L. Effect of conservation tillage on moisture, fertility and yield in Weibei Highland . Agric Res Arid Areas, 2003, 21(1): 54-57 (in Chinese with English abstract) [本文引用:1]
[37] 彭文英, 张雅彬. 免耕对粮食产量及经济效益的影响评述. 干旱地区农业研究, 2006, 24(4): 114-118
Peng W Y, Zhang Y B. Review of impacts of no-tillage on crop yield and economic benefit. Agric Res Arid Areas, 2006, 24(4): 114-118 (in Chinese with English abstract) [本文引用:2]
[38] 蒋斌, 李军, 王学春, 胡伟. 不同施肥水平下宝鸡旱塬连作冬小麦产量与深层土壤水分效应的模拟研究. 西北农林科技大学学报(自然科学版), 2008, 36(8): 73-83
Jiang B, Li J, Wang X C, Hu W. Simulation of yield and deep soil water effects of winter wheat field under different fertilization levels in Baoji dry-land of the Loess Plateau. J Northwest Sci-Tech Univ Agric For (Nat Sci Edn), 2008, 36(8): 73-83 (in Chinese with English abstract) [本文引用:3]
[39] 陈洪松, 邵明安, 王克林. 黄土区深层土壤干燥化与土壤水分循环特征. 生态学报, 2005, 25: 2491-2497
Chen H S, Shao M A, Wang K L. Desiccation of deep soil layer and soil water cycle characteristics on the Loess Plateau. Acta Ecol Sin, 2005, 25: 2491-2497 (in Chinese with English abstract) [本文引用:1]
[40] 樊军, 郝明德, 邵明安. 黄土高原沟壑区草地土壤深层干燥化与氮素消耗. 自然资源学报, 2004, 19: 201-206
Fan J, Hao M D, Shao M A. Soil desiccation and nitrogen consumption of artificial meadow in the Loess Plateau. J Nat Resour, 2004, 19: 201-206 (in Chinese with English abstract) [本文引用:1]
[41] 侯贤清, 贾志宽, 韩清芳, 孙红霞, 王维, 聂俊峰, 杨宝平. 不同轮耕模式对旱地土壤结构及入渗蓄水特性的影响. 农业工程学报, 2012, 28(5): 85-94
Hou X Q, Jia Z K, Han Q F, Sun H X, Wang W, Nie J F, Yang B P. Effects of different rotational tillage patterns on soil structure, infiltration and water storage characteristics in dryland . Trans CSAE, 2012, 28(5): 85-94 (in Chinese with English abstract) [本文引用:1]
[42] Qin H L, Gao W S, Ma Y C, Ma L, Yin C M, Chen Z, Chen C L. Effects of subsoiling on soil moisture under no-tillage for two years. Agric Sci China, 2008, 7: 88-95 [本文引用:1]