基于MaxEnt模型的柑橘木虱在西南地區(qū)潛在分布模擬

王茹琳 王閆利 陳東東 郭翔 李慶 王明田

摘要 ：本研究應(yīng)用生態(tài)位軟件MaxEnt和地理信息系統(tǒng)軟件ArcGIS預(yù)測柑橘木虱在西南地區(qū)的潛在地理分布范圍，篩選影響柑橘木虱分布的關(guān)鍵變量。結(jié)果表明柑橘木虱在西南地區(qū)高適生區(qū)主要位于西藏東南部、四川中東部、重慶大部、貴州大部和云南中北部，總面積為52.03萬km2，占西南地區(qū)總面積的23.89%。影響柑橘木虱分布的關(guān)鍵變量為最冷月最低溫度、最冷季度平均溫度、最濕季度降水量、最暖季度降水量、最暖季度平均溫度和最暖月最高溫度。本研究不僅模擬了柑橘木虱在西南地區(qū)的潛在地理分布，更為該蟲的分布與環(huán)境變量之間關(guān)系的研究提供了理論參考。

關(guān)鍵詞 ：柑橘木虱; 西南地區(qū); MaxEnt模型; 潛在地理分布

中圖分類號：

S 431.2文獻(xiàn)標(biāo)識碼： A

DOI： 10.16688/j.zwbh.2019592

Analysis of the potential distribution of the Asian citrus psyllid， Diaphorina citri Kuwayama in Southwest China using the MaxEnt model

WANG Rulin1，2，3， WANG Yanli2， CHEN Dongdong4， GUO Xiang4， LI Qing3， WANG Mingtian5，6*

（1. Chengdu Institute of Plateau Meteorological， China Meteorology Administration， Heavy Rain and

DroughtFlood Disasters in Plateau and Basin Key Laboratory of Sichuan Province， Chengdu 610072， China;

2. Sichuan Provincial Rural Economic Information Center， Chengdu 610072， China; 3. College of Agronomy，

Sichuan Agricultural University， Chengdu 611130， China; 4. Sichuan Provincial Agrometeorological Center，

Chengdu 610072， China; 5. Sichuan Meteorological Observatory， Chengdu 610072， China; 6. WaterSaving

Agriculture in Southern Hill Area Key Laboratory of Sichuan Province， Chengdu 610066， China）

Abstract ：In this study， the ecological niche modeling software MaxEnt （the maximum entropy model） was combined with ArcGIS （geographic information system） to predict the potential geographic distribution of Diaphorina citri in China. Key environmental factors and the appropriate ranges of their values were also investigated. The results showed that the highly suitable areas of D.citri in Southwest China were mainly located in the Southeast Tibet，the middleeast of Sichuan， most of Chongqing， most of Guizhou and the northcentral part of Yunnan. The total area of the highly suitable area was 52.03×104km2， accounting for 23.89% of the total area of the five provinces.The key variables affecting the distribution of D.citri were the minimum temperature of coldest month， mean temperature of coldest quarter， precipitation of wettest quarter， precipitation of warmest quarter， mean temperature of warmest quarter， maximum temperature of warmest month. The results of this study not only established the suitable distribution areas of D.citri in Southwest China， but also would provide theoretical references for the research of the relationships between the presence of D.citri and environmental variables.

Key words ：Diaphorina citri; Southwest China; MaxEnt model; potential distribution

柑橘在我國熱帶及亞熱帶地區(qū)種植廣泛，是一種重要的經(jīng)濟(jì)作物。西南地區(qū)是我國柑橘的重要產(chǎn)區(qū)，據(jù)統(tǒng)計(jì)，2017年四川、重慶、云南和貴州四?。ㄊ校└涕俜N植面積和產(chǎn)量為64.73萬hm2和781萬t，分別占全國的25.02%和20.46%[1]。西南地區(qū)由于其獨(dú)特的地理特點(diǎn)和氣候特征，五省（市、自治區(qū)）種植品種差異明顯：四川以晚熟品種為主，2017年底種植面積約占全國晚熟柑橘總面積的60%;重慶市38個區(qū)縣中有36個生產(chǎn)柑橘，早中晚熟品種比例為25∶34∶41，栽培面積達(dá)21.33萬 hm2，年產(chǎn)量超過300萬t;云南各市（州）均有種植柑橘，以特早熟品種為主，是我國的特色柑橘生產(chǎn)基地之一;柑橘產(chǎn)業(yè)在貴州農(nóng)業(yè)經(jīng)濟(jì)中地位顯著，品種以甜橙和寬皮柑橘為主，種植面積和年產(chǎn)量均居全省果類之首;西藏自治區(qū)柑橘種植區(qū)，主要集中在雅魯藏布江河谷地帶，年產(chǎn)量約為0.07 萬t。

