王鵬新 劉麗娜 劉峻明 胡 新

(1.中國農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院， 北京 100083； 2.農(nóng)業(yè)部農(nóng)業(yè)災(zāi)害遙感重點實驗室， 北京 100083；3.商丘市農(nóng)林科學(xué)院小麥研究所， 商丘 476000)

0 引言

冠層氣溫是作物冠層內(nèi)空氣本身的溫度，屬于農(nóng)田小氣候中氣象因素,反映農(nóng)作物的生長環(huán)境，對作物生長發(fā)育、產(chǎn)量、凍害和病蟲害的發(fā)生有重要影響[1]。晚霜凍害是一種突發(fā)性的低溫農(nóng)業(yè)氣象災(zāi)害，一般發(fā)生于冬小麥拔節(jié)后，常采用氣溫指標(biāo)法、高光譜分析法等進行監(jiān)測[2-3]。冠層氣溫與發(fā)生凍害時作物的生長狀況緊密聯(lián)系，但由于受太陽輻射、湍流交換、風(fēng)速等影響而頻繁變化，了解其動態(tài)變化具有一定難度。一種相對可行的解決方法是通過綜合考慮農(nóng)田作物生態(tài)系統(tǒng)和水熱運移過程之間相互作用的數(shù)值模型實現(xiàn)對作物生長過程和農(nóng)田小氣候要素的模擬。

FLERCHINGER等[4]建立的SHAW模型，在1991年加入了蒸騰作用和作物冠層，模擬了土壤-植被-大氣系統(tǒng)中能量流動和物質(zhì)循環(huán)過程，在氣象、水文、水土保持等方面得到諸多應(yīng)用[5-6]。此外，SHAW模型也被用于對玉米、小麥等農(nóng)作物冠層內(nèi)的小氣候要素進行模擬，并獲得了較好效果[7-8]。劉峻明等[9]利用SHAW模型模擬冬小麥拔節(jié)期后近地面層不同高度的每小時氣溫變化特征，驗證了模型對麥田氣溫模擬的有效性，并認(rèn)為該模型在遭受晚霜凍時可以準(zhǔn)確地模擬麥田溫度廓線。SHAW模型是一個分層系統(tǒng)，其運行依賴于能準(zhǔn)確描述垂直方向植株結(jié)構(gòu)的參數(shù)。對農(nóng)作物而言，這些參數(shù)主要包括株高、葉寬、葉面積指數(shù)(LAI)、地上干生物量等，但這些參數(shù)較難獲取，通常采用定期采樣的方式獲得，因而限制了SHAW模型的大范圍應(yīng)用。作物生長模型是一種面向過程、機理性的動態(tài)模型，它根據(jù)氣象、土壤、作物品種、作物種植因素，模擬作物生長、發(fā)育過程和產(chǎn)量，其優(yōu)勢主要體現(xiàn)在，以日為時間步長動態(tài)地模擬作物生長發(fā)育過程，能定量地描述作物的基本生理生態(tài)過程和植株結(jié)構(gòu)相關(guān)參數(shù)。目前已經(jīng)建立的作物生長模型多達(dá)百種，其中WOFOST模型可以通過改變作物參數(shù)而用于不同作物種類或品種，其普適性已經(jīng)得到了肯定[10]，廣泛用于作物估產(chǎn)、確定生長關(guān)鍵因子、物候監(jiān)測等方面[11-13]。

本文耦合WOFOST模型和SHAW模型，利用WOFOST模型模擬冬小麥生長過程，以其輸出結(jié)果中的若干生物指標(biāo)作為反映冬小麥結(jié)構(gòu)的直接或間接參數(shù)，并將這些結(jié)構(gòu)參數(shù)作為SHAW模型的輸入，對農(nóng)田小氣候環(huán)境要素廓線進行模擬，以解決SHAW模型植株結(jié)構(gòu)參數(shù)獲取困難的問題，以期為晚霜凍害監(jiān)測指標(biāo)提供參考依據(jù)。

1 實驗與方法

1.1 研究區(qū)域概況

