基于產(chǎn)量歷史豐歉氣象影響指數(shù)的河北省馬鈴薯產(chǎn)量預(yù)報

薛思嘉，魏瑞江，王朋朋，楊 梅，甄曉菊

（1.河北省氣象與生態(tài)環(huán)境重點實驗室，河北 石家莊050021；2.承德市氣象局，河北 承德067000；3.河北省氣象科學研究所，河北 石家莊050021）

馬鈴薯是繼水稻、小麥和玉米之后世界第四大糧食作物，也是河北省主要種植作物之一，其種植面積由20世紀80年代初的1048×105hm2增長到2018年的1631×105hm2，雖然面積增長僅56%，但其單產(chǎn)迅速增長，由1158.3 kg/hm2增長到6503.9 kg/hm2，增長了4.6倍。馬鈴薯生長發(fā)育和產(chǎn)量形成與氣象條件密切相關(guān)，產(chǎn)量年際之間波動較大，準確預(yù)報馬鈴薯的產(chǎn)量對河北省乃至全國糧食安全具有重要意義。

農(nóng)作物產(chǎn)量預(yù)報方法有農(nóng)學方法[1-2]、統(tǒng)計方法[3-5]、遙感方法[6-7]和作物模型方法[8-9]，目前在產(chǎn)量預(yù)報業(yè)務(wù)中應(yīng)用比較廣泛的統(tǒng)計方法中有基于作物氣候適宜度的方法[10-13]、基于作物產(chǎn)量歷史豐歉氣象影響指數(shù)方法[14-15]和基于關(guān)鍵氣象因子方法[16-17]?；谧魑锂a(chǎn)量歷史豐歉氣象影響指數(shù)方法，其原理是利用預(yù)報年作物生育階段氣象條件與歷史上任意一年的相似性，利用大概率法或加權(quán)平均法或單因子分析法或多因子分析法[18]等方法，確定預(yù)報年作物產(chǎn)量豐歉氣象影響指數(shù)來預(yù)報作物的產(chǎn)量[19]。該方法已應(yīng)用于黑龍江、江西、湖南、山東等省的水稻、油菜、早稻、小麥等作物的產(chǎn)量預(yù)報中[15，20-22]，得到了較好的效果。但由于該方法中大概率法、加權(quán)平均法、單因子分析法和多因子分析法等不同方法之間有一定的差異性和適宜性，哪一種方法更適合于河北省馬鈴薯產(chǎn)量預(yù)報，目前尚未見報道。

本文采用產(chǎn)量歷史豐歉氣象影響指數(shù)方法中常用的大概率法和加權(quán)平均法，分別預(yù)報河北省馬鈴薯的單產(chǎn)豐歉氣象影響指數(shù)，并分析比較兩種預(yù)報方法的準確率高低，對于做好河北省馬鈴薯產(chǎn)量預(yù)報業(yè)務(wù)，服務(wù)當?shù)剞r(nóng)業(yè)生產(chǎn)具有重要意義。

1 資料與方法

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

文中所用氣象資料來自河北省氣象局，包括河北省馬鈴薯主產(chǎn)區(qū)24個氣象站1982—2018年逐日最高氣溫、最低氣溫、降水量和日照時數(shù)；馬鈴薯發(fā)育期資料來自于河北省馬鈴薯主產(chǎn)區(qū)的農(nóng)業(yè)氣象觀測站；馬鈴薯產(chǎn)量資料來自于國家統(tǒng)計局，包括1982—2018年河北省馬鈴薯總產(chǎn)和種植面積。

1.2 資料處理

1.2.1 氣象資料處理

