基于響應(yīng)面法的多礫石砂土滴灌春小麥水肥條件優(yōu)選

趙經(jīng)華，楊庭瑞，張 恒，虎膽·吐馬爾白，馬 亮，陳凱麗

(1.新疆農(nóng)業(yè)大學水利與土木工程學院，烏魯木齊 830052；2.新疆農(nóng)業(yè)科學院，烏魯木齊 830091)

0 引 言

【研究意義】新疆地處內(nèi)陸干旱半干旱地區(qū)，高蒸發(fā)、少降水、干旱多風，灌溉方式也由漫灌等常規(guī)灌溉方式向滴灌等節(jié)水灌溉方式變化。阿勒泰地區(qū)的多礫石砂土由于其孔隙結(jié)構(gòu)不同于其他類型土壤，保水保肥性較差，灌溉方式及灌水量不同時會明顯影響其土壤中水肥的分布，且種植不同作物時其根系分布情況不同也會影響土壤中的水肥，共同影響著作物的生長發(fā)育。【前人研究進展】小麥是新疆的主要糧食作物[1]。新疆小麥的灌溉方式由漫灌等逐漸演變成現(xiàn)在的滴灌技術(shù)[2]。滴灌為局部灌溉，棵間蒸發(fā)和深層滲漏損失較常規(guī)灌溉明顯降低，大幅減少輸水過程中的渠系滲漏、蒸發(fā)等現(xiàn)象，并把水分精確集中地輸送至作物根系土壤[3]。楊榮贊等[4]研究認為，20 cm左右耕作層垂向土壤水分差異較大，土壤水分隨著繼續(xù)深度加深趨于穩(wěn)定。滴灌小麥根系主要分布在40 cm以上土層，合理的控制滴灌灌水策略可改變根系生長及分布，且利于小麥高產(chǎn)[5, 6]。肥料也是影響作物生長發(fā)育及最終產(chǎn)量的關(guān)鍵因素，而氮肥尤其重要。適宜的施氮量可以促進小麥的生長及增產(chǎn)[7, 8]；并且作物水肥有較明顯的交互作用，在水分虧缺下增加施肥和肥料缺乏下增加灌水可以緩解虧水或缺肥下對作物產(chǎn)生的抑制效應(yīng)[9-11]。【本研究切入點】合理的水肥施用會直接影響作物的干物質(zhì)積累，繼而影響作物產(chǎn)量，干物質(zhì)積累量反映著作物生長狀況的好壞[12]。需以尿素用量為主要施肥變化并結(jié)合不同灌水量對滴灌春小麥進行試驗，研究不同水、肥條件下對春小麥的生長及發(fā)育的影響。【擬解決的關(guān)鍵問題】利用方差分析并結(jié)合響應(yīng)面法優(yōu)選水肥，分析該區(qū)多礫石砂土土質(zhì)下適宜小麥灌水、施肥水平，為春小麥種植管理提供參考。

1 材料與方法

1.1 材 料

春小麥選用新春37號。試驗于2019年在阿勒泰地區(qū)福海縣(87°35′56″E～87°36′01″E，47°00′56″N～47°01′56″N)進行。試區(qū)海拔平均高度445 m；極端最高氣溫為39.30℃，極端最低氣溫為-41.20℃，多年平均氣溫為4.20℃；多年平均降水量為121.90 mm，多年平均蒸發(fā)量為1 820.10 mm。試區(qū)土壤質(zhì)地為多礫石砂土，土壤保水保肥性差。表層40 cm為土層，40 cm以下為砂礫石層，土壤滲透系數(shù)8.78 mm/d；土壤肥力較差，有機質(zhì)和全氮含量為0.21%(2.1 g/kg)和0.027%(0.27 g/kg)，速效氮、速效磷和速效鉀的含量分別為19.50、9.04和92.4 mg/kg，其中速效磷和速效鉀的含量可以達到中等肥力土壤的水平。項目區(qū)土壤全鹽含量(土水比1∶5)為0.161%(1.61 g/kg)，pH值8.52，呈弱堿性，總鹽含量0.64 g/kg。表1

