趙偉霞，李久生，楊汝苗，栗巖峰

?

基于土壤水分空間變異的變量灌溉作物產量及節(jié)水效果

趙偉霞，李久生※，楊汝苗，栗巖峰

（中國水利水電科學研究院流域水循環(huán)模擬與調控國家重點實驗室，北京 100048）

提高整個田塊作物生長指標和產量的均勻性是實施變量灌溉水分管理的目標之一。該研究基于土壤可利用水量（available water holding capacity，AWC）將試驗區(qū)劃分為4個水分管理區(qū)，利用相同的灌水控制指標（0.45AWC）進行分區(qū)變量灌溉水分管理；作為對照，基于最小AWC區(qū)的土壤水分進行均一灌溉水分管理。對比變量灌溉和均一灌溉條件下冬小麥、夏玉米生長指標（株高、葉面積指數、地上部分干物質質量）、葉片相對葉綠素含量、產量及其均勻性，分析AWC對作物生長和產量的影響。結果表明，與均一灌溉相比，夏玉米變量灌溉節(jié)水14.1%，冬小麥灌水量相同。與均一灌溉相比，變量灌溉對冬小麥、夏玉米生長指標、葉片相對葉綠素含量和產量的影響均未達到顯著水平，而不同AWC管理區(qū)之間作物生長指標和產量的差異均達到了顯著水平。為獲得更高的作物產量，建議不同AWC管理區(qū)內采用不同的灌水控制指標。研究可為大型噴灌機變量灌溉水分管理決策提供依據。

土壤；灌溉；作物；土壤可利用水量；管理分區(qū)；冬小麥；夏玉米；均勻性

0 引 言

精準灌溉是精準農業(yè)的重要組成部分，將效益-成本比最大化是精準灌溉管理的主要目標之一。在較大農田尺度內，精準灌溉管理的實現主要通過實時、實地的提供適宜的灌溉水量，最大限度地優(yōu)化灌水量，以獲得最高產量和最大經濟效益，同時保護生態(tài)環(huán)境[1-3]。大型噴灌機變量灌溉技術為精準灌溉管理提供了有力工具[3]，而管理區(qū)劃分隨之成為決定沿噴灌機行走方向和沿噴灌機桁架方向實施變量灌水的前提。

土壤可利用水量（available water holding capacity，AWC）指田間持水率與凋萎含水率的差值，反映土壤的實際供水能力。因此，基于AWC的靜態(tài)分區(qū)是實施變量灌溉水分管理的一種重要方法[4-6]。與傳統(tǒng)的均一灌溉管理相比，獲得更均勻的作物生長和產量是評估變量灌溉效果的重要指標之一[6]。然而，近期研究表明[7-9]，不同土壤結構和AWC均會影響作物水分利用、產量和冠層溫度。對這種觀點的驗證有助于提高變量灌溉水分管理水平。

因此，本文以傳統(tǒng)的均一灌溉水分管理為對照，在不同管理區(qū)內基于相同的灌水指標進行變量灌溉水分管理，分析變量灌溉對冬小麥和夏玉米生長、產量及其空間均勻性的影響，并研究管理區(qū)之間作物生長指標和產量的差異，為提高大型噴灌機變量灌溉技術水分管理水平提供參考。

1 材料與方法

1.1 試驗區(qū)概況

試驗在河北省涿州市東城坊鎮(zhèn)中國農業(yè)大學教學試驗場進行。該試驗場位于河北省中部（39°27′N，115°51′E，海拔42 m），屬暖溫帶半濕潤季風區(qū)，大陸性季風氣候特點明顯。多年平均降雨量563.3 mm，年平均溫度11.6 ℃，多年平均凍土深度0.4 m，最長連續(xù)凍結122 d。試驗場地質構造屬太行山山洪沖積扇。冬小麥-夏玉米輪作是該地區(qū)主要的種植模式。

1.2 圓形噴灌機變量灌溉系統(tǒng)

試驗用圓形噴灌機變量灌溉系統(tǒng)由中國水利水電科學研究院研制，該系統(tǒng)通過對原有均一噴灑式圓形噴灌機的改造獲得，主要由電磁閥、壓力調節(jié)器、地緣識別器和控制系統(tǒng)組成（圖1），變量灌水通過調節(jié)電磁閥開閉狀態(tài)和噴灌機行走速度實現[10]。圓形噴灌機中心支軸距末端旋轉式噴頭的距離為139.93 m，末端噴頭的射程為6 m。根據國家標準對噴灌機“有效半徑”的定義[11]，該變量灌溉系統(tǒng)控制的圓形面積為6.55 hm2。試驗區(qū)為灌溉區(qū)域的1/4，面積1.64 hm2。

