摘要：【目的】預(yù)測(cè)巴西含羞草和無(wú)刺含羞草2種入侵植物在海南省的潛在地理分布，確定2種入侵植物在海南省的入侵風(fēng)險(xiǎn)等級(jí)，為海南省入侵植物管理和防治工作提供參考?！痉椒ā坷肕axEnt模型和GIS地理信息系統(tǒng)評(píng)估制約巴西含羞草和無(wú)刺含羞草的主要環(huán)境變量，模擬并劃分2種入侵植物在海南省的潛在入侵地理分布；開(kāi)展2種入侵植物定性定量入侵風(fēng)險(xiǎn)分析，確定其危險(xiǎn)級(jí)別?！窘Y(jié)果】模擬結(jié)果顯示，巴西含羞草和無(wú)刺含羞草受試者工作特征曲線面積（AUC）分別為0.956和0.996，結(jié)果遠(yuǎn)高于隨機(jī)預(yù)測(cè)值0.500，表明模型預(yù)測(cè)結(jié)果具有較高的可信度。模型預(yù)測(cè)結(jié)果顯示，2種入侵植物在我國(guó)的適生區(qū)主要分布在南方，尤其以廣東、海南、福建等省份分布最為密集；其中，在海南省的潛在適宜區(qū)域主要集中在東北部，并有向西南擴(kuò)散的趨勢(shì)，與實(shí)際調(diào)查結(jié)果分布區(qū)域基本吻合。緯度分布預(yù)測(cè)結(jié)果發(fā)現(xiàn)，巴西含羞草主要分布在北緯50o～南緯50o，無(wú)刺含羞草則主要分布在北緯10o～30o。根據(jù)風(fēng)險(xiǎn)指標(biāo)體系準(zhǔn)則層賦值，計(jì)算得出巴西含羞草風(fēng)險(xiǎn)值（R）=69.15，高于61.00，為特別危害物種；無(wú)刺含羞草R=58.00，介于29.50～61.00，為高度危害物種?！窘Y(jié)論】巴西含羞草為特別危害物種，應(yīng)嚴(yán)格禁止引入；無(wú)刺含羞草屬于高度危害，必須嚴(yán)格限制引入。由于2種入侵植物已經(jīng)傳入我國(guó)，應(yīng)采取必要的應(yīng)急措施，開(kāi)展檢疫、監(jiān)測(cè)和滅除等工作，防止其進(jìn)一步擴(kuò)散和造成更大的危害。

關(guān)鍵詞：巴西含羞草；無(wú)刺含羞草；MaxEnt模型預(yù)測(cè)；入侵風(fēng)險(xiǎn)分析；海南省

中圖分類(lèi)號(hào)：S45 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：2095-1191（2024）11-3324-09

The distribution patterns and risk analysis of 2 invasive plants in Hainan

LI Xiao-xia， ZENG Li-wang， ZENG An-yi， DONG Ding-chao*

（Institute of Scientific and Technology Information， China Academy of Tropical Agricultural Sciences/Key Laboratory of

Applied Research on Tropical Crop Information Technology of Hainan， Haikou， Hainan 571101， China）

Abstract：【Objective】To predict the potential geographic distribution of Mimosa diplotricha C. Wright and Mimosa diplotricha var. inermis（ Adelbert） Veldkamp in Hainan， and to determine the invasion risk level of the 2 invasive plants in Hainan， so as to provide reference for the management of invasive plants and control of invasive plants in Hainan. 【Method】The MaxEnt model and geographic information system（ GIS） were used to assess the main environmental varia-bles governing Mimosa diplotricha C. Wright and Mimosa diplotricha var. inermis （Adelbert） Veldkamp， and to simu‐late and delineate the potential geographic distribution of the 2 invasive plants in Hainan； qualitative and quantitative inva‐sion risk analyses of the 2 invasive plants were carried out to determine the level of their risk. The simulation results showed that area under a receiver operating characteristic curve（ AUC） of Mimosa diplotricha C. Wright and Mimosa dip‐lotricha var. inermis（ Adelbert） Veldkamp were 0.956 and 0.996 respectively， which was much higher than the random prediction value（ 0.500）， indicating that the prediction results of the model had high confidence.【 Result】The model pre‐diction results showed that in China，the suitable habitats for 2 invasive plants were mainly distributed in the southern re‐gions， particularly in provinces such as Guangdong， Hainan and Fujian. In Hainan， the potential suitable areas were mainly concentrated in the northeastern part and showed a trend of spreading towards the southwest，which was basically consistent with the distribution areas of actual survey results. According to the latitude distribution prediction results， Mi‐mosa diplotricha C. Wright mainly distributed between 50o northern latitude and 50o southern latitude，while Mimosa dip‐lotricha var. inermis （Adelbert） Veldkamp mainly distributed between 10o northern latitude and 30o northern latitude. Based on the risk index system criteria layer assignment， the calculated risk value R for Mimosa diplotricha C. Wright was 69.15， which was higher than 61.00， classifying it as a particularly harmful species； the risk value R for Mimosa dip‐lotricha var. inermis（ Adelbert） Veldkamp was 58.00， falling between 29.50 and 61.00， classifying it as a highly harmful species. 【Conclusion】Mimosa diplotricha C. Wright is a particularly hazardous species and its introduction should be strictly prohibited； Mimosa diplotricha var. inermis（ Adelbert） Veldkamp is highly hazardous and its introduction must be strictly restricted. Since the 2 invasive plants have been introduced into China， necessary emergency measures should be taken to carry out quarantine， monitoring and eradication， to prevent them from further spreading and causing greater harm.

Key words： Mimosa diplotricha C. Wright； Mimosa diplotricha var. inermis（ Adelbert） Veldkamp； MaxEnt model prediction； invasion risk anlysis； Hainan province

Foundation items： Hainan Natural Science Foundation（624MS105）； Hainan Philosophy and Social Science Plan‐ning Project［HNSK（YB）23-75］； Young Scholars Project of Hainan Natural Science Foundation（323QN299）； Open Fund Project of Key Laboratory of Applied Research on Tropical Crop Information Technology of Hainan in 2023（ZDSYS-KFJJ-202315）

0 引言

