徐 敏, 徐經(jīng)緯, 高 蘋, 吳洪顏*, 羅曉春, 任義方

(1. 江蘇省氣象局, 南京 210008; 2. 南京信息工程大學(xué)氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室,氣候與環(huán)境變化國(guó)際合作聯(lián)合實(shí)驗(yàn)室,氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,資料同化研究與應(yīng)用中心, 南京 210044)

?

基于大尺度因子的江蘇稻區(qū)稻瘟病氣象等級(jí)長(zhǎng)期預(yù)測(cè)

徐 敏1, 徐經(jīng)緯2, 高 蘋1, 吳洪顏1*, 羅曉春1, 任義方1

(1. 江蘇省氣象局, 南京 210008; 2. 南京信息工程大學(xué)氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室,氣候與環(huán)境變化國(guó)際合作聯(lián)合實(shí)驗(yàn)室,氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,資料同化研究與應(yīng)用中心, 南京 210044)

為較早地準(zhǔn)確預(yù)測(cè)水稻稻瘟病發(fā)生發(fā)展的氣象等級(jí),利用大氣環(huán)流和太平洋海溫對(duì)氣象條件影響的滯后性,采用最優(yōu)相關(guān)和空間拓?fù)浞治黾夹g(shù),結(jié)合滑動(dòng)平均和主成分識(shí)別法,篩選出對(duì)江蘇稻區(qū)稻瘟病指數(shù)影響最顯著且穩(wěn)定獨(dú)立的大尺度預(yù)報(bào)因子,分別建立了基于大氣環(huán)流因子和基于海溫因子的稻瘟病氣象等級(jí)長(zhǎng)期預(yù)測(cè)模型。經(jīng)歷史擬合和試報(bào)檢驗(yàn),模型效果理想,能提前一個(gè)月預(yù)測(cè)出水稻稻瘟病發(fā)生的氣象等級(jí)。該模型的預(yù)測(cè)結(jié)果對(duì)江蘇稻區(qū)稻瘟病防治具有重要意義。

水稻稻瘟病; 氣象等級(jí)預(yù)測(cè); 海溫; 大氣環(huán)流指數(shù)

水稻稻瘟病是由被稱為灰梨孢或稻梨孢的真菌引起,該真菌屬于半知菌亞門[1]。稻瘟病菌不僅能侵染水稻,還能侵染小麥、大麥等農(nóng)作物的莖葉,以及植物的根部[2]。稻瘟病是我國(guó)稻區(qū)最嚴(yán)重的病害之一,一般流行年份可造成10%～20%的減產(chǎn),嚴(yán)重發(fā)生時(shí)損失高達(dá)40%～50%[3],嚴(yán)重制約我國(guó)水稻生產(chǎn)。近年來(lái)在我國(guó)南北稻區(qū)及部分感病品種上呈偏重流行態(tài)勢(shì),老病區(qū)和感病品種普遍發(fā)病,具有點(diǎn)多面廣、發(fā)病品種多、發(fā)生頻率趨高、面積趨大、為害趨重的特點(diǎn)[4-6]。

