武萬(wàn)里 劉垚 馬菁 張學(xué)儉

摘要：為評(píng)估分析寧夏區(qū)域枸杞癭螨危害的風(fēng)險(xiǎn)特征，探討枸杞癭螨致災(zāi)的氣象因子指標(biāo)，基于2013-2018年寧夏全區(qū)枸杞癭螨田間自動(dòng)采集器、病蟲(chóng)害移動(dòng)采集系統(tǒng)監(jiān)測(cè)資料，應(yīng)用Maxent模型及ArcGIS技術(shù)，構(gòu)建了枸杞癭螨適生分布模型，分析了生態(tài)環(huán)境因子對(duì)枸杞癭螨分布的影響。結(jié)果表明：對(duì)枸杞癭螨分布影響較大的環(huán)境變量主要為6月、7月、8月、9月降水量、10月的平均氣溫和最高地溫;構(gòu)建的枸杞癭螨潛在地理分布模型的AUC值達(dá)到0.950，屬“非常好”的標(biāo)準(zhǔn)，可用于研究枸杞癭螨適生性地理分布研究;枸杞癭螨的適生區(qū)范圍主要在中寧核心枸杞產(chǎn)區(qū)和銀川市賀蘭山沿山枸杞產(chǎn)區(qū);枸杞癭螨高、次高、中等氣候風(fēng)險(xiǎn)區(qū)的面積分別為1323.01 km2、2779.96 km2、4346.14 km2，分別占寧夏國(guó)土總面積的1.99%、4.19%、6.55%。研究結(jié)果對(duì)于制定合理防控措施、控制枸杞癭螨災(zāi)害具有重要的指導(dǎo)意義。

關(guān)鍵詞：枸杞;枸杞癭螨;Maxent模型;適生性分析

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

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1008-0457（2022）02-0066-005

國(guó)際DOI編碼：10.15958/j.cnki.sdnyswxb.2022.02.010

