摘要：發(fā)展針對異常氣候條件下的作物估產(chǎn)模型，對于理解氣候變化對作物產(chǎn)量的影響，提高估產(chǎn)模型的適用性具有重要意義。本文提出了一種綜合氣象災(zāi)害指標(biāo)、遙感植被指數(shù)與趨勢產(chǎn)量的隨機森林回歸估產(chǎn)方法，基于該方法構(gòu)建了中國五大麥區(qū)的小麥單產(chǎn)估算模型，并選擇典型的災(zāi)害年份在縣級尺度和麥區(qū)尺度上分別進(jìn)行了精度驗證，對比分析了不同麥區(qū)輸入變量對估產(chǎn)模型構(gòu)建的重要性。結(jié)果顯示：五大麥區(qū)估產(chǎn)模型擬合精度的R2均在0.95以上，各麥區(qū)縣級實際單產(chǎn)與預(yù)測單產(chǎn)平均相對誤差均低于0.060，區(qū)級均低于0.049;輸入變量在不同麥區(qū)的重要性存在差異。標(biāo)準(zhǔn)化降水蒸散指數(shù)在五大麥區(qū)重要性均較高，干熱風(fēng)指標(biāo)在西北春麥區(qū)與北部春麥區(qū)相比其他區(qū)域更重要，月平均溫度距平與月平均降水距平在各個麥區(qū)的重要性差異不大，小麥拔節(jié)期與抽穗期的標(biāo)準(zhǔn)化差分植被指數(shù)(NDVI)重要性相對其他時間段的NDVI更高。本文構(gòu)建的估產(chǎn)模型能夠滿足異常氣候條件下估產(chǎn)的精度，可以為極端氣候條件下大尺度的估產(chǎn)研究提供參考。

　　關(guān)鍵詞：隨機森林;估產(chǎn);標(biāo)準(zhǔn)化降水蒸散指數(shù);標(biāo)準(zhǔn)化差值植被指數(shù);干熱風(fēng)指標(biāo);溫度距平;降水距平

　　從20世紀(jì)上半葉至今，地球大氣中的CO2濃度不斷升高，變暖趨勢越來越顯著，全球氣候變化，導(dǎo)致異常氣候事件頻發(fā)[1]。相關(guān)研究表明，盡管農(nóng)業(yè)生產(chǎn)受氣候變化的影響在不同區(qū)域有一定的差異，但整體上弊大于利[2]。我國農(nóng)業(yè)受氣候變化的影響十分嚴(yán)重[3]，在各類農(nóng)業(yè)氣象災(zāi)害中，干旱和洪澇是影響我國農(nóng)作物生產(chǎn)最主要的氣象災(zāi)害。1980年以來，我國洪旱災(zāi)害頻發(fā)。1998年的長江和松花江流域特大洪水事件和2016年的長江中下游及太湖洪澇事件直接帶來上千億元經(jīng)濟(jì)損失[4]。

　　相比洪澇，干旱對我國農(nóng)業(yè)生產(chǎn)的影響更為嚴(yán)重[5]。有關(guān)統(tǒng)計數(shù)據(jù)表明，2007至2017年全國農(nóng)作物旱災(zāi)面積年均值為1.53×107hm2，糧食損失年均值為2.14億t，造成直接經(jīng)濟(jì)損失年均值為282.15億元，占GDP均值的0.21%[6]。在全國各種氣象災(zāi)害中，干旱導(dǎo)致的農(nóng)作物受災(zāi)面積占比約56%，洪澇導(dǎo)致的農(nóng)作物受災(zāi)面積占比約24%[7]。此外，低溫冷害、干熱風(fēng)等氣候事件的頻繁發(fā)生也給我國農(nóng)作物安全生產(chǎn)帶來了巨大挑戰(zhàn)。在中國北方麥區(qū)，干熱風(fēng)災(zāi)害的發(fā)生一般會導(dǎo)致小麥減產(chǎn)5%~10%，在受災(zāi)嚴(yán)重的年份減產(chǎn)率能達(dá)到20%~30%[8]。

