2023, Vol. 44
随着经济全球化的加深，我国进口货物和进口国家数量不断增长，同时，许多有害杂草伴随着国际流通进入我国[1]。研究表明，进口额每增长1个百分点，植物疫情在海关截获量约增长2个百分点[2]。近年来有多起南方三棘果Emex australis在海关截获的报告。例如，2012—2015年广东黄埔口岸在主要进境粮食中截获的有害生物中，南方三棘果截获高达90批次排第7位。此外在广东的肇庆和汕头口岸以及江苏口岸也都有截获报道[3-4]。
南方三棘果是蓼科亦模属的冬型一年生杂草，喜生长在湿润和半湿润热带、亚热带和温带气候区，被我国列入进口植物检疫对象。目前在我国内地尚未发现，仅在台湾有分布记录。南方三棘果单株植株可以产生1000粒以上的果实，种子在土壤中休眠多年，仍有繁殖能力，且具有繁殖速度快、传播迅速、危害性大和防治难度高的特点，一旦大面积入侵我国，难以根治，将会给农业生产造成巨大破坏[5-7]。
随着信息化技术的发展，目前用于预测物种潜在分布区的模型有CLIMEX、GARP、DIVA-GIS、BIOCLIM和MaxEnt等。相比于其他模型，MaxEnt生态位模型在分布数据较少的情况下，预测精确度和准确度高[8-9]。其通过将研究区域的环境因子转换为栅格数据，形成多个栅格单元，把物种已知分布区域的栅格单位作为样点，然后以分布点区域的环境变量数据对所有栅格单元进行数理统计，得到当前和未来一段时间内研究区域的每一个栅格单元分布的概率以及环境变量数值与分布概率之间的关系[10-12]。
由于各物种的研究和关注程度不同，物种在世界范围内详实的分布记录存在差异，因此物种的分布数据获得量存在差异[12]。南方三棘果属于数据获得量较少的物种，因此选择MaxEnt作为预测南方三棘果在我国潜在适生区的模型。全球气温呈现升高趋势，会影响植物的潜在分布区域变化[13-14]。气候系统模式能够很好地反映气候变化的趋势，对未来一段时间的气候变化进行准确预测[15-17]。因此本文采用国家气候中心4种气候情景下的未来气候数据对南方三棘果在我国的适生区进行分析预测。
本研究将MaxEnt模型与ArcGIS软件相结合，得到南方三棘果潜在分布区的可视化地图，以期能够为有关部门在潜在分布区的相关地区开展入侵风险评估、检验检疫和监测预警等提供理论参考，对防控南方三棘果传入、定殖和扩散，保障农业生产和生物多样性健康以及生态安全和经济贸易安全等具有重要意义[18]。
1 材料与方法 1.1 技术路线技术路线包括查找和处理数据、模型优化与构建以及制作可视化地图3部分，按照“物种分布点查找、处理→查找影响分布的环境因子以及SPSS分析→MaxEnt模型参数优化→构建MaxEnt模型→ArcGIS制作可视化图”的研究方案进行，数据分析贯穿于整个研究过程。
1.2 相关软件和程序MaxEnt软件3.4.1版、ArcGIS软件10.4.1版、SPSS Statistics软件20.0版、R语言3.6.3版、ENMeval程序包以及Microsoft Excel 2016。
1.3 环境数据及处理从世界气候数据网站(
|
表 1 环境变量因子 Table 1 Environmental variable factors |
参考文献[21-22]方法，用ArcGIS将所有环境变量坐标系及基准面统一定义，并转出为ASCⅡ格式，以待MaxEnt模型预测分布；使用R语言将TIFF格式的环境变量因子统一转换为ASC格式，以待使用R语言ENMeval 程序包对MaxEnt模型进行优化。为了避免变量间的高度共线性导致的过拟合和减少运算量，对环境数据进行相关性分析。将物种分布点与环境变量输入MaxEnt模型运算10次，获得各个环境变量的贡献率，将环境数据和物种分布点数据在ArcGIS处理为“*．txt”格式文件，将文件利用Excel转化为“*．csv”格式文件，最后将分析数据导入SPSS进行双变量相关性分析，作双尾检验。对相关系数(R)>0.8的环境变量，选择贡献率较高的作为主要变量进行下一步分析，舍去贡献率较低项Bio5、Bio6、Bio7、Bio10、Bio12、Bio16和Bio17。舍去贡献率为0的环境变量Bio8。由于南方三棘果为冬型一年生草本植物，在最湿月和最暖季度是种子形态[23-24]，因此舍去Bio13和Bio18。最终得到12个环境变量因子(表1)，分别为Bio1、Bio2、Bio3、Bio4、Bio9、Bio11、Bio14、Bio15、Bio19、Elev、Slo和Asp。
1.4 南方三棘果分布数据及处理在全球生物多样性信息交换平台(
从标准地图服务系统(
使用R语言ENMeval程序包将南方三棘果分布点分为2部分，75%随机样本数据用于训练集，25%用于测试集，对MaxEnt软件中调控倍频(Regularization multiplier，RM)和特征组合(Feature combination，FC)参数进行优化，利用数据包设定RM从0.1到4，每次增加0.1，共40个调控倍频，同时设定29个特征组合，即H、L 、LH、LP、LPH、LPT、LPTH、LQ、LQH、LQP、LQPH、LQPT、LQPTH、LQT、LQTH、LT、P、PH、PT、PTH、Q、QH、QP、QPH、QPT、QT、QTH、T、TH，其中L为线性(Linear)，Q为二次型(Quadratic)，H为片段化(Hinge)，P为乘积型(Product)，T为阈值型(Threshold)，ENMeval 语言包将上述 1160 种参数组合进行检验，最终采用最小信息准则(The minimum information criterion AICc value)模型评估受试者工作特征曲线(Receiver operating curve，ROC)下的面积(Area under curve，AUC)差异检验模型的拟合度与复杂度，当AICc值最低(AICc=0)时模型参数组合最佳[28]，用于MaxEnt模型建模。
1.7 构建 MaxEnt 预测模型参考前人所述方法[29-33]，将处理后的当前环境变量和物种分布点数据输入MaxEnt模型，模型参数设置随机选择25%的样本作为测试集，75%的样本数据为训练集，同时勾选制作响应曲线功能以获得分布概率与气候因子的关系，勾选刀切法检验环境变量的重要程度，RM与FC取优化后的参数，其他参数为默认设置。以AUC评价模型的准确性。AUC越大，表明模型的可靠性越高，0.9以上表示模型结果优秀可靠，为降低模型的随机性和模拟结果的不确定性，模型重复运行10次，以运行的均值作为模型运算结果。
1.8 绘制物种分布图将MaxEnt模型生成的ASC格式物种分布数据导入ArcGIS软件，设置坐标系为地理坐标系GCS_WGS_1984，通过掩膜提取获得南方三棘果在我国的分布图，对预测结果利用自然间断点分级法进行重分类，得到南方三棘果的非适生区、低适生区、中适生区和高适生区。
2 结果与分析 2.1 MaxEnt生态位模型优化结果基于12个环境变量因子和南方三棘果的647个分布点数据，MaxEnt软件初始默认RM=1，FC为LQPH，AICc=142.939251，经ENMeval语言包优化后RM=1，FC为LQPT，AICc=0，表明其为最优参数组合。在优化参数下进行模型预测，采用ROC分析法进行模型准确度检验。经10次重复训练得到的AUC > 0.9( 图1)，模拟结果优秀可靠。
|
图 1 MaxEnt模型ROC预测结果 Fig. 1 The prediction results of ROC by MaxEnt model |
使用12个环境因子建立模型，模型参数为ENMeval语言包优化结果，RM=1，FC为LQPT，采用刀切法评估各环境因子对南方三棘果地理分布的重要性。贡献率超过15%的包括最冷季度平均气温(Bio11，27.7%)、等温性(Bio3，24.9%)、温度季节变化标准差(Bio4，18.2%)，这3个因子的贡献率合计达70.8%。贡献率低于15%的包括年平均温度(Bio1，14.0%)、最干月降水量(Bio14，6.5%)、降水量变异系数(Bio15，6.5%)、最冷季度降水量(Bio19，2.7%)、高程海拔(Elev，1.4%)、坡度(Slo，0.6%)、最干季度平均温度(Bio9，0.3%)、昼夜温差月均值(Bio2，0.2%)、坡向(Asp，0.2%)。
刀切法检验表明，仅使用单一环境变量时，正则化训练增益和测试增益都表明最重要的环境因子为最冷季度平均温度(Bio11)(图2)。
|
图 2 环境变量重要性的刀切法检验结果 Fig. 2 Results of Jackknife for the importance of the environmental variables Asp: 坡向；Slo: 坡度；Elev: 高程海拔；Bio9: 最干季度平均温度；Bio4: 温度季节变化标准差；Bio3: 等温性 ( Bio2/ Bio7) ×100；Bio2: 昼夜温差月均值；Bio19: 最冷季度降水量；Bio15: 降水量变异系数；Bio14: 最干月降水量；Bio11: 最冷季度平均温度；Bio1: 年平均温度 Asp: Aspect; Slo: Slope; Elev: Elevation; Bio9: Mean temperature of the driest quarter; Bio4: Standard deviation of seasonal variation of temperature; Bio3: Isothermality; Bio2: Mean diurnal range; Bio19: Precipitation of the coldest quarter; Bio15: Precipitation seasonality; Bio14: Precipitation of the driest month; Bio11: Mean temperature of the coldest quarter; Bio1: Annual mean temperature |
一般认为存在概率≥0.5对应的环境因子范围是物种适宜分布的阈值。根据响应曲线确定影响南方三棘果分布环境因子的阈值：最冷季度平均温度9.35~12.76 ℃、等温性45.2%~53.72%以及温度季节变化标准差337.06~616.40 ℃(图3)。
|
图 3 环境因子响应曲线 Fig. 3 Environmental factor response curve |
将南方三棘果分布点和环境变量数据导入MaxEnt模型进行模拟，得到南方三棘果在当前(2022年)气候下的潜在分布区(图4)。总适生面积约69.13万km2，高适生区、中适生区、低适生区分别为11.64万、13.37万、44.02万km2。南方三棘果的适生范围主要集中在我国云南、广西、广东、福建以及西藏局部，其中，高适生区和中适生区主要位于云南、海南以及西藏局部，低适生区主要位于广西、广东、台湾、福建和四川。
|
图 4 南方三棘果在中国的适生区(2022年) Fig. 4 Current suitable regions of Emex australis in China (2022) |
如表4所示，SSP1-2.6气候情景下，2040年南方三棘果总适生区面积约63.54万km2，较2022年分布面积减少约5.48万km2，其中，低适生区面积减少5.93万km2，中适生区面积增加 6.57万km2，高适生区面积减少6.12万km2；SSP2-4.5气候情景下，2040年南方三棘果总适生区面积约65.10万km2，较2022年分布面积减少约3.92万km2，其中，低适生区面积减少2.31万km2，中适生区面积增加2.31万km2，高适生区面积减少3.92万km2；SSP3-7.0气候情景下，2040年南方三棘果总适生区面积约63.34万km2，较2022年分布面积减少约5.68万km2，其中，低适生区面积减少9.74万km2，中适生区面积增加6.51万km2，高适生区面积减少2.45万km2；SSP5-8.5情景下，2040年南方三棘果总适生区面积约62.90万km2，较2022年分布面积减少约6.12万km2，其中，低适生区面积减少5.42万km2，中适生区面积增加3.39万km2，高适生区面积减少4.09万km2。
|
|
表 4 不同时期气候情景下南方三棘果在中国潜在适生区的面积 Table 4 Areas of the potential suitable region for Emex australis in China under different climate scenarios in different periods |
如图5所示，对不同气候情景下南方三棘果适生区的空间格局进行对比分析，结果表明，2040年，未来4种情景模式下的适生区都呈现轻微缩小趋势，广东和广西的部分低适生区变为非适生区，云南的高适生区变为中适生区，海南的高、中、低适生区都呈减少趋势。从单个情景模式下看也有一些地区变得更加适应南方三棘果分布，SSP1-2.6气候情景模式下西藏出现局部高适生区，四川的局部低适生区转变为中适生区；SSP2-4.5气候情景模式下福建的一些非适生区变为低适生区；SSP3-7.0气候情景模式下四川部分低适生区转变为中适生区，局部地区变为高适生区，西藏出现局部高适生区；SSP5-8.5气候情景模式下，四川和西藏的一些非适生区变为低适生区，其中，西藏出现局部高适生区。这些变化区域为南方三棘果分布的敏感区。
|
图 5 4种气候情景模式下南方三棘果在中国的适生区 Fig. 5 The suitable growth reqion of Emex australis in China under four climate scenarios |
本模型的AUC测试值 > 0.9，模型结果优秀可靠。最冷季度平均气温(Bio11)、等温性(Bio3)和温度季节变化标准差(Bio4)是影响南方三棘果地理分布的主要环境因子。正则化训练增益和测试增益以及贡献率都表明最冷季度平均气温(Bio11)对南方三棘果分布影响最大，贡献率为27.7%。环境因子响应曲线表明，当最冷季度平均气温为9.35~12.76 ℃时，南方三棘果发生概率≥0.5。南方三棘果的适宜分布区总面积约为69.02万km 2，适生区域主要集中于云南、广东、广西和福建地区，其中高、中适生区主要集中于云南。
通过对2022—2040年国家气候中心气候系统模式(BCC-CSM2-MR)的4种未来气候情景下的南方三棘果潜在分布区研究发现，总体上适生区呈现轻微减少趋势。四川、西藏、广东和广西属于南方三棘果分布的敏感区，易随气候变化而导致适生区等级以及适生面积发生变化。广东、广西、海南和云南的一些适生区在未来一段时间会出现转变为非适生区或降低适生等级的情况。福建部分地区在SSP2-4.5情景模式下会更适宜南方三棘果生存，西藏和四川部分地区在SSP1-2.6、SSP3-7.0和SSP5-8.5情景下会更适宜南方三棘果生存。
其中，西藏的适生区尤为特殊，主要集中在藏南地区的错那县和墨脱县以及朗县最南部。这一部分地区由于地形原因导致气候复杂，高原气候和亚热带、热带季风气候同时存在，加上地处边疆，适生区西边与不丹接壤，南边与印度接壤，导致藏南外来物种入侵严重，值得相关部门进一步关注。中国科学院植物研究所的植物物种信息系统(
南方三棘果是具有严重危害的检疫杂草，对农业生产会造成巨大破坏[35-36]。云南是我国葡萄十大产地之一，南方三棘果在我国的高、中适生区主要集中于云南，南方三棘果一旦定殖于葡萄园，可能会有部分种子在机械化收割时与葡萄一起收割，影响葡萄的品相和食品安全性[37]。藏南是西藏的主要农区，藏南地区生产的小麦是西藏的主要粮食作物之一，研究表明南方三棘果茎的水提取物和植株残留物会抑制小麦种子萌发和早期幼苗的生长，降低小麦产量[38-39]。油菜是广西、广东和福建的主要油料作物，在农田里南方三棘果与油料作物存在种间竞争关系[40]。
外来入侵杂草传入与人员和货物的流动密切相关，作为国与国之间管控人员和货物流通的重要节点，海关和动植物检验检疫局实施的植物检疫是防止外来杂草传入的最重要方式。在适生区域内以及邻近适生区域的海关和动植物检验检疫局应当加大对南方三棘果的检疫力度，防止南方三棘果的传入[41]。尤其值得注意的是，广东、广西和福建的海关每年都有大量的货物进口，加上这些地区存在大面积的南方三棘果适生区，海关和动植物检验检疫局的植物疫情数据也表明在广东和福建的海关有多起南方三棘果的截获记录，因此，这些地区是南方三棘果传入的高风险地区。
有害杂草入侵的定殖和扩散阶段在适生区范围内进行，其中，外部气候环境的变化是影响其扩散速度的重要因素，试验表明最冷季度平均温度是影响南方三棘果分布的最重要环境因子，在适生区气候适宜其生存的重要年份要加大南方三棘果的监测，有效防止其定殖、扩散，以免影响农作物的质量和产量[42]。生物入侵杂草的分布与气候环境、地形因素、土壤类型以及植被类型等多种因素相关，在使用MaxEnt生态位模型预测入侵杂草潜在分布区时，目前尚未有统一的影响因子选择。朱煜建等[43]在预测疣果匙荠Bunias orientalis适生区时选择了气候环境、地形因素以及土壤类型的环境变量因子。郭燕青等[44]在预测假臭草Praxelis clematidea时仅选择了气候因素有关因子。本研究选用与环境和地形有关的变量因子，缺少土壤类型、物种类型等与南方三棘果分布有关的其他变量，这可能对适生区的预测带来一定的偏差。伴随着我国与世界各国日益频繁的贸易往来，生物入侵难以避免。近些年假高粱Sorghum halepense和豚草Ambrosia artemisiifolia等入侵生物对我国农业生产和生态安全造成了巨大破坏。采用MaxEnt生态位模型将信息学技术与生物的生态学特性结合，对各种入侵杂草的潜在分布区进行预测，对我国的口岸检疫和入侵杂草的防控具有指导意义。
| [1] |
陈凤新, 蒙彦良, 杜杰, 等. 中国外来生物入侵与社会经济因素多元回归分析[J]. 植物检疫, 2021, 35(4): 8-14. DOI:10.19662/j.cnki.issn1005-2755.2021.04.002 (  0) |
| [2] |
李瑞法, 柳之光, 楚伟, 等. 我国进境植物疫情截获量的时序特征及预测[J]. 植物检疫, 2016, 30(4): 1-5. DOI:10.19662/j.cnki.issn1005-2755.2016.04.001 ( 0) |
| [3] |
王昕, 戴良英. 黄埔口岸进境粮食截获疫情分析[J]. 植物检疫, 2016, 30(6): 80-82. DOI:10.19662/j.cnki.issn1005-2755.2016.06.023 ( 0) |
| [4] |
伏建国, 高振兴, 李敏. 江苏口岸首次从进境加拿大油菜籽中截获检疫性杂草南方三棘果[J]. 植物检疫, 2010, 24(4): 68. ( 0) |
| [5] |
MENDONCA L P, OLIVERIA E E, ANDERAZZA F, et al. Host potential and adaptive responses of Drosophila suzukii (Diptera: Drosophilidae) to barbados cherries
[J]. Journal of Economic Entomology, 2019, 112(6): 3002-3006. ( 0) |
| [6] |
王瑞, 万方浩. 入侵植物银毛龙葵在中国的适生区预测与早期监测预警[J]. 生态学杂志, 2016, 35(7): 1697-1703. ( 0) |
| [7] |
叶建峰, 唐丽霞, 吴先健, 等. 东莞口岸从进境燕麦粒中截获检疫性杂草南方三棘果[J]. 植物检疫, 2010, 24(2): 42. DOI:10.19662/j.cnki.issn1005-2755.2010.02.015 ( 0) |
| [8] |
HERNANDEZ P A, GRAHAM C H, MASTER L L, et al. The effect of sample size and species characteristics on performance of different species distribution modeling methods[J]. Ecography, 2010, 29(5): 773-785. ( 0) |
| [9] |
KHANGHAH S S, MOAMERI M, GHORBANI A, et al. Modeling potential habitats and predicting habitat connectivity for Leucanthemum vulgare Lam. in northwestern rangelands of Iran
[J]. Environmental Monitoring and Assessment, 2022, 194(2): 1-16. ( 0) |
| [10] |
王运生. 生态位模型在外来入侵物种风险评估中的应用研究[D]. 长沙: 湖南农业大学, 2007.
( 0) |
| [11] |
LUO C, XU W H, ZHOU Z X, et al. Habitat prediction for forest musk deer (Moschus berezovskii) in Qinling mountain range based on niche model
[J]. Acta Ecologica Sinica, 2011, 31(5): 1221-1229. ( 0) |
| [12] |
房锋, 张朝贤, 黄红娟, 等. 基于MaxEnt的麦田恶性杂草节节麦的潜在分布区预测[J]. 草业学报, 2013, 22(2): 62-70. DOI:10.11686/cyxb20130208 ( 0) |
| [13] |
柳晓燕, 李俊生, 赵彩云, 等. 基于MAXENT模型和ArcGIS预测豚草在中国的潜在适生区[J]. 植物保护学报, 2016, 43(6): 1041-1048. ( 0) |
| [14] |
唐兴港, 袁颖丹, 张金池. 气候变化对杉木适生区和生态位的影响[J]. 植物研究, 2022, 42(1): 151-160. DOI:10.7525/j.issn.1673-5102.2022.01.016 ( 0) |
| [15] |
李安, 李良涛, 高萌萌, 等. 基于MaxEnt模型和气候变化情景入侵种黄顶菊在中国的分布区预测[J]. 农学学报, 2020, 10(1): 60-67. DOI:10.11923/j.issn.2095-4050.cjas20190700109 ( 0) |
| [16] |
姜大膀, 富元海. 2 ℃全球变暖背景下中国未来气候变化预估[J]. 大气科学, 2012, 36(2): 234-246. DOI:10.3878/j.issn.1006-9895.2011.11074 ( 0) |
| [17] |
JAMAL Z A, ABOU-SHAARA H F, QAMER S, et al. Future expansion of small hive beetles, Aethina tumida, towards North Africa and South Europe based on temperature factors using maximum entropy algorithm
[J]. Journal of King Saud University-Science, 2021, 33(1): 1-17. DOI:10.1016/j.jksues.2019.09.005 ( 0) |
| [18] |
KONG F, TANG L, HE H, et al. Assessing the impact of climate change on the distribution of Osmanthus fragrans using MaxEnt
[J]. Environmental Science and Pollution Research, 2021, 28(26): 34655-34663. DOI:10.1007/s11356-021-13121-3 ( 0) |
| [19] |
辛晓歌, 吴统文, 张洁, 等. BCC模式及其开展的CMIP6试验介绍[J]. 气候变化研究进展, 2019, 15(5): 533-539. ( 0) |
| [20] |
王翀, 林慧龙, 何兰, 等. 紫茎泽兰潜在分布对气候变化响应的研究[J]. 草业学报, 2014, 23(4): 20-30. DOI:10.11686/cyxb20140403 ( 0) |
| [21] |
SONG R, MA Y, HU Z, et al. MaxEnt modeling of Dermacentor marginatus (Acari: Ixodidae) distribution in Xinjiang, China
[J]. Journal of Medical Entomology, 2020, 57(5): 1659-1667. DOI:10.1093/jme/tjaa063 ( 0) |
| [22] |
ZHANG M G, SLIK J, MA K P. Using species distribution modeling to delineate the botanical richness patterns and phytogeographical regions of China[J]. Scientific Reports, 2016, 6(1): 22400. doi: 10.1038/srep22400.
( 0) |
| [23] |
GRAHAM M H. Confronting multicollinearity in ecological multiple regression[J]. Ecology, 2003, 84(11): 2809-2815. DOI:10.1890/02-3114 ( 0) |
| [24] |
张路. 基于MAXENT模型预测齿裂大戟在中国的潜在分布区[J]. 生物安全学报, 2015, 24(3): 194-200. DOI:10.3969/j.issn.2095-1787.2015.03.003 ( 0) |
| [25] |
任子春. 基于MaxEnt模型在全球变暖条件下菟丝子属全寄生植物及其寄主的潜在分布预测[D]. 临汾: 山西师范大学, 2020.
( 0) |
| [26] |
郭水良, 高平磊, 娄玉霞. 应用MaxEnt模型预测检疫性杂草毒莴苣在我国的潜分布范围[J]. 上海交通大学学报(农业科学版), 2011, 29(5): 15-19. ( 0) |
| [27] |
WILLIAMS P, YA T, JIAN Y, et al. The bumblebees of Sichuan (Hymenoptera: Apidae, Bombini)[J]. Systematics and Biodiversity, 2009, 7(2): 101-187. DOI:10.1017/S1477200008002843 ( 0) |
| [28] |
叶兴状, 张明珠, 赖文峰, 等. 基于MaxEnt优化模型的闽楠潜在适宜分布预测[J]. 生态学报, 2021, 41(20): 8135-8144. ( 0) |
| [29] |
ELITH J, GRAHAM C H, ANDERSON R P, et al. Novel methods improve prediction of species distributions from occurrence data[J]. Ecography, 2010, 29(2): 129-151. ( 0) |
| [30] |
张海娟, 陈勇, 黄烈健, 等. 基于生态位模型的薇甘菊在中国适生区的预测[J]. 农业工程学报, 2011, 27(S1): 413-418. ( 0) |
| [31] |
朱耿平, 范靖宇, 王梦琳, 等. ROC曲线形状在生态位模型评价中的重要性——以美国白蛾为例[J]. 生物安全学报, 2017, 26(3): 184-190. DOI:10.3969/j.issn.2095-1787.2017.03.002 ( 0) |
| [32] |
PEARCE J L, BOYCE M S. Modelling distribution and abundance with presence-only data[J]. Journal of Applied Ecology, 2006, 43(3): 405-412. DOI:10.1111/j.1365-2664.2005.01112.x ( 0) |
| [33] |
KOZAK K H, GRAHAM C H, WIENS J J. Integrating GIS-based environmental data into evolutionary biology[J]. Trends in Ecology & Evolution, 2008, 23(3): 141-148. ( 0) |
| [34] |
陈星, 赵韦, 包海龙. 西藏山南地区采集纪行——错那县[J]. 生命世界, 2018(3): 94-95. ( 0) |
| [35] |
张玉环, 郝建华, 吴海荣, 等. 外来入侵植物胜红蓟的胚胎学观察及繁殖系统研究[J]. 植物科学学报, 2020, 38(2): 162-172. DOI:10.11913/PSJ.2095-0837.2020.20162 ( 0) |
| [36] |
王桔红, 史生晶, 陈文, 等. 鬼针草与含羞草化感作用及其入侵性的研究[J]. 草业学报, 2020, 29(4): 81-91. DOI:10.11686/cyxb2019480 ( 0) |
| [37] |
POHLNER D. Emex australis and dried vine fruit production in Sunraysia
[J]. Plant Protection Quarterly, 1996, 11(4): 150-153. ( 0) |
| [38] |
JAVAID M M, TANVEER A, ALI H H, et al. Wheat yield loss in a two species competition with Emex australis and Emex spinosa
[J]. Planta Daninha, 2016, 34: 35-46. DOI:10.1590/S0100-83582016340100004 ( 0) |
| [39] |
ABBAS R N, TANVEER A, ALI A, et al. Effects of Emex australis Steinh on germination and early seedling growth of wheat (Triticum aestivum L. )
[J]. Allelopathy Journal, 2010, 25(2): 513-520. ( 0) |
| [40] |
SCOTT J K, BOWRAN D G. Workshop to identify research priorities for Emex species
[J]. Plant Protection Quarterly, 1996, 11(4): 175-176. ( 0) |
| [41] |
PUCHAŁKA R, DYDERSKI M K, VÍTKOVÁ M, et al. Black locust (Robinia pseudoacacia L. ) range contraction and expansion in Europe under changing climate
[J]. Global Change Biology, 2021, 27(8): 1587-1600. DOI:10.1111/gcb.15486 ( 0) |
| [42] |
梁忆冰. 植物检疫对外来有害生物入侵的防御作用[J]. 植物保护, 2002(2): 45-47. DOI:10.3969/j.issn.0529-1542.2002.02.018 ( 0) |
| [43] |
朱煜建, 周伟光, 张昊, 等. 刺蒺藜草在我国的适生区预测[J]. 植物检疫, 2022, 36(1): 71-76. DOI:10.19662/j.cnki.issn1005-2755.2021.00.028 ( 0) |
| [44] |
郭燕青, 史梦竹, 李建宇, 等. 基于MaxEnt模型的假臭草潜在分布区预测[J]. 热带亚热带植物学报, 2019, 27(3): 250-260. DOI:10.11926/jtsb.3977 ( 0) |