肖正利，敬順華，劉海玲，周 勇，彭艷紅，王 歡

(遠(yuǎn)安縣林業(yè)局，湖北 宜昌 444200)

松材線蟲病，又名松樹線蟲萎蔫病、松樹萎蔫病、松樹枯萎病，其病原——嗜木傘真滑刃線蟲(Bursaphelenchusxylophius)，是林業(yè)中最具危害性的病害之一，被稱為“松樹癌癥”[1-3]，自1982年被發(fā)現(xiàn)入侵江蘇南京開始[4]，截至2019年底，已迅速擴(kuò)散蔓延到全國(guó)18個(gè)省市666個(gè)縣[5-6]。在溫帶，幾乎所有松樹種類都可以被松材線蟲侵染[2]，造成巨大的經(jīng)濟(jì)損失，對(duì)松林及國(guó)土生態(tài)安全構(gòu)成了嚴(yán)重威脅?？刂七@一重大毀滅性和傳染性病害擴(kuò)散蔓延是林業(yè)有害生物防控亟待解決的重大科學(xué)技術(shù)難題。

松材線蟲自然傳播主要靠媒介昆蟲?？勺鳛槊浇槔ハx的有13種[2]，國(guó)內(nèi)發(fā)現(xiàn)能攜帶松材線蟲的昆蟲有 6 種[7]。該研究區(qū)域自然傳播媒介為昆蟲松褐天牛(Monochamusalternatus，又名松墨天牛[2])，成蟲羽化飛出后在松樹上取食或產(chǎn)卵時(shí)將松材線蟲傳播到鄰近的松樹。人為傳播為以帶有松材線蟲和傳播媒介的病材、病枝及其加工品等通過(guò)人工、運(yùn)輸工具等搬運(yùn)到另一地區(qū)。這種傳播方式不受自然屏障限制，而且速度快，是最主要也最危險(xiǎn)的傳播方式。

松材線蟲病疫情的擴(kuò)散和傳播具有一定的空間分布，獲取疫情侵染源空間位置信息，輔助計(jì)算機(jī)GIS系統(tǒng)，對(duì)疫點(diǎn)空間數(shù)據(jù)信息、位置信息、屬性特征運(yùn)用空間分析技術(shù),解釋事物內(nèi)在的規(guī)律，并在空間查詢、分析、擴(kuò)散模擬、統(tǒng)計(jì)和預(yù)測(cè)等方面也發(fā)揮著巨大的作用。

因此，研究松材線蟲病疫情傳播的空間分布信息，對(duì)疫情傳播擴(kuò)散空間分布數(shù)據(jù)定量化并模擬分析，建立松材線蟲病疫情擴(kuò)散、蔓延、風(fēng)險(xiǎn)評(píng)估模型，可以幫助解決決策、疫情管理、風(fēng)險(xiǎn)預(yù)測(cè)，進(jìn)行疫情監(jiān)測(cè)預(yù)警；在疫情調(diào)查、防治、監(jiān)測(cè)等方面也發(fā)揮著巨大的作用，為松材線蟲防控和決策提供理論依據(jù)，以期遏制松材線蟲病的擴(kuò)散蔓延速度。

1 材料與方法

1.1 研究區(qū)域概況

研究地遠(yuǎn)安縣(111°14′～111°53′E，31°22′～31°53′N)位于湖北省長(zhǎng)江中、上游鄂西北，宜昌之東北,系鄂西北山地向江漢平原過(guò)渡地帶，地貌以中低山為主,地勢(shì)為西北高，東南低。海拔76～1 325 m，平均海拔500 m。屬亞熱帶型大陸季風(fēng)氣候區(qū)，氣候四季分明，自然條件優(yōu)越。日照充足，年均溫12～16℃，極端最高氣溫41.4℃，最低氣溫<-7℃。年太陽(yáng)輻射80～108 kCal/cm，≥10℃活動(dòng)積溫 3 700～4 900℃，年日照時(shí)數(shù) 1 830 h，無(wú)霜期210～275 d，年均相對(duì)濕度79%。降水量充沛，年均降水量多在 1 000～1 100 mm，冬、秋季雨水少，春、夏季降水多，年均蒸發(fā)量 1 300 mm。該區(qū)域人口相對(duì)密集,交通發(fā)達(dá)，經(jīng)濟(jì)活動(dòng)頻繁，人口大規(guī)模流動(dòng)有利于松材線蟲疫情的傳播。