　　隨著遙感技術(shù)的不斷發(fā)展與創(chuàng)新應(yīng)用，許多專家學(xué)者將遙感技術(shù)應(yīng)用到農(nóng)作物估產(chǎn)中[9]，為作物估產(chǎn)研究開辟了新途徑[10]。從模型建立的理論角度 出發(fā)，當(dāng)前主流的作物估產(chǎn)模型可以分為四類：經(jīng)驗統(tǒng)計模型、光能利用率模型、作物生長模擬模型和耦合模型[11]。探究氣候變化對作物產(chǎn)量的影響，大多采用經(jīng)驗統(tǒng)計模型中的氣象估產(chǎn)模型。氣象估產(chǎn)模型發(fā)展較早，始于20世紀(jì)70年代末，基本的氣象估產(chǎn)模型有三種。第一種，直接建立氣象因子和作物產(chǎn)量之間的回歸模型;第二種，首先計算相鄰兩年作物產(chǎn)量差和氣象因子差，然后建立作物產(chǎn)量差和氣象因子差之間的回歸模型，進(jìn)而進(jìn)行作物產(chǎn)量預(yù)測。

　　第三種，將作物產(chǎn)量分解為趨勢產(chǎn)量和氣象產(chǎn)量兩部分，前者反映由技術(shù)進(jìn)步(如灌溉、施肥、新品種等)導(dǎo)致的產(chǎn)量的長期變化，后者反映由自然氣候要素(光照、降水、輻射等)引起的產(chǎn)量的短期波動。常用的氣象估產(chǎn)模型構(gòu)建方法包括線性和非線性回歸兩大類[12-13]。隨著數(shù)據(jù)挖掘技術(shù)的發(fā)展，神經(jīng)網(wǎng)絡(luò)、支持向量機、隨機森林等新的數(shù)據(jù)挖掘手段也應(yīng)用于估產(chǎn)模型的構(gòu)建[14-16]。

　　Breiman[17]在2001年提出了隨機森林算法。它是一種集成學(xué)習(xí)算法，具有穩(wěn)定性好，預(yù)測精度高，不易產(chǎn)生過擬合等優(yōu)點，不需要顧慮多重共線性的問題，與神經(jīng)網(wǎng)絡(luò)和線性回歸方法相比，其表現(xiàn)更穩(wěn)定，具有較強的抗干擾能力。該算法已經(jīng)廣泛應(yīng)用于干旱監(jiān)測[18]、山體滑坡預(yù)測[19]、洪澇[20]和火災(zāi)風(fēng)險[21]分析等方面。但目前應(yīng)用隨機森林算法進(jìn)行異常氣候條件下的作物估產(chǎn)的研究甚少。

　　目前大多研究也并沒有驗證所建立的估產(chǎn)模型在災(zāi)害年份的估產(chǎn)效果。為此，本文基于隨機森林算法，以氣象估產(chǎn)模型為基礎(chǔ)，結(jié)合氣象災(zāi)害指數(shù)和遙感植被指數(shù)，在中國五大小麥種植區(qū)分別建立估產(chǎn)模型，并驗證其在典型災(zāi)害年份的估產(chǎn)精度，為探究氣候變化對小麥產(chǎn)量的影響，以及建立滿足異常氣候下估產(chǎn)精度的模型提供參考。

　　1研究區(qū)概況

　　本文研究范圍包括西北春麥區(qū)、北部冬麥區(qū)、黃淮海冬麥區(qū)、長江中下游冬麥區(qū)與西南冬麥區(qū)[22-23]。西北春麥區(qū)位于黃淮上游三大高原的交匯地帶，該區(qū)由寧夏全區(qū)、甘肅省大部分地區(qū)，內(nèi)蒙古西北部地區(qū)構(gòu)成，冬天寒冷、夏季炎熱，春季與秋季多風(fēng)，天氣較為干燥，白天和晚上溫差很大。降水不足是影響該區(qū)小麥作物生長的重要因素，部分地區(qū)春小麥生長后期有干熱風(fēng)危害。北部冬麥區(qū)位于我國中緯度地帶，由北京全市、天津全市、河北、山西以及陜西省部分地區(qū)組成，屬于暖溫帶季風(fēng)區(qū)域。

