基于SIMETAW模型的不同地區(qū)燕麥需水量研究

王沖，林倩，石曉宇，李碩，揭曉婧，褚慶全*

(中國(guó)農(nóng)業(yè)大學(xué)農(nóng)學(xué)院，農(nóng)業(yè)農(nóng)村部農(nóng)作制度重點(diǎn)實(shí)驗(yàn)室，北京 100193；2. 煙臺(tái)市農(nóng)業(yè)技術(shù)推廣中心， 山東 煙臺(tái) 264001)

作物需水量能夠直觀反映作物需水情況[1]，是作物充分灌溉的理論基礎(chǔ)[2]，研究作物需水的時(shí)空變化特征能為優(yōu)化作物布局提供依據(jù)，因此加強(qiáng)作物需水特征研究刻不容緩[3]。SIMETAW模型是由美國(guó)加州大學(xué)戴維斯分校研究開發(fā)的農(nóng)業(yè)水資源管理與規(guī)劃模型，多數(shù)學(xué)者利用SIMETAW模型對(duì)作物需水量計(jì)算展開了大量研究：劉洋[4]計(jì)算了遼寧阜新地區(qū)不同年型下的玉米(ZeamaysL.)蒸散量、需水量、降水耦合度等指標(biāo)，并分析了氣候變化對(duì)玉米水分需求的影響；基于1960—2006年氣象資料，史金麗[5]對(duì)河西走廊地區(qū)的春小麥、春玉米、馬鈴薯(SolanumtuberosumL.)的蒸散量進(jìn)行了模擬，認(rèn)為不同分區(qū)下作物蒸散量的差異達(dá)到顯著水平；劉洋等[6]基于1951—2009年的標(biāo)準(zhǔn)氣象資料計(jì)算了遼寧阜新地區(qū)花生(Arachishypogaea)的蒸散量、有效降水量及灌溉需水量，表明有效降水量呈下降趨勢(shì)，而灌溉需水量呈升高趨勢(shì)。前人對(duì)水稻(Oryzasativa)、小麥(TriticumaestivumL.)、玉米等大宗作物需水量的研究較多，均認(rèn)為作物生育期需水量發(fā)生了變化，且不同地區(qū)的變化存在空間差異[7-11]。

燕麥(AvenaL.)是世界八大谷類作物中營(yíng)養(yǎng)價(jià)值較高的糧飼兼用作物，具有耐寒、耐旱、耐貧瘠等生長(zhǎng)習(xí)性[12]，包括皮燕麥和裸燕麥兩大類。國(guó)外以種植皮燕麥為主，消費(fèi)方式多為飼用，我國(guó)以種植裸燕麥為主，消費(fèi)方式多為食用。我國(guó)90%的春燕麥種植于年降水量200～400 mm的干旱半干旱區(qū)，其水分需求與所處環(huán)境不相適應(yīng)[13]。與玉米、高粱(SorghumbicolorL.)、小麥等作物相比，燕麥被認(rèn)為是耗水量相對(duì)較大的禾谷類作物[14]。華北燕麥區(qū)、北方燕麥區(qū)、西南燕麥區(qū)是我國(guó)的三個(gè)優(yōu)勢(shì)產(chǎn)區(qū)[15-17]，其中燕麥主要分布在我國(guó)華北地區(qū)，種植面積占我國(guó)燕麥總播種面積的80%以上[18]，對(duì)東北地區(qū)而言，在小麥種植面積逐年下降的前提下，燕麥對(duì)東北地區(qū)種植結(jié)構(gòu)調(diào)整的作用不斷上升[19]。但近年來受氣候變化的影響，溫度升高和降水減少的疊加效應(yīng)導(dǎo)致華北和東北地區(qū)的干旱程度不斷加深[20-21]，水分供應(yīng)不足成為制約我國(guó)華北和東北地區(qū)燕麥生產(chǎn)和布局的重要因素，嚴(yán)重影響著當(dāng)?shù)匮帑湹恼ＩL(zhǎng)，探討燕麥生育期需水情況及空間分布顯得尤為重要。