【研究意義】隨著全球氣候變化和人類(lèi)活動(dòng)不斷擴(kuò)展，外來(lái)物種入侵的數(shù)量呈上升趨勢(shì)，給全球范圍內(nèi)造成了巨大的經(jīng)濟(jì)損失（Paini et al.，2016；陶永祥等，2017；Bacher et al.，2018；Masood，2018）。由于外來(lái)入侵植物治理難度大，對(duì)其進(jìn)行潛在分布區(qū)預(yù)測(cè)以及建立有效的預(yù)警機(jī)制，是控制其種群擴(kuò)張的關(guān)鍵手段和重要基礎(chǔ)（李慧琪等，2015；Zhang et al.，2021）。海南是我國(guó)最大的熱帶島嶼，擁有豐富的生物多樣性和獨(dú)特的生態(tài)系統(tǒng)，是我國(guó)的重要生態(tài)安全屏障和生物資源庫(kù)，同時(shí)也是我國(guó)入侵植物的重災(zāi)區(qū)。當(dāng)前，巴西含羞草（Mimosa diplotricha C. Wright）和無(wú)刺含羞草［Mimosa diplotricha var. inermis （Adelbert） Veldkamp］已在海南省?？凇①僦?、澄邁和三亞等多個(gè)市（縣）分布，并呈迅猛擴(kuò)散態(tài)勢(shì)，對(duì)海南省的農(nóng)業(yè)生產(chǎn)和生態(tài)環(huán)境造成了嚴(yán)重?fù)p害。因此，了解和評(píng)估巴西含羞草和無(wú)刺含羞草2種入侵植物在海南省的分布格局和風(fēng)險(xiǎn)狀況，對(duì)于制定有效的防治策略和保護(hù)海南省的生物多樣性和生態(tài)安全具有重要意義?！厩叭搜芯窟M(jìn)展】目前，生態(tài)位模型［最大熵模型（MaxEnt）、BIOCLIM模型和DOMAIN模型］已被廣泛用于預(yù)測(cè)物種的潛在適生區(qū)（Song et al.，2020；Anand et al.，2021）。與其他生態(tài)模型相比，MaxEnt模型能準(zhǔn)確地模擬大部分外來(lái)入侵物種的潛在分布區(qū)域，是目前具有較高預(yù)測(cè)精度的生態(tài)位模型（王聰?shù)龋?018；Xian et al.，2023）。研究顯示，通過(guò)MaxEnt模型預(yù)測(cè)外來(lái)入侵物種的潛在分布區(qū)域，可以更準(zhǔn)確地評(píng)估其可能帶來(lái)的風(fēng)險(xiǎn)和危害，進(jìn)而制定針對(duì)性的防控措施。Yang等（2018）研究入侵植物銀膠菊（Galinsoga quadrira‐diata）擴(kuò)張對(duì)氣候變化的響應(yīng)，結(jié)果發(fā)現(xiàn)，在氣候變化條件下銀膠菊在中國(guó)的適宜區(qū)可能減少并向北移動(dòng)，東北地區(qū)入侵風(fēng)險(xiǎn)增加，需緊急采取預(yù)防措施。塞依丁·海米提等（2019）運(yùn)用MaxEnt模型預(yù)測(cè)刺蒼耳（Xanthium spinosum）在新疆的潛在分布格局，結(jié)果發(fā)現(xiàn)，在2種不同的氣候變化情景下，至2070年，刺蒼耳分布區(qū)中心可能向新疆伊犁州的奎屯方向移動(dòng)。Arogoundad等（2020）結(jié)合MaxEnt模型、遙感數(shù)據(jù)和環(huán)境變量，對(duì)夸祖魯—納塔爾省的銀膠菊進(jìn)行了模擬，并基于此預(yù)測(cè)了5種潛在的外來(lái)入侵植物在內(nèi)布拉斯加州北普拉特河地區(qū)蔓延的可能性。Zhang等（2021）研究桉樹(shù)枝癭姬小蜂（Leptocybe invasa）在中國(guó)的潛在分布核心，預(yù)測(cè)結(jié)果有助于識(shí)別和遏制其擴(kuò)散趨勢(shì)，為制定有效的控制措施提供理論支持。潘銘心等（2022）利用MaxEnt模型預(yù)測(cè)加拿大一枝黃花（Solidago canadensis L.）在我國(guó)的潛在適生區(qū)，結(jié)果發(fā)現(xiàn)，其高度適生區(qū)集中在多個(gè)省份和沿海地區(qū)，為防控該物種入侵提供了重要依據(jù)。陳思明（2023）、方鏵等（2023）利用MaxEnt模型預(yù)測(cè)互花米草（Spartina alterniflora）在我國(guó)潛在適生區(qū)的研究發(fā)現(xiàn)，福建漳江口紅樹(shù)林等保護(hù)區(qū)面臨較高的入侵風(fēng)險(xiǎn)。Xian等（2023）對(duì)紫莖澤蘭（Ageratina adenophora）入侵中國(guó)后生態(tài)位擴(kuò)張的研究發(fā)現(xiàn)，紫莖澤蘭在墨西哥（原產(chǎn)地）和中國(guó)（入侵地）之間的生態(tài)位重疊度較低，在中國(guó)潛在地理分布主要集中在西南地區(qū)的熱帶和亞熱帶地區(qū)，其地理分布格局由降水和溫度變量共同決定。李曉霞等（2024）運(yùn)用MaxEnt模型對(duì)小花十萬(wàn)錯(cuò)（Asystasia gangetica subsp. micrantha）的研究發(fā)現(xiàn)，小花十萬(wàn)錯(cuò)潛在分布受氣候變暖影響，預(yù)測(cè)其將北遷且集中于華南沿海及高海拔地區(qū)，因此，應(yīng)針對(duì)高溫多雨、溫差小區(qū)域制定小花十萬(wàn)錯(cuò)重點(diǎn)防治措施。林敬梧和桂東偉（2024）運(yùn)用MaxEnt模型預(yù)測(cè)火炬樹(shù)（Rhus typhina）在北美和中國(guó)的空間分布差異，結(jié)果表明，火炬樹(shù)在北美和中國(guó)適宜分布區(qū)集中于北緯30°～40°濕潤(rùn)地帶，入侵中國(guó)后其氣候空間和生態(tài)位發(fā)生較大變化，需加強(qiáng)防疫監(jiān)測(cè)和監(jiān)管以保護(hù)生物多樣性?！颈狙芯壳腥朦c(diǎn)】野外調(diào)查發(fā)現(xiàn)，巴西含羞草和無(wú)刺含羞草具有廣泛的生態(tài)適應(yīng)性和繁殖能力，且能在多種環(huán)境中形成單一群落，對(duì)本地植物的生長(zhǎng)和生態(tài)平衡造成威脅。然而，對(duì)于巴西含羞草和無(wú)刺含羞草在海南省的具體分布格局、風(fēng)險(xiǎn)狀況及其與環(huán)境因子的關(guān)系等方面的研究仍需進(jìn)一步深入。【擬解決的關(guān)鍵問(wèn)題】利用MaxEnt模型和GIS地理信息系統(tǒng)評(píng)估制約巴西含羞草和無(wú)刺含羞草2種入侵植物的主要環(huán)境變量，模擬并劃分其在海南省的潛在入侵地理分布，并開(kāi)展2種入侵植物定性定量入侵風(fēng)險(xiǎn)分析，確定其危險(xiǎn)級(jí)別，為海南省入侵植物管理和防治工作提供參考。