　　該區(qū)冬寒春旱，降水不足且分布不均勻。黃淮海冬麥區(qū)位于黃河中下游，由山東省全部、河南省大部分區(qū)域以及陜西省與山西省南部、安徽省與江蘇省北部、河北南部的少部分地區(qū)組成。該區(qū)北部小麥在低溫年份有遭受低溫冷害的風(fēng)險，南部地區(qū)氣溫較高，冬小麥返青期不明顯。降水量南多北少，東多西少，小麥生育期降雨可基本滿足，但北部仍偶爾有旱災(zāi)發(fā)生，需灌溉。

　　長江中下游冬麥區(qū)由湖北全省、湖南全省、江西全省、浙江全省、上海全市，安徽省與浙江省大部分地區(qū)組成。該區(qū)大部分地區(qū)的自然條件適宜小麥生長，但該區(qū)降水極不平衡，南部降水過多不適宜小麥生長，北部降水較少時有干旱發(fā)生。西南冬麥區(qū)由重慶全市、貴州全省、四川省東部與云南省北部組成。該麥區(qū)地形復(fù)雜，包括山地、盆地、平原等多種地形，其中山地為該區(qū)主要地形，復(fù)雜的地形也導(dǎo)致該區(qū)域氣候環(huán)境差異較大，農(nóng)作物種植、成熟時間差異較大，除平原地區(qū)外，農(nóng)田地塊較為破碎。

　　2技術(shù)路線與研究方法

　　主要包括以下步驟：1)利用原始?xì)庀髷?shù)據(jù)與遙感數(shù)據(jù)建立相應(yīng)氣象災(zāi)害指標(biāo)與遙感植被指數(shù);2)利用直線滑動平均法和HP濾波法擬合趨勢產(chǎn)量，選擇最佳擬合方法;3)以氣象災(zāi)害指標(biāo)、遙感植被指數(shù)為輸入變量，實際產(chǎn)量為輸出變量，建立樣本庫;4)利用隨機森林算法在不同麥區(qū)建立估產(chǎn)模型并進(jìn)行精度驗證;5)篩選出各麥區(qū)不同類型災(zāi)害年份，并進(jìn)行估產(chǎn)驗證;6)對不同麥區(qū)輸入變量進(jìn)行重要性評價。

　　2.1數(shù)據(jù)獲取與預(yù)處理

　　通過查閱資料發(fā)現(xiàn)在我國眾多不同類型氣象災(zāi)害中，對農(nóng)作物生產(chǎn)影響較大的災(zāi)害類型主要包括旱災(zāi)、洪澇、低溫凍害以及干熱風(fēng)四種[24-28]。為此，本文收集了氣象數(shù)據(jù)與遙感數(shù)據(jù)，用于構(gòu)建表示干旱、洪澇、低溫冷害、干熱風(fēng)四種災(zāi)害的指標(biāo)，以及表示小麥不同時期長勢情況的遙感植被指數(shù)指標(biāo)，用于估產(chǎn)模型的建立。具體包括：1)來自全球標(biāo)準(zhǔn)化降水蒸散指數(shù)(StandardizedPrecipitationEvapotranspirationIndex,SPEI)數(shù)據(jù)庫1999—2018年的SPEIbasev2.6數(shù)據(jù)集，時間尺度為1~48月，分辨率為0.5°;2)來自歐空局哥白尼氣象數(shù)據(jù)中心的1999—2018年的月平均降水量、月平均溫度、14時風(fēng)速、14時相對濕度、日最高溫度，除了日最高溫度分辨率為0.1°外，其他數(shù)據(jù)的空間分辨率均為0.25°;3)來自美國航空航天局的2000—2018年的16d合成的MOD13Q1-NDVI數(shù)據(jù)，空間分辨率為250m。

　　2.2輸入?yún)��?shù)計算

　　五大麥區(qū)范圍廣，小麥種植時間差異較大，各 麥區(qū)的小麥生育期間分別為黃淮海冬麥區(qū)9月-次年6月、北部冬麥區(qū)9月-次年6月、西南冬麥區(qū)8月次年7月、長江中下游冬麥區(qū)10月-次年5月、西北春麥區(qū)3月-8月，因此本文以各麥區(qū)小麥的生育期并集為研究時段來對參數(shù)進(jìn)行計算。