為了解早播和晚播燕麥需水特征，解決農(nóng)作物需水問題，本文基于華北和東北地區(qū)15個(gè)氣象臺(tái)站1960—2014年的標(biāo)準(zhǔn)氣象數(shù)據(jù)，采用SIMETAW模型測(cè)算出近55年來燕麥生育期需水量，解析了該區(qū)域燕麥生育期需水量和有效降水量的空間分布特點(diǎn)，并進(jìn)一步比較了早播和晚播燕麥生育期的降水耦合度，以期為華北和東北燕麥種植區(qū)布局結(jié)構(gòu)的調(diào)整和水資源的合理配置提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)域

華北地區(qū)選取內(nèi)蒙古自治區(qū)、河北省、山西省，其中，內(nèi)蒙古自治區(qū)選取烏蘭浩特、烏拉蓋、赤峰、多倫、阿巴嘎旗、察哈爾右翼前旗、呼和浩特7個(gè)站點(diǎn)，河北省選取豐寧、張家口2個(gè)站點(diǎn)，山西省選取右玉、大同2個(gè)站點(diǎn)。燕麥?zhǔn)莾?nèi)蒙古自治區(qū)農(nóng)牧交錯(cuò)帶的傳統(tǒng)優(yōu)勢(shì)作物；河北省獨(dú)特的氣候條件和優(yōu)良的土壤環(huán)境適合燕麥的生長(zhǎng)；山西省的忻州盆地、大同盆地、呂梁及太行山區(qū)是該省燕麥的主要種植區(qū)[18]。東北地區(qū)包括黑龍江省、吉林省、遼寧省，其中，黑龍江省選取嫩江、哈爾濱2個(gè)站點(diǎn)，吉林省選取白城1個(gè)站點(diǎn)，遼寧省選取沈陽(yáng)1個(gè)站點(diǎn)。燕麥主要種植在黑龍江省的西北部地區(qū)[22]，吉林省進(jìn)行燕麥生產(chǎn)不但增加糧食產(chǎn)量，還一定程度地減少春季揚(yáng)塵天氣[23]，遼寧省的中部及西部地區(qū)是燕麥的種植聚集區(qū)[24]。

1.2 數(shù)據(jù)來源

該研究所用的華北和東北地區(qū)具有代表性的15個(gè)標(biāo)準(zhǔn)氣象站點(diǎn)1960—2014年的逐日氣象數(shù)據(jù)來源于網(wǎng)站信息公開的國(guó)家氣象局。氣象要素主要包括降水量、日照時(shí)數(shù)、平均相對(duì)濕度、大氣壓強(qiáng)、平均風(fēng)速、日最低氣溫、日平均氣溫、日最高氣溫等，各氣象臺(tái)站的地理特征數(shù)值主要包括經(jīng)度、緯度、海拔高度等。作物參數(shù)來源于國(guó)家燕麥?zhǔn)w麥產(chǎn)業(yè)技術(shù)體系的不同生態(tài)區(qū)的田間試驗(yàn)數(shù)據(jù)，主要包括播種面積、密度、播種日期、收獲日期等。

1.3 計(jì)算方法

利用SIMETAW模型分別計(jì)算早播和晚播燕麥生育期需水量, SIMETAW模型的校正與正確使用方法見參考文獻(xiàn)[25-26]。

1.3.1參考作物蒸散量 采用FAO推薦的Penman-Monteith公式[27]計(jì)算。

(1)

式中，ET0為參考作物蒸散量(mm)，Δ為水汽壓對(duì)溫度的斜率(kPa·℃-1)，Rn為到達(dá)作物表面的凈輻射(MJ·m-2·d-1)，G為土壤熱通量密度(MJ·m-2·d-1)，γ為干濕球常數(shù)(kPa·℃-1)，T為作物冠層2 m高處的空氣濕度(℃)，u2為作物冠層2米高處的風(fēng)速(m·s-1)，es為飽和水汽壓(kPa)，ea為實(shí)際水汽壓(kPa)。

