Prediction and analysis of the global suitability of Lymantria dispar based on MaxEnt
-
摘要:
目的 舞毒蛾是一种食叶性的国际性检疫害虫,给世界上许多国家和地区造成了严重的经济损失。该研究筛选出限制舞毒蛾分布的环境变量,利用MaxEnt软件预测舞毒蛾当前及未来气候条件下的全球适生区范围,明确舞毒蛾在不同气候条件下的适生区变化。 方法 利用ArcGIS软件设置缓冲区筛选舞毒蛾在全球的分布点数据;利用MaxEnt、SPSS和ArcGIS软件根据环境变量贡献率、刀切法和变量相关性分析对生物气候变量、月总降水量、月平均最高温度和月平均最低温度4种环境变量进行筛选;利用R软件计算调控倍频和特征组合等因子调整MaxEnt模型参数;利用MaxEnt模型预测当前和未来不同情境条件下舞毒蛾全球适生区的分布范围。 结果 经过缓冲区筛选得到734个舞毒蛾的分布点数据;MaxEnt模型结果中,测试遗漏率与理论遗漏率吻合度高,而且模型AUC值为0.940;MaxEnt模型预测当前条件下舞毒蛾在全球的高、中度适生区主要集中在欧洲的大部分地区,北美洲中东部,亚洲的东西部,而非洲、大洋洲和南美洲分布较少。此外,舞毒蛾在未来气候条件下北半球适生区的边界向北偏移,北美洲以及欧亚大陆的高、中度适生区的面积扩增明显。 结论 舞毒蛾的分布受多种环境变量影响,并且温度和降水与舞毒蛾的特定发育阶段相吻合。MaxEnt模型的预测结果能够反映舞毒蛾在全球的分布特征。该研究可为防治舞毒蛾和制定相关检疫措施提供理论依据。 Abstract:Objective Lymantria dispar (Linnaeus) is a leaf-eating international quarantine pest, which has caused serious economic losses to many countries and regions in the world. In this study, environmental variables limiting the distribution of Lymantria dispar were screened out, and MaxEnt software was used to predict the global adaptation range of Lymantria dispar under current and future climatic conditions, so as to clarify the adaptation range of Lymantria dispar under different climatic conditions. Method We used ArcGIS software to set up buffer zones to screen the global distribution data of Lymantria dispar. We used MaxEnt, SPSS and ArcGIS software to screen four kinds of environmental variables, i.e. biological climate variable, monthly total precipitation, monthly mean maximum temperature and monthly mean minimum temperature, according to the contribution rate of environmental variables, knife cutting method and variable correlation analysis. We used R software to calculate the regularization multiplier and feature classes and other factors to adjust the parameters of MaxEnt model. We used MaxEnt model to predict the global distribution range of Lymantria dispar under current and future climatic conditions. Result The distribution data of 734 Lymantria dispar were obtained by buffer screening. In the MaxEnt model results, the test omission rate was in high agreement with the theoretical omission rate, and the AUC value of the model was 0.940. MaxEnt model predicted that Lymantria dispar was mainly distributed in most parts of Europe, central and eastern North America, western and eastern Asia, and less in Africa, Oceania and South America under current conditions. In addition, Lymantria dispar boundaries in the northern hemisphere will shift to the north under future climatic conditions, and the areas of high and moderate adaptability in North America and Eurasia will expand significantly. Conclusion The distribution of Lymantria dispar is influenced by a variety of environmental variables and temperature and precipitation are consistent with the specific development stage of Lymantria dispar. The predicted results of MaxEnt model can reflect the distribution characteristics of Lymantria dispar in the world. This study can provide theoretical basis for prevention and control of Lymantria dispar and establishment of relevant quarantine measures. -
Key words:
- Lymantria dispar /
- MaxEnt /
- adaptive area prediction /
- CMIP6
-
表 1 MaxEnt参数设置值
Table 1. Parameter setting values of MaxEnt
选项
Option参数 Parameter 默认值
Default value设置值
Setting value随机选取测试集比例
Random text percentage0 25 重复训练次数
Number of repetitions1 10 最大重复次数
Maximum number of repetitions500 5000 表 2 影响舞毒蛾分布的环境变量及其贡献率
Table 2. Environmental variables affecting Lymantria dispar distribution and their contribution rates
编号 No. 气候变量 Climatic variable 贡献率 Contribution rate/% Bio2 温差月均值 Monthly mean temperature difference 3.7 Bio9 最干季度平均温度 Mean temperature of driest quarter 1.2 Bio10 最暖季度平均温度 Mean temperature of warmest quarter 37.9 Bio14 最干月降水量 Precipitation of driest month 39.9 Bio16 最湿季度降水量 Precipitation of wettest quarter 3.8 Prec3 3月平均降水量 Average precipitation in March 13.5 表 3 多种气候情景的AUC值
Table 3. AUC values under various climatic scenarios
气候情景
Climate change scenario年份
YearAUC值
AUC value当前
Current当前 Current 0.940 低强迫情景 ssp126
Low compulsion scenario ssp1262021—2040 0.940 2041—2060 0.944 2061—2080 0.942 2081—2100 0.943 中等强迫情景ssp245
Moderate compulsion scenario ssp2452021—2040 0.943 2041—2060 0.943 2061—2080 0.944 2081—2100 0.944 中等至高等强迫情景ssp370
Moderate to high compulsive scenario ssp3702021—2040 0.942 2041—2060 0.941 2061—2080 0.942 2081—2100 0.943 高等强迫情景ssp585
High compulsion scenario ssp5852021—2040 0.945 2041—2060 0.941 2061—2080 0.944 2081—2100 0.943 -
[1] 张丽茹. 舞毒蛾的生物学特性及综合防治技术[J]. 现代农业科技, 2020(6):116−117. doi: 10.3969/j.issn.1007-5739.2020.06.072Zhang L R. Biological characteristics and comprehensive control techniques of Lymantira dispar[J]. Modern Agricultural Science and Technology, 2020(6): 116−117. doi: 10.3969/j.issn.1007-5739.2020.06.072 [2] Keena M A, Côté M, Grinberg P S, et al. World distribution of female flight and genetic variation in Lymantria dispar (Lepidoptera: Lymantriidae)[J]. Environmental Entomology, 2008, 37(3): 636−649. doi: 10.1603/0046-225X(2008)37[636:WDOFFA]2.0.CO;2 [3] 钱路, 安榆林, 徐梅, 等. 舞毒蛾不同地理种群基于AFLP分子标记的遗传分析[J]. 林业科学, 2011, 47(10):104−110. doi: 10.11707/j.1001-7488.20111016Qian L, An Y L, Xu M, et al. AFLP analysis of different geographic populations of the gypsy moth, Lymantria dispar (Lepidoptera: Lymantriidae)[J]. Scientia Silvae Sinicae, 2011, 47(10): 104−110. doi: 10.11707/j.1001-7488.20111016 [4] Luque G M, Bellard C, Bertelsmeier C, et al. The 100th of the world’s worst invasive alien species[J]. Biological Invasions, 2014, 16(5): 981−985. doi: 10.1007/s10530-013-0561-5 [5] 谵运清, 安输林, 刘翔, 等. 由亚洲型舞毒蛾检疫问题引发的几点思考[J]. 植物检疫, 2012, 26(4):76−78.Zhan Y Q, An Y L, Liu X, et al. Some thoughts on quarantine of Gypsy moth in Asia[J]. Plant quarantine, 2012, 26(4): 76−78. [6] Srivastava V, Griess V C, Keena M A. Assessing the potential distribution of Asian gypsy moth in Canada: a comparison of two methodological approaches[J]. Sentific Reports, 2020, 10(22): 1−10. [7] Phillips S J, Anderson R P, Schapire R E. Maximum entropy modeling of species geographic distributions[J]. Ecological Modelling, 2006, 190(3/4): 231−259. [8] Merow C, Smith M J, Silander J A. A practical guide to MaxEnt for modeling species’ distributions: what it does, and why inputs and settings matter[J]. Ecography, 2013, 36(10): 1058−1069. doi: 10.1111/j.1600-0587.2013.07872.x [9] Phillips S J, Dudík M. Modeling of species distributions with MaxEnt: new extensions and a comprehensive evaluation[J]. Ecography, 2008, 31(2): 161−175. doi: 10.1111/j.0906-7590.2008.5203.x [10] Elith J, Phillips S J, Hastie T, et al. A statistical explanation of MaxEnt for ecologists[J]. Diversity & Distributions, 2011, 17(1): 43−57. [11] 王运生, 谢丙炎, 万方浩, 等. ROC曲线分析在评价入侵物种分布模型中的应用[J]. 生物多样性, 2007(4):365−372. doi: 10.3321/j.issn:1005-0094.2007.04.005Wang Y S, Xie B Y, Wan F H, et al. Application of ROC curve analysis in evaluating the performance of alien species’ potential distribution models[J]. Biodiversity Science, 2007(4): 365−372. doi: 10.3321/j.issn:1005-0094.2007.04.005 [12] 李一琳, 丁长青. 基于GIS和MaxEnt技术对濒危物种褐马鸡的保护空缺分析[J]. 北京林业大学学报, 2016, 38(11):34−41.Li Y L, Ding C Q. Reserve gap analysis of endangered brown eared pheasant (Crossoptilon mantchuricum) through GIS and MaxEnt technology[J]. Journal of Beijing Forest University, 2016, 38(11): 34−41. [13] 赵佳强, 石娟. 基于新型最大熵模型预测刺槐叶瘿蚊(双翅目: 瘿蚊科)在中国的适生区[J]. 林业科学, 2019, 55(2):118−127.Zhao J Q, Shi J. Prediction of the potential geographical distribution of Obolodiplosis robiniae (Diptera: Cecidomyiidae) in China based on a novel Maximum Entropy model[J]. Scientia Silvae Sinicae, 2019, 55(2): 118−127. [14] 黄梦伊, 赵佳强, 石娟. 基于MaxEnt对桉树枝瘿姬小蜂在中国发生趋势的预测[J]. 北京林业大学学报, 2020, 42(11):64−71.Huang M Y, Zhao J Q, Shi J. Predicting occurrence tendency of Leptocybe invasa in China based on MaxEnt[J]. Journal of Beijing Forestry University, 2020, 42(11): 64−71. [15] Fang W, Duo W, Ge G, et al. Potential distributions of the invasive barnacle scale Ceroplastes cirripediformis (Hemiptera: Coccidae) under climate change and implications for its management[J]. Journal of Economic Entomology, 2021, 114(1): 82−89. doi: 10.1093/jee/toaa245 [16] 张春华, 和菊, 孙永玉, 等. 基于MaxEnt模型的紫椿适生区预测[J]. 北京林业大学学报, 2017, 39(8):33−41.Zhang C H, He J, Sun Y Y, et al. Distributional change in suitable areas for Toona sureni based on MaxEnt model[J]. Journal of Beijing Forest University, 2017, 39(8): 33−41. [17] 吕汝丹, 何健, 刘慧杰, 等. 羽叶铁线莲的分布区与生态位模型分析[J]. 北京林业大学学报, 2019, 41(2):70−79.Lü R D, He J, Liu H J, et al. Distribution and niche modeling analysis of Clematis pinnata[J]. Journal of Beijing Forestry University, 2019, 41(2): 70−79. [18] Inoue M N, Suzuki-Ohno Y, Haga Y, et al. Population dynamics and geographical distribution of the gypsy moth, Lymantria dispar, in Japan[J]. Forest Ecology and Management, 2019, 434: 154−164. doi: 10.1016/j.foreco.2018.12.022 [19] 姜彤, 吕嫣冉, 黄金龙, 等. CMIP6模式新情景(SSP-RCP)概述及其在淮河流域的应用[J]. 气象科技进展, 2020, 10(5):102−109. doi: 10.3969/j.issn.2095-1973.2020.05.016Jiang T, Lü Y R, Huang J L, et al. New scenarios of CMIP6 model (SSP-RCP) and its application in the Huaihe river basin[J]. Advances in Meteorological Science and Technology, 2020, 10(5): 102−109. doi: 10.3969/j.issn.2095-1973.2020.05.016 [20] 张丽霞, 陈晓龙, 辛晓歌. CMIP6情景模式比较计划(ScenarioMIP)概况与评述[J]. 气候变化研究进展, 2019, 15(5):519−525. doi: 10.12006/j.issn.1673-1719.2019.082Zhang L X, Chen X L, Xin X G. Short commentary on CMIP6 scenario model intercomparison project (ScenarioMIP)[J]. Climate Change Research, 2019, 15(5): 519−525. doi: 10.12006/j.issn.1673-1719.2019.082 [21] Sillero N. What does ecological modelling model? A proposed classification of ecological niche models based on their underlying methods[J]. Ecological Modelling, 2011, 222(8): 1343−1346. doi: 10.1016/j.ecolmodel.2011.01.018 [22] 朱耿平, 刘强, 高玉葆. 提高生态位模型转移能力来模拟入侵物种的潜在分布[J]. 生物多样性, 2014, 22(2):223−230. doi: 10.3724/SP.J.1003.2014.08178Zhu G P, Liu Q, Gao Y B. Improving ecological niche model transferability to predict the potential distribution of invasive exotic species[J]. Biodiversity Science, 2014, 22(2): 223−230. doi: 10.3724/SP.J.1003.2014.08178 [23] Phillips S J, Anderson R P, Dudík M, et al. Opening the black box: an open-source release of Maxent[J]. Ecography, 2017, 40(7): 887−893. doi: 10.1111/ecog.03049 [24] 刘洋, 石娟. 气候变化背景下埃及吹绵蚧在中国的适生区预测[J]. 植物保护, 2020, 46(1):108−117.Liu Y, Shi J. Prediction of potential geographical distribution of Icerya aegyptiaca in China under climate change[J]. Plant Protection, 2020, 46(1): 108−117. [25] Muscarella R, Galante P J, Soley-Guardia M, et al. ENMeval: an R package for conducting spatially independent evaluations and estimating optimal model complexity for Maxent ecological niche models[J]. Methods in Ecology & Evolution, 2015, 5(11): 1198−1205. [26] Shcheglovitova M, Anderson R P. Estimating optimal complexity for ecological niche models: a jackknife approach for species with small sample sizes[J]. Ecological Modelling, 2013, 269: 9−17. [27] Zhong Q, Zhang J E, Ditommaso A, et al. Predicting invasions of Wedelia trilobata (L.) Hitchc. with MaxEnt and GARP models[J]. Journal of Plant Research, 2015, 128: 763−775. doi: 10.1007/s10265-015-0738-3 [28] Pearson R G, Raxworthy C J, Nakamura M, et al. Original article: predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar[J]. Journal of Biogeography, 2007, 34(1): 102−117. [29] Liu Y, Shi J. Predicting the potential global geographical distribution of two icerya species under climate change[J]. Forests, 2020, 11(6): 1−19. [30] 赵晓冏, 巩娟霄, 赵莎莎, 等. 样本量及其空间分布对物种分布模型的影响[J]. 兰州大学学报: 自然科学版, 2018, 54(2):70−77.Zhao X J, Gong J X, Zhao S S, et al. Impact of sample size and spatial distribution on species distribution model[J]. Journal of Lanzhou University: Natural Sciences, 2018, 54(2): 70−77. [31] Morey A C, Venette R C. Minimizing risk and maximizing spatial transferability: challenges in constructing a useful model of potential suitability for an invasive insect[J]. Annals of the Entomological Society of America, 2020, 113(2): 100−113. doi: 10.1093/aesa/saz049 [32] 闻连仁, 赵春芝, 张洪彦, 等. 舞毒蛾生物学特性观察[J]. 吉林林业科技, 1989(2):34−37.Wen L R, Zhao C Z, Zhang H Y, et al. Observation on biological characteristics of Lymantria dispar[J]. Journal of Jilin Forestry Science and Technology, 1989(2): 34−37. [33] 卢小雨, 陈洪俊, 陈乃中, 等. 亚洲型舞毒蛾在北美的适生性[J]. 昆虫知识, 2009, 46(3):398−402, 495.Lu X Y, Chen H J, Chen N Z, et al. Potential geographic distribution of Asian gypsy moth, Lymantria dispar, in north America[J]. Chinese Bulletin of Entomology, 2009, 46(3): 398−402, 495. [34] Heikkinen R K, Luoto M, Araújo M B, et al. Methods and uncertainties in bioclimatic envelope modelling under climate change[J]. Progress in Physical Geography, 2006, 30(6): 751−777. doi: 10.1177/0309133306071957 [35] Sinclair S J, White M D, Newell G R. How useful are species distribution models for managing biodiversity under future climates?[J]. Ecology & Society, 2010, 15(1): 175−183. -