研究區(qū)商丘市位于河南省最東部，地理位置為114°49′～116°39′ E，33°43′～34°52′ N，地處黃河沖積平原東南部，屬于淮河流域。商丘市屬暖溫帶季風(fēng)氣候，四季分明，年均溫在13～15℃，年降水量在650～780 mm，年日照時數(shù)2 205～2 428 h，年均蒸發(fā)量1 461 mm，無霜期221 d[8]。商丘市土壤類型為潮土，呈中性至弱酸性。冬小麥播種時間一般為10月上中旬，12月中旬進入越冬期，次年2月下旬進入返青期，拔節(jié)期在3月中下旬，4月中下旬抽穗開花，5月底或6月初小麥成熟收割。

1.2 田間實驗與數(shù)據(jù)采集

在河南省商丘市農(nóng)林科學(xué)研究所雙八實驗站進行田間實驗，位于北緯34°31′55″，東經(jīng)115°42′37″，海拔高度50.1 m。實驗田冬小麥品種為“豫麥18號”，播種期為2015年10月20日，次年3月8日進入拔節(jié)期。2016年3月12日—5月29日在實驗田中設(shè)置2個樣點進行定位觀測。

(1)冬小麥冠層各高度空氣溫度的觀測。在實驗田土壤表面上方不同高度安裝空氣溫度傳感器，高度分別為5、10、20、30、40、50、60、70、80 cm，每隔1 h觀測1次氣溫。

(2)上墊面氣象數(shù)據(jù)的觀測。在實驗田土壤表面上方2 m高度處安裝傳感器，每隔1 h觀測1次空氣溫度、空氣相對濕度、風(fēng)速、太陽輻射強度和降水量。取1 d內(nèi)空氣溫度觀測值中的最大值和最小值分別作為日最高氣溫和日最低氣溫，取1 d內(nèi)風(fēng)速測量值的平均值作為日平均風(fēng)速，取1 d內(nèi)太陽輻射量的累計值作為日太陽輻射量，取前一日20點至當(dāng)日20點之間降水量的累計值作為當(dāng)日降水量，取每日08:00空氣相對濕度與飽和水汽壓的乘積作為每日08:00水汽壓。另外，2015年冬小麥播種期至2016年田間觀測之前的時間段內(nèi)的每日氣象數(shù)據(jù)均由商丘市氣象站提供。

(3)土壤溫濕度的觀測。在土壤不同深度處埋藏土壤溫度和土壤水分傳感器，深度分別為0、5、10、20、30、40 cm，每隔1 h觀測1次土壤溫度和濕度。

(4)取樣觀測。在觀測時段內(nèi)，每隔4 d對株高、葉寬、生長錐或幼穗所在高度、葉面積指數(shù)(LAI)、地上各部分的干生物量進行測定。在拔節(jié)期和孕穗期，測量植株基部到最上部展開葉的葉尖或者抽穗的距離記為株高；在抽穗期到成熟期過程中，測量從植株基部到穗頂(不包括芒)的距離記為株高。葉寬為葉片中最寬距離，測量主莖所有葉片的葉寬并計算平均值記為植株的葉寬。

(5)冬小麥生育期的調(diào)查界定。田間觀察確定實驗田中冬小麥進入關(guān)鍵物候期(拔節(jié)-抽穗-開花-灌漿-成熟)的時間。

1.3 WOFOST模型本地化標(biāo)定

WOFOST模型根據(jù)氣象和土壤條件，以作物生理生態(tài)過程為基礎(chǔ)對一年生作物的根、葉、莖、穗生物量及土壤水分動態(tài)進行模擬[14]。其模擬過程主要通過氣候、作物、土壤3個模塊協(xié)同完成，其根據(jù)需要分別在潛力模式、水分限制模式、養(yǎng)分限制模式等3種模式下模擬作物生長過程。研究使用實驗田觀測或者商丘市氣象站提供的降水資料作為輸入，不考慮灌溉，開展水分限制模式的冬小麥生長模擬。

為使WOFOST模型準(zhǔn)確模擬冬小麥的生長發(fā)育過程，需要對模型參數(shù)進行標(biāo)定。研究中對冬小麥的關(guān)鍵生長參數(shù)和土壤參數(shù)進行調(diào)整。對關(guān)鍵生長參數(shù)，先通過查閱文獻(xiàn)、WOFOST模型推薦、實驗測定等方式確定可能取值范圍，再運用FSEOPT優(yōu)化程序結(jié)合田間實測值確定最終參數(shù)值[14-16]。土壤參數(shù)主要實測了田間持水量、凋萎系數(shù)、飽和含水率等。主要參數(shù)見表1。

表1 WOFOST模型主要生長參數(shù)Tab.1 Main growth parameters of WOFOST