　　2.2.1干旱指標(biāo)本文選用SPEI8-6、SPEI6-10、SPEI6-10、SPEI5-8、SPEI7-12分別表示西北春麥區(qū)、北部冬麥區(qū)、黃淮海冬麥區(qū)、長江中下游冬麥區(qū)、西南冬麥區(qū)中小麥整個生育期內(nèi)的干旱程度。其中，SPEIm-n為m月份時間尺度為n的SPEI干旱指數(shù)，用來描述m月份過去n個月整體的干旱情況。具體計算流程是首先將SPEI干旱指數(shù)重采樣至1km空間尺度，用全國耕地范圍數(shù)據(jù)進(jìn)行耕地掩膜，然后以縣域范圍內(nèi)SPEI干旱指數(shù)均值作為描述某縣小麥生育期內(nèi)整體干旱程度的指標(biāo)。

　　2.2.2干熱風(fēng)指標(biāo)我國小麥干熱風(fēng)具體可以分為高溫低濕型、雨后青枯型和旱風(fēng)型3種[29]。國家氣象局發(fā)布了2019年干熱風(fēng)國家標(biāo)準(zhǔn)，采用日最高氣溫、14時空氣相對濕度和14時風(fēng)速組合，結(jié)合20cm土壤相對濕度確定了3種類型小麥干熱風(fēng)指標(biāo)。

　　不同地區(qū)干熱風(fēng)類型不同，評定標(biāo)準(zhǔn)也不同，同時干熱風(fēng)災(zāi)害主要發(fā)生在我國北方麥區(qū)，具有局地性特征。為了統(tǒng)一各麥區(qū)干熱風(fēng)指標(biāo)，本文采用了“三三三”標(biāo)準(zhǔn)對各麥區(qū)小麥生育期并集內(nèi)干熱風(fēng)有效天數(shù)進(jìn)行統(tǒng)計，也即當(dāng)小麥生育期并集內(nèi)某天同時滿足日最高氣溫>30℃、14時空氣相對濕度<3014>3m/s三個條件時即為一個干熱風(fēng)有效天。干熱風(fēng)有效天數(shù)指小麥生育期并集內(nèi)干熱風(fēng)有效天的累積。

　　2.2.3洪澇與低溫冷害指標(biāo)洪澇指標(biāo)的計算流程是：首先在像元尺度上計算小麥生育期內(nèi)某月份平均降水的歷年均值，然后用當(dāng)月平均降水實際值與歷年均值作差得到當(dāng)月平均降水距平值，最后取縣域內(nèi)耕地像元距平值的均值作為描述該縣當(dāng)月洪澇情況的指標(biāo)(下文簡稱為月平均降水距平)。低溫冷害指標(biāo)的計算方法與洪澇指標(biāo)類似，計算過程不再贅述，下文簡稱為月平均溫度距平。

　　2.2.4遙感植被指數(shù)計算首先對MOD13Q1-NDVI數(shù)據(jù)利用均值濾波從250m重采樣至1km，然后用全國土地利用分布圖對2000—2018年小麥生育期內(nèi)MOD13Q1-NDVI數(shù)據(jù)進(jìn)行耕地掩膜，最后用全國縣級矢量對各縣耕地范圍內(nèi)NDVI取均值，作為各縣的遙感植被指數(shù)。

　　2.3趨勢產(chǎn)量擬合

　　在氣象估產(chǎn)模型中，通常將小麥實際產(chǎn)量分解為趨勢產(chǎn)量、氣象產(chǎn)量和隨機產(chǎn)量三部分，其中隨機產(chǎn)量往往忽略不計。因此，實際單產(chǎn)可以分解為趨勢單產(chǎn)和氣象單產(chǎn)之和。 本文利用直線滑動平均法(包括3a和5a直線滑動平均法)和HP濾波法[30]對五大麥區(qū)所有縣小麥實際單產(chǎn)進(jìn)行趨勢產(chǎn)量擬合，用實際單產(chǎn)減去趨勢單產(chǎn)得到氣象單產(chǎn)。