圖1 圓形噴灌機變量灌溉系統(tǒng)

1.3 變量灌溉管理分區(qū)

基于AWC對試驗區(qū)進行管理分區(qū)。確定AWC的直接方法是通過田間試驗[12]或室內試驗[13]測量不同土層的田間持水量和凋萎含水量。但是直接測量法具有勞動強度大、耗時長和投資高的缺點。因此，Rawls等[14]通過對大量文獻和數據的分析，提出了基于土壤粒級組成的田間持水率和凋萎含水率間接估算方法，并在變量灌溉管理分區(qū)中應用[6]。為獲得AWC空間分布圖，將試驗地塊按照12 m×12 m的網格劃分為110個單元（圖2），在每個單元內分別以網格中心為圓心，在半徑1 m的圓周上用土鉆取4個剖面，分層混合后用MS2000激光粒度儀（英國馬爾文儀器有限公司）測量0～0.2、＞0.2～0.4、＞0.4～0.6、＞0.6～0.8 m土層的土壤砂粒、粉粒、黏粒體積百分數，并根據美國農業(yè)部推薦方法確定土壤質地[15]。土壤質地及其剖面分布統(tǒng)計結果如表1所示。不同土層的砂粒、粉粒和黏粒含量變異系數（）、土壤質地復雜程度均隨土壤深度增加呈增大趨勢，且黏粒變異系數大于砂粒。

圖2 試驗區(qū)土壤可利用水量面狀分布圖

為獲得所有單元不同土層的田間持水量，共選取13個網格用環(huán)刀法進行剖面取土，并采用威爾克斯法[13]測定13個網格中心不同土層（0～0.2、＞0.2～0.4、＞0.4～0.6、＞0.6～0.8 m）的田間持水量和容重。以田間持水量為因變量（，cm3/cm3）、以黏粒（1，%）、粉粒（2，%）、砂粒體積百分數（3，%）為自變量利用二次多項式逐步回歸法擬合方程式（1）。

=2068.22-312.861-45.843+2.9212-0.2122+

0.2532+2.721·2+3.171·3+0.0452·3

2=0.72，=0.001 （1）

表1 試驗區(qū)土壤顆粒體積分數及質地沿剖面分布

將其他單元不同土層的土壤粒級組成代入式（1）即可獲得所有單元的田間持水量。為獲得所有單元不同土層的凋萎含水量，利用ROSETTA軟件[16]和砂粒、粉粒、黏粒體積百分數分別擬合110個取樣點0～0.2、＞0.2～0.4、＞0.4～0.6 m土層的van Genuchten模型方程（式（2）），將凋萎含水率對應的土壤水勢=1500 kPa[17]代入式（2）即可計算凋萎含水量。

()=θ+(θ-θ)/[1+(·)]1-1/n（2）

式中()為土壤體積含水量，cm3/cm3；θ為土壤飽和含水量，cm3/cm3；為土壤殘余含水量，cm3/cm3；（cm-1）、為待定的特性參數。

利用上述方法獲得的田間持水量和凋萎含水量分別計算所有單元0～0.2、＞0.2～0.4、＞0.4～0.6 m土層的AWC，取0～0.6 m土層的AWC平均值代表作物主要根區(qū)的AWC，然后利用ArcGIS軟件中的克里格插值法和“Jenks”自然斷點法根據作物主要根區(qū)的AWC生成AWC面狀分布圖[6]，如圖2所示。由圖2可知，試驗區(qū)內AWC變化范圍較大，表明在該試驗區(qū)進行變量灌溉管理試驗的可行性。

根據AWC分類等級將試驗區(qū)劃分為4個管理區(qū)，1~4區(qū)AWC變化范圍分別為152.2～161.4、＞161.4～170.9、＞170.9～185.1、＞185.1～204.7 mm。因為噴灌機第一跨控制面積較小，不做變量灌溉試驗處理[18]。每個管理區(qū)對應的面積分別為0.33、0.76、0.18和0.27 hm2（圖3）。

注：VRI和URI分別為變量灌溉和均一灌溉，下同。

1.4 試驗設計及過程

在管理區(qū)1、2、3、4內，按照面積近似相等原則分別設置變量灌溉（V）和均一灌溉（U）2個灌水處理，分別記為1V、1U、2V、2U、3V、3U、4V和4U，共8個試驗處理。在每個處理內，將每個12 m×12 m的網格作為1個重復，1V、1U、2V、2U、3V、3U、4V和4U處理的重復數分別為6、6、17、19、3、3、4和4。