枸杞是寧夏的農(nóng)業(yè)優(yōu)勢(shì)特色產(chǎn)業(yè)之一，目前全區(qū)枸杞種植面積已超過(guò)600 km2，占全國(guó)種植面積一半以上，形成了以中寧為核心，清水河流域和銀川以北為重點(diǎn)布局的枸杞產(chǎn)業(yè)經(jīng)濟(jì)帶。但近些年來(lái)，隨著氣候變化以及枸杞產(chǎn)區(qū)面積的逐年擴(kuò)大，枸杞病蟲(chóng)害的發(fā)生也越來(lái)越嚴(yán)重，對(duì)枸杞的產(chǎn)量和品質(zhì)造成嚴(yán)重影響。枸杞癭螨（拉丁學(xué)名Aceria macrodonis ?Keifer），又稱大瘤癭螨[1]，主要分布于寧夏、內(nèi)蒙古、甘肅、新疆、山西、陜西、青海等地的枸杞產(chǎn)區(qū)，是枸杞的主要害蟲(chóng)之一，主要危害枸杞的葉片、花蕾、幼果、嫩莖、花瓣及花柄[2]，優(yōu)先侵染枸杞幼嫩部位[3]。嫩莖受害時(shí)，在頂端葉芽處會(huì)形成長(zhǎng)3～5 mm的丘狀蟲(chóng)癭，受害嚴(yán)重時(shí)葉片會(huì)扭曲變形，嚴(yán)重的頂端嫩葉卷曲膨大成拳頭狀，變成褐色，提前脫落，造成禿頂枝條，停止生長(zhǎng)，危害嚴(yán)重時(shí)，可能造成50%以上減產(chǎn)。在寧夏地區(qū)，枸杞癭螨為常發(fā)害蟲(chóng)，1年可繁殖十余代。每年4月，越冬的成螨開(kāi)始出蟄活動(dòng)，5-6月份和8-9月份出現(xiàn)2次危害高峰，11月成螨開(kāi)始越冬，并通過(guò)枸杞木虱成蟲(chóng)攜帶越冬[4]。由于枸杞癭螨發(fā)生的時(shí)間跨度長(zhǎng)，大多數(shù)時(shí)間生活于蟲(chóng)癭內(nèi)，防治難度大，防治成本高，因此，科學(xué)有效地進(jìn)行枸杞癭螨監(jiān)測(cè)預(yù)報(bào)工作至關(guān)重要。生態(tài)位理論是生態(tài)學(xué)最重要的理論之一[5]，廣泛應(yīng)用于生物生境分析和預(yù)測(cè)工作中。其中，最大熵（Maxent）模型是目前應(yīng)用較廣的用于模擬給定環(huán)境條件下物種出現(xiàn)概率的物種分布模型[6]。多年的研究證明，Maxent模型在物種分布預(yù)測(cè)研究和應(yīng)用上表現(xiàn)出了很好的預(yù)測(cè)能力，它具有模擬精度高、操作簡(jiǎn)便、結(jié)果穩(wěn)定、所需樣本量小等優(yōu)點(diǎn)[7]，廣泛應(yīng)用在物種生境適宜區(qū)的預(yù)測(cè)和評(píng)價(jià)工作中。實(shí)踐中，Li等[8]利用最大熵模型（Maxent）和地理信息系統(tǒng)（GIS）軟件，根據(jù)三種氣候情景模擬了椰棕扁葉甲的未來(lái)分布。孫佩珊等[9]利用Maxent模型預(yù)測(cè)了地中海實(shí)蠅在中國(guó)的適生區(qū)域，評(píng)估了其入侵中國(guó)的風(fēng)險(xiǎn)。王蕾等[10]基于Maxent模型分析了新疆特色林果區(qū)春尺蠖發(fā)生的風(fēng)險(xiǎn)，ROC評(píng)價(jià)顯示，對(duì)新疆特色林果區(qū)春尺蠖潛在風(fēng)險(xiǎn)區(qū)模擬效果較好。劉博等[11]利用Maxent模型預(yù)測(cè)了小巢粉虱在中國(guó)的潛在地理分布。Li等[12]利用Maxent生態(tài)位模擬軟件，模擬了云斑天牛的分布，預(yù)測(cè)了氣候變化情景下的潛在變化范圍。馬菁等[13]利用Maxent模型預(yù)測(cè)了寧夏枸杞蚜蟲(chóng)適生分布。本研究擬應(yīng)用Maxent模型，分析寧夏區(qū)域枸杞癭螨的危害特征，探討枸杞癭螨致災(zāi)的氣象因子等指標(biāo)參數(shù)，構(gòu)建枸杞癭螨適生地理分布預(yù)測(cè)模型，分析寧夏全區(qū)枸杞癭螨適生分布區(qū)域，研究寧夏枸杞癭螨的適生地理分布特征，旨在為及時(shí)掌握寧夏枸杞癭螨地理分布特征及發(fā)生范圍、為制定合理防控措施、控制癭螨災(zāi)害提供參考和理論依據(jù)。

1數(shù)據(jù)與方法

1.1枸杞癭螨分布數(shù)據(jù)采集