河北省馬鈴薯種植品種更替緩慢，據(jù)2019年統(tǒng)計，河北省馬鈴薯產(chǎn)區(qū)主要種植品種有“冀張薯”系列、“荷蘭15”、“費烏瑞它”、“克新一號”等，畝產(chǎn)分別為1700～2000、1500～2000、2000～2500、1600～2100 kg，種植面積占比分別為52.3%、22.9%、12.0%、1.1%。根據(jù)實際觀測表明，河北省馬鈴薯全生育期在4—9月，其中發(fā)芽期通常在4月下旬—5月中旬，幼苗期通常在5月中旬—6月中旬，塊莖形成期通常在6月中旬—7月中旬，塊莖增長期通常在7月中旬—8月上旬，淀粉積累期通常在8月上旬—9月中旬，成熟期通常在9月中旬—9月下旬。

基于河北省主產(chǎn)區(qū)24個氣象站歷史年與預(yù)報年作物播種后的日平均氣溫、降水量和日照時數(shù)等要素，按照參考文獻[23]的資料處理方法，分別計算河北省1982—2018年歷年馬鈴薯播種后全省平均積溫、全省標準化降水量和全省累積日照時數(shù)。

1.2.2 馬鈴薯產(chǎn)量資料的處理

為了分析氣象條件變化對產(chǎn)量的影響，對馬鈴薯單產(chǎn)進行處理[14]：

式中，ΔYi為第i年和第i-1年馬鈴薯單產(chǎn)豐歉值，稱作馬鈴薯產(chǎn)量豐歉氣象影響指數(shù)，Yi和Yi-1分別為第i年和第i-1年河北省馬鈴薯單產(chǎn)值。

1.3 預(yù)報方法

由于相同地區(qū)的土壤肥力、土壤類型和農(nóng)業(yè)投入等條件在相鄰兩年相對穩(wěn)定，所以相鄰兩年作物單產(chǎn)的波動主要由氣象條件的差異引起的[24]。即：

式中，ΔX為相鄰兩年氣象要素的變化量，ΔX=Xi-Xi-1，Xi為當年氣象要素（平均積溫、標準化降水量和累積日照時數(shù)等）值，Xi-1為前一年對應(yīng)的氣象要素值。

如果兩個相鄰年份之間氣象條件的差異變化不大，即認為這兩個年份同一作物的產(chǎn)量變化也相近，所以可以通過利用氣象條件相似年的方法來預(yù)報作物產(chǎn)量。相似年篩選采用氣象因子綜合聚類方法，即計算預(yù)報年（k年）馬鈴薯播種至產(chǎn)量預(yù)報時間段某氣象要素的變化量與歷史上任意一年（i年）同時段該氣象要素變化量的歐氏距離（dik）和相關(guān)系數(shù)（rik），再計算預(yù)報年氣象因子與歷史上任意一年的相似度Cik）：

Cik數(shù)值越高，則相似度越高。

為了減少預(yù)報的偶然性，選取全省平均積溫、全省標準化降水量和全省累積日照時數(shù)3個氣象要素的前3個Cik最大值，分別作為溫度、降水和日照3個因子的最佳相似年、第二相似年和第三相似年，每年就得到9個相似年，再分別采用大概率法和加權(quán)平均法計算預(yù)報年馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)。

針對環(huán)境衛(wèi)生學中大氣部分知識的測試情況，實驗班和對照班的成績平均分分別為（9.07±0.59）分和（7.65±1.35）分。對兩班的成績采用“有序分類變量兩組獨立樣本的秩和檢驗”方法進行分析，發(fā)現(xiàn)兩班成績相比，差異具有統(tǒng)計學意義（由于n值較大，采用正態(tài)近似檢驗，Zc=3.951 6，大于臨界值3.290 5，P＜0.001）。實驗班的成績高于對照班。實驗班和對照班學生的各分數(shù)段人數(shù)情況見（表1）所示。

1.3.1 大概率法

大概率法是在所得到的9個相似年的基礎(chǔ)上，根據(jù)9個相似年中產(chǎn)量為增產(chǎn)年或減產(chǎn)年的出現(xiàn)幾率，以出現(xiàn)幾率較多的增產(chǎn)年或減產(chǎn)年的產(chǎn)量豐歉氣象影響指數(shù)平均值為預(yù)報值，即：