供試作物為冬小麥（濟麥22）和夏玉米（強生928）。冬小麥2013年10月16日播種，行距15 cm，播種量375 kg/hm2，2014年6月12日收獲，收獲時秸稈同時粉碎并覆蓋在土壤表面。夏玉米2014年6月13日播種，株距26 cm，行距55 cm，10月2日收獲。

根據土壤含水量下限值確定灌水時間。結合FAO-56關于冬小麥、夏玉米的推薦值[19]和試驗區(qū)土壤含沙量較大的特點，土壤水消耗比率取為0.45，即當土壤水分虧缺量達到0.45AWC時進行灌溉。對于變量灌溉的4個管理區(qū)（1V、2V、3V和4V），進行分區(qū)灌溉管理，即當各處理內0～0.6 m土層的土壤水分虧缺量達到0.45AWC時該區(qū)進行灌溉。對于均一灌溉的4個管理區(qū)（1U、2U、3U和4U），灌溉進行統(tǒng)一管理，即當最小AWC區(qū)（1U）0～0.6 m土層的土壤水分虧缺量達到該區(qū)0.45AWC時4個管理區(qū)同時用噴灌機變量灌溉系統(tǒng)進行灌溉，以保證所有管理區(qū)均不缺水，以期獲得最高產量[4,20]。變量灌溉處理和均一灌溉處理采用相同的灌水定額，拔節(jié)期結束前為20 mm，拔節(jié)期結束后為30 mm。為防止凍傷，所有處理均于2013年11月18日進行冬灌，灌水量42 mm。

為判斷灌水時間并觀測作物生育期內的土壤含水量動態(tài)變化，在所有處理的每個重復內分別埋設1根Trime-T3探管（德國IMKO公司），共53根，埋設位置在網格中心（圖3），埋設深度0.8 m，每隔7 d測量1次。同時，在1U、1V、2V、3V和4V處理分別埋設1組EM50土壤水分傳感器（美國Decagon Devices公司），埋設深度為0.12，0.24，0.36，0.48和0.6 m，每隔30 min采集1次土壤含水量數據，通過Trime-T3與EM50數據的對比精確判斷灌水時間。因為4V和4U區(qū)土壤為南水北調輸水管線鋪設后的回填土，0.4 m土層以下為較大礫石，Trime-T3探管和EM50土壤水分傳感器埋設深度為0.4 m，灌水時間由0～0.4 m土層的土壤含水量進行判斷。

所有處理采用相同的施肥管理。根據測土配方施肥結果，冬小麥施N量為175 kg/hm2，P2O5為138 kg/hm2，K2O為90 kg/hm2。其中，54 kg/hm2N和全部的P2O5、K2O均以底肥施入土壤，所用肥料為磷酸二銨（N、P2O5質量分數分別為18%、46%）和氯化鉀（K2O質量百分數為60%），利用播種機播種的同時進行施肥；剩余施N量分別在冬小麥返青期和拔節(jié)期用尿素進行追施，追肥比例為返青期:拔節(jié)期=1:1.75，施肥方式為撒施。夏玉米施N量為228.6 kg/hm2，P2O5為55.2 kg/hm2，K2O為63.0 kg/hm2，其中21.6 kg/hm2N和全部的P2O5、K2O均以底肥施入土壤，所用肥料類型和施肥方式與冬小麥相同；剩余施N量分別在夏玉米拔節(jié)期和大喇叭口期用尿素進行追施，追肥比例為拔節(jié)期：大喇叭口期=1:2。

1.5 指標測定與計算

在冬小麥和夏玉米生育期內，測量的作物生長指標主要包括株高、葉面積指數（leaf area index，LAI）和地上部分干物質質量，測量的生理指標為葉片相對葉綠素含量（soil and plant analyzer development，SPAD）。測量開始前，在每個處理的每個重復內隨機選取長勢均勻的冬小麥4株和夏玉米3株進行掛牌標記。利用鋼卷尺和葉綠素儀（SPAD502，日本柯尼卡美能達公司）分別對標記植株的株高和SPAD值進行測量，每隔10 d測量1次。測量SPAD時，冬小麥生長前期測定新的完全展開葉，后期測定旗葉；夏玉米苗期至抽雄前測定最上部完全展開葉，抽雄后測定穗位葉。冬小麥LAI采用比葉重法[21]測量，然后將所有樣品烘干后測量地上部分干物質質量，每隔16 d左右測量1次。夏玉米LAI采用長寬系數法通過對標記植株所有綠葉的長寬測量進行計算，系數為0.75。夏玉米地上部分干物質質量測量時，在每個處理的所有重復內隨機選取3個重復作為取樣點，然后在每個取樣點內分別選取代表性植株3株，烘干后測量。