柑橘黃龍病是柑橘的毀滅性病害，尚無有效防治辦法，在亞洲、非洲和美洲均造成嚴(yán)重的經(jīng)濟(jì)損失，嚴(yán)重威脅柑橘產(chǎn)業(yè)的發(fā)展[2]。柑橘黃龍病在我國柑橘主要產(chǎn)區(qū)均有分布，西南地區(qū)除重慶和西藏外均有發(fā)生。在四川主要分布于涼山州、攀枝花市和宜賓市;在貴州主要分布于黔西南州、黔東南州和黔南州;在云南分布于紅河州、文山州、曲靖市、楚雄市、玉溪市、麗江市和昭通市[34]。柑橘黃龍病病原為Candidatus Liberibacter spp.，包括亞洲種、非洲種和巴西種，在中國分布的為亞洲種[5]。亞洲柑橘木虱Diaphorina citri（Kuwayama）（以下稱柑橘木虱）是柑橘黃龍病病菌的傳播介體，也是其田間自然傳播的重要途徑之一[6]。柑橘木虱以增殖式體內(nèi)循環(huán)傳播的方式傳播黃龍病，一旦獲菌即可終生帶菌并傳菌[78]。由于目前尚無監(jiān)測黃龍病早期侵染的方法，因此開展柑橘木虱潛在地理分布研究，對柑橘黃龍病的預(yù)測和防控意義重大。

對病蟲害進(jìn)行潛在地理分布預(yù)測是有害生物定量風(fēng)險(xiǎn)評估的重要研究內(nèi)容之一，常用預(yù)測軟件有CLIMEX、BIOCLIM、GARP和MaxEnt等。利用上述軟件，國內(nèi)外學(xué)者開展了大量研究，為評估病蟲害入侵、擴(kuò)散和定殖可能性并制定相應(yīng)防控策略提供了科學(xué)依據(jù)[9]。根據(jù)最大熵理論開發(fā)的MaxEnt模型是非商業(yè)軟件，2004年問世以來在生態(tài)學(xué)、物種進(jìn)化、生物安全和保護(hù)生物學(xué)等領(lǐng)域應(yīng)用廣泛。我國植保學(xué)者借助MaxEnt進(jìn)行了模擬病蟲鼠草地理分布、預(yù)測外來入侵生物定殖可能和分析病蟲害分布與氣候變化關(guān)系等方面的研究。例如趙文娟等[10]用MaxEnt預(yù)測玉米霜霉病5種病原菌在中國的適生區(qū);齊國君等[11]利用MaxEnt分析和預(yù)測了稻水象甲Lissorhoptrus oryzophilus Kuschel在中國的入侵?jǐn)U散動態(tài);唐繼洪等[12]應(yīng)用MaxEnt闡明了氣候變化背景下草地螟Loxostege sticticalis在中國越冬區(qū)的變化規(guī)律。

目前對柑橘木虱生物學(xué)特性[13]、寄主選擇[14]、防治措施[1516]、對黃龍病的傳播特性[17]等方面研究較深入，而對其潛在分布的研究較少[18]。因此本研究選用MaxEnt模型模擬柑橘木虱在西南地區(qū)的適生性，并從生態(tài)位角度篩選影響其分布的關(guān)鍵環(huán)境變量并分析二者之間的關(guān)系，期望為柑橘木虱的預(yù)測預(yù)報(bào)和有效防控提供參考。

1 材料與方法

1.1 地圖數(shù)據(jù)

本文所用地圖底圖由國家氣象信息中心提供。

1.2 環(huán)境變量

為分析柑橘木虱在西南地區(qū)的適生性，選擇19個生物氣候變量和海拔數(shù)據(jù)作為初始模型的變量。上述數(shù)據(jù)均下載自WorldClim官網(wǎng)，為1970年-2000年的平均值數(shù)據(jù)，分辨率為2.5 arcmin[19]。MaxEnt是一種基于氣候相似性原理的數(shù)學(xué)模型，用于探討物種地理分布與環(huán)境變量之間的關(guān)系[20]。環(huán)境變量的選擇是決定模擬準(zhǔn)確性的關(guān)鍵，本文參照Wang等[21]的方法并略加改進(jìn)篩選變量，具體方法為：首先使用全部環(huán)境變量建立初始模型，保留對建模貢獻(xiàn)率較高的環(huán)境變量。其次利用ArcGIS軟件提取20個變量的值，利用Pearson相關(guān)系數(shù)計(jì)算變量之間相關(guān)性，對于相關(guān)系數(shù)高于0.8的變量，根據(jù)柑橘木虱的生物學(xué)意義和初始模型中各變量的貢獻(xiàn)率篩選變量。經(jīng)過上述程序，最終保留了9個變量用于建立最終模型（表1）。

1.3 柑橘木虱的分布數(shù)據(jù)

