基于InVEST模型的昆明市生境质量时空演变分析

曾欣怡宋钰红

(1.西南林业大学园林园艺学院，云南 昆明 650224；2.国家林业与草原局西南风景园林工程技术研究中心，云南 昆明 650224)

引言

随着城市化进程的加快，全球正在经历前所未有的生境质量退化，导致生物多样性丧失严重[1]，生物多样性保护面临着巨大挑战。生境质量指生态系统为生物的生存与发展提供适宜条件的能力，在一定程度上可以反映区域生物多样性的情况，是生态系统服务功能的基础，更是影响人类福祉的关键因素[2]。人类活动与土地利用变化是导致生境质量下降的重要原因[3]，是生境质量的一个重要风险因素，基于土地利用对生境质量变化进行研究，对区域生物多样性保护与可持续发展有重要价值。

目前，关于生境质量评估的研究主要集中在2点。建立指标体系。如，杨建强等[4]通过生境质量和生态响应2个方面构建指标来评价莱州湾的生境质量，此方法主要基于实地调研数据构建评价体系，结果准确但消耗人力物力较大、时间较长，不利于广泛应用于大尺度研究。构建模型模拟。如，刘汉仪等[5]以土地利用为基础运用InVEST模型分析了粤港澳大湾区生境质量时空演变特征；王志强等[6]通过HSI模型对扎龙国家级自然保护区丹顶鹤繁殖生境适宜性进行了定量评价；高浩翔等[7]采用最大熵模型分析了野生马麝夏季生境适宜性。其中，InVEST模型因其获取数据相对容易、结果可视化强，评估精度较高被国内外学者广泛应用于自然保护区[8]、流域[9]、省域[10]、市域[11]等尺度。

云南省生物多样性极其丰富[12]，昆明作为云南省省会，滇中城市群中心城市，其城市发展与生态环境保护面临较大冲突。因此，本研究以昆明市为研究区，使用InVEST模型中Habitat Quality模块，基于2000年、2010年、2020年3期的土地利用数据分析2000—2020年昆明市生境质量以及变化情况，以期为昆明市权衡区域生态环境保护与发展，打造高原湖滨生态宜居城市提供理论依据。

1 研究区概况

昆明市位于我国西南地区，云贵高原中部(E102°10′～103°40′，N24°23′～26°22′)，下辖14个区县，总面积21060.5km2。地势整体呈北高南低的格局，以湖盆岩溶高原地貌形态为主，市中心海拔约1891m。属于亚热带高原季风气候，年平均气温15℃，四季如春，有“春城”之称。生态本地优越，森林覆盖率达到52.62%，生物多样性总体水平上物种较丰富，生态系统类型较多，局部地区生物多样性高度丰富。2021年，昆明市常住人口850.20万人，地区生产总值7222.5亿元。

2 数据来源与研究方法

2.1 数据来源