稻瘟病的發(fā)生流行受菌源、氣象因子及寄主植物水稻三者影響,因此可根據(jù)菌量、氣象條件、栽培條件或寄主生育狀況、侵染狀況進(jìn)行預(yù)測(cè)[2]。稻瘟病屬于再侵染頻繁的病害,初始菌量的多少?zèng)Q定流行程度,氣象條件也對(duì)其影響很大,溫度主要影響稻瘟病菌的生長(zhǎng)發(fā)育,濕度則影響病原菌孢子的形成、萌發(fā)和侵入,雨水飛濺易致氣傳孢子的釋放和傳播。王建忠等[7]指出稻瘟病的發(fā)生、流行及危害與氣象條件密切相關(guān),屬于典型的“氣象型”病害,在稻瘟病流行年份,水稻生長(zhǎng)前期的氣候特點(diǎn)主要表現(xiàn)為溫暖、多雨寡照;水稻生長(zhǎng)中后期的氣候特點(diǎn)表現(xiàn)為低溫、多雨寡照[8-9]。因此,利用氣象條件與稻瘟病間的關(guān)系,建立稻瘟病氣象等級(jí)預(yù)測(cè)模型是可行的。我國(guó)學(xué)者從20世紀(jì)80年代起,在農(nóng)作物病害發(fā)生發(fā)展與氣象條件的關(guān)系、發(fā)生的氣象環(huán)境成因、氣象預(yù)測(cè)預(yù)報(bào)等研究方面取得了重要進(jìn)展,何永坤等[10]從稻瘟病發(fā)病機(jī)理出發(fā)，建立了其發(fā)生發(fā)展的氣象條件促病指數(shù)和氣象條件等級(jí)預(yù)報(bào)指標(biāo),利用72 h的短期天氣預(yù)報(bào)開(kāi)展稻瘟病發(fā)生發(fā)展的氣象等級(jí)預(yù)報(bào);郭瑞鴿等[11]針對(duì)江西早稻利用統(tǒng)計(jì)方法建立了稻瘟病發(fā)生的氣象條件監(jiān)測(cè)預(yù)警方法,開(kāi)展預(yù)警時(shí)同樣采用的是短期天氣預(yù)報(bào);黃春艷等[12]應(yīng)用逐步回歸法建立了黑龍江稻區(qū)稻瘟病流行趨勢(shì)預(yù)測(cè)數(shù)學(xué)模型,預(yù)報(bào)因子為6月下旬至7月中旬的各旬關(guān)鍵氣象因子;胡毅等[13]利用貝葉斯判別法和神經(jīng)網(wǎng)絡(luò)建立了四川省稻區(qū)稻瘟病發(fā)病率及其等級(jí)預(yù)報(bào)模式,預(yù)報(bào)因子為3月上旬至6月上旬各旬的主要?dú)庀笠亍?/p>

從已有研究來(lái)看,稻瘟病氣象等級(jí)的預(yù)測(cè)時(shí)效基本上只有3 d左右,預(yù)測(cè)精度依賴于短期天氣預(yù)報(bào)的準(zhǔn)確度,預(yù)測(cè)時(shí)效短不利于稻瘟病防治工作的開(kāi)展。為了延長(zhǎng)預(yù)測(cè)時(shí)效,本研究將以江蘇稻區(qū)為例,在綜合稻瘟病指數(shù)和稻瘟病氣象等級(jí)指標(biāo)的研究成果基礎(chǔ)上,根據(jù)長(zhǎng)期預(yù)報(bào)方法,利用最優(yōu)相關(guān)和空間拓?fù)浞治黾夹g(shù),篩選出對(duì)綜合稻瘟病指數(shù)影響顯著的海溫預(yù)報(bào)因子和大氣環(huán)流預(yù)報(bào)因子,分別建立其長(zhǎng)期預(yù)報(bào)模型,并進(jìn)行模型檢驗(yàn)。若能提前一個(gè)月且準(zhǔn)確地預(yù)測(cè)出稻瘟病氣象等級(jí),則可為提前采取有效防治措施提供科學(xué)依據(jù),同時(shí)對(duì)保障糧食安全具有十分重要的現(xiàn)實(shí)意義。

1 資料與方法

1.1 研究資料

(1) 氣象資料:1978-2015年,江蘇全省67個(gè)氣象站逐日氣象資料,日降水量、日平均溫度、日均相對(duì)濕度、日照時(shí)數(shù)。

(2) 環(huán)流指數(shù):1978-2015年,來(lái)自國(guó)家氣候中心的74項(xiàng)逐月大氣環(huán)流指數(shù),包括不同區(qū)域副熱帶高壓面積和強(qiáng)度指數(shù)、極渦面積指數(shù)、緯向環(huán)流指數(shù)、經(jīng)向環(huán)流指數(shù)、南方濤動(dòng)指數(shù)、東亞槽強(qiáng)度和位置等。

(3) 海溫資料:1978-2015年,來(lái)自美國(guó)大氣與海洋管理局網(wǎng)站氣候診斷中心的太平洋海區(qū)(10°S～50°N、120°E～80°W)逐月海表面溫度,水平分辨率是5°×5°(圖1)。