本文通過查詢物種數(shù)據(jù)庫GBIF （https：∥www.gbif.org/）和EPPO（https：∥www.eppo.int/），共獲得了187個分布點(diǎn)。通過查閱大量公開發(fā)表的期刊論文獲得了104個分布點(diǎn)。分布數(shù)據(jù)中有經(jīng)緯度的可直接使用，只有地名（精確到鄉(xiāng)鎮(zhèn)）的分布點(diǎn)則通過百度坐標(biāo)拾取系統(tǒng)獲得經(jīng)緯度。參照江宏燕等[22]的方法去除上述291個分布點(diǎn)中的重復(fù)數(shù)據(jù)，由于環(huán)境數(shù)據(jù)分辨率為2.5 arcmin，即每一格網(wǎng)中僅保留離中心點(diǎn)最近的一個數(shù)據(jù)。通過上述程序共保留了135個柑橘木虱在中國的分布數(shù)據(jù)進(jìn)行建模（圖1）。

1.4 模型構(gòu)建和檢驗(yàn)

MaxEnt軟件（版本 3.4.1） 為開源軟件，可從美國自然歷史博物館網(wǎng)站（http：∥biodiversityinformatics.amnh.org/open_source/maxent/）免費(fèi)下載[23]。模型構(gòu)建步驟如下：1）輸入柑橘木虱的分布數(shù)據(jù)和19個環(huán)境變量建立初始模型。初始模型中設(shè)置“Random test percentage”為25，選擇“Do jackknife to measure variable importance”計(jì)算變量的貢獻(xiàn)率，設(shè)定輸出路徑，其他均設(shè)為默認(rèn)值。2） 參照1.2中方法篩選建模變量。3） 將柑橘木虱的分布數(shù)據(jù)和篩選的變量導(dǎo)入MaxEnt軟件重建模型。選擇“Random seed”，設(shè)置replicate為10， 選擇“Make pictures of predictions”“Do jackknife to measure variable importance”和創(chuàng)建響應(yīng)曲線，其他設(shè)置與初始模型一致[24]。將MaxEnt軟件輸出的ASCII文件導(dǎo)入 ArcGIS軟件中轉(zhuǎn)換為柵格文件，將其與西南地區(qū)行政區(qū)劃圖進(jìn)行疊加，得到柑橘木虱在該區(qū)域的適生區(qū)劃圖。根據(jù)IPCC關(guān)于物種分布可能性（P）的解釋并結(jié)合實(shí)際情況，將適生等級重新分為4類，并使用不同顏色顯示：高適生區(qū)（P≥0.8，紅色），中適生區(qū)（0.6≤P<0.8，橙色），低適生區(qū)（0.4≤P<0.6，黃色），不適生區(qū)（P<0.4，白色）[25]。

模型的檢驗(yàn)：ROC曲線（receiver operating characteristic curve） 是一種檢驗(yàn)生態(tài)位模型準(zhǔn)確性的有效指標(biāo)。ROC曲線以曲線下面積值即AUC值為判別指標(biāo)，AUC值取值范圍為0.5～1，越接近于1模型越準(zhǔn)確[26]。

2 結(jié)果與分析

2.1 柑橘木虱在西南地區(qū)的潛在地理分布

西南地區(qū)柑橘木虱高適生區(qū)位于西藏東南部、四川中東部、重慶大部、貴州大部和云南中北部（圖2），總面積達(dá)52.03萬km2，占西南地區(qū)總面積的23.89%。高適生區(qū)在西藏主要位于山南地區(qū)和林芝地區(qū)，總面積為2.93萬km2，占高適生區(qū)總面積的5.62%;四川省除甘孜州和阿壩州，其他市（州）均有高適生區(qū)分布，總面積為13.48萬km2;重慶市的35個區(qū)（縣）存在高適生區(qū)，總面積占全市總面積的75.69%;貴州省所有市（州）均有高適生區(qū)分布，總面積占全省總面積的89.12%;云南省的高適生區(qū)分布在德宏州、保山市、大理州、楚雄市、昆明市、文山州、紅河州和玉溪州等，總面積占高適生區(qū)總面積的29.9%。

西南地區(qū)柑橘木虱中適生區(qū)總面積為14.51萬km2，占西南地區(qū)總面積的6.67%。云南省中適生區(qū)分布最廣，集中在南部的臨滄市、普洱市、西雙版納州、紅河州和東北部的昭通市，總面積達(dá)11.17萬km2。四川、貴州、西藏和重慶的中適生區(qū)呈現(xiàn)零散分布，面積分別為1.75萬、0.71萬、0.61萬km2和0.27萬km2。

2.2 環(huán)境變量的重要性