1.3.2作物需水量 某時(shí)段作物供水充分情況下，作物系數(shù)與參考作物蒸散量的乘積稱為作物需水量，公式如下。

ETc=Kc×ET0

(2)

式中，ETc為逐日作物需水量(mm)，Kc為作物系數(shù)，ET0為參考作物蒸散量(mm)。

1.3.3有效降水量 作物旱作條件下，除去地表滲漏和徑流到作物根部以下以及淋洗鹽分深層滲漏部分的降雨量外，其余用來滿足作物蒸散所需的降雨量[28]稱為有效降水量，本文利用聯(lián)合國(guó)糧農(nóng)組織推薦的比較作物需水量和降水量值大小的方法[27]確定有效降水量，公式如下。

(3)

式中，ETri為日有效降水量(mm)，ET為日降水量(mm)，ETr為作物全生育期有效降水量(mm)，n為生育期按旬分組的數(shù)量。

1.3.4缺水量 相同生長(zhǎng)階段下，作物需水量與有效降水量之差即為缺水量。若結(jié)果為正值，表明有效降水量沒有達(dá)到作物需水要求，若結(jié)果為負(fù)值，表明有效降水量滿足作物需水要求，若結(jié)果為零值，表明水分供需平衡。

ETaw=ETc-ETr

(4)

式中，ETaw為缺水量(mm)。

1.3.5水分盈虧指數(shù)

不考慮灌溉條件下，水分盈虧指數(shù)(I)為作物生育期有效降水量與需水量的差值占作物需水量的比值[29]，高于0表示水分盈余，低于 0表示水分虧缺，等于 0表示水分供需平衡。

(5)

1.3.6降水耦合度 了解作物需水與自然降水的耦合度，能夠更好地利用自然降水滿足作物需水量，對(duì)高產(chǎn)高效節(jié)水種植制度有重要意義。降水耦合度指在一定時(shí)間、一定地區(qū)內(nèi)，自然降水對(duì)某作物生育期需水的滿足程度，取值介于0～1之間，計(jì)算公式如下。

(6)

式中，λi代表第i階段(月)的作物需水與自然降水的耦合度，Pi代表第i階段(月)內(nèi)自然降水量(mm)，ETci代表第i階段(月)內(nèi)作物需水量(mm)。

全生育期作物需水與自然降水的耦合度(λ)等于各階段(各月)耦合度的以其需水模數(shù)(ETci/ETc)為權(quán)重的加權(quán)平均值。

(7)

2 結(jié)果與分析

2.1 燕麥生育期需水量和有效降水量分布特征

華北和東北地區(qū)的早(晚)播燕麥1960—2014年生育期需水量和有效降水量平均值如圖1所示，研究區(qū)域各站點(diǎn)早播燕麥生育期需水量在305.13～362.88 mm之間，平均值為334.73 mm，晚播燕麥在279.60～348.92 mm之間，平均值為319.13 mm。早播和晚播燕麥生育期需水量空間分布不同，絕大多數(shù)站點(diǎn)早播燕麥生育期需水量明顯高于晚播燕麥(阿巴嘎旗、呼和浩特和張家口除外)。早播燕麥生育期需水量的高值區(qū)(>340 mm)主要分布在白城、多倫、嫩江、察哈爾右翼前旗、赤峰，低值區(qū)(<320 mm)主要分布在豐寧、右玉、烏蘭浩特、張家口；而嫩江、察哈爾右翼前旗、阿巴嘎旗、大同、哈爾濱(>330 mm)為晚播燕麥生育期需水量的高值區(qū)，沈陽(yáng)、右玉、豐寧、白城為晚播燕麥生育期需水量的低值區(qū)(<310 mm)。

