閔凡祥，王貴江，楊 帥，胡林雙，王文重，魏 琪，李文楓

（1黑龍江省農(nóng)業(yè)科學(xué)院馬鈴薯研究所，哈爾濱 150086；2黑龍江省農(nóng)業(yè)科學(xué)院，哈爾濱 150086；3黑龍江省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)遙感與信息研究所，哈爾濱 150086）

0 引言

馬鈴薯營(yíng)養(yǎng)豐富，糧菜飼兼用，增產(chǎn)增收潛力大，對(duì)保障國(guó)家糧食安全有著重要的意義[1]。數(shù)據(jù)統(tǒng)計(jì)顯示，黑龍江省馬鈴薯平均單產(chǎn)僅為16.5 t/hm2，僅為發(fā)達(dá)國(guó)家單產(chǎn)的1/3[2]，導(dǎo)致單產(chǎn)水平較低主要原因之一就是馬鈴薯晚疫病爆發(fā)頻繁。馬鈴薯晚疫病是由致病疫霉菌引起，導(dǎo)致馬鈴薯莖葉死亡和塊莖腐爛的一種毀滅性卵菌病害[3]，雖然，化學(xué)藥劑防治可以有效控制馬鈴薯晚疫病發(fā)生與流行，但化學(xué)防治技術(shù)普遍存在用藥不及時(shí)、過(guò)度施藥和盲目施藥等問(wèn)題，造成防治效果不理想、病菌抗藥性增強(qiáng)、次生有害生物再猖獗和環(huán)境污染等一系列嚴(yán)重后果，嚴(yán)重影響馬鈴薯單產(chǎn)水平提高，已成為馬鈴薯產(chǎn)業(yè)發(fā)展重要制約瓶頸之一[4]。

為了提升馬鈴薯晚疫病防治效率，避免在防治過(guò)程中用藥不及時(shí)和盲目過(guò)度施藥，達(dá)到精準(zhǔn)防治的目的，世界各國(guó)均在馬鈴薯晚疫病監(jiān)測(cè)預(yù)警技術(shù)方面開(kāi)展大量工作，目前國(guó)際上已報(bào)道的預(yù)測(cè)預(yù)報(bào)模型有17種之多，均是根據(jù)當(dāng)?shù)貧夂蛱攸c(diǎn)及馬鈴薯晚疫病發(fā)生規(guī)律等相關(guān)參數(shù)建立模型，并已在各地區(qū)廣泛應(yīng)用，如在丹麥NEGFRY模型、美國(guó)BLITECAST模型和比利時(shí)CARAH模型等[5]。中國(guó)在馬鈴薯晚疫病監(jiān)測(cè)預(yù)警技術(shù)研究方面起步較晚，目前，應(yīng)用最廣泛是由謝開(kāi)云等[6]引入了比利時(shí)CARAH模型，已經(jīng)在10余省進(jìn)行應(yīng)用，在指導(dǎo)馬鈴薯晚疫病防控和促進(jìn)馬鈴薯產(chǎn)業(yè)健康發(fā)展中發(fā)揮了重要作用[7]。以往報(bào)道比利時(shí)CARAH模型應(yīng)用，僅僅是單一模型防效分析，參數(shù)校驗(yàn)，缺少系統(tǒng)比較，本試驗(yàn)利用比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型，針對(duì)黑龍江省氣象數(shù)據(jù)、馬鈴薯晚疫病發(fā)生情況和殺菌機(jī)特點(diǎn)等，通過(guò)田間試驗(yàn)，與丹麥NEGFRY模型、每7天施藥處理和每10天施藥處理方法進(jìn)行評(píng)價(jià)，比較分析不同處理的防治效果，明確CARAH監(jiān)測(cè)預(yù)警模型減藥次數(shù)和經(jīng)濟(jì)效益分析，進(jìn)而建立黑龍江省馬鈴薯晚疫病CARAH監(jiān)測(cè)預(yù)警模型體系，指導(dǎo)黑龍江省馬鈴薯晚疫病綜合防治工作。

1 材料與方法

1.1 供試材料

