陳超飛，柳雙環(huán)，郭大辛，徐芳平，馬孝義

(1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，陜西 楊凌 712100;2.西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室，陜西 楊凌 712100)

玉米是中國第一大糧食作物。2015年中國玉米種植面積達(dá)3 811.9萬hm2，總產(chǎn)量達(dá)2.246億t[1]。關(guān)中地區(qū)是陜西糧食主產(chǎn)區(qū)，主要實(shí)行冬小麥-夏玉米輪作制[2]。該區(qū)夏季容易出現(xiàn)極端天氣，降雨分配不均導(dǎo)致伏旱時有發(fā)生[3]，這對夏玉米的生產(chǎn)極為不利，氣候變化、水資源短缺和田間管理不當(dāng)更加劇了生產(chǎn)的不穩(wěn)定性[4]。因此，研究夏玉米節(jié)水灌溉對保障關(guān)中地區(qū)糧食穩(wěn)產(chǎn)高產(chǎn)、提高水資源利用率有重要意義。

田間試驗(yàn)往往耗時耗力，且易受氣候變化和病蟲害等環(huán)境因素影響，結(jié)合作物模型進(jìn)行農(nóng)田或區(qū)域尺度的作物生產(chǎn)模擬、預(yù)測和評估成為國際上的研究熱點(diǎn)[5]。作物模型克服了傳統(tǒng)試驗(yàn)研究對象單一、試驗(yàn)周期長、時間和空間有限制的缺點(diǎn)，有助于田間管理決策[6]。相較其他作物模型而言，AquaCrop模型以良好的適用性、較少的參數(shù)和較好的模擬效果吸引了大批科研工作者的關(guān)注和應(yīng)用[7]。世界各地學(xué)者對該模型適用性展開的大量研究表明，利用該模型能較好模擬不同水分處理多種農(nóng)作物的產(chǎn)量形成機(jī)制[8-10]。國內(nèi)學(xué)者應(yīng)用AquaCrop模型對冬小麥[11-12]，玉米[13-14]，棉花[15]，大豆[16]，大蔥[17]和西北胡麻[18]等多種作物進(jìn)行模型校準(zhǔn)和驗(yàn)證，評價了生物量、產(chǎn)量和水分利用效率等方面的精度和適用性。目前，AquaCrop模型主要用于模擬不同作物產(chǎn)量和水分利用效率方面的研究。國內(nèi)較缺乏對夏玉米不同水氮組合情境下的模擬研究，且AquaCrop作物模型在關(guān)中地區(qū)針對夏玉米模擬的研究更是少見報道。

于2016年和2017年6-10月在陜西武功縣高新農(nóng)業(yè)示范園針對夏玉米設(shè)置不同梯度灌水和施氮處理的大田試驗(yàn)。利用試驗(yàn)數(shù)據(jù)對AquaCrop模型進(jìn)行校驗(yàn)，評價其適用性，最后模擬不同降水年型下不同灌溉策略對夏玉米產(chǎn)量和水分利用效率的影響，以期為模型在關(guān)中地區(qū)的應(yīng)用和夏玉米不同典型年的灌溉制度提供參考。

1 材料與方法

1.1 研究區(qū)概況

田間試驗(yàn)分別于2016年和2017年6-10月在陜西武功縣高新農(nóng)業(yè)示范園(34°21′N，108°03′E)進(jìn)行，冬小麥-夏玉米輪作是該區(qū)主要的糧食種植制度。試驗(yàn)區(qū)位于關(guān)中平原中西部，地處渭河一、二階地。該區(qū)多年平均降水量596.8 mm，夏玉米生長季6-9月的多年平均降水量為360.2 mm，有效降水量316.5 mm；平均氣溫為13.3℃，平均日照時數(shù)2 095 h，年輻射總量480.90 kJ，無霜期221 d；土壤為典型土，0～20 cm土層容重為1.32 g·cm-3，pH值為8.14，有機(jī)質(zhì)含量10.20 g·kg-1，全氮含量0.96 g·kg-1，堿解氮質(zhì)量分?jǐn)?shù)為60.3 mg·kg-1，速效磷質(zhì)量分?jǐn)?shù)為10.28 mg·kg-1，速效鉀質(zhì)量分?jǐn)?shù)為128.5 mg·kg-1。