圖1 1978-2015年北太平洋海溫場(chǎng)286個(gè)格點(diǎn)編號(hào)示意圖Fig.1 The north Pacific SST field 286 lattice code diagram

1.2 研究方法

1.2.1 綜合稻瘟病指數(shù)的算法

適宜的氣象條件是水稻稻瘟病大發(fā)生的直接誘發(fā)因素,當(dāng)氣溫在20～30℃、空氣相對(duì)濕度90%以上、稻株表面水膜保持6～10 h,稻瘟病就容易發(fā)生[10]。對(duì)于江蘇稻區(qū),6—8月為稻瘟病發(fā)生流行的主要時(shí)期,該時(shí)段覆蓋了水稻的葉瘟和穗頸瘟發(fā)病期,因此在構(gòu)建稻瘟病指數(shù)時(shí),綜合考慮了葉瘟指數(shù)和穗頸瘟指數(shù),并對(duì)綜合稻瘟病指數(shù)所對(duì)應(yīng)的氣象等級(jí)進(jìn)行了劃分,具體計(jì)算公式和等級(jí)劃分見(jiàn)文獻(xiàn)[14]。

1.2.2 最優(yōu)化算法

為了實(shí)現(xiàn)提前一個(gè)月預(yù)測(cè)出稻瘟病氣象等級(jí)的目標(biāo),充分利用稻瘟病發(fā)生的氣象條件對(duì)大氣環(huán)流與海溫響應(yīng)的滯后性[15-16],采用相關(guān)分析篩選出與綜合稻瘟病指數(shù)相關(guān)的大氣環(huán)流指數(shù)和海溫格點(diǎn)。為了取得更加顯著的相關(guān)效果,應(yīng)用最優(yōu)化相關(guān)處理技術(shù)[17]進(jìn)行膨化處理,經(jīng)過(guò)最優(yōu)化處理可尋找到相關(guān)顯著的因子。

最優(yōu)化處理計(jì)算公式如下:

H=(|x-b|/B+0.5)a

(1)

其中,x為大尺度預(yù)報(bào)因子(大氣環(huán)流指數(shù)或太平洋海溫),H為經(jīng)過(guò)變換的大尺度預(yù)報(bào)因子,a、b為待定參數(shù),根據(jù)經(jīng)驗(yàn),a、b、B取值范圍如下:

a∈(-10,-1/10)或a∈(1/10,10)

xmin+(xmax-xmin)/4≤b≤xmax-(xmax-xmin)/4

B=max(xmax-b,b-xmin)

待定參數(shù)a、b用最優(yōu)化技術(shù)求出,令目標(biāo)函數(shù)為:

f(a,b)=1-R2→min

(2)

其中,R為變換后的大尺度預(yù)報(bào)因子H與綜合稻瘟病指數(shù)的相關(guān)系數(shù),經(jīng)過(guò)最優(yōu)化處理后,可得到一系列與綜合稻瘟病指數(shù)相關(guān)最顯著的大尺度預(yù)報(bào)因子。

2 結(jié)果與分析

2.1 基于大尺度氣象因子的綜合稻瘟病指數(shù)預(yù)測(cè)模型

2.1.1 大氣環(huán)流預(yù)報(bào)因子篩選與預(yù)測(cè)模型的建立

由于江蘇地區(qū)稻瘟病的易發(fā)關(guān)鍵期是8月份,因此預(yù)報(bào)因子截止到7月份,通過(guò)因子膨化,從上一年1月至當(dāng)年7月共19個(gè)月,對(duì)74項(xiàng)大氣環(huán)流特征量進(jìn)行所有時(shí)段組合,并計(jì)算各組合時(shí)段的平均值,經(jīng)過(guò)篩選,得到21個(gè)兼具穩(wěn)定性和獨(dú)立性的大氣環(huán)流影響因子,具體為:不同時(shí)間段北非大西洋北美副高面積指數(shù)(110°W～60°E)、東太平洋副高面積指數(shù)(175°W～115°W)、西太平洋副高強(qiáng)度指數(shù)(110°W～180°)、南海副高強(qiáng)度指數(shù)(100°E～120°E)、太平洋副高強(qiáng)度指數(shù)(110°E～115°W)、不同時(shí)間段北半球副高脊線(5°E～360°)、北半球副高脊線(5°E～360°)、北半球副高北界(5°E～360°)、不同時(shí)間段太平洋副高北界(110°E～115°W)、太平洋區(qū)極渦面積指數(shù)(2區(qū)150°E～120°W)、北美區(qū)極渦面積指數(shù)(3區(qū)120°W～30°W)、北半球區(qū)極渦強(qiáng)度指數(shù)(5區(qū)0～360°)、不同時(shí)間段大西洋歐洲環(huán)流型C、亞洲緯向和經(jīng)向環(huán)流指數(shù)(60°E～150°E)、印緬槽(15°N～20°N,80°E～100°E)、編號(hào)臺(tái)風(fēng)。