表1 土壤容重和田間持水量Table 1 Soil bulk density and field water holding capacity

1.2 方 法

1.2.1 試驗設(shè)計

播種量為420 kg/hm2。試驗采用1管4行種植方式，小麥行距0.15 m。滴灌帶間距0.6 m，滴頭間距0.25 m，滴頭流量2.6 L/h，額定工作壓力0.1 MPa。小麥生育期劃分：苗期4月22日～5月14日，分蘗期5月15日～5月27日，拔節(jié)期5月28日～6月4日，孕穗期6月5日～6月9日，抽穗揚花期6月10日～6月25日，灌漿期6月26日～7月5日，成熟期7月6日～7月22日。

灌水定額設(shè)置3個水平，分別為30 mm(W1)、45 mm(W2)、60 mm(W3)，在4月22日各處試驗理灌出苗水45mm。施肥設(shè)置3個水平，低(N0)、中(N1)、高(N2)。N0處理不施底肥及氮肥，N1、N2處理底肥種類及用量均相同，僅尿素用量變化。通過施肥量可得N1處理下施氮量為179 kg/hm2、N2處理下施氮量為317 kg/hm2。水肥試驗共設(shè)9個處理，即W1N0、W2N0、W3N0、W1N1、W2N1、W3N1、W1N2、W2N2、W3N2，3次重復，共27個處理。表2，表3

表2 春小麥灌溉制度Table 2 Irrigation system of spring wheat

1.2.2 測定指標

株高：各處理小區(qū)各生育期分別取5株小麥并標記后測量，取均值得株高，小麥株高在抽穗前為地表到最高葉尖的高度，抽穗后為從地表至麥穗的頂端，不連芒。

干物質(zhì)：各處理小區(qū)各生育期內(nèi)取25 cm行長樣本，取地上部分，先在105℃ 烘干30 min，滅活樣本使其停止呼吸消耗，后在80℃下烘干至恒重。

有效穗數(shù)[13]：每穗結(jié)實粒數(shù)在5粒以上的穗數(shù)(被病、蟲危害造成的空穗亦作有效穗計算)。

每穗粒數(shù)[13]：從樣點取20個麥穗混合脫粒，數(shù)清總粒數(shù)，然后用穗數(shù)除總粒數(shù)。

千粒重[13]：每千粒的重量，重復3次取平均值。

產(chǎn)量[13]：收獲曬干揚凈后的實際產(chǎn)量。

表3 春小麥各處理施肥Table 3 Fertilization status of each treatment for spring wheat

1.3 數(shù)據(jù)處理

采用Excel 2019數(shù)據(jù)處理、Design Expert10.0軟件優(yōu)化參數(shù)，Origin 2017繪圖，SPSS 25.0統(tǒng)計學軟件進行方差分析和LSD分析。通過Design Expert10.0軟件進行響應(yīng)曲面分析，利用Design Expert10.0響應(yīng)面軟件中Miscellaneous design(MD)模塊對工藝參數(shù)進行優(yōu)化[14]。

2 結(jié)果與分析

2.1 單因素及多因素方差

2.1.1 不同水肥處理對小麥株高的影響

研究表明，同一施氮水平下，隨著灌水量的增加，小麥株高也隨之增加；但同一灌水水平下，隨著施肥量的增加，小麥株高大致表現(xiàn)為增長趨勢，N0處理與N1、N2處理最后株高有較大差幅，N1、N2之間株高差距不明顯。拔節(jié)孕穗期小麥的株高增長最快，抽穗揚花期以后株高增長趨于緩慢。在拔節(jié)孕穗期， N2施肥水平下株高均與N0施肥水平下株高差異顯著(P<0.05)，與N1施肥水平下株高差異不顯著(P>0.05)。抽穗揚花期W3N0、W2N1、W3N1、W1N2、W2N2、W3N2與W1N0、W2N0株高差異顯著(P<0.05)。在抽穗揚花期至成熟期，N0、N1施肥水平下，W1的小麥株高較小，W1灌水水平不利于小麥的生長。在灌漿期與成熟期W3N1、W3N2的株高均高于其他處理下小麥株高，與W1N0、W2N0、W1N1的株高差異顯著(P<0.05)，但兩者之間差異不顯著(P>0.05)。灌水量與施肥量對株高的影響顯著，少灌水和低施肥處理阻礙了株高的生長，適宜的灌水量及施肥量可以有效地促進小麥株高的增長。表4