注：各器官的維持呼吸作用速率以CH2O計。

1.4 SHAW模型描述與本地化標(biāo)定

SHAW模型用于模擬從植物冠層、雪被、枯落物、地表到土壤一定深度的一維垂直廓線[17-19]。模型對大氣-植物-土壤系統(tǒng)中物質(zhì)能量傳輸過程有著清晰的數(shù)學(xué)表達(dá)，綜合考慮冠層、雪、枯落物和土壤中的熱、水、溶質(zhì)通量而得到一個同步解作為系統(tǒng)的解[20-21]。該系統(tǒng)在垂直方向進行分層，系統(tǒng)中能量和水分輸入分別由上邊界的氣象條件和下邊界土壤條件決定。

SHAW模型的作物層基于能量平衡方程計算水分的吸收和蒸發(fā)，基于K理論求解層內(nèi)溫濕度廓線。假定葉片熱容量可忽略，作物層每個節(jié)點需滿足以下的葉能量平衡方程

Rn+Hl+LvEl=0

(1)

其中

(2)

(3)

式中Rn——作物層吸收的全波段凈輻射，W/m2Hl——作物層間從葉片傳輸?shù)娇諝鉄崃浚琖/m2El——作物層間水汽通量，kg/(m2·s)LAI——葉面積指數(shù)ρa——空氣密度，kg/m3ca——空氣比熱容，J/(kg·K)Ta、Tl——氣溫、葉溫，℃ρva、ρvs——空氣的水汽密度、葉片氣孔腔內(nèi)水汽密度，kg/m3

Lv——蒸發(fā)潛熱，J/kg

rh、rv、rs——熱量傳輸阻力、水汽傳輸阻力、氣孔阻力，s/m

模型中植物生物物理特征參數(shù)、土壤理化參數(shù)和水力特性參數(shù)主要來自文獻(xiàn)[8-9]或者基于實測數(shù)據(jù)計算獲得。具體參數(shù)見表2。

表2 SHAW模型站點主要特征參數(shù)Tab.2 Main characteristic parameters of SHAW at experimental site

注：臨界葉水勢(PLEAF0)指氣孔阻力為最小值的2倍時的葉水勢。Ks、ψe和b通過Campbell提出的利用土壤結(jié)構(gòu)、容重、顆粒組成等建立的經(jīng)驗方程計算獲得。

1.5 WOFOST模型與SHAW模型的耦合

WOFOST模型以日為時間步長運行，輸入每日氣象數(shù)據(jù)(氣溫、風(fēng)速、太陽輻射、水汽壓和降水等)模擬冬小麥生長過程，獲得每日各器官干物質(zhì)積累量、LAI等數(shù)據(jù)，并將這些數(shù)據(jù)作為SHAW模型的輸入?yún)?shù)。SHAW模型所需要的株高與葉寬數(shù)據(jù)無法直接通過WOFOST模型獲得，本研究分別基于單位面積的莖干質(zhì)量與LAI建立轉(zhuǎn)換模型來獲得。

冬小麥拔節(jié)之后莖稈伸長，莖稈壁增厚且直徑變粗，從而株高和莖粗不斷增大，而抽穗開花之后設(shè)定株高保持不變。在此過程中冬小麥的單位面積莖干質(zhì)量也不斷增大。應(yīng)用2016年冬小麥拔節(jié)期-開花期實測的H(株高，m)和Wst(單位面積莖干質(zhì)量，kg/m2)，構(gòu)建的兩者間的回歸方程為

(4)

方程達(dá)到極顯著水平。應(yīng)用此方程將WOFOST模型輸出的單位面積莖干質(zhì)量數(shù)據(jù)轉(zhuǎn)換為株高，并作為SHAW的輸入。

在拔節(jié)期和孕穗期，冬小麥葉寬、葉長和LAI均不斷增大。抽穗期葉寬繼續(xù)增大，而LAI逐漸增至最大值后減小。抽穗期之后葉寬基本保持不變。將在2016年冬小麥拔節(jié)期-抽穗期實測LAI與w(葉寬，cm)建立相關(guān)方程為

LAI=-24.625w2+77.073w-53.783

(R2=0.92，P<0.001)

(5)

方程達(dá)到極顯著水平。應(yīng)用此方程將WOFOST模型輸出的LAI數(shù)據(jù)轉(zhuǎn)換為葉寬，并作為在拔節(jié)-抽穗期間SHAW的輸入。

