發(fā)布時間：所屬分類：農(nóng)業(yè)論文瀏覽：1次
摘 要： 摘要：【目的】當(dāng)前湖北省松材線蟲Bursaphelenchusxylophilus疫情形勢嚴(yán)峻。通過研究湖北省松材線蟲病的入侵風(fēng)險，分析松材線蟲病在湖北省的危害程度，為當(dāng)前疫情的防治工作提供建議和參考。【方法】結(jié)合3S技術(shù)，應(yīng)用生態(tài)位因子模型(ENFA)，選取影響松材線蟲
摘要：【目的】當(dāng)前湖北省松材線蟲Bursaphelenchusxylophilus疫情形勢嚴(yán)峻。通過研究湖北省松材線蟲病的入侵風(fēng)險，分析松材線蟲病在湖北省的危害程度，為當(dāng)前疫情的防治工作提供建議和參考。【方法】結(jié)合“3S”技術(shù)，應(yīng)用生態(tài)位因子模型(ENFA)，選取影響松材線蟲定殖和傳播的4類影響因素(氣候、植被、地形、人類活動)，對松材線蟲病在湖北的入侵風(fēng)險進(jìn)行了預(yù)測和評價。【結(jié)果】湖北省松材線蟲病高風(fēng)險區(qū)面積38884.62km2，占湖北省總面積的20.92%，高風(fēng)險區(qū)主要集中于海拔較低，人類活動頻繁的的中部和東部地區(qū)，中風(fēng)險區(qū)面積66501.84km2，占總面積的35.77%，低風(fēng)險區(qū)面積80513.54km2，占總面積的43.31%，中低風(fēng)險地區(qū)主要分布在林地稀少的江漢平原和西部的高海拔山地。松材線蟲偏好分布在溫度較高，降水豐富，海拔較低，離人類居住點較近且人類活動頻繁的針葉林地區(qū)。通過交叉驗證(cross-validation)對模型的預(yù)測進(jìn)行檢驗，得到P/E曲線，曲線呈單調(diào)遞增且Boyce指數(shù)很高，說明模型精度很高。
【結(jié)論】ENFA模型能很好模擬松材線蟲病的風(fēng)險區(qū)域，模型結(jié)果可為湖北省各縣市的松材線蟲病防治決策提供參考。圖2表3參25關(guān)鍵詞：森林保護(hù)學(xué);松材線蟲;生態(tài)位因子模型;風(fēng)險評估
松材線蟲Bursaphelenchusxylophilus是危害性極高的外來有害生物，也是中國林業(yè)檢疫性有害生物[1]。它主要通過松墨天牛Monochamusalternatus等媒介昆蟲在松樹體內(nèi)傳播，引發(fā)破壞性極強(qiáng)的松材線蟲病[2]。自1982年在南京中山陵首次發(fā)現(xiàn)被松材線蟲感染的黑松Pinusthunbergia以來，松材線蟲病已經(jīng)給中國造成直接和間接損失上千億元[3]。截至2020年4月，松材線蟲病已擴(kuò)散至湖北省的86個市縣地區(qū)，湖北省的林業(yè)生態(tài)安全面臨著巨大的威脅[4]。運用生態(tài)位模型對外來有害生物的適宜生境和風(fēng)險進(jìn)行評估是當(dāng)前研究熱點之一[5−6]，比較廣泛的生態(tài)位模型有BIOCLIM、GARP、GLM、MAXENT、ENFA等，這些模型都建立在HUTCHINSON[7]生態(tài)位原理的基礎(chǔ)上，即將物種的分布視為與環(huán)境因素異質(zhì)性密切相關(guān)的因變量，通過計算物種點的環(huán)境因素變量來確定物種的生態(tài)位，然后利用物種的生態(tài)位模擬研究該物種在研究區(qū)域內(nèi)的分布適宜性。生態(tài)位因子分析(ecological-nichefactoranalysis，ENFA)模型的優(yōu)點在于模型計算只需要已出現(xiàn)的物種點數(shù)據(jù)，而不需要非出現(xiàn)點的數(shù)據(jù)[8]。并且相較于其他生態(tài)位模型，ENFA模型結(jié)果還包括各指標(biāo)因子與物種的相關(guān)關(guān)系，可對物種的生態(tài)位直接進(jìn)行解釋。ENFA模型已成功模擬了大熊貓Ailuropodamelanoleuca[8]和亞洲象Elephasmaximus的適宜生境[9]、人象沖突的風(fēng)險評估[10]、蝴蝶Rhopalocera的異地適生區(qū)[11]、蕨類Pteridophyta植物的空間分布預(yù)測[12]、林麝Moschusberezovskii的生境預(yù)測[13]、歐洲松雞Tetraourogallus的保護(hù)[14]等。本研究結(jié)合“3S”技術(shù)，應(yīng)用ENFA模型，對湖北省松材線蟲病的潛在風(fēng)險區(qū)進(jìn)行評估，研究影響湖北省松材線蟲定殖和傳播的主要影響因子，分析風(fēng)險區(qū)的面積大小和分布狀況，以期為湖北省松材線蟲病的防治提供科學(xué)依據(jù)。
1數(shù)據(jù)來源與研究方法
1.1數(shù)據(jù)來源
物種分布點數(shù)據(jù)來源于國家林業(yè)和草原局公告(2020年第4號)[15]，共得到松材線蟲病在湖北省(29°01′53″~33°06′47″N，108°21′42″~116°07′50″E)的分布點數(shù)據(jù)86條(圖1)。數(shù)據(jù)包括：①氣候數(shù)據(jù)。主要為2008−2017年湖北省34個氣象站點的數(shù)據(jù)，來源于國家氣象科學(xué)數(shù)據(jù)共享服務(wù)平臺(http://data.cma.cn)。②地形數(shù)據(jù)。主要為湖北省的數(shù)字高程模型(DEM)數(shù)據(jù)(2009年)，來自于地理空間數(shù)據(jù)云(http://www.gscloud.cn)。③人類活動干擾數(shù)據(jù)。包括道路數(shù)據(jù)(2015年)、土地利用/覆蓋變化數(shù)據(jù)(2015年)、人口分布密度數(shù)據(jù)(2015年)，均來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn)。④植被因素。主要包括植被覆蓋類型(2015年)和歸一化植被指數(shù)(NDVI)數(shù)據(jù)(2017年)，來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn)。
1.2生態(tài)位因子分析模型(ENFA)
ENFA模型是研究物種分布的多變量分析方法[16]。該模型采用主成分分析方法，利用物種分布點的數(shù)據(jù)和生態(tài)地理變量(EGV)數(shù)據(jù)，將物種的生態(tài)位和生態(tài)幅度與整個研究區(qū)的環(huán)境因子平均狀態(tài)和標(biāo)準(zhǔn)差相比較，從EGV中提取主要信息組合為互不相關(guān)的特征矩陣，特征矩陣和特征向量用來生成生境適宜性圖[17]。物種的分布與EVG之間的相互關(guān)系通過計算邊際性、特殊性和耐受性3個指標(biāo)得到[18]。邊際性系數(shù)(M)為物種分布區(qū)EGV平均值和整個研究區(qū)內(nèi)EGV平均值的差異。M>0時，表明該物種偏好EGV平均水平以上的環(huán)境;M<0時，表明物種偏好EGV平均水平以下的環(huán)境。M的絕對值越大，物種對EGV的偏好程度越高，物種的分布也越不隨機(jī)。
1.3ENFA模型指標(biāo)體系
松材線蟲病的定殖和傳播是一個復(fù)雜的病害系統(tǒng)。溫度、降水量、日照等氣候條件影響松材線蟲及其寄主松墨天牛的生存和定殖[19];林分狀況、樹齡、樹種的豐富度、森林郁閉度等因子影響著媒介昆蟲松墨天牛的種群密度[20];海拔、坡度和坡向等地形因子通過影響氣候和土壤從而間接影響植物的生長和分布;松墨天牛的自然擴(kuò)散能力不強(qiáng)[2]，因此松材線蟲主要借助人為傳播，人類活動直接影響松材線蟲病擴(kuò)散的速度和規(guī)模。本研究從氣候、植被、地形、人類活動干擾等4個影響因素中選取15個EGV指標(biāo)，用于松材線蟲病潛在風(fēng)險區(qū)的評估(表1)。
1.4數(shù)據(jù)處理
所有生態(tài)地理變量數(shù)據(jù)均預(yù)先在ArcGIS10.2中進(jìn)行處理。氣候數(shù)據(jù)進(jìn)行克里金(Kriging)插值處理，地形數(shù)據(jù)從DEM提取，植被和人類活動干擾數(shù)據(jù)需轉(zhuǎn)換為距離數(shù)據(jù)或者頻率數(shù)據(jù)[19]，植被覆蓋度利用像元二分模型計算NDVI得到[21]。所有數(shù)據(jù)轉(zhuǎn)換為250m×250m空間分辨率和相同的投影，并將ArcGIS柵格數(shù)據(jù)轉(zhuǎn)換為ENFA所能使用的柵格數(shù)據(jù)格式，然后利用ENFA對數(shù)據(jù)進(jìn)行Box-Cox標(biāo)準(zhǔn)化后進(jìn)行分析。
2結(jié)果與分析
2.1ENFA模型結(jié)果
通過ENFA模型得到邊際性系數(shù)M=1.639，M>1，說明松材線蟲對各種環(huán)境條件的選擇不是隨機(jī)的。特殊性系數(shù)S=8.822，耐受性系數(shù)T=0.113，S較大，T小，說明松材線蟲的定殖和傳播是受特定環(huán)境條件影響，在環(huán)境條件適宜的情況下入侵風(fēng)險較高。
相關(guān)期刊推薦：《浙江農(nóng)林大學(xué)學(xué)報》是中國林業(yè)類和綜合性農(nóng)業(yè)科學(xué)類核心期刊之一，主要報道林學(xué)基礎(chǔ)學(xué)科、森林培育學(xué)、森林經(jīng)理學(xué)、經(jīng)濟(jì)林學(xué)、林業(yè)工程、森林保護(hù)學(xué)、林木遺傳育種學(xué)、森林生物學(xué)、生態(tài)學(xué)、生物技術(shù)、園林學(xué)和園藝學(xué)等學(xué)科的學(xué)術(shù)論文、問題討論和研究簡報，適當(dāng)刊登與農(nóng)林相關(guān)的其他學(xué)科的稿件。
從表2可知：第1列為邊際因子，特殊性解釋比例為98.00%，其他為特殊因子，包括特異因子1和特異因子2，解釋比例分別為1.11%、0.49%，3個因子的累積貢獻(xiàn)率達(dá)99.60%，即解釋了100%的邊際性和99.60%的特殊性，由這3個因子生成最終的風(fēng)險等級圖。用交叉驗證對模型的預(yù)測進(jìn)行檢驗，得到呈單調(diào)遞增趨勢的P/E曲線，Boyce平均值為0.90，標(biāo)準(zhǔn)差為0.20，說明模型預(yù)測能力很好。根據(jù)P/E曲線將湖北省松材線蟲病風(fēng)險區(qū)分為3個等級：低風(fēng)險區(qū)(0
2.3松材線蟲的潛在風(fēng)險區(qū)分析
利用ArcGIS10.2對湖北省松材線蟲的入侵風(fēng)險等級區(qū)域進(jìn)行計算，發(fā)現(xiàn)湖北省松材線蟲病高風(fēng)險區(qū)面積38884.62km2，占湖北省總面積的20.92%，中風(fēng)險區(qū)面積66501.84km2，占總面積的35.77%，低風(fēng)險區(qū)面積80513.54km2，占總面積的43.31%(表3)。
湖北省松材線蟲病入侵風(fēng)險較高的地區(qū)主要有：荊門市東部和中西部大部分地區(qū);荊州市的西南部和南部地區(qū);宜昌市的東部和中部地區(qū);天門市的北部地區(qū);襄樊市的東北部和中部大部分地區(qū);孝感市的北部大部分地區(qū);黃岡市的東南部和北部的大部分地區(qū);黃石市的東部和中部大部分地區(qū);咸寧市的大部分地區(qū);十堰市的中部地區(qū);恩施市的中部和東北部的少部分地區(qū);隨州市的東部和中部地區(qū);武漢市的西北部。這些地區(qū)基本都是海拔較低的山地和平原地區(qū)，森林覆蓋率高，人類活動頻繁，導(dǎo)致松材線蟲入侵的風(fēng)險也相對較高。
中風(fēng)險區(qū)和低風(fēng)險區(qū)主要分布在江漢平原的大部分地區(qū)和西部的高海拔山地，主要包括鄂州市、仙桃市、潛江市、天門市、武漢市的大部分地區(qū)、荊州市的東部和中部大部分地區(qū)、黃岡市東南部的小部分地區(qū)、宜昌市的南部和西部大部分地區(qū)、恩施州的大部分地區(qū)、襄樊市的西部大部分地區(qū)和十堰市的西南部大部分地區(qū)以及與之接壤的神農(nóng)架林區(qū)。江漢平原地勢平坦，人口密集，但森林很少，大部分低風(fēng)險區(qū)則分布于此;湖北省西部大部分地區(qū)都是海拔較高的山地，特別是神農(nóng)架林區(qū)的平均海拔就有1700m，一部分低風(fēng)險區(qū)就位于這些海拔很高的地區(qū);西部地區(qū)人類活動較少，但一些低海拔地區(qū)仍然有森林覆蓋，因此有大部分中度風(fēng)險區(qū)分布于此。
3結(jié)論與討論
較前人研究[5−6,22]，本研究在模型指標(biāo)因子的選擇上更加全面，不局限于氣候因素，因此模型結(jié)果可以更好地了解松材線蟲病在湖北省的傳播規(guī)律。模型通過交叉驗證，Boyce平均值為0.90，標(biāo)準(zhǔn)差為0.20，很好地模擬了湖北省松材線蟲病的潛在風(fēng)險區(qū)域。
本研究表明：年均氣溫、年降水量、夏季氣溫及降水均是影響松材線蟲定殖的主要氣候因素。湖北省適宜的氣候給松材線蟲及其寄主松墨天牛提供了良好的生存環(huán)境。地形因子的邊際性系數(shù)表明：松材線蟲喜歡分布于海拔較低，坡度平緩，并位于西坡和南坡的地區(qū)，因此松材線蟲病普遍發(fā)生于海拔低于700m的地區(qū)。植被因子的邊際性系數(shù)表明：松材線蟲偏好于植被覆蓋度較低的針葉林。人口密度、居民點、鐵路和農(nóng)田是影響松材線蟲病在湖北省內(nèi)傳播的幾個主要人為因素。
根據(jù)模型預(yù)測結(jié)果，湖北省松材線蟲病入侵風(fēng)險較高的區(qū)域廣泛分布于湖北省的12個地級市，部分地區(qū)已經(jīng)暴發(fā)了比較嚴(yán)重的松材線蟲病，這些區(qū)域都是海拔較低的山地或者平原，氣候適宜，人類活動頻繁，松材線蟲入侵的風(fēng)險特別高。因此在風(fēng)險較高的地區(qū)必須加強(qiáng)松材線蟲病的防范和治理力度。首先，必須加大高風(fēng)險區(qū)內(nèi)針葉林的疫情監(jiān)測力度，在夏季運用遙感和無人機(jī)監(jiān)測等相關(guān)技術(shù)，實現(xiàn)松材線蟲病疫情的大面積監(jiān)測;其次，針對已經(jīng)感染松材線蟲病的林區(qū)，可采取目前較為成熟的生物防治、引誘劑和化學(xué)防治來降低媒介昆蟲的蟲口數(shù)量，從而減小松材線蟲病疫區(qū)的范圍;松材線蟲病防治最重要的環(huán)節(jié)就是染病松木的處理，松材線蟲病疫區(qū)必須徹底清理染病瀕死和枯死的松木，在受災(zāi)程度嚴(yán)重的地區(qū)則應(yīng)建立疫木集中處理場地，對砍伐的疫木進(jìn)行集中處理[23];同時，政府和相關(guān)部門也要做好與松材線蟲病害相關(guān)的宣傳和管理工作，嚴(yán)格控制松木及其加工制品的流通，防止松材線蟲借助松木制品再次擴(kuò)散到其他地區(qū)[24]。
中風(fēng)險區(qū)和高風(fēng)險區(qū)相鄰，大部分針葉林也都分布在中風(fēng)險地區(qū)，因此在這些地區(qū)主要做好松材線蟲病的監(jiān)測和檢疫工作，對于新發(fā)病的地區(qū)要及時處理疫木，砍伐的疫木盡快就地焚燒或粉碎處理。
低風(fēng)險區(qū)主要分布在江漢平原中部和西北部的山區(qū)和神農(nóng)架林區(qū)。江漢平原中部地勢平緩，氣候適宜，但是松材線蟲的寄主植物松樹卻很少，西北部山區(qū)和神農(nóng)架林區(qū)海拔較高，大部分松林都分布于人跡罕至的地區(qū)，因此感染松材線蟲的風(fēng)險相對較低。但隨著城鎮(zhèn)的擴(kuò)張，人類活動的加劇，氣候變暖和松材線蟲的遺傳多樣性提高[25]，其感染風(fēng)險也會相對提高，因此也要做好松材線蟲病的檢疫工作，并對松林進(jìn)行定期巡查，阻止松材線蟲的進(jìn)一步蔓延。——論文作者：沈鵬，李功權(quán)