云南思茅土地利用变化方式对亚洲象栖息地的影响.pdf

No.3FORESTRESOURCESMANAGEMENT第3期June20232023年6 月林业资源管理云南思茅土地利用变化方式对亚洲象栖息地的影响金瑜，王娟，杨玥3,杨雄伟4,马惠敏4（1.思茅区六顺镇林业服务中心，云南思茅6 6 50 0 0；2.西南林业大学，昆明6 50 2 2 4；3.普洱市林业和草原局，云南思茅665000；4.思茅区林业和草原局，云南思茅6 6 50 0 0）摘要：2 0 1 7 一2 0 2 0 年，对思茅区亚洲象进行追踪监测，以第三次全国国土调查为基础数据，开展外业调查获得2020年土地利用类型数据，分析土地利用变化对亚洲象种群分布的关联性；对思茅区亚洲象栖息地适宜性进行预测和等级分析。结果表明：1）亚洲象喜欢在海拔6 0 0 1 1 50 m、坡度0 2 0、东北和东南方位、距离硝塘和河流较近、外界干扰较少的区域活动；适宜亚洲象的栖息地面积仅占思茅区总面积的1 4%，生境破碎化严重。2）2 0 0 9 一2 0 1 9 年，思茅区土地利用变化主要表现为林地、耕地和草地面积减小并转化为其他地类；种植园地、建设用地与水域面积增加，种植园地年均增长率达到3.9%，土地利用变化是导致亚洲象的适宜栖息地不断减少和破碎化的主要原因。3）2 0 1 7 一2 0 2 0 年，亚洲象在林地的活动占总监测次数的54.9%，这表明林地是亚洲象栖息活动的重要土地利用类型。关键词：亚洲象；土地利用变化；栖息地适宜性；最大滴模型中图分类号：S718.53文献标识码：A1文章编号：1 0 0 2-6 6 2 2（2 0 2 3)0 3-0 0 46-1 0D0I:10.13466/ki.lyzygl.2023.03.007Impacts of Land Use Changes on the Habitat ofAsian Elephants in Simao District of YunnanJIN Yu,WANG Juan,YANG Yue,YANG Xiongwei,MA Huimin4(1.Forestry Service Center of Liushun Town,Simao District,Puer,Yunnan 665000,China;2.Southwest Forestry University,Kunming 650224,China;3.Forestry and Grassland Bureau of Simao District,Puer,Yunnan 665000,China;4.Puer Forestry and Grassland Bureau,Puer,Yunnan665000,China)Abstract:Based on the tracking and monitoring of Asian Elephants in Simao District of Yunnan Provinceduring 2017 to 2020,and field survey for the land-utilization types in 2020,then combining with thebackground data of the third national-land survey in 2009 and 2019,this paper analyzed the correlationbetween land-utilization and Asian Elephant population distribution.The habitat suitability for AsianElephant in Simao District were predicted and evaluated.The main results were as follows:1)AsianElephant favoured to live in the regions where were 600 1 150 m of elevation,0 20of slope,northeastand southeast,close to the pond and river,less external interference.suitable habitat regions were onlyaccounted for 14%with severe habitat fragmentation.2)From 2009 to 2019,the characteristics of landuse type in Simao District were changed as follows:the areas for woodland,farmland and grassland were收稿日期：2 0 2 3-0 3-2 0；修回日期：2 0 2 3-0 5-1 7作者简介：金瑜（1 9 8 7），女，云南墨江人，工程师，硕士研究生，主要研究方向：森林资源保护和生物多样性保护与利用。Email:通讯作者：王娟（1 9 6 5），女，云南昆明人，教授，博士生导师，主要研究方向：生物多样性保护与利用。Email:47金瑜等土地利用变化方式对亚洲象栖息地的影响第3期decreased and transfered into other types.While the land areas for plantation,construction and waterstorage were increased,plantation with an annual growth rate by 3.9%.Land use change was the mainreason for decreased and fragmentation of the suitable habitat for Asian Elephant.3)The resaults of AsianElephant monitoring during 2017-2020 showed that the activity frequency of Asian Elephant in thewoodlands accounted for 54.9%.It indicated that woodland was an important place for the habitation ofAsian Elephant.Key words:Asian Elephant;land use change;habitat suitability;Maximum Entropy Model近年来，亚洲象（Elephasmaximus）肇事给社区群众的农业和林业生产甚至人身安全造成了极大的危害，引起社会各界的广泛关注。开展亚洲象种群防控与保护研究，对缓解人象冲突和实现人象和谐相处具有重要意义。亚洲象是我国I级重点保护野生动物 )，被自然资源保护联盟（IUCN）列为濒危物种 2 ；喜欢群居，单独活动较少，活动范围广 3;喜好在植被相对稀疏的竹阔混交林和灌木林地活动 4;阴?湿润、隐蔽性较好且距离水源较近（不超过1 0 0 m）的栖息地是亚洲象的最适宜栖息地 5。人象争地最终造成亚洲象栖息地逐渐缩小甚至丧失 6 ，导致人类和亚洲象产生了最大程度的空间交集 7 ，从而引发人象矛盾冲突 8 。翟佳豪等 9 对比分析了1 9 8 7 2 0 1 8 年西双版纳橡胶林沿道路、河流和边境线的变化，发现期间西双版纳橡胶种植面积在道路、河流沿线不断扩大。刘鹏等1 0 利用地理信息系统（Geographic Information System，GIS）与生态位因子分析（Ecological NicheFactorAnalysis,ENFA）模型分析得出，天然林是普洱市亚洲象适宜栖息的主要植被类型之一，但适宜栖息地范围较小且破碎化严重，连通度也较低，受人类活动干扰程度较大。魏建华 1 1 对南滚河流域近30 年土地利用变化对亚洲象栖息地的影响进行了研究，结果表明，南滚河流域土地利用变化使亚洲象丧失了部分生境，改善和提高适宜栖息地间的连通度是保护亚洲象和减少人象冲突的重要途径。Naha等 1 2 在印度东北部地区开展了人象冲突研究，从土地利用变化特征和环境因素等多方面对人象冲突的原因进行了剖析。目前，对亚洲象的研究多集中在道路等工程项目建设对亚洲象栖息地的影响及保护对策 1 3-1 4、亚洲象肇事原因分析及缓解对策探讨 1 5-1 6 、亚洲象种群数量结构分析以及栖息地适宜性预测17-18、亚洲象历史演变过程 1 9 、人象冲突发生原因的分析 2 0-2 1 等方面。云南省思茅区土地利用方式的变化不可避免会导致生境破碎化程度突出、森林面积减少，从而对亚洲象生存和繁衍活动产生巨大影响，致使亚洲象的最适宜生境急速退化、野外食物迅速减少。本研究通过分析思茅区的土地利用特征、亚洲象活动与土地利用变化的相关性，进行亚洲象栖息地适宜性预测，为合理转变土地利用方式和修复、改善亚洲象栖息生境提供重要的参考依据。1研究区概况1.1地理位置思茅区（2 2 2 7 2 30 6 N,1001910127E）位于云南省南部，地势西北高东南低，最高海拔2154.8m，最低海拔57 8 m，属南亚热带高原性季风气候，年平均气温1 9.4；年均降雨量1 340.9 mm。全区属于高原亚热带南部季风常绿阔叶林地带，森林植被种类丰富，针阔叶林地域分布明显。1.2中国亚洲象分布及其数量变迁亚洲象曾分布在我国从北到南的许多区域，其生境以0.5km/a的速度南移 2-2 3。2 0 世纪50 80年代，亚洲象在我国除藏南区域外，主要分布在云南南部地区一西双版纳州的景洪市和勐腊县、普洱市的西盟县和思茅区、临沧市的沧源县和德宏州的盈江县 2 4。2 0 世纪6 0 7 0 年代，由于当地的生态环境被严重破坏，思茅区和德宏的亚洲象逐渐消失 3。2 0 世纪8 0 年代，生活在普洱西盟、德宏盈江的亚洲象消失 2 5】,1 9 9 2 年起，野生亚洲象又重新在思茅市境内出现 3目前，野生亚洲象在我国境内仅分布于云南省48第3期林业资源管理普洱市思茅区、宁洱县、景谷县、澜沧县和江城县，西双版纳州景洪市、勐腊县和勐海县以及临沧市沧源县 2 6 。人类活动的开发与利用导致亚洲象分布核心区域的栖息地生境退化 2 7 我国亚洲象种群数量从2 0 世纪8 0 年代的1 40余头增加到2 0 2 2 年的30 0 多头 2 8 。其中临沧市沧源县南滚河亚洲象种群数量大约为1 2 头 2 9 ，分布的种群与普洱市和西双版纳州的种群处于地理隔离。西双版纳州勐腊亚洲象7 4头 2 6 ，勐海县与澜沧县1 9 头 30 ,西双版纳州勐养亚种群7 0 头 31 ，常在思茅区六顺镇、思茅港镇和倚象镇往返活动。普洱市境内亚洲象主要在思茅区、江城县、澜沧县活动，常年活动的象群约1 0 群1 2 7 头（数据来源：普洱市林业和草原局），其中，思茅区常年活动的有1 1 9 头，是普洱市亚洲象活动的重点区域。2研究方法2.1亚洲象数量及分布区调查2.1.1亚洲象野外监测和跟踪调查2017年开始对思茅区的亚洲象活动进行追踪观测，共聘请1 7 名亚洲象监测员。要求监测员采用直接跟踪的方法，每天至少选取耕地等亚洲象活动的土地利用类型之一进行记录，记录观测点和实际活动点的GPS定位点、象群数量等信息2.1.2访问调查2019年1 2 月至2 0 2 2 年6 月，在亚洲象经常出没的村庄，向经验丰富的居民和亚洲象监测员、护林员了解亚洲象活动规律及活动范围、取食习性等生活习性，以及村民对缓解人象冲突的观点与见解。2.1.3亚洲象野外监测数据处理用Excel汇总2 0 1 7 2 0 2 0 年思茅区亚洲象活动野外追踪监测数据，将其导人ArcGIS10.7软件，建立思茅区亚洲象分布点矢量数据，分析亚洲象的活动轨迹，确定亚洲象的活动范围，并结合收集的相关材料和数据，确定本研究的重点区域。由于在外业调查中采集的亚洲象分布点数据过于集中，分布点间的相关性较强，将会影响MaxEnt（最大熵）模型预测结果，因此选取1 0 0 m距离对数据进行去相关处理，最终确定348 6 个亚洲象活动分布点。将确定的亚洲象分布点属性表导出到Excel中，属性表中仅保留物种名称、经度、纬度3种信息，将其另存为.csv格式。2.1.4环境变量选择高程数据来源于美国太空总署（NASA）和国防部国家测绘局（NIMA）联合测量的数据（精度1 2.5m，http:/srtm.csi.cgiar.org/）。运用ArcGIS10.7软件计算坡度和坡向。结合研究区的实际情况和多次环境变量检验，最终选择海拔、坡度、坡向、距公路距离、距居民点距离、距河流距离、硝塘、林地面积、种植园地面积和耕地面积等1 0 个环境变量因子。利用ArcCIS10.7将计算整理好的环境变量因子转换成ASCI（美国信息交换标准代码）格式，用于MaxEnt模型运算。2.2最大熵模型预测方法2.2.1最大熵模型（MaxEnt）目前，很多学者采用生态位模型、机理模型和回归模型等3种模型进行栖息地适宜性分析 32】,其中，生态位模型中的MaxEnt模型采用物种分布点数据和环境变量数据相结合的方式来评估物种的栖息地适宜性，其优势在于即使在物种分布点数据较少时，也能得到比其他方法更为准确的预测结果 3。该方法已被广泛应用于物种栖息地适宜性、物种分布模拟等研究，如采用MaxEnt模型分析亚洲象生境适宜性及潜在分布区 342.2.2模型运行使用MaxEnt3.3.4建立模型，导人亚洲象分布点（.csv数据）和环境变量（ASCII）数据，随机选取亚洲象分布点数据的7 5%用作模型训练，其余的25%用于模型的验证，模型重复运行1 0 次，其他参数采用默认值，以1 0 次重复模拟的平均值为最终结果。2.2.3模型性能评价利用MaxEnt模型预测结果ROC曲线（ReceiverOperating Characteristic Curve，受试者操作特征曲线）进行模型性能分析，采用AUC(Area Under Curve，曲线下面积）为评判依据。AUC取值范围为0 1，一般认为，当0.5 AUC0.6时，预测结果失败；0.6AUC0.7时,预测效果较差;0.7 AUC0.8时,预测效果一般；0.8 0.9时,预测精度较高 3 49金瑜等土地利用变化方式对亚洲象栖息地的影响A第3期2.2.4栖息地适宜性等级划分通常依据模型敏感度和特异性数值相等的方法，选择MaxEnt模型的阈值，选择同等或者最大训练敏感性和特异性值作为模型阈值 36 ，然后根据阈值用ArcGIS10.7进行重分类，将思茅区亚洲象栖息地划分为4个等级：非生境 0,0.2 5）、边际生境0.25,0.50）、较适宜生境0.50,0.7 5）、最适宜生境 0.7 5,1.0 ，依据等级划分绘制思茅区亚洲象栖息地适宜性分布图。2.3思茅区土地利用类型分布特征调查2.3.1土地利用类型数据来源以2 0 0 9 年第二次全国国土调查和2 0 1 9 年第三次全国国土调查数据（以下简称国土三调）为基础数据，思茅区第三次全国国土调查全面采用0.5m分辨率的卫星遥感影像制作调查底图，结合现地调查举证获取成果数据。通过逐年对林地地类或林相变化遥感判读、森林资源档案核实、现地核查等手段，对2 0 1 7 一2 0 2 0 年林地“一张图”进行数据更新和完善信息，2.3.2土地利用类型野外外业调查2021年1 2 月至2 0 2 2 年6 月，以2 0 1 9 年国土三调为本底数据，采用样线调查法，在思茅区亚洲象重点活动区开展土地利用类型外业调查,共选择2 6 条亚洲象经常活动的路线，以样线能覆盖各种植被类型为前提，根据具体情况选取样线长度（一般为几公里至十几公里），分别在有亚洲象活动痕迹的地块、样线方向发生改变的地点、地类发生变化的地点，进行土地利用类型调查，并记录硝塘点及周边环境条件。外业调查采集了30 2 6 个点，并结合国土三调本底数据和林地“一张图”、森林督查等数据进行分析整理，同时还走访了熟悉当地情况的护林员和村民了解和证实2 0 1 7 年以来的土地利用现状及变化方式，从而获得2 0 2 0 年土地利用类型的相关数据。2.3.3土地利用类型分析方法参照土地利用现状分类 37 建立思茅区土地利用分类系统，利用ArcGIS10.7软件分析得到2009一2 0 1 9 年思茅区土地利用分布图和土地利用转移矩阵，并进一步分析思茅区土地利用类型的变化特征和分布格局；结合2 0 1 7 2 0 2 0 年土地利用类型和亚洲象野外监测点数据，分析思茅区土地利用变化与亚洲象活动的相关性。3研究结果3.1思茅区亚洲象分布状况1976年记录到7 头亚洲象在云仙乡、南屏镇一带活动 38 ；随后直到1 9 9 5年5头亚洲象才从西双版纳州重新迁入思茅区南屏镇。2 0 1 8 年是思茅区亚洲象数量增长最快的时期，从6 4头增加至9 6 头，各乡镇象群数量均增长，增长的主要原因除了相邻的江城县和西双版纳州象群的迁入外，还包括新象的诞生。2 0 1 9 年调查显示,思茅区常年活动的亚洲象约1 1 9 头，分4个象群和3头独象，其中，倚象镇44头、六顺镇45头、思茅港镇2 6 头、云仙乡4头（图1）。除了常驻象群外，还时常监测到西双版纳州象群间歇性的觅食迁入。2 0 2 2 年9 月1 6 日，六顺镇统计到1 0 4头亚洲象，亚洲象数量的增加导致六顺镇亚洲象肇事频率和范围增加。N思茅街道龙源乡思港锁简象镇洪市洪市10km洲象话动范围业洲象分板点图1思茅区亚洲象分布活动示意图Fig.1 Distribution of Asian Elephants in Simao district3.2MaxEnt模型栖息地适宜性分析3.2.1模型预测结果评估和环境变量重要性分析MaxEnt模型运行1 0 次后的AUC平均值为0.815，表明模型预测效果较好，可用于研究思茅区亚洲象栖息地适宜性。由表1 和图2 可知，硝塘对模型预测的贡献率最大，而海拔对模型预测的重要性最大，因此硝塘和海拔是影响亚洲象栖息地适宜性的重要环境因子。其次是耕地面积、种植园地面积、公路、林地面积和河流等因子，其余因子重要性相对较小。50第3期林业资源管理表1 各各环境变量对模型预测的相对贡献率和重要性排名Tab.1 Relative contribution rate and importance ranking ofeach environmental variable to the model prediction%环境变量相对贡献率重要性硝塘55.321.8海拔27.343.2耕地面积4.28.1种植园地面积3.69.0公路3.46.5林地面积2.74.3河流1.64.6居民点0.71.1坡度0.70.8坡向0.30.7公路坡向坡度居民点林地面积河流海拔硝塘种植园地面积耕地面积0.500.550.600.650.700.750.80AUC有变量无变量所有变量图2基于刀切法的环境因子重要性分析Fig.2 Importance analysis of environmental factor based onJackknifetest3.2.2亚洲象分布与栖息地环境变量关系MaxEnt模型分析结果显示，各环境变量对亚洲象活动的影响程度不同，对栖息地的适宜性随环境变量的改变而变化（图3）。亚洲象喜欢的生境特征为：海拔6 0 0 1 1 50 m；坡度在0 2 0 之间活动较频繁，在2 0 40 之间也经常活动，对坡度的适宜范围较广；对于坡向，主要偏向在东北和东南方位活动；硝塘、河流对亚洲象有积极的吸引作用；亚洲象喜欢在公路较少、外界干扰较少的区域活动，也喜欢在居民点附近、人为活动干扰较少的种植园地、耕地以及周边集中连片的林地进行取食活动。3.2.3栖息地适宜性分析本研究选择同等训练灵敏度和特异性来确定闵值，阈值为0.456 3，从模型预测思茅区亚洲象栖息地适宜性分布（图4）可知，亚洲象栖息的最适宜区主要位于思茅港镇的橄榄坝村和那澜村，六顺镇的官房村、团结村、南邦河村、炮掌山村，云仙乡的骂木村、挖令村，倚象镇的纳吉村和蚌弄村等地；较适宜区主要集中在六顺镇的那棵落村、南屏镇的曼歇坝村和南岛河村、云仙乡的坝塘村和黄竹林等地。10010095909080%/率谢用883R70%/率谢旺甲60504065306020551050000.00.20.40.70.81.0 1.21.404006008001000120014001600180020002200硝塘/（个/hm)海拔/m(a)硝塘(b)海拔9090807080%/率旺用%/率谢旺甲60705040503020401030020-100.01020304050 607080-1000100200300400500600700800耕地面积/hm种植园地面积/hm(c)耕地面积(d）种植园地面积51金瑜等土地利用变化方式对亚洲象栖息地的影响第3期901008590808075%/率谢雅甲%/率潮班甲288702885540503045204035103000.00.20.40.60.81.01.2-2024681012 141618202224林地面积/hmx106距公路距离/km(e)林地面积（f）距公路距离100959090%/率谢旺甲8885807075657060-10123456780123456距河流距离/km距居民点距离/km(g）距河流距离(h）距居民点距离9592909185809088757065876086-10010203040506070-50050100150200250300350400坡度/（）坡向的aspect数值/（）(i)坡度()坡向图3环境因子对预测概率贡献值的反应曲线Fig.3 Response curve of environmental factor to probabilityN非生境边际生境010较适生境km最适生境图4思茅区亚洲象栖息地适宜分布图Fig.4 Suitable habitat distribution map of AsianElephantinSimaoarea利用ArcCIS10.7统计得出，非生境面积2 2 0 4km占思茅区实际管辖总面积38 9 9 km的56.5%；边际生境面积1 1 52 km，占2 9.5%；较适宜生境面积370km，占9.5%；最适宜生境面积1 7 3km，占4.5%。由上述分析可知，思茅区适宜亚洲象栖息的生境面积较少，最适宜生境和较适宜生境仅占总面积的1 4%，生境破碎化严重，大面积的连续适宜生境较少。3思茅区亚洲象活动与土地利用变化相关性研究3.3.1思茅区土地利用变化特征分析1)思茅区土地利用类型分布。从2 0 0 9 年和52第3 期林业资源管理2019年土地利用现状图（图5）可知，2 0 0 9 年和2 0 1 9年思茅区土地利用类型以林地为主，且连片覆盖范围广泛，通常与种植园地相连；种植园地分布面积仅次于林地，靠近耕地，与林地相连；耕地与建设用地和种植园地紧密相连，主要分布在地势较平缓的山区、平坝或水域周围；其他土地类型分布相对较少。N林地草地10km耕地水域0种植园地其他土地建设用地(a)2009年N林地草地10耕地0水域km种植园地其他土地建设用地(b)2019年图52 0 0 9 年和2 0 1 9 年思茅区土地利用现状Fig.5 Land use status in 2009 and 2019 in Sima0 area2）思茅区土地利用变化特征。2 0 0 9 2 0 1 9 年思茅区林地、耕地和草地的面积减小，主要转化为其他地类，其中草地面积减少最显著；种植园地、建设用地与水域面积明显增大，其中种植园地扩张最显著；其他土地变化不明显。由于种植园地等土地扩张致使林地破碎化愈加严重，亚洲象栖息地遭到严重破坏，栖息地间的连通性也大幅度降低，适宜亚洲象栖息的生境不断缩减。从思茅区土地利用类型转移矩阵（表2）可知，20092019年，林地面积减少了9 6.6 6 km，并向其他地类转化。其中，林地向种植园地转化最明显，转化了2 1 1.0 1 km；向建设用地转化也较明显，转化了46.6 4km；其次是向水域和耕地转化，分别转化了44.1 4km和34.6 3km，这主要是因思茅区城乡扩建和水利设施等的建设占用了大面积林地。种植园地面积扩张了1 56.44km，主要种植咖啡、茶叶和橡胶，而种植园地向林地转化最明显，转化了7 3.0 1 km，主要是部分农户替换种植桉树以及对一些地块进行退耕还林。耕地面积减少了8 6.1 9 km，向林地转化最明显，转化了7 1.2 3km，向种植园地转化次之，主要是由退耕还林以及种植咖啡等经济树种所致。草地面积减少了1 0 0.2 5km，向林地转化了93.44km。湿地和其他土地变化较小但面积均有所增加。表2 思茅区土地利用类型转移矩阵Tab.2 Land use type transfer matrix of Simao districtkm2019年增长量林地种植园地耕地草地建设用地水域湿地其他土地总计林地2.609.68211.0134.637.9746.6444.144.851.092.960.01-96.66种植园地73.01270.3123.614.4916.686.900.310.84396.15156.44耕地71.2357.39115.270.3712.7111.591.000.87270.4386.19草地93.448.305.0218.601.886.190.520.09134.04100.252009年建设用地12.284.383.941.4653.971.220.020.5777.8457.74水域2.480.911.480.713.4722.870.320.1032.3464.59湿地0.440.020.220.110.043.870.2904.992.32其他土地0.790.270.070.080.190.150.091.642.01总计2.863.35552.59184.2433.79135.5896.937.313.653.877.4453金瑜等土地利用变化方式对亚洲象栖息地的影响第3期由表2、表3可知，林地、耕地和草地面积减少，并向其他地类转化了2 8 3.1 km，年均变化率为负值，分别为-0.3%，3.2%，7.5%，年均分别减少了9.6 7,8.6 2,1 0.0 2 km，其中草地年均减少的变化率最大。种植园地、建设用地和水域呈增长趋势，种植园地年均增长量最大，为1 5.6 4km，水域年均增长率最大，为1 9.9%。湿地和其他土地变化较小。表32 0 1 7 2 0 2 0 年思茅区亚洲象土地利用状况调查表Tab.3 Land use status survey of Asian Elephants in Simao district from 2017 to 2020林地耕地种植园地水域建设用地草地合计年份次数/次头次/头次数/次头次/头次数/次头次/头次数/次头次/头次数/次头次/头次数/次头次/头次数/次头次/头20177866.49541230783793326169237158151381 6451318720181416109295363 4344463865113592436201882.52118 88720191 3626 4965052.6534193.447491025796100239813 46420201.72511 1986453.4823633.366836252.82918763合计5289351182.09812647160713904311363141798414269.393643013.3.2思茅区亚洲象活动与土地利用的分布格局随着亚洲象种群数量的增加，林地成为亚洲象栖息活动的重点场所（图6），耕地和种植园地由于食源相对丰富且易取食也成为亚洲象取食活动频繁的场所。尤其是思茅港镇橄榄坝村和那澜村、南屏镇曼歇坝村和大开河、倚象镇的纳吉村和蚌弄村，由于耕地和种植园地分布范围广、地势较平坦加之水源丰富，而成为亚洲象经常光顾取食的场所。N林地水域耕地其他土地种植园地010亚洲象分布点km建设用地草地土地利用类型野外调查点图6思茅区亚洲象的土地利用现状图Fig.6 Land use status of Asian Elephants in Simao district由表3可知，在各类土地利用类型中，亚洲象对林地的利用最高，其次是对耕地和种植园地的利用。在林地累计监测次数和头次（监测次数和头次均指监测到亚洲象的次数和头次）分别占总监测次数9 39 3次的54.9%和总监测头次6 430 1 头次的54.6%，表明分布广泛，天然食源丰富、隐蔽性好的林地是亚洲象栖息活动的最佳选择地；耕地监测占总监测次数的21.6%和总监测头次的1 9.7%，亚洲象在农地的活动多数是取食农作物；种植园地监测占总监测次数的17.1%和总监测头次的2 1.6%，亚洲象在种植园地除了取食套种的玉米、芒果等农作物外，有时还休息玩要；建设用地监测占总监测次数的4.1%和总监测头次的3.1%，亚洲象经常进人村寨和家中觅食、活动；在其它类型的地类监测次数也相对较少。这与MaxEnt模型分析的适宜区结果一致。亚洲象分布数量和活动频率的变化与思茅区土地利用变化密切相关，种植园地面积的扩大和部分区域耕地面积的增加以及森林的采伐破坏了原本适宜亚洲象的栖息地（表2、表3），导致亚洲象在耕地和种植园地活动的频率和分布的数量增加，致使人象矛盾加剧。不合理的土地利用方式成为阻碍亚洲象种群数量增长的重要因素，同时也是引发人象矛盾的关键因素。4丝结论与建议4.1结论本研究通过分析2 0 1 7 2 0 2 0 年思茅区亚洲象追踪监测调查获取的分布点，结合国土三调、林地“一张图”等开展思茅区土地利用类型调查，最终选取1 0 个环境变量，并运用ArcCIS10.7和MaxEnt模54第3 期林业资源管理型分析思茅区土地利用变化对亚洲象活动的影响，得到如下结论：1)MaxEnt模型预测结果显示，硝塘和海拔是影响亚洲象栖息地适宜性的最重要因子，亚洲象喜欢在海拔6 0 0 1 1 50 m、坡度0 2 0 东北和东南方位、距离硝塘和河流较近、外界干扰较少的区域活动；适宜亚洲象的栖息地面积仅占1 4%，生境破碎化严重。22009年以来，思茅区土地利用变化特征主要表现为林地、耕地和草地面积减小并转化为其他地类；种植园地、建设用地与水域面积增加，种植园地年均增长率达到3.9%，林地利用的变化是导致亚洲象的适宜栖息地不断减少和破碎化的主要原因。3)20172020年的监测显示，亚洲象在林地活动占总监测次数的54.9%，表明林地是亚洲象栖息活动的重要土地利用类型，人类活动以及不合理的土地利用方式是造成适宜栖息地破坏的主要因素。4.2建议1）推进生态廊道建设，提高适宜栖息地间的连通度，减少栖息地“岛屿化”，促进亚洲象种间交流。2)鼓励转变土地利用方式，建立亚洲象森林公园。转变土地利用方式，丰富天然食源，修复和改善亚洲象栖息环境，加快推进普洱亚洲象公园的建设。3）建立生物防护带，减少人象空间交集区。在临近亚洲象经常出人地的村寨，在其耕地边缘种植一定宽度的亚洲象不喜欢且树干结实的树种，如茶叶、咖啡、花椒等，改变传统的种植模式，积极发展坚果、油茶等经济林产业，探索新的发展模式。参考文献：1 汪松，解炎.中国物种红色名录（第二卷）：脊椎动物M.北京：高等教育出版社，2 0 0 9.2 Smallidge P J,Leopold D J.Vegetation management for the mainte-nance and conservation of butterfly habitats in temperate human-dominated landscapes J.Landscape and Urban Planning,1997,38:259-280.3张立，王宁，王宇宁，等.云南思茅亚洲象对栖息地的选择与利用 J.兽类学报，2 0 0 3,2 3（3）：1 52-1 9 5.【4杨正斌，陈明勇，董永华，等.西双版纳国家级自然保护区勐养子保护区亚洲象生境现状分析 J.林业调查规划，2 0 0 6（3）：49-52.5李芝喜.亚洲象生境评价J.环境遥感，1 9 9 6，1 1（2）：108-115.6代云川.气候变暖背景下橡胶林适宜区的空间扩张及其对亚洲象栖息地的影响 D.昆明：云南师范大学,2 0 1 7.7靳莉.中国亚洲象肇事原因和对策研究 J.野生动物，2 0 0 8,2 9(4):220-223.8 Lee P C,Graham M D.African elephants and human-elephantinteractions:Implications for conservation J.International ZooYearbook,2006,40(1):9-19.9翟佳豪，刘影，肖池伟.1 9 8 7 2 0 1 8 年西双版纳橡胶林时空变化及其线状特征 J.热带地理，2 0 2 2，42（8）：1 37 6-1 38 5.【1 0 刘鹏，代娟，曹大藩，等.普洱市亚洲象栖息地适宜度评价 J.生态学报，2 0 1 6,36（1 3）：41 6 3-41 7 0.11】魏建华.南滚河流域土地利用变化对亚洲象生境的影响D.昆明：云南大学，2 0 1 9.12 Naha D,Sathyakumar S,Dash S,et al.Assessment and predictionof spatial patterms of human-elephant conflicts in changing landcover scenarios of a human-dominated landscape in North BengalJ.PloS one,2019,14(2):e0210580.13黄元，谭德安，白梁芳，等.思茅至江城拟建高速公路沿线亚洲象栖息地适宜性分析 J.林业建设，2 0 1 9,2 1 0（6）：6 7-7 3.【1 4】王冀，王云，关磊，等.热带雨林公路建设对野生动物的影响及保护研究进展 J.生态学杂志，2 0 1 9,38（1 0）：31 8 3-31 8 8.15肖洒.普洱地区亚洲象肇事特点及损失情况分析D.北京：北京林业大学，2 0 1 5.【1 6 谢文芳，宋军平，苏海萍，等，西双版纳国家级自然保护区野生动物肇事补偿现状及缓解对策 J.林业调查规划，2 0 2 0，45(2):183 186.17吕婷，郑璇，朱高凡，等.西双版纳勐腊亚洲象种群数量动态分析 J.林业建设，2 0 1 9,2 1 0（6）：9 1-9 6.18杨子诚.中国亚洲象分布及局域种群数量评估D.北京,北京林业大学，2 0 1 7.19】何成馨.近50 年西双版纳人象关系演变研究 D.昆明：云南大学,2 0 1 3.【2 0 彭勇，吴顺福.人象冲突原因分析及缓解措施J.山东林业科技,2 0 1 6,2 2 6(5):9 4-9 6.21】唐勤.西双版纳人象冲突与缓解对策D.昆明：昆明理工大学,2 0 0 7.22】孙刚，许青，金昆，等.野象在中国的历史性消退及与人口压力关系的初步研究 J.东北林业大学学报，1 9 9 8，2 6（4）：47-50.【2 3刘雪华.3S技术与中国野生动物生境评价M.北京：中国林业出版社,2 0 1 1.24张洁.历史时期中国境内亚洲象相关问题研究【D.西安：陕西师范大学,2 0 0 8.25张立.中国亚洲象保护研究M.北京：科学出版社，2 0 1 8.55