　　以氣象單產(chǎn)為輸出變量，以四類氣象災(zāi)害指標(biāo)為輸入變量，分別在五大麥區(qū)構(gòu)建隨機森林回歸模型。模型擬合精度越高說明對氣象產(chǎn)量擬合效果越好，也即對趨勢產(chǎn)量擬合效果越好，以此來篩選最適合的趨勢產(chǎn)量擬合方法。模型擬合精度用驗證樣本的決定系數(shù)R2，均方根誤差(RMSE)和平均絕對誤差(MAE)進(jìn)行評價。

　　2.4隨機森林回歸與驗證

　　隨機森林算法的主要思想是從原始訓(xùn)練集中有放回地抽取k個樣本，且每個樣本的樣本容量均與原始訓(xùn)練集的大小一致;然后對每個樣本分別進(jìn)行回歸樹建模，得到k個建模結(jié)果，最后以每一棵回歸樹預(yù)測結(jié)果的平均值作為最終預(yù)測結(jié)果[31]。

　　隨機森林估產(chǎn)模型的輸入變量包括趨勢單產(chǎn)1個、干熱風(fēng)有效天數(shù)1個、SPEI干旱指數(shù)1個、生育期內(nèi)月平均降水距平N個、生育期內(nèi)月平均溫度距平N個、生育期內(nèi)每16d的NDVI植被指數(shù)M個，合計2N+M+3個。需要注意的是，不同麥區(qū)小麥的生育期不同，因此生育期內(nèi)月平均溫度距平和月平均降水距平的個數(shù)N在不同麥區(qū)將發(fā)生改變。以樣本集中3/4樣本子集作為訓(xùn)練樣本建立回歸模型，1/4樣本子集作為驗證樣本進(jìn)行估產(chǎn)模型擬合精度的驗證。隨機森林回歸模型擬合精度用模型決定系數(shù)R2、驗證樣本的均方根誤差、平均絕對誤差以及平均相對誤差進(jìn)行評價。

　　2.5災(zāi)害年份估產(chǎn)精度驗證

　　采用麥區(qū)內(nèi)主要省(直轄市)農(nóng)作物受災(zāi)面積累計的大小來表示各麥區(qū)整體受干旱、洪澇、低溫冷害三類氣象災(zāi)害導(dǎo)致的農(nóng)作物受災(zāi)情況。

　　另外，由于年鑒統(tǒng)計數(shù)據(jù)中不包含干熱風(fēng)災(zāi)害統(tǒng)計數(shù)據(jù)，因此對于不同麥區(qū)干熱風(fēng)災(zāi)害年份的確定，主要參考了其他專家學(xué)者的相關(guān)研究。災(zāi)害年份確定之后，將五大麥區(qū)災(zāi)害年份對應(yīng)的樣本進(jìn)行剔除，對各麥區(qū)重新建立隨機森林回歸模型，然后將各麥區(qū)旱災(zāi)、洪澇災(zāi)害、低溫冷害、干熱風(fēng)災(zāi)害年份對應(yīng)的樣本的自變量代入重新訓(xùn)練好的模型進(jìn)行災(zāi)害年份產(chǎn)量預(yù)測。

　　最后在縣級尺度和麥區(qū)尺度分別進(jìn)行了精度驗證�？h級尺度驗證用各麥區(qū)中所有縣小麥單產(chǎn)真實值和預(yù)測值的相對誤差的平均值表示。麥區(qū)尺度驗證用各麥區(qū)實際單產(chǎn)和預(yù)測單產(chǎn)的相對誤差進(jìn)行評價。其中，各麥區(qū)實際單產(chǎn)是先用各縣實際單產(chǎn)與各縣實際種植面積計算得到麥區(qū)實際總產(chǎn)量，再用實際總產(chǎn)量除以總種植面積得到。同理，各麥區(qū)預(yù)測單產(chǎn)是先用各縣預(yù)測單產(chǎn)與各縣實際種植面積計算得到麥區(qū)總預(yù)測產(chǎn)量，再用麥區(qū)總預(yù)測產(chǎn)量除以總種植面積得到。