將以上各參數(shù)值、實驗田上墊面的實測氣象數(shù)據(jù)和土壤溫濕度數(shù)據(jù)等輸入SHAW模型，驅(qū)動SHAW模型運行。SHAW模型以小時為時間步長，輸出冠層垂直方向各高度層間的能量、水分和溶質(zhì)通量以及描述作物冠層小氣候環(huán)境的冠層空氣溫度等數(shù)據(jù)。

1.6 模擬結(jié)果的精度評價指標(biāo)

1.6.1 模型效率

模型效率(Model efficiency，ME)用于表達(dá)氣溫整體模擬效果。ME的值越接近1時，表明模擬效果越好。ME的計算公式為

(6)

1.6.2 均方根誤差

均方根誤差(Root mean square error, RMSE)用來反映氣溫模擬值與實測值之間的總體差異。RMSE越接近于0，說明模擬誤差越小，精度越高。RMSE的計算公式為

(7)

1.6.3 平均偏差

平均偏差(Mean bias error, MBE)用來反映氣溫模擬值與實測值之間的平均偏差，負(fù)值表示氣溫被低估，正值表示氣溫被高估。MBE越接近于0，說明模擬精度越高。MBE的計算公式為

(8)

2 結(jié)果與分析

2.1 冬小麥冠層各高度的氣溫模擬

圖1 冬小麥不同高度層氣溫模擬值和實測值對比(2016年3月12日—5月29日)(*** 表示0.001顯著水平)Fig.1 Comparisons of simulated and measured air temperatures at different heights within winter wheat canopy from March 12 to May 29, 2016

利用WOFOST-SHAW模型對實驗田冬小麥2016年3月12日(拔節(jié)期)—5月29日(成熟期)(儒略日72-150)冬小麥冠層5、10、20、30、40、50、60、70、80 cm高度的每小時氣溫進行模擬，并通過氣溫實測值驗證模型的模擬結(jié)果(圖1)?？梢钥闯觯趯痈鞲叨葰鉁啬M值與實測值均集中分布在1∶1線附近，兩者具有顯著的線性相關(guān)關(guān)系，各高度層的R2均不小于0.93(P<0.001)，表明WOFOST-SHAW模型的氣溫模擬值和實測值較一致。分析圖1中氣溫模擬值和實測值的統(tǒng)計指標(biāo)，發(fā)現(xiàn)在除80 cm的高度上，WOFOST-SHAW模型的ME均大于0.90，RMSE均在1.41～2.11℃之間。在5～40 cm高度上，WOFOST-SHAW模型的MBE值均為正值，模擬值偏大0.51～0.93℃。在50～80 cm高度上，WOFOST-SHAW模型的MBE值均為負(fù)值，模擬值偏小0.21～1.38℃。這些結(jié)果表明，WOFOST-SHAW模型的模擬誤差在合理誤差范圍內(nèi)，整體模擬精度較高。另外，冠層上部(60～80 cm高度)和冠層下部(5～20 cm高度)的ME分別介于0.85～0.94、0.91～0.93之間，而冠層中部(30～50 cm高度)ME為0.95，冠層中部的ME值高于冠層上部和冠層下部；冠層上部、中部、下部的RMSE分別介于1.77～3.09℃、1.41～1.53℃、1.86～1.89℃之間，MBE分別介于-1.38～-0.21℃、-0.25～0.55℃、0.78～0.93℃之間，表明冠層中部的模擬誤差最小，模擬效果優(yōu)于冠層上部和冠層下部。

圖2 冬小麥冠層不同高度氣溫模擬值的絕對誤差頻率分布Fig.2 Frequency distribution of absolute errors of simulated air temperature at different heights within winter wheat canopy