該區(qū)域森林面積13.87萬(wàn)hm2，森林植被類型屬于亞熱帶落葉闊葉林，森林覆蓋率75.6%，地段性森林植被類型為落葉闊葉混交林，以殼斗科(Fagaceae)物種為主體的森林群落居多。松科植物主要為馬尾松，占比26.5%，馬尾松林呈純林和混交分布；闊葉混交林占比67%，其它林占比6.5%,馬尾松在該區(qū)域森林植被類型上占有一定的生態(tài)空間。

該病2006年傳入宜昌市[8],2012年傳入遠(yuǎn)安，2019年遠(yuǎn)安成為疫區(qū)[9]。松材線蟲病入侵馬尾松林致病力強(qiáng)，死亡速度快、數(shù)量大、擴(kuò)散快，治理難度大，猝不及防，給該區(qū)域造成了嚴(yán)重的經(jīng)濟(jì)和生態(tài)損失，對(duì)松林資源構(gòu)成嚴(yán)重威脅。

1.2 數(shù)據(jù)收集和處理

收集該區(qū)域松材線蟲疫情管理小班數(shù)據(jù)庫(kù)(*.SHP格式)，基于ArcGIS軟件平臺(tái)，以林地變更調(diào)查暨森林資源動(dòng)態(tài)監(jiān)測(cè)數(shù)據(jù)庫(kù)為基礎(chǔ)底圖，建立具有空間位置、管理因子、立地因子、生態(tài)因子、林木因子等的松材線蟲病管理數(shù)據(jù)庫(kù)，能與安卓系統(tǒng)GIS軟件平臺(tái)數(shù)據(jù)交互，用于野外作業(yè)，包括普查、除治現(xiàn)場(chǎng)監(jiān)督、疫區(qū)區(qū)劃、檢查驗(yàn)收等；所匯集的數(shù)據(jù)庫(kù)點(diǎn)、面圖層，每個(gè)點(diǎn)視為一個(gè)侵染源(或疫源、疫點(diǎn))。

1.3 研究方法

為便于研究和分析，將松材線蟲疫情數(shù)據(jù)庫(kù)面圖層轉(zhuǎn)成點(diǎn)圖層，采用空間點(diǎn)模式分析；空間點(diǎn)模式分析是對(duì)點(diǎn)狀地理實(shí)體或事件的空間分布和模式進(jìn)行研究的方法，在相關(guān)領(lǐng)域得到廣泛應(yīng)用[10]；它是以空間坐標(biāo)為主導(dǎo),以屬性數(shù)據(jù)作為權(quán)重進(jìn)行統(tǒng)一描述[11]，每個(gè)疫點(diǎn)事件都可以抽象化為空間上的一個(gè)位置點(diǎn)，用來(lái)描述任何類型的事件數(shù)據(jù)[12]。研究方法采用空間可視化、三維趨勢(shì)、空間自相關(guān)、多環(huán)緩沖區(qū)、克里金插值5種方法[13-16]。所有數(shù)據(jù)基于ArcGIS 10.8軟件平臺(tái)分析、統(tǒng)計(jì)。

1.3.1空間可視化分析