對(duì)上述篩選出的顯著相關(guān)因子,進(jìn)行逐步回歸,由其中7個(gè)因子建立了基于大氣環(huán)流因子的綜合稻瘟病指數(shù)的預(yù)測(cè)模型,方程如下:

Zhl=-0.539H1+0.685H2-0.490H3+0.364H4+

0.231H5-0.914H6-3.150H7+14.267

式中,Zhl為水稻綜合稻瘟病指數(shù),H1……H7為經(jīng)過(guò)線性和非線性處理的大氣環(huán)流預(yù)報(bào)因子,計(jì)算公式和含義見(jiàn)表1。

表1 綜合稻瘟病指數(shù)預(yù)測(cè)模型中的大氣環(huán)流預(yù)報(bào)因子1)

Table 1 The mathematical expression of atmospheric circulation predictor in the long-term prediction model of composite index of rice blast

預(yù)報(bào)因子Predictor表達(dá)式Expression影響時(shí)段Periodaffected含義Symbolicmeaning相關(guān)系數(shù)CorrelationcoefficientH1(|x1-14.64|/3.36+0.5)2.97上年5月Maylastyear北半球副高脊線Subtropicalhighridgelineofthenorthernhemisphere-0.421H2(|x2-10.19|/4.81+0.5)-2.46上年1月-上年3月January-Marchlastyear北半球副高脊線Subtropicalhighridgelineofthenorthernhemisphere0.403H3(|x3-216.56|/37.44+0.5)-2.27上年12月Decemberlastyear太平洋區(qū)極渦面積指數(shù)PolarvortexareaindexofthePacificarea-0.404H4(|x4-7.12|/7.88+0.5)-2.77當(dāng)年5月-當(dāng)年7月May-Julythepresentyear大西洋歐洲環(huán)流型EuropeanAtlanticcirculationtype0.453

續(xù)表1 Table 1(Continued)

預(yù)報(bào)因子Predictor表達(dá)式Expression影響時(shí)段Periodaffected含義Symbolicmeaning相關(guān)系數(shù)CorrelationcoefficientH5(|x5-95.58|/27.25+0.5)-2.61上年5月-上年7月May-Julylastyear亞洲緯向環(huán)流指數(shù)Asianzonalcirculationindex0.506H6(|x6-41.50|/11+0.5)-1.50當(dāng)年6月-當(dāng)年7月June-Julythepresentyear亞洲經(jīng)向環(huán)流指數(shù)Asianmeridionalcirculationindex-0.547H7(|x7-16|/26+0.5)0.75當(dāng)年3月Marchthepresentyear印緬槽Indian-Myanmartrough-0.506

1) 通過(guò)0.01顯著性檢驗(yàn)的F值是0.418,通過(guò)0.001顯著性檢驗(yàn)的F值是0.519,x1……x7是對(duì)應(yīng)的大氣環(huán)流指數(shù)值。Fvalue is 0.418 by significance test of 0.01.Fvalue is 0.519 by significance test of 0.001.x1……x7is the atmospheric circulation index value.

2.1.2 海溫預(yù)報(bào)因子篩選與預(yù)測(cè)模型的建立