為進一步進行模擬結(jié)果的精度評價，計算WOFOST-SHAW模型對冬小麥冠層各高度氣溫模擬值與實測值的差值的頻率分布(圖2)?？煽闯觯鞲叨饶M結(jié)果絕對誤差頻率分布圖形狀相似，各高度絕對誤差分布寬度基本一致，均在[-9.5,7.5]℃范圍內(nèi)，絕對誤差頻率分布直方圖峰值基本上均落在[-0.5,0.5]℃之間。在[-9.5, -0.5]℃誤差范圍內(nèi)的冬小麥冠層上部的頻率高于冠層下部，而在[0.5, 7.5]℃誤差范圍內(nèi)，則結(jié)果相反，說明冠層上部模擬值整體偏低，冠層下部模擬值整體偏高，這與之前所述的不同高度MBE值反映的模擬偏差的結(jié)果一致。WOFOST-SHAW模型在冠層各高度的氣溫模擬誤差主要集中分布于[-1.5, 2.5]℃之間，在[-1.5, -0.5]℃、[-0.5, 0.5]℃、[0.5, 1.5]℃和[1.5, 2.5]℃范圍內(nèi)的比例分別為15%、26%、23%、12%左右。說明超過75%的氣溫模擬值誤差在[-1.5, 2.5]℃范圍內(nèi)，WOFOST-SHAW模型在冬小麥冠層各高度的氣溫模擬值精度較好。

圖3 冬小麥不同高度日最低氣溫模擬值與實測值對比Fig.3 Comparisons of simulated and measured daily minimum air temperatures at different heights within winter wheat canopy

2.2 冠層日最低氣溫的模擬

冬小麥在拔節(jié)到抽穗期間遇到0℃以下低溫時易造成晚霜凍害，日最低氣溫常作為霜凍害監(jiān)測指標(biāo)[22]。相對于氣象站1.5 m百葉箱內(nèi)氣溫，冠層不同高度的氣溫能更真實地反映冬小麥生長環(huán)境，因此分析WOFOST-SHAW模型對冬小麥冠層不同高度日最低氣溫的模擬情況有利于了解霜凍害發(fā)生的冬小麥冠層低溫環(huán)境。在冬小麥冠層各高度分別取1 d內(nèi)氣溫模擬值和實測值的最小值作為冠層日最低氣溫的模擬值與實測值，并進行對比(圖3，由于實驗田斷電，儒略日84、85、124、146、147的氣象數(shù)據(jù)等采集不完整，造成SHAW模型在當(dāng)天無法正常運行模擬，從而模擬值和實測值均缺失)。可以發(fā)現(xiàn)，在觀測時間段內(nèi)冠層各高度日最低氣溫的模擬值與實測值基本一致。模型對3月14日、3月27日、4月4日、4月29日、5月16日(儒略日74、87、95、120、137)出現(xiàn)的日最低氣溫極小值和3月17日、4月2日、4月22日、5月5日、5月22日(儒略日77、93、113、126、143)出現(xiàn)的極大值都較好地體現(xiàn)了，說明模擬值能較好地體現(xiàn)出冬小麥冠層日最低氣溫的波動。整體看冠層各高度日最低氣溫的模擬值與實測值擬合度較高，且能反映實測值的變化趨勢。

為進行冠層日最低氣溫模擬結(jié)果的精度評價，計算其模擬值和實測值的各個統(tǒng)計指標(biāo)(表3)?？梢钥闯?，模型在冠層各高度層日最低氣溫的ME平均值為0.86，RMSE和MBE的平均值分別為1.44℃、

表3 冬小麥冠層不同高度日最低氣溫模擬結(jié)果統(tǒng)計分析Tab.3 Statistical analysis of simulated daily minimum air temperature at different heights within winter wheat canopy

0.53℃，表明模型對冠層日最低氣溫整體模擬效果比較合理。WOFOST-SHAW模型在5、10、20、30 cm高度的ME均不小于0.92，40、50 cm高度的ME在0.83～0.86之間，而60、70、80 cm高度的ME下降至0.74～0.79之間，說明模型的ME隨高度增加而下降；同時RMSE總體上隨高度上升而增加，距離冠層頂部越近，模擬誤差越大；模型在各高度的MBE基本上均大于0，各層模擬值稍高于實測值?？傊?，WOFOST-SHAW模型可以較好地模擬冠層日最低溫度，存在隨冠層高度增加模擬精度下降、日最低氣溫模擬值偏高的問題。

進一步分析WOFOST-SHAW模型對零下冠層日最低氣溫的模擬效果。在3月12日—3月16日、3月27日(處于拔節(jié)期)在實驗田上方2 m高度處均觀測到0℃以下氣溫，而在3月12日—27日時間段內(nèi)對低溫敏感的冬小麥生長錐的所在高度由3 cm增至9 cm。分別以5 cm和10 cm高度為例，分析WOFOST-SHAW模型在低溫天氣出現(xiàn)時冬小麥易受凍高度層的氣溫模擬效果(圖4)，發(fā)現(xiàn)在3月12日—3月16日、3月27日這6天模型的整體氣溫模擬值與實測值比較接近，且變化趨勢一致。在3月14日、16日和27日模型對0℃以下冠層日最低氣溫的模擬值與實測值基本一致，而3月12日、13日和15日模型對日最低氣溫的模擬值高于實測值，但仍然可以反映冬小麥生長錐所在位置低于0℃的環(huán)境氣溫。