將采集獲取的松材線蟲疫點(diǎn)地面屬性數(shù)據(jù)和空間數(shù)據(jù)信息，以疫點(diǎn)坐標(biāo)經(jīng)度、緯度和發(fā)生株數(shù)等建立地理信息數(shù)據(jù)庫(kù)，疊加基礎(chǔ)地理地圖為工作底圖，通過(guò)地理信息方法進(jìn)行可視化，直觀地顯示出疫情的空間分布狀況，揭示疫情與某些事物之間的關(guān)聯(lián)，繪制松材線蟲病疫情空間分布圖，并做趨勢(shì)分析(Trend analysis)。

1.3.2多環(huán)緩沖區(qū)(Multiple Ring Buffer)分析

松材線蟲自然擴(kuò)散需依靠松褐天牛成蟲攜帶遷飛擴(kuò)散來(lái)實(shí)現(xiàn)。松褐天牛自然擴(kuò)散能力不強(qiáng)，成蟲羽化后主要飛翔目的是取食、產(chǎn)卵、尋找配偶，自然遷飛擴(kuò)散能力有限。松褐天牛飛行潛能被激發(fā)時(shí)，在野外一次性飛行的距離也僅能達(dá)到0.8～1 km[17]，以此規(guī)律，通過(guò)疫點(diǎn)病原的自然擴(kuò)散，可建立空間擴(kuò)散模型，以對(duì)象疫點(diǎn)為基本核心，根據(jù)指定鄰域距離，在其周圍建立不同半徑寬度的緩沖區(qū)域；運(yùn)用多環(huán)緩沖區(qū)工具，對(duì)所有要素點(diǎn)(疫點(diǎn))按設(shè)定的距離條件圍繞其要素中心形成一定緩沖區(qū)多邊形實(shí)體，實(shí)現(xiàn)數(shù)據(jù)在二維空間得以擴(kuò)展的信息，計(jì)算得到疫區(qū)自然擴(kuò)散指數(shù)的空間模型圖。公式為：

B={x││d(x,A)≤R}

(1)

式中：A疫點(diǎn)半徑為R的緩沖區(qū)，以對(duì)象A為核心，建立半徑為d的多邊形區(qū)域范圍內(nèi)所有點(diǎn)的集合。d取值為最小歐氏距離，R為鄰域半徑即緩沖區(qū)建立的條件。

1.3.3克里金插值(Kriging)分析

克里金插值法是以變異函數(shù)理論和結(jié)構(gòu)分析為基礎(chǔ)，在有限區(qū)域內(nèi)對(duì)區(qū)域化變量進(jìn)行無(wú)偏最優(yōu)估計(jì)的一種方法，是對(duì)未知點(diǎn)特定鄰域范圍的測(cè)量點(diǎn)或者特定數(shù)量的相鄰測(cè)量點(diǎn)數(shù)值進(jìn)行加權(quán)相加，以求得未知點(diǎn)的數(shù)值，實(shí)現(xiàn)對(duì)未知點(diǎn)的預(yù)測(cè)。其中周圍測(cè)量點(diǎn)的權(quán)重根據(jù)半方差函數(shù)來(lái)確定[16,18]。利用材線蟲病疫點(diǎn)已知空間分布和疫點(diǎn)松林受害屬性值，估計(jì)空間上任意一點(diǎn)疫情的屬性值?？死锝鸩逯倒綖椋?/p>

(2)

(3)

式中：Z(x0)為x0處的估計(jì)值，Z(xi)為xi處的屬性值，λi為克里金權(quán)重系數(shù)，n為點(diǎn)個(gè)數(shù)。

1.3.4空間自相關(guān)(Global Moran′s I)分析