1 材料與方法

1. 1 數(shù)據(jù)來(lái)源及處理

1. 1. 1 巴西含羞草和無(wú)刺含羞草分布數(shù)據(jù)的收集與處理 通過(guò)檢索Inaturalist、全球生物多樣性信息共享平臺(tái)以及野外實(shí)地調(diào)查，分別搜集到巴西含羞草標(biāo)本記錄1502條、無(wú)刺含羞草標(biāo)本記錄45條。將收集的數(shù)據(jù)排布至Excel表格中，去除無(wú)效樣本（相同采集地點(diǎn)和時(shí)間的重復(fù)樣本、地理位置和時(shí)間不清楚的樣本）。將坐標(biāo)通過(guò)ArcGIS 10.8映射在世界地圖上，并通過(guò)SDMtoolbox工具的減小空間自變量功能處理，最終得到巴西含羞草位置記錄1225條、無(wú)刺含羞草位置記錄35條，并轉(zhuǎn)化成csv格式進(jìn)行保存。

1. 1. 2 環(huán)境變量的篩選和數(shù)據(jù)處理 環(huán)境氣候數(shù)據(jù)來(lái)源于WorldClimate全球氣候數(shù)據(jù)庫(kù)（https：//www.worldclim.org/），包含歷史條件下的生物氣候變量和未來(lái)（2060—2080年）的氣候數(shù)據(jù)變量，氣候數(shù)據(jù)精度為30″。通過(guò)聯(lián)合國(guó)糧農(nóng)組織土壤門(mén)戶網(wǎng)站（https：//www.fao.org/）下載海拔（GloElev）和土壤質(zhì)量（Sq1～Sq7）等數(shù)據(jù)。19個(gè)生物氣候變量環(huán)境因子之間有一定的相關(guān)性，7個(gè)土壤質(zhì)量因子與植被分布存在一定相關(guān)性。利用ArcGIS 10.8中的SDMtoolbox插件Conversion Tools將環(huán)境圖層轉(zhuǎn)換成MaxEnt可用的ASCII格式。

1. 2 風(fēng)險(xiǎn)評(píng)估方法

本研究基于層次分析法（Analytic hierarchy pro‐cess，AHP）原理，參照張瑞海等（2019）、張國(guó)良等（2020）提出的外來(lái)植物風(fēng)險(xiǎn)指標(biāo)體系和判斷標(biāo)準(zhǔn)，包括4個(gè)主要指標(biāo)層、17個(gè)具體指標(biāo)和32個(gè)問(wèn)題進(jìn)行量化賦值（張?jiān)赖龋?022），并計(jì)算出入侵植物的風(fēng)險(xiǎn)值（R），即R=0.40R1+0.15R2+0.20R3+0.25R4，式中，R1為入侵性、R2為適應(yīng)性、R3為擴(kuò)散性、R4為危害性。根據(jù)R值確定危害級(jí)別。

2 結(jié)果與分析

2. 1 巴西含羞草和無(wú)刺含羞草適生區(qū)分析

2. 1. 1 當(dāng)前氣候條件下巴西含羞草和無(wú)刺含羞草分布 巴西含羞草主要分布于美國(guó)、墨西哥、哥倫比亞、委內(nèi)瑞拉、巴西、印度以及東南亞國(guó)家等（圖1），我國(guó)已確認(rèn)的分布地有臺(tái)灣、香港、廣州、東莞、三亞和?？诘?；無(wú)刺含羞草則主要分布于我國(guó)澳門(mén)、深圳、惠州、廣州等珠三角城市及海南?。▓D2）。通過(guò)對(duì)2種入侵植物生長(zhǎng)區(qū)域的對(duì)比分析，發(fā)現(xiàn)巴西含羞草主要分布在北緯50o～南緯50o，無(wú)刺含羞草則主要分布在北緯10o～30o；2種入侵植物在我國(guó)的適生區(qū)主要分布在南方，尤其以廣東、海南、福建等省份分布最為密集。

2. 1. 2 巴西含羞草和無(wú)刺含羞草在海南省的適生區(qū)預(yù)測(cè) 為進(jìn)一步提高預(yù)測(cè)結(jié)果的精確性，本研究使用MaxEnt模型進(jìn)行了5次環(huán)境因子的重要性分析，結(jié)果顯示，最冷月最低溫度、最暖季度降水量、最濕季度降水量、最冷季度平均溫度和最暖月最高溫度等環(huán)境因子對(duì)2種入侵植物的生長(zhǎng)影響較大。2種入侵植物受試者工作曲線面積（AUC）分別為0.956和0.996，結(jié)果遠(yuǎn)高于隨機(jī)預(yù)測(cè)值0.500，表明模型預(yù)測(cè)結(jié)果具有較高的可信度。