1.1.1 供試品種 本研究采用品種為國(guó)家種質(zhì)克山馬鈴薯種質(zhì)資源庫(kù)提供的荷蘭品種‘費(fèi)烏瑞它’，英文名Favorita，種薯級(jí)別為G3，該品種對(duì)晚疫病高度敏感。

1.1.2 藥劑 72%霜脲氰錳鋅（商品名：克露），杜邦中國(guó)集團(tuán)有限公司生產(chǎn)；50%氟啶胺（商品名：福帥得），日本石原產(chǎn)業(yè)株式會(huì)社生產(chǎn)；68.75%氟菌·霜霉威SC（商品名：銀法利），拜耳作物科學(xué)(中國(guó))有限公司生產(chǎn)；23.4%雙炔酰菌胺（商品名：瑞凡），先正達(dá)生物科技(中國(guó))有限公司提供；75%代森錳鋅，江蘇龍燈化學(xué)有限公司生產(chǎn)。

1.1.3 氣象站 本試驗(yàn)利用Vantage Pro2型美國(guó)DAVIS便攜式6要素氣象站，可采集田間的溫度、相對(duì)濕度、降雨量、光照、風(fēng)速和風(fēng)向等氣象因子[6]，儀器憑借太陽(yáng)能板供電，數(shù)據(jù)實(shí)時(shí)采集后通過(guò)GPRS進(jìn)行遠(yuǎn)程數(shù)據(jù)傳輸。

1.2 試驗(yàn)方法

1.2.1 試驗(yàn)站點(diǎn)及處理 2018年，田間試驗(yàn)設(shè)置于黑龍江現(xiàn)代農(nóng)業(yè)示范區(qū)，該區(qū)域氣候濕潤(rùn)，溫度適宜，適合晚疫病發(fā)生。本試驗(yàn)設(shè)置為5個(gè)處理，分別為比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型（簡(jiǎn)稱(chēng)CM）、每7天施藥處理（簡(jiǎn)稱(chēng)WM）、每10天施藥處理（簡(jiǎn)稱(chēng)TM）施藥和丹麥NegFry模型處理（簡(jiǎn)稱(chēng)NM）和空白對(duì)照（簡(jiǎn)稱(chēng)CK）。每處理重復(fù)4次，小區(qū)面積24 m2，完全隨機(jī)區(qū)組排列，四周設(shè)隔離行和接種行。

1.2.2 數(shù)據(jù)的采集方法 本研究待田間出苗80%后，利用便攜式6要素氣象站收集氣象數(shù)據(jù)，同時(shí)調(diào)查病害發(fā)生情況。

1.2.3 確定馬鈴薯晚疫病潛在侵染嚴(yán)重程度的方法 馬鈴薯晚疫病潛在侵染嚴(yán)重程度主要是根據(jù)比利時(shí)Guntz-Divoux基本原理進(jìn)行判斷，原理為當(dāng)每日的平均溫度和相對(duì)濕度大于90%時(shí)，持續(xù)小時(shí)數(shù)越多，馬鈴薯晚疫病潛在侵染嚴(yán)重程度就越大[8]。

1.2.4 確定馬鈴薯晚疫病侵染循環(huán)和施藥劑時(shí)間方法 馬鈴薯晚疫病侵染循環(huán)曲線是由每天平均溫度決定，每天平均溫度對(duì)應(yīng)著一定分值（見(jiàn)表1）。在晚疫病菌存在情況下，病原菌從侵染到萌發(fā)，即為晚疫病侵染循環(huán)。用曲線表示侵染循環(huán)，為每日平均溫度得分值累加為7分，對(duì)應(yīng)為病害循環(huán)曲線。如果沒(méi)有采取預(yù)防措施，則田間出現(xiàn)晚疫病初期癥狀。比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型建議，當(dāng)累加得分值在5～6分時(shí)，是最佳施藥時(shí)間。下一次噴藥時(shí)間確定是根據(jù)殺菌劑持效期、新的侵染循環(huán)、降雨和品種抗性判斷。

