叢 鑫 張立志 徐征和 龐桂斌 馮嚴明 趙登良
(1.濟南大學水利與環境學院,濟南 250022; 2.山東省水利科學研究院,濟南 250014)
小麥作為世界第二大糧食作物,其產量和品質對保障國家糧食安全具有重要作用。灌溉和施肥是影響小麥生長發育和最終產量的兩個相互耦合的重要因子[1-2]。在農業生產中,農民試圖依靠水肥來提高作物產量,但過量的水肥并不能顯著提高作物產量,而且造成水肥利用效率低、經濟效益降低,甚至污染環境[3]。因此,綜合考慮灌溉和氮肥的協同影響、優化灌溉量和施氮量是實現農業可持續發展的關鍵。
研究者通過大量試驗研究了水氮運籌下的作物產量、水肥利用效率與經濟效益效應,探索最佳的水氮配比。薛麗華等[4]在新疆地區研究了水氮耦合對冬小麥的影響,結果表明,灌溉量為3 900 m3/hm2、施氮量為270 kg/hm2是其試驗條件下冬小麥產量最高的水氮高效運籌模式;宋明丹等[5]在陜西地區研究不同水氮水平冬小麥干物質積累特征及產量效應,結果表明,施氮量為210 kg/hm2時冬小麥可以獲得較高的干物質量和產量;周加森等[6]在河北省實施水肥一體化微噴灌下冬小麥田間試驗,推薦減氮20%(減后施氮量192 kg/hm2)下冬小麥水氮利用效率最高;鄭春風等[7]在河南省通過大田試驗研究了減施氮對冬小麥經濟效益的影響,結果表明,將普通尿素習慣施用量減少10%~20%(減后施用量180 kg/hm2)及等量控釋尿素減少20%不僅不會造成冬小麥大幅度減產,反而提高了冬小麥經濟效益。
上述研究只對個別指標進行了水氮模式下的尋優,在最優模式探索中,很難做到平衡水肥效率、產量和經濟效益的多重目標。基于冬小麥生產、經濟效益和環境效益來確定水肥管理多目標優化的研究較為匱乏。本文研究不同灌溉和施氮肥制度對冬小麥產量、水分利用效率、氮肥偏生產力和經濟效益的影響,描述和量化產量、水分利用效率、氮肥偏生產力和經濟效益對灌溉量和施氮量的響應,以冬小麥的多重目標為基礎,通過空間分析確定適宜的灌溉量和施氮量,以期達到高產、節水、節肥的目的。
2018—2019年在山東省濟南市長清區山東省灌溉試驗中心站(36°34′N,116°50′E)利用水肥滲漏研究池進行試驗。滲漏池下的廊道裝有滲漏液收集裝置。研究區域氣象參數每5 min由試驗中心站附近的自動氣象站收集一次,年降雨量654.7 mm,夏季降雨量占全年的62%,年平均潛在蒸散量1 399.70 mm,年平均氣溫14.3℃(2008—2018年),生長季的降雨量與氣溫如圖1所示。研究區土壤質地為壤土,土壤容重在1.40~1.67 g/cm3之間,總孔隙度在48.7%~54.2%之間。種植制度為冬小麥-夏玉米一年兩熟制。每個滲漏池的長、寬、高分別為2.0、3.33、1.2 m,池內土壤取自附近農田。2018年播種前,各池內土壤理化性質基本一致。0~20 cm耕層土壤養分含量為有機質質量比10.89 g/kg,全氮質量比1.07 g/kg,速效磷質量比22.26 mg/kg,速效鉀質量比83.98 mg/kg。
試驗設計3個微噴灌溉水平與3個施氮水平。3個灌溉水平基于當地冬小麥長期計算的參考作物需水量(ETc),設計整個生育期灌溉水量為I1(60%ETc,300.0 mm)、I2(75%ETc,370.0 mm)、I3(ETc,495.0 mm),小麥生育期降雨量100 mm,具體灌溉量見表1。每個滲漏池安裝相同規格的微噴灌管道與水表來控制灌溉量。3個施氮水平在當地農民普施氮肥(330.0 kg/hm2)的基礎上,設計N1(180.0 kg/hm2,優化45%)、N2(255.0 kg/hm2,優化23%)、N3(330.0 kg/hm2)。磷、鉀肥施用量均為150.0 kg/hm2,均底施;氮肥底施50%,追施50%。供試冬小麥品種為“山農28號”,于2018年10月16日,采用播種機播種,播種量130.0 kg/hm2,每個小區種植5行,行距20 cm。冬小麥于2019年6月1日收獲。其他田間管理采用與當地田間相同的標準,以確保小麥在整個生長季生長良好。