將MaxEnt得出的預(yù)測(cè)結(jié)果導(dǎo)入到ArcGIS 10.8中，并對(duì)柵格數(shù)據(jù)進(jìn)行可視化處理。具體如下：0～1的連續(xù)柵格數(shù)據(jù)中，0代表最低生存概率，1代表最高生存概率。生存概率越接近1，表示生存適宜程度越高。采用自然間斷點(diǎn)分級(jí)法將預(yù)測(cè)結(jié)果劃分為多個(gè)等級(jí)，設(shè)定入侵植物生存概率在［0，0.3）區(qū)間內(nèi)的地區(qū)為低適生區(qū)，生存概率在［0.3，0.7）區(qū)間內(nèi)的地區(qū)為中適生區(qū)，生存概率在［0.7，1］區(qū)間內(nèi)的地區(qū)為高適生區(qū)。同時(shí)，將生存概率≥0.7作為高適生區(qū)的判別閾值，對(duì)2種入侵植物的適生性等級(jí)圖進(jìn)行二值化處理，并將已入侵區(qū)域和潛在入侵區(qū)域進(jìn)行比較分析，結(jié)合實(shí)地踏查數(shù)據(jù)，劃定2種入侵植物的管控區(qū)、高風(fēng)險(xiǎn)入侵區(qū)和監(jiān)測(cè)區(qū)域。

從圖3可看出，當(dāng)前氣候背景下，除樂(lè)東、三亞、瓊中黎族苗族自治縣、保亭黎族苗族自治縣、五指山市、陵水黎族自治縣、白沙黎族自治縣和昌江黎族自治縣等地2種入侵植物的適生等級(jí)存在差異外，其余地區(qū)的適生區(qū)基本相同。數(shù)據(jù)庫(kù)以及實(shí)際調(diào)研所得到巴西含羞草和無(wú)刺含羞草的分布數(shù)據(jù)基本集中在高熱、高濕地區(qū)，證明模型預(yù)測(cè)與2種植物的實(shí)際分布吻合，預(yù)測(cè)結(jié)果有較高的可信度。從圖3還可看出，巴西含羞草和無(wú)刺含羞草入侵風(fēng)險(xiǎn)較大的地區(qū)均集中在海南省的東北部地區(qū)，但無(wú)刺含羞草的入侵風(fēng)險(xiǎn)區(qū)域面積更大。

2. 1. 3 影響巴西含羞草和無(wú)刺含羞草分布格局的氣候因子 為更全面地了解巴西含羞草和無(wú)刺含羞草的生存環(huán)境需求，本研究分析了2種入侵植物原產(chǎn)地和入侵地的分布數(shù)據(jù)以及全球環(huán)境數(shù)據(jù)。同時(shí)，根據(jù)MaxEnt模型運(yùn)算結(jié)果，確定其高適生區(qū)4個(gè)主要環(huán)境因子的變化范圍，并計(jì)算其在模型中的貢獻(xiàn)率。結(jié)果（表1）顯示，2種入侵植物生長(zhǎng)環(huán)境因子有較高的相似性，即包含最冷月最低溫度、最暖季度降水量、最濕季度降水量、最冷季度平均溫度等4個(gè)高貢獻(xiàn)率環(huán)境因子，且4個(gè)高貢獻(xiàn)率環(huán)境因子對(duì)MaxEnt模型的總貢獻(xiàn)率在70.0%以上，說(shuō)明上述4個(gè)環(huán)境因子是影響巴西含羞草和無(wú)刺含羞草分布的主要依賴環(huán)境因子，且影響因子存在高度相似；同時(shí)也表明兩者的生長(zhǎng)環(huán)境相似，因此兩者有極大可能生長(zhǎng)在同一片區(qū)域。

2. 1. 4 巴西含羞草和無(wú)刺含羞草在海南省的預(yù)警分析 在MaxEnt分析結(jié)果基礎(chǔ)上對(duì)巴西含羞草和無(wú)刺含羞草的適宜生境等級(jí)圖進(jìn)行疊加處理，識(shí)別出高度適宜的生境區(qū)域。這一過(guò)程產(chǎn)生了代表高適宜區(qū)域的柵格數(shù)據(jù)集，隨后將柵格數(shù)據(jù)集與行政區(qū)域圖層結(jié)合，生成巴西含羞草和無(wú)刺含羞草適生性風(fēng)險(xiǎn)預(yù)測(cè)圖（圖4和圖5）。從圖4和圖5可看出，巴西含羞草和無(wú)刺含羞草已擴(kuò)散到海南省東部高適生區(qū)的邊緣，繼續(xù)向適生度較低的西南方向蔓延。建議在瓊海市、萬(wàn)寧市、定安縣和澄邁縣等市（縣）建立阻截帶，同時(shí)在瓊中黎族苗族自治縣和陵水黎族自治縣等建立入侵檢測(cè)機(jī)制，并且在其他未被入侵地區(qū)開(kāi)展檢測(cè)調(diào)查。

2. 2 巴西含羞草和無(wú)刺含羞草入侵風(fēng)險(xiǎn)分析結(jié)果

2. 2. 1 定性分析

2. 2. 1. 1 全球范圍內(nèi)的入侵與影響 （1）巴西含羞草原產(chǎn)巴西，為東南亞和太平洋島嶼上常見(jiàn)的嚴(yán)重雜草，廣布于非洲中部、太平洋地區(qū)多個(gè)島嶼、加勒比海地區(qū)以及毛里求斯、越南、柬埔寨、馬來(lái)西亞、印度尼西亞、菲律賓和澳大利亞等地（車(chē)晉滇，2010），是太平洋沿岸的一種嚴(yán)重危害雜草。在薩摩亞西部，巴西含羞草已成為農(nóng)業(yè)、牧場(chǎng)、荒地和路邊等地的主要雜草；在烏波盧島，有高達(dá)85%的村莊受到巴西含羞草侵襲；在澳大利亞，巴西含羞草已成為熱帶牧場(chǎng)和北昆士蘭甘蔗園的主要雜草；巴西含羞草最初作為花卉被引入我國(guó)，隨后逸生至野外，其最早的記載可追溯到1956年的《廣州植物志》（侯寬昭，1956）和《海南植物志 第二卷》（陳煥鏞，1965）。目前，巴西含羞草在我國(guó)的廣東、廣西、海南、臺(tái)灣、香港和云南等地均有分布。（2）無(wú)刺含羞草原產(chǎn)印度尼西亞，現(xiàn)廣布于全世界熱帶地區(qū)?！稄V西植物名錄 第二冊(cè)》（廣西植物研究所，1971）《廣東植物志 第五卷》（吳德鄰，2003）《廣西植物志 第二卷》（李樹(shù)剛，2005）均收載本變種，其標(biāo)本采集于海南萬(wàn)寧。目前，無(wú)刺含羞草在我國(guó)澳門(mén)、香港、福建、廣東、廣西、海南和香港均有分布，因其對(duì)農(nóng)業(yè)生產(chǎn)和生物多樣性的嚴(yán)重影響而被列為我國(guó)二級(jí)嚴(yán)重入侵物種。