空間自相關(guān)分析是一種空間統(tǒng)計(jì)方法，可揭示空間變量的區(qū)域結(jié)構(gòu)形態(tài),是檢驗(yàn)?zāi)匙兞吭谔囟臻g位置上的屬性值是否與其相鄰空間位置上的屬性值顯著相關(guān)[15-16]?？臻g自相關(guān)分為全局和局部空間自相關(guān)。前者是檢測(cè)變量的屬性值在空間上是否具有相關(guān)性和規(guī)律性，后者是表示變量的屬性值在空間上的聚集或者離散區(qū)域。Moran′s I指數(shù)從整體上反映研究區(qū)域的空間相關(guān)性，取值范圍是[-1，1]。當(dāng)I>0時(shí)，表示正相關(guān)，且越接近于1 聚集性越大；當(dāng)I<0，表示負(fù)相關(guān)；I=0，是隨機(jī)現(xiàn)象[12，15]。局部空間自相關(guān)，是在全局分析基礎(chǔ)上沒有聚集出現(xiàn)時(shí)，用來(lái)檢查是否有小的、可能被忽視的聚集，可反映子區(qū)域間的相關(guān)程度。

利用該方法對(duì)松材線蟲疫點(diǎn)區(qū)域空間傳播的空間自相關(guān)性進(jìn)行實(shí)證分析，描述松材線蟲病疫點(diǎn)空間整體分布狀況，檢驗(yàn)要素(疫點(diǎn))屬性值是否與其相鄰空間點(diǎn)上的屬性值相關(guān)聯(lián)，其結(jié)果為正相關(guān)和負(fù)相關(guān)。正相關(guān)表明屬性值變化與其相鄰空間單元具有相同的變化趨勢(shì)，代表空間現(xiàn)象有集聚性的存在；負(fù)相關(guān)則為隨機(jī)現(xiàn)象。全局空間自相關(guān)的Moran＇s I 的計(jì)算公式為：

(4)

(5)

(6)

E[I]=-1/(n-1)

(7)

V[I]=E[I2]-E[I]2

(8)

2 結(jié)果與分析

2.1 松材線蟲病發(fā)生概況

2019年度，該區(qū)域松材線蟲病發(fā)生范圍分布于7個(gè)鄉(xiāng)鎮(zhèn)83個(gè)村，受害松林面積 2 345.6 hm2，疫點(diǎn)(小班) 1 976個(gè)。疫點(diǎn)松林受害1～800株，年度清理枯死、萎蔫病松樹 124 408株(表1)。危害程度處于輕微，處置方式為現(xiàn)場(chǎng)焚燒，傳播途徑為自然傳播和人為傳播。疫情發(fā)生的區(qū)域存在較高的自然傳播擴(kuò)散風(fēng)險(xiǎn)，當(dāng)前，雖然成災(zāi)率發(fā)生在3‰以下，加之采取了防治措施，除治率達(dá)到100%,疫情可控可防。

表1 松材線蟲病疫點(diǎn)概況

2.2 空間可視化分析

利用ArcGIS 軟件平臺(tái)對(duì)疫點(diǎn)圖層進(jìn)行符號(hào)的設(shè)置，以疫點(diǎn)字段株數(shù)屬性值，采用自然間斷點(diǎn)分級(jí)法，將松林受害株數(shù)從低到高的順序分成8個(gè)等級(jí)，每個(gè)等級(jí)100株，由點(diǎn)的大小以地圖可視化形式來(lái)表示疫情空間分布，得到疫點(diǎn)散點(diǎn)分布圖(圖1)，其分布特點(diǎn)從聚集到稀疏，小點(diǎn)居多。

圖1 疫點(diǎn)分布Fig.1 Distribution of epidemic points

由圖1可知，所有鄉(xiāng)鎮(zhèn)都有疫情發(fā)生，且分布較為分散，受害松林嚴(yán)重區(qū)域?yàn)榛炙骆?zhèn)、茅坪場(chǎng)鎮(zhèn)，疫情高發(fā)地區(qū)主要集中在東南地區(qū)，為高風(fēng)險(xiǎn)區(qū)，且疫情存在隔年或翌年發(fā)病現(xiàn)象[19]，疫點(diǎn)短期內(nèi)難以根除。

2.3 三維趨勢(shì)分析