圖4 冬小麥冠層5 cm和10 cm高度氣溫模擬值與實測值對比Fig.4 Comparisons of simulated and measured air temperatures at 5 cm and 10 cm height within winter wheat canopy

3 討論

SHAW模型為分層結(jié)構(gòu)，可以獲得垂直方向的溫度分布，因此植株結(jié)構(gòu)參數(shù)對于SHAW模型有重要影響；而由于作物在不斷生長，所以植株結(jié)構(gòu)參數(shù)處于動態(tài)變化中，而傳統(tǒng)的田間采樣方式存在周期長、破壞性大、費時費力的缺點，這成為SHAW模型的限制因素。WOFOST模型的優(yōu)勢在于為SHAW模型提供動態(tài)的作物生長信息，基于這些信息轉(zhuǎn)換為植株結(jié)構(gòu)參數(shù)，使得SHAW模型可以獲得動態(tài)的植株結(jié)構(gòu)信息。WOFOST-SHAW模型在冬小麥拔節(jié)之后對冠層各高度層氣溫的整體模擬效果較好，表明在無實測生物量數(shù)據(jù)的情況下，WOFOST-SHAW模型可以合理模擬田間冠層氣溫整體動態(tài)變化。傳統(tǒng)人工采樣方式獲取的生物量將某一次采樣值代表某一時段，周期長，破壞性大、且費時費力。而WOFOST模型能較合理地模擬每日冬小麥的生物量變化，它提供動態(tài)植株結(jié)構(gòu)參數(shù)后，使SHAW模型適應(yīng)性增強。同時，WOFOST模型本身的也具有局限性。在水分限制生長模式下，它僅考慮水分而假設(shè)養(yǎng)分、病蟲害和耕作措施等其他條件均處于最優(yōu)狀況，致使開花期之后WOFOST模型模擬生物量的偏差大于拔節(jié)期和抽穗期。這可能導(dǎo)致WOFOST-SHAW模型在冬小麥生長后期氣溫模擬值偏差變大，出現(xiàn)整體看冠層上部與冠層下部各高度層誤差變化趨勢不一致的現(xiàn)象。雖然WOFOST-SHAW模型在冠層各高度模擬效果存在一定差異，但整體看WOFOST-SHAW模型的模擬結(jié)果與實測值擬合度高，解決了SHAW模型參數(shù)獲取困難的問題。

WOFOST-SHAW模型可以較好地模擬冠層日最低氣溫，可以為冬小麥晚霜凍害監(jiān)測提供比較準(zhǔn)確的監(jiān)測指標(biāo)。同時隨冠層高度增加模型模擬精度下降，這與前人研究結(jié)果一致[5]。從極端低溫的模擬情況來看，WOFOST-SHAW模型可以模擬冬小麥低溫敏感生育期內(nèi)易受凍高度層的低于0℃的低溫環(huán)境，并合理模擬零下日最低氣溫，同時模型存在低溫模擬值偏高的問題，模擬精度有待進一步提高。

4 結(jié)論

(1)提出了一種通過WOFOST模型直接或間接獲取SHAW模型運行所需的生物指標(biāo)的方法，該方法不僅避免了傳統(tǒng)采樣方式費時費力的不足，有效地解決了SHAW模型植株結(jié)構(gòu)參數(shù)獲取困難的問題，而且可以獲取連續(xù)變化的生物量參數(shù)，提高了SHAW模型的適用性。

(2)WOFOST-SHAW模型對冬小麥冠層各高度氣溫模擬值與實測值的R2均不小于0.93，ME基本上均大于0.90，表明模型能較好地模擬冠層氣溫，整體模擬精度比較高；除80 cm高度外，各高度層的RMSE均在1.41～2.11℃之間，超過75%的氣溫模擬誤差在-1.5～2.5℃范圍內(nèi)，表明冠層各高度層模擬誤差在合理誤差范圍內(nèi)。