2. 2. 1. 2 傳播性與繁殖能力 巴西含羞草和無(wú)刺含羞草因具有種子數(shù)量多、壽命持久、傳播途徑多樣和繁殖能力強(qiáng)等特性而被廣泛傳播（韋春強(qiáng)等，2013）。巴西含羞草莢果上有刺毛，易黏附在動(dòng)物的皮毛上隨動(dòng)物移動(dòng)而傳播；其種子?；祀s于各種器物或土壤中，并隨交通工具擴(kuò)散至各個(gè)角落；也可借助水流進(jìn)行傳播；此外，人為無(wú)意識(shí)的攜帶行為也是其傳播的一條途徑。無(wú)刺含羞草種子可隨苗木、草皮交易傳播；也可通過(guò)人為引種栽培，逸生。巴西含羞草和無(wú)刺含羞草在熱帶、亞熱帶地區(qū)生態(tài)適應(yīng)性強(qiáng)，生長(zhǎng)迅速，能在多種土壤和環(huán)境條件下生存。根據(jù)MaxEnt模型分析結(jié)果可知，巴西含羞草和無(wú)刺含羞草的適生區(qū)包含了海南省的大部分農(nóng)業(yè)產(chǎn)區(qū)、山林、生態(tài)保護(hù)區(qū)，故2種入侵植物均可在海南省形成完整的傳播鏈并進(jìn)行大規(guī)模繁殖。

2. 2. 1. 3 經(jīng)濟(jì)與生態(tài)影響 巴西含羞草和無(wú)刺含羞草與農(nóng)作物競(jìng)爭(zhēng)光、水、養(yǎng)分和空間，降低農(nóng)業(yè)生產(chǎn)的效率和質(zhì)量，且在大量生長(zhǎng)后影響農(nóng)業(yè)機(jī)械的正常運(yùn)行，增加農(nóng)業(yè)生產(chǎn)的成本和難度；同時(shí)，也會(huì)破壞土壤結(jié)構(gòu)和肥力，導(dǎo)致土壤侵蝕和退化；此外，巴西含羞草和無(wú)刺含羞草還會(huì)改變土壤的微生物群落，影響土壤的生物多樣性和功能。隨著2種植物入侵的加劇，會(huì)排擠和取代本地的植物物種，降低植被的多樣性和穩(wěn)定性，甚至影響其他生物，如昆蟲(chóng)、鳥(niǎo)類(lèi)和野生動(dòng)物等的生存和繁衍，破壞生態(tài)系統(tǒng)的平衡。

2. 2. 1. 4 檢疫與防治難易度 巴西含羞草和無(wú)刺含羞草由于種子數(shù)量多、壽命長(zhǎng)、傳播途徑多、生態(tài)適應(yīng)性強(qiáng)、抗藥性強(qiáng)，容易形成單一的群落且難以根除；其種子還可以通過(guò)水流、風(fēng)力、動(dòng)物、車(chē)輛、衣物等方式進(jìn)行遠(yuǎn)距離傳播，因而難以有效控制（韋春強(qiáng)等，2013）。此外，巴西含羞草和無(wú)刺含羞草的植株具有很強(qiáng)的適應(yīng)性和再生能力，可以在各種環(huán)境條件下生長(zhǎng)，難以殺滅。目前，巴西含羞草和無(wú)刺含羞草的防治方法主要包括物理、化學(xué)和生物防治等（馬金雙，2014）。

2. 2. 2 定量分析 根據(jù)外來(lái)植物風(fēng)險(xiǎn)指標(biāo)體系和判斷標(biāo)準(zhǔn)，并參照相關(guān)外來(lái)物種風(fēng)險(xiǎn)評(píng)估方法（丁暉等2006；張國(guó)良等，2020），提出巴西含羞草和無(wú)刺含羞草風(fēng)險(xiǎn)指標(biāo)體系和判斷標(biāo)準(zhǔn)，進(jìn)行定量分析。該定量分析指標(biāo)體系準(zhǔn)則層包括4個(gè)指標(biāo)層（R1、R2、R3和R4）的17個(gè)具體指標(biāo)和32個(gè)問(wèn)題，根據(jù)評(píng)價(jià)指標(biāo)的賦分值表（張?jiān)赖龋?022），設(shè)定2種入侵植物R為100，其中R1賦值40、R2賦值15、R3賦值20、R4賦值25。按R=R1×0.40+R2×0.15+R3×0.20+R4×0.25計(jì)算風(fēng)險(xiǎn)值，確定危險(xiǎn)級(jí)別（表2）。通過(guò)計(jì)算得出巴西含羞草R=69.15，超過(guò)了閾值61.00，為特別危害物種；無(wú)刺含羞草R=58.00，介于閾值29.50～61.00，為高度危害物種。針對(duì)特別危害的巴西含羞草應(yīng)堅(jiān)決禁止其引入，以防對(duì)本地生態(tài)系統(tǒng)造成不可逆的損害；而對(duì)于高度危害的無(wú)刺含羞草，雖不完全禁止其引入，但應(yīng)嚴(yán)格限制其數(shù)量，并采取足夠的監(jiān)管和防范措施，以控制其逃逸和擴(kuò)散，確保本地生態(tài)安全。

3 討論

本研究利用MaxEnt模型結(jié)合GIS技術(shù)分析了巴西含羞草和無(wú)刺含羞草在海南省的分布格局，發(fā)現(xiàn)2種入侵植物的潛在適宜區(qū)域主要集中在海南東北部地區(qū)，并有向西南方向擴(kuò)散的趨勢(shì)，與其現(xiàn)有的分布區(qū)域基本一致。說(shuō)明MaxEnt模型能較好地預(yù)測(cè)入侵植物的分布范圍，與其他研究者的結(jié)論相符（岳茂峰等，2016；李璇等，2018）。MaxEnt模型預(yù)測(cè)結(jié)果還顯示，2種入侵植物的分布受到多個(gè)環(huán)境因子影響，其中最重要的是最冷月最低溫度、最暖季度降水量、最濕季度降水量和最冷季度平均溫度，這些環(huán)境因子與入侵植物的生態(tài)適應(yīng)性有關(guān)，也與海南的