利用空間探索數(shù)據(jù)分析(ESDA)將疫點(diǎn)數(shù)據(jù)做三維趨勢(shì)分析，將區(qū)域內(nèi)每一個(gè)疫點(diǎn)松林受害株數(shù)值看作是該單元格幾何中心處的值,用三維數(shù)據(jù)點(diǎn)(X,Y,Z)來(lái)表示,其中：X和Y分別代表各區(qū)域幾何中心處的經(jīng)度和緯度方向,(i)表示疫點(diǎn)受害松林株數(shù)值，i分別在平面XY和YZ上進(jìn)行投影散點(diǎn)圖,通過(guò)散點(diǎn)做多項(xiàng)式擬合,得到一條曲線,用此曲線模擬疫點(diǎn)松林受害株數(shù)值在經(jīng)緯度方向上的某種變化趨勢(shì)。計(jì)算所得該區(qū)域疫情三維趨勢(shì)圖(圖2)，X軸為方向東，Y軸為方向北，Z軸代表疫點(diǎn)受害松林株數(shù)值，XY面為疫點(diǎn)，疫點(diǎn)桿的高低代表受害株數(shù)值。從桿的高低分析疫情趨勢(shì)為東南高，西北低；從投影趨勢(shì)線值分析，趨勢(shì)線近似于直線，總體較平、低，表明未造成大的受災(zāi)或嚴(yán)重區(qū)域，危害程度處于較輕。

圖2 疫點(diǎn)三維趨勢(shì)Fig.2 Three-dimensional trend of epidemic sites

2.4 空間自相關(guān)分析

對(duì)所有疫點(diǎn)作全局整體分析，所涉距離閥值以媒介昆蟲遷飛距離，取值1 km,即每個(gè)疫點(diǎn)搜索周圍半徑1 km范圍疫點(diǎn)，檢驗(yàn)空間位置上的屬性值是否與其相鄰空間位置上的屬性值相關(guān)。全局自相關(guān)計(jì)算結(jié)果(表2)：Moran′s I=-0.222 6，Z=2.5578，P=0.010 5，即莫蘭指數(shù)(Moran′s I)取值范圍I<0，為負(fù)相關(guān)關(guān)系，為隨機(jī)現(xiàn)象，聚類模式的可能性小于5%，表明疫源傳播是隨機(jī)擴(kuò)散現(xiàn)象，傳播范圍內(nèi)任何一株松樹均存在侵染的風(fēng)險(xiǎn)。

表2 空間自相關(guān)指數(shù)

2.5 多環(huán)緩沖區(qū)分析

利用緩沖區(qū)分析可預(yù)測(cè)疫情傳播范圍，并依據(jù)其范圍進(jìn)行合理的監(jiān)測(cè)設(shè)防。松材線蟲病疫情的傳播擴(kuò)散與媒介昆蟲松褐天牛成蟲遷飛距離相關(guān)，通過(guò)距離可建立疫源的空間擴(kuò)散模型。因此，根據(jù)自然擴(kuò)散傳播距離和疫點(diǎn)的分布狀況，按松褐天牛成蟲遷飛距離潛能被激發(fā)飛行距離1 km傳播擴(kuò)散發(fā)展趨勢(shì),以等間距對(duì)疫點(diǎn)為圓心分成10個(gè)等級(jí)多環(huán)緩沖區(qū)分析，每個(gè)緩沖區(qū)100 m,計(jì)算得到疫區(qū)自然擴(kuò)散指數(shù)的空間模型圖(圖3)。由圖3可知,模型圖體現(xiàn)了空間連續(xù)影響和空間疊加效應(yīng)，由環(huán)內(nèi)向外，風(fēng)險(xiǎn)依次降低，連續(xù)的效果反映了疫點(diǎn)傳播的影響因素隨距離的漸變效果，青色區(qū)域?yàn)榫嚯x疊加效應(yīng)，風(fēng)險(xiǎn)最高。對(duì)緩沖區(qū)距離分析結(jié)果分布狀況進(jìn)行評(píng)估，整個(gè)環(huán)內(nèi)存在傳播的高風(fēng)險(xiǎn)區(qū)及傳播的風(fēng)險(xiǎn)隨著距離的漸遠(yuǎn)而降低。

