田開慧陳怡帆 周宏威

(1.湘西自治州林業(yè)局森防站,湖南 湘西土家苗族自治州 416300;2.國(guó)家林業(yè)和草原局生物災(zāi)害防控中心/林草有害生物監(jiān)測(cè)預(yù)警國(guó)家林業(yè)和草原局重點(diǎn)實(shí)驗(yàn)室,沈陽(yáng) 110034;3.東北林業(yè)大學(xué) 機(jī)電工程學(xué)院,哈爾濱 150040)

0 引言

馬尾松毛蟲(Dendrolimuspunctatus)屬鱗翅目,枯葉蛾科[1]。分布于我國(guó)秦嶺至淮河以南各省,主要危害馬尾松(Pinusmassoniana)、濕地松(PinusselliottiiEngelm)、火炬松(Pinustaeda)和黑松(Pinusthunbergii)等松屬植物,松樹受害輕則松針枯黃卷曲,重則大片枯死。2021年馬尾松毛蟲在我國(guó)發(fā)生面積35.63×104hm2,在華北北部、長(zhǎng)江中下游局部地區(qū)偏重發(fā)生[2]。馬尾松毛蟲造成危害后,林木易受到其他病蟲害的入侵,在湖南湘西州松材線蟲病發(fā)生面積僅次于馬尾松毛蟲,最易引起暴發(fā)危害。松材線蟲病主要是由松材線蟲(Bursaphelenchusxylophilus)所引起的病害,在我國(guó)松墨天牛(Monochamusalternatus)是其主要傳媒昆蟲,主要危害馬尾松、赤松(Pinusdensiflora)、黑松和火炬松等松屬植物,致病力強(qiáng),傳播快,寄主死亡速度快[3]。而湘西州松林面積28.38×104hm2,占全州林分總面積的31.55%;蓄積量占全州森林總蓄積量的41.30%[4]。馬尾松毛蟲和松材線蟲病一旦在湘西州大規(guī)模暴發(fā),將會(huì)給森林資源帶來巨大的損失。本研究旨在分析湘西州2種主要病蟲害發(fā)生現(xiàn)狀和未來發(fā)生情況,為其防治提供指導(dǎo)。

1 材料與方法

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

本文采用的湘西州2007—2021年林業(yè)有害生物發(fā)生面積統(tǒng)計(jì)數(shù)據(jù)來源于國(guó)家林業(yè)和草原局生物災(zāi)害防控中心林業(yè)有害生物防治信息管理系統(tǒng),見表1。

表1 2007—2021年湘西州部分病蟲害發(fā)生面積Tab.1 The occurrence area of some diseases and insect pests in Xiangxi from 2007 to 2021 hm2

1.2 發(fā)生面積預(yù)測(cè)分析方法

根據(jù)《林業(yè)有害生物監(jiān)測(cè)預(yù)報(bào)管理規(guī)范》(LY/T 2517—2015)[5]和《林業(yè)有害生物監(jiān)測(cè)預(yù)報(bào)技術(shù)規(guī)范》(LY/T 2516—2015)[6]要求,對(duì)2007—2021年間湘西州馬尾松毛蟲和松材線蟲病發(fā)生數(shù)據(jù)進(jìn)行收集、整理和篩選,由于湘西州馬尾松毛蟲和松材線蟲病發(fā)生面積數(shù)據(jù)是一系列隨時(shí)間推移而形成的數(shù)據(jù)序列,且該序列具有非線性、不平均性、多樣性和規(guī)律性等特點(diǎn),因此用非線性自回歸移動(dòng)平均模型(ARIMA)與支持向量機(jī)(SVM)組合模型近似描述這個(gè)序列[7]。在使用組合模型時(shí),需將病蟲害發(fā)生數(shù)據(jù)Yt視為由線性部分?jǐn)?shù)據(jù)Lt和非線性部分?jǐn)?shù)據(jù)Nt組成,即:

在本文中選用的ARIMA非線性模型表達(dá)公式為

式中:L是滯后算子;p為原始數(shù)據(jù)本身的滯后數(shù);q為差分的次數(shù);d為預(yù)測(cè)誤差的滯后數(shù);φi和θi是待估計(jì)的參數(shù)值;εt為隨機(jī)誤差,均值0。

在使用ARIMA非線性模型進(jìn)行預(yù)測(cè)后,得到的預(yù)測(cè)結(jié)果與原數(shù)據(jù)序列得到的殘差進(jìn)行樣本重構(gòu),組成新的非線性樣本集作為SVM樣本集用于SVM進(jìn)行預(yù)測(cè),最后將2次預(yù)測(cè)結(jié)果相加,得到最終的預(yù)測(cè)結(jié)果。原理如圖1所示。