氣候和地形特征相吻合。

外來(lái)植物風(fēng)險(xiǎn)評(píng)估包括分析外來(lái)入侵植物特性、環(huán)境和人為因素，評(píng)價(jià)其入侵可能性和影響，從社會(huì)、經(jīng)濟(jì)和政治等方面確定可接受的危害程度，制定和實(shí)施降低風(fēng)險(xiǎn)措施（張國(guó)良等，2020；張?jiān)赖龋?022）。本研究對(duì)巴西含羞草和無(wú)刺含羞草進(jìn)行了分布、危害及防治等分析，建立合適的風(fēng)險(xiǎn)評(píng)估體系，通過(guò)定量分析的辦法確定其危險(xiǎn)級(jí)別。研究結(jié)果表明，巴西含羞草和無(wú)刺含羞草在海南省的入侵風(fēng)險(xiǎn)均較高，其中巴西含羞草的風(fēng)險(xiǎn)略高于無(wú)刺含羞草，與2種入侵植物的生物學(xué)特性和危害現(xiàn)狀等相一致。巴西含羞草具有較強(qiáng)的競(jìng)爭(zhēng)力、繁殖力和適應(yīng)力等，當(dāng)前已在海南省形成了大面積的單一群落，對(duì)本地生物多樣性和生態(tài)系統(tǒng)功能帶來(lái)嚴(yán)重的威脅。無(wú)刺含羞草入侵性雖然沒(méi)有巴西含羞草強(qiáng)，但也具有較高的生長(zhǎng)速度、擴(kuò)散能力和耐干旱性等特性，且在海南省的分布范圍不斷擴(kuò)大，對(duì)農(nóng)業(yè)生產(chǎn)和自然保護(hù)區(qū)造成了不良影響。因此，建議在關(guān)鍵區(qū)域建立生物阻截帶，并設(shè)立長(zhǎng)期監(jiān)測(cè)站點(diǎn)；同時(shí)，加強(qiáng)公眾教育和跨學(xué)科合作，政府出臺(tái)政策法規(guī)支持對(duì)巴西含羞草和無(wú)刺含羞草的管理，控制其擴(kuò)散，保護(hù)本地生物多樣性。

本研究使用的入侵植物分布數(shù)據(jù)是基于文獻(xiàn)和網(wǎng)絡(luò)記錄，可能存在一定的誤差和不完整性，未能反映2種入侵植物的實(shí)際分布狀況。因此，有必要通過(guò)實(shí)地調(diào)查和監(jiān)測(cè)來(lái)獲取更準(zhǔn)確和更新的入侵植物分布數(shù)據(jù)，以提高模型的預(yù)測(cè)精度和可靠性。此外，本研究使用的環(huán)境因子是基于全球尺度的數(shù)據(jù)集，可能存在一定的粗糙度和不適用性，未能充分考慮海南省的地域特征和微觀環(huán)境變化。因此，有必要使用更高分辨率和更具代表性的環(huán)境因子數(shù)據(jù)，以提高模型的預(yù)測(cè)靈敏度和適應(yīng)性。

本研究使用的定性或定量指標(biāo)風(fēng)險(xiǎn)評(píng)估方法可能存在一定的主觀性和模糊性，未能充分反映入侵植物的風(fēng)險(xiǎn)程度和影響范圍。因此，有必要使用更客觀和綜合的風(fēng)險(xiǎn)評(píng)估方法，以提高風(fēng)險(xiǎn)評(píng)估的科學(xué)性和有效性。

4 結(jié)論

MaxEnt模型分析結(jié)果，巴西含羞草和無(wú)刺含羞草在海南省的潛在適宜區(qū)域主要集中在海南省東北部地區(qū)，與其現(xiàn)有的分布區(qū)域基本一致。海南省氣候總體上較適合巴西含羞草和無(wú)刺含羞草生長(zhǎng)，且隨著農(nóng)業(yè)生產(chǎn)、動(dòng)植物活動(dòng)等因素，其擴(kuò)散范圍仍在不斷擴(kuò)大。根據(jù)2種入侵植物的入侵風(fēng)險(xiǎn)分析結(jié)果，巴西含羞草和無(wú)刺含羞草分別為特別危害物種和高度危害物種，對(duì)海南省原生生態(tài)環(huán)境和農(nóng)業(yè)生產(chǎn)均可能造成不良影響，因此，迫切需要采取措施，開(kāi)展檢疫、監(jiān)測(cè)和滅除等，防止其進(jìn)一步擴(kuò)散和造成更大的危害。

參考文獻(xiàn)（（References））：

車(chē)晉滇. 2010. 中國(guó)外來(lái)雜草原色圖鑒［M］. 北京：化學(xué)工業(yè)出版社. ［Che J D. 2010. A color atlas of alien weeds in China［M］. Beijing：Chemical Industry Press.］

陳煥鏞. 1965. 海南植物志：第2卷［M］. 北京：科學(xué)出版社：204-215. ［Chen H Y. 1965. Flora of Hainan，Vol. 2［M］. Bejing：Science Press：204-215.］

陳思明. 2023. 互花米草（Sparina alernilora）潛在分布格局的空間尺度效應(yīng)［J］. 生態(tài)學(xué)報(bào)，43（14）：6058-6068.［ Chen S M. 2023. Spalial seale elect of potential distribution pat‐tern of Sparina alermifora［J］. Acta Eeologica Sinica，43（14）：6058-6068.］ doi：10.5846/stxb202207162032.

丁暉，石碧清，徐海根. 2006. 外來(lái)物種風(fēng)險(xiǎn)評(píng)估指標(biāo)體系和評(píng)估方法［J］. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，22（2）：92-96.［ Ding H，Shi B Q，Xu H G. 2006. Index system and methodology for risk assessment of alien species［J］. Joumal of Ecology and Rural Environment，22（2）：92-96.］ doi：10.3969/j.issn.1673-4831.2006.02.020.