圖3 疫點(diǎn)多環(huán)緩沖區(qū)分布Fig.3 Epidemic point multi-ring buffer distribution

2.6 克里金插值預(yù)測(cè)分析

以松褐天牛攜帶松材線蟲自然遷飛擴(kuò)散距離1 km作為疫源點(diǎn)傳播的范圍邊界值，即擴(kuò)散距離范圍閾值取半徑R=1 km。對(duì)疫點(diǎn)使用克里金插值法分析計(jì)算并作圖，采用自然間斷點(diǎn)分級(jí)法將受災(zāi)區(qū)域分成8個(gè)等級(jí)來(lái)比較，預(yù)測(cè)疫點(diǎn)不同受害株數(shù)擴(kuò)散的范圍，得到克里金插值結(jié)果(圖4)。如圖4所示，克里金插值預(yù)測(cè)結(jié)果與多環(huán)緩沖區(qū)范圍分布圖大致相同，結(jié)果均能大致反映出該區(qū)域的松林受害株數(shù)空間分布特征。松林受害風(fēng)險(xiǎn)的分布具有相對(duì)集中的特點(diǎn),分別以55株為倍數(shù)遞增的8種不同顏色區(qū)域估計(jì)了松林受害程度，區(qū)間值在0～440，其結(jié)果特征為疫點(diǎn)聚集區(qū)域向外圍逐漸遞減，松林受害株數(shù)呈現(xiàn)聚集分布。因此，其結(jié)果可以對(duì)其進(jìn)行風(fēng)險(xiǎn)分區(qū)的劃分。

圖4 疫點(diǎn)克里金插值Fig.4 Kriging interpolation map of the epidemic point

3 結(jié)論與討論

本文基于GIS空間分析功能，將空間數(shù)據(jù)管理功能、空間分析建模方法、空間可視化同松材線蟲病疫情相結(jié)合，闡述松材線蟲病疫情空間分布特征，以及風(fēng)險(xiǎn)性估計(jì)建模理論、方法。

3.1 結(jié)論

1)空間可視化直觀地反映了疫情發(fā)生分布區(qū)域、受害松林?jǐn)?shù)據(jù)值空間關(guān)系，該區(qū)域2019年度所有鄉(xiāng)鎮(zhèn)都為疫區(qū)，且疫點(diǎn)分布較廣、較多，疫點(diǎn)主要集中在縣域東南地區(qū)，自然傳播擴(kuò)散風(fēng)險(xiǎn)加劇。

2)三維空間趨勢(shì)圖反映了疫情在經(jīng)緯度方向上的整體變化情況，該區(qū)域受害松林株數(shù)桿值趨勢(shì)東南高，西北低，疫點(diǎn)圖投影多項(xiàng)式擬合曲線近似于直線，總體較平、低，疫情在經(jīng)緯度方向上變化趨勢(shì)較低，表明未造成大的受災(zāi)或嚴(yán)重區(qū)域，危害程度較輕。

3)探索了疫情屬性值在空間位置間的相關(guān)關(guān)系，運(yùn)用全局空間自相關(guān)分析整個(gè)區(qū)域疫情的空間分布模式,搜索每個(gè)疫點(diǎn)周圍半徑1 km范圍，其結(jié)果Moran′s I<0，為負(fù)相關(guān)關(guān)系，表明疫點(diǎn)傳播擴(kuò)散是隨機(jī)分布，疫情的傳播是由疫點(diǎn)向四周隨機(jī)擴(kuò)散。