本研究共設(shè)立了21個(gè)枸杞癭螨蟲(chóng)害監(jiān)測(cè)點(diǎn)，分別位于中寧核心產(chǎn)區(qū)、清水河流域、銀川及以北三大枸杞種植帶的21個(gè)規(guī)模化枸杞種植基地，其中，中寧核心產(chǎn)區(qū)有：中寧百瑞源、鑫陽(yáng)、紅寶、大地生態(tài)、中杞、杞泰、早康、舟塔、杞源祥、不老子、吳忠杞愛(ài)、金沙灣枸杞種植基地;銀川及以北植產(chǎn)區(qū)有：銀川百瑞源、南梁農(nóng)場(chǎng)、園林場(chǎng)、暖泉、惠農(nóng)西夏堂枸杞種植基地;清水河流域種植產(chǎn)區(qū)有：紅寺堡百瑞源、同心潤(rùn)德、易捷、固原正杞紅枸杞種植基地。采用田間自動(dòng)采集器和病蟲(chóng)害移動(dòng)采集系統(tǒng)，結(jié)合人工調(diào)查的方式采集蟲(chóng)口數(shù)量和密度。具體做法是：（1）利用田間自動(dòng)采集器對(duì)枸杞癭螨進(jìn)行區(qū)域化遠(yuǎn)程監(jiān)測(cè)，將信息素黃板上監(jiān)控到的枸杞癭螨等蟲(chóng)害的蟲(chóng)口數(shù)量和密度圖像信息，通過(guò)3G或4G等無(wú)線通信模式傳輸至后臺(tái)數(shù)據(jù)服務(wù)器，并進(jìn)行圖像識(shí)別處理入庫(kù)，遠(yuǎn)程獲取各枸杞種植基地病蟲(chóng)監(jiān)測(cè)點(diǎn)的枸杞癭螨監(jiān)測(cè)信息;（2）利用枸杞病蟲(chóng)害移動(dòng)采集系統(tǒng)，采取隨機(jī)監(jiān)測(cè)點(diǎn)監(jiān)測(cè)和人工客戶端填報(bào)等方式，將監(jiān)測(cè)信息傳輸至后臺(tái)服務(wù)器。采集時(shí)間為2013—2018年。5年間共計(jì)獲取點(diǎn)位數(shù)據(jù)34 300余條。

1.2氣象數(shù)據(jù)獲取

用于枸杞癭螨風(fēng)險(xiǎn)分布預(yù)測(cè)分析的氣象數(shù)據(jù)，源自寧夏氣象檔案館，包含近30年的氣溫、降水、輻射、風(fēng)速月平均數(shù)據(jù)，以及19個(gè)能反映監(jiān)測(cè)區(qū)域水熱變化的氣候變量，來(lái)源于Worldclim氣候數(shù)據(jù)庫(kù)（https：//worldclim.org/data/worldclim21.html）。

1.3建模方法與評(píng)價(jià)標(biāo)準(zhǔn)

使用Maxent模型對(duì)寧夏枸杞癭螨風(fēng)險(xiǎn)分布區(qū)域進(jìn)行預(yù)測(cè)分析，并利用ArcGIS進(jìn)行空間特征分析，具體作法為：（1）將實(shí)際調(diào)查的枸杞癭螨地理信息數(shù)據(jù)統(tǒng)一整理為Maxent模型指定的CSV格式數(shù)據(jù)，作為樣本輸入數(shù)據(jù);利用ArcGIS提取與枸杞癭螨地理信息數(shù)據(jù)所對(duì)應(yīng)的各變量數(shù)據(jù)值，作為模型輸入的變量數(shù)據(jù);（2）使用刀切法評(píng)判各變量的貢獻(xiàn)率，剔除貢獻(xiàn)率小的變量。將枸杞癭螨分布數(shù)據(jù)和4-10月降水量、氣溫以及19個(gè)生物氣候變量等，再次輸入模型，其中，隨機(jī)選取75%的監(jiān)測(cè)數(shù)據(jù)作為訓(xùn)練樣本用于建模，剩余的25%監(jiān)測(cè)數(shù)據(jù)用于模型檢驗(yàn);（3）利用AUC值評(píng)價(jià)模型的模擬精度。AUC值是ROC曲線下的面積，取值范圍在0～1之間，越接近1，表示模型預(yù)測(cè)精度越高[14]。AUC值的評(píng)價(jià)按五級(jí)劃分，即0.5≤AUC<0.6、0.6≤AUC<0.7、0.7≤AUC<0.8、0.8≤AUC<0.9、0.9≤AUC<1.0共5個(gè)區(qū)間。對(duì)應(yīng)的評(píng)價(jià)等級(jí)為失敗、較差、一般、好、非常好5個(gè)評(píng)價(jià)等級(jí)。通常，當(dāng)AUC值>0.75時(shí)，可判定模型模擬結(jié)果準(zhǔn)確度高;（4）Maxent模型的預(yù)測(cè)輸出結(jié)果為枸杞癭螨的存在概率，利用ArcGIS軟件插值，可得到枸杞癭螨的存在概率分布值。將存在概率劃分為高風(fēng)險(xiǎn)區(qū)（存在概率>0.8）、次高風(fēng)險(xiǎn)區(qū)（0.6≤存在概率<0.8）、中風(fēng)險(xiǎn)區(qū)（0.4≤存在概率<0.6）、次低風(fēng)險(xiǎn)區(qū)（0.2≤存在概率<0.4）和低風(fēng)險(xiǎn)區(qū)（存在概率<0.2）5個(gè)等級(jí)，并繪制枸杞癭螨風(fēng)險(xiǎn)分布預(yù)測(cè)圖。

