APSIM模型在河北山前平原農業耗水研究中的應用

趙彥茜， 安 塞， 齊永青， 朱 驥

(1.河北地質大學，河北石家莊 050031； 2.河北科技大學，河北石家莊 050018；3.中國科學院遺傳與發育生物學研究所農業資源研究中心/中國科學院農業水資源重點實驗室，河北石家莊 050021)

河北山前平原主要是指河北省境內京津以南的河北山前平原，屬暖溫帶半干旱半濕潤大陸性季風氣候，雨熱同期，降水主要集中于夏季，適合多種作物生長。

冬小麥—夏玉米一年兩熟制是河北山前平原主要的種植制度[1]，但該地區大部分區域的降水量無法滿足一年兩熟制下作物對水分的要求，原因在于冬小麥生育期內水分虧缺非常嚴重[2]。通常情況下，該地區多年平均降水量為400～550 mm，冬小麥生育期內降水量為90～130 mm，夏玉米生育期內降水量為300～360 mm，而冬小麥年平均耗水量為380～450 mm，夏玉米年平均耗水量則為270～330 mm[3-5]，冬小麥生育期水分虧缺為250～320 mm，須要進行多次灌溉來保證產量，而夏玉米只要適時播種，自然降水可滿足其對水分的要求，生育期內基本不須要灌溉或者僅須要灌溉1次出苗水，只有在極少的尤其干旱的年份須要進行二次灌溉，以緩解“卡脖旱”的影響。河北山前平原作物高產主要是靠提取地下水灌溉來保證，而長期的灌溉已經造成河北山前平原地下水資源的嚴重虧損，導致河北山前平原地下水位持續下降[6]。所以，農業生產耗水研究對于供水問題突出的河北山前平原來說十分重要，而應用于這一研究的主要就是進行田間定位試驗和利用作物模型模擬這2種方法，本研究主要是在已有田間試驗資料的基礎上應用APSIM模型對河北山前平原的農業耗水狀況進行研究分析。

作物模型可以在已有研究資料基礎上對在多種模擬條件下的作物生長情況進行模擬，預測未知風險，優化農業資源管理措施，為農業生產的管理經營提供科學指導，APSIM模型就是作物模型中的典型代表。APSIM模型是一種可用于模擬農業生產系統中各主要組成部分的機制模型，它是由澳大利亞的聯邦科工組織以及昆士蘭州政府的農業生產系統組(agricultural production systems research unit，簡稱APSRU)所開發建立的優秀作物模型[7-8]。APSIM模型建立最初的目的是在農業系統里進行長期資源管理試驗時，對在氣候變化、作物的遺傳特征、土壤環境以及管理措施等因子影響下的作物生產力提供一個更準確地預測[9]。在國外，APSIM模型關于水資源-產量綜合效應和水分利用的研究較多，涉及面較廣[10-13]。Acua等利用APSIM模型對澳大利亞塔斯馬尼亞島上農業用地的地表水平衡、水分利用效率和小麥產量進行模擬，分析了灌溉、施肥等因素對產量和水分利用效率的影響[10]。而在國內，關于APSIM模型的應用研究也漸趨增多[14-24]。Sun等利用APSIM模型在華北平原量化灌溉對于地下水和產量的影響[17]。van Oort等利用APSIM模型模擬了華北平原種植制度變化下地下水的變化情況[18]。

為研究河北山前平原的農業耗水狀況，以中國科學院欒城農業生態系統試驗站2009—2013年的冬小麥和夏玉米實測資料對APSIM模型進行率定和驗證，對河北山前平原水資源-產量綜合效應進行研究分析。

1 材料與方法

1.1 試驗區概況

中國科學院欒城農業生態系統試驗站位于河北省石家莊市欒城區，地處37°53′N、114°41′E，海拔為50.1 m。該試驗站位于河北平原北部的太行山山前平原，屬暖溫帶半濕潤季風氣候，年降水量為400～550 mm，其中7—9月為雨季，降水量為300～360 mm，占年降水量60%以上，年均溫為10～13 ℃，無霜期為180～230 d，主要作物種植模式為冬小麥—夏玉米一年兩熟。

1.2 試驗方法