圖1 ARIMA-SVM組合模型原理Fig.1 ARIMA-SVM model principles

2 結(jié)果與分析

2.1 湘西州林業(yè)有害生物發(fā)生面積及主要種類歷史變化

據(jù)統(tǒng)計(jì),2007—2021年間湘西州累計(jì)發(fā)生病蟲害32.60×104hm2,約占全州森林總面積的三分之一。利用“雙周期移動(dòng)平均值法”[8]對(duì)2007—2021年間湘西州有害生物發(fā)生面積進(jìn)行的統(tǒng)計(jì)分析,如圖2所示;選取連續(xù)發(fā)生年份大于10 a的馬尾松毛蟲、云南松毛蟲、松茸毒蛾、松墨天牛和油茶炭疽病等病蟲害,加上重點(diǎn)關(guān)注的松材線蟲病,這6種病蟲害的發(fā)生面積統(tǒng)計(jì)分析如圖3所示。

圖2 近15 a湘西州林業(yè)有害生物發(fā)生總面積及趨勢(shì)線圖Fig.2 The occurrence area and trend line diagram of forest pests in Xiangxi Prefecture in the past 15 years

圖3 近15 a湘西州6種主要林業(yè)有害生物發(fā)生面積統(tǒng)計(jì)圖Fig.3 Statistical graph of the occurrence area of 6 major forestry pests in Xiangxi in the past 15 years

分析結(jié)果顯示,湘西州林業(yè)有害生物發(fā)生面積呈周期性變化,一般3～5 a為1個(gè)周期;其中2010年、2015年、2018年為波峰年,且2018年達(dá)到歷史峰值56 591.83 hm2。2019年處于波谷,當(dāng)前正處于顯性上升周期,且湘西州當(dāng)前馬尾松毛蟲發(fā)生數(shù)據(jù)逐漸攀升,松材線蟲病發(fā)生逐漸加重。因此,本文對(duì)湘西州的主要林業(yè)有害生物馬尾松毛蟲、松材線蟲病進(jìn)行發(fā)生數(shù)據(jù)分析及趨勢(shì)預(yù)測(cè)。

2.2 湘西州馬尾松毛蟲歷史發(fā)生數(shù)據(jù)分析

湘西州馬尾松毛蟲在2007—2021年連續(xù)發(fā)生,2016年起波及至全州7縣1市,如圖4所示。調(diào)查數(shù)據(jù)顯示,古丈縣發(fā)生面積占比最高,平均占比26.51%;花垣縣在2016年首次發(fā)現(xiàn)馬尾松毛蟲成災(zāi),占當(dāng)年疫情發(fā)生總面積31.2%。馬尾松毛蟲周期性暴發(fā)明顯,第1個(gè)暴發(fā)周期為2008—2012年,第2個(gè)暴發(fā)周期為2013—2016年,第3個(gè)暴發(fā)周期為2017—2019年,周期縮短;峰值的增幅分別為157.4%與90.6%。該地松毛蟲暴發(fā)呈周期縮短、暴發(fā)趨勢(shì)明顯、世代重疊現(xiàn)象顯著和疫情態(tài)勢(shì)嚴(yán)峻特點(diǎn)。

圖4 湘西州各縣松毛蟲發(fā)生面積統(tǒng)計(jì)圖Fig.4 Statistical map of the occurrence area of Dendrolimus punctatus in counties of Xiangxi Prefecture

2.3 湘西州松材線蟲病歷史發(fā)生數(shù)據(jù)分析

2017年,湘西州首次在永順、保靖2縣發(fā)現(xiàn)松材線蟲病災(zāi)情,發(fā)生面積186.68 hm2,2018年疫情擴(kuò)散至6個(gè)縣,發(fā)生面積合計(jì)4 473.56 hm2,發(fā)生面積增長(zhǎng)近24倍,疫情從湘西州東北方向入侵,向西南方向擴(kuò)散蔓延,如圖5所示。2017年松毛蟲在湘西州暴發(fā),削弱樹勢(shì),弱勢(shì)的松樹易致松墨天牛侵入并增殖,加快松樹死亡;次年大量天牛從“松毛蟲危害松木”上羽化,侵染健康松林,導(dǎo)致松材線蟲病進(jìn)一步擴(kuò)散和流行[9]。湘西州松材線蟲病疫情在區(qū)域上高度集中,永順縣2021年疫情發(fā)生面積占全州82.8%,瀘溪縣為唯一未發(fā)現(xiàn)疫情縣。

2.4 發(fā)生面積預(yù)測(cè)結(jié)果

