煙草馬鈴薯Y 病毒病與氣象因子的相關(guān)性

付修廷，游堂貴，李曉燕，黃，曾廳余，倪 霞，陸文林，何湘華

1.云南省煙草公司昭通市公司，云南省昭通市昭陽(yáng)區(qū) 657000

2.云南省昭通市氣象局，云南省昭通市昭陽(yáng)區(qū) 657000

煙草馬鈴薯Y 病毒?。?］是由馬鈴薯Y 病毒引起的煙草上的一種病毒病害［2］，在全國(guó)各煙區(qū)普遍發(fā)生［3］；發(fā)病早、病歷長(zhǎng)、范圍廣、損失大，該病的發(fā)生除引起煙草產(chǎn)量降低外，更為嚴(yán)重的是影響煙葉品質(zhì)，帶有煙草馬鈴薯Y 病毒病的烤煙病葉烤后色澤和品質(zhì)顯著降低，已成為影響昭通烤煙生產(chǎn)的主要病害之一。在煙草馬鈴薯Y 病毒病發(fā)生規(guī)律及防治方面，前人已有研究，余春英等［2］提出由于氣候條件的變化，馬鈴薯種植面積增大以及抗病品種單一導(dǎo)致煙草馬鈴薯Y 病毒病由次要病害逐漸上升為主要病害；孫宏宇等［3］認(rèn)為蚜蟲的發(fā)生量是影響煙草馬鈴薯Y 病毒病發(fā)生程度的重要因素。但針對(duì)煙草馬鈴薯Y 病毒病在昭通發(fā)生流行規(guī)律及其與當(dāng)?shù)貧夂蛞蛩氐南嚓P(guān)性分析方面尚未見報(bào)道。因此，依據(jù)作物病蟲害氣象條件預(yù)測(cè)分析方面的研究方法［4-7］，分析了昭通氣候因素對(duì)煙草馬鈴薯Y 病毒病發(fā)生的影響，旨在明確馬鈴薯Y 病毒病感染煙草的生育時(shí)段和發(fā)生流行的適宜氣候條件，為煙草馬鈴薯Y 病毒病的監(jiān)測(cè)、預(yù)測(cè)以及有效防治提供依據(jù)。

1 材料和方法

1.1 材料

2009—2011年4—8 月逐日氣溫、降水、日照資料來(lái)源于昭陽(yáng)區(qū)氣象局；煙草發(fā)病株率及病情指數(shù)數(shù)據(jù)來(lái)源于云南省煙草公司昭通市公司煙葉生產(chǎn)技術(shù)中心，2009—2011年5 月25 日—8 月30 日每年在昭陽(yáng)區(qū)太平辦事處水平村社區(qū)煙田的調(diào)查觀測(cè)數(shù)據(jù)見表1。

表1 2009—2011年煙草馬鈴薯Y 病毒病調(diào)查觀測(cè)結(jié)果（%）

1.2 調(diào)查方法

烤煙團(tuán)棵期開始到煙葉成熟采收結(jié)束，田間固定5點(diǎn)取樣，每點(diǎn)不少于30 株，每5 天調(diào)查1 次，根據(jù)國(guó)家標(biāo)準(zhǔn)方法［8］，計(jì)算發(fā)病株率及病情指數(shù)。

煙草馬鈴薯Y 病毒病的分級(jí)［8］：0 級(jí)全株無(wú)?。? 級(jí)心葉脈明或輕微花葉，病株無(wú)明顯矮化；3 級(jí)1/3 葉片花葉但不變形，或病株老化為正常株高的3/4 以上；5 級(jí)1/3 至1/2 葉片花葉，或少數(shù)葉變形，或主脈邊黑，或病株矮化為正常株高 的2/3 或3/4；7 級(jí)1/2 至1/3 葉片 花葉，或變形或主側(cè)脈壞死，或病株矮化為正常株高的1/2至2/3。

1.3 數(shù)據(jù)處理

采用多因子逐步回歸方法進(jìn)行數(shù)據(jù)的統(tǒng)計(jì)分析。將發(fā)病株率（Yi）、病情指數(shù)（Yj）設(shè)為預(yù)測(cè)變量，各時(shí)段平均氣溫（Ti）、累計(jì)雨量（Ri）、累計(jì)日照（Si）設(shè)為因變量，其中Ti，Ri，Si（i=n-1）表示調(diào)查當(dāng)天前1～5 d,1～10 d,1～15 d,1～20 d 和1～30 d 的平均氣溫、累計(jì)雨量、累計(jì)日照時(shí)數(shù)。利用PASW Statistics 18 軟件，進(jìn)行影響因子的相關(guān)性分析。