海洋對(duì)我國(guó)氣候的影響存在后延效應(yīng),因此我們計(jì)算了上年1月至當(dāng)年7月海溫與綜合稻瘟病指數(shù)的相關(guān)系數(shù),受篇幅限制,文中以相關(guān)性較好的4個(gè)月份為例(圖2),從海溫場(chǎng)的相關(guān)系數(shù)空間分布圖上可以看出,綜合稻瘟病指數(shù)與上年5月(圖2a)、6月(圖2b)、7月(圖2c)、9月(圖2d)太平洋海溫均存在高相關(guān)區(qū)(藍(lán)色區(qū)域,通過(guò)0.05顯著性檢驗(yàn)),上年5月的顯著負(fù)相關(guān)區(qū)位于(23°N～27°N,150°E～180°),上年6月和7月的顯著負(fù)相關(guān)區(qū)域較為接近,均位于25°N附近的東太平洋區(qū)域,上年9月的顯著負(fù)相關(guān)區(qū)已經(jīng)轉(zhuǎn)移到赤道以南的西太平洋區(qū)域。由此可見(jiàn),太平洋海溫對(duì)綜合稻瘟病指數(shù)具有可預(yù)報(bào)性,而且海溫的后延效應(yīng)可大幅提高預(yù)測(cè)模型的預(yù)測(cè)時(shí)效。

圖2 綜合稻瘟病指數(shù)與太平洋海溫的相關(guān)系數(shù)Fig.2 Correlation coefficient between composite index of rice blast and sea surface temperature

由于太平洋海溫是一個(gè)場(chǎng),若是獨(dú)立的格點(diǎn)具有相關(guān)性并沒(méi)有實(shí)際意義,所以在尋找綜合稻瘟病指數(shù)的海溫相關(guān)因子時(shí)需要考慮空間相關(guān)格點(diǎn)數(shù),在此利用空間拓?fù)浞治鲈?篩選出空間相連格點(diǎn)≥4個(gè)的海區(qū),將此海區(qū)影響時(shí)段的海溫平均值作為一個(gè)海溫預(yù)報(bào)因子,然后再將所選因子進(jìn)行最優(yōu)化處理。經(jīng)過(guò)篩選,共有22個(gè)海溫因子入選,經(jīng)過(guò)逐步回歸后,由其中7個(gè)因子構(gòu)建了預(yù)測(cè)模型,方程如下:

Zsst=-1.935S1-0.284S2+0.003S3+15.276S4-

0.770S5+0.133S6-0.338S7-0.544

式中,Zsst為水稻綜合稻瘟病指數(shù),S1……S7為經(jīng)過(guò)線性和非線性處理的海溫預(yù)報(bào)因子,計(jì)算公式和相關(guān)信息見(jiàn)表2。

表2 綜合稻瘟病指數(shù)預(yù)測(cè)模型中的海溫預(yù)報(bào)因子1)

Table 2 The mathematical expression of sea surface temperature predictor in the long-term prediction model of composite index of rice blast

預(yù)報(bào)因子Predictor表達(dá)式Expression影響時(shí)段Periodaffected高相關(guān)海區(qū)格點(diǎn)數(shù)Latticesitesinhighlycorrelatedarea相關(guān)系數(shù)CorrelationcoefficientS1(|y1-136.75|/30+0.5)0.82上年5月Maylastyear4-0.479S2(|y2-299.45|/17.45+0.5)-3.54上年7月Julylastyear5-0.543S3(|y3-246.98|/12.13+0.5)-10上年9月Septemberlastyear90.417S4(|y4-232.13|/10.27+0.5)0.125上年3月Marchlastyear70.559S5(|y5-279.11|/18.31+0.5)-2.12上年4月Aprillastyear5-0.516S6(|y6-274.29|/20.44+0.5)-4.24上年6月Junelastyear110.525S7(|y7-167.59|/37.88+0.5)4.56上年7月Julylastyear170.514

1) 通過(guò)0.01顯著性檢驗(yàn)的F值是0.418,通過(guò)0.001顯著性檢驗(yàn)的F值是0.519,y1……y7是對(duì)應(yīng)海區(qū)海溫格點(diǎn)的平均值。Fvalue is 0.418 by significance test of 0.01.Fvalue is 0.519 by significance test of 0.001.y1……y7is the atmospheric circulation index value.

2.2 模型的擬合檢驗(yàn)和試報(bào)檢驗(yàn)

