李 濤,李澤華,余仲東,李 斌
(1四川省森林病蟲防治檢疫總站,四川 成都 610082;2西北農林科技大學 林學院,陜西 楊凌 712100)
美國白蛾(Hyphantriacunea),又稱秋幕毛蟲、網幕毛蟲,是一種世界性檢疫害蟲。它原產北美洲,1922年在加拿大最早發(fā)現,此后相繼在美國、墨西哥局部地區(qū)發(fā)生。第二次世界大戰(zhàn)期間(1940年之后),隨著地區(qū)間的物資運輸從美洲擴散到歐洲的許多國家,先后入侵匈牙利、羅馬尼亞、保加利亞、奧地利、原蘇聯(lián)、波蘭和法國等國家。1945年傳入日本,1957年傳入韓國,1961年傳入朝鮮,1979年傳入中國。最早在我國與朝鮮交界處的遼寧省丹東地區(qū)發(fā)現,之后陸續(xù)傳播擴散到山東(1981年)、陜西(1984年)、河北(1989年)、上海(1994年)、天津(1995年)、北京(2003年)、河南(2008年)、吉林(2010年)、江蘇(2011年)、安徽(2012年)等省(市)[1]。李淑賢等[2]預測,美國白蛾可分布于我國21.20°-46.33°N,97.80°-132.11°E和36.81°-41.85°N,76.00°-94.66°E的廣大地區(qū)。
美國白蛾食性非常雜,寄主植物幾乎包括疫情區(qū)所有林木、果樹、蔬菜、農作物、園林花卉以及雜草等,能夠危害175~300種農林植物及野生植物[3-4];適應性強,不取食15 d后仍然可以正常繁殖;繁殖量大,成蟲產卵量高達1 800粒[5];傳播能力強,每年可自然擴散35~40 km,更可隨貨物、交通工具等進行遠距離傳播;缺乏原產地天敵的有效制約,常在新傳入地暴發(fā)成災,造成的經濟損失重大[6]。我國政府高度重視美國白蛾的防控工作,自1979年發(fā)現以來,林業(yè)部或國家林業(yè)局分別于1984年、1996年、2004年、2013年連續(xù)4次將其列入全國森林植物檢疫對象或全國林業(yè)檢疫性有害生物名單,并作為重大林業(yè)有害生物在全國范圍內組織開展了防治工作[1]。1998年,國家林業(yè)局啟動了美國白蛾國家級治理工程,工程范圍覆蓋北京、天津、河北3省(市)的67個縣(區(qū)),工程治理取得明顯成效。但近年來,美國白蛾出現了新的危害苗頭,新疫點不斷出現,發(fā)生面積明顯反彈,危害程度不斷加劇,擴散速度明顯加快,因此做好美國白蛾在中國的分布預警和危險性評估,對于美國白蛾的防治工作具有重要意義。本研究基于MaxEnt模型,結合我國105個美國白蛾已發(fā)生區(qū)環(huán)境變量,按1∶400萬尺度勾勒出四川省行政區(qū)劃面要素,并計算美國白蛾在四川省的分布范圍和分布面積,在此基礎上,根據傳統(tǒng)風險評價指標體系,評估了美國白蛾入侵四川省的風險,為盡早建立美國白蛾防控體系和應急預案提供依據。
MaxEnt軟件是從其官方網站(http://www.cs. princeton.edu/~schapire/maxent/)下載的。與MaxEent配合使用的GIS軟件是美國Esri公司ArcGis軟件。
氣候環(huán)境變量數據下載于WorldClim(http ://www.worldcfim.org/),包括1950-2000年的溫度、降水和海拔等20個變量因子。利用ArcGis提取中國行政區(qū)內環(huán)境變量并轉化為MaxEnt模型需要的ASCII格式數據,根據計算機的模型運行能力和物種分布數據本身的準確度,選擇WorldClim環(huán)境數據的空間分辨率為2.5 arc-min。根據MaxEnt模型自身附帶的Jackknife刀切法對初始變量的重要性進行權重分析。
按照MaxEnt模型要求格式,整理所收集美國白蛾在我國的105個已知分布地點(來自當前文獻資料及相關研究調查網站,包含了國內所有美國白蛾發(fā)生過的省區(qū), 參考文獻[7]和國家林業(yè)局2016年8號公告)。所得分布按照物種、經度、緯度(經度、緯度來自Google earth)的順序在Excel中列出,并另存為MaxEnt可以識別的csv格式。
根據已知分布區(qū),選擇MaxEnt模型中maxentResults.csv文件最適合已知分布區(qū)的閾值規(guī)則進行驗證,并將預測結果ASC格式通過ArcMap軟件進行轉化,采用ROC曲線分析法[8],計算AUC值,對預測模型的準確性進行檢驗[9-10]。
用MaxEnt軟件初步模擬出美國白蛾在中國的適生區(qū),將結果導入ArcGis,之后根據比例尺為1∶400萬的中國地圖做底圖進行校正,并單獨勾勒出四川省的分布,利用ArcGis軟件中的拓撲功能生成面要素,求出四川省適生區(qū)面積。中國行政區(qū)劃數據來自國家基礎地理信息網(http://nfgis.nsdigov.cn/),下載獲得比例尺為1∶400萬的中國國界及省界圖。
參照蔣青等[11]的評估體系,按目標層、準則層、指標層、評判指標進行分層賦值和綜合評估,全面評估美國白蛾在四川省的風險狀況。
美國白蛾適生區(qū)環(huán)境變量權重分析表明,20個環(huán)境變量因子中,alt(海拔)、bio1(年平均溫度)、bio5(最暖月最高溫度)、bio10(最暖季度平均溫度)、bio11(最冷季度平均溫度)、bio6(最冷月最低溫度)、bio8(最濕季度平均溫度)、bio9(最干季度平均溫度)、bio12(年均降水量) 9個環(huán)境因子為美國白蛾適生區(qū)范圍的主要影響因子,其中最冷季度平均溫度、海拔、最干季度平均溫度對美國白蛾分布范圍貢獻最大(圖1)。
利用2.1節(jié)中得出的9個環(huán)境變量因子對美國白蛾在中國的適生區(qū)預測表明,美國白蛾在我國渤海灣、膠州半島及東部沿海省區(qū)發(fā)生較為嚴重,并向西、向南逐漸深入擴散,新疆喀什、京津塘地區(qū)為發(fā)生區(qū)域(圖2),這與已報道的情況[7,12-13]基本吻合。
alt.海拔;bio1.年平均溫度;bio2.晝夜溫差月均值;bio3.等溫性;bio4.溫度季節(jié)變化標準差;bio5.最暖月最高溫度;bio6.最冷月最低溫度;bio7.年均溫度變化范圍;bio8.最濕季度平均溫度;bio9.最干季度平均溫度;bio10.最暖季度平均溫度;bio11.最冷季度平均溫度;bio12.年均降水量;bio13.最濕月降水量;bio14.最干月降水量;bio15.降水量變異系數;bio16.最濕季度降水量;bio17.最干季度降水量;bio18.最暖季度降水量;bio19.最冷季度降水量alt.Altitude;bio1.Average annual temperature;bio2.Average monthly temperature difference between night and day;bio3.Isothermal;bio4.Seasonal temperature standard deviation;bio5.Highest temperature of the warmest month;bio6.Lowest temperature of the coldest month;bio7.Average annual temperature range;bio8.Average temperature of the wettest season;bio9.Average temperature of the driest season;bio10.Average temperature of the warmest season;bio11.Average temperature of the coldest season;bio12.Average annual precipitation;bio13.Precipitation of the wettest month;bio14.Precipitation of the driest month;bio15.Precipitation variation coefficient;bio16.Precipitation of the wettest season;bio17.Precipitation of the driest season;bio18.Precipitation of the warmest season;bio19.Precipitation of the coldest season
紅色表示已有分布,值越大越適生Red color means distributed areas,and bigger value means more adaptive
利用ROC曲線分析法對美國白蛾在中國的適生區(qū)預測結果進行準確性評價,結果表明,MaxEnt模型獲得AUC平均值為0.949,明顯高于隨機概率值0.5,且AUC值處于0.9~1.0,表明預測結果準確可靠,環(huán)境變量與預測的物種地理分布模型之間相關性大(圖3)。
圖3 美國白蛾在中國適生區(qū)預測結果的AUC檢測Fig.3 AUC test of Hyphantria cunea prediction in China
通過圖像校正后,得到四川省行政區(qū)劃面要素(圖4)。利用ArcGis軟件求得四川省面積為48.3萬km2,與實際面積48.6萬km2相差0.3萬km2,相對誤差為0.6%,符合要求,并求得美國白蛾在四川省的適生區(qū)面積為11.52萬km2,占全省計算總面積的23.85%。
圖4 圖像校正后的四川省行政區(qū)劃圖Fig.4 Verified map of Sichuan
由圖5可知,除攀枝花市、甘孜藏族自治州、阿壩藏族羌族自治州、涼山彝族自治州4個市(州)外,四川省其余各市幾乎都是美國白蛾適生區(qū),其中以成都市最為適宜美國白蛾生存。
綠色深度表示美國白蛾的適生程度Deep of color represents the adaptive of Hyphantria cunea
根據美國白蛾的生物學特性、傳入可能性、危害與控制等因素進行綜合評估,各評估因子詳見表1。參照蔣青等[11]風險評估值計算公式:
P1= P11,
P3=0.4·P31+0.4·P32+0.2·P33,
P4=Max(P41,P42,P43),
P5=(P51+P52+P53)/3,
表1 美國白蛾在四川的風險評估體系Table 1 Risk assessment index for Hyphantria cunea in Sichuan
美國白蛾在中國的分布,很多學者采用不同的方法做了預測[2,13-14],均得到相似的結果。李淑賢等[2]利用GARP適生性軟件對美國白蛾在我國的適生區(qū)進行了預測,認為適生范圍分布在21.20°-46.33°N, 76°-132.11°E的區(qū)域內。本研究預測結果涵蓋了李淑賢等[2]的預測結果,同時范圍還要更廣泛,擴散趨勢更明顯,主要體現在發(fā)生的西南界上,增加了四川、海南、臺灣地區(qū),將新疆部分地區(qū)的分布等級進一步提高。這種差異可能與分析軟件不同和分析中使用的變量因子有關。本研究采用最大熵模型MaxEnt,即自然條件下有害生物會盡可能大地尋找自己的適生范圍[15-16],利用9個影響美國白蛾自然分布的環(huán)境變量因子進行了適生區(qū)預測,擴大了目前國內該蟲實際分布范圍,提高了預測的準確性,研究結果還進一步預測了美國白蛾在四川不同縣區(qū)的分布。雖然目前未在四川省內發(fā)現有美國白蛾分布,但四川省具備美國白蛾適生的自然地理條件、寄主成分和傳播條件,因此極有可能遭到美國白蛾的入侵。本研究評估結果表明,美國白蛾在四川省屬于中等危險性有害生物,其風險綜合評價值為1.80。MaxEnt預測表明,除攀枝花市、甘孜藏族自治州、阿壩藏族羌族自治州、涼山彝族自治州4個市(州)外,四川省其余各市幾乎都是美國白蛾適生區(qū),其中成都市最為適宜美國白蛾生存。
目前,美國白蛾在我國呈快速擴散蔓延之勢,截至2015年底,已有北京、天津、遼寧、吉林、河北、河南、山東、陜西、安徽、江蘇等10個省(直轄市、自治區(qū))共500多個縣(區(qū))(國家林業(yè)局公告2016第8號)遭到美國白蛾的入侵危害,造成了巨大損失。隨著我國經濟建設的加快和物流貿易的頻繁,美國白蛾極有可能向內陸省份擴散。與四川省毗鄰的陜西省在1984年就發(fā)現美國白蛾危害,經過多年艱苦的防治得到有效控制,并沒有翻越秦嶺侵入四川[17]。美國白蛾具有繁殖能力強、寄主廣、食性雜、傳播快、耐餓性強、危害嚴重、天敵跟進慢等特點[18-21],已引起各級政府部門的高度重視,做好美國白蛾的監(jiān)測和預警分析勢在必行。四川地處長江上游,是全國森林資源大省、生態(tài)建設大省、林業(yè)經濟大省和物種保護大省。全省林地面積2 400萬hm2,居全國第3位;森林蓄積17.3億m3,居全國第3位;森林面積1 733 萬hm2,居全國第4位;森林覆蓋率36%,是全球生態(tài)多樣性的熱點地區(qū)之一,也是長江上游生態(tài)安全屏障,生態(tài)區(qū)位性極其重要。一旦遭到美國白蛾入侵,很難根除,后果會很嚴重。因此,四川省應提前做好美國白蛾的風險評估和預警分析,在美國白蛾適生區(qū)加大檢疫檢驗力度,嚴防死守,確保該蟲不侵入、不擴散、不危害,盡早建立美國白蛾防控體系和應急預案,并加強美國白蛾宣傳力度,做好技術儲備,打贏這場人與白蛾的防控戰(zhàn)爭。
[1] 宋玉雙.美國白蛾的綜合管理 [M].哈爾濱:東北林業(yè)大學出版社,2015.
Song Y S.Integrate management ofHyphantriacunea[M].Harbin:Forestry Publication Press Housing of Northeast Forestry University,2015.
[2] 李淑賢,高寶嘉,張東風.美國白蛾危險性評估研究 [J].中國農學通報,2009,25(10):202-206.
Li S X,Gao B J,Zhang D F.Studies of risk assessment ofHypanthiacunea(Drury) [J].Chinese Agriculture Bulletin,2009,25(10):202-206.
[3] 金瑞華,魏淑秋,梁憶冰.黑頭型美國白蛾在我國適生地初探 [J].植物檢疫,1991(4):241-246.
Jin R H,Wei S Q,Liang Y B.Adaptability of American white moth with black head in China [J].Plant Qurantine,1991(4):241-246.
[4] 趙鐵珍,高 嵐,柯水發(fā),等.美國白蛾入侵部分非經濟損失評估案例分析 [J].北京林業(yè)大學學報,2007,29(3):171-177.
Zhao T Z,Gao L,Ke S F,et al.Case study on the evaluation of non-economic loss afterHyphantriacuneainvading China [J].Journal of Beijing Forestry University,2007,29(3):171-177.
[5] 魏曉棠,肖海軍,楊 東,等.美國白蛾生物學上的幾個重要特性 [J].植物檢疫,2006,20(1):14-17.
Wei X T,Xiao H J,Yang D,et al.Several important biological characters in the fall webworm,Hyphantriacunea[J].Plant Quarantine,2006,20(1):14-17.
[6] 劉海軍,駱有慶,溫俊寶,等.京津地區(qū)紅脂大小蠢、美國白蛾和銹色粒肩天牛風險階評價 [J].北京林業(yè)大學學報,2005,27(2):81-87.
Liu H J,Luo Y Q,Wen J B,et al.Pest risk assessment ofDendroctonusvalens,HyphantriacuneaandAprionaswainsoniin Beijing area [J].Journal of Beijing Forestry University,2005,27(2):81-87.
[7] 李素霞,張 杰,張 斌,等.美國白蛾(Hyphandriacunea)分布、危害及防治研究進展 [J].世界農藥,2013,35(3):41-46.
Li S X,Zhang J,Zhang B,et al.The distribution,damage and control of fall webworm(Hyphantriacunea) [J].World Pesticides,2013,35(3):41-46.
[8] 王運生,謝丙炎,萬方浩.ROC 曲線分析在評價入侵物種分布模型中的應用 [J].生物多樣性,2007,15(4):365-372.
Wang Y S,Xie B Y,Wan F H.Application of ROC curve analysis in evaluating the performance of alien species’potential distribution models [J].Biodiversity Science,2007,15(4):365-372.
[9] Phillips S J,Anderson R P,Schapire R E.Maximum entropy modeling of species geographic distributions [J].Ecological Modeling,2006,190(3/4):231-259.
[10] Metz C E.Basic principles of ROC analysis [J].Seminars in Nuclear Medicine,1978,8(4):283-298.
[11] 蔣 青,梁憶冰,王乃楊,等.有害生物危險性評價的定量分析方法研究 [J].植物檢疫,1995(4):208-211.
Jiang Q,Liang Y B,Wang N Y,et al.Quantitative analysis of pest risk management [J].Plant Quarantine,1995(4):208-211.
[12] 申衛(wèi)星,郭惠玲,遲延凱,等.美國白蛾在泰山的適生性分析 [J].林業(yè)科學,2012,48(6):165-169.
Shen W X,Guo H L,Chi Y K,et al.Adaptability analysis of American white moth in the Mount Tai [J].Scientia Silvae Sinicae,2012,48(6):165-169.
[13] 林 偉.美國白蛾在中國適生性的初步研究 [D].北京: 北京農業(yè)大學,1991.
Lin W.Preliminary study of adaptability of America white moth in China [D].Beijing:Beijing Agriculture University,1991.
[14] 季英超,秦 萌,吳立峰,等.二種危害小麥的草螟在中國的適生區(qū) [J].應用昆蟲學報,2014,51(1):271-277.
Ji Y C,Qin M,Wu L F,et al.Potential distribution of two wheat pests,AgriphilaaeneociliellaandA.vulgivagella,in China [J].Chinese Journal of Applied Entomology,2014,51(1):271-277.
[15] Beaumont L J,Hughes L,Poulsen M.Predicting species distributions: use of climatic parameters in BIOCLIM and its impact on predictions of species,current and future distributions [J].Ecological Modeling,2005,186(2):251-270.
[16] 韓陽陽,王 眾,項 楊,等.基于MaxEnt生態(tài)位模型的松材線蟲在中國的適生區(qū)預測分析 [J].南京林業(yè)大學學報(自然科學版),2015,39(1):6-10.
Han Y Y,Wang Z,Xiang Y,et al.Prediction of potential distribution ofBursaphelenchusxylophilusin China based on MaxEnt ecological niche model [J].Journal of Nanjing Forestry University(Natural Sciences Edition),2015,39(1):6-10.
[17] 劉 興,李高華,馬友信.美國白蛾在陜西的發(fā)生與檢疫防治 [J].西北大學學報,1998,28(3):258-260.
Liu X,Li G H,Ma Y X.Quarantine and control of fall webworm in Shaanxi province [J].Journal of Northwest University,1998,28(3):258-260.
[18] 鞠 珍,趙 靜,丁福波,等.饑餓程度對美國白蛾生長發(fā)育和繁殖的影響 [J].昆蟲知識,2008,45(3):437-440.
Ju Z,Zhao J,Ding F B,et al.Effects of starvation on the development and reproduction ofHyphantriacunea[J].Chinese Bulletin of Entomology,2008,45(3):437-440.
[19] 魏建榮,楊忠岐,王傳真,等.天敵昆蟲對美國白蛾的生物控制研究 [J].林業(yè)科學,2004,40(2):90-97.
Wei J R,Yang Z Q,Wang C Z,et al.Study on biological control of fall webworm by natural insect enemy [J].Scientia Silvae Sinicae,2004,40(2):90-97.
[20] 楊忠岐.利用天敵昆蟲控制我國重大林木害蟲研究進展 [J].中國生物防治,2004,20(4):221-227.
Yang Z Q.Advance in bio-control researches of the important forest insect pests with natural enemies in China [J].Chinese Journal of Biological Control,2004,20(4):221-227.
[21] 馮術快,盧緒利.北京市昌平區(qū)美國自蛾生物學特性觀察及綜合防治 [J].植物保護,2009,35(5):168-169.
Feng S K,Lu X L.Observation and integrate control of American white moth in Beijing Changping [J].Plant Protection,2009,35(5):168-169.