表1 病害侵染循環(huán)開(kāi)始后每日分值的計(jì)算標(biāo)準(zhǔn)

1.2.5 晚疫病調(diào)查及數(shù)據(jù)處理方法 每次噴施藥劑前，進(jìn)行晚疫病田間調(diào)查。發(fā)病后，需1周調(diào)查2次，采用James分級(jí)計(jì)數(shù)法[9]，計(jì)算病害發(fā)病百分率。利用病害發(fā)病百分率和調(diào)查時(shí)間，得出AUDPC(Area under disease progress curve)值。AUDPC值越大，則表示病害發(fā)病越嚴(yán)重，病害的防治效果就越差，反之，防治效果越顯著。對(duì)不同處理AUDPC值進(jìn)行平均值多重比較，采用LSD法進(jìn)行顯著性分析[10]。

1.2.6 田間產(chǎn)量和晚疫病病薯率調(diào)查方法 對(duì)小區(qū)試驗(yàn)進(jìn)行測(cè)產(chǎn)，主要調(diào)查植株數(shù)、大中小薯比例、晚疫病爛薯率和單位產(chǎn)量[11]。計(jì)算單株產(chǎn)量及公頃產(chǎn)量和晚疫病病薯率[12]。

2 結(jié)果與分析

2.1 不同晚疫病監(jiān)測(cè)預(yù)警模型的結(jié)果分析

2018年5月9日試驗(yàn)地播種，出苗齊全時(shí)間為6月2日。所以，本年CARAH監(jiān)測(cè)預(yù)警模型起止時(shí)間從6月2日開(kāi)始，截止時(shí)間為9月30日，生育期結(jié)束，整個(gè)馬鈴薯生育期溫度曲線（藍(lán)色曲線）、相對(duì)濕度曲線（綠色曲線）、降雨（藍(lán)色柱狀）、病害嚴(yán)重度（黃色柱形）和病程曲線（紅色線）情況見(jiàn)圖2，由圖2可知，本試驗(yàn)區(qū)晚疫病潛在發(fā)生時(shí)間段共計(jì)44次，其中極重侵染為20次、重度侵染有10次、病害中度侵染有11次、病害輕度侵染為3次，圖2中黃色柱形高矮寬度表示為潛在發(fā)生晚疫病嚴(yán)重程度，柱形越高，寬度越大，表明晚疫病發(fā)生越嚴(yán)重。

圖1 2018年馬鈴薯晚疫病潛在發(fā)生嚴(yán)重程度及侵染循環(huán)情況

圖2 NegFry測(cè)報(bào)負(fù)值分析模型對(duì)首次施藥分析

2018年，利用丹麥NegFry監(jiān)測(cè)預(yù)警模型分析田間氣象數(shù)據(jù)，結(jié)果顯示從6月16日起，至6月25日，每日風(fēng)險(xiǎn)值累加達(dá)到143.71，6月25日的日風(fēng)險(xiǎn)值為16，所以，6月25日為丹麥NegFry監(jiān)測(cè)預(yù)警模型第一次噴施殺菌劑時(shí)間（見(jiàn)圖3）。第二次施藥時(shí)間，需通過(guò)每日風(fēng)險(xiǎn)值進(jìn)行累加，每日風(fēng)險(xiǎn)值累加大于等于40時(shí)，需噴施殺菌劑。之后，將每日風(fēng)險(xiǎn)值歸0，重新累加，進(jìn)而確定下一次施藥時(shí)間[13]。

圖3 丹麥NegFry監(jiān)測(cè)預(yù)警模型施藥次數(shù)情況

2.2 不同處理施藥情況結(jié)果分析

