崔紹朋 王澤鑫 吉晟男

摘 要：研究有害物種分布格局，對制定行之有效的預防和檢疫措施至關(guān)重要。為明確雙條杉天牛在山西省境內(nèi)的適生區(qū)及其對環(huán)境因素的響應，通過文獻回顧并提取20個記錄發(fā)生點，應用MaxEnt模型分析了雙條杉天牛的潛在地理分布。結(jié)果表明：山西省南部臨汾和運城盆地為雙條杉天牛侵害最嚴重地區(qū)，中度適生區(qū)集中在中部及呂梁山西部盆地。當溫度在18～33℃范圍內(nèi)，適生度呈“S”形激增，而當NDVI處于0.2～0.6區(qū)間時，適生度急劇下降。本文的研究結(jié)果有助于理解雙條杉天牛的生物學習性，可為其防治對策提供科學依據(jù)。

關(guān)鍵詞：雙條杉天牛;林業(yè)有害生物;MaxEnt模型;適生區(qū)

中圖分類號：S7

文獻標識碼：A

DOI：10.19754/j.nyyjs.20200915025

收稿日期：2020-07-30

基金項目：山西省高等學?？萍紕?chuàng)新項目（項目編號：2019L0368）;山西農(nóng)業(yè)大學科技創(chuàng)新基金（項目編號：2017YJ18）

作者簡介：崔紹朋（1987-），男，博士，講師。研究方向：動物生態(tài)與資源利用;通訊作者吉晟男。

有害生物，是指在一定條件下，對社會經(jīng)濟和人身安全帶來危害的生物[1]，往往由于其大量群居且生存條件適宜無天敵干擾，而嚴重影響人類的生產(chǎn)生活，一般包括動物、植物、微生物和病毒等，因此有害生物的適生區(qū)分析對其預防和治理至關(guān)重要[2]。

雙條杉天牛（Semanotusbifasciatus）是一種常見的林業(yè)有害生物，隸屬于鞘翅目（Coleoptera）天牛科（Cerambycidae），是一種鉆蛀性害蟲[3，4]。該物種主要危害側(cè)柏（Platycladusorientalis）等柏科植物，初孵幼蟲蛀道穿過韌皮部造成粒狀流脂，進而在韌皮部與邊材之間鉆蛀，啃食為害，形成扁圓形不規(guī)則的蟲道，幼蟲將糞屑排在蟲道內(nèi)，流脂現(xiàn)象加劇;隨后幼蟲在木質(zhì)中向下蛀食，蛀道變粗，糞屑不排出，前蛀后填，充塞堅實，蛀道長可達20～40cm，在蛀道末端做蛹室，危害在離樹干表皮約2cm處，表現(xiàn)為上方填充滿木屑雜質(zhì)，封口處呈橢圓形[4]。被害林輕者長勢衰退，林木木質(zhì)降低，重者林木枝干中空，在外界條件作用下易折斷，嚴重時會造成整株或整枝樹木死亡[3，4]。目前，雙條杉天牛是國家明令確定的35種檢疫對象之一，同時也被列入全國林業(yè)危險性有害生物名單[5]。

MaxEnt（Maximum Entropy）模型是目前使用范圍最廣、精準度較高的生態(tài)位模型之一[6]，因其操作簡單，表現(xiàn)穩(wěn)健，被廣泛應用于動植物潛在適生區(qū)預測、環(huán)境因子檢測以及未來氣候變化影響等研究[7]。該模型以最大熵理論為基礎，其原理是在服從于一定限定條件的情況下，通過計算目標地區(qū)概率分布的最大熵來估計物種分布概率[6]。MaxEnt具有樣本量小、預測準確、只需現(xiàn)存分布點、結(jié)果易于解釋、環(huán)境變量權(quán)重可用刀切法測定等優(yōu)點，目前已經(jīng)被廣泛應用于入侵生物的潛在地理分布預測，成為最流行的預測模型之一[6-8]。

本文基于雙條杉天牛的發(fā)生地點和環(huán)境變量，運用MaxEnt模型預測雙條杉天牛的潛在適生區(qū)，旨在預測雙條杉天牛的潛在發(fā)生區(qū);探究環(huán)境影響變量對雙條杉天牛的生存影響;為防治雙條杉天牛的發(fā)生提供理論依據(jù)。

1?數(shù)據(jù)來源與處理

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

通過廣泛查閱文獻資料[3，4，9，10]，同時檢索山西省林業(yè)和草原有害生物防治檢疫局官網(wǎng)（http：//lcj.?shanxi.gov.cn/fjj/）發(fā)布的歷年有害生物發(fā)生報告，使用谷歌地球（Google Earth）提取經(jīng)緯度坐標，利用ArcGIS 10.6篩選并去除相似點，最終共保留20個分布點，位于運城、晉中、呂梁、臨汾和太原5個地級市的12個縣域內(nèi)，進一步將數(shù)據(jù)匯總保存為CSV格式以便后續(xù)分析。