表4 各水肥處理下不同生育期春小麥的株高變化Table 4 Plant height of spring wheat at different growth stages under different water and fertilizer treatments(cm)

2.1.2 不同水肥處理對小麥干物質(zhì)積累的影響

研究表明，在同一施肥水平下，小麥生育前期干物質(zhì)積累隨灌水的增加無明顯規(guī)律，在生育后期表現(xiàn)為隨著灌水量的增加，小麥干物質(zhì)積累也增加；在同一灌水水平下，小麥生育前期干物質(zhì)積累隨施肥的增加也無明顯規(guī)律，在生育后期表現(xiàn)為隨施肥量的增加小麥干物質(zhì)積累也增加。小麥全生育期在分蘗期與抽穗揚花期各處理干物質(zhì)積累差異均不顯著(P>0.05)，在灌漿期W3N1、W2N2、W3N2與W1N0、W2N0干物質(zhì)積累差異顯著(P<0.05)。在成熟期W2N2、W3N2與W1N0、W1N1干物質(zhì)積累差異顯著(P<0.05)。小麥的干物質(zhì)積累情況在W1和N0處理下較差，灌水量與施肥量顯著影響著小麥干物質(zhì)積累，干物質(zhì)積累會隨著灌水量或施肥量的增加而增加。表5

表5 各水肥處理下不同生育期春小麥的干物質(zhì)積累變化Table 5 Dry matter accumulation of spring wheat at different growth stages under different water and fertilizer treatments(kg/hm2)

2.1.3 不同水肥處理對小麥產(chǎn)量及產(chǎn)量構(gòu)成的影響

研究表明，灌水對每穗粒數(shù)和產(chǎn)量有極顯著影響(P<0.01)，對有效穗數(shù)有顯著影響(P<0.05)，對千粒重無顯著影響(P>0.05)。施肥對每穗粒數(shù)和產(chǎn)量有極顯著影響(P<0.01)，對千粒重和有效穗數(shù)無顯著影響(P>0.05)。水肥交互作用對產(chǎn)量有極顯著影響(P<0.01)，對產(chǎn)量構(gòu)成要素均無顯著影響(P>0.05)。在N0、N1施肥水平下，千粒重大致表現(xiàn)為隨灌水量增加而增加；在同一灌水水平下，千粒重隨施肥量的增加表現(xiàn)為先增后降的趨勢；W3N1與W1N0千粒重差異顯著(P<0.05)，與其他處理差異不顯著(P>0.05)。在同一施肥水平下，隨灌水量增加每穗粒數(shù)也增加；在同一灌水水平下，隨施肥量的增加每穗粒數(shù)也增加；W3N1、W3N2與W1N0、W2N0、W3N0、W1N1、W1N2的每穗粒數(shù)差異顯著(P<0.05)。在同一施肥水平下，有效穗數(shù)隨灌水量的增加而增加；在同一灌水水平下，有效穗數(shù)隨施肥量的增大呈現(xiàn)先增后降的趨勢；有效穗數(shù)除W3N0外其余處理間差異不顯著(P>0.05)。在同一施肥水平下，產(chǎn)量隨灌水量的增多而增長；在同一灌水水平下，產(chǎn)量隨施肥量的增加呈先增后降的趨勢；W3N1處理與其他處理間產(chǎn)量均差異顯著(P<0.05)。W1N0、W2N0、W3N0處理之間產(chǎn)量兩兩差異顯著，W1N1、W2N1、W3N1處理間兩兩差異顯著， W1N2與W2N2、W3N2差異顯著，在一定情況下灌水可以促進小麥增產(chǎn)。W2N0、W2N1、W2N2處理之間產(chǎn)量兩兩差異顯著(P<0.05)，W3N0、W3N1、W3N2處理間兩兩差異也顯著(P<0.05)，施肥量存在一個閾值，當施肥量超過該閾值時，不會促進小麥增產(chǎn)，甚至產(chǎn)生抑制的效果。在小麥種植過程中要控制灌水與施肥量。表6

