2050年3種氣候模式下稻縱卷葉螟對(duì)我國(guó)風(fēng)險(xiǎn)分析

王閆利,王茹琳,姜 淦,沈沾紅,林 姍

(1.四川省農(nóng)村經(jīng)濟(jì)綜合信息中心，四川 成都 610072;2.高原與盆地暴雨旱澇災(zāi)害四川省重點(diǎn)實(shí)驗(yàn)室, 四川 成都 610072)

?

2050年3種氣候模式下稻縱卷葉螟對(duì)我國(guó)風(fēng)險(xiǎn)分析

王閆利1,2,王茹琳1,姜 淦1,沈沾紅1,林 姍1

(1.四川省農(nóng)村經(jīng)濟(jì)綜合信息中心，四川 成都 610072;2.高原與盆地暴雨旱澇災(zāi)害四川省重點(diǎn)實(shí)驗(yàn)室, 四川 成都 610072)

摘要：采用最大熵生態(tài)位模型(MaxEnt)作為分析軟件，結(jié)合ArcGIS，預(yù)測(cè)了2050年時(shí)間段稻縱卷葉螟在我國(guó)的風(fēng)險(xiǎn)區(qū)變化。結(jié)果表明:稻縱卷葉螟在中國(guó)的極高風(fēng)險(xiǎn)區(qū)為上海、浙江大部、江蘇南部、安徽大部、湖北中東部、湖南東部、四川西南部、廣西中部等地；高風(fēng)險(xiǎn)區(qū)為貴州、廣東、福建、江西、重慶、廣西大部。2050年，在3種排放情景下，稻縱卷葉螟在中國(guó)的極高風(fēng)險(xiǎn)區(qū)面積均將減少，B2a情景下高風(fēng)險(xiǎn)區(qū)面積將增加。

關(guān)鍵詞：稻縱卷葉螟；MaxEnt；風(fēng)險(xiǎn)分析；預(yù)測(cè)

0引言

稻縱卷葉螟是一種熱帶和亞熱帶的重要水稻害蟲(chóng)，隸屬于昆蟲(chóng)綱(Isecta)、鱗翅目(Lepidoptera)、螟蛾科(Pyralidae)[1]。稻縱卷葉螟廣泛分布于亞洲、澳大利亞、東非等國(guó)家和地區(qū)。在我國(guó)該蟲(chóng)主要危害南方稻區(qū)和長(zhǎng)江中下游地區(qū)[2]，而在20世紀(jì)90年代后半期，稻縱卷葉螟發(fā)生范圍逐漸擴(kuò)增且強(qiáng)度明顯加重，目前已成為我國(guó)水稻主產(chǎn)區(qū)上常發(fā)性重要害蟲(chóng)。21世紀(jì)以來(lái)曾多次造成我國(guó)水稻大幅度減產(chǎn)，造成嚴(yán)重的經(jīng)濟(jì)損失[3]。

生態(tài)位模型的工作原理是：利用物種的地理分布數(shù)據(jù)(分布點(diǎn)經(jīng)、緯度信息)和環(huán)境數(shù)據(jù)，以不同生物所需的特殊生存環(huán)境(生態(tài)位要求)為依據(jù)，利用特定數(shù)學(xué)算法，模擬或歸納特定物種的生態(tài)位需求，將其投射到目標(biāo)地區(qū)，因此而獲得的結(jié)果即為目標(biāo)物種的豐富度、出現(xiàn)概率或生境適宜度[4-7]。近年來(lái)，多種生態(tài)位模型(MaxEnt、GARP、Bioclim、ENFA和Domain等)在物種潛在地理分布預(yù)測(cè)方面取得了良好的效果，同時(shí)在生態(tài)學(xué)、考古學(xué)、生物地理學(xué)、進(jìn)化和保護(hù)生物學(xué)等領(lǐng)域得到了廣泛應(yīng)用[8-13]。MaxEnt模型是一種基于最大熵原理，在模擬計(jì)算過(guò)程中，結(jié)合只出現(xiàn)型數(shù)據(jù)模型和生態(tài)位原理的對(duì)物種地理分布進(jìn)行預(yù)測(cè)的方案，是目前預(yù)測(cè)物種潛在分布的最新的生態(tài)位模型。MaxEnt模型同其他生態(tài)位模型相比具有僅需少量數(shù)據(jù)即可進(jìn)行預(yù)測(cè)且較其它方法更精確的優(yōu)點(diǎn)[14-16]。