(3)WOFOST-SHAW模型對冠層日最低氣溫的模擬時，ME平均值為0.86，RMSE與MBE的平均值分別為1.44℃、0.53℃，說明模型對冠層日最低氣溫的擬合度較高。同時，模型在冬小麥易受凍高度層比較準(zhǔn)確地反映低于0℃的低溫環(huán)境和日最低氣溫。

1 劉云, 宇振榮, 孫丹峰. 冬小麥冠氣溫差及其相關(guān)影響因素關(guān)系研究[J]. 灌溉排水學(xué)報, 2004, 23(1):30-35. LIU Yun, YU Zhenrong, SUN Danfeng. Difference of canopy-air temperature and its affecting factors[J]. Journal of Irrigation and Drainage, 2004, 23(1):30-35. (in Chinese)

2 楊邦杰, 王茂新, 裴志遠(yuǎn). 冬小麥凍害遙感監(jiān)測[J]. 農(nóng)業(yè)工程學(xué)報, 2002, 18(2): 136-140. YANG Bangjie, WANG Maoxin, PEI Zhiyuan. Monitoring freeze injury to winter wheat using remote sensing[J]. Transactions of the CSAE, 2002, 18(2): 136-140. (in Chinese)

3 王慧芳, 王紀(jì)華, 董瑩瑩, 等. 冬小麥凍害脅迫高光譜分析與凍害嚴(yán)重度反演[J]. 光譜學(xué)與光譜分析, 2014, 34(5): 1357-1361. WANG Huifang, WANG Jihua, DONG Yingying, et al. Monitoring freeze stress levels on winter wheat from hyperspectral reflectance data using principal component analysis[J]. Spectroscopy and Spectral Analysis, 2014, 34(5): 1357-1361. (in Chinese)

4 FLERCHINGER G N, PIERSON F B. Modelling plant canopy effects on variability of soil temperature and water: model calibration and validation[J]. Journal of Arid Environments, 1997, 35(4): 641-653.

5 LI R, SHI H, FLERCHINGER G N, et al. Modeling the effect of antecedent soil water storage on water and heat status in seasonally freezing and thawing agricultural soils[J]. Geoderma, 2013, 206(9): 70-74.

6 ZHANG Y, CHENG G, XIN L, et al. Coupling of a simultaneous heat and water model with a distributed hydrological model and evaluation of the combined model in a cold region watershed[J]. Hydrological Processes, 2013, 27(25): 3762-3776.

7 肖薇, 鄭有飛, 于強. 基于SHAW模型對農(nóng)田小氣候要素的模擬[J]. 生態(tài)學(xué)報, 2005, 25(7): 1626-1634. XIAO Wei, ZHENG Youfei, YU Qiang. Evaluation of SHAW model in simulating energy balance, leaf temperature and micrometeorological variables within a maize canopy[J]. Acta Ecologica Sinica, 2005,25(7):1626-1634. (in Chinese)

8 劉峻明, 汪念, 王鵬新, 等. 基于SHAW模型的冬小麥近地面層氣溫模擬[J/OL]. 農(nóng)業(yè)機械學(xué)報, 2015, 46(增刊): 274-282.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=2015s044&flag=1. DOI：10.6041/j.issn. 1000-1298.2015.S0.044. LIU Junming, WANG Nian, WANG Pengxin, et al. Simulation of air temperature within winter wheat near-ground layer based on SHAW model[J/OL].Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(Supp.): 274-282.(in Chinese)