方鏵，陳星彤，劉明月，張永彬，苗正紅，滿衛(wèi)東，張清文，寇財(cái)垚，李想. 2023. 基于MaxEnt 模型的互花米草潛在生境適宜性分析［J］. 草地學(xué)報(bào)，31（11）：3514-3524.［ Fang H，Chen X T，Liu M Y，Zhang Y B，Miao Z H，Man W D，Zhang Q W，Kou C Y，Li X. 2023. Analysis on the poten‐tial habitat suitability of Spartina alterniflora based on MaxEnt model［J］. Acta Agrestia Sinica，31（11）：3514-3524.］ doi：10.11733/i.issn.1007-0435.2023.11.030.

廣西植物研究所. 1971. 廣西植物名錄 第二冊(cè)［M］. 桂林：廣西植物研究所：280.［ Guangxi Institute of Botany. 1971. A checklist of plants in Guangxi，Vol. 2［M］. Guilin：Guangxi Institute of Botany：280.］

侯寬昭. 1956. 廣州植物志［M］. 北京：科學(xué)出版社：311.［ Hou K Z. 1956. Flora of Guangzhou［M］. Beijing：Science Press：311.］

李慧琪，趙力，祝培文，張嫣，朱耿平，劉強(qiáng). 2015. 入侵植物長(zhǎng)芒莧在中國(guó)的潛在分布［J］. 天津師范大學(xué)學(xué)報(bào)（自然科學(xué)版），35（4）：57-61.［ Li H Q，Zhao L，Zhu P W，Zhang Y，Zhu G P，Liu Q. 2015. Potential distribution of invasive palmer amaranth in China（ Amaranthus palmeri）［J］. Jour‐nal of Tianjin Normal University（ Natural Science Edition），35（4）：57-61.］ doi：10.3969/j.issn.1671-1114.2015.04.013.

李樹(shù)剛. 2005. 含羞草科［M］//中國(guó)科學(xué)院廣西植物研究所. 廣西植物志 第二卷. 南寧：廣西科學(xué)技術(shù)出版社：416-441. ［Li S G. 2005. Mimosaceae［M］//Flora of Guangxi，Vol. 2. Guangxi Institute of Botany，Chinese Academy of Sciences. Nanning：Guangxi Science and Technology

Press：416-441.］

李璇，李垚，方炎明. 2018. 基于優(yōu)化的MaxEnt模型預(yù)測(cè)白櫟在中國(guó)的潛在分布區(qū)［J］. 林業(yè)科學(xué)，54（8）：153-164.［ Li X，Li Y，F(xiàn)ang Y M. 2018. Prediction of potential suitable distribution areas of Quercus fabri in China based on an optimized MaxEnt model［J］. Scientia Silvae Sinicae，54（8）：153-164.］ doi：10.11707/j.1001-7488.20180817.

李曉霞，胡寬義，曾安逸，董定超. 2024. 基于MaxEnt生態(tài)位模型的小花十萬(wàn)錯(cuò)在中國(guó)的潛在分布［J］. 西南農(nóng)業(yè)學(xué)報(bào)，37（4）：860-868.［ Li X X，Hu K Y，Zeng A Y，Dong D C. 2024. Potential distribution of Asystasa gangetiea ssp. micrantha in China based on the MaxEnt ecological miche model［J］. Southwest China Journal of Agricultural Seien-ces，37（4）：860-868.］ doi：10.16213/j.enki.scjas.2024.4.020.

林敬梧，桂東偉. 2024. 火炬樹(shù)在北美和中國(guó)空間分布差異對(duì)比及防范入侵措施［J］. 生態(tài)學(xué)報(bào)，44（4）：1692-1699.［ Lin J W，Gui D W. 2024. Comparing spatial distribution of torch tree in North America and China and its preventing measures［J］. Acta Ecologica Sinica，44（4）：1692-1699.］ doi：10.20103/j.stxb.202204100939.

馬金雙. 2014. 中國(guó)外來(lái)入侵植物調(diào)研報(bào)告（下卷）［M］. 北京：高等教育出版社：592-593.［ Ma J S. 2014. The survey reports on Chinese alien invasive plants（Vol. 2）［M］. Bei‐jing：Higher Education Press：592-593.］

潘銘心，朱思睿，張震. 2022. 外來(lái)入侵植物加拿大一枝黃花在中國(guó)的適生區(qū)預(yù)測(cè)［J］. 西安文理學(xué)院學(xué)報(bào)（自然科學(xué)版），25（1）：90-96. ［Pan M X，Zhu S R，Zhang Z. 2022. Prediction of suitable habitat of invasive plant of Solidago canadensis L. in China［J］. Journal of Xi’an University （Natural Science Edition），25（1）：90-96.］ doi：10.3969/j.issn.1008-5564.2022.01.016.

塞依丁·海米提，努爾巴依·阿布都沙力克，許仲林，阿爾曼·解思斯，邵華，維尼拉·伊利哈爾. 2019. 氣候變化情景下外來(lái)入侵植物刺蒼耳在新疆的潛在分布格局模擬［J］. 生態(tài)學(xué)報(bào)，39（5）：1551-1559.［ Sayit H，Nurbay A，Xu Z L，Arman J，Shao H，Vinira Y. 2019. Simulation of potential distribution patterns of the invasive plant species Xanthium spinosum L. （Bathurst burr） in Xinjiang under climate change［J］. Acta Ecologica Sinica，39（5）：1551-1559.］ doi：10.5846/stxb201802040299.

陶永祥，趙建偉，王蘭新，郭賢明. 2017. 西雙版納自然保護(hù)區(qū)外來(lái)入侵植物現(xiàn)狀調(diào)查［J］. 山東林業(yè)科技，（1）：58-61. ［Tao Y X，Zhao J W，Wang L X，Guo X M. 2017. Investi‐gation on the present situation of alien invasive plants in Xishuangbanna Nature Reserve［J］. Shandong Forestry Scinence and Technology，（1）：58-61.］ doi：10.3969/j.issn. 1002-2724.2017.01.014.

王聰，蔡普默，儀傳冬，陳家驊. 2018. 2007—2017年外來(lái)入侵害蟲(chóng)風(fēng)險(xiǎn)評(píng)估的文獻(xiàn)分析及常用風(fēng)險(xiǎn)評(píng)估模型介紹［J］. 中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，23（8）：225-238.［ Wang C，Cai P M，Yi C D，Chen J H. 2018. Bibliometric based analysis of the risk assessment of alien invasive pest from 2007 to 2017 and introduction of several risk assessment models［J］. Journal of China Agricultural University，23（8）：225-238.］ doi：10.11841/i.issn.1007-4333.2018.08.24.