為了檢驗(yàn)基于大氣環(huán)流因子和基于海溫因子的預(yù)測(cè)模型的擬合效果,將綜合稻瘟病指數(shù)的實(shí)際值和模擬值進(jìn)行了對(duì)比(圖3),可以看出兩類預(yù)測(cè)模型的模擬結(jié)果與綜合稻瘟病指數(shù)的實(shí)際值都非常接近,波動(dòng)特征基本一致,尤其是病害重的極端年份,1980、1986、1987、2003、2011年的綜合稻瘟病指數(shù)均超過(guò)12,對(duì)照表1,這些年份稻瘟病發(fā)生的氣象條件等級(jí)為強(qiáng)適宜,兩類模型均能準(zhǔn)確模擬出稻瘟病發(fā)生程度的氣象等級(jí)?；诖髿猸h(huán)流因子的預(yù)測(cè)模型的擬合值平均相對(duì)誤差是7.8%,基于海溫因子的預(yù)測(cè)模型的擬合值平均相對(duì)誤差是7.5%,海溫模型的誤差相對(duì)小一些。兩類模型均通過(guò)了0.01的顯著性檢驗(yàn),擬合效果好。

圖3 1978-2012年綜合稻瘟病指數(shù)兩類長(zhǎng)期預(yù)測(cè)模型的擬合檢驗(yàn)Fig.3 Fitting test of two long-term prediction models of composite index of rice blast during 1978-2012

利用2013-2015年的資料進(jìn)行試報(bào)檢驗(yàn),從表4可以看出,利用大氣環(huán)流因子構(gòu)建的預(yù)測(cè)模型預(yù)測(cè)出的近3年綜合稻瘟病指數(shù)與實(shí)際值的誤差范圍是0.18～0.67,其中2014年預(yù)測(cè)出的氣象條件對(duì)病害的適宜程度為“弱適宜”,低于實(shí)際的氣象等級(jí)“適宜”,其余兩年的氣象等級(jí)預(yù)報(bào)均為正確;利用海溫因子構(gòu)建的預(yù)測(cè)模型預(yù)測(cè)出的近3年綜合稻瘟病指數(shù)與實(shí)際值的誤差范圍是-0.34～0.28,預(yù)測(cè)值與實(shí)際值非常接近,預(yù)測(cè)的氣象條件對(duì)稻瘟病病害的適宜程度均正確。從近3年試報(bào)效果來(lái)看,基于海溫因子的模型預(yù)測(cè)效果要好于基于環(huán)流因子的模型預(yù)測(cè)效果。

表3 綜合稻瘟病指數(shù)預(yù)測(cè)模型的試報(bào)檢驗(yàn)(綜合稻瘟病指數(shù)/氣象等級(jí))

Table 3 Test of two long-term prediction models of composite index of rice blast (composite index of rice blast / meteorology grade)

年份Year基于環(huán)流因子的預(yù)報(bào)值Forecastvalueofpredictionmodelbasedonthecirculationfactors基于海溫因子的預(yù)報(bào)值Forecastvalueofpredictionmodelbasedonseasurfacetemperaturefactors實(shí)際值A(chǔ)ctualvalue20138.89/(1)8.79/(1)9.07/(1)20149.97/(1)10.68/(2)10.35/(2)201511.14/(2)11.74/(2)11.81/(2)

3 結(jié)論與討論

本研究為稻瘟病氣象等級(jí)的預(yù)測(cè)研究提供了新思路,通過(guò)找尋稻瘟病的大尺度氣象預(yù)報(bào)因子,構(gòu)建了稻瘟病氣象等級(jí)的長(zhǎng)期預(yù)測(cè)模型,并對(duì)模型進(jìn)行了驗(yàn)證,通過(guò)了0.01顯著性水平檢驗(yàn),說(shuō)明預(yù)測(cè)模型切實(shí)可行,可提前一個(gè)月預(yù)測(cè)出稻瘟病的氣象等級(jí),為提前準(zhǔn)備稻瘟病的科學(xué)防治提供了充足時(shí)間,彌補(bǔ)了原有稻瘟病氣象等級(jí)預(yù)測(cè)時(shí)效短的不足。大尺度氣象預(yù)報(bào)因子的找尋,一是根據(jù)中長(zhǎng)期天氣預(yù)報(bào)原理,海溫的變化會(huì)引起大氣環(huán)流系統(tǒng)的改變[18],從而影響到各地的氣象條件,而氣象條件的變化會(huì)影響稻瘟病的發(fā)生發(fā)展;二是充分利用了海溫和大氣環(huán)流對(duì)局地氣象條件影響的滯后性,使得預(yù)測(cè)模型的預(yù)測(cè)時(shí)效延伸到了一個(gè)月。