本次試驗研究模擬驗證時間段為2009—2013年，熟制是冬小麥—夏玉米一年兩熟制，冬小麥品種為科農199，夏玉米品種為鄭單958。冬小麥播種日期為10月5—15日，夏玉米播種日期為6月9—17日，冬小麥生育期內進行3～4次灌溉，夏玉米生育期內進行0～2次灌溉。

1.3 模型參數

1.3.1 氣象參數 本研究所需的氣象參數包括最高氣溫、最低氣溫、相對濕度、降水量、日照時數等，均來源于中國科學院欒城農業生態系統試驗站田間Milos520自動氣象站。

1.3.2 土壤參數 APSIM模型土壤參數主要包括各土層容重(BD)、飽和體積含水量(SAT)、萎蔫系數(LL15)、風干系數(AirDry)等參數(表1)。

表1 土壤水分參數

1.3.3 作物參數 作物參數主要是根據田間觀測數據來確定，研究中的APSIM模型作物屬性模塊主要包括該研究區域冬小麥和夏玉米品種遺傳特性參數、作物生長發育進程、植株形態以及產量形成等相關參數(表2)。

表2 冬小麥和夏玉米作物參數

1.4 降水年型分類

為研究河北山前平原降水對農業生產的影響，根據欒城站試驗區的降水資料對該地區的降水年型進行分類，采用的標準為國內常用的降水保證率分類標準[25-28]。降水保證率表示降水量某一界限值出現的可靠程度，根據降水保證率將降水年份劃分為不同的年型：25%保證率的降水年份為濕潤年，50%保證率的降水年份作為平水年，75%保證率的降水年份為干旱年。

1.5 模型檢驗方法

模型檢驗方法主要是以下2個國際上統一的指標：確定系數R2和均方根誤差(root mean square error，簡稱RMSE)。計算公式如下：

(1)

(2)

2 結果與分析

2.1 APSIM模型驗證

利用APSIM模型對試驗區0～1.8 m土壤含水量SWC、產量、蒸散量ET等進行模擬，模擬結果表明，SWC模擬結果的決定系數為0.79(圖1)，冬小麥和夏玉米產量模擬結果的均方根誤差和確定系數分別為502 kg/hm2和0.83(圖2)，作物蒸散量ET模擬結果的均方根誤差和決定系數分別為32.7 mm 和0.88(圖2)，每月蒸散量ET模擬結果的均方根誤差和確定系數分別為13.2 mm和0.925(圖3)。這些結果表明，APSIM模型適合應用于研究分析河北山前平原水資源消耗和糧食生產的綜合效應。

APSIM模型對0.0～1.8 m的土壤含水量的模擬結果較好。根據實測資料，可以觀察到土壤含水量在1年中的7、8月達到全年最大值，超過600 mm，之后隨著小麥生育期的到來，土壤含水量維持在450～550 mm這個略低的水平，APSIM模型的模擬結果可以將土壤含水量的這一動態變化趨勢很清楚地展現出來，對研究農業生產和土壤水動態變化的關系有重要意義。

APSIM模型在一定程度上可以模擬出氣候風險對作物生長的影響以及相應的產量波動，相關性分析擬合精度達0.83。冬小麥和夏玉米生育期內蒸散量ET的年際變化較為明顯，在APSIM模型的模擬結果中，這一動態變化模擬得也較為準確，擬合精度達到0.91。

通過2009—2013年的每月蒸散量ET實測數據可以發現，冬小麥生育期內月蒸散量ET的峰值是在5月，夏玉米生育期內月ET的峰值是在8月，這分別對應的是冬小麥和夏玉米的灌漿期，而月ET全年谷值則是在6月和冬季月，全年呈現出“兩峰一谷”的特點。APSIM模型的模擬結果非常明顯地把這一特點表現出來，確定系數達到0.89，模擬效果較好。

2.2 河北山前平原1986—2015年不同降水年型情況下冬小麥—夏玉米一年兩熟制的作物生產耗水分析

根據研究區域1986—2015年的降水資料，可以計算出年降水豐水年的分類界限為570 mm，枯水年分類界限為 400 mm，小麥季豐水年的分類界限為170 mm，枯水年分類界限為105 mm，玉米季豐水年的分類界限為400 mm，枯水年分類界限為250 mm。根據上述分類標準，可將1986—2015年這30年分成豐水年、平水年、枯水年等3個類型(圖4)。