1.2 試驗(yàn)方案設(shè)計

2016-2017年田間試驗(yàn)在同一地塊進(jìn)行，設(shè)置水、氮2個因素。水分設(shè)4個水平：I0，雨養(yǎng)(不灌水)；I1，低水30 mm；I2，適水60 mm；I3高水90 mm。氮素4個水平：N0，不施氮；N1，為推薦施氮量的67%，120 kg·hm-2；N2，推薦施氮量，180 kg·hm-2；N3，推薦施氮量的133%，240 kg·hm-2，共16個處理，3次重復(fù)。試驗(yàn)小區(qū)為東西走向，寬5.7 m，長28 m，面積159.6 m2。裂區(qū)設(shè)計，水分為主區(qū)，肥料為副區(qū)，小區(qū)間由壟隔開，壟寬0.5 m，高0.2 m。灌水方式為畦灌，灌水時間依據(jù)旱情決定，2016年灌水日期為8月13日，2017年灌水日期為7月11日。供試氮肥為尿素(含N≥46%)，磷肥為過磷酸鈣(含P2O5≥16%)，鉀肥為硫酸鉀(含K2O≥50%)。各處理施磷肥(60 kg·hm-2)和鉀肥(60 kg·hm-2)，均作為基肥在播前一次性施入。玉米供試品種為鄭單958，行距65 cm，株距30 cm，密度為51 300株·hm-2。2016年于6月14日播種，10月1日成熟；2017年于6月15日播種，10月7日成熟。其他管理措施如雜草控制、病蟲害管理等按照一般大田管理規(guī)范進(jìn)行。

1.3 觀測指標(biāo)及處理

1.3.1 葉面積指數(shù) 每隔10～15 d測量標(biāo)記選擇供試植株的所有綠葉葉片的長度(從葉領(lǐng)到葉尖)及寬度(葉片最寬處)，葉面積指數(shù)(leaf area index，LAI)采用下式計算[19]：

(1)

式中，ρ為種植密度(株·hm-2)；n為第j株玉米的總?cè)~片數(shù)；m為測定株數(shù)；Lij為葉片葉領(lǐng)到葉尖的長度(m)；Bij為葉片最寬處寬度(m)；0.75為與葉形有關(guān)的葉面積回歸系數(shù)。

1.3.2 冠層覆蓋度 冠層覆蓋度(canopy cover，CC)是指土壤表面綠色冠層覆蓋的面積百分比，由LAI計算獲得，其計算公式如下[20]：

CC=1.005[1-exp(-0.6×LAI)]1.2

(2)

1.3.3 土壤含水量 本試驗(yàn)使用德國IMKO公司制造的TRIME-TDR測定土壤含水量(soil water content，SWC)，測量深度為1m，間隔為20 cm。每次播種前和收獲后，用土鉆取土測量土壤含水量，并校準(zhǔn)TRIME-TDR儀器。

1.3.4 地上部生物量及產(chǎn)量 玉米進(jìn)入拔節(jié)期后，每隔10～15 d測定地上部生物量(aboveground biomass，B)，每小區(qū)取3株能夠代表作物平均長勢的植株，剪去地下部分后于105℃下殺青0.5 h,75℃烘至恒重，測其干重，乘以種植密度即為地上部生物量。玉米成熟后，選取能夠代表作物整體長勢的20株成熟玉米，果穗風(fēng)干后經(jīng)人工脫粒，75℃烘至恒重后測算產(chǎn)量。

1.3.5 水分利用效率 水分利用效率(water use efficiency，WUE)是指單位水分消耗獲得的生物量或產(chǎn)量，產(chǎn)量水分利用效率WUE可按下式計算：

(3)

式中，Y為籽粒產(chǎn)量(kg·hm-2)；ET為作物耗水量，(mm)，由AquaCrop模型模擬計算獲得。

1.4 構(gòu)建AquaCrop模型數(shù)據(jù)庫

AquaCrop模型需要輸入的數(shù)據(jù)包括氣象、作物、土壤、田間管理等。