2 結(jié)果與分析

2.1 關(guān)鍵氣象因子篩選

根據(jù)煙草馬鈴薯Y 病毒病發(fā)病規(guī)律與預(yù)測(cè)需求，對(duì)2009—2011年每年5—8 月各時(shí)段的光熱水三氣象因子與煙草馬鈴薯Y 病毒病發(fā)病株率、病情指數(shù)間的關(guān)系進(jìn)行相關(guān)性分析，結(jié)果見表2。煙草馬鈴薯Y 病毒病發(fā)病株率、病情指數(shù)與平均氣溫、累計(jì)降水量均呈正相關(guān)，與累計(jì)日照時(shí)數(shù)均呈負(fù)相關(guān)。表明氣溫高、降水多、日照少，容易引發(fā)煙草馬鈴薯Y 病毒病的發(fā)生。通過F 檢驗(yàn)結(jié)果表明，與發(fā)病株率和病情指數(shù)達(dá)到極顯著水平的氣象因子分別有10 個(gè)和7 個(gè)。

表2 煙草馬鈴薯Y 病毒病發(fā)病株率及病情指數(shù)與氣象因子的相關(guān)性① （n=360）

2.2 模型的建立

利用多元回歸分析方法，選取各時(shí)段煙草馬鈴薯Y病毒病發(fā)病株率、病情指數(shù)與光熱水三要素復(fù)相關(guān)系數(shù)較大的因子建立模型，即選取通過F 檢驗(yàn)達(dá)到極顯著的因子組合建立回歸模型，得到煙草馬鈴薯Y 病毒病與三因子權(quán)重組合回歸模型，見表3 和表4。

2.3 模型檢驗(yàn)

選取復(fù)相關(guān)系數(shù)0.60 以上，同時(shí)達(dá)到極顯著的模型進(jìn)行檢驗(yàn)，如沒有達(dá)到極顯著的回歸方程，則選取相關(guān)性較好的回歸方程進(jìn)行模擬檢驗(yàn)，并根據(jù)模擬檢驗(yàn)結(jié)果與實(shí)測(cè)值進(jìn)行馬鈴薯Y 病毒病發(fā)病株率、病情指數(shù)與氣象要素的相關(guān)模擬分析。

（1）發(fā)病株率與光熱水因子的模擬分析。圖1 表明，2009—2011年馬鈴薯Y 病毒病發(fā)病株率實(shí)測(cè)值與調(diào)查前10，15，20 和30 d 的氣象因子模擬趨勢(shì)基本一致；其復(fù)相關(guān)系數(shù)分別為0.66，0.71，0.74 和0.78；標(biāo)準(zhǔn)誤差分別為2.85，2.68，2.54 和2.35，最大誤差4.4～7.4，出現(xiàn)在2009年7 月20 日—9 月10 日。

表3 煙草馬鈴薯Y 病毒病發(fā)病株率與三氣象因子權(quán)重組合回歸模型

表4 煙草馬鈴薯Y 病毒病病情指數(shù)與三氣象因子權(quán)重組合回歸模型

圖1 病株率與三氣象因子的模擬

分析顯示：馬鈴薯Y 病毒病發(fā)病初期在6 月上旬,其后發(fā)病率逐漸上升，至6 月下旬末期開始進(jìn)入病情發(fā)展期，7 月進(jìn)入發(fā)病盛期。其中2009年病情發(fā)展期至盛期持續(xù)時(shí)間57 d，時(shí)段在7 月5 日—8 月30 日；2010年持續(xù)時(shí)間52 d，時(shí)段在6 月15 日—8 月15 日；2011年持續(xù)時(shí)間僅27 d，時(shí)段在7 月15 日—8 月10 日，該年份煙草馬鈴薯Y 病毒病發(fā)病相對(duì)偏輕，其原因主要是該年份昭陽(yáng)區(qū)降水量很少，出現(xiàn)歷史上最重的春夏秋連旱現(xiàn)象所致。