CM處理時(shí)間依次為6月28日、7月10日、7月24日、8月5日、8月13日、8月29日和9月7日，共計(jì)噴藥7次，每次施藥依次為75%代森錳鋅、68.75%氟菌·霜霉威SC、72%霜脲氰錳鋅、23.4%雙炔酰菌胺、68.75%氟菌·霜霉威SC、72%霜脲氰錳鋅、50%氟啶胺（見(jiàn)表2）；NM處理噴藥時(shí)間依次分別為6月25日、7月5日、7月16日、7月21日、8月4日、8月14日、8月24日、8月30日和9月4日（圖3），噴施藥劑順序依次為75%代森錳鋅、68.75%氟菌·霜霉威SC、72%霜脲氰錳鋅、23.4%雙炔酰菌胺、68.75%氟菌·霜霉威SC、72%霜脲氰錳鋅、23.4%雙炔酰菌胺、75%代森錳鋅、50%氟啶胺；WM和TM處理，首次噴藥時(shí)間選擇丹麥NegFry模型時(shí)間，為6月25日，WM處理共計(jì)噴藥10次，TM處理噴藥8次，具體噴藥時(shí)間見(jiàn)表2。通過(guò)比較，WM處理、TM處理和NM處理和與CM處理相比，減少施藥次數(shù)分別為3次、1次和2次。

表2 不同處理噴施藥劑時(shí)間及施藥類(lèi)型情況

2.3 CARAH監(jiān)測(cè)預(yù)警模型與病害實(shí)際發(fā)生擬合情況分析

2018年，采用CARAH監(jiān)測(cè)預(yù)警模型，分析出馬鈴薯整個(gè)生育期內(nèi)出現(xiàn)44次病害病程曲線（見(jiàn)圖2）。實(shí)際調(diào)查顯示馬鈴薯晚疫病各個(gè)級(jí)別發(fā)生次數(shù)為41次，CARAH監(jiān)測(cè)預(yù)警模型分析晚疫病潛在發(fā)生時(shí)間段共計(jì)44次，與田間實(shí)際病害發(fā)生數(shù)擬合度達(dá)到93%，具體情況如下：6月26—7月14日，田間初始出現(xiàn)發(fā)病中心病株，模型與實(shí)際病害發(fā)生情況擬合度為67%。7月14—8月5日，田間病害進(jìn)入指數(shù)期，葉片發(fā)病嚴(yán)重程度達(dá)中度侵染，模型與實(shí)際病害發(fā)生情況擬合度達(dá)91%。8月5日—8月31日，病害開(kāi)始大流行，模型與實(shí)際病害發(fā)生情況擬合度達(dá)100%。9月1日—9月10日，晚疫病發(fā)病達(dá)高峰，之后，病害進(jìn)入衰退期。通過(guò)調(diào)查表明，CARAH監(jiān)測(cè)預(yù)警模型的結(jié)果與實(shí)際發(fā)生期高度接近（見(jiàn)表3），測(cè)報(bào)準(zhǔn)確率達(dá)到93%。

表3 CARAH監(jiān)測(cè)預(yù)警模型的驗(yàn)證結(jié)果

2.4 不同處理防治效果比較

通過(guò)田間病害發(fā)生調(diào)查，結(jié)果表明：6月30日空白對(duì)照(CK)首先發(fā)現(xiàn)病斑，僅零星葉片發(fā)病，其他試驗(yàn)處理均沒(méi)有發(fā)現(xiàn)病斑。7月14日CM處理、WM處理和NM處理陸續(xù)發(fā)病，計(jì)入病害指數(shù)增長(zhǎng)期。截至生育期結(jié)束，各處理TM、CM、NM和WM病害發(fā)病百分率分別為37.50%、26.25%、22.25%和22.00%，9月1日CK發(fā)病率為100%，全部死亡（圖4）。計(jì)算不同處理AUDPC值，CK、TM、CM、NM和WM分別為12.25、1.85、1.27、1.11和1.08（圖5）。對(duì)上述AUDPC值進(jìn)行數(shù)據(jù)差異顯著性分析，分析顯示，各處理與CK比較達(dá)顯著水平，WM、CM和NM處理之間比較差異并不顯著。CM處理與WM和NM處理達(dá)到同等效果。

圖4 不同處理的晚疫病病害發(fā)展曲線

圖5 不同處理AUDPC值

2.5 不同處理對(duì)田間馬鈴薯產(chǎn)量及商品薯率的影響