　　2.6特征重要性評價方法

　　用隨機森林中某個樹節(jié)點的樣本方差與分裂后兩個葉子節(jié)點樣本方差進(jìn)行差值運算，用該差值來評價某個特征的重要性程度，差值降低越高表示該特征重要性越大。計算得到模型所有特征重要性后，對所有特征重要性值作歸一化處理[32]。主要說明的是：對于某一模型，所有特征重要性值累積為1，特征重要性具體的值本身并無意義，因此分析時關(guān)注的是特征重要性值之間的相對差異。

　　3研究結(jié)果

　　以2000—2019年《中國農(nóng)村統(tǒng)計年鑒》全國各省農(nóng)作物自然氣候災(zāi)害損失數(shù)據(jù)為基礎(chǔ)，以各麥區(qū)主要省(直轄市)農(nóng)作物受災(zāi)面積累計的大小為依據(jù)，對五大麥區(qū)不同類型災(zāi)害的典型年份進(jìn)行確定。例如，西北春麥區(qū)主要包括甘肅、寧夏以及內(nèi)蒙古西北部地區(qū)，內(nèi)蒙古西北部地區(qū)多為戈壁沙漠地帶，故用寧夏與甘肅農(nóng)作物受災(zāi)面積的累計來表示該區(qū)整體農(nóng)作物受災(zāi)狀況;西南冬麥區(qū)主要包括重慶與貴州全部、四川中東部、云南中北部以 及陜西南部少部分地區(qū)，故用重慶與貴州農(nóng)作物受災(zāi)面積累計并參考四川、云南農(nóng)作物受災(zāi)面積來表示該區(qū)整體農(nóng)作物受災(zāi)狀況，其他三個麥區(qū)進(jìn)行類似分析來確定不同類型災(zāi)害的典型年份。

　　低溫冷害、干旱和洪澇災(zāi)害嚴(yán)重年份在西北春麥區(qū)分別為2004年、2000年和2013年，在北部冬麥區(qū)分別為2006年、2001年和2007年，在黃淮海冬麥區(qū)分別為2005年、2001年和2003年，在長江中下游冬麥區(qū)分別為2008年、2001年和2003年，在西南冬麥區(qū)分別為2008年、2001年和2002年。

　　干熱風(fēng)災(zāi)害有較強的區(qū)域性，其主要集中在我國河西地區(qū)、黃淮海平原等地區(qū)[33]，即西北春麥區(qū)、北方冬麥區(qū)以及黃淮海冬麥區(qū)。已有研究表明[26,29,34-35]，1960年至2017年河西地區(qū)干熱風(fēng)發(fā)生總次數(shù)和日數(shù)先緩慢減少后迅速增多，其中2004年的干熱風(fēng)發(fā)生總天數(shù)最高，持續(xù)天數(shù)也達(dá)到了歷年最長(12d);1961年至2015年黃淮海平原地區(qū)干熱風(fēng)日數(shù)總體呈下降趨勢，但是在2001年黃淮海平原地區(qū)仍發(fā)生了一次嚴(yán)重的干熱風(fēng)災(zāi)害，區(qū)域平均干熱風(fēng)日數(shù)達(dá)8.1d。因此，本文將　2004年作為西北春麥區(qū)干熱風(fēng)災(zāi)害年份，2001年作為北部冬麥區(qū)與黃淮海冬麥區(qū)的干熱風(fēng)災(zāi)害年份。長江中下游冬麥區(qū)與西南冬麥區(qū)未進(jìn)行干熱風(fēng)災(zāi)害年份的確定。

　　4討論

　　4.1本研究的主要科學(xué)貢獻(xiàn)

