裴宏偉,張紅娟,王飛梟,王 鑫,沈彥俊,馬 宏
(1.河北建筑工程學(xué)院,河北 張家口 075000; 2. 河北省水質(zhì)工程與水資源綜合利用重點(diǎn)實(shí)驗(yàn)室,河北 張家口 075000;3. 張家口市農(nóng)業(yè)高效節(jié)水研究所,河北 張家口 075000; 4. 中國(guó)科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心,河北 石家莊 050021)
水資源匱乏是河北壩上等半干旱地區(qū)作物生長(zhǎng)發(fā)育的主要限制因素[1],同時(shí)也是制約農(nóng)牧業(yè)可持續(xù)發(fā)展的關(guān)鍵資源要素[2]。蒸散(Evapotranspiration,ET)是水在“土壤-作物-大氣”系統(tǒng)中循環(huán)的重要途徑,也是農(nóng)田水分平衡的關(guān)鍵環(huán)節(jié)[3-5]。河北壩上位于河北西北部,處于中國(guó)北方農(nóng)牧交錯(cuò)帶中部核心區(qū)。該地區(qū)年降雨量普遍低于400 mm,且年內(nèi)分配不均勻、年際變化較大[6],水資源不足成為限制本地區(qū)農(nóng)牧業(yè)健康發(fā)展的重要資源因素。隨著北方農(nóng)牧交錯(cuò)帶農(nóng)業(yè)結(jié)構(gòu)調(diào)整工作的不斷深入,特別是國(guó)家層面對(duì)于本地區(qū)發(fā)展高效農(nóng)牧業(yè)的重大戰(zhàn)略需求,裸燕麥作為一種適宜當(dāng)?shù)貧夂?、具有較高經(jīng)濟(jì)價(jià)值的糧飼一體作物,成為本地區(qū)發(fā)展高效農(nóng)牧業(yè)的理想作物種類(lèi)[7-10]。然而,水是農(nóng)業(yè)的命脈,針對(duì)北方農(nóng)牧交錯(cuò)帶寒冷而干旱的氣候條件,如何利用有限的水資源實(shí)現(xiàn)最高作物產(chǎn)量及水分利用效率,已經(jīng)成為本地區(qū)亟待解決的問(wèn)題。
雨養(yǎng)作物的田間水分消耗主要通過(guò)蒸散耗水來(lái)完成,因此,精確估算雨養(yǎng)作物田間蒸散量并明析其變化規(guī)律對(duì)于作物需水管理、旱情預(yù)測(cè)等方面具有重要的理論意義及應(yīng)用價(jià)值。當(dāng)前,國(guó)內(nèi)外學(xué)者對(duì)測(cè)定作物蒸散量的方法進(jìn)行了大量研究,如波文比量平衡法[11]、渦度相關(guān)法[12]、遙感法[13]、大孔徑激光閃爍儀法[14-15]和大型稱(chēng)重式蒸滲儀法[16-17]等。其中,通過(guò)大型稱(chēng)重式蒸滲儀對(duì)作物蒸散量的測(cè)定深受廣大學(xué)者的關(guān)注,且大型稱(chēng)重式蒸滲儀法因具有較高的精度,常被用作檢驗(yàn)其他方法所測(cè)結(jié)果的標(biāo)準(zhǔn)方法[18-20]。自20世紀(jì)50年代以來(lái),已有許多學(xué)者利用大型稱(chēng)重式蒸滲儀對(duì)不同管理方式下的小麥、玉米、水稻、馬鈴薯等作物的田間蒸散耗水規(guī)律等進(jìn)行了大量的研究[21-26]。但是,目前針對(duì)河北壩上地區(qū)雨養(yǎng)條件下裸燕麥的田間蒸散規(guī)律及模擬分析的研究卻鮮有報(bào)道。本文以河北壩上地區(qū)雨養(yǎng)裸燕麥為研究對(duì)象,基于參數(shù)優(yōu)化后的RZWQM2模型,對(duì)2018、2019年裸燕麥生育期內(nèi)田間蒸散動(dòng)態(tài)過(guò)程進(jìn)行模擬,并通過(guò)大型稱(chēng)重式蒸滲儀測(cè)定的試驗(yàn)數(shù)據(jù)對(duì)裸燕麥生育期內(nèi)各時(shí)段的田間蒸散量進(jìn)行深入研究,分析土壤水分及氣象因素對(duì)蒸散的影響,以期為河北壩上地區(qū)農(nóng)田水分高效管理提供理論依據(jù)。
試驗(yàn)于2018年和2019年在張家口市農(nóng)業(yè)高效節(jié)水研究所開(kāi)展。該研究所位于河北省張家口市張北縣(41°15′N(xiāo),114°7′E,海拔約1 393.3 m),多年(1960—2018年)平均降雨量385 mm左右,年平均氣溫為3.2℃,屬于中溫帶大陸性季風(fēng)氣候。試驗(yàn)所需的氣象數(shù)據(jù)由該試驗(yàn)區(qū)所安裝的小型氣象站測(cè)定,主要包括平均風(fēng)速(km·d-1)、最高及最低溫度(℃)、平均相對(duì)濕度(%)、日照時(shí)數(shù)(h)及降雨量(mm)。土壤含水量由智能測(cè)墑系統(tǒng)Insentek Sensor(簡(jiǎn)稱(chēng)“智墑”,東方智感(浙江)科技股份有限公司生產(chǎn))測(cè)定,監(jiān)測(cè)間隔1 h。
RZWQM2模型是由美國(guó)農(nóng)業(yè)部大平原系統(tǒng)研究所開(kāi)發(fā),該模型包括土壤化學(xué)、土壤物理及土壤微生物的過(guò)程,不僅可以較好地模擬作物的生長(zhǎng)過(guò)程,也能模擬日尺度的土壤蒸發(fā)和作物蒸騰。RZWQM2的土壤水平衡模塊采用Green-Ampt方程進(jìn)行土壤水滲透過(guò)程的模擬,同時(shí)使用Richards方程在不同的土層間重新分配水分。而潛在蒸散發(fā)是使用擴(kuò)展的Shuttleworth-Wallace(S-W)模型來(lái)量化裸露或殘余覆蓋的土壤蒸發(fā)和作物蒸騰的每日潛在速率,該方程經(jīng)過(guò)修正,包括空氣動(dòng)力學(xué)和能量通量上的地表作物殘留動(dòng)力學(xué)。
利用裸燕麥2個(gè)完整生育期(2018年和2019年,一年一熟)的實(shí)測(cè)數(shù)據(jù)運(yùn)行S-W模型。所用實(shí)測(cè)數(shù)據(jù)包括氣象數(shù)據(jù)(日均最低及最高氣溫、相對(duì)濕度、降雨量、風(fēng)速等)、作物生長(zhǎng)及田間管理數(shù)據(jù)(株高、葉面積指數(shù)、生物量、產(chǎn)量、播種及收獲時(shí)間、施肥量等)、土壤水分?jǐn)?shù)據(jù)、裸燕麥每日田間蒸散量數(shù)據(jù)以及土壤理化特性參數(shù)。
供試裸燕麥品種為壩莜1號(hào),種植方式為條播,裸燕麥分別于2018年6月15日和2019年6 月5日播種, 2018年9月23日和2019年9月11日收獲,播種密度為450萬(wàn)株·hm-2,種植行距25 cm,播種深度4 cm。裸燕麥在種植前施入復(fù)合肥90 kg·hm-2作為底肥,作物生長(zhǎng)期間沒(méi)有追肥,其他田間管理方式與當(dāng)?shù)乇3忠恢?,裸燕?018—2019年生育期見(jiàn)表1。整個(gè)生育期內(nèi)的降雨量分別為146.5 mm和252.8 mm。
表1 裸燕麥生育期(m-d)
裸燕麥種植區(qū)域安裝有大型稱(chēng)重式蒸滲儀,該蒸滲儀由一個(gè)不銹鋼鐵箱(長(zhǎng)×寬×高為2 m×2 m×2.3 m)組成,精度為32 g(約為0.01 mm水層的深度),其填土土壤物理參數(shù)見(jiàn)表2。裸燕麥生育期內(nèi),在不進(jìn)行灌溉處理的條件下每間隔30 min稱(chēng)取填土重量1次。此外,由于蒸滲儀底部無(wú)滲漏液出流,故可根據(jù)特定時(shí)間段內(nèi)的土體質(zhì)量變化計(jì)算出田間蒸散量。樣方內(nèi)土壤含水率依據(jù)在蒸滲儀內(nèi)安裝的智能測(cè)墑系統(tǒng)進(jìn)行物聯(lián)網(wǎng)自動(dòng)實(shí)時(shí)監(jiān)測(cè)[27],并每間隔7 d左右對(duì)裸燕麥的葉面積、株高及生物量進(jìn)行一次人工測(cè)定。
表2 裸燕麥試驗(yàn)區(qū)土壤物理參數(shù)
對(duì)RZWQM2模型模擬結(jié)果采用絕對(duì)平均誤差(MAE)、決定系數(shù)(R2)和一致性指數(shù)(d)進(jìn)行評(píng)價(jià),計(jì)算公式:
(1)
(2)
(3)
雨養(yǎng)條件下,2018年和2019年裸燕麥田間蒸散量總體擬合趨勢(shì)呈現(xiàn)“單峰型”曲線(圖1),且實(shí)測(cè)值與模擬值均在孕穗期達(dá)到最大,實(shí)測(cè)最大值分別為7.2 mm和7.7 mm,模擬值的最大值為6.6 mm和7.0 mm,偏低了8.33%、9.09%(圖1);從表3可以看出,2018年和2019年裸燕麥田間蒸散量模擬結(jié)果較好,決定系數(shù)(R2)分別為0.71和0.68,一致性指數(shù)(d)分別為0.65和0.63,絕對(duì)平均誤差(MAE)分別為1.11 mm·d-1和1.19 mm·d-1; 2018年裸燕麥生育期內(nèi)總蒸散量的模擬值為358.73 mm,相比實(shí)測(cè)總蒸散量高了10.19 mm,而2019年模擬總蒸散量相對(duì)實(shí)測(cè)總蒸散量低了69.5 mm(表4)。
圖1 RZWQM2模型對(duì)裸燕麥田間蒸散量模擬分析Fig.1 Analysis of evapotranspiration quantity simulation in naked oat field by RZWQM2 model
表3 裸燕麥田間蒸散模擬值與實(shí)測(cè)值的模擬結(jié)果統(tǒng)計(jì)分析
表4 裸燕麥生育期內(nèi)田間蒸散、降雨量及土壤水分變化
2018年和2019年都在生育期前期和后期擬合較好,而在中期擬合偏差相對(duì)較大,可能是因?yàn)樵诼阊帑溕诘那昂笃冢镩g蒸散量都較小,而裸燕麥中期為作物生長(zhǎng)旺盛階段,日均蒸散量較大,在同樣大小的相對(duì)誤差情況下,使得在數(shù)值上體現(xiàn)出裸燕麥生育期前后期絕對(duì)偏差較小,而在生育期中期出現(xiàn)絕對(duì)偏差增大的情況。
由圖2可以看出,裸燕麥整個(gè)生育期內(nèi)孕穗期的日蒸散量最大,而在播種~出苗期的日蒸散量最小,出苗期后日蒸散量呈現(xiàn)逐漸增加的趨勢(shì),在孕穗期達(dá)到峰值,在開(kāi)花期~成熟期日蒸散量又逐步下降。
圖2 裸燕麥生育期內(nèi)田間蒸散量與降雨變化特征Fig.2 Characteristics of evapotranspiration and rainfall during the growth period of naked oat
2018年和2019年裸燕麥生育期內(nèi)逐日蒸散量的變化特征(表5)表明,拔節(jié)期~開(kāi)花期出現(xiàn)日蒸散量最大值,分別為7.20 mm·d-1和7.67 mm·d-1,這一時(shí)段日蒸散量平均值分別為4.50 mm·d-1和4.35 mm·d-1;播種期~出苗期日均蒸散量最小值分別為0.10 mm·d-1和0.25 mm·d-1,該時(shí)段的日蒸散量平均值分別為0.53 mm·d-1和0.52 mm·d-1;當(dāng)裸燕麥在開(kāi)花期~成熟期時(shí),由于作物達(dá)到成熟而吸收水分較少,使得日均蒸散量也隨之減小,其逐日平均蒸散量分別為1.98 mm·d-1和 2.54 mm·d-1。
表5 裸燕麥各生育期田間蒸散情況
選取了裸燕麥孕穗期的兩種典型天氣(多云7月9日和晴天7月15日)對(duì)逐時(shí)蒸散量的變化規(guī)律進(jìn)行分析。裸燕麥田不同天氣條件下蒸散量的變化整體呈現(xiàn)早晚變低、中午左右增高的“單峰型”曲線(圖3)。在晴天和多云都表現(xiàn)為06∶00左右蒸散量開(kāi)始逐漸增大,而兩種不同天氣達(dá)到峰值的時(shí)間不同,晴天峰值出現(xiàn)在15∶00左右,蒸散量達(dá)到了0.73 mm;多云天氣則在12∶00左右蒸散量達(dá)到最大,最大為0.58 mm,可能是因?yàn)樵诙嘣铺鞖庀抡鐣r(shí)分不會(huì)出現(xiàn)裸燕麥氣孔關(guān)閉,與晴天相比出現(xiàn)峰值較早,在峰值過(guò)后蒸散量開(kāi)始降低,到20∶00左右蒸散量開(kāi)始接近于零,其多云夜間平均蒸散量為0.03 mm;晴天由于正午時(shí)刻氣溫較高使得裸燕麥出現(xiàn)“午休”現(xiàn)象,導(dǎo)致氣孔關(guān)閉, 降低了蒸散量。多云和晴天兩種不同天氣的逐時(shí)蒸散量變化過(guò)程表明,裸燕麥田間蒸散量主要在白天,而在夜間蒸散量較低,趨近于零且相對(duì)穩(wěn)定,由于夜間溫度和濕度較低,因此蒸散量較小。
圖3 裸燕麥田典型日蒸散量變化Fig.3 Variation of typical daily evapotranspirationin naked oat field
裸燕麥各生育期的日均蒸散量(圖4)表現(xiàn)為:孕穗期>開(kāi)花期>成熟期>拔節(jié)期>出苗期。出苗期裸燕麥蒸散量小于其他生育期,在本研究中因出苗期葉面積較小而未對(duì)其進(jìn)行觀測(cè)。在拔節(jié)期隨著氣溫和葉面積指數(shù)的逐漸增大,日均蒸散量也較出苗期增大,日均蒸散量為1.65 mm·d-1。隨著裸燕麥的生長(zhǎng)發(fā)育,土壤的蒸發(fā)增大,土壤水分也較拔節(jié)期降低,而裸燕麥在孕穗期的葉面積指數(shù)大于拔節(jié)期,葉面積指數(shù)為4.66,因而蒸散量較大,孕穗期內(nèi)裸燕麥日均蒸散量達(dá)到5.41 mm·d-1,是整個(gè)生育期中蒸散量最大的階段;經(jīng)過(guò)孕穗期的蒸散及大量耗水,土壤含水率降低,因而土壤蒸發(fā)也降低,但因葉面積指數(shù)的增大(LAI=5.07),裸燕麥的蒸騰量增加,對(duì)土壤蒸發(fā)的降低進(jìn)行了彌補(bǔ),因此其日均蒸散量也相對(duì)較大,日均蒸散量為4.91 mm·d-1,相比孕穗期減少了9.24%;而在成熟期葉面積指數(shù)相比孕穗期減少了45%,且該生育期作物逐漸成熟,作物所需水分逐漸減小,因而日均蒸散量減小。在裸燕麥整個(gè)生育期日均蒸散量變化趨勢(shì)與指示作物生長(zhǎng)狀況的葉面積指數(shù)高度一致,同時(shí)受到氣溫等環(huán)境因子的影響(圖4)。在裸燕麥生育初期,由于植株較小(葉面積指數(shù)不足1.5)且氣溫較低,日均蒸散量較低;在裸燕麥生長(zhǎng)旺盛的孕穗期和開(kāi)花期,裸燕麥生長(zhǎng)加快,葉面積指數(shù)迅速增加(4.66~5.07),日均蒸散量均顯著高于前期;到裸燕麥成熟階段,葉片的蒸騰能力減弱,加上土壤濕度和氣溫逐漸降低,因此日均蒸散量相比孕穗期和開(kāi)花期有明顯下降。
圖4 裸燕麥各生育期日均蒸散量、葉面積指數(shù)、土壤含水率及氣溫Fig.4 Daily evapotranspiration, leaf area index, soil moisturecontent and temperature of naked oat at each growth stage
研究基于大型稱(chēng)重式蒸滲儀及RZWQM2模型對(duì)河北壩上地區(qū)2018年和2019年雨養(yǎng)裸燕麥田間蒸散進(jìn)行測(cè)定及模擬,探討了河北壩上地區(qū)裸燕麥整個(gè)生育期和典型日田間蒸散特征,并結(jié)合其土壤水分、作物產(chǎn)量、葉面積指數(shù)、氣溫及降雨等因素進(jìn)行分析,得出以下結(jié)論:
1)RZWQM2模型在通過(guò)擴(kuò)展的S-W ET模型對(duì)2018年和2019年裸燕麥生育期內(nèi)前后期田間蒸散量的擬合程度表現(xiàn)較好,而在中期擬合相差較大;2018、2019年決定系數(shù)(R2)分別為0.71和0.68,平均絕對(duì)誤差(MAE)分別為1.11 mm·d-1和1.19 mm·d-1, 且整體的模擬結(jié)果及趨勢(shì)與實(shí)測(cè)結(jié)果較為一致,均表現(xiàn)為先增大后減小的態(tài)勢(shì)。
2)2018年和2019年裸燕麥生育期內(nèi),田間蒸散量總體呈現(xiàn)先升高后降低的“單峰型”特征曲線,其變化趨勢(shì)與葉面積指數(shù)變化表現(xiàn)為一致的態(tài)勢(shì),日均蒸散量0.52~5.41 mm·d-1,孕穗期日均蒸散量達(dá)到最大,分別為5.02 mm·d-1和5.41 mm·d-1。
3)不同天氣條件下逐時(shí)蒸散量分析表明,裸燕麥田間逐時(shí)蒸散量基本呈現(xiàn)為早晚較低、中午較高的“單峰型”曲線,多云和晴天下的逐時(shí)蒸散變化趨勢(shì)一致,蒸散都主要集中在白天,但多云天氣下蒸散量較低。