本研究选取2000年、2010年和2020年共3期的昆明市土地利用数据，数据来源于30m全球地表覆盖数据GlobeLand30(http：//www.globallandcover.com/)，N48_20和N48_25图幅号，其中，将昆明市土地利用分为耕地、林地、草地、灌木地、湿地、水体、建设用地7类。该数据总体精度83%，是目前为止能公开获取到的分辨率较高、分类精度较好的全球土地覆被数据集[13]。

2.2 生境质量评估方法

本研究采用生态系统服务评估和权衡(InVEST)模型中生境质量(Habitat Quality)模块对昆明市2000年、2010年和2020年生境质量进行评估。生境质量模块综合分析土地利用和威胁生物多样性的因素来得到生境质量地图[14]，其数值大小代表生境质量水平，数值越高表示该区域生境质量越好，物种生存条件更适宜，生物多样性水平就越高[15]。该模块通过威胁因子对生境的影响计算生境类型的退化程度，通过生境适宜度与生境退化度来计算生境质量[11]。其中，生境退化度计算公式：

式中，Dxj代表土地利用类型j中栅格x受胁迫程度(生境退化度)；R代表威胁因子个数；r代表单个威胁因子；Yr代表威胁因子r所占栅格数；y代表威胁因子r所在的栅格；Wr代表威胁因子r的权重；ry代表威胁因子强度；irxy代表威胁因子对每个栅格产生的影响；Sjr代表土地利用类型对不同威胁因子的相对敏感程度。

生境质量的计算公式：

式中，Qxj代表土地利用类型j中栅格x的生境质量水平，范围为0～1；Hj代表土地利用类型j的生境适宜度；k为半饱和常数；z采用模型默认参数[16]。

结合前人研究与研究区实际情况[15，17]，本研究将土地利用中耕地与人造地表作为威胁因子。具体参数设置主要参考InVEST模型手册以及相关案例研究[14，16]，威胁因子最大影响距离、权重、衰减方式如表1所示，各用地类型的生境适宜度及其对威胁因子的相对敏感度如表2所示。

表1 威胁源最大影响距离及权重

表2 土地利用对各威胁源的敏感度

3 结果与分析

3.1 生境质量指数变化

生境质量指数范围为0～1，数值越大表示生境质量越高，采用分区统计方法得到昆明市生境质量指数，2000年、2010年、2020年的生境质量指数平均值分别为0.7、0.7、0.68，整体呈减少的变化趋势；2000—2010年基本不变(增幅1%)，2010—2020年降幅4%。如表3所示，2000—2010年，昆明市14个区县中半数的区县生境质量指数不变，5个区县生境质量指数增加，增幅均在1.5%左右，生境质量指数呈降低趋势的区县为官渡区和禄劝彝族苗族自治县，降幅分别为3.2%和1.2%；2010—2020年昆明市14个区县中12个区县的生境质量指数呈降低趋势，降幅各异，其中降幅最大的为官渡区，达到25%，呈贡区次之，降幅达20.3%；富民县和禄劝彝族苗族自治县生境质量指数不变；2000—2020年，14个区县中仅富民县的生境质量指数不变，其余13个区县均呈降低趋势，其中官渡区降幅最大，达到27.4%。

表3 2000—2020年昆明市各区县生境质量指数

3.2 生境质量等级的面积变化

结合研究区实际情况与相关案例，将生境质量结果分为低[0，0.4]、较低[0.4，0.6]、中[0.6，0.8]、较高[0.8，0.9]和高[0.9，1.0]5个等级，如表4所示。

表4 2000—2020年昆明市各等级生境质量面积

2000年、2010年和2020年昆明市生境质量处于中等级、较高等级和高等级区域的面积相对平稳；低等级区域面积呈增加趋势，面积占比分别为1.7%、2.1%和6.1%；较低等级生境质量面积呈降低趋势，面积占比分别为38.4%、36.7%和33.9%。从时间尺度上看，2000—2010年，低、中、高等级面积增加，增幅分别为22.4%、4.5%和1.1%；较低、较高等级面积减少，降幅分别为4.5%、11.8%。2010—2020年，低、较高等级面积增加，增幅分别为191.1%、7.3%；较低、中、高等级面积减少，降幅分别为7.4%、3.3%、1.6%。2000—2020年，低、中等级面积增加，增幅分别为256.5%和1%；较低、较高、高等级面积减少，降幅分别为11.5%、5.3%和0.5%。整体上看，昆明市低等级生境质量面积呈增加趋势，但速率不快，而且整体占比不高；较低等级生境质量面积呈减少趋势，但整体占比最高；中、较高和高等级生境质量面积大致维持稳定，整体占研究区面积约60%，表明研究区生态本底优越，这也得益于昆明市近年来生态保护策略。

3.3 生境质量等级空间分布特征

从生境质量等级空间分布来看，如图1a、1b、1c所示，2000年、2010年和2020年昆明市生境质量整体呈现“北高南低，西高东低”的格局。3期生境质量等级中，高等级区域主要分布于研究区西部，包括禄劝彝族苗族自治县中西部、西山区大部分、安宁市西部、晋宁区南部，其他区县零散分布，用地类型以林地为主；较高等级区域呈片式分布，集中分布于滇池、阳宗海、松华坝水库、上游水库、清水海、云龙水库等水域；中等级区域主要分布于东川区以及寻甸回族彝族自治县东部，其他区县零散分布，富民县北部区域中等级面积显著增多；较低等级区域主要分布于嵩明县中部、石林彝族自治县西北部，在寻甸回族彝族自治县和宜良县大面积分布，其他区县零散分布；低等级区域分布变化最为明显，随着社会发展与人类活动影响，向周围扩张，分布范围具体来讲，自五华区东南部向西北部扩张，自盘龙区西南部向东北部扩张，自西山区东北向西南扩张，自官渡区西部向南部扩张，并向东北延伸至长水机场方向、向西南延伸至呈贡区，自安宁市中部向西部扩张，晋宁区环滇池向南扩张，其他各区县分布较为分散，在各自低等级分布区域周围缓慢扩张，主要用地类型为建设用地。

注：a、b、c为生境质量空间分布；d为等级转移格局。图1 2000—2020年昆明市生境质量空间分布及等级转移格局

3.4 生境质量等级转移分析

整体上看，如图1d所示，研究区生境质量等级转移面积变化最大为较低等级转为低等级，其面积为734km2。研究区生境质量等级转为高等级的区域以中等级转高等级为主，面积为503.9km2。生境质量等级转为中等级区域以较低等级转中等级为主，面积为569km2。生境质量等级转为较低等级区域以中等级转较低等级为主，面积为381.6km2。研究区生境质量等级不变面积17660km2，占研究区总面积83.9%。生境质量等级降低面积2068.8km2，占研究区总面积9.8%，生境质量等级升高面积1331.7km2，占研究区总面积6.3%。

如表5所示，2000年高等级生境质量区域中有0.7%和2.5%的面积分别转化为低和较低等级生境质量，有8.5%和0.1%的面积分别转化为中和较高等级，88.2%的面积仍然保持高等级。由高等级转化为低和较低等级的区域零散分布于主城5区、其他区县也有少量分布，其主要原因是用地类型由林地转为建设用地。

表5 2000—2020年昆明市生境质量等级面积占比变化转移矩阵

较高等级生境质量区域中有2.1%和5.3%的面积分别转化为低和较低等级生境质量，有6%和3.2%的面积分别转化为中和高等级，83.3%的面积仍然保持较高等级。由较高等级转化为低和较低等级的区域主要零散分布于西山区东部、官渡区西部、嵩明县中部，其主要原因是用地类型由水体转为耕地和建设用地。

中等级生境质量区域有2.6%和6.4%的面积分别转化为低和较低等级生境质量，有0.3%和8.5%的面积分别转化为较高和高等级，82.1%的面积仍然保持中等级。由中等级转化为低和较低等级的区域主要分布于东川区东部、官渡区中部，其他区县少量分布，其主要原因是用地类型由草地转为建设用地和耕地。中等级转化为高等级区域整体分散，少部分集中分布于富民县昂大本山，其主要原因是用地类型由草地转为林地

较低等级生境质量区域中有9.1%的面积转化为低等级生境质量，有7%、0.3%和2.2%的面积分别转化为中、较高和高等级，81.3%的面积仍然保持较低等级。由较低等级转化为低等级的区域主要分布于滇池东北部、东南部和西南部、盘龙区和五华区城郊、官渡区东北部以及其他区县建成区域，其主要原因是这类区域用地类型由耕地转化为建设用地。较低等级转化为中等级区域主要分布于富民县东北部、寻甸回族彝族自治县中部，其他区县零散分布，其主要原因是用地类型由耕地转为草地。

低等级生境质量区域中有3.7%的面积转化为较低等级生境质量，有1%、0.1%和0.4%的面积分别转化为中、较高和高等级生境质量，94.8%的面积仍然保持低等级。由低等级转化为较低等级及以上的区域面积占比低、零散分布于各区县。

4 讨论

生境质量是生态文明建设的重要方面，对维持区域生物多样性起着关键作用。近年来，昆明经济建设的推进使得其社会经济高速发展，城镇化水平不断提高，建设用地、耕地增加明显，使得昆明市域生境质量显著下滑。根据统计，2000年昆明市常住人口481万人，GDP636.13亿元，到2020年昆明市常住人口达到846万人，GDP达到6733.79亿元，可见经济与人口增长带来的城市扩张是影响昆明市生境质量的重要因素，造成中心城区，即五华区、盘龙区、西山区、官渡区、呈贡区生境质量下降尤其明显。最近几年随着生态文明建设的推进，昆明部分区域开展的退耕还林、换草、还湖活动使得部分区域的生境质量有所提高。

昆明市作为滇中最大城市，其生态安全问题不容小视。在发展经济的同时也应注重生态环境的保持，协调人与自然的关系，统筹发展与保护之间的关系。因此在城市建设时应减少对林地、草地等自然环境的破坏，加强城市建设的集约化发展。并针对需要优先保护的区域划定生态保护区，确定好发展与保护的界限，实现绿色发展。

5 结论

研究基于土地利用数据，采用InVEST模型对昆明市20年来生境质量时空演变进行了分析，研究结果表明，昆明市整体生境质量处于中等级且呈下滑趋势，其中官渡区下滑最为严重；整体上昆明市各生境质量等级面积由高到低为较低、高、中、低、较低，处于中等级、较高等级和高等级区域的面积相对平稳，占研究区面积60%，低等级区域面积呈增加趋势，较低等级生境质量面积呈降低趋势；昆明市生境质量呈“北高南低，西高东低”的分布格局，最高点位于西部山区，最低点位于滇池东岸的城区以及各县的县政府所在地；20年来，昆明市生境质量等级以不变为主，生境质量等级转移面积变化最大为较低等级转为低等级。生境质量降低主要为城镇区域，是城市扩张的结果，生境质量增加主要为远离城镇的偏远山区。本文研究结果为协调昆明市城市发展与生态保护，以及生境保护视角下的城市规划提供了理论参考。