1.4.1 氣象數(shù)據(jù) 武功氣象站1961-2017年的逐日降水、氣溫等氣象數(shù)據(jù)來源于國家氣象科學(xué)數(shù)據(jù)共享服務(wù)平臺(http://data.cma.cn/site/index.html)中國地面氣候資料日值數(shù)據(jù)集(V3.0)。參考作物蒸發(fā)蒸騰量ET0采用Penman-Montieth方程[21]進(jìn)行計算。2016年和2017年夏玉米生育期內(nèi)逐日最高和最低氣溫、太陽輻射、參考蒸發(fā)蒸騰量、最小相對濕度及降水量見圖1。

1.4.2 土壤數(shù)據(jù) 模型輸入的土壤數(shù)據(jù)包括土層深度及層數(shù)、土壤質(zhì)地、容重、各土層的凋萎含水量、田間持水量和飽和含水量等參數(shù)。試驗(yàn)區(qū)土壤為典型土，用土鉆分別在0～20、20～40、40～60、60～80 cm和80～100 cm 5個深度取土，風(fēng)干過2 mm篩后采用英產(chǎn)Mastersizer 2000型激光粒度儀測定土壤質(zhì)地。各土層容重采用環(huán)刀法測定，采用離心法測定土壤水分特征曲線，用RETC軟件擬合得到凋萎含水量、田間持水量和飽和含水量，輸入AquaCrop模型建立土壤數(shù)據(jù)庫文件，詳細(xì)土壤參數(shù)如表1所示。

圖1 2016-2017年武功縣夏玉米生育期氣象數(shù)據(jù)Fig.1 Meteorological data of summer maize growth period in Wugong County from 2016 to 2017

1.4.3 作物參數(shù)數(shù)據(jù) 作物參數(shù)中的植株密度和物候期等參數(shù)根據(jù)試驗(yàn)記錄直接輸入模型。歸一化水分生產(chǎn)力、參考收獲指數(shù)，冠層增長衰減系數(shù)和水分脅迫響應(yīng)系數(shù)等重要參數(shù)根據(jù)模型手冊和文獻(xiàn)提供的取值范圍采用“試錯法”進(jìn)行修正，部分默認(rèn)參數(shù)如基底溫度、上限溫度等采用模型手冊上的推薦值。參照Vanuytrecht等[22]介紹的校準(zhǔn)順序：第一步，校準(zhǔn)冠層覆蓋度；第二步，校準(zhǔn)地上部生物量；第三步，校準(zhǔn)產(chǎn)量。利用2017年田間試驗(yàn)數(shù)據(jù)進(jìn)行模型調(diào)試，使模型模擬結(jié)果與實(shí)測結(jié)果相符合，經(jīng)過校準(zhǔn)，得到‘鄭單958’的品種參數(shù)，部分參數(shù)詳見表2。

AquaCrop模型沒有明確考慮養(yǎng)分的循環(huán)與平衡，而是采用半定量的評估方法來描述養(yǎng)分脅迫對生物量的影響[23]。模型引入相對于最大地上部生物量的比例系數(shù)Brel，公式如下：

(4)

式中，Bstress為無肥料脅迫和水分脅迫條件下得到的地上部生物量(t·hm-2)；Bref為有肥料脅迫但無水分脅迫條件下得到的地上部生物量(t·hm-2)；Brel的范圍為0%～100%，0%表示有養(yǎng)分脅迫時作物生長失敗，100%表示沒有養(yǎng)分脅迫，Bstress和Bref參照當(dāng)?shù)厣a(chǎn)現(xiàn)狀、前人研究成果[24-26]以及本次試驗(yàn)數(shù)據(jù)確定。模型內(nèi)嵌自動校準(zhǔn)程序，使用4個脅迫參數(shù)來表現(xiàn)肥料脅迫的響應(yīng)：冠層擴(kuò)張(Ksexp,f)，最大冠層覆蓋度(KsCCx)，生物量水分生產(chǎn)力(KsWP)和冠層衰減(fCDecline)[27]。通過校準(zhǔn)，結(jié)果如表3。

表1 試驗(yàn)區(qū)土壤數(shù)據(jù)

表2 AquaCrop模型玉米校準(zhǔn)參數(shù)

1.5 模型評價指標(biāo)