表1 冬小麥生育期各階段灌溉量Tab.1 Irrigation for each stage of winter wheat growing seasons mm
1.3.1土壤含水率
在土壤深度為20、40、60、80、100 cm處安裝中子探針來實時監測冬小麥全生育期土壤含水率變化。計算100 cm內的體積含水率,并將其轉換為質量含水率。
1.3.2地上干物質累積量
在小麥每次澆水前后與收獲時,隨機取代表性的植株3株,用剪刀沿土壤表面將冬小麥整個植株剪下,相應裝袋編好序號,帶回實驗室后將所有樣品用清水洗凈、晾干,對應編號裝入紙袋,在105℃干燥箱中殺青30 min,降溫至70℃干燥至質量恒定后用電子天平稱量,即為地上干物質累積量(GB,kg/hm2)。
1.3.3產量與構成因素
收獲時,從每個地塊隨機抽取1個1 m2樣本由人工剪穗。記錄每平方米的有效穗數,然后從中隨機選取20穗,量取有效穗長(EL,cm),手工脫粒得到單穗粒數。最后用固定脫粒機對樣方小麥進行脫粒,脫粒風干后隨機數出1 000粒小麥,以獲得小麥千粒質量。則小麥理論產量為
Y=EPSGTKN×10-2
(1)
式中Y——小麥理論產量,kg/hm2
EP——每平方米的有效穗數,穗/m2
SG——單穗粒數,粒
TKN——千粒質量,g
1.3.4農田耗水量和水分利用效率
土壤蓄水量計算公式為
W=hρω×10/100
(2)
式中W——土壤蓄水量,mm
h——土層深度,cm
ρ——土壤容重,g/cm3
ω——土壤含水率,%
農田耗水量的計算基于土壤水分平衡公式[8],即
ET=P+U+I-D-R-ΔW
(3)
式中ET——農田耗水量,mm
P——降雨量,mm
U——與地下水的交換量,mm
I——灌溉量,mm
D——深層滲漏量,mm
R——地表徑流量,mm
ΔW——生育期前后土壤蓄水量變化量,mm
試驗滲漏池為有底側坑且試驗期間并未產生地表徑流,因此U和R可以忽略,方程簡化為
ET=P+I-D-ΔW
(4)
階段耗水量計算公式為
ETi=Pi+Ii-Di-(SWSi+1-SWSi)
(5)
式中ETi——階段耗水量,mm
Pi——某個生育期降雨量,mm
Ii——某個生育期灌溉量,mm
Di——某個生育期深層滲漏量,mm
SWSi+1——某個生育期末的土壤蓄水量,mm
SWSi——某生育期初的土壤蓄水量,mm
耗水模系數(Ecp)計算公式為
Ecp=ETi/ET×100%
(6)
耗水強度計算公式為
ETd=ETi/d
(7)
式中ETd——耗水強度,mm/d
d——生育階段時間,d
水分利用效率(WUE)計算公式[9]為
WUE=Y/ET
(8)
1.3.5氮肥偏生產力
氮肥偏生產力計算公式為
NPFP=Y/N
(9)
式中NPFP——氮肥偏生產力,kg/kg
N——施氮量,kg/hm2
1.3.6經濟效益
經濟效益計算公式為
EB=YR-IW-NW-OW
(10)
式中EB——經濟效益,元/hm2
YR——冬小麥產量收入,元/hm2
IW——灌溉水分投入,元/hm2
NW——施氮投入,元/hm2
OW——其他投入,元/hm2
采用Excel 2010與SPSS進行數據分析,表中每個指標值均為3個處理的均值,采用ANOVA與LSD(0.05顯著水平)對數據進行方差分析和多重檢驗。采用Matlab進行數據擬合與繪圖。
由圖2(圖中不同小寫字母表示處理間差異顯著(P<0.05),下同)可知,水氮互作下冬小麥ET和ETi之間存在明顯的差異(P<0.05)。I1耗水量在358~369 mm范圍內,I2耗水量在396~406 mm范圍內,I3耗水量在439~450 mm范圍內。說明同一施肥處理下,作物ET隨灌溉量的增加而增加,與I3比較,I1和I2分別節水18.5%和9.8%左右。而同一灌溉處理下,受施肥影響不顯著,這說明在冬小麥生長發育過程中,灌水對作物ET的影響較大。以抽穗-收獲期為例:在N1施肥條件下,I3比I2和I1處理的ET分別高11.21%和48.93%;在N2施肥條件下,I3比I2和I1處理的ET分別高10.50%和58.34%;在N3施肥條件下,I3比I2和I1處理的ET分別高9.87%和64.28%。
從各生育階段看,9種水氮處理耗水量的差異出現在中后期,主要是抽穗-收獲期差異顯著。這是因為在小麥拔節期后開始實施不同灌溉量,可見加大灌溉定額會降低生育前期(播種-拔節期)和中期(拔節-抽穗期)耗水比例,加大生育后期(抽穗-收獲期)耗水比例,這有利于提高作物產量。
表2為各處理3個生育階段的階段Ecp與ETd??梢钥闯?,所有處理的階段Ecp均表現為播種-拔節期最高,抽穗-收獲期最低;所有處理的ETd均表現為拔節-抽穗期最高,播種-拔節期最低。
播種-拔節期,小麥株體小,水分消耗以地表蒸發為主,ETd遠低于生育中期和后期的ETd;但該生育階段占全生育期的70%,Ecp顯著高于其他生育期。在該生育階段只存在氮的處理,不存在水分處理,從表2可看出,氮對Ecp與ETd的影響并不顯著。