2結(jié)果與分析

2.1模型檢驗(yàn)

ROC曲線法（即AUC法）是Maxent模型精度評(píng)價(jià)上應(yīng)用最廣的方法，用于模擬結(jié)果的驗(yàn)證。從ROC曲線圖（圖1）可以看出，基于所有變量構(gòu)建的枸杞癭螨潛在地理分布模型的AUC值為0.950，表明所構(gòu)建模型的預(yù)測(cè)準(zhǔn)確性達(dá)到“非常好”的標(biāo)準(zhǔn)，可用于枸杞癭螨潛在地理分布研究。

2.2各變量對(duì)模型預(yù)測(cè)的貢獻(xiàn)率

通過(guò)刀切法對(duì)氣象數(shù)據(jù)變量進(jìn)行相關(guān)性分析，得到各變量對(duì)枸杞癭螨分布影響的貢獻(xiàn)率（圖2），AUC值排名前六的變量分別為：6月降水量、10月氣溫、10月最高地溫、7月降水量、9月降水量、8月降水量。其中，6月降水量和10月氣溫AUC值>0.9，達(dá)到“非常好”的標(biāo)準(zhǔn)，10月最高地溫、7月降水量、9月降水量、8月降水量AUC值介于0.8～0.9之間，屬于“好”的標(biāo)準(zhǔn)。

2.3各氣象變量與癭螨潛在分布的關(guān)系

枸杞癭螨潛在地理分布模型預(yù)測(cè)結(jié)果顯示（圖3），6月降水量、10月平均氣溫、10月最高地溫、7月降水量、9月降水量、8月降水量對(duì)枸杞癭螨發(fā)生發(fā)展影響較大。6月、7月和8月降水量在20～40 mm適宜枸杞癭螨發(fā)生發(fā)展，基本在30 mm左右蟲(chóng)口密度達(dá)到最大，之后蟲(chóng)口密度隨著降水量的增加逐漸下降。9月降水量在15 mm以下和30～45 mm范圍較適宜枸杞癭螨發(fā)生發(fā)展。10月平均氣溫和10月最高地溫對(duì)枸杞癭螨發(fā)生影響較為明顯，適宜的10月平均氣溫范圍在6.0 ℃以上，適宜的10月最高地溫范圍在38 ℃以上。

2.4枸杞癭螨風(fēng)險(xiǎn)分析

基于Maxent模型的寧夏枸杞癭螨風(fēng)險(xiǎn)分布預(yù)測(cè)結(jié)果，利用ArcGIS軟件得到枸杞癭螨的存在概率分布。寧夏枸杞癭螨蟲(chóng)害風(fēng)險(xiǎn)區(qū)主要集中在銀川市賀蘭山沿山枸杞產(chǎn)區(qū)和中寧枸杞產(chǎn)區(qū)（圖4）。其中，高和次高氣候風(fēng)險(xiǎn)區(qū)位于中寧枸杞核心產(chǎn)區(qū);中等氣候風(fēng)險(xiǎn)區(qū)位于銀川市賀蘭山沿山枸杞產(chǎn)區(qū)及清水河流域北部的枸杞產(chǎn)區(qū)。枸杞癭螨的高、次高、中等氣候風(fēng)險(xiǎn)區(qū)的面積分別為1323.01 km2、2779.96 km2、4346.14 km2，分別占寧夏國(guó)土面積的1.99%、4.19%、6.55%。