表6 不同水肥水平對小麥的產(chǎn)量及產(chǎn)量構(gòu)成方差變化Table 6 Variance analysis of different water and fertilizer levels on wheat yield and yield composition

2.2 響應(yīng)面模型結(jié)果及方差

研究表明，處理9灌水定額為45 mm時，尿素施用300 kg/hm2，小麥產(chǎn)量為6 301.77 kg/hm2，相比高灌水施肥處理7，小麥產(chǎn)量更高；與處理8相比，在同一尿素施用水平下，小麥產(chǎn)量隨灌水定額的增加有所增長，灌水定額為60 mm時，尿素施用300 kg/hm2，小麥產(chǎn)量最優(yōu)為7 040.93 kg/hm2。表7

表7 模型參數(shù)的設(shè)計及運行數(shù)據(jù)Table 7 Design and operation data of model parameters

Y1為產(chǎn)量響應(yīng)值，分析灌水量X1及施肥量X2對小麥產(chǎn)量的影響，以Y1產(chǎn)量為響應(yīng)值建立的二次回歸方程如式(1)所示：

Y1=6 289.37+875.28X1+469.75X2+195.34X1X2-339.25X12-1 131.97X22

(1)

響應(yīng)面產(chǎn)量模型的P<0.05，該模型顯著；但在擬合時，采樣點有時會發(fā)生不全落在響應(yīng)曲面上的情形，響應(yīng)面模擬多元相關(guān)系數(shù)R2為0.954，調(diào)整相關(guān)系數(shù)R2為0.877，該模型的相關(guān)性良好，擬合效果較好；F(X1)>F(X2)，灌水量對小麥產(chǎn)量的影響較顯著。表8

灌水量與施肥量的交互作是顯著的。當灌水量一定時，小麥產(chǎn)量隨著施肥量增加呈先增后降的趨勢，并且在300～400 kg/hm2的施肥量產(chǎn)量達到最大，之后小麥產(chǎn)量有降低趨勢。當施肥量一定時，隨著灌水量的提高小麥產(chǎn)量也逐漸提高，在灌水量60 mm時產(chǎn)量表現(xiàn)較好，水分對小麥生長有較好的促進作用。灌水量與施肥量對小麥產(chǎn)量影響最優(yōu)的為灌水量為60 mm，施肥量為388 kg/hm2，在此灌水施肥水平下，小麥產(chǎn)量為6 923.1 kg/hm2。圖1，圖2

表 8 產(chǎn)量響應(yīng)值回歸方程方差Table 8 Variance analysis of regression equation for yield response value

圖1 灌水量與施氮量處理下小麥產(chǎn)量等高線分布Fig.1 Contour distribution of the effect of irrigation amount and urea amount on wheat yield

3 討 論