圖3是模型輸出的各環(huán)境變量對柑橘木虱分布貢獻(xiàn)率的統(tǒng)計(jì)，結(jié)果表明10個環(huán)境變量對模型的貢獻(xiàn)率均大于0，說明此次模擬并未選入無關(guān)變量。貢獻(xiàn)率大于5%的變量為最冷月最低溫度（BIO6，35.9%）、最冷季平均溫度（BIO11，19.8%）、最暖季度降水量（BIO18，14.2%）、最濕潤季度降水量（BIO16，9.9%）、最暖季度平均溫度（BIO10，7.4%）和最暖月最高溫度（BIO5，6.2%），累積貢獻(xiàn)率高達(dá)93.4%，說明上述6個環(huán)境變量含有柑橘木虱最適生境的有效信息，是模擬柑橘木虱潛在地理分布的關(guān)鍵。

刀切法通過分別計(jì)算僅此變量、除此變量和所有變量的規(guī)則化訓(xùn)練增益反映變量對模型的重要性。圖4是刀切法輸出的結(jié)果：最冷月最低溫度的訓(xùn)練增益最高，超過1.7，是最關(guān)鍵的關(guān)鍵變量。最冷季度平均溫度、最濕季度降水量和最暖季度降水量的訓(xùn)練增益均超過1.6，均為較重要的環(huán)境變量。最干月降水量和降水季節(jié)性變異系數(shù)的訓(xùn)練增益較低。

2.3 環(huán)境變量對柑橘木虱的影響

環(huán)境變量與物種存在概率之間的響應(yīng)曲線可反映二者間的關(guān)系。由圖5看出，當(dāng)最冷月最低溫度低于-3.28℃時，柑橘木虱的存在概率低于0.33，隨著溫度的升高存在概率逐漸升高，高于22℃達(dá)到最高值。最冷季度平均溫度響應(yīng)曲線變化趨勢為：溫度在6.03～17.88℃時，柑橘木虱存在概率高于0.33，10.26℃時存在概率最高。最濕季度降水量低于548.66 mm時柑橘木虱的存在概率低于0.33，隨著降雨量的增高存在概率迅速升高，在619.02 mm時達(dá)到峰值，當(dāng)降雨量超過1 189.75 mm時存在概率再次降到0.5之下。其他4種關(guān)鍵變量響應(yīng)曲線的適宜范圍分別為：最暖季度平均溫度24.46～34.27℃，最暖季度降水量503.73～1 533.58 mm，最暖月最高溫度28.6～40.91℃和溫度季節(jié)性變化標(biāo)準(zhǔn)差56.83～818.03（圖5）。

2.4 模型的評價

生態(tài)位模型常用的評價指標(biāo)包括準(zhǔn)確度、敏感度、特異度、Kappa指標(biāo)、TSS指標(biāo)和ROC曲線等。ROC曲線由于不受閾值影響而被認(rèn)為是目前最好的評價指標(biāo)之一。MaxEnt軟件可直接繪制ROC曲線并計(jì)算AUC值，所以此方法的應(yīng)用最廣泛。例如，Wang等[27]利用MaxEnt分析馬鈴薯甲蟲在全球分布范圍時使用ROC曲線對模擬精度進(jìn)行了評價;唐繼洪等[12]應(yīng)用MaxEnt模型預(yù)測了草地螟在我國的越冬區(qū)，并選用ROC曲線評價了預(yù)測的準(zhǔn)確性。因此，本研究選用ROC曲線來評價預(yù)測精度。圖6a為初始模型的ROC曲線，訓(xùn)練數(shù)據(jù)和測試數(shù)據(jù)的AUC值分別為0.966和0.956。 由判別標(biāo)準(zhǔn)可知，初始模型的準(zhǔn)確性級別為極好。圖6b為最終模型的ROC曲線，10次重復(fù)的AUC平均值為0.988，說明本研究的模擬效果“很好”，研究結(jié)果可用于研究柑橘木虱在我國西南地區(qū)的潛在分布。

3 討論

3.1 柑橘木虱在西南地區(qū)潛在地理分布模擬

本研究基于MaxEnt模型模擬了西南地區(qū)柑橘木虱的潛在地理分布，結(jié)果表明高適生區(qū)主要位于西藏東南部、四川中東部、重慶大部、貴州大部和云南中北部，與汪善勤等[18]的模擬結(jié)果相比，本研究預(yù)測的高適生區(qū)分布更廣，這是由于不同的分布數(shù)據(jù)、不同的環(huán)境變量和不同的模擬理論所致。MaxEnt基于最大熵理論量化物種分布，此法對未知分布信息不加限制，較其他同類模型優(yōu)勢更明顯，被國內(nèi)外學(xué)者廣泛應(yīng)用。Elith等[20]對比了多種生態(tài)位模型的預(yù)測能力，結(jié)論是MaxEnt在16種模型中預(yù)測能力最高;張海濤等[28]比較了4種常用的生態(tài)位模型的預(yù)測準(zhǔn)確性，結(jié)果表明MaxEnt模型的預(yù)測效果更好。在植保領(lǐng)域我國學(xué)者應(yīng)用該模型對地中海實(shí)蠅Ceratitis capitata[29]、中國飛蝗Locusta migratoria[30]、菊方翅網(wǎng)蝽Corythucha marmorata[31]、草地貪夜蛾Spodoptera frugiperda[32]、刺槐突瓣細(xì)蛾Chrysaster ostensackenella[33]、中對長小蠹Platypus parallelus[34]、紅脂大小蠹Dendroctonus valens[35]、木瓜秀粉蚧Paracoccus marginatus[36]和斑翅果蠅Drosophila suzukii[37]等病蟲害成功地進(jìn)行了潛在地理分布預(yù)測?；谏鲜鲈?，本文選擇MaxEnt作為模擬軟件。