不同處理由于防治晚疫病的效果不同，導(dǎo)致產(chǎn)量差異各不相同。產(chǎn)量最高是CARAH監(jiān)測(cè)預(yù)警模型處理(CM)，其次是WM、NM、TM和CK，分別是43.24、42.61、42.27、39.74和29.12t/hm2。采用LSD法對(duì)數(shù)據(jù)進(jìn)行平均值多重比較，結(jié)果顯示W(wǎng)M、CM和NM與CK處理相比產(chǎn)量差異達(dá)0.05顯著水平，增產(chǎn)效果較明顯。所以，利用CARAH監(jiān)測(cè)預(yù)警模型可以有效控制馬鈴薯晚疫病的流行，減少因此病害造成損失，增產(chǎn)效果較顯著，CARAH監(jiān)測(cè)預(yù)警模型與空白對(duì)照相比增產(chǎn)高達(dá)48.47%。對(duì)不同處理馬鈴薯晚疫病病薯率調(diào)查發(fā)現(xiàn)，病薯率從高到低依次為CK、NM、TM、CM和WM（表4），比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型處理與對(duì)照相比病薯率降低4.43個(gè)百分點(diǎn)，所以，CARAH監(jiān)測(cè)預(yù)警模型處理能夠顯著降低馬鈴薯塊莖感染病原菌的概率，從而減少次年初侵染源。

表4 產(chǎn)量和病薯率情況及差異分析

商品薯率高低是由大中薯率多少所決定，因此，本研究對(duì)不同處理大中小薯率進(jìn)行統(tǒng)計(jì)，明確不同處理對(duì)商品薯率影響。研究結(jié)果表明：病害防控效果越好，相應(yīng)小薯率就越低，對(duì)應(yīng)商品薯率就越高（見(jiàn)圖6）。其中，商品薯率最高是WM，為95.06%，其次為NM為94.14%，CM為93.98%，TM為89.64%，CK為81.35%。與CK相比，采用CARAH監(jiān)測(cè)預(yù)警模型處理提高大中薯率12.63個(gè)百分點(diǎn)，小薯率降低12.63個(gè)百分點(diǎn)，所以，本實(shí)驗(yàn)CARAH監(jiān)測(cè)預(yù)警模型、每7天噴試藥劑處理和丹麥NegFry監(jiān)測(cè)預(yù)警模型處理可有效提高馬鈴薯商品薯率。

圖6 不同處理大中小薯率

2.6 經(jīng)濟(jì)效益分析

調(diào)查結(jié)果表明（表5），比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型處理(CM)與空白對(duì)照處理(CK)相比較，CM處理新增產(chǎn)值為14116元/hm2，不計(jì)算人工費(fèi)用，扣除用藥成本3150元/hm2，新增純收入為10966元/hm2；每7天噴藥處理(WM)與空白對(duì)照(CM)相比，WM處理新增產(chǎn)值為13491元/hm2，人工費(fèi)用除外，去除藥劑成本4500元/hm2，WM處理新增純收入8991元/hm2?？瞻讓?duì)照處理(CM)與丹麥NEGFER測(cè)報(bào)模型處理(NM)比較，NM處理新增產(chǎn)值為13147元/hm2，扣除用藥成本4050元/hm2，NM處理新增純收入9097元/hm2?？瞻讓?duì)照處理(CK)與每10天噴藥處理(TM)相比，TM處理新增產(chǎn)值為10615元/hm2，減去藥劑成本3600元/hm2，TM處理新增純收入7015元/hm2。由此得出，采用CM處理較CK增收10966元/hm2，用藥成本比WM處理節(jié)約1350元/hm2，新增純收入比WM處理增加1975元/hm2，較NM處理增收1869元/hm2。所以，利用CARAH監(jiān)測(cè)預(yù)警模型防治馬鈴薯晚疫病，即可增加經(jīng)濟(jì)效益，又能夠降低農(nóng)藥使用量，達(dá)到增效和減藥的目的。

表5 不同馬鈴薯晚疫病防治處理的經(jīng)濟(jì)效益評(píng)價(jià) 元/hm2

3 討論與結(jié)論