小麥生長受到氣候、土壤、水分、營養(yǎng)、管理等影響，其中水分和氮肥對小麥生長獲產(chǎn)是極為重要的。小麥生長過程中會施用含氮、磷、鉀的等多種肥料，其中氮素施用量與小麥生長發(fā)育、籽粒產(chǎn)量和品質(zhì)緊密相關(guān)[15]。試驗發(fā)現(xiàn)在同一灌水水平下，不施肥處理與施肥處理下的小麥株高、干物質(zhì)積累、有效穗數(shù)、千粒重和每穗粒數(shù)均差幅較大，也與不施肥處理中未施底肥有關(guān)；2個施肥處理進行對比發(fā)現(xiàn)，施肥量對小麥株高的影響較小，2個施肥處理生育末期株高差異較小；N2施肥處理的小麥干物質(zhì)積累量相比N1處理的要高；N2施肥處理的小麥千粒重、有效穗數(shù)和產(chǎn)量比N1處理的要低。施肥量影響春小麥的生長發(fā)育，且存在施肥量適宜的臨界閾值。王玲玲等[16]研究認為氮素用量超過一定水平時，小麥氮素利用效率和氮肥生產(chǎn)效率均下降，小麥產(chǎn)量及品質(zhì)也受到影響。叢鑫等[17]也認為小麥灌水與施氮需要適宜的水平，才可以達到較好的經(jīng)濟效益。響應(yīng)面方法對試驗結(jié)果進行模擬分析后，在灌水60 mm，施尿素388 kg/hm2的情況下可取得較好的小麥產(chǎn)量；對試驗結(jié)果進行方差分析后發(fā)現(xiàn)在灌水60 mm，施尿素300 kg/hm2對小麥生長及獲產(chǎn)有較好的影響；試驗所得產(chǎn)量與響應(yīng)面優(yōu)選所得產(chǎn)量兩者相差較小，產(chǎn)量僅只有1%差距，但在施尿素量有一定的變化幅度，可能是由于在小麥種植過程中的肥料施用的不均勻存在隨機性，以及在模擬過程中存在采樣點不全落在響應(yīng)面上的誤差。通常響應(yīng)曲面模型評價結(jié)果的優(yōu)劣是由多元相關(guān)系數(shù)R2來反映，R2是完全擬合的度量值，反映了響應(yīng)面模擬符合給定數(shù)據(jù)的程度，通常要求R2值在0.9以上[18]。多元相關(guān)系數(shù)R2為0.954，模型擬合效果較好；調(diào)整相關(guān)系數(shù)R2為0.877，有87.7%的響應(yīng)值的變化可以通過此試驗方法建立的模型進行闡釋[19]。響應(yīng)面模型評價后模擬效果是較好的。較好的小麥種植效果需要適宜的灌水與施肥水平，并最終將反映在獲產(chǎn)上。響應(yīng)面法在灌溉制度、農(nóng)田施肥等方面優(yōu)化有較好的分析表現(xiàn)，模擬優(yōu)選所得結(jié)論與試驗結(jié)果相近，后續(xù)可利用響應(yīng)面法進行試驗設(shè)計優(yōu)化并優(yōu)選方案。

4 結(jié) 論

4.1不施肥處理相比較施肥處理，小麥生長、產(chǎn)量構(gòu)成及產(chǎn)量狀況較差；在試驗設(shè)計的水肥處理下，千粒重和有效穗數(shù)隨著施肥量的增加表現(xiàn)出了先增后降的趨勢，每穗粒數(shù)為單增趨勢；隨灌水量的增加千粒重、每穗粒數(shù)、有效穗數(shù)也有所增加；僅在小麥生育后期，增加灌水量與施肥量小麥干物質(zhì)積累隨之增加，300 kg/hm2施肥水平、60 mm灌水水平下，春小麥生長發(fā)育狀況較好。

4.2灌水量60 mm，施肥量300 kg/hm2下小麥的干物質(zhì)積累、產(chǎn)量構(gòu)成及產(chǎn)量表現(xiàn)較好，最終產(chǎn)量為7 040.93 kg/hm2；響應(yīng)面法優(yōu)化后的最佳水肥條件為灌水60 mm，施肥388 kg/hm2，該條件下，模擬小麥產(chǎn)量為6 923.1 kg/hm2，與試驗較好水肥條件接近，產(chǎn)量變化僅有1%，該響應(yīng)面模型與試驗結(jié)果相吻合。選取W3N1為較好的水肥處理。