冬小麥和夏玉米收獲前進行人工采樣測產。測量冬小麥產量時，分別在每個處理的每個重復內選取具有代表性的長勢均勻的1 m×1.05 m的樣方1個，人工脫粒后烘干計產。夏玉米產量測量時，樣方為2行×1.56 m的區(qū)域，人工脫粒后烘干計產。

1.6 數據分析

冬小麥、夏玉米生長指標（株高、LAI、地上部分干物質質量）、SPAD和產量的空間均勻性利用變異系數進行定量描述，為樣本數值的標準差相對于平均數大小的相對量。為了分析變量灌溉水分管理和AWC對作物生長和產量的影響，利用SAS軟件中的PROC GLM方差分析過程進行組內觀察值數目不等的單向分組資料的方差分析[22]。

2 結果與分析

2.1 均一灌溉和變量灌溉處理灌溉制度差異

冬小麥、夏玉米生育期內不同灌水處理的累積灌水量和降雨量如表2所示。冬小麥生育期內總降雨量61.2 mm，小于冬小麥生育期內的多年平均降雨量126 mm，屬于干旱年。夏玉米生育期內總降雨量311.6 mm，小于有秸稈覆蓋條件下太行山山前平原夏玉米總耗水量350～360 mm[23]，說明試驗年份冬小麥和夏玉米均需要進行灌溉。

為提高肥料利用效率，冬小麥和夏玉米進行追肥后在沒有降雨時均進行灌溉。4月3日和4月13日施冬小麥分蘗肥和拔節(jié)肥時，各處理分別灌水20和16 mm；在7月15日施夏玉米拔節(jié)期肥時各處理分別灌水20 mm，7月29日的抽穗肥隨降雨施入。冬小麥和夏玉米的土壤含水量分別在播種后6個月和1.5個月達到灌水下限。從4月22日到6月5日內灌水，加上冬灌水和施肥水，冬小麥每個處理整個生育期總灌水量248 mm，均一灌溉處理和所有變量灌溉管理區(qū)的總灌水量相等。但是由于各處理土壤持水能力不同，土壤含水量達到灌水下限的時間不等，最大相差3 d，3V處理先進行灌水，然后是均一灌溉處理進行灌水。夏玉米均一灌溉總灌水量為70 mm，與變量灌溉2V和4V處理的灌水量相等，但灌水時期不盡相同；1V和3V處理只灌水2次，總灌水量40 mm。與冬小麥相比，由于降雨量增加，夏玉米各處理土壤含水量達到灌水下限的時間差異減小，最大相差1 d。按各處理面積與其總灌水量的乘積計算每個處理的灌水量，與均一灌溉管理相比，夏玉米生育期內變量灌溉管理共節(jié)水14.1%。

表2 冬小麥和夏玉米不同灌水處理的灌水量

注：1V～4V分別表示1～4區(qū)的變量灌溉處理；U表示均一灌溉；下同。

Note: 1V-4V are variable rate irrigation management in zone 1, 2, 3, and 4, respectively; U is uniform rate irrigation management; Same as below.

2.2 作物生育期內土壤含水量動態(tài)變化

冬小麥和夏玉米生育期內0～0.4和＞0.4～0.6 m土層各處理平均土壤含水量動態(tài)變化如圖4所示。由于4區(qū)0.4 m以下土壤含水量無法測定，圖4b和圖4d只給出3個區(qū)土壤含水量測定結果。冬小麥生育期內，0～0.4和＞0.4～0.6 m土層的土壤含水量變化趨勢相似，在冬小麥快速生長期，各處理土壤含水量較高且基本相等，進入生長中期實施變量灌溉水分管理后，各處理土壤含水量差異增大，除成熟期因為降雨使土壤含水量出現了明顯的增大外，各處理土壤含水量整體上處于下降趨勢。在所有均一灌溉處理（1U～4U）內，雖然灌水時間相同，但因為土壤持水能力不同，3U處理的土壤含水量一直處于最低水平，4U處理的土壤含水量則處于最高水平。進行變量灌溉水分管理后，1V和3V處理的土壤含水量分別較1U和3U處理有所提高，2V和4V處理則分別較2U和4U處理有所降低。