近年來(lái)，隨著全球氣候變暖的不斷加劇，出現(xiàn)了病蟲(chóng)害繁殖代數(shù)增多、越冬死亡率降低和發(fā)生期提前等一系列問(wèn)題，因此有必要采取不同的研究方法評(píng)估氣候變化對(duì)病蟲(chóng)害未來(lái)分布的影響。預(yù)測(cè)氣候變化背景下的病蟲(chóng)害適生區(qū)的變化，對(duì)及時(shí)掌握病蟲(chóng)害的發(fā)生發(fā)展規(guī)律、建立健全病蟲(chóng)害防控機(jī)制具有重要參考意義。本研究以水稻的重要害蟲(chóng)稻縱卷葉螟為研究對(duì)象，借助MaxEnt模型和ArcGIS評(píng)估稻縱卷葉螟“當(dāng)前時(shí)段”的分布狀態(tài)和“2050年時(shí)段”的變化，得到不同時(shí)期稻縱卷葉螟種群分布面積的變化，以此揭示氣候變化對(duì)該蟲(chóng)害帶來(lái)的深刻影響，以期為制定合理的管理和保護(hù)策略提供科學(xué)依據(jù)。

1材料與方法

1.1稻縱卷葉螟地理分布信息的獲取

實(shí)地調(diào)查、國(guó)內(nèi)外發(fā)表論文、查詢數(shù)據(jù)庫(kù)及標(biāo)本館記錄為獲取物種地理分布數(shù)據(jù)的4種主要途徑[17-18]。在本研究中，通過(guò)全球物種多樣性信息庫(kù)(GBIF，http://www.gbif.org/)、高校教學(xué)標(biāo)本資源共享平臺(tái)(http://mnh.scu.edu.cn/)和檢索國(guó)內(nèi)外公開(kāi)發(fā)表的有關(guān)稻縱卷葉螟的論文[19,33]，共獲得稻縱卷葉螟在全球的分布點(diǎn)238個(gè)。在238個(gè)分布點(diǎn)中，有準(zhǔn)確經(jīng)緯度信息的直接使用；對(duì)于只有地名的分布點(diǎn)，通過(guò)地名查詢系統(tǒng)GeoName(http://www.geonames.org/)查詢具體經(jīng)緯度。MaxEnt軟件運(yùn)行時(shí)，將稻縱卷葉螟的地理信息數(shù)據(jù)保存為“*.CSV”格式文件，保存順序?yàn)椋何锓N名(species)、經(jīng)度(dd long)和緯度(dd lat)[34-35]。

1.2環(huán)境層數(shù)據(jù)和地圖數(shù)據(jù)的獲取與處理

1.2.1環(huán)境層數(shù)據(jù)的獲取與處理環(huán)境數(shù)據(jù)包括年平均氣溫、最暖月最高溫度、最濕月降水量等19個(gè)全球環(huán)境生物氣候變量，變量代碼、描述以及指標(biāo)代表的意義見(jiàn)表1，數(shù)據(jù)空間分辨率為5 arcmin，“當(dāng)前時(shí)段”(1950～2000年)氣候數(shù)據(jù)可在WORLDCLIM網(wǎng)站(http://www.worldclim.org/)查詢并免費(fèi)下載[36]。IPCC AR4針對(duì)未來(lái)的氣候變化，提供了多種氣候模式，本文選擇了3種具有代表性的氣候模式，分別為：A1b模式(能源需求平衡)、A2a(能源需求較高)和B2a(能源需求較低)，時(shí)間區(qū)間選擇為2050s(2041～2050年)，該時(shí)段上述3種氣候模式數(shù)據(jù)可從國(guó)際熱帶農(nóng)業(yè)中心網(wǎng)站的未來(lái)氣候數(shù)據(jù)庫(kù)免費(fèi)下載。

1.2.2ArcGIS所需地圖數(shù)據(jù)的獲取與處理MaxEnt軟件輸出的結(jié)果需借助ArcGIS進(jìn)行分級(jí)顯示，本研究中中國(guó)分析的底圖(1∶400萬(wàn))下載自國(guó)家基礎(chǔ)地理信息系統(tǒng)官方網(wǎng)站，世界行政區(qū)劃圖下載自Blue Marble Geographics[37]。

1.3MaxEnt模型預(yù)測(cè)與檢驗(yàn)

參照程軍[38]、郭彥龍[39]、張熙驁等[40]從全球238個(gè)稻縱卷葉螟分布點(diǎn)中隨機(jī)選取25%作為MaxEnt模擬的測(cè)試集數(shù)據(jù)，剩余的75%分布點(diǎn)則作為訓(xùn)練集，輸出格式為“ASCII”文件，其它設(shè)定均選擇“默認(rèn)”。

研究證實(shí)[41-43]，ROC接受曲線是目前最為廣泛使用的、針對(duì)物種分布模擬結(jié)果準(zhǔn)確性驗(yàn)證最為有效的指標(biāo)。ROC曲線以1-特異度(假陽(yáng)性率)為橫坐標(biāo)，1-遺漏率(真陽(yáng)性率)為縱坐標(biāo)繪制而成，通過(guò)計(jì)算曲線下方面積(area under curve，AUC)而獲得評(píng)價(jià)指標(biāo)。評(píng)價(jià)標(biāo)準(zhǔn)為：AUC>0.9為“極好”，0.8

表1　用于物種分布模擬的 19個(gè)生物氣象因子

2結(jié)果與分析

2.1MaxEnt模擬結(jié)果驗(yàn)證

表2為“當(dāng)前時(shí)段”和“2050年時(shí)段”MaxEnt模型預(yù)測(cè)的AUC值，“當(dāng)前時(shí)段”模型AUC值為0.982，預(yù)測(cè)效果為“極好”，說(shuō)明模擬結(jié)果可用于稻縱卷葉螟風(fēng)險(xiǎn)等級(jí)預(yù)測(cè)。“2050年時(shí)段”A1b、A2a和B2a 3種情景下，AUC值分別為0.982、0.985和0.984，預(yù)測(cè)效果均為“極好”，表明3種模擬結(jié)果均適用于本研究。

表2　稻縱卷葉螟風(fēng)險(xiǎn)等級(jí)預(yù)測(cè)模型的AUC值

2.2稻縱卷葉螟在中國(guó)的風(fēng)險(xiǎn)等級(jí)分析

根據(jù)稻縱卷葉螟在全球的分布記錄數(shù)據(jù)，利用MaxEnt生態(tài)位模型和ArcGIS軟件，稻縱卷葉螟在中國(guó)的風(fēng)險(xiǎn)等級(jí)區(qū)劃結(jié)果如圖1所示。采用“均分等級(jí)法”對(duì)預(yù)測(cè)結(jié)果進(jìn)行分類，按風(fēng)險(xiǎn)等級(jí)從低到高順序依次分為5個(gè)等級(jí)。在當(dāng)前氣候條件下，稻縱卷葉螟在中國(guó)的適宜生長(zhǎng)區(qū)(極高風(fēng)險(xiǎn)區(qū))為上海、浙江大部、江蘇南部、安徽大部、湖北中東部、湖南東部、四川西南部、廣西中部等地，總面積約78.41萬(wàn)km2。高風(fēng)險(xiǎn)區(qū)為貴州、廣東、福建、江西、重慶、廣西大部，總面積約80.2萬(wàn)km2。

浙江金融職業(yè)學(xué)院智慧圖書(shū)館，空間預(yù)約、電腦登錄、活動(dòng)報(bào)名與考勤都已經(jīng)實(shí)現(xiàn)了人臉識(shí)別，是全國(guó)高校首家將人臉識(shí)別技術(shù)應(yīng)用到圖書(shū)館所有業(yè)務(wù)流程，實(shí)現(xiàn)一站式服務(wù)的圖書(shū)館。

2.32050年3種情景下稻縱卷葉螟在中國(guó)的適生性分析

2050年，在3種排放情景下，稻縱卷葉螟在中國(guó)的極高風(fēng)險(xiǎn)區(qū)面積均將減少，A1b情景下分布范圍以目前為中心向內(nèi)減少，總面積約為69.5萬(wàn)km2；A2a和B2a情景下極高風(fēng)險(xiǎn)區(qū)面積均將明顯減少，總面積分別為58.93萬(wàn)km2和58.24萬(wàn)km2，其中江西、浙江、湖北和四川等地減少較為明顯，B2a情景下高風(fēng)險(xiǎn)區(qū)面積將增加，總面積約為85.66萬(wàn)km2(圖2)。

2.4稻縱卷葉螟在我國(guó)各省區(qū)的氣候風(fēng)險(xiǎn)綜合指數(shù)分析

表3為稻縱卷葉螟在中國(guó)各省區(qū)的氣候風(fēng)險(xiǎn)等級(jí)綜合指數(shù)(CRRI)，計(jì)算方法參照孫兵等[44]的研究方法，公式如下：

上式中 ，Bi代表i等級(jí)的氣候風(fēng)險(xiǎn)指數(shù)，APi是Bi地區(qū)對(duì)應(yīng)的氣候風(fēng)險(xiǎn)指數(shù)等級(jí)所對(duì)應(yīng)的面積百分比。

結(jié)果表明：在“當(dāng)前時(shí)段”氣候條件下，在我國(guó)范圍內(nèi)，從氣候適應(yīng)性綜合指數(shù)分析，稻縱卷葉螟在上海、浙江、香港、江西、安徽和湖南等省市區(qū)風(fēng)險(xiǎn)等級(jí)較高，氣候風(fēng)險(xiǎn)綜合指數(shù)為50. 01～65.60。稻縱卷葉螟在甘肅、吉林、黑龍江、內(nèi)蒙古、新疆、山西、寧夏和青海的風(fēng)險(xiǎn)等級(jí)較低，綜合指數(shù)為2.75～2.98。

圖2　2050年稻縱卷葉螟在中國(guó)的適生區(qū)域區(qū)劃圖

2050年，A1b情景下，上海、湖南、江蘇、臺(tái)灣、浙江、湖北、廣西、重慶、貴州和安徽等地風(fēng)險(xiǎn)等級(jí)較高，氣候風(fēng)險(xiǎn)綜合指數(shù)為56.87～71.80。天津、北京、遼寧、河北、西藏、甘肅和青海等地風(fēng)險(xiǎn)等級(jí)較低，綜合指數(shù)低于7.0。

B2a情景下，上海、香港、江蘇、浙江、安徽、福建、湖南、臺(tái)灣、貴州、廣西、重慶和廣東等地風(fēng)險(xiǎn)等級(jí)較高，氣候風(fēng)險(xiǎn)綜合指數(shù)為46.56～74.02。西藏、甘肅、河北、山西、吉林、黑龍江、內(nèi)蒙古、新疆、寧夏和青海等地風(fēng)險(xiǎn)等級(jí)較低，綜合指數(shù)為4.18～6.04。

3結(jié)論與討論

MaxEnt模型是目前國(guó)內(nèi)外在物種適生區(qū)分析領(lǐng)域常用的一種預(yù)測(cè)軟件，它基于機(jī)器學(xué)習(xí)和數(shù)學(xué)統(tǒng)計(jì)，相對(duì)于其他一些常用預(yù)測(cè)模型有其獨(dú)特的優(yōu)勢(shì)[45]。本研究所建立模型的AUC值均大于0.982，模型預(yù)測(cè)準(zhǔn)確性標(biāo)準(zhǔn)為“極好”，結(jié)果適用于稻縱卷葉螟適生區(qū)預(yù)測(cè)的研究。

本文使用MaxEnt和ArcGIS軟件相結(jié)合的方法，分析了當(dāng)前氣候條件下稻縱卷葉螟在我國(guó)各省區(qū)的風(fēng)險(xiǎn)等級(jí)區(qū)劃。結(jié)果表明，稻縱卷葉螟在中國(guó)極高風(fēng)險(xiǎn)區(qū)為上海、浙江大部、江蘇南部、安徽大部、湖北中東部、湖南東部、四川西南部、廣西中部等地，總面積約78.41萬(wàn)km2。高風(fēng)險(xiǎn)區(qū)為貴州、廣東、福建、江西、重慶、廣西大部，總面積約80.2萬(wàn)km2。本文研究了2050年3種排放情景下稻縱卷葉螟在我國(guó)的潛在地理分布，3種情景雖是在特殊情況下的氣候模擬，但從一定程度上反映了在未來(lái)氣候變化的發(fā)展趨勢(shì)。結(jié)果表明，2050年，3種排放情景下，稻縱卷葉螟在中國(guó)的極高風(fēng)險(xiǎn)區(qū)面積均將減少，B2a情景下高風(fēng)險(xiǎn)區(qū)面積將增加。

生態(tài)位模型對(duì)物種生境選擇的分析、預(yù)測(cè)物種潛在分布和物種在環(huán)境變化下的空間分布具有重要意義，目前已被廣泛應(yīng)用于上述領(lǐng)域的研究中。MaxEnt模型實(shí)質(zhì)上是基于貝葉斯定理的一種推斷方法，該模型的工作思路為：根據(jù)物種的分布數(shù)據(jù)(僅需分布區(qū)數(shù)據(jù))和環(huán)境數(shù)據(jù)(包括氣候、海拔、坡度和植被等)，基于目標(biāo)物種的生態(tài)位需求，以最大熵原理作為統(tǒng)計(jì)推斷工具，研究目標(biāo)物種在目標(biāo)區(qū)域的分布幾率及生物多樣性，因此，MaxEnt模型可應(yīng)用于模擬重點(diǎn)監(jiān)測(cè)、檢疫和防控對(duì)象在大尺度環(huán)境背景下的分布及其變化。MaxEnt模型同其它模型相比是一種較新的研究方法，且具有以下優(yōu)點(diǎn)：首先，在基礎(chǔ)數(shù)據(jù)要求方面，MaxEnt模型運(yùn)行僅需要分布區(qū)數(shù)據(jù)，而不需要非分布區(qū)數(shù)據(jù)，現(xiàn)實(shí)中非分布區(qū)數(shù)據(jù)往往獲取難度更大；其次，國(guó)內(nèi)外學(xué)者針對(duì)多種生態(tài)位模型的預(yù)測(cè)結(jié)果的比較分析表明，MaxEnt模型較其它一些模型(包括GARP，BIOCLIM，CLIMEX，DOMAIN等)的AUC值最大，預(yù)測(cè)效果最佳；第三，即使樣本量較小，MaxEnt模型的規(guī)則化程序仍可阻止模擬過(guò)擬合，因此更適合模擬生態(tài)位較窄、分布數(shù)據(jù)有限的物種。基于以上幾點(diǎn)原因，本研究選擇了MaxEnt作為分析軟件對(duì)稻縱卷葉螟的適生性分布進(jìn)行評(píng)估[46-48]。

作為一種基于最大熵的生態(tài)位模型，MaxEnt模型有其優(yōu)點(diǎn)，但同樣尚有不甚滿意的地方。首先，由于模型進(jìn)行模擬時(shí)只使用了有限的“發(fā)生數(shù)據(jù)”，因此預(yù)測(cè)結(jié)果與病蟲(chóng)害實(shí)際的分布有一定偏差，僅代表與分布區(qū)相似的環(huán)境條件；第二，最大熵模型使用模擬時(shí)段的生態(tài)位參數(shù)的方式預(yù)測(cè)物種在未來(lái)的適生性，沒(méi)能以約束形式考慮物種對(duì)氣候變化的適應(yīng)性，這可能導(dǎo)致預(yù)測(cè)結(jié)果不夠客觀甚至有些偏頗；第三，由于模型使用的環(huán)境變量時(shí)間跨度大，以文本為例，“當(dāng)前時(shí)段”跨度達(dá)50年，加之物種對(duì)環(huán)境的適應(yīng)性和響應(yīng)滯后的累積，則很有可能導(dǎo)致最初可被忽視的因素最后變?yōu)橹饕蛩亍ａ槍?duì)以上幾點(diǎn)，要求在應(yīng)用MaxEnt模型進(jìn)行適生分布預(yù)測(cè)時(shí)，首先應(yīng)充分學(xué)習(xí)并掌握目標(biāo)物種的生物學(xué)和生態(tài)學(xué)特性，據(jù)此選擇影響分布的關(guān)鍵環(huán)境因子和精確又有代表性的地理分布點(diǎn)數(shù)據(jù)，提高預(yù)測(cè)結(jié)果的準(zhǔn)確性，建立相對(duì)準(zhǔn)確的模型。在本文分析稻縱卷葉螟在我國(guó)的風(fēng)險(xiǎn)區(qū)時(shí)，選取的19個(gè)因子為與溫度和濕度相關(guān)的生物氣候因子，沒(méi)有考慮影響稻縱卷葉螟分布的其它環(huán)境因素，例如寄主范圍、栽培類型、植被類型、天敵分布等，這必然會(huì)導(dǎo)致預(yù)測(cè)結(jié)果的偏差。在下一步的工作中，應(yīng)綜合考慮包括氣候因素在內(nèi)的多種復(fù)合環(huán)境因素對(duì)稻縱卷葉螟分布的影響，這樣才有可能使預(yù)測(cè)結(jié)果更加準(zhǔn)確。