全國農(nóng)業(yè)技術(shù)推廣服務(wù)中心2014年專項(xiàng)調(diào)查顯示，柑橘木虱在四川、貴州、云南的發(fā)生北界分別為涼山州雷波縣（28°24′N，103°59′E），黔東南州從江縣（25°52′N， 108°06′E）和麗江市華坪縣（26°38′N，101°14′E）[4]。2017年云南和四川植物檢疫部門在兩省交界的昭通市綏江縣和宜賓市屏山縣發(fā)現(xiàn)柑橘木虱[6]，說明柑橘木虱有逐步向北擴(kuò)散的趨勢。汪善勤等應(yīng)用生物氣候模型CLIMEX比較了1990年到2050年柑橘木虱適生區(qū)的發(fā)展趨勢，發(fā)現(xiàn)柑橘木虱在四川和云南的適生區(qū)均有北移趨勢[18]。本研究結(jié)果指出，四川盆地中東部、云南中北部和貴州大部為柑橘木虱的高適生區(qū)，表明上述地區(qū)具備適宜柑橘木虱定殖的氣候條件，未來有向上述地區(qū)擴(kuò)散的可能，與上述研究結(jié)果一致。因此預(yù)測為高適生區(qū)且柑橘木虱已發(fā)生的地區(qū)，應(yīng)著重開展防治工作;尚未發(fā)現(xiàn)柑橘木虱而模擬為高適生區(qū)的地區(qū)，農(nóng)業(yè)部門應(yīng)高度重視，加強(qiáng)調(diào)查和檢疫工作，防止柑橘木虱傳入。2018年四川省農(nóng)業(yè)廳開展的柑橘木虱專項(xiàng)調(diào)查顯示，柑橘木虱尚未擴(kuò)散到成都平原。考慮到柑橘黃龍病對柑橘產(chǎn)業(yè)的嚴(yán)重威脅，四川省已在宜賓地區(qū)著手建立柑橘木虱阻截帶，防止其沿長江流域蔓延擴(kuò)散。本研究顯示，重慶市35個區(qū)（縣）有柑橘木虱高適生區(qū)分布，總面積占全市面積75.71%。且重慶西南部與云南昭通和四川宜賓交界，柑橘木虱傳入風(fēng)險(xiǎn)極高，因此該市雖尚未發(fā)現(xiàn)柑橘木虱，但仍不可放松警惕。

3.2 環(huán)境變量對柑橘木虱潛在地理分布的影響

病蟲害的發(fā)生、繁殖和傳播與環(huán)境、寄主、農(nóng)業(yè)生態(tài)系統(tǒng)和管理水平等緊密相關(guān)。氣象要素是環(huán)境變量的重要組成，在其他因素相對固定的情況下，氣象因素可能會成為決定病蟲害大尺度分布的決定性因素[3839]。研究表明，溫度、濕度和大氣壓力等因素對柑橘木虱的潛在地理分布影響顯著[40]。