早播和晚播燕麥生育期有效降水量波動(dòng)范圍分別為84.24～224.43 mm和120.02～236.30 mm，平均值分別為132.70 mm和173.07 mm，且早(晚)播燕麥生育期有效降水量大于平均值的站點(diǎn)比重均為47%。早(晚)播燕麥生育期有效降水量空間分布特征相似，均呈現(xiàn)明顯的經(jīng)向和緯向分布，表現(xiàn)為由北向南遞減，由東向西遞減，高值區(qū)主要位于嫩江、沈陽(yáng)、哈爾濱、豐寧等地，低值區(qū)主要集中在察哈爾右翼前旗、張家口、阿巴嘎旗等地。上述研究結(jié)果表明，在不考慮灌溉條件的前提下，燕麥生育期有效降水量均不能滿足燕麥的需水要求，從降水量與需水量的空間分布來看，降水量低值區(qū)大多是燕麥需水量的高值區(qū)，因此，各地區(qū)應(yīng)結(jié)合當(dāng)?shù)匮帑湹男杷徒邓卣?，積極采取相應(yīng)措施以充分利用自然降水。

2.2 燕麥生育期水分盈虧指數(shù)比較分析

研究區(qū)域各站點(diǎn)均表現(xiàn)出不同程度的水分虧缺(表1)，早播燕麥生育期缺水量在130.67～258.48 mm之間，平均值為202.03 mm；水分盈虧指數(shù)波動(dòng)范圍在-0.75～-0.37之間，平均值為 -0.60，察哈爾右翼前旗、阿巴嘎旗、張家口水分虧缺最為嚴(yán)重，水分盈虧指數(shù)分別為-0.75、-0.75、-0.72，而嫩江、沈陽(yáng)、哈爾濱水分虧缺程度較小，其值分別為-0.37、-0.46、-0.49；晚播燕麥生育期缺水量在89.84～219.54 mm之間，平均值為144.49 mm，水分盈虧指數(shù)在-0.65～-0.30之間，平均值為-0.46，與早播燕麥相比，晚播燕麥水分虧缺程度較小，說明晚播燕麥能充分利用有效降水量，提高降水資源的利用率。在空間分布上，與早播燕麥分布類似，晚播燕麥水分虧缺最嚴(yán)重的地區(qū)在阿巴嘎旗、察哈爾右翼前旗、張家口，水分虧缺指數(shù)分別為-0.65、-0.60、-0.59，豐寧、嫩江、哈爾濱、沈陽(yáng)水分虧缺程度較小，其值均不小于-0.35。通過水分盈虧指數(shù)的分析得出，各區(qū)域站點(diǎn)燕麥的水分盈虧指數(shù)均小于0，表現(xiàn)出燕麥生長(zhǎng)期間所需的水分不足，遭遇不同程度的干旱脅迫，且不同地點(diǎn)的燕麥?zhǔn)艿降乃痔澣贝笮〔煌?/p>

表1 1960—2014年華北和東北地區(qū)燕麥生育期缺水量與水分盈虧指數(shù)Table 1 Water deficit and moisture loss index of oat during growing season in North and Northeast China from 1960 to 2014

2.3 燕麥生育期降水耦合度分析