影響病蟲(chóng)害生長(zhǎng)及分布的主要生態(tài)因子有氣候、生物、土壤等幾大類，其中氣候條件與病蟲(chóng)害生命活動(dòng)關(guān)系密切，是決定病蟲(chóng)害分布的主要因素之一。氣候變化可能對(duì)生態(tài)系統(tǒng)和經(jīng)濟(jì)生產(chǎn)產(chǎn)生不可逆轉(zhuǎn)的影響，國(guó)內(nèi)外研究證實(shí)，氣候變化特別是氣候變暖增加了各地的農(nóng)業(yè)熱量資源，提高了作物的復(fù)種指數(shù)[49-50]，同時(shí)氣候變暖亦可促使病蟲(chóng)害的發(fā)生發(fā)展、分布范圍及危害程度發(fā)生顯著變化[51]。年均氣溫的增高，增加了病蟲(chóng)害發(fā)育所需的可利用積溫，造成害蟲(chóng)成蟲(chóng)發(fā)生期延長(zhǎng)、發(fā)育代數(shù)增多，最終種群數(shù)量大幅度上升；氣溫上升還可導(dǎo)致病蟲(chóng)害分布區(qū)擴(kuò)大，促使害蟲(chóng)向北向西遷移，適宜害蟲(chóng)棲息的地理范圍擴(kuò)大[52-59]。