考慮到生態(tài)因子的綜合性和復雜性，選用21個能反映物種生態(tài)位的環(huán)境變量（表1）。其中，植被覆蓋指數(shù)（NDVI）來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心（www.resdc.cn），19個生物氣候變量和海拔均從WorldClim數(shù)據(jù)庫（Version 2.0）下載[11]。將所有環(huán)境變量數(shù)據(jù)導入ArcGIS，裁剪并統(tǒng)一坐標系，將柵格大小重采樣為1km×1km，轉(zhuǎn)化為ASC格式儲存。

1.2?模型構(gòu)建與評估

選擇MaxEnt軟件（Version 3.4.1）模擬雙條杉天牛的適生區(qū)，創(chuàng)建響應曲線（Response curves）以分析不同變量的影響，輸出邏輯斯蒂值（Logistic）計算結(jié)果，進行10次重復，重復運行模式選擇交叉驗證（Cross-validate），隨機檢驗百分比（Random test percentage）設置為30%，其余參數(shù)選擇默認。將雙條杉天牛的適生區(qū)分為4類：非適生區(qū)，0≤P<0.2;低度適生區(qū)，0.2≤P<0.5;中度適生區(qū)，0.5≤P<0.8;高度適生區(qū)，0.8≤P≤1;其中P為物種分布概率[6-8]。

采用ROC曲線（受試者工作特征曲線）下面積（Area under ROCcurve，AUC）對模型分析結(jié)果進行精度檢驗。AUC值越大，表示與隨機分布相距越遠，環(huán)境變量與預測的物種地理分布模型之間相關(guān)性越大，即模型預測效果越好[6]。ROC曲線的評估標準為：AUC在0～0.6，失敗;0.6～0.7，較差;0.7～0.8，一般;0.8～0.9，良好;0.9～1.0，極好[8]。AUC值不受閾值影響，其結(jié)果更加客觀精準，近年來ROC曲線分析法在潛在適生區(qū)預測模型評價中的應用越來越廣泛[8]。

2?結(jié)果與分析

2.1?模型表現(xiàn)

依據(jù)雙條杉天牛適生區(qū)分布模型運行結(jié)果，AUC平均值為0.959，說明預測達到非常好的效果，表明基于MaxEnt模型預測出的分布區(qū)域不具有隨機性，模型得到的雙條杉天牛在山西省境內(nèi)潛在適生區(qū)預測結(jié)果可靠且精度較高。

2.2?山西省雙條杉天牛適生區(qū)

模型結(jié)果表明，雙條杉天牛在山西省境內(nèi)的潛在分布面積共32377km2，占全省總面積的20.7%，其中低度、中度和高度適生區(qū)分別為22119km2（14.1%）、8695km2（5.6%）和1563km2（1.0%），大部分集中在山西省中部和南部，中條山山麓向北至呂梁山南段區(qū)域尤其明顯。高適生區(qū)僅出現(xiàn)在南部臨汾和運城盆地，具體縣域主要為運城市鹽湖區(qū)、聞喜中部、河津中西部、稷山中部和新絳東部，臨汾市侯馬地區(qū)、曲沃西部和堯都中部，此外北部晉城和呂梁地區(qū)也有零星分布。中度適生區(qū)主要出現(xiàn)在山西省中部及呂梁山西部盆地區(qū)域，并被低適生區(qū)包圍。

2.3?環(huán)境影響變量

刀切法檢驗結(jié)果顯示，最熱月最高氣溫（Bio5）和植被覆蓋指數(shù)（NDVI）對雙條杉天牛在山西省的潛在分布影響最大，2個環(huán)境變量對模型的貢獻率分別為41.1%和35.7%，合計貢獻近76.8%（表1）。響應曲線分析表明，隨著氣溫的上升，雙條杉天牛適生度不斷升高，在33℃左右達到頂峰。在12～18℃范圍內(nèi)，其適生度趨于穩(wěn)定，沒有太過劇烈的起伏，但在18～33℃區(qū)間內(nèi)，適生度呈“S”形激增，隨后趨于穩(wěn)定。與此相反，隨著植被覆蓋指數(shù)的增加，雙條杉天牛適生度呈“S”形降低，NDVI處于0.2～0.6范圍時，降速最大（圖2）。

3?討論