在本研究中，采用以下統(tǒng)計指標(biāo)來評價模型校準(zhǔn)和驗(yàn)證的精度：決定系數(shù)R2，均方根誤差RMSE，標(biāo)準(zhǔn)均方根誤差NRMSE，一致性指數(shù)d和納什效率系數(shù)ENS。

(5)

(6)

(7)

(8)

(9)

表3 土壤肥料脅迫參數(shù)校準(zhǔn)表

1.6 不同降水年型灌溉情景模擬

利用Pearson-Ⅲ型分布模型對1961-2017年夏玉米生育期降水量進(jìn)行分析，得到濕潤年1990年、平水年1989年和干旱年1999年降水量分別為377.5、283.1 mm和207.3 mm。選擇苗期、拔節(jié)期、抽雄期和灌漿期4個關(guān)鍵生育期制定模擬灌溉方案,土壤初始含水量設(shè)置為田間持水量；肥料脅迫設(shè)置為100%，即無肥料脅迫；灌溉模擬情景設(shè)置為雨養(yǎng)(P0)、灌1水(P1～P4)和灌2水(P5～P10)，灌水定額分別為60 mm和120 mm，按全排列共11種處理，各灌溉方案見表4。

2 結(jié)果與分析

2.1 模型校準(zhǔn)結(jié)果

利用2017年的試驗(yàn)數(shù)據(jù)進(jìn)行校準(zhǔn)，試驗(yàn)共有4個水分處理，限于篇幅，主要分析充分灌水處理I3與雨養(yǎng)處理I0。不同施氮水平下夏玉米冠層覆蓋度校準(zhǔn)結(jié)果如圖2，地上部生物量校準(zhǔn)結(jié)果如圖3，土壤含水量校準(zhǔn)結(jié)果如圖4。

2017年各處理的冠層覆蓋度和地上部生物量的校準(zhǔn)誤差統(tǒng)計指標(biāo)見表5。所有處理下冠層覆蓋度的R2、RMSE和ENS分別為0.868、11.83%、0.743，表明AquaCrop模型能夠較為準(zhǔn)確地模擬出不同水氮組合情景下冠層覆蓋度的動態(tài)變化過程。所有處理地上部生物量的R2、RMSE和ENS分別為0.975、0.970 t·hm-2、0.968，表明AquaCrop模型能夠準(zhǔn)確地模擬出不同水氮組合情景地上生物量的動態(tài)變化過程。所有處理下土壤含水量的R2、RMSE和ENS分別為0.743、20.1 mm、0.490,表明AquaCrop模型能夠較為準(zhǔn)確地模擬出不同水氮組合情景下土壤含水量的動態(tài)變化過程。

表4 模擬灌溉方案/mm

注：N0，N1，N2，N3：0%，67%，100%和133%的推薦施氮量；I0，I1，I2，I3：雨養(yǎng)，虧水，適水，充分灌水。下同。Note: N0, N1, N2, N3: Application of 0%, 67%, 100% and 133% of the recommended nitrogen rate; I0, I1, I2, I3:Rainfed, deficit, suitable and adequate irrigation. The same below.圖2 充分灌水與雨養(yǎng)條件下夏玉米冠層覆蓋度校準(zhǔn)結(jié)果Fig.2 Calibration results of summer maize canopy cover under full irrigation and rainfed conditions.

圖3 充分灌水與雨養(yǎng)條件下夏玉米地上部生物量動態(tài)模擬Fig.3 Calibration results of summer maize aboveground biomass under full irrigation and rainfed conditions

圖4 充分灌水與雨養(yǎng)條件下夏玉米土壤含水量校準(zhǔn)結(jié)果Fig.4 Calibration results of summer maize soil water content under full irrigation and rainfed conditions

如圖5(a)所示，所有試驗(yàn)處理下夏玉米產(chǎn)量的R2、RMSE、NRMSE、d和ENS分別為0.879、0.297 t·hm-2、4.390、0.959和0.812。如圖5(b)所示，最終生物量的R2、RMSE、NRMSE、d和ENS分別為0.873、1.108 t·hm-2、6.977、0.904和0.502。各處理的相對誤差見表6，產(chǎn)量的相對誤差最小為0.60%，對應(yīng)處理為I3N2，最大為7.63%，對應(yīng)處理為I1N3；最終生物量的相對誤差最小為1.27%，對應(yīng)處理為I1N0，最大為11.04%，對應(yīng)處理為I0N2。