韋春強(qiáng)，趙志國(guó)，丁莉，潘玉梅. 2013. 廣西新記錄入侵植物［J］. 廣西植物，33（2）：275-278. ［Wei C Q，Zhao Z G，Ding L，Pan Y M. 2013. Newly recorded invasive plants in Guangxi［J］. Guihaia，33（2）：275-278.］ doi：10.3969/i.issn.1000-3142.2013.02.026.

吳德鄰. 2003. 含羞草科［M］//中國(guó)科學(xué)院華南植物研究所. 廣東植物志 第五卷. 廣州：廣東科技出版社：141-163. ［Wu D L. 2003. Mimosaceae［M］//South China Institute of Botany，Chinese Academy of Sciences. Flora of Guang‐dong，Vol. 5. Guangzhou：Guangdong Science and Tech‐nology Press：141-163.］

岳茂峰，馮莉，崔燁，張純，田興山. 2016. 基于MaxEnt模型的入侵植物白花鬼針草的分布預(yù)測(cè)及適生性分析［J］. 生物安全學(xué)報(bào)，25（3）：222-228. ［Yue M F，F(xiàn)eng L，Cui Y，Zhang C，Tian X S. 2016. Prediction of the potential distri‐bution and suitability analysis of the invasive weed，Bidens alba（ L.） DC［J］. Jouranl of Biosafety，25（3）：222-228.］ doi：10.3969/j.issn.2095-1787.2016.03.013.

張國(guó)良，付衛(wèi)東，孫玉芳，王忠輝，張瑞海. 2020. 外來(lái)入侵物種監(jiān)測(cè)與控制［M］. 北京：中國(guó)農(nóng)業(yè)出版社.［ Zhang G L，F(xiàn)u W D，Sun Y F，Wang Z H，Zhang R H. 2020. Monito-ring and control of alien invasive species［M］. Beijing：China Agriculture Press.］

張瑞海，張國(guó)良，宋振，王忠輝，付衛(wèi)東. 2019. 印加孔雀草入侵生態(tài)風(fēng)險(xiǎn)評(píng)估及其管理措施［J］. 生物安全學(xué)報(bào)，28（1）：71-75.［ Zhang R H，Zhang G L，Song Z，Wang Z H，F(xiàn)u W D. 2019. The invasive assessment and management measure of Tagetes minuta L［. J］. Journal of Biosafety，28（1）：71-75.］ doi：103969/i.issn.2095-1787.2019.01.013.

張?jiān)溃瑥垏?guó)良，袁至立，宋振，王忠輝，付衛(wèi)東. 2022. 速生槐葉蘋(píng)入侵風(fēng)險(xiǎn)評(píng)估及管理措施分析［J］. 生物安全學(xué)報(bào)，31（4） ：364-374.［ Zhang Y，Zhang G L，Yuan Z L，Song Z，Wang Z H，F(xiàn)u W D. 2022. Invasive risk assessment and management measure analysis of Salvinia molesta D. S. Mitchel［l J］. Journal of Biosafety，31（4）：364-374.］ doi：10.3969/j.issn.2095-1787.2022.04.010.

Anand V，Oinam B，Singh I H. 2021. Predicting the current and future potential spatial distribution of endangered Rucer‐vus eldii eldii（ Sangai） using MaxEnt mode［l J］. Environ‐mental Monitoring and Assessment，193：147-165. doi：10.1007/s10661-021-08950-1.

Arogoundade A M，Odindi J，Mutanga O. 2020. Modelling Parthenium hysterophorus invasion in KwaZulu-Natal Province using remotely sensed data and environmental variables［J］. Geocarto International，35（13）：1450-1465. doi：10.1080/10106049.2019.1581268.

Bacher S，Blackburn T M，Essl F，Cenovesi P，Heikkila J，Jeschke J M，Jones C，Keller R，Kenis M，Kueller C，Marti‐nou A F. Nentwig W，Perl J，Py?ek P，Rabitsch W，Richard‐son D M，Roy H E，Saul W C，Scalera R，Vila M，Wilson J R U，Kumschick S. 2018. Socio-economic impact classifi‐cation of alien taxa （SEICAT）［J］. Methods in Ecology and Evolution，9（1）：159-168. doi：10.1007/10265-015-0738-3.

Masood E. 2018. The battle for the soul of biodiversity［J］. Nature，560（7719）：423-425. doi：10.1038/d41586-018-05984-3.

Paini D R，Sheppard A W，Cook D C，De Barro P J，Worner S P，Thomas M B. 2016. Global threat to agriculture from invasive species［J］. Proceedings of the National Academy of Sciences，113（27）：7575-7579. doi：10.1073/pnas.1602 205113.

Song X Y，Milne R I，F(xiàn)an X Y，Xie S Y，Zhang L，Zheng H L，F(xiàn)an L Q，Chung J M，Chung M G，Ma T Xu X T，Wang J，M K. 2020. Blow to the Northeast？Intraspecific differen‐tiation of Populus davidiana suggests a north-eastward skew of a phylogeographic break in East Asia［J］. Journal of Biogeography，48（1）：187-201. doi：10.1111/jbi.13992.

Xian X G，Zhao H X，Wang R，Zhang H B ，Chen B X，Liu W X，Wan F H. 2023. Evidence of the niche expansion of crofton weed followinginvasion in China［J］. Ecology and Evolution，13（1）：e9708. doi：10.1002/ece3.9708.

Yang Y Bo，Liu G，Xiong S，Zhang W G，Cai X W. Ren Z L . Yao N N，Zhu Z H，Nie H. 2018. Where will invasive plants colonize in response to climate change：Predicting the invasion of Galinsoga quadriradiata in China［J］. Inter‐national Journal of Environmental Research，12：929-938. doi：10.1007/s41742-018-0146-3.

Zhang H，Song J Y，Zhao H X，Li M，Han W H. 2021. Predic-ting the distribution of the invasive species Leptocybe invasa：Combining MaxEnt and geodetector models［J］. Insects，12（2）：92-109. doi：10.3390/insects12020092.

（責(zé)任編輯 麻小燕）