分析早(晚)播燕麥全生育期平均降水耦合度發(fā)現(xiàn)，晚播燕麥的降水耦合度高于早播燕麥，早播燕麥介于0.26～0.68之間，晚播燕麥介于0.38～0.78之間(圖2)。從生長(zhǎng)季內(nèi)不同月份看，早播燕麥生育期在4—7月，各月份站點(diǎn)平均降水耦合度7月> 4月> 6月> 5月，表明在5—6月時(shí)自然降水對(duì)燕麥需水的滿足程度較小，在基本成熟期時(shí)(7月)自然降水對(duì)燕麥需水的滿足程度達(dá)到最大，晚播燕麥生育期基本在5—8月，雖然生長(zhǎng)前期階段(5—6月)燕麥需水少，但降水耦合度低，自然降水不能滿足燕麥的需水，在生長(zhǎng)后期階段(7—8月)降水耦合度升高，一定程度上可滿足燕麥對(duì)水分的需求。從區(qū)域上看，東北地區(qū)燕麥的降水耦合度高于華北地區(qū)。其中，早播燕麥降水耦合度，東北地區(qū)在0.41～0.68之間，華北在0.26～0.50之間，而晚播燕麥降水耦合度，東北地區(qū)在0.50～0.78之間，華北在0.38～0.72之間，說明東北地區(qū)的燕麥對(duì)自然降水利用率較高，再通過調(diào)整作物布局結(jié)構(gòu)，可達(dá)到對(duì)水資源的優(yōu)化利用。從各站點(diǎn)看，15個(gè)站點(diǎn)早播和晚播燕麥的降水耦合度相差不大，表現(xiàn)出東高西低、北高南低的空間特征，如高值區(qū)的哈爾濱、嫩江、沈陽(yáng)的早(晚)播燕麥降水耦合度均達(dá)到了0.50以上，而低值區(qū)的張家口、察哈爾右翼前旗、阿巴嘎旗的降水耦合度在0.40左右。上述結(jié)果表明，由于降水時(shí)期與燕麥需水階段的不匹配造成了生長(zhǎng)前期階段的降水耦合度較低，使部分降水不能充分利用，在生長(zhǎng)后期階段自然降水量升高和需水量降低的雙重效應(yīng)導(dǎo)致了燕麥降水耦合度升高。雖然燕麥整體降水耦合度不高，但晚播燕麥全生育期降水耦合度大于早播燕麥，能更大限度地利用降水資源，因此華北和東北地區(qū)可適當(dāng)推遲燕麥播期，提高降水耦合度，充分利用自然降水。

圖2 1960—2014年華北和東北地區(qū)燕麥生育期多年逐月平均降水耦合度Fig.3 Average coupled rainfall degree of oat during growing season in North and Northeast China from 1960 to 2014

3 討論

燕麥主要分布在我國(guó)東北、華北、西北、西南等高寒冷涼、氣候干旱、土壤貧瘠的地區(qū)，其中農(nóng)業(yè)氣候因子的異常變化加劇了華北和東北地區(qū)的生態(tài)脆弱性，過去60年間東北地區(qū)的平均氣溫以每10年增加0.33 ℃的速率不斷上升[30]，可利用降水量表現(xiàn)出下降的趨勢(shì)[31]，而華北地區(qū)是我國(guó)降水量減少最多的區(qū)域[32]。燕麥被認(rèn)為是水分利用率較低的作物[33]，華北和東北地區(qū)一系列的干旱問題成為限制燕麥生長(zhǎng)的重要因素，如何有效利用當(dāng)?shù)氐淖匀唤邓畞頋M足燕麥的需水要求成為本研究的重點(diǎn)內(nèi)容。本文對(duì)比了研究區(qū)域燕麥生育期需水量和有效降水量的空間分布，降水量低值區(qū)多數(shù)是需水量高值區(qū)，表明該區(qū)域水分虧缺比較嚴(yán)重，進(jìn)一步利用缺水量和水分盈虧指數(shù)分析了研究區(qū)域有效降水量對(duì)燕麥生育期需水的盈虧狀況，各站點(diǎn)均受到干旱脅迫的影響，生育期降水量均不能滿足燕麥的用水需求，該指標(biāo)可定量評(píng)估降水對(duì)作物需水的滿足程度[34]。燕麥大多種植在水資源不足或生態(tài)環(huán)境脆弱的條件下，在華北和東北地區(qū)積極發(fā)展燕麥產(chǎn)業(yè)的過程中，應(yīng)充分考慮當(dāng)?shù)氐墓?、溫、水等氣候資源，加強(qiáng)節(jié)水農(nóng)作制的推廣，規(guī)劃燕麥布局以合理分配水資源。