　　1)傳統(tǒng)的氣象估產(chǎn)模型的輸入變量往往是原始的氣象要素，如降水量、溫度等，而非氣候指標(biāo)。氣候指標(biāo)是一定氣候條件下的單項氣候要素或多種氣候要素綜合的特征量，如干旱指數(shù)，熱浪指數(shù)，高溫指數(shù)等。氣候指標(biāo)能更加綜合的反映氣候條件的變化，特別是考慮作物生長條件的農(nóng)業(yè)氣候指標(biāo)，能更好的衡量農(nóng)業(yè)氣象條件的利弊。本文從影響小麥生長的主要農(nóng)業(yè)自然災(zāi)害(干旱、洪澇、干熱風(fēng)和低溫凍害)出發(fā)，計算了SPEI干旱指數(shù)、降水距平、干熱風(fēng)指數(shù)和溫度距平來反映研究區(qū)中四個主要農(nóng)業(yè)自然災(zāi)害的情況，據(jù)此構(gòu)建的估產(chǎn)模型經(jīng)過驗證可以有效的進(jìn)行災(zāi)害年份小麥的產(chǎn)量估計。

　　2)NDVI指數(shù)(包括NDVI的最大值、生長季的NDVI累計值等)是遙感估產(chǎn)中最常用的植被指數(shù)。傳統(tǒng)的遙感統(tǒng)計估產(chǎn)模型往往只依托于遙感指數(shù)，而不考慮非遙感數(shù)據(jù)。本研究綜合使用遙感數(shù)據(jù)和氣象數(shù)據(jù)進(jìn)行估產(chǎn)模型的建立。相比單純的遙感數(shù)據(jù)，加入氣候數(shù)據(jù)能夠增加模型估產(chǎn)結(jié)果的穩(wěn)定性，減少由傳感器衰變，大氣條件差異等造成的植被指數(shù)本身年際間的變化而產(chǎn)生的估產(chǎn)誤差。相比單純使用氣候指數(shù)，遙感數(shù)據(jù)能夠在一定程度上反映人類活動的影響，例如人工灌溉對干旱的抵御作用，從而避免僅僅考慮氣候災(zāi)害指數(shù)而帶來的產(chǎn)量低估。

　　3)本文利用隨機森林算法對各個麥區(qū)估產(chǎn)模型的輸入變量的重要性進(jìn)行了分析，可以幫助我們更好的理解影響各個麥區(qū)小麥產(chǎn)量的主要因素，為不同區(qū)域的小麥防災(zāi)減災(zāi)措施的制定提供參考，也可以為其他區(qū)域的小麥估產(chǎn)研究和其他作物估產(chǎn)模型的建立提供輸入變量篩選的方法借鑒。

　　4.2本研究的不足之處

　　1)計算2000—2018年小麥種植區(qū)各縣四類災(zāi)害指標(biāo)與遙感植被指數(shù)時，考慮到工作量和試驗可行性，以及每年對全國小麥空間分布提取也同樣存在一定的誤差，最終選擇用2010年全國土地利用空間分類圖中耕地范圍來代替各縣小麥種植區(qū)進(jìn)行參數(shù)的計算，這會在一定程度上帶來誤差。例如對于NDVI變量，長江中下游冬麥區(qū)的5月25日至6月9日的NDVI重要性顯示很高，其次是4月7日至4月22日的NDVI，前者小麥已經(jīng)進(jìn)入收獲期，后者在抽穗期，這一結(jié)果和理論預(yù)期不符合。

　　究其原因可能是長江中下游的小麥種植面積總體小于水稻的種植面積，小麥?zhǔn)斋@期時以耕地范圍統(tǒng)計出來的NDVI指數(shù)很大程度上是水稻的信息，5月25日至6月9日大致在早稻的拔節(jié)孕穗期和中稻的移栽期，此時發(fā)生干旱對水稻的產(chǎn)量影響很大。

　　2)本文根據(jù)各省(直轄市)農(nóng)作物自然災(zāi)害損失數(shù)據(jù)進(jìn)行各麥區(qū)災(zāi)害年份的確定，一方面災(zāi)損數(shù)據(jù)并不只針對小麥，另一方面災(zāi)損數(shù)據(jù)是按照行政單元統(tǒng)計的，而麥區(qū)是跨行政單元的。我們以麥區(qū)包括的主要省(直轄市)的災(zāi)損總量來反映整個麥區(qū)的災(zāi)損情況，由此確定各麥區(qū)不同災(zāi)害類型的典型災(zāi)害年份，可能會產(chǎn)生一定誤差。在后續(xù)的研究中可以考慮根據(jù)時間序列氣象數(shù)據(jù)進(jìn)行不同類型災(zāi)害年份的確定。