表3　稻縱卷葉螟在中國(guó)的氣候風(fēng)險(xiǎn)等級(jí)綜合指數(shù)

本研究結(jié)果表明2050時(shí)段，B2a排放情景下，稻縱卷葉螟在我國(guó)的高風(fēng)險(xiǎn)區(qū)面積將增加，總面積達(dá)85.66萬(wàn)km2，說(shuō)明未來(lái)溫室氣體排放模式若為“能源利用較低”情景，則很有可能增加稻縱卷葉螟在我國(guó)的適生度，增加該水稻害蟲(chóng)在我國(guó)的風(fēng)險(xiǎn)范圍；研究還表明，2050年影響稻縱卷葉螟分布的主導(dǎo)氣候因子與目前情景相比有較大變化，這預(yù)示現(xiàn)在影響較小的氣候因子有可能在未來(lái)上升為主要影響因子。這就警示我們?cè)谘芯坎∠x(chóng)害風(fēng)險(xiǎn)區(qū)區(qū)劃時(shí)不僅應(yīng)重視目前的情況，還應(yīng)關(guān)注氣候變化對(duì)病蟲(chóng)害分布的影響。

參考文獻(xiàn)：

[1] 李照會(huì).農(nóng)業(yè)昆蟲(chóng)與鑒定[M].北京:中國(guó)農(nóng)業(yè)出版社,2002:205.