式中，ΔY預(yù)報為采用大概率法得到的馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)預(yù)報值，是9個相似年中馬鈴薯單產(chǎn)為增產(chǎn)年的氣象影響指數(shù)累加，∑ΔYi（-）為9個相似年中馬鈴薯單產(chǎn)為減產(chǎn)年的氣象影響指數(shù)累加，L、M分別為9個相似年中馬鈴薯單產(chǎn)為增產(chǎn)年或減產(chǎn)年的個數(shù)。

1.3.2 加權(quán)平均法

加權(quán)平均法是在所得到的9個相似年的基礎(chǔ)上，根據(jù)9個相似年中馬鈴薯單產(chǎn)為增產(chǎn)年或減產(chǎn)年出現(xiàn)幾率為權(quán)重，采用加權(quán)平均得到預(yù)報年的單產(chǎn)豐歉∑ΔYi（+）氣象影響指數(shù)，即：

1.3.3 準確率

利用上述大概率與加權(quán)平均法對馬鈴薯單產(chǎn)豐歉趨勢ΔY進行預(yù)報，如果預(yù)報的ΔY與實際符號一致，則預(yù)報所得馬鈴薯豐歉趨勢準確，反之不準確。根據(jù)公式（1）可以得到馬鈴薯單產(chǎn)預(yù)報值，有定量預(yù)報準確率=×100%。

2 結(jié)果與分析

2.1 相似年的確定

以2018年作為預(yù)報年為例，2018年不同起始預(yù)報時間的溫度、降水和日照3個氣象因子的最佳相似年、第二相似年和第三相似年見表1。同一起始預(yù)報時間不同的氣象因子的相似年既有相同也有差異，如起始預(yù)報時間為6月1日時，溫度因子相似年中有2015年、降水因子相似年中也有2015年，但相似程度不同，溫度因子中2015年為第三相似年，降水因子中2015年為最佳相似年；不同起始預(yù)報時間，相同氣象因子的相似年份也既有相同又有差異，如以6月1日、7月1日和8月1日為起報時間，溫度因子的相似年中均有1991年，但6月1日起報時為最佳相似年，7月1日起報時為第二相似年，8月1日起報時為第三相似年，相似程度有差異。

表1 2018年不同時間起始預(yù)報時各氣象因子的最佳、第二和第三相似年

2.2 基于大概率法的馬鈴薯產(chǎn)量預(yù)報試驗

表2列出了2001—2015年基于大概率法的河北省馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)趨勢預(yù)報結(jié)果。馬鈴薯歷年單產(chǎn)豐歉氣象影響指數(shù)實際值（ΔY實際）為-34.3%～49.1%，不同起報時間馬鈴薯歷年單產(chǎn)豐歉氣象影響指數(shù)預(yù)報值（ΔY預(yù)報）在-30.2%～56.3%。6月1日起報時，近15 a中馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)增減趨勢預(yù)報結(jié)果與實際變化趨勢一致的有6 a，分別為2015、2014、2008、2006、2005和2003年，占40.0%；7月1日起報時，增減趨勢預(yù)報結(jié)果與實際變化趨勢一致的有9 a，分別為2015、2009—2013、2007、2006和2001年，占60.0%；8月1日起報時，增減趨勢預(yù)報結(jié)果與實際變化趨勢一致的有13 a，分別為2014、2001—2012年，占86.7%。可見，隨著預(yù)報時間的后移，預(yù)報效果越好。2006年預(yù)報效果最好，3次預(yù)報中馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)增減趨勢預(yù)報結(jié)果與實際變化趨勢均一致。2014—2015、2007—2012年、2005、2003和2001年的預(yù)報效果次之，3次預(yù)報中有2次預(yù)報與實際一致。2013、2004、2002年預(yù)報效果最差，只有1次預(yù)報與實際一致。

表2 基于大概率法的近15 a河北省馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)趨勢預(yù)報結(jié)果 %

