用微信扫码二维码

分享至好友和朋友圈

WE ARE COMMITTED TO REPORTING THE LATEST FORESTRY ACADEMIC ACHIEVEMENTS

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于MaxEnt的松墨天牛在中国的潜在分布及其对气候变化的响应

刘勇 李祥乾 李阳娣 杨伟 杨桦

刘勇, 李祥乾, 李阳娣, 等. 基于MaxEnt的松墨天牛在中国的潜在分布及其对气候变化的响应[J/OL]. 四川林业科技, 2024, 45(2)[2023-11-30] doi: 10.12172/202306120001
引用本文: 刘勇, 李祥乾, 李阳娣, 等. 基于MaxEnt的松墨天牛在中国的潜在分布及其对气候变化的响应[J/OL]. 四川林业科技, 2024, 45(2)[2023-11-30] doi: 10.12172/202306120001
LIU Y, LI X Q, LI Y D, et al. MaxEnt modeling for predicting the potential distribution of monochamus alternatus hope (coleoptera: cerambycidae) in china[J/OL]. Journal of Sichuan Forestry Science and Technology, 2024, 45(2)[2023-11-30] doi: 10.12172/202306120001
Citation: LIU Y, LI X Q, LI Y D, et al. MaxEnt modeling for predicting the potential distribution of monochamus alternatus hope (coleoptera: cerambycidae) in china[J/OL]. Journal of Sichuan Forestry Science and Technology, 2024, 45(2)[2023-11-30] doi: 10.12172/202306120001

基于MaxEnt的松墨天牛在中国的潜在分布及其对气候变化的响应


doi: 10.12172/202306120001
详细信息
    作者简介:

    刘勇(1971—),男,四川兴文人,工程师,学士,主要从事森林保护,E-mail:695189352@qq.com

    通讯作者: 杨桦(1982—),男,副教授,博士,主要从事林业害虫综合防治研究,E-mail: yanghua151017@163.com
  • 基金项目:  四川省科技厅国际合作项目(2021YFH0160)

MaxEnt modeling for predicting the potential distribution of Monochamus alternatus Hope (Coleoptera: Cerambycidae) in China