表2 不同水肥處理下各生育期耗水模數和耗水強度Tab.2 Ecp and ETd under different nitrogen and irrigation treatments
拔節-抽穗期,小麥進入需水關鍵期,ETd較播種-拔節期顯著增大,I3N2最高可達5.30 mm/d,最低I1N1也有3.93 mm/d;在相同施氮條件下,不同灌溉之間Ecp與ETd表現出顯著差異,I3大于I1和I2,說明該階段是冬小麥灌水施肥關鍵時期,該階段灌溉定額對冬小麥的用水需求具有顯著的影響。
抽穗-收獲期是冬小麥產量形成的關鍵時期,ETd也較大,同樣表現出隨灌溉定額的增加而顯著增加。
冬小麥生育期耗水主要來源于降雨、灌溉水和農田土壤蓄水。從表3可以看到,降雨量占總耗水量的比例隨灌水量的增加逐漸降低,其中I1N3降雨量占總耗水量的比例最高,為27.93%,I3N2處理降雨量占總耗水量的比例最低,為22.26%。灌溉水占總耗水量的比例隨灌水量的增加呈遞增趨勢,I3處理下比例高達89.93%。土壤蓄水消耗量所占比例隨著灌水量的增加呈先增后減的趨勢,而相同灌水量條件下均以N2處理最高,且與N1處理差異達到顯著水平。說明適量的增加灌溉定額與增施氮肥有助于冬小麥對土壤蓄水的吸收利用。滲漏量所占比例絕對值隨灌溉定額的增加逐漸增加,I3較I2和I1所占比例顯著增加,造成水分大量流失。