2.2 模型驗(yàn)證結(jié)果

采用2016年的試驗(yàn)數(shù)據(jù)進(jìn)行模型驗(yàn)證，土壤含水量的驗(yàn)證結(jié)果如圖6所示，誤差統(tǒng)計見表7，所有試驗(yàn)處理下土壤含水量的R2、RMSE、ENS分別為0.691、15.7 mm和0.671。如圖7(a)所示，所有試驗(yàn)處理下夏玉米產(chǎn)量的R2、RMSE、NRMSE、d和ENS分別為0.919、0.249 t·hm-2、4.112、0.977和0.915。如圖7(b)，最終生物量的R2、RMSE、NRMSE、d和ENS分別為0.860、0.977 t·hm-2、6.407、0.933和0.694。各處理的相對誤差見表8，產(chǎn)量的相對誤差最小為1.51%，對應(yīng)處理為I3N0，最大為6.65%，對應(yīng)處理為I3N1；最終生物量的相對誤差最小為1.24%，對應(yīng)處理為I3N2，最大為12.05%，對應(yīng)處理為I0N3。

表5 2017年不同施氮及灌溉水平下冠層覆蓋度、地上部生物量和土壤含水量校準(zhǔn)誤差統(tǒng)計

圖5 2017年所有灌溉施氮處理下產(chǎn)量(a)和最終生物量(b)校準(zhǔn)結(jié)果Fig.5 Calibration results of yield (a) and final biomass (b) under all differentnitrogen and irrigation levels in 2017

表6 2017年所有灌溉施氮處理下產(chǎn)量、最終生物量校準(zhǔn)結(jié)果

圖6 充分灌水與雨養(yǎng)條件下夏玉米土壤含水量驗(yàn)證結(jié)果Fig.6 Validation results of summer maize soil water content under adequate irrigation and rainfed conditions

表7 2016年不同施氮及灌溉水平下土壤含水量驗(yàn)證誤差統(tǒng)計

2.3 不同降水年型灌溉模擬結(jié)果

利用校準(zhǔn)后的模型對3種不同降水年型進(jìn)行不同灌溉方案的模擬，夏玉米產(chǎn)量與水分利用效率的模擬結(jié)果如圖8、表9所示。不同灌溉方案下，濕潤年產(chǎn)量在7.899～7.998 t·hm-2，平水年產(chǎn)量在7.724～8.066 t·hm-2，干旱年產(chǎn)量在1.984～7.349 t·hm-2。濕潤年由于降水量豐沛，雨養(yǎng)和灌溉條件下的水分利用效率基本不變，其模擬值為2.19 kg·m-3或2.21 kg·m-3；平水年水分利用效率模擬值為2.50～2.72 kg·m-3；干旱年水分利用效率模擬值為0.86～2.33 kg·m-3。

3 討 論

3.1 AquaCrop模型的適用性

本研究利用AquaCrop模型對關(guān)中地區(qū)夏玉米栽培進(jìn)行校準(zhǔn)和驗(yàn)證，夏玉米產(chǎn)量和地上部生物量的模擬值與實(shí)測值一致性高，模擬誤差在可接受范圍內(nèi)，且能準(zhǔn)確模擬土壤含水量的變化趨勢，與國內(nèi)外的研究成果基本吻合。如倪玲等[6]評價了AquaCrop模型在黃土塬區(qū)的適用性，結(jié)果顯示夏玉米模擬產(chǎn)量與實(shí)測產(chǎn)量間的決定系數(shù)R2為0.927，相對誤差在-2.479%～11.182%之間；模擬地上部生物量與實(shí)測地上部生物量間的R2為0.784，模型對產(chǎn)量的模擬效果優(yōu)于對生物量的模擬，也與本研究結(jié)果相似。Abedinpour等[27]模擬了不同水氮組合情景下的玉米產(chǎn)量、生物量及水分生產(chǎn)力的變化特征，其中產(chǎn)量R2為0.96，RMSE為0.1 t·hm-2，ENS為0.98；生物量的R2為0.9，RMSE為0.75 t·hm-2，ENS為0.95。Ran等[14]模擬不同虧缺條件下制種玉米冠層覆蓋度的模擬值與實(shí)測值的R2為0.818，RMSE為12.9%，ENS為0.811；生物量R2為0.929，RMSE為1.972 t·hm-2，ENS為0.903；土壤含水量R2為0.736，RMSE為33.1%，ENS為0.364。這些都表明AquaCrop模型可以準(zhǔn)確模擬出不同田間管理?xiàng)l件下玉米產(chǎn)量和地上部生物量，說明AquaCrop模型可以為關(guān)中地區(qū)夏玉米田間灌溉和決策提供指導(dǎo)。