本文研究表明，山西省南部臨汾和運城盆地是雙條杉天牛侵害最嚴重地區(qū)，中部及呂梁山西部盆地中度適生區(qū)比較集中。盡管山西省北部適生區(qū)面積小且零星分布，并且現(xiàn)有發(fā)生地均在山麓或盆地，但在全球變暖以及人類活動影響的背景下，雙條杉天牛的活動范圍可能會逐漸入侵氣候適宜的中高海拔地區(qū)，同時也應防止南部發(fā)生區(qū)向中部滲透，在相鄰發(fā)生區(qū)的盆地應特別警惕雙條杉天牛的入侵，未來城市行道樹和人造林所受危害程度可能會加劇[2，5]。

不同地區(qū)應做好相應的防范措施。在已發(fā)生嚴重雙條杉天牛蟲害地區(qū)，應開展進一步的專項檢測和防治，可采用理化和生防等方法，結(jié)合當?shù)貙嶋H情況，制定出科學合理、切實有效的措施以控制危害，并及時清除感染病害的疫木，防止向其它地區(qū)擴散蔓延[12，13]。對于潛在適生但還未發(fā)生蟲害地區(qū)，須建立完善的檢查機制，密切監(jiān)測、及時發(fā)現(xiàn)、及早消除。對于非適生區(qū)，在保持警惕的同時，可通過選育對雙條杉天牛抵抗力較強的樹木品種和改善單一的人工林結(jié)構(gòu)等方法，加強小型生境自身抗性，以減少被侵染的概率[4，5]。

雙條杉天牛在氣溫適中偏高，植被覆蓋率較低且結(jié)構(gòu)單一的環(huán)境中適生性最高，也最易成災，這與該物種的生物學特性密切相關(guān)。隨著氣溫的升高，雙條杉天牛適生度也隨之增高，趙源吉[14]通過監(jiān)測聲信號發(fā)現(xiàn)雙條杉天牛在20～28℃最為活躍，此結(jié)果與本文結(jié)論基本一致。山西省最熱月為在7—8月，氣溫一般在25～35℃，此時正是雙條杉天牛幼蟲取食量最大的時段，隨著溫度的升高，雙條杉天?；顒痈宇l繁，適生性相應增強[14]。植被覆蓋指數(shù)的升高會使雙條杉天牛的適生度急劇下降，這可能與生態(tài)系統(tǒng)抗性有關(guān)。植被覆蓋率越高，一定程度上表示生態(tài)系統(tǒng)穩(wěn)定完整，豐富的生物多樣性會顯著增強雙條杉天牛被捕食或抑制的程度。

參考文獻

[1] 國家林業(yè)局植樹造林司.中國林業(yè)檢疫性有害生物及檢疫技術(shù)操作辦法[M].北京：中國林業(yè)出版社，2005.

[2]陳燕婷，黃離，王瑞，等.近10年外來入侵昆蟲風險分析的研究現(xiàn)狀及主要進展[J].生物安全學報，2013，22（04）：217-228.

[3]高蘇嵐，許志春，弓獻詞.雙條杉天牛研究進展[J].中國森林病蟲，2007，26（03）：19-22.

[4]劉麗珍，劉晏冰，陳天林，毛亞琴.雙條杉天牛的危害與防治[J].森林病蟲通訊，1996（04）：17-18.

[5]崔永三，趙宇翔，胡學兵.我國外來林業(yè)有害生物入侵現(xiàn)狀與防控對策[J].中國森林病蟲，2009，28（03）：40-42.

[6]Phillips S J， Anderson R P， Schapire R E.Maximum entropy modeling of species geographic distributions[J]. Ecological Modelling，2006，190（3）：231-259.

[7]Phillips S J， Dudík M. Modeling of species distributions with Maxent： New extensions and a comprehensive evaluation[J]. Ecography，2008，31（2）：161-175.

[8]樊婷婷.7種中國林業(yè)檢疫性害蟲的適生區(qū)預測[D].泰安：山東農(nóng)業(yè)大學，2018.

[9]曹天文，王瑞.山西省天?？疲ㄇ食崮浚┟沎J].山西農(nóng)業(yè)大學學報（自然科學版），2007，27（04）：344-347.

[10]周維民.山西省森林植物檢疫對象分布調(diào)查[J].山西林業(yè)科技，2004（01）：33-34.

[11]FickSE， HijmansR. WorldClim 2： new 1-km spatial resolution climate surfaces for global land areas[J].International Journal of Climatology，2017，37（12）：4302-4315.

[12]李彥民，李君，郭文霞.雙條杉天牛的發(fā)生與防治[J].中國果菜，2009（04）：36.

[13]馬立芹.雙條杉天牛有效積溫測定及生物防治技術(shù)研究[D].北京：北京林業(yè)大學，2009.

[14]趙源吉，韋雪青，溫俊寶，等.雙條杉天牛（semanotusbifasciatus）幼蟲聲音探測技術(shù)初報[J].生態(tài)科學，2009，28（03）：242-246.

（責任編輯?周康）