表3 不同水肥處理下冬小麥耗水結構Tab.3 Composition of water consumption from different water sources under different nitrogen and irrigation treatments
不同生育期下小麥GB在不同灌溉和施氮水平之間存在差異(圖3,圖中ns表示無顯著性差異(P>0.05),*表示顯著性差異(P<0.05),**表示極顯著性差異(P<0.01),下同)。在返青期(圖3a),就氮肥施用量而言,N2和N3下的GB相差不大,但顯著高于N1下的GB,這是因為本試驗是在抽穗后進行水分處理,此時水分對GB未產生顯著性影響。在抽穗期(圖3b),GB隨灌溉水平的變化極顯著。I1和I2處理下的GB與相同氮處理的I3變化相似,但高于I3處理。說明此時灌溉定額過大,影響了干物質的累積。在灌漿期(圖3c),在I1灌溉水平下,N1、N2和N3處理的GB顯著低于其他兩個灌溉水平處理;在I2灌溉水平下,N1、N2和N3的GB差異顯著,由大到小依次為N2、N3、N1。在I3灌溉水平下,N2和N3的GB沒有顯著差異,但顯著高于N1處理(圖3c)。最終,N1、N2和N3在I2下的總GB分別比I1與I3高6.3%、22.9%、15.0%和16.9%、28.6%、22.8%。N2處理較N1使GB增加4.3%~20.6%,N3處理較N1使GB增加5.6%~14.3%。綜上,在一定范圍內,適當提高水氮施用量有助于地上干物質量累積,過量則會產生抑制作用。
灌溉與施氮量水平對所有產量構成因素都有顯著影響(表4),但灌溉與施氮量的交互作用對SG并未構成顯著影響。從產量構成因素來看,灌溉量相同時,EP由大到小均表現為N2、N3、N1,但在I1、I3灌水水平下N1與N3處理中EP未表現出顯著性差異(P>0.05);SG與EL由大到小也表現為N2、N3、N1,I1、I3下的N2與N3處理間差異不顯著,I2下N1與N3處理間差異不顯著;TKN由大到小表現為N2、N3、N1,在3個施氮量處理間有著顯著性差異。3個施氮量水平下,I2下各指標顯著高于其他兩個處理,I1與I3處理間各指標差異不顯著,表明適量灌溉可以提高產量因素,過量灌溉起到抑制作用。

表4 冬小麥產量組分受灌溉、施肥的影響Tab.4 Yield component factors of winter wheat as affected by irrigation and nitrogen
產量(Y)受灌溉與施氮量影響顯著(P<0.01),在同等灌溉量下各施氮處理均隨施氮量的增加先增加后降低,N2處理下產量最高,比N1與N3處理分別高3.6%、2.3%(I1),14.8%、9.3%(I2),12.1%、8.6%(I3),這說明產量與施氮量并不是一直呈現正相關,過量施氮反而會造成產量減少;在相同施氮量下,I2下的產量較I1下產量有明顯的提高,N2下提高比例高達25%,但I3下產量反而較I2有所下降??傮w分析,小麥產量受單一變量影響時,隨灌水量或施肥量的增加,表現為先增加后降低的趨勢,這也與GB的變化一致。因此,只有適當的灌溉和施氮才有利于小麥產量的提高。
冬小麥水分利用效率、氮肥偏生產力受到灌溉與施肥的顯著性影響(圖4、5)。水分利用效率在15.53~20.52 kg/(hm2·mm)之間,I2N2處理最高,I3N1處理最低。在同一灌溉水平下,I1下3種施氮水平下的WUE沒有明顯差異,可能因為灌溉水平偏低無法協同氮肥利用;I2、I3下3種施氮水平的WUE隨著施氮量的增加呈現先增加后降低的趨勢,這說明適量的增施氮肥可以提高產量從而提高WUE。同一施氮水平下,I2和I3由于灌溉使得產量高于I1,但相關研究表明高灌溉使得作物耗水效率提高,這成為限制I3水分利用效率提高的關鍵因素。總體上,相同施氮不同灌水水平下WUE由大到小依次為I2、I1、I3。
氮肥偏生產力范圍為19.81~40.43 kg/kg,I2N1處理最高,I1N3處理最低。在同一灌溉水平下,NPFP隨施氮量的增加而降低,由大到小依次為N1、N2、N3;在同一施氮水平下,不同灌溉下NPFP表現出明顯的差異性,由大到小均為I2、I3、I1,這說明適當的灌溉可以提高氮肥的利用效率。
根據各處理的投入與收益,分析了各處理的經濟效益(表5)。從表5可以看出,I2N2處理經濟效益最高,較同等灌溉量下N1與N3收益分別高2 073.09、1 856.13元/hm2,較同等施氮量下的I1與I3收益分別高3 316.10、2 150.7元/hm2,說明I2N2水肥配施可以提高小麥產最、增加其純收益。I3N3處理的投入最高,水與肥的投入較I2N2高917.5元/hm2,但是水肥投入的提高并沒有提高產量,并導致其經濟效益顯著降低。綜上說明在試驗設計的水肥3個水平下,適量提高水肥投入有利于經濟效益的提高,實際生產中,農民可降低水肥的投入,會帶來更可觀的經濟效益。