4)以疫點(diǎn)為圓心，等間距100 m的10個(gè)等級(jí)多環(huán)緩沖區(qū)自然擴(kuò)散指數(shù)的空間模型圖，體現(xiàn)了空間連續(xù)影響和空間疊加效應(yīng)，連續(xù)的效果較形象地反映了疫情擴(kuò)散的影響因素隨距離的漸變效果。整個(gè)環(huán)內(nèi)存在傳播的高風(fēng)險(xiǎn)區(qū)，由環(huán)內(nèi)向外，風(fēng)險(xiǎn)依次降低，傳播的風(fēng)險(xiǎn)隨著距離的漸遠(yuǎn)而降低。

5)采用克里金插值結(jié)果預(yù)測(cè)疫點(diǎn)不同等級(jí)受害株數(shù)傳播擴(kuò)散的范圍、受害株數(shù)程度，其預(yù)測(cè)范圍與疫點(diǎn)多環(huán)緩沖區(qū)分布范圍大致相同，均能大致反映出該區(qū)域的松林受害株數(shù)區(qū)域空間分布特征，其特征為疫點(diǎn)聚集區(qū)域向外圍逐漸遞減。松林受害風(fēng)險(xiǎn)的分布具有相對(duì)集中的特點(diǎn)，松林受害株數(shù)呈現(xiàn)聚集分布，其結(jié)果可以對(duì)其進(jìn)行風(fēng)險(xiǎn)分區(qū)的劃分，因此克里金插值方法較優(yōu)。

3.2 討論

1)松材線蟲病疫情的發(fā)生、流行與地理環(huán)境之間有著密切的關(guān)系，研究疫情地理分布規(guī)律，是科學(xué)監(jiān)測(cè)和防治的前期手段，在疫情防控中得到了廣泛的應(yīng)用；本文以遠(yuǎn)安縣區(qū)域年度松材線蟲病疫點(diǎn)數(shù)據(jù)庫(kù)為基礎(chǔ)，運(yùn)用疫點(diǎn)空間可視化、三維趨勢(shì)分析、空間自相關(guān)、多環(huán)緩沖區(qū)、克里金插值5種方法對(duì)松線線蟲疫情傳播風(fēng)險(xiǎn)性估計(jì)建模，在擴(kuò)散蔓延的預(yù)測(cè)研究方面都有重要意義。

2)影響松材線蟲病疫情發(fā)生、擴(kuò)散、蔓延的因子很多，絕大部分探討疫點(diǎn)或影響因素的空間位置，忽略了空間信息；本研究存在一定局限性，松材線蟲病疫情的傳播因素與媒介昆蟲的密度、環(huán)境因子、氣候、海拔、松林的分布及密度、人為活動(dòng)干擾等息息相關(guān)，這些因素需要進(jìn)一步疊加研究，其分析估計(jì)結(jié)果不可能精確到具體的山頭地塊和劃分風(fēng)險(xiǎn)等級(jí)；其數(shù)據(jù)量化研究與模擬方法可適用于有害生物發(fā)生、發(fā)展的預(yù)測(cè)預(yù)報(bào)。

3)該區(qū)域松材線蟲的入侵?jǐn)U展迅速，侵染面積逐年擴(kuò)大，嚴(yán)重影響了縣域的生態(tài)環(huán)境建設(shè)和經(jīng)濟(jì)發(fā)展，防治上耗費(fèi)了大量人力、物力、財(cái)力，目前還未取得階段性根除進(jìn)展。該疫情不僅給國(guó)民經(jīng)濟(jì)造成巨大損失，也破壞了自然景觀及生態(tài)環(huán)境，對(duì)豐富的松林資源構(gòu)成嚴(yán)重威脅。研究松材線蟲病疫情空間分布特征，基于風(fēng)險(xiǎn)性估計(jì)建模、可視化分析，倡導(dǎo)在空間格局上采取自然隔離、根除侵染源、阻擋疫源擴(kuò)散，逐步減小范圍，最終達(dá)到移除疫區(qū)。