（2）病情指數(shù)與光熱水因子的模擬分析。圖2 表明，除2010年7 月15 日降水、日照異常偏多，模擬值異常外，馬鈴薯Y 病毒病病情指數(shù)實(shí)測(cè)值與前15，20 和30 d 的光熱水三因素模擬趨勢(shì)基本一致，其中2009年病情指數(shù)實(shí)測(cè)值與模擬值趨勢(shì)較為一致，2009年病情指數(shù)實(shí)測(cè)值較大而模擬值也較大；2010年和2011年因連續(xù)兩年干旱降水量少，病情指數(shù)較小，模擬值相對(duì)實(shí)測(cè)值偏高；復(fù)相關(guān)系數(shù)分別為0.64，0.70 和0.72；標(biāo)準(zhǔn)誤差分別為1.90，1.76 和1.72；最大誤差為5.0，出現(xiàn)在2009年7 月15 日。

圖2 病情指數(shù)與三氣象因子的模擬

3 小結(jié)

（1）煙草馬鈴薯Y 病毒病發(fā)生初期、盛期及其發(fā)生病株率、病情指數(shù)與氣象因子中平均氣溫、累計(jì)降水量、累計(jì)日照時(shí)數(shù)的關(guān)系極為密切，一般情況下煙草馬鈴薯Y 病毒病發(fā)病株率、病情指數(shù)與平均氣溫、累計(jì)降水量呈正相關(guān)，與日照呈負(fù)相關(guān)。在氣溫高、降水多、日照少的氣候環(huán)境條件下容易導(dǎo)致煙草馬鈴薯Y 病毒病發(fā)病,烤煙團(tuán)棵期為煙草馬鈴薯Y 病毒病發(fā)生初期到發(fā)展期，烤煙旺長(zhǎng)期到中部煙葉成熟期為煙草馬鈴薯Y病毒病發(fā)展期至發(fā)病盛期，以后病情穩(wěn)定直到頂葉成熟。

（2）本試驗(yàn)中煙草馬鈴薯Y 病毒病與氣象因子的關(guān)系符合正態(tài)分布，復(fù)相關(guān)系數(shù)達(dá)0.64～0.78，且各因子均通過0.01 或0.05 的水平顯著性檢驗(yàn)；模型擬合檢驗(yàn),平均誤差較小。

（3）煙草馬鈴薯Y 病毒病與氣象因子回歸方程模型的建立需要較長(zhǎng)時(shí)間的數(shù)據(jù)資料,由于技術(shù)和分析手段上的局限性，回歸方程模擬模型有待進(jìn)一步優(yōu)化和逐步完善。

（4）雖然氣象因子與煙草馬鈴薯Y 病毒病有較好相關(guān)性，但為了進(jìn)一步提高煙草馬鈴薯Y 病毒病預(yù)測(cè)模型的準(zhǔn)確率，需加強(qiáng)煙草病害影響因子的研究，同時(shí)還應(yīng)考慮越冬病害基數(shù)、抗病性、病源、傳毒媒介及栽培管理措施等因素的影響。

［1］周本國(guó)，高正良，劉小平.安徽省煙草馬鈴薯Y 病毒?。≒VY）的發(fā)生趨勢(shì)及防治措施［J］.安徽農(nóng)學(xué)通報(bào)，1998,4(1)：34-35.

［2］余春英，張西仲，王定福，等.黔南煙區(qū)煙草馬鈴薯Y 病毒病發(fā)病原因及防治措施［J］.安徽農(nóng)業(yè)科學(xué)，2010，38(10)：5110-5112.

［3］孫宏宇，孫劍萍，欒雙，等.煙草馬鈴薯Y 病毒病流行規(guī)律的初步研究［J］.煙草科技，2005（12）：39-42.

［4］陳懷亮，張宏，李有.農(nóng)作物病蟲害發(fā)生發(fā)展氣象條件與預(yù)報(bào)方法研究綜述［J］.中國(guó)農(nóng)業(yè)氣象，2007，28(2)：212-216.

［5］淡建兵，孔德胤，劉雙平，等.河套取向日葵核菌病發(fā)生程度預(yù)測(cè)預(yù)報(bào)［J］.中國(guó)農(nóng)業(yè)氣象，2012，33(1)：142-147.

［6］賈金明，張明捷，徐巧真，等.黃河中下游小麥條銹病氣象模型預(yù)測(cè)［J］.中國(guó)農(nóng)業(yè)氣象，2008，29（3）：365-370.

［7］賈金明.黃河小麥赤霉病氣象指數(shù)的建立與應(yīng)用［J］.氣象，2002，28（3）：50-53.

［8］GB/T 23222—2008 煙草病蟲害分級(jí)及調(diào)查方法［S］.