燕麥的不同生長(zhǎng)階段對(duì)水分需求不同，分蘗期至抽穗期耗水量最大，約占全生育階段需水量的70%，其次是灌漿期至成熟期，占總需水量的20%左右[35]。燕麥對(duì)干旱較為敏感，在需水量較大的生育階段必需具有足夠的水分，才得以避免干旱缺水對(duì)燕麥可能造成的脅迫，保證燕麥的正常發(fā)育。本文對(duì)燕麥的降水耦合度進(jìn)行了探討，深入了解降水對(duì)燕麥不同時(shí)期需水的滿足情況，為研究區(qū)域燕麥合理利用水資源提供參考。華北和東北地區(qū)7、8月份自然降水較多，早播燕麥的生育期基本在4—7月，晚播燕麥在5—8月。對(duì)早播燕麥而言，生育期內(nèi)5月份需水量較大，但降水遠(yuǎn)遠(yuǎn)不能滿足需水要求，導(dǎo)致降水耦合度較低，各站點(diǎn)平均降水耦合度約為0.28；7月份基本達(dá)到成熟期時(shí)，降水量增加，而生長(zhǎng)需水減少，致使耦合度較高，基本在0.90以上。對(duì)晚播燕麥而言，6月份是生育期需水較大的時(shí)期，但降水耦合度較低，各站點(diǎn)平均降水耦合度約為0.45；8月份達(dá)到成熟期時(shí)，需水量較少，但降水量大，導(dǎo)致降水耦合度較高，接近于1.00。上述研究表明，當(dāng)燕麥處于拔節(jié)期或孕穗期時(shí)，需水較多，自然降水卻少，降水耦合度低；當(dāng)處于基本成熟期時(shí)，需水較少，自然降水卻多，造成水資源不能充分利用。為避免燕麥需水較多的生長(zhǎng)階段與當(dāng)?shù)亟邓袝r(shí)間的不匹配，華北和東北地區(qū)可適當(dāng)推遲燕麥播期，提高降水耦合度，充分利用自然降水。

燕麥的生長(zhǎng)狀況和發(fā)育進(jìn)程受播期的影響，適時(shí)播種可充分利用當(dāng)?shù)氐淖匀唤邓凸鉄豳Y源，減輕低溫干旱對(duì)燕麥生長(zhǎng)前期的危害[36]，燕麥的適宜播期應(yīng)由需水較大的灌漿期與此時(shí)的氣候條件共同決定[37]，氣候溫暖的地區(qū)適宜早播，而氣候冷涼的地區(qū)適宜晚播，有助于提高燕麥的結(jié)實(shí)率。本文中研究區(qū)域的氣候條件以冷涼為主，近年來干旱程度不斷加劇，且晚播燕麥的降水耦合度高于早播燕麥，為了適應(yīng)當(dāng)?shù)氐臍夂驐l件，應(yīng)適當(dāng)推遲燕麥播期，達(dá)到需水量與降水量時(shí)期的一致性，減輕干旱脅迫的影響，從而有效利用自然降水，提高降水耦合度。這一結(jié)論與李鳳霞等[38]的結(jié)論一致，其研究表明在沒有灌溉條件的地區(qū)燕麥不宜早播，若適時(shí)晚播可及時(shí)遇到較大的降水，對(duì)燕麥生長(zhǎng)發(fā)育起到積極的作用。在氣候干旱的華北和東北地區(qū)，為了使燕麥生長(zhǎng)過程中得到充足的水分，在土壤墑情允許的條件下，應(yīng)適當(dāng)推遲燕麥播期，更好地利用自然降水，為當(dāng)?shù)毓?jié)水農(nóng)業(yè)的推廣提供參考依據(jù)。