　　3)本文選擇的SPEI干旱指數(shù)和干熱風(fēng)有效天數(shù)是對小麥生育期整體災(zāi)害情況的描述，然而發(fā)生在不同生育期的相同災(zāi)害對作物產(chǎn)量的影響是不同的，在后續(xù)的研究中需要在小麥不同生育期對災(zāi)害指標(biāo)進(jìn)行計算，以便分析在小麥不同生長階段發(fā)生的氣象災(zāi)害對產(chǎn)量的影響差異。

　　氣象工程論文： 氣象因子對赤霞珠果實GLVs香氣的影響

　　5結(jié)論

　　本文基于隨機森林算法，以氣象災(zāi)害指標(biāo)、遙感植被指數(shù)、趨勢產(chǎn)量為輸入變量，構(gòu)建了五大麥區(qū)的估產(chǎn)模型，并在典型災(zāi)害年份進(jìn)行了估產(chǎn)精度的驗證，得到如下主要結(jié)論：

　　1)精度驗證結(jié)果顯示，各麥區(qū)回歸模型的決定系數(shù)R2均達(dá)到了0.95以上，驗證樣本的平均相對誤差均低于0.073，均方根誤差均低于200.835kg/hm2，平均絕對誤差均低于131.655kg/hm2。各麥區(qū)所有縣災(zāi)害年份實際單產(chǎn)與預(yù)測單產(chǎn)平均相對誤差均低于0.060，各麥區(qū)區(qū)級實際單產(chǎn)與預(yù)測單產(chǎn)的相對誤差均低于0.049。

　　2)輸入特征的重要性在不同麥區(qū)具有一定的差異性，與地區(qū)氣候差異情況有關(guān)。SPEI干旱指數(shù)在五大麥區(qū)重要性均較高，表明干旱仍然是五大麥區(qū)影響小麥生產(chǎn)的重要因素;干熱風(fēng)指標(biāo)重要性在西北春麥區(qū)與北部春麥區(qū)較高，符合干熱風(fēng)主要發(fā)生在河西與黃淮海地區(qū)的局地性特點;各麥區(qū)月平均溫度距平重要性與月平均降水距平重要性整體差異比較小;總體來說拔節(jié)期與抽穗期的小麥NDVI相比其他時間段的NDVI更重要。

　　參考文獻(xiàn)：

　　[1]Masson-DelmotteV,ZhaiP,PiraniA,etal.ClimateChange2021:ThePhysicalScienceBasis,theWorkingGroupIContributiontotheSixthAssessmentReportoftheIntergovernmentalPanelonClimateChange[M].CambridgeUniversityPress:CambridgeandNewYork,2021.

　　[2]姜彤,李修倉,巢清塵,等.《氣候變化2014:影響、適應(yīng)和脆弱性》的主要結(jié)論和新認(rèn)知[J].氣候變化研究進(jìn)展,2014,10(3):157-166.JiangT,LiXC,ChaoQC,etal.Highlightsandunderstandingofclimatechange2014:Impacts,adaptation,andvulnerability[J].ProgressusInquisitionesDEMutationeClimatis,2014,10(3):157-166.

　　[3]權(quán)暢,景元書,譚凱炎.氣候變化對三大糧食作物產(chǎn)量影響研究進(jìn)展[J].中國農(nóng)學(xué)通報,2013,29(32):361-367.QuanC,JingYS,TanKY.Progressinresearchofimpactsofclimatechangeonyieldofthreemajorgraincrops[J].ChineseAgriculturalScienceBulletin,2013,29(32):361-367.

　　[4]賀大興.極端氣候?qū)χ袊�糧食產(chǎn)量影響的定量分析[J].中國農(nóng)業(yè)資源與區(qū)劃,2017,38(4):28-34.HeDX.QuantitativeanalysisontheeffectsofextremeclimateongainproductioninChina[J].ChineseJournalofAgriculturalResourcesandRegionalPlanning,2017,38(4):28-34.

　　作者：朱秀芳1,2,4，李石波3，石振君3，任笑3，程昌秀4