比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型應(yīng)用在中國(guó)擁有20余年歷史，已經(jīng)在10余省進(jìn)行應(yīng)用[14-18]，但對(duì)于模型的評(píng)價(jià)，僅僅是通過(guò)病情指數(shù)調(diào)查，對(duì)模型防治效果進(jìn)行分析，存在著對(duì)模型原理和應(yīng)用技術(shù)掌握不夠而理解出現(xiàn)偏差的情況[7]。本研究以丹麥NEGFRY模型和常規(guī)施藥防治措施作為對(duì)照，引入病程曲線下面積即AUDPC值，明確CARAH監(jiān)測(cè)預(yù)警模型在病害發(fā)展過(guò)程中防治效果。有效解決數(shù)據(jù)分析偏差和缺乏有效對(duì)照的問(wèn)題，提高黑龍江省馬鈴薯晚疫病監(jiān)測(cè)預(yù)警水平。

比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型的準(zhǔn)確性主要受品種抗性、病菌種群結(jié)構(gòu)、區(qū)域環(huán)境和殺菌劑特性等因素會(huì)影響[7]，本試驗(yàn)為避免品種影響，選用感病品種‘費(fèi)烏瑞它’，但在調(diào)查中發(fā)現(xiàn)，相同品種在不同氣候環(huán)境、不同菌源結(jié)構(gòu)下，不同藥劑防治下，晚疫病的流行規(guī)律仍然存在一定的差異。所以，為了保證試驗(yàn)一致性，試驗(yàn)過(guò)程中采用自然發(fā)病和人工接種同時(shí)進(jìn)行方法，接種病原菌為優(yōu)勢(shì)菌株。在調(diào)查數(shù)據(jù)處理上采用AUDPC值過(guò)程量的分析方法，避免單次偏差問(wèn)題出現(xiàn)，以保證整體評(píng)估的準(zhǔn)確性。

比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型原理是根據(jù)感病品種和病原菌特性，即溫度、濕度和降雨情況，分析出可能潛在病害發(fā)生時(shí)間段，利用殺菌劑，無(wú)死角全方位覆蓋防治。但是，“病害三角”中，寄主抗性一樣非常重要，比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型只是針對(duì)感病品種，而對(duì)于耐病品種，如果采用CARAH監(jiān)測(cè)預(yù)警模型防治，勢(shì)必產(chǎn)生不必要的浪費(fèi)。因此，首先需要開(kāi)展對(duì)區(qū)域主要馬鈴薯栽培品種抗性級(jí)別鑒定，明確抗性級(jí)別，對(duì)于感病品種，直接可以采用測(cè)報(bào)模型進(jìn)行防治。對(duì)于中感或中抗品種，需開(kāi)展病原菌在與不同品種侵染曲線研究，即不同品種病程曲線調(diào)查試驗(yàn)，明確不同耐病品種病程發(fā)展曲線，進(jìn)而完善比利時(shí)CARAH監(jiān)測(cè)預(yù)警模型，全方位為黑龍江省馬鈴薯晚疫病防治服務(wù)。

本研究結(jié)果表明CARAH監(jiān)測(cè)預(yù)警模型可用于黑龍江省馬鈴薯晚疫病綜合防治工作中，與高頻繁施藥處理和國(guó)際先進(jìn)的測(cè)報(bào)模型處理相比較，取得同等防治效果，而且CARAH監(jiān)測(cè)預(yù)警模型與丹麥NEGFRY模型處理，每7天和10天施藥處理相比，減少殺菌劑使用次數(shù)1～3次。與CK相比，增產(chǎn)48.47%小薯率降低12.63個(gè)百分點(diǎn)，病薯率降低4.43個(gè)百分點(diǎn)，增加純收入10 966元/hm2。所以，CARAH監(jiān)測(cè)預(yù)警模型可用于指導(dǎo)黑龍江省馬鈴薯晚疫病防治工作，進(jìn)而減少不必要?dú)⒕鷦﹪娛?，達(dá)到精準(zhǔn)治的效果，完成減藥目標(biāo)。