選育和種植抗病品種雖然是控制稻瘟病發(fā)生流行的最經(jīng)濟(jì)、環(huán)保的措施,但在一定年限后,選育出的抗病品種仍可能變成感病品種。因此,當(dāng)?shù)疚敛“l(fā)生流行時(shí),速效的化學(xué)防治仍然是最佳手段,有效防治稻瘟病的前提是進(jìn)行有效的預(yù)測(cè),只有及時(shí)準(zhǔn)確地掌握了稻瘟病發(fā)生發(fā)展的趨勢(shì),才能抓住用藥最佳時(shí)機(jī),實(shí)現(xiàn)有效防治[19]。值得注意的是,通過(guò)本研究構(gòu)建的稻瘟病氣象等級(jí)長(zhǎng)期預(yù)測(cè)模型,預(yù)測(cè)出的氣象等級(jí)與稻瘟病發(fā)生流行實(shí)際等級(jí)可能會(huì)存在不一致的情形,因?yàn)榈疚敛庀蟮燃?jí)是表示氣象因子對(duì)稻瘟病發(fā)生發(fā)展的適宜程度,而實(shí)際田間的稻瘟病除了受氣象條件影響外,還與植株抗病能力、病原數(shù)量、氮肥施用量和施用時(shí)期等有關(guān)。因此,在開(kāi)展稻瘟病的實(shí)際預(yù)測(cè)時(shí),需要綜合考慮氣象因素和管理水平以及水稻本身的抗病能力。

[1] 任鄄勝,肖陪村,陳勇,等.水稻稻瘟病病菌研究進(jìn)展[J].現(xiàn)代農(nóng)業(yè)科學(xué),2008,15(1):19-23.

[2] 曾士邁,楊演.植物病害流行學(xué)[M].北京:農(nóng)業(yè)出版社,1986:184-190.

[3] 孫國(guó)昌,杜新法,陶榮祥.水稻稻瘟病防治研究進(jìn)展和21世紀(jì)研究設(shè)想[J].植物保護(hù),2000,26(1):33-35.

[4] 王艷青.近年來(lái)中國(guó)水稻病蟲(chóng)害發(fā)生及趨勢(shì)分析[J]. 中國(guó)農(nóng)學(xué)通報(bào),2006,22(2):343-347.

[5] 趙夢(mèng),歐陽(yáng)芳,張永生,等.2000-2010年我國(guó)水稻病蟲(chóng)害發(fā)生與為害特征分析[J]. 生物災(zāi)害科學(xué),2014,37(4):275-280.

[6] 舒暢,潘華. 江西水稻病蟲(chóng)害發(fā)生與防治淺議[J]. 氣象與減災(zāi)研究,2006,29(1):61-64.

[7] 王建忠,鄭文靜,張燕之. 水稻稻瘟病的發(fā)病原因及綜合防治[J]. 農(nóng)業(yè)科技通訊, 2006(8):24-25.

[8] 劉靜,王健林,宋迎波,等. 寧夏水稻稻瘟病發(fā)生程度的氣象等級(jí)預(yù)報(bào)[J]. 安徽農(nóng)業(yè)科學(xué),2009, 37(1):5021-5023.

[9] 劉志恒,王世維,魏松紅,等. 2011—2012年遼寧省稻瘟病菌種群動(dòng)態(tài)分析[J]. 沈陽(yáng)農(nóng)業(yè)大學(xué)學(xué)報(bào),2014,45(4):393-397.

[10]何永坤,陽(yáng)園燕,羅孳孳. 稻瘟病發(fā)生發(fā)展氣象條件等級(jí)業(yè)務(wù)預(yù)報(bào)技術(shù)研究[J]. 氣象,2008,34(12):110-113.