表3為近15 a基于大概率法的馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)定量預(yù)報準確率。由表3可見，6月1日、7月1日和8月1日3次預(yù)報時馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)預(yù)報準確率分別在54.5%～96.8%、68.5%～98.5%、67.0%～98.7%，3次平均預(yù)報準確率分別為79.3%、84.1%和89.2%，可見隨著馬鈴薯生育進程的推進，預(yù)報時間越后移，預(yù)報準確率越高。

表3 基于大概率法的近15 a河北省馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)定量預(yù)報準確率/%

2.3 基于加權(quán)平均法的馬鈴薯產(chǎn)量預(yù)報試驗

表4列取了近15 a基于加權(quán)平均法的河北省馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)趨勢預(yù)報結(jié)果。由表4可見，不同預(yù)報時間馬鈴薯歷年單產(chǎn)豐歉氣象影響指數(shù)預(yù)報值（ΔY預(yù)報）在-22.1%～65.2%。6月1日起報時，除2010、2004和2002年外，其余年份馬鈴薯ΔY預(yù)報增減趨勢預(yù)報結(jié)果與實際變化趨勢一致，預(yù)報準確率為80.0%；7月1日起報時，除2011年外，其余年份馬鈴薯 增減趨勢預(yù)報結(jié)果與實際變化趨勢一致，預(yù)報準確率為93.3%；8月1日起報時，除2008年外，其余年份馬鈴薯增減趨勢預(yù)報結(jié)果與實際變化趨勢一致，預(yù)報準確率為93.3%?？梢?月1日和8月1日起報時預(yù)報效果好于6月1日的。從年際間比較看，除2011、2010、2008、2004和2002年分別在8月1日、6月1日、7月1日、6月1日和6月1日起報時馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)增減趨勢預(yù)報結(jié)果與實際變化趨勢不一致外，其余10 a的每年3次預(yù)報均與實際一致，預(yù)報效果良好。

表4 基于加權(quán)平均法的近15年河北省馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)趨勢預(yù)報結(jié)果 %

表5為近15 a基于加權(quán)平均法的馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)定量預(yù)報準確率。6月1日、7月1日和8月1日3次預(yù)報時馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)定量預(yù)報準確率分別在68.5%～82.8%、82.6%～98.9%、77.5%～99.4%，3次預(yù)報平均準確率分別為98.4%、90.5%和91.8%，可見8月1日預(yù)報準確率高于7月1日，7月1日準確率高于6月1日。

表5 基于加權(quán)平均法的近15 a河北省馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)定量預(yù)報準確率/%

2.4 基于大概率法和基于加權(quán)平均法預(yù)報準確率的比較

從馬鈴薯歷年單產(chǎn)豐歉氣象影響指數(shù)增減趨勢預(yù)報準確率看，6月1日預(yù)報時，大概率法平均準確率為40.0%，加權(quán)平均法平均為80.0%；7月1日預(yù)報時，大概率法平均準確率為60.0%，加權(quán)平均法平均為93.3%；8月1日預(yù)報時，大概率法平均準確率為86.7%，加權(quán)平均法平均為93.3%。加權(quán)平均法在3次預(yù)報中平均準確率均≥80%，而大概率法僅在8月1日預(yù)報時平均準確率＞80%（圖1a），在15 a每年3次預(yù)報中，大概率法預(yù)報準確率平均為62.2%，加權(quán)平均法平均為88.9%。

從馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)定量預(yù)報準確率看，6月1日預(yù)報時，大概率法平均預(yù)報準確率為79.3%，加權(quán)平均法為88.4%；7月1日預(yù)報時，大概率法平均預(yù)報準確率為84.1%，加權(quán)平均法為90.5%；8月1日預(yù)報時，大概率法平均預(yù)報準確率為89.2%，加權(quán)平均法為91.8%（圖1b）。在15 a的每年3次預(yù)報中，大概率法預(yù)報準確率平均為84.2%，加權(quán)平均法平均為90.2%。