[2] 李傳明,徐健,楊亞軍,等.人工飼料飼養(yǎng)稻縱卷葉螟的生長(zhǎng)發(fā)育與繁殖[J].中國(guó)水稻科學(xué),2011,25(3):321-325.

[3] 白先達(dá),黃超艷,唐廣田,等.氣象條件對(duì)稻縱卷葉螟遷飛的影響分析[J].中國(guó)農(nóng)學(xué)通報(bào),2010,26(21):262-267.

[4] Batty M. Space, scale, and scaling in entropy maximizing[J]. Geographical Analysis, 2010, 42: 395-421.

[5] Wilson A. Entropy in urban and regional modeling: retrospect and prospect[J]. Geographical Analysis, 2010, 42(4): 364-394.

[6] Elith J, Graham H C, Anderson R P, et al. Novel methods improve prediction of species’ distributions from occurrence data[J]. Ecography, 2006, 29(2): 129-151.

[7] 胡忠俊,張鐿鋰,于海彬.基于MaxEnt模型和GIS的青藏高原紫花針茅分布格局模擬[J].應(yīng)用生態(tài)學(xué)報(bào),2015,26(2):505-511.

[8] Phillips S J, Dudík M , Schapire R E. A maximum entropy approach to species distribution modeling[C]//Proceedings of the Twenty-First International Conference on Machine Learning. 2004.