夏玉米生育期內由于降雨頻繁，0～0.4 m和＞0.4～0.6 m土層的土壤含水量均處于較高水平；0～0.4 m土層土壤含水量變化幅度較＞0.4～0.6 m土層大。進行變量灌溉水分管理后，因為灌水量相對較少，所以由土壤持水能力不同造成的各管理區(qū)土壤含水量差異雖然存在，但均一灌溉和變量灌溉處理之間的土壤含水量差異較冬小麥明顯減小。說明作物生育期內土壤剖面的含水量主要由其持水能力決定，變量灌溉水分管理對土壤剖面土壤含水量的調節(jié)作用與作物生育期內的降水量有關。

2.3 變量灌溉水分管理對作物生長及其均勻性影響

冬小麥、夏玉米生育期內均一灌溉和變量灌溉處理的作物株高、LAI、SPAD、地上部分干物質質量隨時間的變化曲線如圖5所示。冬小麥各處理作物生長過程相似，株高隨時間快速增加，并于5月13日達到最大值（56 cm），然后由于穗發(fā)育而稍微降低；從拔節(jié)期到抽穗期，冬小麥LAI呈緩慢增加趨勢（5.7～6.3），在成熟期快速減小；地上部分干物質質量在冬小麥整個生育期內一直呈增加趨勢；SPAD除了在成熟期因葉片變黃而降低外，其他時間基本保持恒定（52～61）。夏玉米株高、地上部分干物質質量和SPAD的變化過程與冬小麥相似，但LAI從拔節(jié)期到抽穗期呈快速增加趨勢（1.1～4.6）。由冬小麥和夏玉米生育期內作物株高、LAI、SPAD、地上部分干物質質量在所有測量日期的平均值方差分析結果（表3）可知，均一灌溉和變量灌溉處理的冬小麥、夏玉米生長指標和生理指標均未產生顯著差異。

由冬小麥、夏玉米生育期內作物株高、LAI、SPAD、地上部分干物質質量在所有測量日期的平均值（表3）可知，冬小麥和夏玉米株高的值均小于0.1，呈弱變異[16]，表現出較好的空間均勻性。冬小麥SPAD的值稍大于0.1，呈中等變異，而夏玉米SPAD的值小于0.1，呈弱變異。與株高和SPAD相比，LAI和地上部分干物質質量的值相對較大，空間均勻性較差，且夏玉米LAI和地上部分干物質質量的值小于冬小麥。這表明作物生長空間均勻性與作物類型有關，夏玉米生育期內由于較多的降雨減輕了作物生長對灌溉的響應，表現出比冬小麥更好的空間均勻性。與均一灌溉管理相比，變量灌溉處理的冬小麥株高、LAI、地上部分干物質質量和SPAD的值分別增加了0.02、0.04、0.04和0.02，而變量灌溉處理的夏玉米株高、LAI和SPAD的值與均一灌溉處理相等，地上部分干物質質量的值降低了0.05。與作物生長指標和生理指標的方差分析結果相似，均一灌溉和變量灌溉處理的冬小麥、夏玉米產量均無顯著差異（表3）。由值可知，夏玉米產量變異系數為0.16～0.18，小于冬小麥的0.23～0.27。與均一灌溉處理相比，冬小麥值增加了0.04，而夏玉米值降低了0.02。

表3 均一灌溉和變量灌溉處理冬小麥和夏玉米生長指標均值和變異系數

注：不同小寫字母表示在0.05水平上處理間差異顯著，下同。

Note: Different letters after values indicate significant difference among treatments at 0.05 level; same as below.

由上述結果可知，與均一灌溉相比，變量灌溉對作物生長指標、產量及其空間均勻性的影響較小，增產效益不明顯，這與其他研究者的結果[5,20,24-28]相似。因此，為了充分發(fā)揮變量灌溉技術潛在的節(jié)水、增產優(yōu)勢，需進一步提高變量灌溉技術水分管理水平。

2.4 土壤可利用水量對作物生長的影響

因為均一灌溉處理和變量灌溉處理的作物生長指標和產量之間無顯著差異（表3），因此，將同一AWC管理區(qū)所有樣本在所有測量日期的測量結果取平均值后分析AWC對作物生長和產量的影響（表4）。在冬小麥生育期內，除SPAD外，與不同AWC管理區(qū)之間產量的顯著差異結果相一致，不同管理區(qū)之間的作物株高、LAI和地上部分干物質質量均存在顯著差異，且最小的冬小麥株高、LAI和地上部分干物質質量均在產量最小的1區(qū)和3區(qū)獲得。這與李久生等[29]關于干旱區(qū)春小麥的研究結果相同，即AWC的離散程度對小麥產量有明顯影響。在夏玉米生育期內，除不同管理區(qū)之間的作物株高產生了顯著差異外，LAI、地上部分干物質質量、SPAD和產量均無顯著差異。但是在2015年夏玉米生育期內的降雨量為286 mm時，不同管理區(qū)之間的夏玉米產量差異達到了顯著水平[30]。說明同一塊田內土壤持水能力差異對作物生長和產量的影響依賴于作物生育期內的降雨量。