圖7 2016年不同灌溉施氮處理下產(chǎn)量(a)和最終生物量(b)驗(yàn)證結(jié)果Fig.7 Validation results of yield (a) and final biomass (b) under all different nitrogen and irrigation levels in 2016

表8 2016年不同灌溉施氮產(chǎn)量、最終生物量驗(yàn)證結(jié)果

圖8 不同降水年型不同灌溉方案的產(chǎn)量(a)和水分利用效率(b)模擬結(jié)果Fig.8 Simulated yield (a) and WUE (b) in different irrigation schemes under wet, normal and dry year

表9 不同降水年型不同灌溉方案的產(chǎn)量與水分利用效率模擬結(jié)果

與此同時，本文也反映出模型存在一定不足，模型輕微高估了雨養(yǎng)I0條件下的冠層覆蓋度和地上部生物量(如圖2(e)～2(h)與圖3(e)～3(h))，即水分脅迫越嚴(yán)重，模型的模擬誤差越大，這與李玥等[18]研究胡麻的結(jié)果類似；模型整體高估了土壤含水量，尤其是2017年的所有處理(如圖4(a)～4(h))，這與Ran等[14]的研究結(jié)果相似；同時模型低估了衰落期的冠層覆蓋度，這是由于2017年玉米生育后期降水量過多，日照時數(shù)較少，延緩了玉米的衰老，這可能是模型對于衰老期的水分脅迫響應(yīng)過于敏感[19]，模擬的衰老速度過快造成。因此模型對于不同生育期的水分脅迫響應(yīng)仍需改進(jìn)[14]。

3.2 針對不同降水年型灌溉方案的制定

不同生育期受旱會對夏玉米的生長發(fā)育及產(chǎn)量造成不同程度的影響[28]。生長初期地上部生長緩慢，比較耐旱，一般無需灌水，但若長時間遭遇伏旱，則容易造成弱苗，最終導(dǎo)致大幅減產(chǎn)。玉米進(jìn)入拔節(jié)期后，植株生長旺盛，此時缺水會引起營養(yǎng)體生長不良，明顯制約夏玉米株高和葉面積的增長及干物質(zhì)的積累[29]。抽雄期受旱會導(dǎo)致穗粒數(shù)減少，灌漿期則明顯減少百粒重，均會造成減產(chǎn)[19]。因此，本研究按照夏玉米的生育期實(shí)際劃分情況，分別統(tǒng)計3種降水年型在苗期、拔節(jié)期、抽雄期和灌漿期的有效降水量。濕潤年4個生育期有效降水量分別為206.8、57.6、86.9 mm和26.2 mm，各生育期無明顯水分脅迫；平水年4個生育期有效降水量分別為105.7、14.4、131 mm和32 mm，拔節(jié)期出現(xiàn)水分脅迫；干旱年4個生育期有效降水量分別為129.7、6.5、22.9 mm和48.2 mm，拔節(jié)期和抽雄期出現(xiàn)水分脅迫。

不同降水年型灌溉模擬結(jié)果顯示：雨養(yǎng)條件下，濕潤年產(chǎn)量最高，但水分利用效率低于平水年，這是由于濕潤年的降水在各生育期分布相對較均勻，各生育期未發(fā)生明顯的水分脅迫，所以產(chǎn)量最高；濕潤年的耗水量為363 mm，大于平水年297.2 mm，而產(chǎn)量相差不大，因此導(dǎo)致水分利用效率較低。