3結(jié)論與討論

本研究基于枸杞癭螨地理分布信息、環(huán)境及氣象變量，應(yīng)用Maxent模型對(duì)寧夏枸杞癭螨進(jìn)行了適生性分析。結(jié)果表明，寧夏枸杞癭螨蟲(chóng)害風(fēng)險(xiǎn)區(qū)主要集中在銀川市賀蘭山沿山枸杞產(chǎn)區(qū)和中寧枸杞產(chǎn)區(qū)，其中，高和次高氣候風(fēng)險(xiǎn)區(qū)位于中寧枸杞核心產(chǎn)區(qū);中等氣候風(fēng)險(xiǎn)區(qū)位于銀川市賀蘭山沿山枸杞產(chǎn)區(qū)及清水河流域北部的枸杞產(chǎn)區(qū);對(duì)枸杞癭螨蟲(chóng)害分布影響較大的環(huán)境變量主要為6月降水量、10月氣溫、10月最高地溫、7月降水量、9月降水量、8月降水量，枸杞癭螨發(fā)生的適宜降水量范圍為6月降水量在20 mm以上，7月、8月和9月降水量適宜范圍在25 ?mm以上，10月平均氣溫范圍在6.0 ℃以上，最高地溫范圍在38 ℃以上。

枸杞癭螨體型較小，活動(dòng)非常隱蔽，遷徙距離短，在枸杞的生長(zhǎng)季終生在蟲(chóng)癭內(nèi)繁殖，1年中可繁殖多達(dá)10代，世代重疊嚴(yán)重，因此，防治難度較大。另一方面，由于枸杞癭螨可通過(guò)枸杞木虱攜帶在本地越冬，因此，枸杞癭螨發(fā)生的區(qū)域特征比較明顯。

Maxent模型用于寧夏枸杞癭螨風(fēng)險(xiǎn)分布的研究，具有較高的準(zhǔn)確性，但也存在一定的局限性，一是模型所選用的變量都是氣候因子，未涉及自然、植被等其他環(huán)境因子，其結(jié)果偏重于氣候因子對(duì)枸杞癭螨的風(fēng)險(xiǎn)分布影響。二是Maxent模型是基于基礎(chǔ)生態(tài)位理論構(gòu)建的模型，影響因素只考慮了非生物因素而沒(méi)有考慮生物因素對(duì)物種分布產(chǎn)生的影響，如物種代際變化、天敵等，特別是枸杞癭螨與枸杞木虱的越冬攜帶關(guān)系，實(shí)際情況要復(fù)雜的多，所以得到的預(yù)測(cè)結(jié)果比現(xiàn)實(shí)的生態(tài)位要寬。

參考文獻(xiàn)：

[1]張潤(rùn)志.枸杞癭螨Aceria macrodonis Keifer蟲(chóng)癭及若螨[J].應(yīng)用昆蟲(chóng)學(xué)報(bào)，2016，53（4）：705.

[2]劉賽，李建領(lǐng)，郭昆，等.枸杞癭螨對(duì)枸杞侵染的初步研究[J].中藥材，2019，42（6）：1237-1241.

[3]李建領(lǐng)，劉賽，徐常青，等.寧夏枸杞主要害蟲(chóng)發(fā)生規(guī)律與防治策略[J].中國(guó)現(xiàn)代中藥，2017，19（11）：1599-1604.

[4]劉賽，楊孟可，李建領(lǐng)，等.我國(guó)枸杞主產(chǎn)區(qū)癭螨鑒定及其越冬調(diào)查[J].中國(guó)中藥雜志，2019，44（11）：2208-2212.

[5]曹嶸，陳娟，白光潤(rùn).生態(tài)位理論在我國(guó)城市發(fā)展中的應(yīng)用[J].地理與地理信息科學(xué)，2003，19（1）：62-65.