表4 不同變量灌溉管理區(qū)間冬小麥和夏玉米生長指標和產量

3 結論與討論

基于AWC進行管理分區(qū)，并利用圓形噴灌機變量灌溉系統(tǒng)進行變量灌溉水分管理后，以傳統(tǒng)的均一水分管理為對照，對比分析了變量灌溉水分管理對作物株高、LAI、SPAD、地上部分干物質質量、產量，及其空間均勻性的影響，并研究AWC對作物生長和產量的影響。主要結論如下：

1）當不同AWC管理區(qū)采用相同的土壤水分下限值時，與均一灌溉相比，變量灌溉水分管理的節(jié)水效果與作物生育期內的降雨量有關。在本試驗的氣象條件下，冬小麥生育期內變量灌溉與均一灌溉處理灌水量相等，夏玉米生育期內變量灌溉處理節(jié)水14.1%。

2）變量灌溉水分管理對作物生長和產量均勻性的影響較小，與均一灌溉相比，冬小麥變異系數值增加了0.02～0.04，夏玉米除地上部分干物質質量和產量分別降低了0.05和0.02外，其他指標的變異系數值相等。

3）與均一灌溉相比，變量灌溉水分管理對冬小麥和夏玉米株高、LAI、SPAD、地上部分干物質質量和產量的影響均未達到顯著水平。不同AWC管理區(qū)之間的作物生長指標和產量差異達到了顯著水平，且這種差異在作物生育期內降雨量較少時更為明顯。說明不同AWC管理區(qū)具有不同的水分生產潛力，為了充分發(fā)揮變量灌溉節(jié)水增產效益，需要進一步研究不同管理區(qū)適宜的灌水下限值。