[11]郭瑞鴿,杜筱玲,劉文英,等. 江西早稻瘟病發(fā)生的氣象條件監(jiān)測(cè)預(yù)警方法[J]. 氣象與減災(zāi)研究,2008, 31(3):48-51.

[12]黃春艷,張勻華,商世吉,等. 黑龍江稻區(qū)稻瘟病流行趨勢(shì)預(yù)測(cè)數(shù)學(xué)模型研究[J]. 植物保護(hù),1995,21(2):8-11.

[13]胡毅,向衛(wèi)國(guó),趙純波,等. 四川省稻瘟病發(fā)病率及其等級(jí)預(yù)報(bào)初探[J]. 成都?xì)庀髮W(xué)院學(xué)報(bào),1994,9(4):60-67.

[14]任義方,高蘋,朱鳳. 江蘇省稻瘟病發(fā)生程度的氣象等級(jí)預(yù)報(bào)[J]. 江蘇農(nóng)業(yè)科學(xué),2016,44(8):151-154.

[15]高蘋,武金崗,楊榮明,等.江蘇省稻縱卷葉螟遷入期蟲(chóng)情指標(biāo)與西太平洋海溫的遙相關(guān)及其長(zhǎng)期預(yù)報(bào)模型的研究[J].應(yīng)用生態(tài)學(xué)報(bào),2008,19(8):2056-2066.

[16]陳烈庭.太平洋海氣相互作用時(shí)空變化[J].氣象學(xué)報(bào),1983,41:296-304.

[17]湯志成, 孫涵.最優(yōu)化因子處理及加權(quán)多重回歸模型[J]. 氣象學(xué)報(bào),1992, 50(4):514-517.

[18]章基嘉,葛玲.中長(zhǎng)期天氣預(yù)報(bào)基礎(chǔ)[M].北京:氣象出版社,1983:254-282.

[19]齊龍,馬旭,梁柏,等.稻瘟病監(jiān)測(cè)預(yù)測(cè)方法研究現(xiàn)狀及流行風(fēng)險(xiǎn)評(píng)估體系構(gòu)建初探[J].中國(guó)農(nóng)學(xué)通報(bào),2011,27(33):213-216.

(責(zé)任編輯：田 喆)

Long-term prediction models based on large-scale factors for meteorological grade of the rice blast in Jiangsu rice area

Xu Min1, Xu Jingwei2, Gao Ping1, Wu Hongyan1, Luo Xiaochun1, Ren Yifang1

(1.MeteorologicalBureauofJiangsuProvince,Nanjing210008,China; 2.KeyLaboratoryofMeteorologicalDisaster,MinistryofEducation(KLME),JointInternationalResearchLaboratoryofClimateandEnvironmentChange(ILECE),CollaborativeInnovationCenteronForecastandEvaluationofMeteorologicalDisasters(CIC-FEMD),CenterforDataAssimilationResearchandApplication(CDARA),NanjingUniversityofInformationSciences&Technology(NUIST),Nanjing210044,China)

In order to predict the degree of meteorological grade of rice blast early and accurately, two long-term prediction models of meteorological grade for rice blast were established based on the atmospheric circulation factors and based on SST factors by the optimal correlation and spatial topological analysis technology and the impact of atmospheric circulation and Pacific SST on meteorological conditions (tread effect). In addition, the methods of moving average and principal component were used for the forecast factor having stability and independence to the comprehensive index of rice blast. Furthermore, through the historical fitting and forecast testing, the results of the prediction models were satisfactory. The meteorological grade of rice blast could be predicted by the models at least one month in advance. Therefore, the prediction results have a significant guiding meaning for agricultural prevention and control of rice blast in Jiangsu rice area.

rice blast; prediction model of meteorological grade; sea surface temperature; atmospheric circulation index

2016-09-02

2016-12-07

公益性行業(yè)(氣象)科研專項(xiàng)(GYHY201306035);江蘇省氣象局科研基金面上項(xiàng)目(KM201504,KM201707)

S431.2, S 435.111.41

A

10.3969/j.issn.0529-1542.2017.04.006

* 通信作者 E-mail: jsqx_zqzb@126.com