灌1水情況下，濕潤年和平水年均在拔節(jié)期灌水60 mm時產(chǎn)量達(dá)到最大，分別為7.998 t·hm-2和8.066 t·hm-2，對應(yīng)的水分利用效率也最大，分別為2.21 kg·m-3和2.72 kg·m-3，這是因?yàn)槠剿暝诎喂?jié)期出現(xiàn)水分脅迫，此時灌水可以提高產(chǎn)量與水分利用效率。對于干旱年，抽雄期灌水60 mm產(chǎn)量最大，為6.618 t·hm-2，對應(yīng)水分利用效率為2.23 kg·m-3；拔節(jié)期灌水60 mm的產(chǎn)量略低，為5.749 t·hm-2，對應(yīng)水分利用效率為2.05 kg·m-3。這是因?yàn)楦珊的暝诎喂?jié)期和抽雄期缺水嚴(yán)重，在這2個生育期灌水均可提高產(chǎn)量，模型模擬結(jié)果顯示抽雄水獲得的產(chǎn)量要高于拔節(jié)水，表明模型對于營養(yǎng)生長階段的水分脅迫響應(yīng)較為敏感。

灌2水情況下，苗期和拔節(jié)期各灌水60 mm，即方案P5，此時濕潤年和平水年的產(chǎn)量最大，分別為7.998、8.170 t·hm-2，且灌水日期越延后產(chǎn)量越低，但如果存在拔節(jié)期灌水的組合，如方案P8和P9，產(chǎn)量均較其他生育期的灌水組合大。干旱年在拔節(jié)期和抽雄期各灌水60 mm時(方案P8)產(chǎn)量最大，為7.349 t·hm-2，且水分利用效率也達(dá)到最大，為2.33 kg·m-3，這些模擬結(jié)果都表明在玉米拔節(jié)期或抽雄期灌水可提高夏玉米產(chǎn)量，與前人研究成果極為吻合[28]，基本說明了不同灌溉時間和灌溉次數(shù)對夏玉米產(chǎn)量的影響。

濕潤年灌水后的產(chǎn)量和水分利用效率增加幅度較小，如方案P2的產(chǎn)量較雨養(yǎng)條件下的產(chǎn)量僅提高1.2%，水分利用效率提高不足1%，因此若能保證出苗整齊，濕潤年不灌水也能夠獲得較高產(chǎn)量和水分利用效率，同時也達(dá)到節(jié)約水資源的目的。平水年雖然方案P5可獲得最高產(chǎn)量8.170 t·hm-2，但相較方案P2的增幅較小，為1.3%，且水分利用效率低于P2，增加灌水量并沒有顯著提升產(chǎn)量，反而會降低水分利用效率，造成水資源浪費(fèi)，因此平水年可在拔節(jié)期灌水60 mm，推薦灌溉方案P2，可節(jié)約50%的灌溉水資源，產(chǎn)量同樣能維持在較高水平。干旱年采用灌溉方案P8能夠獲得最高的產(chǎn)量和水分利用效率，相較于方案P3，產(chǎn)量和水分利用效率增幅分別為11%和4.5%，考慮到干旱年整體降水較少，推薦灌溉方案P8，即在拔節(jié)期和抽雄期各灌水60 mm，可實(shí)現(xiàn)水資源充分利用和作物高產(chǎn)。

4 結(jié) 論

AquaCrop模型可以較好地模擬關(guān)中地區(qū)不同水氮條件下夏玉米冠層覆蓋度、生物量、產(chǎn)量和土壤含水率。通過不同降水年型不同灌溉方案的模擬，依據(jù)模擬結(jié)果及夏玉米生理需水規(guī)律，同時為了實(shí)現(xiàn)高產(chǎn)和高水分利用效率，得出不同年型的優(yōu)化灌溉制度為：若能保證出苗整齊，濕潤年可不灌水，同樣能夠獲得較高產(chǎn)量和水分利用效率，同時也達(dá)到節(jié)約水資源的目的；平水年推薦在拔節(jié)期灌水60 mm，能節(jié)約50%水資源，亦能保證穩(wěn)產(chǎn)；干旱年推薦在拔節(jié)期和抽雄期各灌水60 mm，可獲得高產(chǎn)和高水分利用效率。