[9] 孫敬松,周廣勝.利用最大熵法(MaxEnt)模擬中國(guó)冬小麥分布區(qū)的年代際動(dòng)態(tài)變化[J].中國(guó)農(nóng)業(yè)氣象,2012,33(4):481-487.

[10] 何奇瑾.我國(guó)玉米種植分布與氣候關(guān)系研究[D].南京:南京信息工程大學(xué),2012.

[11] 朱耿平,劉晨,李敏,等.基于Maxent和GARP模型的日本雙棘長(zhǎng)蠹在中國(guó)的潛在地理分布分析[J].昆蟲(chóng)學(xué)報(bào),2014,57(5):581-586.

[12] 李雙成,高江波.基于GARP模型的紫莖澤蘭空間分布預(yù)測(cè):以云南縱向嶺谷為例[J].生態(tài)學(xué)雜志,2008,27(9):1531-1536.

[13] 楊瑞,張雅林,馮紀(jì)年.利用ENFA生態(tài)位模型分析玉帶鳳蝶和箭環(huán)蝶異地放飛的適生性[J].昆蟲(chóng)學(xué)報(bào),2008,51(3):290-297.

[14] 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.

[15] 王運(yùn)生,謝丙炎,萬(wàn)方浩,等.ROC曲線分析在評(píng)價(jià)入侵物種分布模型中的應(yīng)用[J].生物多樣性,2007,15(4):365-372.

[16] 趙文娟,陳林,丁克堅(jiān),等.利用MaxEnt預(yù)測(cè)玉米霜霉病在中國(guó)的適生區(qū)[J].植物保護(hù),2009,35(2):32-38.

[17] 房鋒,張朝賢,黃紅娟,等.基于Maxent的麥田惡性雜草節(jié)節(jié)麥的潛在分布區(qū)預(yù)測(cè)[J].草業(yè)學(xué)報(bào),2013,229(2):62-70.

[18] 王穎,章桂明,楊偉東,等.基于MaxEnt的大豆南北方莖潰瘍病菌在中國(guó)適生區(qū)的預(yù)測(cè)[J].植物檢疫,2009,23(4):14-17.

[19] 齊會(huì)會(huì),張?jiān)苹?王健,等.稻縱卷葉螟在探照燈下的撲燈節(jié)律[J].中國(guó)農(nóng)業(yè)科學(xué),2014,47(22):4436-4444.

[20] 齊國(guó)君,王政,藍(lán)日青,等.粵北稻區(qū)稻縱卷葉螟的發(fā)生規(guī)律及蟲(chóng)源性質(zhì)分析[J].中國(guó)水稻科學(xué),2013,27(2):177-183.

[21] 廖懷建,黃建榮,方源松,等.水稻和玉米苗上飼養(yǎng)的稻縱卷葉螟對(duì)溫度的反應(yīng)[J].生態(tài)學(xué)報(bào),2013,33(2):409-415.

[22] 蔣春先,楊秀麗,齊會(huì)會(huì),等.中國(guó)華南地區(qū)稻縱卷葉螟遷飛的一次雷達(dá)觀測(cè)[J].中國(guó)農(nóng)業(yè)科學(xué),2012,45(23):4808-4817.

[23] 王鳳英,張孝羲,翟保平.稻縱卷葉螟的飛行和再遷飛能力[J].昆蟲(chóng)學(xué)報(bào),2010,53(11):1265-1272.

[24] 劉芳,奚本貴,包善微,等.氯蟲(chóng)苯甲酰胺對(duì)稻縱卷葉螟的防效及對(duì)稻田有益節(jié)肢動(dòng)物的安全性評(píng)價(jià)[J].植物保護(hù),2009,35(5):139-144.

[25] 王鳳英,胡高,陳曉,等.近年來(lái)廣西南寧稻縱卷葉螟大發(fā)生原因分析[J].中國(guó)水稻科學(xué),2009,23(5):537-545.

[26] 徐德進(jìn),顧中言,徐廣春,等.霧滴密度及大小對(duì)氯蟲(chóng)苯甲酰胺防治稻縱卷葉螟效果的影響[J].中國(guó)農(nóng)業(yè)科學(xué),2012,45(4):666-674.