表5 冬小麥經濟效益分析Tab.5 Economic benefit analysis of winter wheat as affected by irrigation and nitrogen 元/hm2
根據獲得的數據,分別建立了4個模型來預測產量、水分利用效率、偏肥生產力和經濟效益隨灌溉量和施氮量的變化。圖6為灌溉量和施氮量與冬小麥產量、水分利用效率、氮肥偏生產力和經濟效益的關系,其回歸方程分別為
Y=-17 554.18+90.97I+54.69N+0.003 2IN-
0.111 4I2-0.106 7N2(R2=0.946)
WUE=-23.07+0.154 3I+0.110 4N-
0.000 009 437IN-0.000 204 4I2-
0.000 204N2(R2=0.995)
NPFP=-15.21+0.372 9I-0.114 2N-
0.000 011 91IN-0.000 449 1I2+
0.000 015 70N2(R2=0.930)
EB=-46 441.54+195.14I+116.41N+0.007 049IN-
0.245 2I2-0.234 8N2(R2=0.942)
預測值與實測值間顯著相關,因此,4個模型能夠預測各指標在水氮互作下的變化。
由各指標與水氮施用量建立的二元二次回歸方程計算出各指標最大化下的相應水氮變量值(表6)。Y、WUE、NPFP和EB達到最大值時的灌溉量和施氮量分別為410.0 mm、260.0 kg/hm2,70.0 mm、260.0 kg/hm2,410.0 mm、180.0 kg/hm2,400.0 mm、250.0 kg/hm2。由表6可以看出,NPFP實現最大值時,氮肥的施用量較少,而WUE若想達到最大,灌溉量需較??;水肥投入量偏少時,無法實現Y的最大化與EB的最優化。由此可以看出,生產中很難同時獲得所有指標的最大化。