本研究使用刀切法測定了變量的重要性，最冷月最低溫度是最重要的環(huán)境變量，說明該變量與柑橘木虱潛在地理分布的關(guān)系最為緊密。響應(yīng)曲線表明當(dāng)最冷月最低溫度低于-3.28℃時，柑橘木虱的存在概率極低，說明極端低溫是限制柑橘木虱分布的關(guān)鍵因子。在印度，Atwal等[41]研究發(fā)現(xiàn)0℃以下低溫對壓低柑橘木虱的種群數(shù)量非常關(guān)鍵，Yang等[42]在中國也發(fā)現(xiàn)了同樣的規(guī)律。白先進(jìn)等[43]的研究表明，日平均低溫是影響柑橘木虱分布的主要因素;Hall等[44]研究指出冬季低溫是限制柑橘木虱種群增長、地理分布和傳播的主要?dú)庀笠蜃印Ｉ鲜鲅芯拷Y(jié)果均表明冬季低溫對于柑橘木虱的關(guān)鍵作用。研究表明，高溫對昆蟲產(chǎn)生的傷害具有積累效應(yīng)，高溫條件持續(xù)一段時間，會導(dǎo)致昆蟲體內(nèi)水分含量明顯下降造成昆蟲死亡[4546]。本研究結(jié)論指出最暖月最高溫度適合柑橘木虱發(fā)生適宜值范圍為28.6～40.91℃，溫度高于40.91℃時，柑橘木虱存在的可能性很低，這與早期一些研究結(jié)論一致。鄺凡等[47]研究了高溫下柑橘木虱死亡率，指出40℃以上高溫條件下，死亡時間隨溫度升高而降低，高溫是高溫季節(jié)木虱種群消退的主要原因。Hall等[44]研究證實(shí)溫度高于41.6℃時，柑橘木虱將無法產(chǎn)卵。上述結(jié)論說明40℃高溫是制約柑橘木虱分布的重要因素。NaroueiKhandan等[48]研究表明最暖季度的平均溫度高于33℃時，柑橘木虱的存在受到限制，本文中最暖季度平均溫度適宜的上限溫度為34.27℃，這可能是因?yàn)楦涕倌臼粝x和卵對高溫更為敏感[49]。MaxEnt軟件繪制的響應(yīng)曲線反映的是單一環(huán)境變量對柑橘存在概率的影響，忽略了變量間的相互作用，具有一定的局限性。因此本文所得結(jié)論不能完全闡述環(huán)境變量對柑橘木虱的作用，但可作為判斷二者關(guān)系的度量。

不同寄主植物對柑橘木虱發(fā)育和繁殖的影響不同。任素麗等[50]測定了柑橘木虱在柚、酸橘、黃皮、九里香和砂糖橘上的存活情況，結(jié)果表明，若蟲在柚上的存活率最高，在黃皮上最低。成蟲在柚上壽命顯著長于黃皮上的壽命。Tsai等[51]研究了柑橘木虱在多個柑橘品種上的存活率，發(fā)現(xiàn)在西柚上存活率最高，酸橙上的最低，不同品種柑橘上的存活率差異顯著。陳建利等[52]研究了柑橘木虱對13個柑橘品種的取食偏好，結(jié)果表明取食‘紐荷爾臍橙‘福橘和‘佩奇甜橙的數(shù)量顯著高于其他柑橘品種。上述研究證明了寄主對柑橘木虱分布的重要性，分析二者之間的關(guān)系有助于提高模擬的準(zhǔn)確性。本文未對柑橘木虱各種寄主在西南地區(qū)的適生性進(jìn)行分析，會對預(yù)測結(jié)果造成一定影響。因此進(jìn)一步分析時，需考慮寄主特別是柑橘的種植分布，使預(yù)測結(jié)果更有意義。

參考文獻(xiàn)

[1] 吳園. 中國柑橘種植業(yè)技術(shù)效率評估及影響因素分析[J]. 中國農(nóng)業(yè)資源與區(qū)劃， 2018， 39（9）： 94102.

[2] 程春振， 曾繼吾， 鐘云， 等. 柑橘黃龍病研究進(jìn)展[J]. 園藝學(xué)報(bào)， 2018， 40（9）： 16561668.

[3] 范國成， 劉波， 吳如健， 等. 中國柑橘黃龍病研究30年[J]. 福建農(nóng)業(yè)學(xué)報(bào)， 2009， 24（2）： 183190.

[4] 王曉亮， 李瀟楠， 馮曉東， 等. 柑橘黃龍病與柑橘木虱在我國發(fā)生情況調(diào)查[J]. 植物檢疫， 2016， 30（2）： 44.

[5] 胡文召， 周常勇. 柑橘黃龍病病原研究進(jìn)展[J]. 植物保護(hù)， 2010， 36（3）： 3033.

[6] 江宏燕， 吳豐年， 王妍晶， 等. 亞洲柑橘木虱的起源、分布和擴(kuò)散能力研究進(jìn)展[J]. 環(huán)境昆蟲學(xué)報(bào)， 2018， 40（5）： 10141020.

[7] 陳麗芬， 徐昭煥， 王建國. 柑橘木虱的研究進(jìn)展[J]. 貴州農(nóng)業(yè)科學(xué)， 2016， 44（6）： 4247.

[8] 任素麗， 郭長飛， 歐達(dá)， 等. 黃龍病病菌在柑橘木虱體內(nèi)的分布及感染動態(tài)[J]. 應(yīng)用昆蟲學(xué)報(bào)， 2018， 55（4）： 595601.

[9] 李志紅，秦譽(yù)嘉.有害生物風(fēng)險(xiǎn)分析定量評估模型及其比較[J].植物保護(hù)，2018，44（5）：134145.

[10]趙文娟，陳林，丁克堅(jiān)， 等.利用MAXENT預(yù)測玉米霜霉病在中國的適生區(qū)[J].植物保護(hù)， 2009， 35（2）：3238.