圖1 近15 a基于兩種方法的河北省馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)增減趨勢與定量預(yù)報準確率比較

加權(quán)平均法在馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)趨勢預(yù)報和定量預(yù)報上均優(yōu)于大概率法，兩種方法均是隨著預(yù)報時間的后移，預(yù)報準確率越高。

2.5 預(yù)報檢驗

對2016—2018年河北省馬鈴薯產(chǎn)量進行動態(tài)預(yù)報檢驗，檢驗結(jié)果見表6。從產(chǎn)量豐歉氣象影響指數(shù)預(yù)報值上看，大概率法在3 a的9次預(yù)報中，產(chǎn)量增減趨勢與實際趨勢不一致的有4次預(yù)報，分別為2018年7月1日、2017年6月1日和7月1日、2016年6月1日預(yù)報，而加權(quán)平均法在3 a的9次預(yù)報中，僅在2017年6月1日預(yù)報中產(chǎn)量增減趨勢與實際趨勢相反。從定量預(yù)報準確率上看，大概率法和加權(quán)平均法在3 a的9次預(yù)報中平均準確率分別為86.9%和97.4%，采用大概率法準確率≥90%的概率為44.4%，而加權(quán)平均法準確率≥90%的概率為100%，準確率≥95%的概率為88.9%。

表6 兩種方法動態(tài)預(yù)報檢驗結(jié)果 %

3 結(jié)論與討論

（1）利用大概率法和加權(quán)平均法對河北省2001—2015年馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)進行動態(tài)預(yù)報，在15 a內(nèi)每年3次預(yù)報中，大概率法僅8月1日預(yù)報時馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)趨勢預(yù)報平均準確率高于80%，而加權(quán)平均法則在每年的3次預(yù)報中平均準確率均高于80%。采用大概率法的馬鈴薯單產(chǎn)豐歉氣象影響指數(shù)增減趨勢預(yù)報準確率平均為62.2%，定量預(yù)報平均準確率為84.2%，而加權(quán)平均法則分別為88.9%和90.2%，加權(quán)平均法可滿足業(yè)務(wù)應(yīng)用。

（2）對2016—2018年河北省馬鈴薯單產(chǎn)動態(tài)預(yù)報檢驗結(jié)果表明，大概率法在3 a內(nèi)的9次預(yù)報中，馬鈴薯單產(chǎn)增減趨勢與實際不一致的有4次，加權(quán)平均法僅有1次。采用大概率法準確率≥90%的概率為44.4%，而加權(quán)平均法準確率≥90%的概率為100%。

（3）大概率法和加權(quán)平均法兩種方法均隨著預(yù)報時間的向后推移，預(yù)報準確率越來越高，且加權(quán)平均法在產(chǎn)量豐歉氣象影響指數(shù)增減趨勢、定量預(yù)報準確率上均高于大概率法。

以往的研究中，大多基于氣候適宜度、關(guān)鍵氣象因子等方法進行內(nèi)蒙古[25]、山西[26]、寧夏[27]、甘肅[28]等地的馬鈴薯產(chǎn)量預(yù)報，基于產(chǎn)量歷史豐歉氣象影響指數(shù)的馬鈴薯產(chǎn)量預(yù)報尚未見報道。本文基于產(chǎn)量歷史豐歉氣象影響指數(shù)，研究并對比了大概率法和加權(quán)平均法在河北省馬鈴薯產(chǎn)量預(yù)報中的效果，加權(quán)平均法預(yù)報準確率較高，基本能滿足業(yè)務(wù)預(yù)報的需要。本文僅考慮了氣象因子對馬鈴薯產(chǎn)量的影響，未考慮氣象因子的協(xié)同作用，且由于馬鈴薯產(chǎn)量的高低除受氣象條件影響外，還與土壤、肥力、作物品種和外界干預(yù)等因素有關(guān)，所以產(chǎn)量預(yù)報結(jié)果有一定的局限性，還應(yīng)在今后的工作中深入研究，以進一步提高預(yù)報準確率。