[6]Phillips S J，Dudik M.Modeling of species distributions with Maxent：new extensions and a comprehensive evaluation[J].Ecography，2008，31（2）：16l-175.

[7]Steven J P，Miroslav D，Robert E S.A maximum entropy approach to species distribution modeling//Proceedings of the twenty-first international conference on machine learning[C].Canada：Banff Alberta，2004：655-662.

[8]Li X Y，Xu D P，Jin Y.Predicting the current and future distributions of Brontispa longissima （Coleoptera：Chrysomelidae） under climate change in China[J].Global Ecology and Conservation，2020，25：e01444.

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

[10]王蕾，羅磊，劉平，等.基于Maxent模型分析新疆特色林果區(qū)春尺蠖發(fā)生風(fēng)險(xiǎn)[J].新疆農(nóng)業(yè)科學(xué)，2019，56（9）：1691-1700.

[11]劉博，覃偉權(quán)，閻偉.基于Maxent模型的小巢粉虱在中國(guó)的潛在地理分布[J].環(huán)境昆蟲(chóng)學(xué)報(bào)，2019，41（6）：1276-1286.

[12]Li A N，Wang J W，Wang R L，et al.Maxent modeling to predict current and future distributions of Batocera lineolata（Coleoptera：Cerambycidae） under climate change in China[J].coscience，2020，27（1）：23-31.

[13]馬菁，劉垚，張學(xué)儉.基于Maxent模型的寧夏枸杞蚜蟲(chóng)適生分布預(yù)測(cè)研究[J].地理與地理信息科學(xué)，2019，35（5）：40-45.

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

Adaptability Analysis of Aceria macrodonis Keifer in Ningxia Based on Maxent Model

Wu Wanli1，2，Liu Yao1，2*，Ma Jing3，Zhang Xuejian3

（1.Key Laboratory for Meteorological Disaster Monitoring and Early Warning and Risk Management of Characteristic Agriculture in Arid Regions，CMA，Yinchuan，Ningxia 750002，China;2.Ningxia Key Lab of Meteorological Disaster Prevention and Reduction，Yinchuan，Ningxia 750002，China;3.Agricultural Economy and Information Technology Research Institute of Ningxia Academy of Agriculture and Forestry，Yinchuan，Ningxia 750002，China）

Abstract：In order to evaluate the damage risk of Aceria macrodonis Keifer in Ningxia，and to explore impact of meteorological factors on distribution of A.macrodonis，a suitable distribution model of A.macrodonis was constructed by Maxent model and ArcGIS，based on the monitoring data of A.macrodonis in Ningxia from 2013 to 2018.The effects of ecological factors on the distribution of A.macrodonis were analyzed.The results showed that precipitation in June，July，August and September，mean temperature and maximum ground temperature in October were the main environment variables that affected the distribution of A.macrodonis.The AUC value of potential distribution model of A.macrodonis reached 0.950，which means the model can be used to study ecological suitability of A.macrodonis.The suitable areas of A.macrodonis were mainly in the core producing areas of lycium barbarum in Zhongning and areas along Helanshan，Yinchuan.The area of the high，second high and medium climate risk areas of A.macrodonis were 1323.01 km2，2779.96 km2 and 4346.14 km2，respectively，accounting for 1.99%，4.19% and 6.55% of the total land area of Ningxia.The results have important guiding significance for formulating reasonable prevention and control measures and controlling the A.macrodonis.

Keywords：Chinese wolfberry;Aceria macrodonis Keifer;Maxent model;adaptability analysis

收稿日期：2021-06-25;

修回日期：2021-08-04

基金項(xiàng)目：第五批“寧夏青年科技人才托舉工程”項(xiàng)目;寧夏農(nóng)林科學(xué)院全產(chǎn)業(yè)鏈創(chuàng)新示范項(xiàng)目（QCYL-2018-1102）

通訊作者：劉垚（1987—），女，博士，高級(jí)工程師，主要從事農(nóng)業(yè)病蟲(chóng)害氣象相關(guān)研究，Email：liuyao314@163.com.

2264501705363