[11]齊國君， 高燕， 黃德超， 等. 基于MAXENT的稻水象甲在中國的入侵?jǐn)U散動態(tài)及適生性分析[J]. 植物保護(hù)學(xué)報(bào)， 2012， 39（2）： 129136.

[12]唐繼洪， 程云霞， 羅禮智， 等. 基于Maxent模型的不同氣候變化情景下我國草地螟越冬區(qū)預(yù)測[J]. 生態(tài)學(xué)報(bào)， 2017， 37（14）： 48524863.

[13]余繼華， 黃振東， 張敏榮， 等. 亞洲柑橘木虱的生物學(xué)研究進(jìn)展及防治[J]. 植物檢疫， 2018， 32（5）： 813.

[14]胡燕， 王雪峰， 周常勇.柑橘黃龍病菌亞洲種、蟲媒及植物寄主互作研究進(jìn)展[J]. 園藝學(xué)報(bào)， 2016， 43（9）： 16881698.

[15]余繼華， 汪恩國. 柑橘黃龍病發(fā)生危害與防治指標(biāo)研究[J]. 浙江農(nóng)業(yè)學(xué)報(bào)， 2009， 21（4）： 370374.

[16]姚廷山， 周彥， 周常勇. 亞洲柑橘木虱的發(fā)生與防治研究進(jìn)展[J]. 果樹學(xué)報(bào)， 2018， 35（11）： 14131421.

[17]馬曉芳， 張學(xué)潮， 徐海君. 柑橘木虱兩地理種群的內(nèi)共生菌群落組成及傳菌能力[J]. 昆蟲學(xué)報(bào)， 2012， 55（10）： 11491153.

[18]汪善勤， 肖云麗， 張宏宇. 我國柑橘木虱潛在適生區(qū)分布及趨勢分析[J]. 應(yīng)用昆蟲學(xué)報(bào)， 2015， 52（5）： 11401148.

[19]FICK S E， HIJMANS R J. WorldClim 2： new 1km spatial resolution climate surfaces for global land areas [J]. International Journal of Climatology， 2017， 37（12）： 43024315.

[20]ELITH J， GRAHAM C H， ANDERSON R P， et al. Novel methods improve prediction of species distributions from occurrence data [J]. Ecography， 2006， 29（2）： 129151.

[21]WANG Rulin， LI Qing， HE Shisong， et al. Modeling and mapping the current and future distribution of Pseudomonas syringae pv. actinidiae under climate change in China [J/OL]. PLoS ONE， 2018， 13（2）： e0192153. DOI：10.1371/journal.pone.0192153.

[22]江宏燕，吳豐年，王妍晶，等.亞洲柑橘木虱的起源、分布和擴(kuò)散能力研究進(jìn)展[J].環(huán)境昆蟲報(bào)，2018，40（5）：10141020.

[23]PHILLIPS S J， ANDERSON R P， SCHAPIRE R E. Maximum entropy modeling of species geographic distributions [J]. Ecological Modelling， 2006， 190（34）： 231259.

[24]KUMAR S， STOHLGREN T J. Maxent modeling for predicting suitable habitat for threatened and endangered tree Canacomyrica monticola in New Caledonia [J]. Journal of Ecology and Natural Environment， 2009， 1（4）： 9498.

[25]ZHU Yuanjun， WEI Wei， LI Hao， et al. Modelling the potential distribution and shifts of three varieties of Stipa tianschanica in the eastern Eurasian Steppe under multiple climate change scenarios [J]. Global Ecology and Conservation， 2018， 16： e00501.

[26]王運(yùn)生， 謝丙炎， 萬方浩， 等. ROC曲線分析在評價入侵物種分布模型中的應(yīng)用[J]. 生物多樣性， 2007， 15（4）：365372.

[27]WANG Cong， HAWTHORNE D， Qin Yujia， et al. Impact of climate and host availability on future distribution of Colorado potato beetle [J/OL]. Scientific Reports， 2017， 7（1）： 4489.DOI：10.1038/s41598017046077.

[28]張海濤， 羅渡， 牟希東， 等. 應(yīng)用多個生態(tài)位模型預(yù)測福壽螺在中國的潛在適生區(qū)[J]. 應(yīng)用生態(tài)學(xué)報(bào)， 2016， 27（4）：12771284.

[29]孫佩珊，姜帆，張祥林，等.地中海實(shí)蠅入侵中國的風(fēng)險(xiǎn)評估[J].植物保護(hù)學(xué)報(bào)，2017，44（3）：436444.

[30]張杰，張旸，趙振勇，等.中國飛蝗（Locusta migratoria）災(zāi)害地理分布模擬及其生物氣候因子分析[J].干旱區(qū)地理，2019，42（3）：590598.

[31]王志華，于靜亞，沈錦，等.入侵害蟲菊方翅網(wǎng)蝽在中國的潛在分布預(yù)測[J].環(huán)境昆蟲學(xué)報(bào)，2019，41（3）：626633.