[1] 汪懋華．“精細農業(yè)”發(fā)展與工程技術創(chuàng)新[J]．農業(yè)工程學報，1999，15(1)：l－8． Wang Maohua. Development of precision agriculture and innovation of engineering technologies[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 1999, 15(1): 1－8. (in Chinese with English abstract)

[2] 趙春江，薛緒掌，王秀，等. 精準農業(yè)技術體系的研究進展與展望[J]. 農業(yè)工程學報，2003，19(4)：7－12. Zhao Chunjiang, Xue Xuzhang, Wang Xiu, et al. Advance and prospects of precision agriculture technology system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2003, 19(4): 7－12. (in Chinese with English abstract)

[3] O′Shaughness S A, Rush C. Precision agriculture: Irrigation [M]//Encyclopedia of Agriculture and Food Systems. Manhattan, New York: Academic Press, 2014: 521－535.

[4] Ritchie J T, Amato M. Field evaluation of plant extractable soil water for irrigation scheduling[J]. Acta Horticulturae, 1990, 278: 595－615.

[5] Nijbroke R, Hoogenboom G, Jones W J. Optimizing irrigation management for a spatially variable soybean field[J]. Agricultural Systems, 2003, 76(1): 359－377.

[6] King B A, Stark J C, Wall R W. Comparison of site-specific and conventional uniform irrigation management for potatoes[J]. Applied Engineering in Agriculture, 2006, 22(5): 677－688.

[7] Duncan H A. Locating the Variability of Soil Water Holding Capacity and Understanding its Effects on Deficit Irrigation and Cotton Lint Yield[D]. Tennessee: University of Tennessee, 2012.

[8] Tolk J A, Evett S R. Do soil textural properties affect water use efficiency? [EB/OL]. (2015-07-26) [2016-01-26]. http:// elibrary.asabe.org/abstract. asp?aid=46432.

[9] DeJonge K C, Taghvaeian S, Trout T J. Soil variability effects on canopy temperature in a limited irrigation environment [EB/OL]. (2015-07-26) [2016-01-26]. http://elibrary.asabe.org/ abstract.asp? aid=46425.

[10] 趙偉霞，李久生，楊汝苗，等．田間試驗評估圓形噴灌機變量灌溉系統(tǒng)水量分布特性[J]．農業(yè)工程學報，2014，30(22)：53－62． Zhao Weixia, Li Jiusheng, Yang Rumiao, et al. Field evaluation of water distribution characteristics of variable rate center pivot irrigation system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(22): 53－62. (in Chinese with English abstract)

[11] 農業(yè)灌溉設備中心支軸式和平移式噴灌機水量分布均勻度的測定：GB/T 19797-2005[S].北京：中國標準出版社，2005.

[12] Cassel D K, Nielsen D R. Chapter 36: Field capacity and water availability[M]//Klute A. Methods of Soil Analysis Part 1: Physical and Mineralogical Methods. Madison, Wis: ASA, 1986: 902－926.

[13] Wilcox J C. Rate of soil drainage following an irrigation Ⅲ: A new concept of the upper limit of available moisture[J]. Canadian Journal of Soil Science, 1962, 122－128.

[14] Rawls W J, Brakensiek D L, Saxtonn K E. Estimation of soil water properties[J]. Transactions of the ASAE, 1982, 25 (5): 1316－1320.

[15] Shirazi M A, Boersma L. A unifying quantitative analysis of soil texture[J]. Soil Science Society of America Journal, 1984, 48: 142－147.

[16] Schaap M G, Leij F J, van Genuchten M T. Rosetta: a computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions[J]. Journal of Hydrology, 2001, 251(3/4): 163－176.

[17] 雷志棟，楊詩秀，謝森傳. 土壤水動力學[M]. 北京：清華大學出版社，1987.

[18] Dukes M D, Perry C. Uniformity testing of a variable rate center pivot irrigation control systems[J]. Precision Agriculture, 2006, 7(3): 205－218.

[19] Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration: Guidelines for computing crop water requirements[M]//FAO Irrigation and Drainage Paper 56. Rome, Italy: FAO, 1998.

[20] Hedley C B, Bradbury S, Ekanayake J. Spatial irrigation scheduling for variable rate irrigation[J]. Proceedings of the New Zealand Grassland Association, 2010, 72: 97-102.

[21] 劉镕源，王紀華，楊貴軍，等．冬小麥葉面積指數地面測量方法的比較[J]．農業(yè)工程學報，2011，27(3)：220－224．Liu Rongyuan, Wang Jihua, Yang Guijun, et al. Comparison of ground-based LAI measuring methods on winter wheat[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(3): 220－224. (in Chinese with English abstract)

[22] 胡小平，王長發(fā). SAS基礎及統(tǒng)計實例教程[M]. 西安：西安地圖出版社，2001.

[23] 張喜英，裴冬，胡春勝. 太行山山前平原冬小麥和夏玉米灌溉指標研究[J]. 農業(yè)工程學報，2002，18(6)：36－41. Zhang Xiying, Pei Dong, Hu Chunsheng. Index system for irrigation scheduling of winter wheat and maize in the piedmont of mountain Taihang[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2002, 18(6): 36－41. (in Chinese with English abstract)

[24] LaRue J L. Variable rate irrigation 2010 field results[EB/OL]. (2011-07-26)[2016-01-01]. http://elibrary.asabe.org/abstract. asp?aid=37311.

[25] Hedley C B, Yule I J. Soil water status mapping and two variable-rate irrigation scenarios[J]. Precision Agriculture, 2009, 10: 342－355.

[26] Stone K C, Bauer P J, Busscher W J, et al. Variable-rate

irrigation management using an expert system in the eastern coastal plain[J]. Irrigation Science, 2015, 33: 167－175.

[27] Ritchie J T, Amato M. Field evaluation of plant extractable soil water for irrigation scheduling[J]. Acta Horticulturae, 1990, 278: 595－615.

[28] DeJonge K C, Kaleita A L, Thorp K R. Simulating the effects of spatially variable irrigation on corn yields, costs, and revenue in Iowa[J]. Agricultural Water Management, 2007, 92(1): 99－109.

[29] 李久生，饒敏杰，張建君. 土壤及噴灌水量不均勻性對干旱區(qū)春小麥產量影響的試驗研究[J]. 農業(yè)工程學報，2002，18(3)：15－21. Li Jiusheng, Rao Minjie, Zhang Jianjun. Influence of spatial variation of soil and nonuniform sprinkler irrigation on spring wheat yield in arid regions[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2002, 18(3): 15－21. (in Chinese with English abstract)

[30] Zhao Weixia, Li Jiusheng, Li Yanfeng.Comparison of irrigation scheduling with water balance and sensor-based methods for a variable rate irrigation system of center pivot [EB/OL].(2016-07-17)[2016-01-01].http://elibrary.asabe.org

Yields and water-saving effects of crops as affected by variable rate irrigation management based on soil water spatial variation

Zhao Weixia, Li Jiusheng※, Yang Rumiao, Li Yanfeng

(100048,)

Improving the uniformity of crop growth parameters and yield across the field is one of the objectives in utilizing the variable rate irrigation (VRI) technology. To improve the water management level of VRI, the crop growth parameters, including plant height, leaf area index (LAI), aboveground dry matter, leave relative chlorophyll content (SPAD), yield, and their spatial uniformity were compared between VRI and uniform rate irrigation (URI) managements, and the influences of soil available water holding capacity (AWC) on crop growth parameters and yield were also analyzed. This study was conducted at the experimental station of China Agricultural University in Zhuozhou, Hebei Province (39.45°N and 115.85°E) in 2014. This site belongs to the Taihang mountain alluvial flood plain and experiences a warm and semi-humid monsoon climate with an annual mean temperature of 11.6°C and an annual mean precipitation of 563.3 mm.One quadrant of the center-pivot controlled area (1.64 hm2) was used in the experiment. The main soil types were loam and sandy loam, and both the coefficients of variation for sand and clay percentiles increased with depth in the site. During the growing seasons of winter wheat and summer maize, the seasonal rainfall for winter wheat (61.2 mm) was substantially lower than the long-term average, while it (311.6 mm) was ample for summer maize. According to the relationships between field capacity, wilting point, and clay, silt, and sand percentiles measured from the upper 0.6 m of the profile, the experimental area was delineated into 4 management zones with AWC varying from 152.2 mm to 161.4 mm for zone 1, from 161.4 mm to 170.9 mm for zone 2, from 170.9 mm to 185.1 mm for zone 3, and from 185.1 to 204.7 mm for zone 4. Each zone was then equally divided into 2 sub-zones to represent the VRI and URI treatments. Each VRI treatment was individually managed with an equal irrigation trigger point of 0.45AWC. For the URI treatments, irrigation was triggered when soil water content in the management zone having minimum AWC values (zone 1) depleted to 0.45AWC. The soil water content was measured with Trime access tubes and Decagon soil moisture sensors. A center-pivot irrigation machine that had been modified to enable to apply variable rate in different plots was used to implement the experiments. The results indicated that the soil water content was different among management zones with different AWC values, and the difference increased during the growing season of winter wheat as the less rainfall than that for summer maize. Correspondently, the total irrigation amount was same for winter wheat between URI and VRI treatments, while it was 14.1% lower for summer maize than that for the URI treatment since the frequent precipitation decreased the dependence on irrigation. Compared with URI management, there was minor effect of VRI management on spatial uniformity of winter wheat and summer maize; the variation coefficient values for all parameters of winter wheat increased from 0.02 to 0.04, while it decreased by 0.05 and 0.02 for aboveground dry matter and yield of summer maize, respectively. No significant influences of VRI management on crop growth parameters and yield were detected, while their differences were significant among the management zones. The experimental results suggested that the management zones with different AWC should be managed individually with different irrigation trigger points to improve the uniformity of crop growth and to maximize the crop yield.

soils; irrigation; crops; soil available water holding capacity; management zone; winter wheat; summer maize; uniformity

10.11975/j.issn.1002-6819.2017.02.001

S274.3; S275.5

A

1002-6819(2017)-02-0001-07

2016-07-07

2016-12-19

國家自然科學基金項目（51309251）；中國水科院科研專項（2016TS05）

趙偉霞，女，河南長葛人，高級工程師，博士，主要從事節(jié)水灌溉理論與技術研究。北京中國水利水電科學研究院流域水循環(huán)模擬與調控國家重點實驗室，100048。Email：zhaowx@iwhr.com

李久生，男，河北邢臺人，研究員，博士生導師，主要從事節(jié)水灌溉理論與技術研究。北京中國水利水電科學研究院流域水循環(huán)模擬與調控國家重點實驗室，100048。Email：lijs@iwhr.com

趙偉霞，李久生，楊汝苗，栗巖峰.基于土壤水分空間變異的變量灌溉作物產量及節(jié)水效果[J]. 農業(yè)工程學報，2017，33(2)：1－7. doi：10.11975/j.issn.1002-6819.2017.02.001 http://www.tcsae.org

Zhao Weixia, Li Jiusheng, Yang Rumiao, Li Yanfeng. Yields and water-saving effects of crops as affected by variable rate irrigation management based on soil water spatial variation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(2): 1－7. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.02.001 http://www.tcsae.org