表6 最大產量、水分利用效率、偏肥生產力和經濟效益及其相應的灌溉量和施氮量Tab.6 Maximum yield, WUE, NPFP, EB and their corresponding irrigation and N fertilizer application rates
進一步分析發現,冬小麥Y、WUE和EB的可接受值分別有90%、85%和80%的重疊,但重疊值與最大值偏離太遠。因此,本研究制定大于等于95%的最大值認為是可以接受的。但NPFP與其他3個指標呈現相反的趨勢,與其他指標的重疊不超過80%??紤]95%置信區間,得出灌溉量為359.8~428.9 mm,施氮量為225.4~280.9 kg/hm2,此時Y、WUE和EB可同時達到其最大值的95%以上。
華北平原是地下水嚴重超采區,冬小麥是華北平原的主要高耗水作物,在當前華北平原缺水嚴重現狀下,探索冬小麥最佳的水肥管理制度,可為華北平原冬小麥節水節肥生產提供理論與技術參考。
從耗水特性來看,冬小麥ETi及ET在在同一施肥處理下隨灌溉定額的增加而明顯增加,受灌水影響顯著,這與彭致功等[10]和栗麗等[11]的研究結果一致。本研究結果表明,隨著灌溉定額的增加,降雨量和土壤蓄水量占冬小麥總耗水量的比例呈降低趨勢,灌溉量和滲漏量占總耗水量的比例呈增加趨勢,與張笑培等[12]的研究結果基本一致。且試驗結果表明,通過減少灌溉量可有效地調節冬小麥對土壤蓄水的吸收利用,而生育前期控水、后期適當地增加灌溉量可以有效地減少冬小麥生育前期的耗水,增加后期耗水,這有助于提高產量,達到節水增產的效果。
小麥產量、水肥利用效率與經濟效益是衡量小麥高產性和高效性的主要指標,適宜的水氮運籌能提高作物水肥利用效率,增加產量與經濟效益。本文研究了灌溉-氮素供給對冬小麥產量、水分利用效率、氮肥偏肥生產力和經濟效益的影響。通過灌溉量-施氮量投入優化模型,確定不同水氮條件下冬小麥產量、水分利用效率、氮肥偏肥生產力和經濟效益綜合效益最大化的灌溉量和施氮量,為冬小麥種植田間管理提供科學依據。
提高產量是灌溉施肥管理的主要目的之一[13]。冬小麥干物質和產量的形成密切相關,本試驗中冬小麥成熟期干物質累積量和產量隨水氮變化相一致,最高的處理均是I2N2,分別達到22 693.5、8 351.29 kg/hm2,較其他處理提高5.5%~41.1%和9.1%~29.5%。因此,適宜的水氮配施更有助于產量的提高,這與李晶晶等[14]研究發現“過量的水與肥都不利于小麥產量的增加”的結果相一致。同時,本研究發現,I2處理較I1處理在抽穗灌漿期的灌溉更為重要,使得干物質累積量與產量有了大幅的提高,但高定額的灌溉(I3)反而降低了干物質積累量與產量,這在秦姍姍等[15]的研究中也得到了證明,抽穗灌漿水對小麥產量有重要作用,當降水不足會嚴重影響冬小麥產量,適當灌水可增加冬小麥產量,過大灌水量無助于小麥產量的提高。
本文研究結果表明,I3下冬小麥水分利用效率較I1、I2水分利用效率明顯降低,這與已有研究[16]結果“在一定的范圍內,過量的灌水會降低作物水分利用效率”相一致。隨著施氮量的增加,在I1條件下,水分利用效率未表現出明顯差異;在I2、I3條件下,水分利用效率隨施氮量的增加先增加后降低。說明在本試驗條件下,一定范圍內增施氮肥能夠提高冬小麥的水分利用效率,這也與文獻[17-18]的缺氮會降低水分利用效率的規律相一致。這是因為當氮肥超過合理范圍時,分配到生殖器官的光合產物比例下降,導致產量下降,進而降低水分利用效率[19]。
氮肥偏生產力隨著施氮量的增加而減小,而不同灌水定額之間由大到小依次為I2、I3、I1,這與文獻[20-21]結論相近。氮肥偏生產力均以I2最高,其他作物也發現了類似的結果,包括冬棗[22]、馬鈴薯[23]和番茄[24]。綜上,適宜的灌溉量才可提高氮肥偏生產力,若灌溉量過多,土壤水分和空氣的比例逐漸變小,影響到作物根系的生長特性,進而影響偏肥生產力、水分利用效率與產量[25]。
諸多研究人員通過多變量回歸和空間分析結合,建立了水和肥料投入與作物產量和水肥利用效率之間的關系[26-27]。研究表明,當灌溉量和施氮量高于最適水平并不能提高各指標值,在某些情況下反而會產生相反的效果。當產量達到最大值,灌溉量和施氮量均高于最佳的灌溉施氮水平,這說明過量的水肥不僅會影響經濟收益,而且會造成水肥淋失污染地下水。因此,確定灌溉施肥策略十分必要。本研究考慮95%置信區間,建立了水氮投入與冬小麥產量、水分利用效率和經濟效益的關系。灌溉量359.8~428.9 mm、施氮量225.4~280.9 kg/hm2能使冬小麥的產量、水分利用效率和經濟效益均達到大于等于最大值的95%。
(1)水氮互作效應對冬小麥耗水規律存在影響,其中灌水對耗水量ET、耗水模數Ecp和耗水強度ETd的影響顯著,施肥對其影響不顯著。耗水量隨灌溉定額的增加而增加,增加灌溉定額會增大作物生育后期耗水比例,有利于提高作物產量;同時,適量增加灌溉定額與增施氮肥有助于冬小麥對土壤蓄水的吸收利用。
(2)地上干物質累積量與產量隨水氮的變化表現出一致性,均隨灌溉量或施氮量的增加呈先增加后降低的趨勢。因此,只有適當的灌溉和施氮才有利于冬小麥最終產量的提高。
(3)冬小麥水分利用效率隨著灌溉量與施氮量的增加呈先增后減的趨勢;氮肥偏生產力隨著施氮量的增加而減小,隨著灌溉量的增加呈先增后降的趨勢;經濟效益表現出隨灌溉量與施氮量的增加先增加、后減小的趨勢。
(4)在試驗條件下,灌溉量為359.8~428.9 mm(72.6%ETc~87.5%ETc)與施氮量為225.4~280.9 kg/hm2的水氮耦合模式可在保證高產量、高經濟效益的同時提高水肥利用效率,從而實現高產、節水、節肥的目標。