本文使用的非線性自回歸移動(dòng)平均模型(ARIMA)與支持向量機(jī)(SVM)組合模型在MATLAB2018a環(huán)境下運(yùn)行,通過試驗(yàn),實(shí)現(xiàn)對(duì)湘西州馬尾松毛蟲和松材線蟲病未來發(fā)生面積的預(yù)測(cè),見表2,如圖6和圖7所示。該組合模型可以在時(shí)間序列數(shù)據(jù)中捕獲一組不同標(biāo)準(zhǔn)時(shí)間結(jié)構(gòu),在擬合周期性波動(dòng)的非平穩(wěn)數(shù)據(jù)時(shí)表現(xiàn)良好,適用于湘西州馬尾松毛蟲和松材線蟲病發(fā)生面積預(yù)測(cè)。

預(yù)測(cè)結(jié)果顯示,馬尾松毛蟲在湘西州發(fā)生面積仍將延續(xù)周期性暴發(fā)趨勢(shì),預(yù)計(jì)將在2023—2025年間迎來新的暴發(fā)期,且峰值持續(xù)升高,預(yù)計(jì)最高發(fā)病面積可達(dá)50 002.50 hm2,2026年開始回落;松材線蟲病自2017年在湘西州發(fā)現(xiàn)以來,始終呈現(xiàn)較緩慢的上升趨勢(shì),但從預(yù)測(cè)結(jié)果可以看出,未來松材線蟲病增長(zhǎng)速率將大幅增高,發(fā)生面積將呈指數(shù)型上升,至2026年發(fā)生面積可達(dá)16 360.82 hm2,年同比增長(zhǎng)32.21%。

2.5 發(fā)生面積預(yù)測(cè)結(jié)果分析

預(yù)測(cè)結(jié)果顯示,馬尾松毛蟲的發(fā)生將繼續(xù)延續(xù)周期性暴發(fā)趨勢(shì),與其歷史發(fā)生規(guī)律相符,暴發(fā)周期逐漸縮短,是由湘西州自然及地理因素造成的發(fā)生態(tài)勢(shì)逐漸加劇,主要原因包括:

(1)高溫高濕,利于松毛蟲繁衍;馬尾松占比高,食物基礎(chǔ)好。湘西州冬春交替時(shí)多為陰雨寡照天氣,5月的高溫高濕天氣增多,利于松毛蟲繁衍[10]。馬尾松林面積占湘西州松林總面積的70.68%,豐富的資源給馬尾松毛蟲暴發(fā)成災(zāi)提供了食物基礎(chǔ)。

圖5 湘西州各縣松材線蟲病發(fā)生面積統(tǒng)計(jì)圖Fig.5 Statistical map of the occurrence area of pine wilt disease in counties of Xiangxi Prefecture

表2 湘西州馬尾松毛蟲、松材線蟲病發(fā)生趨勢(shì)預(yù)測(cè)結(jié)果Tab.2 Prediction results of occurrence trend of Dendrolimus punctatus and pine wilt disease in Xiangxi Prefecture

圖6 湘西州馬尾松毛蟲病發(fā)生面積預(yù)測(cè)結(jié)果Fig.6 Prediction results of D. punctatus in Xiangxi Prefecture

圖7 湘西州松材線蟲病發(fā)生面積預(yù)測(cè)結(jié)果Fig.7 Prediction results of pine wilt disease in Xiangxi Prefecture

(2)林分簡(jiǎn)單,森林抵抗力差。湘西州人工營(yíng)造的純林面積大,其中以馬尾松為主,分布范圍廣,林分結(jié)構(gòu)相對(duì)簡(jiǎn)單,抵抗林業(yè)生物災(zāi)害的能力不強(qiáng),一旦發(fā)生危害,易出現(xiàn)快速蔓延,使湘西州全域成災(zāi)[11]。

(3)世代重疊,疫情惡化。經(jīng)過近年來持續(xù)重點(diǎn)防治,該地馬尾松毛蟲蟲口密度、發(fā)生面積在持續(xù)人為干擾下呈下降態(tài)勢(shì),但發(fā)生周期的峰谷差值變小,仍處于多發(fā)態(tài)勢(shì),世代重疊現(xiàn)象增多[12]。

同樣地,預(yù)測(cè)結(jié)果顯示松材線蟲病將持續(xù)暴發(fā),且疫情逐年加重。松材線蟲病最初在永順縣發(fā)生,而后逐年由東北向西南方向擴(kuò)散,一直未得到有效的阻隔防治,仍有大片松林存在潛在風(fēng)險(xiǎn),根據(jù)預(yù)測(cè)結(jié)果若維持現(xiàn)有防治力量,將呈現(xiàn)暴發(fā)增長(zhǎng)相符。而松材線蟲病在湘西州成災(zāi)還存在人為和生境因素:

(1)疫木運(yùn)輸,導(dǎo)致疫情傳入及擴(kuò)散。近年隨著交通、電力和電信等基建項(xiàng)目在湘西州落地,物流頻繁、運(yùn)輸業(yè)興起,疫情隨疫木及其制品和包裝材料入侵及擴(kuò)散蔓延[14]。

(2)松毛蟲災(zāi)害嚴(yán)重,松木樹勢(shì)弱。湘西州松毛蟲災(zāi)害嚴(yán)重,削弱樹勢(shì);弱勢(shì)的松樹更易感染松材線蟲病[6]。

(3)除治能力不足,未能及時(shí)遏制疫情蔓延,社會(huì)對(duì)林業(yè)病蟲害危害認(rèn)知度不高、防控意識(shí)不強(qiáng);當(dāng)?shù)亓謽I(yè)部門防控資金緊張、防控技術(shù)落后。

3 結(jié)論與討論

湘西州森林占比高,尤以松林為重,警惕其主要蟲害馬尾松毛蟲、松材線蟲病的發(fā)生在生態(tài)保護(hù)及經(jīng)濟(jì)發(fā)展上都具有重大意義。湘西州歷年發(fā)生數(shù)據(jù)顯示:馬尾松毛蟲具有明顯周期性暴發(fā)趨勢(shì),峰值顯著激增,世代重疊態(tài)勢(shì)明顯;松材線蟲病從東北向西南擴(kuò)散蔓延,疫情逐年加重。究其原因,湘西州松毛蟲、松材線蟲病災(zāi)害疫情是由于氣候、林分結(jié)構(gòu)、生態(tài)條件以及人為活動(dòng)等因子共同導(dǎo)致。

ARIMA-SVM組合模型能綜合反映森林病蟲害發(fā)生面積時(shí)間序列的趨勢(shì)變化、周期性變化以及隨機(jī)性變化,經(jīng)過反復(fù)測(cè)試修正各項(xiàng)參數(shù)以確定最優(yōu)模型,對(duì)時(shí)間序列分析,以實(shí)現(xiàn)短期預(yù)測(cè)[15]。利用ARIMA-SVM組合模型的預(yù)測(cè)結(jié)果顯示,若保持現(xiàn)有防治能力不變,馬尾松毛蟲將在2023—2025年迎來又一暴發(fā)期;松材線蟲病發(fā)生面積將呈指數(shù)型上升。湘西州應(yīng)該加強(qiáng)監(jiān)測(cè)、檢疫、防治等工作力度,警惕馬尾松毛蟲、松材線蟲病在湘西州暴發(fā)。但是ARIMA-SVM組合模型也有一些不足之處,只能用于短期預(yù)測(cè)[16],在實(shí)際運(yùn)用中還應(yīng)考慮受到的其他未知因素影響,需在下一步工作中進(jìn)一步探索影響因素變量,修正模型,以期更準(zhǔn)確預(yù)測(cè)森林病蟲害發(fā)生趨勢(shì)。

通過分析 2007—2021 年湘西州林業(yè)有害生物發(fā)生面積數(shù)據(jù),對(duì)湘西州馬尾松毛蟲和松材線蟲病發(fā)病面積的短期預(yù)測(cè),能為林業(yè)有害生物防治提供科學(xué)依據(jù)。因此,對(duì)湘西縣的林業(yè)有害生物防治提出幾點(diǎn)建議。

(1)應(yīng)提升監(jiān)測(cè)預(yù)報(bào)預(yù)警水平,以湘西州的3個(gè)國(guó)家級(jí)中心測(cè)報(bào)點(diǎn)為核心,帶動(dòng)全州林業(yè)有害生物監(jiān)測(cè)工作規(guī)范化、信息化、科學(xué)化,實(shí)時(shí)監(jiān)測(cè)災(zāi)情動(dòng)態(tài)變化[17-19],特別是加強(qiáng)對(duì)湘西州唯一未發(fā)生松材線蟲病的瀘溪縣的疫情監(jiān)測(cè),警惕花垣縣新發(fā)馬尾松毛蟲成災(zāi)事故重演。

(2)要強(qiáng)化檢疫工作,加大檢疫執(zhí)法力度,嚴(yán)防疫木運(yùn)輸造成的疫情擴(kuò)散。

(3)要及時(shí)做好防治工作,對(duì)已發(fā)生馬尾松毛蟲、松材線蟲病的疫區(qū)做好制定防治預(yù)案,做好防治藥劑藥械準(zhǔn)備及防治資金籌措[20]。對(duì)僅發(fā)生馬尾松毛蟲的疫區(qū),抓住防控關(guān)鍵時(shí)間節(jié)點(diǎn),做好松材線蟲寄主松墨天牛的防治工作。