[27] 黃學(xué)飛,張孝羲,翟保平.交配對(duì)稻縱卷葉螟飛行能力及再遷飛能力的影響[J].南京農(nóng)業(yè)大學(xué)學(xué)報(bào),2010,33(5):23-28.

[28] 張巍,張志罡,付秀芹,等.轉(zhuǎn)Bt基因水稻對(duì)稻縱卷葉螟幼蟲(chóng)體內(nèi)三種保護(hù)酶活性的影響[J].昆蟲(chóng)學(xué)報(bào),2008,51(10):1022-1027.

[29] 高忠文,陶嶺梅,蘇建亞,等.防治稻縱卷葉螟高毒農(nóng)藥替代藥劑的室內(nèi)篩選[J].中國(guó)水稻科學(xué),2008,22(6):631-636.

[30] 齊國(guó)君,秦冉冉,肖滿開(kāi),等.安徽安慶混作稻區(qū)稻縱卷葉螟第三、四代發(fā)生規(guī)律研究[J].中國(guó)水稻科學(xué),2008,22(5):513-518.

[31] 王亓翔,許路,吳進(jìn)才.水稻品種對(duì)稻縱卷葉螟抗性的物理及生化機(jī)制[J].昆蟲(chóng)學(xué)報(bào),2008,51(12):1265-1270.

[32] 高月波,陳曉,陳鐘榮,等.稻縱卷葉螟(Cnaphalocrocismedinalis)遷飛的多普勒昆蟲(chóng)雷達(dá)觀測(cè)及動(dòng)態(tài)[J].生態(tài)學(xué)報(bào),2008,28(11):5238-5247.

[33] 王翠花,包云軒,王建強(qiáng),等.2003年稻縱卷葉螟重大遷入過(guò)程的大氣動(dòng)力機(jī)制分析[J].昆蟲(chóng)學(xué)報(bào),2006,49(4):604-612.

[34] 盧輝,鐘義海,陳青,等.基于MaxEnt模型的對(duì)粒材小蠹的適生性分析[J].熱帶作物學(xué)報(bào),2013,34(11):2239-2245.