9 劉峻明, 汪念, 王鵬新, 等. SHAW模型在冬小麥晚霜凍害監(jiān)測中的適用性研究[J/OL]. 農(nóng)業(yè)機械學(xué)報, 2016, 47(6): 265-274.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20160635&flag=1. DOI：10.6041/j.issn. 1000-1298.2016.06.035. LIU Junming,WANG Nian, WANG Pengxin, et al. Applicability of simultaneous heat and water model for monitoring late frost injury of winter wheat[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(6): 265-274. (in Chinese)

10 何亮, 侯英雨, 趙剛, 等. 基于全局敏感性分析和貝葉斯方法的WOFOST作物模型參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報, 2016, 32(2): 169-179. HE Liang, HOU Yingyu, ZHAO Gang,et al. Parameters optimization of WOFOST model by integration of global sensitivity analysis and Bayesian calibration method[J]. Transactions of the CSAE, 2016, 32(2): 169-179.(in Chinese)

11 劉峻明, 李曼曼, 王鵬新, 等. 基于LAI時間序列重構(gòu)數(shù)據(jù)的冬小麥物候監(jiān)測[J]. 農(nóng)業(yè)工程學(xué)報, 2013, 29(19): 120-129. LIU Junming, LI Manman, WANG Pengxin, et al. Monitoring of phenology by reconstructing LAI time series data for winter wheat[J]. Transactions of the CSAE, 2013, 29(19): 120-129.(in Chinese)

12 黃健熙, 李昕璐, 劉帝佑, 等. 順序同化不同時空分辨率LAI的冬小麥估產(chǎn)對比研究[J/OL]. 農(nóng)業(yè)機械學(xué)報, 2015,46(1): 240-248.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20150134&flag=1.DOI：10.6041/j.issn.1000-1298.2015.01.034. HUANG Jianxi, LI Xinlu, LIU Diyou, et al. Comparison of winter wheat yield estimation by sequential assimilation of different spatio-temporal resolution remotely sensed LAI datasets[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(1): 240-248.(in Chinese)

13 CHENG Z, MENG J, WANG Y. Improving spring maize yield estimation at field scale by assimilating time-series HJ-1 CCD data into the WOFOST model using a new method with fast algorithms[J]. Remote Sensing, 2016, 8(4): 303.

14 張建平, 趙艷霞, 王春乙, 等. 基于WOFOST作物生長模型的冬小麥干旱影響評估技術(shù)[J]. 生態(tài)學(xué)報, 2013, 33(6):1762-1769. ZHANG Jianping, ZHAO Yanxia, WANG Chunyi, et al. Evaluation technology on drought disaster to yields of winter wheat based on WOFOST crop growth model[J]. Acta Ecologica Sinica, 2013, 33(6) : 1762-1769. (in Chinese)

15 鄔定榮, 歐陽竹, 趙小敏, 等. 作物生長模型WOFOST在華北平原的適用性研究[J]. 植物生態(tài)學(xué)報, 2003, 27(5): 594-602. WU Dingrong, OUYANG Zhu, ZHAO Xiaomin, et al. The applicability research of WOFOST model in North China Plain[J]. Chinese Journal of Plant Ecology, 2003, 27(5): 594-602.(in Chinese)

16 張雪芬, 余衛(wèi)東, 王春乙. 基于作物模型災(zāi)損識別的黃淮區(qū)域冬小麥晚霜凍風(fēng)險評估[J]. 高原氣象, 2012, 31(1): 277-284. ZHANG Xuefen, YU Weidong, WANG Chunyi.Risk evaluation for spring frost disaster of winter wheat in Yellow River-Huai River regions based on crop model[J]. Plateau Meteorology, 2012, 31(1): 277-284.(in Chinese)

17 FLERCHINGER G N, HANSO C L, WIGHT J R. Modeling evapotranspiration and surface energy budgets across a watershed[J]. Water Resources Research, 1996, 32(8): 2539-2548.

18 FANG Q X, MA L, FLERCHINGER G N, et al. Modeling evapotranspiration and energy balance in a wheat-maize cropping system using the revised RZ-SHAW model[J]. Agricultural and Forest Meteorology, 2014, 194(3): 218-229.

19 YU Q, FLERCHINGER G N, XU S, et al. Energy balance simulation of a wheat canopy using the RZ-SHAW (RZWQM-SHAW) model[J]. Transactions of the ASABE, 2007, 50(5): 1507-1516.

20 FLERCHINGER G N, REBA M L, LINK T E, et al. Modeling temperature and humidity profiles within forest canopies[J]. Agricultural and Forest Meteorology, 2015, 213: 251-262.

21 QI Z, MA L, BAUSCH W C, et al. Simulating maize production, water and surface energy balance, canopy temperature, and water stress under full and deficit irrigation[J]. Transactions of the ASABE, 2016, 59(2): 623-633.

22 鐘秀麗, 王道龍, 趙鵬，等.黃淮麥區(qū)冬小麥拔節(jié)后霜凍溫度出現(xiàn)規(guī)律研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報， 2007, 15(5): 17-20. ZHONG Xiuli, WANG Daolong, ZHAO Peng, et al. Occurrence of frost temperature in Huanghuai wheat production zone after winter wheat elongation[J]. Chinese Journal of Eco-Agriculture, 2007, 15(5): 17-20. (in Chinese)