根據調整驗證后的模型參數，本研究運用APSIM模型對河北山前平原1986—2015年的農業生產耗水進行模擬(圖5)。小麥生長階段為每年10月15日至次年6月14日，生長階段內灌溉4次，每次灌溉量為80 mm；玉米生長階段為每年6月15日至10月5日，生長階段內灌溉2次，每次灌溉量為60 mm。通過模擬研究降水對農業生產的影響，并從全年、小麥生育期、玉米生育期等3個角度分析河北山前平原近30年的農業生產耗水規律。

根據模擬結果，可計算得出豐水年平均年產量為 13 622 kg/hm2，平水年平均年產量為13 698.8 kg/hm2，枯水年平均年產量為12 814 kg/hm2，30年總平均年產量為 13 444.9 kg/hm2。根據上述驗證所用實測數據，這個模擬結果是較為符合實際的。模擬結果表明，平水年和豐水年年產量水平均比較高，枯水年產量水平最低，且豐水年年產量水平略低于平水年，這主要是因為豐水年部分年份降水分配較差或者夏季暴雨集中對作物生產產生了不利影響， 例如根據驗證所用實測數據，2009年年降水量為609.5 mm，但年產量僅為 11 163.9 kg/hm2，這主要是因為夏季“9.20暴雨”使得玉米產量大幅減產。而且雖然降水量大小對農業生產有一定程度的影響，但由于灌溉量大，所以產量水平差異不大。

為了解河北山前平原水分虧缺嚴重的詳細情況，首先必須要根據該地區種植制度來進行具體分析。圖6模擬結果表明，小麥年平均產量為6 613.8 kg/hm2，玉米年平均產量為 6 831.1 kg/hm2，玉米年平均產量高于小麥年平均產量，但小麥年平均蒸散量為402.7 mm，多于玉米年平均蒸散量 321.0 mm。小麥生育期內平均水分入滲量為39.5 mm，玉米生育期內平均水分入滲量為132.7 mm，這主要是由于河北山前平原玉米生育期平均降水量為354.7 mm，小麥生育期平均降水量為136.4 mm，河北山前平原降水主要集中在玉米生育期。水分虧損量(水分虧損量=蒸散量-降水量)是衡量作物耗水的重要指標，小麥生育期內平均水分虧損量遠大于玉米生育期內平均水分虧損量(表3)。

由表3可知，玉米生育期內土壤水分大致處于盈余狀態，年平均虧損量為-33.7 mm，只有在枯水年有少量虧損。而無論枯水年還是豐水年，小麥生育期內水分虧損量都比較大，年平均虧損量為266.3 mm，生育期內降水根本無法滿足作物生長需求。小麥耗水量極大，生育期內必須要大量灌溉才能滿足其生長需求，這不利于該地區水資源的可持續利用。因此，調整冬小麥—夏玉米一年兩熟制這一耗水嚴重的種植制度對于河北山前平原地區就顯得十分重要。

表3 小麥玉米不同年型水分虧損量

3 結論

APSIM模型可以較好地模擬河北山前平原冬小麥和夏玉米的產量、蒸散量ET以及土壤水儲量的動態變化，適用于研究分析河北山前平原水資源消耗和產量水平的綜合效應。本研究利用APSIM模型模擬不同年型下的小麥和玉米生產耗水狀況，模擬結果表明，在河北山前平原地區，冬小麥生育期內的降水量遠遠低于冬小麥的耗水量，冬小麥—夏玉米一年兩熟制這一耗水嚴重的種植制度不利于河北山前平原水資源的持續利用，河北山前平原地區種植制度的調整勢在必行。

此外，APSIM模型無法考慮病蟲害的影響，同時由于一些人為因素和自然災害的緣故，實測數據資料出現了偏差，這些均在一定程度上降低了APSIM模型的模擬精度。但從總體上看，這對APSIM模型的整體模擬效果影響不大，APSIM模型在河北山前平原產量水平、水資源消耗的模擬研究上仍然有較好的適應性。