More Information
    Corresponding author: yanghua151017@163.com
  • 摘要: 由于世界气候变暖而导致的害虫空间分布格局的巨大变化,增加了害虫预防难度。通过MaxEnt模型,结合GIS方法,模拟松墨天牛在中国的潜在分布区,预测、对比、分析松墨天牛在未来3种气候情景(RCP2.6、RCP4.5和RCP8.5)下分布范围和空间布局,及其中心节点移动的规律。结果表明:所有AUC值均超过0.95,准确性也处于“极好”的状态。在当前气候环境下,松墨天牛在中国的适生区范围为90°~125°E,20°~41°N,由西向东呈连续分布,集中分布于我国东南部;在未来气候条件下,中适生区面积减少,低适生区、高适生区面积增加,并向高纬度、高海拔地区扩散,RCP8.5扩散范围最大。适生区中心点向西北、东北方向迁移,迁移趋势在高浓度排放情景下 (RCP 8.5) 响应最为剧烈。为松墨天牛的合理区划提供了有效手段, 为应对气候变暖提供了理论依据。
  • 图  1  环境因子刀切法检验结果

    Fig.  1  Environmental factor test results by knife method

    图  2  最暖季降雨量(bio18)

    Fig.  2  Rainfall in the warmest season(bio 18)

    图  3  当前气候条件下的ROC曲线图

    Fig.  3  ROC curve for current climatic conditions

    图  4  未来不同气候情境下MaxEnt模拟结果

    Fig.  4  MaxEnt simulation results under different climate scenarios in the future

    图  5  松墨天牛高适生区分布

    Fig.  5  Distribution of high suitability area of Monochamus alternatus

    图  6  气候变化情景(RCP2.6、RCP4.5和RCP8.5)下松墨天牛高适生区中心点位置及其转移

    Fig.  6  The location and transfer of mean center of highly suitable area of Monochamus alternatus under climate change scenarios

    图  7  气候变化情景(RCP2.6、RCP4.5和RCP8.5)下松墨天牛总适生区中心点位置及其转移

    Fig.  7  The location and transfer of mean center of totally suitable area of Monochamus alternatus under climate change scenarios

    表  1  生物气候特征表

    Tab.  1  List of bioclimatic characteristics

    生物气候特征bioclimatic characteristics 缩写abbreviation
    年平均气温 bio1
    平均气温日较差 bio2
    等温性 bio3
    气温季节性变动系数(标准差*100) bio4
    最热月份最高温度 bio5
    最冷月份最低温度 bio6
    气温年较差 bio7
    最湿季度平均温度 bio8
    最干季度平均温度 bio9
    最暖季度平均温度 bio10
    最冷季度平均温度 bio11
    年降水量 bio12
    最湿月份降水量 bio13
    最干月份降水量 bio14
    降水量季节性变化 bio15
    最干季度降水量 bio16
    最湿季度降水量 bio17
    最暖季度降水量 bio18
    最冷季度降水量 bio19
    最低温度 tmin
    最高温度 tmax
    平均温度 tavg
    降水量 prec
    下载: 导出CSV

    表  2  AUG评价标准表

    Tab.  2  Evaluation criterion of AUC

    AUC值范围Range of AUC value评价标准Evaluation criterion
    0.5≤AUC<0.6失败
    0.6≤AUC<0.7不佳
    0.7≤AUC<0.8较好
    0.8≤AUC<0.9
    ≥0.9极好
    下载: 导出CSV

    表  3  当前气候条件下松墨天牛的潜在适生区预测

    Tab.  3  Prediction of potential habitat area of Monochamus alternatus under current climatic conditions

    省份名称
    Province
    非适生区
    面积
    unsuitable area(km2)<5%
    低适生区
    面积
    lowly suitable
    area (km2)
    5-33%
    中适生区
    面积
    moderately suitable area
    (km2)33-67%
    高适生区
    面积
    highly suitable
    area (km2)
    ≥67%
    总面积
    Total(km2)
    总适生
    区面积
    Total suitable area
    (km2)
    占全国高适生
    区百分比
    Percentage of highly
    suitable areas in China
    占全国总适生
    区百分比
    Percentage of
    total suitable
    areas in China
    黑龙江Heilongjiang543697.920.000.000.005436980.000.00000.000
    内蒙古Inner Mongolia1290607.660.000.000.0012906080.000.00000.000
    新疆Xinjiang1755208.360.000.000.0017552080.000.00000.000
    吉林Jilin210277.781857.64329.860.0021246521880.00000.010
    辽宁Liaoning121267.360.0030746.533541.670.00155556342880.00000.220
    甘肃Gansu408628.484861.111510.42243.0641524366150.000590.016
    河北Hebei158333.3437864.58277.780.00196476381420.00000.194
    北京Beijing9565.977673.610.000.001724076740.00000.445
    山西Shanxi155503.474236.110.000.0015974042360.00000.267
    天津Tianjin3315.978767.360.000.001208387670.00000.726
    陕西Shaanxi149392.3630659.7210468.7513194.44203715543230.06480.267
    宁夏Ningxia52760.420.000.000.00527600.000.00000.000
    青海Qinghai713211.8217.360.000.0071322917.000.00000.00002
    山东Shandong37430.5669878.4737569.458298.611531771157470.05420.756
    河南Henan66458.3351927.0821232.6421597.22161215947570.13401.658
    江苏Jiangsu399.310.0026406.2570364.5897170967710.72410.996
    安徽Anhui260.426197.9239253.4788090.281338021335420.65840.998
    四川Sichuan148194.4558628.4764809.03183663.204552953071010.40330.675
    湖北Hubei2517.3617135.4219913.19135885.421754511729340.77450.986
    重庆Chongqing0.00468.75694.4476041.6777205772050.98491
    上海Shanghai0.00312.503020.832309.03584256420.40921
    浙江Zhejiang260.42937.50277.7892378.4793854935940.98430.997
    湖南Hunan0.000.006961.80186961.811939421939420.96411
    江西Jiangxi1215.282395.833645.83145503.471527601515450.95250.992
    云南Yunnan25086.8170538.20190364.5955729.173147193166320.16310.927
    贵州Guizhou0.000.00260.42159479.171597491597400.99841
    福建Fujian0.00972.223628.47104097.221086981086980.95771
    广西Guangxi0.000.006267.36202569.452088372088370.97001
    台湾Taiwan16111.116440.978819.44104.1731476153650.00330.488
    广东Guangdong17.3686.8135902.78118593.751546011545830.7671
    香港Hongkong0.000.00937.500.009389380.00001
    海南Hainan16701.3911423.61520.830.0028646119440.00000.417
    总(中国)Total (China)5886423.71393281.24486614.581665104.19843539825757670.19730.3054
    下载: 导出CSV

    表  4  当前和未来气候条件下松墨天牛的预测面积

    Tab.  4  Predicted areas for Monochamus alternatus under current and future climatic conditions

    时间Decade气候情景Climate scenarios非适生区

    Unsuitable area
    低适生区

    Lowly suitable area
    中适生区

    Moderately suitable area
    高适生区

    Highly suitable area
    预测面积Predicted area(km2)占总面积的比例Proportion of total area预测面积Predicted area(km2)占总面积的比例Proportion of total area预测面积Predicted area(km2)占总面积的比例Proportion of total area预测面积Predicted area(km2)占总面积的比例Proportion of total area
    1950-2000Current58864240.69564240280.05014866150.057516651040.1968
    2050sRCP2.652701220.62337778820.09204111280.048619956770.2360
    RCP4.553264930.63007914930.09363339760.039520028470.2369
    RCP8.550147570.593110819440.12804488890.053119092190.2258
    2070sRCP2.652619970.62247436630.08804828130.057119663370.2326
    RCP4.551998610.61508204170.09704550350.053819794970.2341
    RCP8.552396880.61977176390.08494272050.050520702780.2449
    下载: 导出CSV
  • [1] IPCC(Intergovernmental Panel On Climate Change). Climate Change 2007: Impacts, adaptations and vulnerability. fourth assessment report of working group II [R]. Cambridge, UK: University Press.
    [2] 丁一汇,任国玉,石广玉,等. 气候变化国家评估报告(Ⅰ):中国气候变化的历史和未来趋势[J]. 气候变化研究进展,2006,2(1):3−8.
    [3] Chen IC, Hill JK, Ohlemüller R, et al. Rapid rangeshifts of species associated with high levels of climatewarming[J]. Science, 2011, 333: 1024−1026. doi: 10.1126/science.1206432
    [4] Pio DV, Engler R, Linder HP, et al. Climate change effects on animal and plant phylogenetic diversity in southern Africa[J]. Global Change Biology, 2014, 20: 1538−1549. doi: 10.1111/gcb.12524
    [5] Dieleman CM, Branfireun BA, Mclaughlin JW, et al. Climate change drives a shift in peatland ecosystem plant community: implications for ecosystem function and stability[J]. Global Change Biology, 2015, 21(1): 388−395. doi: 10.1111/gcb.12643
    [6] 杨春平,赵霞,王嘉雯等. 基于MaxEnt模型的长足大竹象在中国潜在分布区及其对气候变化的响应[J]. 四川农业大学学报,2020,38(06):755−763.
    [7] Raza MM, Khan MA, Arshad M, et al. Impact of global warming on insects[J]. Archives of Phytopathology and Plant Protection, 2015, 48(1): 84−94. doi: 10.1080/03235408.2014.882132
    [8] Porter JH, Parry ML, Carter TR. The potential effects of climatic change on agricultural insect pests[J]. Agricultural and Forest Meteorology, 1991, 57(3): 221−240.
    [9] Emerson KJ, Merz CR, Catchen JM, et al. Resolving postglacial phylogeography using high-throughput sequencing[J]. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(37): 16196−16200.
    [10] Raes N, Roos MC, Slik JWF, et al. Botanical richness and endemicity patterns of Borneo derived from species distribution models[J]. Ecography, 2009, 32: 180−192. doi: 10.1111/j.1600-0587.2009.05800.x
    [11] Bateman BL, Vanderwal J, Williams SE, et al. Biotic interactions influence the projected distribution of a specialist mammal under climate change[J]. Diversity and Distributions, 2012, 18: 861−872. doi: 10.1111/j.1472-4642.2012.00922.x
    [12] Zhang MG, Zhou ZK, Chen WY, et al. Major declinesof woody plant species ranges under climate change in Yunnan, China[J]. Diversity and Distributions, 2013, 20: 405−415.
    [13] Anderson RP, Raza A. The effect of the extent of the study region on GIS models of species geographic distribution and estimates of niche evolution: preliminary tests with motane rodents(genus Nephelomys) in Venezuela[J]. Journal of Biogeography, 2010, 37(7): 1378−1393. doi: 10.1111/j.1365-2699.2010.02290.x
    [14] 朱耿平,刘国卿,卜文俊,等. 生态位模型的基本原理及其在生物多样性保护中的应用[J]. 生物多样性,2013,21(1):90−98.
    [15] 许仲林,彭焕华,彭守璋. 物种分布模型的发展及评价方法[J]. 生态学报,2015,35(2):557−567.
    [16] Yang XQ, S. P. S. Kushwaha, S. Saran, et al. Maxent modeling for predicting the potential distribution of medicinal plant, Justicia adhatoda L. in Lesser Himalayan foothills[J]. Ecological Engineering, 2013, 51: 83−87. doi: 10.1016/j.ecoleng.2012.12.004
    [17] 齐国君,高燕,黄德超,等. 基于MAXENT的稻水象甲在中国的入侵扩散动态及适生性分析[J]. 植物保护学报,2012,39(2):129−136.
    [18] 冯益明,刘洪霞. 基于Maxent与GIS的锈色棕榈象在中国潜在的适生性分析[J]. 华中农业大学学报,2010,29(5):552−556.
    [19] 王茹琳,李庆,封传红,等. 西藏飞蝗潜在分布对气候变化响应研究[J]. 东北农业大学学报,2017,48(08):60−71. doi: 10.3969/j.issn.1005-9369.2017.08.008
    [20] 朱耿平,王晓静,刘国卿,等. 悬铃木方翅网蝽在我国的潜在分布分析[J]. 应用昆虫学报,2012,49(6):1652−1658.
    [21] 朱诚棋,王博,沈婧等. 松墨天牛综合防治进展[J]. 中国植保导刊,2017,37(02):19−24. doi: 10.3969/j.issn.1672-6820.2017.02.004
    [22] 闫冠华,李巧萍,邢超. 不同温室气体排放情景下未来中国地面气温变化特征[J]. 南京信息工程大学学报,2011,03(1):36−46. doi: 10.3969/j.issn.1674-7070.2011.01.004
    [23] 徐汝梅,叶万辉. 生物入侵--理论与实践 [M]. 北京:科学出版社,2003:159−185.
    [24] 李淑华. 气候变化与害虫的生长繁殖、越冬和迁飞[J]. 华北农学报,1994,9(2):110−114. doi: 10.3321/j.issn:1000-7091.1994.02.022
    [25] 杨忠武,黄吉平,杨春生,等. 桂林市松褐天牛生物学特性研究[J]. 广西林业科学,2010,39(01):5−7+10. doi: 10.3969/j.issn.1006-1126.2010.01.002
    [26] 王玲萍,陈顺立,武福华,等. 松墨天牛幼虫空间格局的研究[J]. 福建林学院学报,2002,22(1):78−81.
    [27] 周书永,陈绘画,徐卫民,等. 松墨天牛成虫发生期预测[J]. 东北林业大学学报,2013,41(07):100−103.
    [28] 周樟庭,徐真旺,戴黎瑶,等. 冰雪冻灾后松墨天牛成虫种群数量动态研究[J]. 中国森林病虫,2011,30(03):18−19+24. doi: 10.3969/j.issn.1671-0886.2011.03.006
    [29] 徐正会,段艳,史胜利,等. 云南省松墨天牛地理分布及危害程度调查[J]. 西南林学院学报,2010,30(02):37−4.
    [30] Ha H, Heumann BW, Liesch M, Wang X. Modelling potential conservation easement locations using physical and socio-economic factors: A case-study from south-east Michigan[J]. Applied Geography, 2016, 75: 104−115. doi: 10.1016/j.apgeog.2016.08.009
    [31] 郭杰,刘小平,张琴,等. 基于Maxent模型的党参全球潜在分布区预测[J]. 应用生态学报,2017,28(3):992−1000.
    [32] 张天蛟,刘刚. 提高生态位模型时间转移能力的方法研究[J]. 中国农业大学学报,2017,22(2):98−105. doi: 10.11841/j.issn.1007-4333.2017.02.012
    [33] 王茹琳,李庆,封传红,等. 基于MaxEnt的西藏飞蝗在中国的适生区预测[J]. 生态学报,2017,37(24):8556−8566.
    [34] 沈阳,于晶,郭水良. 不同气候变化情境下中国木灵藓属和蓑藓属植物的潜在分布格局[J]. 生态学报,2015,35(19):6449−6459.
    [35] 张晓华,高云,祁悦,等. IPCC第五次评估报告第一工作组主要结论对《联合国气候变化框架公约》进程的影响分析[J]. 气候变化研究进展,2014,10(1):14−19.
    [36] Li RQ, Xu M, Wong MHG, et al. Climate change threatens giant panda protection in the 21st century[J]. Biological Conservation, 2015, 182: 93−101. doi: 10.1016/j.biocon.2014.11.037
    [37] Worthington TA, Zhang T, Logue DR, et al. Landscape and flow metrics affecting the distribution of a federally-threatened fish: Improving management, model fit, and model transferability[J]. Ecological Modelling, 2016, 342: 1−18. doi: 10.1016/j.ecolmodel.2016.09.016
    [38] 王茹琳,李庆,何仕松等. 中华猕猴桃在中国潜在分布及其对气候变化响应的研究[J]. 中国生态农业学报,2018,26(01):27−37.
    [39] Yue TX, Fan ZM, Sun XF, et al. Surface modelling of global terrestrial ecosystems under three climate change scenarios[J]. Ecological Modelling, 2011, 222(14): 2342−2361. doi: 10.1016/j.ecolmodel.2010.11.026
    [40] 赵晓冏,巩娟霄,赵莎莎,等. 样本量及其空间分布对物种分布模型的影响[J]. 兰州大学学报(自然科学版),2018,54(02):208−215.
    [41] 朱耿平,刘强,高玉葆. 提高生态位模型转移能力来模拟入侵物种的潜在分布[J]. 生物多样性,2014,22(2):223−230.
    [42] Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions[J]. Ecological Modelling, 2006, 190(3): 231−259.
    [43] 涂正彬,戴强,张文. 最大熵理论Maxent模型在林业中的应用[J]. 四川林勘设计,2016(02):10−14+25.
    [44] 孔维娜,王慧,李捷,等. 温湿度对松墨天牛越冬幼虫寿命的影响[J]. 山西农业大学学报(自然科学版),2006,26(3):294−295.
    [45] 泽桑梓,闫争亮,赵涛,等. 我国松墨天牛防治及引诱技术研究的现状及前景分析[J]. 西部林业科学,2010,39(3):93−96. doi: 10.3969/j.issn.1672-8246.2010.03.021
    [46] 张星耀. 森林病理学研究的生态数学方法 [M]. 北京:中国林业出版社,1999:76−98.
    [47] 张星耀,骆有庆. 中国森林重大生物灾害 [M]. 北京:中国林业出版社,2003.
    [48] 陈瑜,马春森. 气候变暖对昆虫影响研究进展[J]. 生态学报,2010,30(08):2159−2172.
    [49] 常晓娜,高慧璟,陈法君,等. 环境湿度和降雨对昆虫的影响[J]. 生态学杂志,2008,27(4):619−625.
    [50] Parmesan C. Ecological and evolutionary responses to recent climate change[J]. Annual Review of Ecology, Evolution, and Systematics, 2006, 37: 637−669. doi: 10.1146/annurev.ecolsys.37.091305.110100
    [51] 张殷波,刘彦岚,秦浩,等. 气候变化条件下山西翅果油树适宜分布区的空间迁移预测[J]. 应用生态学报,2019,30(02):496−502.
    [52] Porter JH, Parry ML, Carter TR. The potential effects of climatic change on agricultural insect pests[J]. Agricultural and Forest Meteorology, 1991, 57(3): 221−240.
    [53] 王庆. 极端气候对三种森林虫害发生趋势的影响 [D]. 中国林业科学研究院,2016.
  • [1] 张仕芳.  基于MaxEnt模拟樟科无根藤属物种的潜在适生区 . 四川林业科技, doi: 10.12172/202309250002
    [2] 李旭, 彭雪峰, 彭培好.  基于MaxEnt模型的四川省黄杉潜在适生区模拟及主导环境因子探究 . 四川林业科技, doi: 10.12172/202309190002
    [3] 黄雪梅, 黄怡, 贾泽旭, 陈德朝.  气候变化背景下药用植物牡丹的潜在适生区分析 . 四川林业科技, doi: 10.12172/202302050001
    [4] 雷蕾, 李雨, 卢林, 杨墅冰, 杨桦.  利用川硬皮肿腿蜂携带高毒力白僵菌微胶囊防治松墨天牛的研究 . 四川林业科技, doi: 10.12172/202205240001
    [5] 杨梅, 石艳, 石博文, 陈玲, 杨桦.  球孢白僵菌及酶增效剂对松墨天牛幼虫的毒力测定 . 四川林业科技, doi: 10.12172/202209230001
    [6] 吴明, 李锋, 王茹琳, 赵金鹏.  四川省野核桃生境适应性及适生区划研究 . 四川林业科技, doi: 10.12172/202104150002
    [7] 胡超, 于静.  基于MaxEnt和GIS预测四川省杉木良种在湖北省同一适宜引种生态区 . 四川林业科技, doi: 10.12172/202108130002
    [8] 陈德朝, 达郎周, 鄢武先, 邓东周, 贺丽, 吴世磊, 余凌帆, 杨靖宇, 王嘉智, 张利.  基于AHP分析法的川西北高寒沙地适生治沙灌木筛选评价研究 . 四川林业科技, doi: 10.12172/201905210001
    [9] 任金, 于正伦, 吴让伟, 李国耀, 陈小平, 李山琦, 杨兴伟, 付俊雪.  四川达州市松墨天牛生物学特性及防控措施最佳时机的研究 . 四川林业科技, doi: 10.12172/202108170002
    [10] 梁勤彪, 麻文建, 马甲强, 李波, 王建平, 贺军花.  广元市松褐天牛危害情况调查 . 四川林业科技, doi: 10.12172/201911180002
    [11] 余吉, 付明霞, 宋心强, 高飞, 杨彪, 李生强.  基于MaxEnt模型的四川大相岭保护区藏酋猴(Macaca thibetana)生境适宜性评价 . 四川林业科技, doi: 10.12172/202003170002
    [12] 吴桂康, 陈章铭, 杨桦, 杨伟.  云南松林松墨天牛发生规律及生物学特性 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2019.03.017
    [13] 陈国熙, 马正发.  采用花绒寄甲防治松褐天牛的初步研究 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2017.01.021
    [14] 乔麦菊, 唐卓, 施小刚, 程跃红, 胡强, 李文静, 张和民.  基于MaxEnt模型的卧龙国家级自然保护区雪豹(Panthera uncia)适宜栖息地预测 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2017.06.001
    [15] 陶园媛, 向茂榕, 刘昊泽, 王理顺, 何勇, 杨伟, 杨桦, 杨春平.  4株天牛致病菌对川硬皮肿腿蜂的毒力测定 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2017.04.007
    [16] 李灿, 刘贤安, 王娟, 彭培好, 邵怀勇.  基于MaxEnt的四川省红豆杉潜在分布区分析及适宜性评价 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2017.05.001
    [17] 康小龙, 李晓冬, 李华清, 杨春琳, 康钉荣, 徐勇, 钟金, 王小梅, 马文俊.  APF-1型诱捕剂防治“松褐天牛”试验研究 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2016.05.009
    [18] 王茹琳, 姜淦, 王闫利, 林姗, 沈沾红.  气候变暖情境下华山松大小蠹在中国的潜在分布区预测 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2015.01.018
    [19] 唐中芬, 华启尧, 杨佐忠, 王吉斌.  松诱木防治松墨天牛试验初探 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2014.03.027
    [20] 陈家德.  加拿大应对气候变化的政策机制及其林业碳计量模型——赴加拿大太平洋林业中心考察报告 . 四川林业科技, doi: 10.16779/j.cnki.1003-5508.2013.02.024
  • 加载中
  • 图(7) / 表(4)
    计量
    • 文章访问数:  95
    • HTML全文浏览量:  36
    • PDF下载量:  10
    • 被引次数: 0
    出版历程
    • 收稿日期:  2023-06-12
    • 网络出版日期:  2023-11-30

    基于MaxEnt的松墨天牛在中国的潜在分布及其对气候变化的响应

    doi: 10.12172/202306120001
      作者简介:

      刘勇(1971—),男,四川兴文人,工程师,学士,主要从事森林保护,E-mail:695189352@qq.com

      通讯作者: 杨桦(1982—),男,副教授,博士,主要从事林业害虫综合防治研究,E-mail: yanghua151017@163.com
    基金项目:  四川省科技厅国际合作项目(2021YFH0160)

    摘要: 由于世界气候变暖而导致的害虫空间分布格局的巨大变化,增加了害虫预防难度。通过MaxEnt模型,结合GIS方法,模拟松墨天牛在中国的潜在分布区,预测、对比、分析松墨天牛在未来3种气候情景(RCP2.6、RCP4.5和RCP8.5)下分布范围和空间布局,及其中心节点移动的规律。结果表明:所有AUC值均超过0.95,准确性也处于“极好”的状态。在当前气候环境下,松墨天牛在中国的适生区范围为90°~125°E,20°~41°N,由西向东呈连续分布,集中分布于我国东南部;在未来气候条件下,中适生区面积减少,低适生区、高适生区面积增加,并向高纬度、高海拔地区扩散,RCP8.5扩散范围最大。适生区中心点向西北、东北方向迁移,迁移趋势在高浓度排放情景下 (RCP 8.5) 响应最为剧烈。为松墨天牛的合理区划提供了有效手段, 为应对气候变暖提供了理论依据。

    English Abstract

    • IPCC 第四次气候评估报告指出,全球气候变化已成事实[1],在过去的100年中(1906年至2005年),中国的年平均地表气温明显上升[2],气候变暖对生态系统的结构、功能,群落组成及物种分布等都会产生影响[3-6]。昆虫是动物界种类最多,数量最大的变温动物,气候变化直接作用于昆虫的行为、发育、分布、生殖和生存[7],对农林业生产和人类健康存在重要影响[8-9]

      物种分布模型(species distribution models, SDM)是一种重要的生物地理学手段,已被广泛用于对可能存在的生物栖息地进行预测,并对其进行地理分布规律研究[10-12]。已有14种被广泛使用 [13-16],而 MaxEnt模型是最准确的一个,它具有操作简单,计算时间短,结果稳定,用量少等优点。近年来,利用Max Ent模型,已经成功预测稻水象甲Lissorhoptrus oryzophilus Kuschel [17]、锈色棕榈象Rhynchophorus ferrugineus Oliver[18]、西藏飞蝗Locusta migratoria tibetensis Chen[19]和悬铃木方翅网蝽Corythucha ciliate Say[20]等昆虫在国内适生区。

      松褐天牛Monochamus alternatus Hope,又称松墨天牛,是鞘翅目Coleoptera 天牛科 Cerambycidae的一种昆虫[21],主要为害马尾松Pinus massoniana Lamb.,其次为害云杉Picea asperata Mast.、冷杉Abies fabri Craib等,其成虫是松材线虫Bursaphelenchus xylophilus的主要携带者,体内所携带的扩散型4龄松材线虫在松墨天牛补充营养和取食过程中,通过伤口感染健康松树。被感染的健康松树在短期内因脱水死亡,从而引起松树枯死[21-23]。目前尚无可行的方法来预防和治理松材线虫病,因而通过控制松墨天牛这一传播媒介,从而控制松材线虫病的扩散显得尤为关键。随着气候变暖,昆虫发生期提前,其地理分布朝着更高纬度和更高海拔地区扩散[1]。目前,我国在害虫与气候变暖关系方面的研究主要集中在农业害虫上[24]。对松墨天牛的研究多集中于生物学特性、防治方法等方面[25-29],而对于松墨天牛在全国范围的预测预警研究较少。为了探索松墨天牛合理区划的有效手段,通过结合MaxEnt模型和GIS技术,模拟其在中国的潜在分布区。预测、对比、分析未来3种气候情景(RCP2.6、RCP4.5和RCP8.5)下松墨天牛的分布范围、空间格局和中心点位移的变化趋势,为未来应对气候变化提供基础数据支持。

      • 选取MaxEnt模型软件,预测松墨天牛在中国的适生区分布。登录普林斯顿大学Robert Schapire计算机科学研究中心主页,可免费下载,最新版本为3.3.3k版[30]。通过ArcGIS预测松墨天牛实际及潜在分布区。

      • 运用生态位模型模拟物种适生区时[31-32],首先通过访问四川省林业有害生物综合信息系统(sc.30120.org/fpgis/secure/Login.aspx?ReturnUrl=%2ffpgis%2f)获取在2012—2016年间四川境内松墨天牛的分布数据[6];通过查询物种分布数据库、参考已发表的文献资料,获取其余地区松墨天牛分布数据。

        查询的数据库包括:“国际农业与生物科学中心数据库(CABI, http://www.cabi.org/) ”、“全球物种多样性信息库(GBIF, http://www.gbif.org/)”、“中国西南地区动物资源数据库(http://www.swanimal.csdb.cn)”以及教学标本资源共享平台 (http://mnh.scu.edu.cn/)。利用Google Earth查询各分布点的经纬度,根据要求,剔除重复、模糊及相邻记录的分布点,得到全省680个分布点,其他地区371个,最后用 Excel对得到的经纬度数进行处理,以*. CSV的格式保存 [33]

      • 环境变量为1970-2000年间的67个环境变量,其中包括19个生物气候变量和月平均气候数据(表1)。环境变量数据从worldclim网站(http://www.worldclim.org/version2)中下载。IPCC第五次评估中公布了4条典型污染路径(RCP2.6, RCP4.5, RCP6.0, RCP8.5)[22]。考虑到 RCP4.5和 RCP6.0都是中等水平的CO2排放情景,故选择优先级更高RCP6.0。所以我们选取了最低(RCP2.6)、中等(RCP4.5)、最高(RCP8.5)共3个排放情景作为未来气候数据。未来时间范围包括2050s (2041—2060年) 和2070s (2061—2080年)[34],数据可以从国际热带农业中心 (CIAT) 网站下载[35]

        表 1  生物气候特征表

        Table 1.  List of bioclimatic characteristics

        生物气候特征bioclimatic characteristics 缩写abbreviation
        年平均气温 bio1
        平均气温日较差 bio2
        等温性 bio3
        气温季节性变动系数(标准差*100) bio4
        最热月份最高温度 bio5
        最冷月份最低温度 bio6
        气温年较差 bio7
        最湿季度平均温度 bio8
        最干季度平均温度 bio9
        最暖季度平均温度 bio10
        最冷季度平均温度 bio11
        年降水量 bio12
        最湿月份降水量 bio13
        最干月份降水量 bio14
        降水量季节性变化 bio15
        最干季度降水量 bio16
        最湿季度降水量 bio17
        最暖季度降水量 bio18
        最冷季度降水量 bio19
        最低温度 tmin
        最高温度 tmax
        平均温度 tavg
        降水量 prec
      • 选取比例尺为1:4 000 000的中国国界、省界和县界行政区划图作为分析的底图。数据可从国家基础地理信息中心(http://www.ngcc.cn/)获取。

      • 为避免模型过度拟合、提高模拟精度,分析预测物种分布通用的19个生物气候变量和月平均气候数据之间的相关性,去掉相关度高但对MaxEnt模型预测结果贡献不大的变量[6,36]。参考Worthington等方法[37]筛选潜在环境变量:首先利用刀切法 (Jackknife)测定环境变量对模型预测的贡献大小,去除贡献小的变量。其次,利用SPSS软件进行spearman相关分析进一步筛选关键的限制因子,并对相关系数|r|≥0.8的变量进行评估。对比初始模型中的贡献率,剔除贡献率小的变量,以提高模型模拟精度。通过这些分析方法,我们能够精确地确定那些高度相关但对模型预测结果贡献不大的变量,并将其排除在分析范围之外,从而得到更可靠和准确的结果。

      • 使用MaxEnt模型对其进行分析,将松墨天牛分布点数据和环境因子数据作为输入。为了建立预测模型,随机选取75%的松墨天牛分布点作为训练集(training data)。剩余25%松墨天牛分布点作为测试集(test data)验证模型。选取刀切法来测定各变量的权重、选择创建环境变量响应曲线,其余参数为模型的默认值。通过 ArcGIS工具箱中的格式化转换工具,将模式中的数据转换成 Raster格式,从而可以在 ArcGIS中进行显示。使用 MaxEnt软件进行计算,通过ROC (Receiver Operating Characteristic)曲线分析法,评价模拟结果的精度。在ROC曲线分析法中,数值的范围介于0到1之间,数值越接近1,表示物种存在的可能性越大。曲线下面积即AUC值反映预测精度,理论取值范围为0.5到1,值越接近1,表示精度越高。评价标准见(表2) [6,17-19]

        表 2  AUG评价标准表

        Table 2.  Evaluation criterion of AUC

        AUC值范围Range of AUC value评价标准Evaluation criterion
        0.5≤AUC<0.6失败
        0.6≤AUC<0.7不佳
        0.7≤AUC<0.8较好
        0.8≤AUC<0.9
        ≥0.9极好
      • 通过刀切法评估不同环境变量对松墨天牛潜在适生区预测的贡献。计算“仅该变量”、“除此变量”和“所有变量”模拟时的训练得分,“仅该变量”得分较高,说明该因子预测能力较高,对物种分布作用大;“除此变量”得分显著降低,说明其携带更多的独特信息,对物种分布较为重要[33,38]

      • 根据松墨天牛的实际情况,参考IPCC报告中有关可能性评估的划分方法,利用“重新分类”功能[38],设置不同等级的分布值,标准为:概率<0.05为非适生区;0.05≤概率<0.33为低适生区;0.33≤概率<0.66为中适生区;概率≥0.66为高适生区。

      • 参考Yue等人的计算方法[38-39],统计高适生区不同时段的面积变化情况及其质心的位移,计算公式如下:

        $$ \left\{ \begin{gathered} {{x}}({{t}}) = \sum\limits_{{{i}} = 1}^I {\frac{{{{si}}({{t}}) \cdot X{{i}}({{t}})}}{{S(t)}}} \\ {{y}}({{t}}) = \sum\limits_{{{i}} = 1}^I {\frac{{{{si}}({{t}}) \cdot Y{{i}}({{t}})}}{{S(t)}}} \end{gathered} \right. $$ (1)

        式中:t表示不同时段;I为高适生区的单位栅格数目;si(t)为时段t的单位栅格面积;s(t)为时段t中高适生区的总面积;(Xi(t),Yi(t))为时段t中高适生区单位栅格的质心坐标; (X(t),Y(t))为时段t中高适生区的质心坐标。

        $$ D = \sqrt {\mathop {\left( {x\left( {t + 1} \right) - x\left( t \right)} \right)}\nolimits^2 + \mathop {\left( {y\left( {t + 1} \right) - y\left( t \right)} \right)}\nolimits^2 } $$
        $$ \theta = {\mathrm{arctg}}\left( {\frac{{y\left( {t + {\text{1}}} \right) - y\left( t \right)}}{{x\left( {t + {\text{1}}} \right) - x\left( t \right)}}} \right) $$ (2)

        式中:D为时段tt+1的高适生区位移距离;θ 为时段tt+1的高适生区位移方向,0°<θ<90°表示位移方向为东北,90°<θ<180°表示位移方向为西北,180°<θ<270°表示位移方向为西南,270°<θ<360°表示位移方向为东南。

      • 选择19个生物气候变量和月平均气候数据作为初始环境变量,研究表明[40-41],变量之间的自相关等问题会在预测过程中引入冗余信息,影响预测结果。经过环境变量的筛选后,保留了9个主要环境变量,在此基础上重建松墨天牛在中国分布的最大熵模型。刀切法检验结果显示(见图1) ,对松墨天牛分布贡献前三的环境变量依次为:最暖季降雨量(bio18)、2月最高温(tmax2)、1月最高温(tmax1),根据变量贡献标准,bio18贡献率最大,最适最暖季降雨量为550 mm(见图2),松墨天牛分布随降雨量增多而增加,当超过最大降雨量时,分布减少。

        图  1  环境因子刀切法检验结果

        Figure 1.  Environmental factor test results by knife method

        图  2  最暖季降雨量(bio18)

        Figure 2.  Rainfall in the warmest season(bio 18)

      • 图3中基于主导环境变量构建的松墨天牛在中国分布模型的AUC值为0.954,为“极好”标准;图4中所有预测结果AUC值也达到“极好”标准。表明预测 可信度较高,可用于分析气候变化对松墨天牛在中国分布的影响。

        图  3  当前气候条件下的ROC曲线图

        Figure 3.  ROC curve for current climatic conditions

        图  4  未来不同气候情境下MaxEnt模拟结果

        Figure 4.  MaxEnt simulation results under different climate scenarios in the future

      • 由统计分析(见表3)可知,当前松墨天牛在中国的适生区范围为90°~125°E,20°~41°N,位于亚热带、暖温带之间,由西向东呈连续分布,集中于我国东南部,总面积2.58×106 km2,占国土面积的26.80%。其中,高适生区面积为1.67×106 km2,占适生区总面积的64.70%,其主要分布于长江以南的绝大多数省份;中适生区面积为4.87×105 km2,占适生区总面积的18.87%,其环绕于高适生区周边;低适生区面积为4.24×105 km2,占适生区总面积的16.43%,包括辽宁、山东、河北等省份。

        表 3  当前气候条件下松墨天牛的潜在适生区预测

        Table 3.  Prediction of potential habitat area of Monochamus alternatus under current climatic conditions

        省份名称
        Province
        非适生区
        面积
        unsuitable area(km2)<5%
        低适生区
        面积
        lowly suitable
        area (km2)
        5-33%
        中适生区
        面积
        moderately suitable area
        (km2)33-67%
        高适生区
        面积
        highly suitable
        area (km2)
        ≥67%
        总面积
        Total(km2)
        总适生
        区面积
        Total suitable area
        (km2)
        占全国高适生
        区百分比
        Percentage of highly
        suitable areas in China
        占全国总适生
        区百分比
        Percentage of
        total suitable
        areas in China
        黑龙江Heilongjiang543697.920.000.000.005436980.000.00000.000
        内蒙古Inner Mongolia1290607.660.000.000.0012906080.000.00000.000
        新疆Xinjiang1755208.360.000.000.0017552080.000.00000.000
        吉林Jilin210277.781857.64329.860.0021246521880.00000.010
        辽宁Liaoning121267.360.0030746.533541.670.00155556342880.00000.220
        甘肃Gansu408628.484861.111510.42243.0641524366150.000590.016
        河北Hebei158333.3437864.58277.780.00196476381420.00000.194
        北京Beijing9565.977673.610.000.001724076740.00000.445
        山西Shanxi155503.474236.110.000.0015974042360.00000.267
        天津Tianjin3315.978767.360.000.001208387670.00000.726
        陕西Shaanxi149392.3630659.7210468.7513194.44203715543230.06480.267
        宁夏Ningxia52760.420.000.000.00527600.000.00000.000
        青海Qinghai713211.8217.360.000.0071322917.000.00000.00002
        山东Shandong37430.5669878.4737569.458298.611531771157470.05420.756
        河南Henan66458.3351927.0821232.6421597.22161215947570.13401.658
        江苏Jiangsu399.310.0026406.2570364.5897170967710.72410.996
        安徽Anhui260.426197.9239253.4788090.281338021335420.65840.998
        四川Sichuan148194.4558628.4764809.03183663.204552953071010.40330.675
        湖北Hubei2517.3617135.4219913.19135885.421754511729340.77450.986
        重庆Chongqing0.00468.75694.4476041.6777205772050.98491
        上海Shanghai0.00312.503020.832309.03584256420.40921
        浙江Zhejiang260.42937.50277.7892378.4793854935940.98430.997
        湖南Hunan0.000.006961.80186961.811939421939420.96411
        江西Jiangxi1215.282395.833645.83145503.471527601515450.95250.992
        云南Yunnan25086.8170538.20190364.5955729.173147193166320.16310.927
        贵州Guizhou0.000.00260.42159479.171597491597400.99841
        福建Fujian0.00972.223628.47104097.221086981086980.95771
        广西Guangxi0.000.006267.36202569.452088372088370.97001
        台湾Taiwan16111.116440.978819.44104.1731476153650.00330.488
        广东Guangdong17.3686.8135902.78118593.751546011545830.7671
        香港Hongkong0.000.00937.500.009389380.00001
        海南Hainan16701.3911423.61520.830.0028646119440.00000.417
        总(中国)Total (China)5886423.71393281.24486614.581665104.19843539825757670.19730.3054

        选取2050s和2070s两个未来时段为, RCP2.6、RCP4.5和RCP8.5三种气候模式,上述情景下松墨天牛在研究区的高适生区仍主要集中在湖北、重庆、浙江、湖南、贵州等省,中适生区减少,逐渐发展为高适生区,其中包括四川、河南等省高适生区已覆盖大面积中适生区,西藏自治区东南部也开始逐渐出现高适生区。低适生区面积增加,向高纬度、高海拔地区扩增。在不同的时段和气候情景下,松墨天牛的分布范围和面积会有所差异。2050年RCP8.5总适生面积最大,低适生区面积增加最多,2070年RCP8.5高适生区面积最大(见表4)。

        表 4  当前和未来气候条件下松墨天牛的预测面积

        Table 4.  Predicted areas for Monochamus alternatus under current and future climatic conditions

        时间Decade气候情景Climate scenarios非适生区

        Unsuitable area
        低适生区

        Lowly suitable area
        中适生区

        Moderately suitable area
        高适生区

        Highly suitable area
        预测面积Predicted area(km2)占总面积的比例Proportion of total area预测面积Predicted area(km2)占总面积的比例Proportion of total area预测面积Predicted area(km2)占总面积的比例Proportion of total area预测面积Predicted area(km2)占总面积的比例Proportion of total area
        1950-2000Current58864240.69564240280.05014866150.057516651040.1968
        2050sRCP2.652701220.62337778820.09204111280.048619956770.2360
        RCP4.553264930.63007914930.09363339760.039520028470.2369
        RCP8.550147570.593110819440.12804488890.053119092190.2258
        2070sRCP2.652619970.62247436630.08804828130.057119663370.2326
        RCP4.551998610.61508204170.09704550350.053819794970.2341
        RCP8.552396880.61977176390.08494272050.050520702780.2449

        图  5  松墨天牛高适生区分布

        Figure 5.  Distribution of high suitability area of Monochamus alternatus

      • 根据研究结果显示,在当前及RCP2.6情景下,松墨天牛高适生区面积呈先增后减的趋势。其面积从当前的1.67×106 km2增加至2070s的1.97×106 km2;总适生面积持续增加,从2.58×106 km2增加到2070s的3.19×106 km2(见表4) 。当前至2070s,高适生区中心点由湖南新化(当前) 经湖南怀化(2050) 位移回湖南新化 (2070s) (见图6) ;总适生区中心点由湖南溆浦(当前) 位移至湖南沅陵(2050、2070s) (见图7)。

        图  6  气候变化情景(RCP2.6、RCP4.5和RCP8.5)下松墨天牛高适生区中心点位置及其转移

        Figure 6.  The location and transfer of mean center of highly suitable area of Monochamus alternatus under climate change scenarios

        图  7  气候变化情景(RCP2.6、RCP4.5和RCP8.5)下松墨天牛总适生区中心点位置及其转移

        Figure 7.  The location and transfer of mean center of totally suitable area of Monochamus alternatus under climate change scenarios

        研究结果显示,在当前及RCP4.5情景下,松墨天牛高适生区面积呈先增后减的趋势。面积从当前的1.67×106 km2增加至2070s的1.98×106 km2;总适生面积持续增加,从2.58×106 km2增加至2070s的3.25×106 km2(见表4) 。当前至2070s, 高适生区中心点由湖南新化(当前)位移至湖南溆浦(2050、2070s)(见图6) ;总适生区中心点由湖南溆浦(当前) 经湖南泸溪(2050s)位移至湖南沅陵(2070s)(见图7)。

        当前及RCP8.5情景下,其高适生区面积呈持续增加趋势,由当前的1.67×106 km2增加至2070s的2.07×106 km2;总适生面积呈现先增后减趋势,由2.58×106 km2增加至2070s的3.22×106 km2(见表4)。当前至2070s, 高适生区中心点由湖南新化(当前)经湖南涟源(2050s)位移至湖南安化(2070s) (见图6);总适生区中心点由湖南溆浦(当前,2050s) 位移至湖南沅陵(2070s)(见图7)。

      • MaxEnt软件基于最大熵原理,以物种分布变量和环境变量为基础,对约束条件下可能的最大熵分布进行统计分析[42]。AUC值表示了不同试验的准确度,并且不受发病率和诊断阈值的影响[43]。根据研究结果,无论是在当前的气候模式下还是未来的气候模式下,对松墨天牛适生区的预测AUC值均超过0.95,这意味着松墨天牛的实际分布与模型预测的地理分布之间具有较高的拟合度,可以利用这些预测结果来进行松墨天牛在中国地理分布与气候关系的研究。根据Maxent结果显示,当前松墨天牛在中国的适生区在90°~125°E ,20°~41°N的亚热带、暖温带之间,分布于华东、华南地区。刀切法检验结果显示,影响松墨天牛适生区分布的决定条件是最暖季降雨量(bio18),适生范围500~1000 mm,这与夏季( 尤其是6—7月) 为松墨天牛普遍的盛发期有关,当最暖季降雨量小于550 mm时,湿度可能最适宜,松墨天牛的羽化率大[44-45];当降雨量超过550 mm时,成虫出孔数减少[46-47]。气候变暖对昆虫的发育、繁殖、物候发生期、种群丰富度、寄主植物等核心生命活动产生直接影响[3-6,48]。气象因子中温度和降水是决定病虫害空间格局的主要因素,对松墨天牛的发生与分布有显著的影响,其变化范围与病虫害潜在的空间格局有十分密切的关系[49]。高适生区和总适生区在中国的中心点均分布于湖南省北部,且在未来气候变化情境下的迁移趋势一致。位移方向为东北、西北方向,即向高纬度、高海拔方向位移,与Maxent模型预测得到的适宜分布图直观观测的变化结果相一致,同时与其他相关研究得到的全球气候变暖会引起物种向高海拔和高纬度迁移的结论相一致[50-51]。迁移趋势均在高浓度排放情景下 (RCP 8.5) 响应最为剧烈,这是因为气候变暖导致昆虫地理分布向高纬度、高海拔地区扩散[52]。王庆运用低温冰冬雪指数K法研究发现,极端的低温冰冻雨雪气候对松墨天牛的发生有促进作用[53]

        气候是决定地球上物种分布的最主要因素,而物种分布格局的变化则是对气候变化最明确和直接的反映。在全球气候变化日益加剧的背景下,利用MaxEnt模型,结合GIS技术,预测物种适生分布区的空间变化和迁移趋势已成为虫害防治的关键问题,将为虫害综合治理提供重要的科学支持。

    参考文献 (53)

    目录

      /

      返回文章
      返回