[32]秦譽(yù)嘉，藍(lán)帥，趙紫華，等.遷飛性害蟲草地貪夜蛾在我國的潛在地理分布[J].植物保護(hù)，2019，45（4）：4347.

[33]樊婷婷，高尚坤，孟凡玲，等.外來入侵新害蟲刺槐突瓣細(xì)蛾在中國的適生區(qū)預(yù)測[J].林業(yè)科學(xué)，2019，55（6）：8695.

[34]范靖宇，吳哿，朱耿平，等.檢疫性害蟲中對長小蠹在中國的適生區(qū)[J].林業(yè)科學(xué)，2019，55（6）：8185.

[35]崔驍芃，石娟，王海香，等.基于MaxEnt模型的紅脂大小蠹在中國適生區(qū)的預(yù)測[J].植物保護(hù)學(xué)報(bào)，2019，46（4）：925926.

[36]宋子驕，秦譽(yù)嘉，馬福歡，等.基于MaxEnt的木瓜秀粉蚧適生性分析[J].植物檢疫，2019，33（5）：7378.

[37]孫宏禹，李志紅，楊普云，等.基于最大熵模型的斑翅果蠅適生性研究[J].中國植保導(dǎo)刊，2016，36（12）：6673.

[38]MONTEITH J L. Agricultural meteorology： evolution and application [J]. Agricultural and Forest Meteorology， 2000， 103（1）： 59.

[39]QIN Jianglin， YANG Xiuhao， YANG Zhongwu， et al. New technology for using meteorological information in forest insect pest forecast and warning systems [J]. Pest Management Science， 2017， 73（12）： 25092518.

[40]MARTINI X， STELINSKI L L. Influence of abiotic factors on flight initiation by Asian citrus psyllid （Hemiptera： Liviidae） [J]. Environmental Entomology， 2017， 46（2）： 369375.

[41]ATWAL A S， CHAUDHARY J P， RAMZAN M. Studies on the development and field population of citrus psylla， Diaphorina citri Kuwayama （Psyllidae： Homoptera） [J]. Journal of ResearchPunjab Agricultural University， 1970， 7（3）： 333338.

[42]YANG Yueping， HUANG Mingdu， ANDREW G， et al. Distribution， biology， ecology and control of the psyllid Diaphorina citri Kuwayama， a major pest of citrus： a status report for China [J]. International Journal of Pest Management， 2006， 52（4）： 343352.

[43]白先進(jìn)， 鄧崇嶺， 陸國保， 等. 柑桔木虱耐寒性調(diào)查研究[J]. 中國南方果樹， 2008， 37（6）： 2224.

[44]HALL D G， WENNINGER E J， HENTZ M G. Temperature studies with the Asian citrus psyllid， Diaphorina citri： Cold hardiness and temperature thresholds for oviposition [J/OL]. Journal of Insect Science， 2011， 11（1）： 83.DOI：10.1673/031.011.8301.

[45]鐘景輝. 松突圓蚧花角蚜小蜂成蟲的耐熱性[J]. 福建農(nóng)林大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2010， 39（5）： 465470.

[46]張會娜， 呂建華， 亞川川， 等. 高溫處理對玉米象成蟲死亡率和水分含量的影響[J]. 河南工業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2014， 35（6）： 4548.

[47]鄺凡， 林夢琪， 袁水秀， 等. 高溫對亞洲柑橘木虱死亡率和活動行為的影響[J]. 南方農(nóng)業(yè)學(xué)報(bào)， 2017， 48（9）： 16001604.

[48]NAROUEIKHANDAN H A， HALBERT S E，WORNER S P， et al. Global climate suitability of citrus huanglongbing and its vector， the Asian citrus psyllid， using two correlative species distribution modeling approaches， with emphasis on the USA [J]. European Journal of Plant Pathology， 2016， 144（3）： 655670.

[49]LIU Yinghong， TSAI J H. Effects of temperature on biology and life table parameters of the Asian citrus psyllid， Diaphorina citri Kuwayama Homoptera： Psyllidae [J]. Annals of Applied Biology， 2000， 137（3）： 201206.

[50]任素麗， 歐達(dá)， 張利荷， 等. 不同寄主植物對柑橘木虱發(fā)育和繁殖的影響[J]. 應(yīng)用昆蟲學(xué)報(bào)， 2018， 55（4）： 5661.

[51]TSAI J H， LIU Yinghong. Biology of Diaphorina citri （Homoptera： Psyllidae） on four host plants [J]. Journal of Economic Entomology， 2000， 93（6）：17211725.

[52]陳建利， 阮傳清， 劉波， 等. 柑橘木虱對柑橘不同品種的趨性[J]. 福建農(nóng)業(yè)學(xué)報(bào)， 2011， 26（2）： 280283.

（責(zé)任編輯：田 喆）