[35] 雷軍成,徐海根.基于MaxEnt的加拿大一枝黃花在中國(guó)的潛在分布區(qū)預(yù)測(cè)[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2010,26(2):137-141.

[36] 齊國(guó)君,陳婷,高燕,等.基于Maxent的大洋臀紋粉蚧和南洋臀紋粉蚧在中國(guó)的適生區(qū)分析[J].環(huán)境昆蟲(chóng)學(xué)報(bào),2015,37(2):219-223.

[37] 趙力,朱耿平,李敏,等.入侵害蟲(chóng)西部喙緣蝽和紅肩美姬緣蝽在中國(guó)的潛在分布[J].天津師范大學(xué)學(xué)報(bào):自然科學(xué)版,2015,35(1):75-78.

[38] 程軍,沈陽(yáng),郭水良.氣候變暖對(duì)花葉蘚屬植物在中國(guó)的潛在分布范圍影響的預(yù)測(cè)[J].生態(tài)科學(xué),2015,34(1):9-16.

[39] 郭彥龍,衛(wèi)海燕,路春燕,等.氣候變化下桃兒七潛在地理分布的預(yù)測(cè)[J].植物生態(tài)學(xué)報(bào),2014,38(3):249-261.

[40] 張熙驁,隋曉云,呂植,等.基于Maxent的兩種入侵性魚(yú)類(麥穗魚(yú)和鯽)的全球適生區(qū)預(yù)測(cè)[J].生物多樣性,2014,22(2):182-188.

[41] 張穎,李君,林蔚,等.基于最大熵生態(tài)位元模型的入侵雜草春飛蓬在中國(guó)潛在分布區(qū)的預(yù)測(cè)[J].應(yīng)用生態(tài)學(xué)報(bào),2011,22(11):2970-2976.

[42] 張曉峰,佟友貴,楊建東,等.基于Maxent生態(tài)位模型預(yù)測(cè)球果角脛象在云南的潛在地理分布[J].環(huán)境昆蟲(chóng)學(xué)報(bào),2013,35(1):28-32.

[43] 韓陽(yáng)陽(yáng),王焱,項(xiàng)楊,等.基于Maxent生態(tài)位模型的松材線蟲(chóng)在中國(guó)的適生區(qū)預(yù)測(cè)分析[J].南京林業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版,2015,39(1):6-10.

[44] 孫兵,婁玉霞,郭水良.野生蔬菜鴨兒芹在中國(guó)栽培的氣候適應(yīng)性研究[J].上海農(nóng)業(yè)學(xué)報(bào),2014,30(1):85-89.

[45] 王運(yùn)生,謝丙炎,萬(wàn)方浩,等.相似穿孔線蟲(chóng)在中國(guó)的適生區(qū)預(yù)測(cè)[J].中國(guó)農(nóng)業(yè)科學(xué),2007,40(11):2502-2506.

[46] Brotons L, Thuiller W, Araujo M B, et al. Presence-absence versus presence-only modeling methods for predicting bird habitat suitability[J]. Ecography, 2004, 27(4): 437-448.

[47] 徐燁,王金龍,劉雅婷,等.鳳仙花壞死斑病毒在中國(guó)的適生性風(fēng)險(xiǎn)評(píng)估[J].植物保護(hù),2013,39(6)：26-31.

[48] 曲偉偉,李志紅,黃貴修,等.利用MAXENT預(yù)測(cè)橡膠樹(shù)棒孢霉落葉病在中國(guó)的適生區(qū)[J].植物保護(hù),2011,37(4):52-57.

[49] 趙宗慈.為IPCC第五次評(píng)估報(bào)告提供的全球氣候模式預(yù)估[J].氣候變化研究進(jìn)展,2009,5(4):241-243.

[50] 張國(guó)慶.氣候變化對(duì)生物災(zāi)害發(fā)生的影響及對(duì)策[J].現(xiàn)代農(nóng)業(yè)科技,2011(1):318-321.

[51] 李劍泉,李智勇,易浩若.森林與全球氣候變化的關(guān)系[J].災(zāi)害學(xué),2011,26(2):35-38.

[52] 王鴻斌,張真,孔祥波,等.入侵害蟲(chóng)紅脂大小蠹的適生區(qū)和適生寄主分析[J].林業(yè)科學(xué),2007,43(10):71-76.

[53] 馬菲,何友元,趙林忠,等.桃蛀螟在新疆的適生區(qū)預(yù)測(cè)[J].植物保護(hù),2012,38(2)：136-138.

[54] 王培新,杜薇薇,李孟樓.氣候變暖與林木害蟲(chóng)的適應(yīng)性變化[J].西北林學(xué)院學(xué)報(bào),2011,26(3):124-128.

[55] 吳志剛,曲偉偉,張澤華,等.基于CLIMEX的苜蓿籽象甲在中國(guó)的適生區(qū)分析[J].植物保護(hù),2012,38(3):63-66.

[56] 何善勇,溫俊寶,駱有慶,等.氣候變暖情景下松材線蟲(chóng)在我國(guó)適生區(qū)范圍[J].應(yīng)用昆蟲(chóng)學(xué)報(bào),2012(1):236-243.

[57] 張立敏,陳斌,李正躍,等.氣候變化對(duì)水稻三化螟物候及多度的影響(英文)[J].西南農(nóng)業(yè)學(xué)報(bào),2012,25(3):885-893.

[58] 劉明春,蔣菊芳,郭小芹,等.氣候變化背景下玉米棉鈴蟲(chóng)消長(zhǎng)動(dòng)態(tài)預(yù)測(cè)及影響因素研究[J].干旱地區(qū)農(nóng)業(yè)研究,2014,32(2):114-118,141.

[59] 馬春森,馬罡,趙飛.氣候變暖對(duì)麥蚜的影響[J].應(yīng)用昆蟲(chóng)學(xué)報(bào),2014,51(6):1435-1443.

(責(zé)任編輯：許晶晶)

Risk Analysis ofCnaphalocrocismedinalisin China under Three Emissions Scenarios in 2050

WANG Yan-li1,2,WANG Ru-lin1, JIANG Gan1, SHEN Zhan-hong1, LIN Shan1

(1. Rural Economic Information Centre of Sichuan Province, Chengdu 610072, China; 2. Sichuan Provincial Key Laboratory of Heavy Rain and Drought-Flood Disasters in Plateau and Basin, Chengdu 610072, China)

Abstract：The potential distribution ofCnaphalocrocismedinalisGuenee in 2050 of China was predicted by using MaxEnt and ArcGIS model. The results showed that, Shanghai, Zhejiang, Anhui, south of Jiangsu, east of Hubei, east of Hunan, southeast of Sichuan, middle of Guangxi were the extremely high risk area forCnaphalocrocismedinalisin China. The high risk areas were in Guizhou, Guangdong, Fujian, Jiangxi, Chongqing, and most of Guagnxi. In 2050, the areas ofCnaphalocrocismedinaliswill tend to decrease in extremely high risk areas, while the high risk areas will tend to expand under B2a emissions scenario.

Key words：CnaphalocrocismedinalisGuenee; MaxEnt; Risk analysis; Prediction

收稿日期：2015-10-27

基金項(xiàng)目：四川省農(nóng)村經(jīng)濟(jì)綜合信息中心業(yè)務(wù)技術(shù)攻關(guān)課題(農(nóng)信課題201403)；中國(guó)氣象局“2012年業(yè)務(wù)專項(xiàng)經(jīng)費(fèi)”資助。

作者簡(jiǎn)介：王閆利(1978─)，女，工程師，從事農(nóng)業(yè)信息化管理、病蟲(chóng)害與氣象關(guān)系研究、為農(nóng)氣象信息管理工作。

中圖分類號(hào)：S431.9

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1001-8581(2016)06-0076-06