考虑地形与生物多样性的InVEST模型及其在县域生物多样性安全格局分析中的应用.pdf

考虑地形与生物多样性的 InVEST 模型及其在县域生物多样性安全格局分析中的应用陈朝1，杨贤房3,4，陈进栋2*，王琦1，刘婷1，邓南荣11.华南土壤污染控制与修复国家地方联合工程研究中心,广东省农业环境综合治理重点实验室,广东省科学院生态环境与土壤研究所2.福建师范大学地理科学学院3.江西师范大学地理与环境学院4.赣南师范大学地理与环境工程学院摘要县域尺度生物多样性安全格局精细评估可为区域生物多样性保护规划、实施与管理提供科学参考。以粤北山区翁源县为例，采用 InVEST 生境质量模型，结合区域地形与生物多样性特征，以生境质量指数、地形与生物多样性综合指数为评价指标，开展区域生物多样性安全格局综合评估。结果表明：翁源县整体生境质量指数（HQI）较高，尤其在人为扰动较弱的水系和林区最为明显，生境质量指数较低区集中在建设用地、耕地等人类活动影响强烈区域；地形多样性综合指数（TDCI）与地形起伏度和地形坡度位置指数分布格局较为一致，而与地形湿度指数相对差值分布格局相反；生物多样性指数（BCI）分布格局受到香农-威纳指数物种多样性本底及生物多样性丰富度、人类足迹指数、生物多样性重要度与植被总覆盖度的交互影响；生物多样性综合安全格局指数与 HQI、BCI 及 TDCI 分布有较好的一致性，低风险/低脆弱（高安全）区集中在县域内城镇建设用地、交通水利和耕地等人类活动强烈区，高风险/高脆弱（低安全）区主要分布在县域北部、西北部、东南部、东北部自然保护和饮用水源等人类活动干扰微弱区。通过与区域生态保护红线等相关资料进行验证与分析，改进优化后的评估结果更能精细和客观地反映研究区生物多样性安全格局。关键词InVEST 生境质量模型；生物多样性综合指数；地形多样性综合指数；生物多样性安全格局；翁源县中图分类号：X176 文章编号：1674-991X（2023）04-1345-09doi：10.12153/j.issn.1674-991X.20220783InVEST Model considering terrain and biodiversity and its application in theanalysis of county biodiversity security pattern:a case study of Wengyuan CountyCHEN Zhao1,YANG Xianfang3,4,CHEN Jindong2*,WANG Qi1,LIU Ting1,DENG Nanrong11.National-Regional Joint Engineering Research Center for Soil Pollution Control and Remediation in South China,Guangdong Key Laboratory of Integrated Agro-environmental Pollution Control and Management,Institute of Eco-environmental and Soil Sciences,Guangdong Academy of Sciences2.School of Geographic Sciences,Fujian Normal University3.School of Geography and Environment,Jiangxi Normal University4.School of Geography and Environmental Engineering,Gannan Normal UniversityAbstractDetailed assessment of the county-level biodiversity security pattern can provide scientific reference forthe planning,implementation and management of regional biodiversity conservation.Taking Wengyuan County inthe mountainous area of Northern Guangdong Province as an example,InVEST habitat quality model was used toachieve a comprehensive assessment of the regional biodiversity security pattern,by combining with regionaltopography and biodiversity characteristics and applying the habitat quality index and the composite topography andbiodiversity index as evaluation indicators.The results showed that the overall habitat quality index (HQI)of 收稿日期：2022-08-05基金项目：广东省科学院发展专项资金项目（2019GDASYL-0105042）；国家自然科学基金项目（41877514，42277479）作者简介：陈朝（1978），男，高级工程师，博士，主要从事土地利用、国土空间规划与生态安全相关研究，*通信作者：陈进栋（1986），男，工程师，博士，主要从事土地资源管理研究， Vol.13，No.4环境工程技术学报第 13 卷，第 4 期Jul.，2023Journal of Environmental Engineering Technology2023 年 7 月陈朝，杨贤房，陈进栋，等.考虑地形与生物多样性的 InVEST 模型及其在县域生物多样性安全格局分析中的应用 J.环境工程技术学报，2023，13（4）：1345-1353.CHEN Z,YANG X F,CHEN J D,et al.InVEST Model considering terrain and biodiversity and its application in the analysis of county biodiversity securitypattern:a case study of Wengyuan CountyJ.Journal of Environmental Engineering Technology，2023，13（4）：1345-1353.Wengyuan County was good,especially in the river system and forest areas with weak anthropogenic disturbance,and the areas with low HQI were concentrated in the those strongly influenced by human activities such asconstruction land and cultivated land.The distribution pattern of topographic diversity composite index(TDCI)wasconsistent with that of relief amplitude(RA)and slope position index(SPI),but opposite to that of the relativedifference of topographic wetness index(RDW).The distribution pattern of biodiversity composite index(BCI)wasaffected by the interaction between species diversity background in Shannon-Wiener index(H)and biodiversityrichness (BR),human footprint index (HF),biodiversity significance (BS)and total vegetation coverage.Thedistribution of biodiversity security pattern index(BSPI)was well consistent with that of HQI,BCI and TDCI.Low-risk/low-vulnerability areas were concentrated in the county areas with strong human activities,such as urbanconstruction land,transportation and water conservancy,and cultivated land.High-risk/high-vulnerability areas weremainly distributed in the areas with weak anthropogenic disturbances,such as nature protection and drinking watersources in the north,northwest,southeast and northeast of Wengyuan County.Through the verification and analysisof the relevant data such as the ecological protection redline,it was shown that the improved and optimizedevaluation results could reflect the biodiversity security pattern of the study area more precisely and objectively.Key wordsInVEST habitat quality model;biodiversity composite index;topographic diversity composite index;biodiversity security pattern;Wengyuan County 生物多样性是生物及其与环境形成的生态复合体以及与此相关的各种生态过程的总和，通常包含遗传多样性、物种多样性和生态系统多样性1-2，其中保护物种多样性和生态系统多样性是为了实现对生物多样性和自然栖息地的整体保护，属于区域或景观尺度生态安全格局及环保策略评估的重要目标3-4。自 1992 年联合国环境与发展大会（UNCED）正式签署生物多样性公约以来，众多生物多样性保护国际组织如世界自然保护联盟（IUCN）、国际鸟盟（Bird Life International）、联合国环境保护世界监测中心（UNEP-WCMC）、世界自然基金会（WWF）、全球森林观察（GFW）、大自然保护协会（TNC）等对全球和区域物种、生态系统层次的生物多样性安全格局越来越重视2,5-6。较多学者提出需要系统考虑物种丰富度以及物种多样性的重要度（范围稀有性）与人类活动扰动等因子效应，以综合反映物种多样性对区域生物多样性安全格局的影响7-8。地形是影响区域生物多样性安全格局的重要因素9。已有研究表明海拔高度10、坡位11、坡向12等地形因子通过影响区域水、热能量分配与物质迁移转换，进而影响到局域气候环境下的生物多样性分布格局，而中海拔、坡谷地区对应的中复杂程度地形被较多研究证明有较高物种丰富度13-14。常用的地形指数如地形位置指数（TPI）、地形湿度指数（TWI）、地形起伏度（RA）可以较好地反映地形坡位与分形类别、流域土壤水分饱和状态及地表起伏程度等方面的空间分异特征，因而被众多研究者采用以描述地形因子对区域生物多样性格局的影响15-16。但综合量化地形空间变化特征指标以反映对区域生物安全格局的研究尚且少见。生境多样性是生态系统多样性形成的基本条件，也是保护生物多样性的重要基础17。生境多样性一般用生境质量衡量，生境质量是生态系统为个体或种群的生存提供适宜的生产条件的能力，反映着生物栖息地质量状况，是生态系统和景观生物多样性的直接体现18，也被视为区域生物多样性和生态服务水平的重要表征和反映，同时也是区域生态安全保障和人类福祉提升的关键环节19。目前，国内外开展的生境质量研究主要集中在不同区域尺度下的土地利用重建20、河湖21、自然保护区22和特定物种23的生境质量变化及生态脆弱性评估24等方面，研究方法多采用生境适宜性评价模型，如机理和生态位类模型 HIS、MAXENT、GARP、ENFA等，经典统计学模型 LRM、GLM、MARS 等以及机器学习类模型 ANN、CART、RF 等25。在众多的生境质量量化评估模型中，由美国斯坦福大学、世界自然基金会和大自然保护协会共同开发的 InVEST生境质量模型26，充分考虑了土地利用变化对生态系统的影响，而且具有数据较易获取、模型参数调整灵活、评估准确度高、空间可视化能力强等优点27，因而在生态系统生境质量及服务评估中得到广泛应用。较多学者20,28-29已应用此模型评估了特定流域和区域生境质量时空格局变化。InVEST 生境质量模型侧重于从土地利用角度系统考虑生境质量优劣，但忽略了物种多样性本身及地形、水热等其他重要生境因子影响，而综合考虑生境多样性和物种多样性对于环境管理和基础研究具有重要价值3,17。鉴于此，笔者采用 InVEST 生境 1346 环境工程技术学报第 13 卷质量模型评估研究区生物多样性安全格局，再采用考虑地形与物种多样性方法改进优化上述格局，以期多视角综合反映研究区生物多样性安全格局，并为区域生物多样性保护与可持续利用、国土空间格局优化、生态安全及环境质量评估与治理等相关领域研究提供理论支撑与决策依据。1研究区概况翁源县（11330E11418E，2407N2440N），地处粤北山区“三山两谷”地形中的中、南部山系间的滃江盆地，流水和岩溶作用形成了典型的中低山、丘陵岗地与平原、阶地、水系纵横交错的复杂地形地貌（图 1）30，是北江上游重要的水源地、水源涵养区和南岭山地生物多样性维护区。翁源县属中亚热带季风气候，冬短夏长、雨热同期，光、热、水资源十分丰富。年均气温 20.3，年均降水量为1 778 mm，境内主要河流滃江为珠江水系北江的一级支流31。县域内重要野生动、植物保护资源包括云豹（Neofelis nebulosa）、豹（Leopard）、蟒 蛇（Pythonbivittatus）、黄 腹 角 雉（Tragopan caboti）、金 猫（Catopuma temminckii）、穿山甲（Manis）、水獭（Lutralutra）、大灵虎纹蛙（Hoplobatrachus）、水鹿（Cervusequinus）、斑 林 狸（Prionodon pardicolor）、鼬 獾（Melogale moschata）、赤麂（Muntiacus vaginalis）、野猪（Sus scrofa）、中华竹鼠（Rhizomys sinensis）、大白鹭（Ardea alba）；杉木（Cunninghamia lanceolata）、木荷（Schima superba）、罗浮柿（Diospyros morrisiana）、毛棉杜鹃花（Rhododendron moulmainense）、黧蒴锥（Castanopsis fissa）、栲（C.fargesii）、罗浮锥（C.faberi）、米槠（C.carlesii）、华润楠（Machilus chinensis）、鼠刺（Iteachinensis）、毛竹（Phyllostachys edulis）和甜槠（C.eyrei）等32-33。2研究方法 2.1数据来源及预处理采用的数据主要包括 IUCN 濒危物种红色名录物种多样性 2021 年空间分布矢量数据以及生物多样性丰富度（BR）1 km 栅格数据；UNEP-WCMC 提供的 2020 年人类足迹（HF）1 km 栅格数据；GFW 提供的 2018 年生物多样性重要度（范围稀有度，BS）1 km栅格数据；地理空间数据云（http:/ m 分辨率数字高程模型（DEM）数据；韶关市林业局提供的 2019 年森林资源二类调查成果中的翁源县植被总覆盖度（ZFGD）矢量数据；韶关市翁源县生态环境局提供的生态保护红线划定成果；韶关市翁源县自然资源局提供的翁源县 2019 年第三次国土调查统一时点更新数据。为便于计算县域尺度地形、生物多样性综合指数与生境质量指数，将不同尺度栅格大小统一重采样为 30 m 栅格。2.2地形多样性综合指数考虑到地形因素对生物多样性的复杂影响机理，参考相关研究文献及研究区地形特点，构建融合基于 TPI 的坡度位置指数（SPI）11、基于 TWI15的地形湿度指数相对差值（RDW），与 RA34的地形多样性综合指数（TDCI），以反映研究区生物多样性的地形因子影响程度。具体计算公式如下：TDCI=(SPI+RDW+RA)/3（1）TPI=z01nni=1zi（2）SPI=6TPISD50.5SDTPISD40.5SDTPI 530.5SDTPI 0.5SD,slope 52SDTPI 0.5SD1TPI SD（3）RDW=TWImaxTWImin;TWI=ln(SCAtan)（4）RA=HmaxHmin（5）TI=(TIiTImin)/(TImaxTImin)（6）式中：TDCI 为地形多样性综合指数；SPI、RDW、RA分别为 SPI、RDW 与 RA 归一化处理之后的地形指数；TIi为第 i 个地形指数栅格值；TImax和TImin分别为地形指数最大和最小值；SPI 为坡度位置指数，该指数综合了坡度和 TPI 所蕴含的地形特征信息，将地貌划分为 6 类11，根据研究区地形地貌特点并参考已有研究成果35-36，将山脊、上坡、中坡、平坡、下坡和山谷对应的生物多样性从低到高分别赋值 16；SD 为 TPI 计算结果的标准差；图 1 研究区位置及高程Fig.1 Location and DEM of the study area 第 4 期陈朝等：考虑地形与生物多样性的 InVEST 模型及其在县域生物多样性安全格局分析中的应用 1347 TPI 为地形位置指数，表示预定移动卷积窗口半径（MR）内的中心点（z0）高程与平均高程之差，m；n 为MR 半径内 DEM 栅格数量；zi为第 i 个栅格高程，m；RDW 为给定邻域卷积窗口内的地形湿度指数最大值 TWImax与最小值 TWImin的相对差值；TWI 为地形湿度指数，表示单位长度等高线上地表水所流经的上游区域汇水面积（SCA）与地表坡度正切值（tan）之比的自然对数；RA 为地形起伏度，m；Hmax和 Hmin分别为单位邻域面积内最大和最小高程，m。采用均值突变法37确定最佳邻域分析窗口，通过 Python 编程批量计算 N 个递增分析窗口（33，150150）下的平均 TPI、RDW 与 RA 最佳窗口，计算结果分别为 13、8、7 个 DEM 单元（其中 TPI 采用环形邻域，RDW 与 RA 采用矩形邻域），即 390 m 环形、240 和 210 m 矩形卷积窗口。2.3生物多样性综合指数构建生物多样性综合指数（BCI）以反映受人类活动扰动下的重要物种、群落层次下的生物多样性状态，计算公式如下：BCI=(BR+H)(1HF)ZFGDBS/2（7）H=si=1PilnPi（8）式中：BR 为 IUCN 生物多样性丰富度，根据每个 1km 栅格中可能出现的生态系统类型中物种数量求和计算；HF 为人类足迹数据，反映生态系统受人类活动扰动的压力状态，包括生物多样性、水、碳压力指数；ZFGD 为植被总覆盖度；BS 为生物多样性重要度（范围稀有性），根据 1 km 栅格像素面积除以每种物种的范围面积计算得出，即给定像素包含物种范围比例，将所有物种的计算值相加，反映每个像素对那里出现的物种的总体重要性；H为香农-威纳（Shannon-Wiener）指数，反映物种的多样性指数；s 为总的物种数；Pi为第 i 个物种个体数占总个体数的比例。当物种数为 1 时，H达最小值 0；当物种数为2 以上，且每个物种个体数量相等时，H达到最大值38。根据翁源县界提取 IUCN 所有红色名录两栖动物、哺乳动物、爬行动物、鱼类、淡水类群、鸟类和植物类别评估物种的空间范围对应的翁源县域内物种分布范围，在 ArcGIS 软件中分层提取每一物种分布范围计算 H指数并转成 30 m 栅格，最后叠加求和计算翁源县的 H指数。2.4InVEST 生境质量模型InVEST 模型生境质量模块对于生境质量的评价是通过分析生境斑块在所处基质中受到的各种威胁的综合影响来进行的27。衡量生境质量有 4 个要素：每种威胁的相对影响、每种生境类型对于每种威胁源的相对敏感性、自然栖息地（生境斑块）与威胁来源的距离以及土地受法律保护的程度。假设威胁源对生境的干扰强度是线性衰退，栅格 x 处的生境类型为 j，则生境类型在位置的退化程度如下所示26：Dxj=Rr=1Yry=1wr/Rr=1wrryirxyxSjr（9）irxy=1(dxy/drmax)（10）式中：Dxj为生境类型 j 中第 x 个生境像元的生境退化程度；r 为生境威胁源；wr为威胁源 r 的权重；y 为威胁源 r 中的栅格数；ry为土地利用图中每个栅格上威胁源的个数；Sjr为土地利用类型 j 对威胁源 r 的敏感性；x为生境保护程度；irxy为栅格 y 中的威胁源r 对栅格 x 的影响；dxy为栅格 x 与栅格 y 的距离，km；drmax为威胁源 r 的影响范围，km。生境质量指数的计算公式如下：Qxj=Hj1(DzxjDzxj+kz)（11）式中：Qxj为土地利用类型 j 中栅格 x 的生境质量；Hj为土地利用类型 j 的生境适宜度；k 为半饱和常数，一般取生境退化程度最大值的 1/2；z 为模型默认参数，可设置为 2.5。生境质量得分域为 01（生境适宜度最高时得分为 1，最低时得分为 0）。参考模型推荐值及已有研究成果26-28,39，威胁源及其相对权重和最大影响距离、生境适宜度及其对威胁源的相对敏感度赋值见表 1 和表 2。表 1 威胁源及其相对权重和最大影响距离Table 1 Threat sources and their relative weights andmaximum impact distances威胁源相对权重最大影响距离/km城镇用地0.82农村用地0.61耕地0.30.5铁路0.40.6主要公路0.40.7工矿用地12.5 2.5顾及地形与生物多样性的 InVEST 模型生物多样性安全格局基于生境质量指数（Qxj）、生物多样性综合指数（BCI）和地形多样性综合指数（TDCI）构建和计算区域水平的生物多样性安全格局指数（BSPI），计算公式如下：1348 环境工程技术学报第 13 卷BSPIx=Qxj1+BCIx2+TDCIx3（12）式中：1、2、3为采用层次分析法（AHP）计算的权重系数，综合考虑各指标对研究区生物多样性影响程度并参考前人研究成果28,40构建比较矩阵，计算比较矩阵得到的最大特征值（3.00 92）所对应归一化特征向量，获取 Qxj、BCIx和 TDCIx权重，CR 为0.008 8（0.1），通过一致性检验。按照自然断点法将生境质量、生物多样性及地形多样性综合指数及生物多样性安全格局指数评价结果从低至高分成 15 级。生物多样性安全格局按生物多样性安全格局指数高低分级，分别对应高安全、较高安全、中安全、较低安全、低安全 5 个等级。3结果与分析 3.1生境质量指数分布格局通过 InVEST 生境质量模型对研究区生境质量进行评价，得到 HQI 分布格局（图 2）。从评价结果来看，翁源县生境质量高值区主要分布在森林生态系统、草地生态系统和湿地生态系统，尤其是湿地生态系统中的河湖水系区是生境质量最高区，这些区域受人为扰动较小，生态系统相对稳定，为区域生物提供了较好的生态环境；低值区主要分布在人工生态系统和农田生态系统，这些区域人类活动强烈，生态系统受人为扰动较大，生境质量较低。总体来看翁源县大部分区域生境质量良好，尤其在人为扰动较弱的水系和林地区最为明显，生境质量较低区集中在建设用地、耕地和其他农用地等人类活动强烈区域。图 2 生境质量指数分布格局Fig.2 Habitat quality index distribution pattern 3.2地形与生物多样性指数分布格局SPI、RDW、RA 和 TDCI 分析结果如图 3 所示。总体来看，TDCI 与 RA 和 SPI 的分布格局类似，有较好的一致性，而与 RDW 分布格局相反。RA 较大区在中尺度效应加持下，坡上、山脊与山谷区生物生境更趋多样化，有利于不同生境物种的共 表 2 生境适宜度及其对威胁源的相对敏感度Table 2 Habitat suitability and its relative sensitivity to threat sources生态系统类型地类生境适宜度生境对各威胁源的敏感性城镇用地农村用地耕地铁路主要公路工矿用地农田生态系统耕地0.50.50.300.40.30.8森林生态系统林地10.80.650.50.60.70.9绿地0.80.70.50.350.550.50.8陆域自然保留地10.80.650.50.60.70.9园地0.70.60.450.350.60.650.8草地生态系统牧草地0.70.70.50.350.550.50.8湿地生态系统陆地水域0.80.750.60.40.550.50.8湿地10.750.550.40.60.70.9人工生态系统居住用地0.1000000广场用地0.1000000商业服务业设施用地0000000公共管理与公共服务设施用地0000000区域基础设施用地0000000仓储用地0000000道路与交通设施用地0000000工业用地0000000采矿盐田用地0000000公共设施用地0000000其他其他农用地0.20.2000.10.10.3特殊用地0.10.2000.10.10.3第 4 期陈朝等：考虑地形与生物多样性的 InVEST 模型及其在县域生物多样性安全格局分析中的应用 1349 存，物种丰富度更高；地形起伏度小、坡度位置平坦区，RDW 较大，尤其是河湖水系区，这些区域具有更大的坡面汇流面积或较低的水力坡降，土壤更容易达到饱和而产流，因此蕴藏有更丰富的生物多样性。BCI 高值区主要分布在县域西北部、东南与东北部，这些区域分别是翁北、西北森林自然公园、青云山自然保护区与贵东饮用水源-半溪自然保护区所在地，自然生境质量良好，物种多样性丰富，HF 低，但 ZFGD 与 BS 较高，因此，BCI 较高；低值区与TDCI、HQI 低值分布格局类似，主要为受人类活动扰动压力较大的建设用地与耕地区，BR 与 H本身较低，HF 高，ZFGD 与 BS 下降，导致 BCI 较低。3.3生物多样性综合安全格局基于生物多样性综合安全格局指数（BSPI）的分级结果见图 4 及表 3。从总体分布来看，翁源县BSPI 与 HQI、BCI 及 TDCI 分布有较好的一致性，且具有更精细的格局分布。低风险/低脆弱（高安全）区域面积为 37 873.71 hm2，占全域总面积的 17.44%，主要分布在县域内城镇建设用地、交通水利和耕地等人类活动干扰频繁而强烈区，物种相对贫乏，生物多样性相对较低，生态系统类型单一、脆弱，生境质量较差；中风险/中脆弱（中安全）区域面积 27 643.77hm2，占全域总面积的 12.73%，主要分布在滃江及支流河湖水系周边缓冲带、人工和农田生态系统周边与山麓结合带，为人类活动中度干扰区，物种丰富度一般，局部生物多样性较丰富，但生物多样性总体水平一般；高风险/高脆弱（较低安全及低安全）区域面积达到 151 638.30 hm2，占全域总面积的 69.83%，主要分布在翁北、西北森林公园、自然保护区与翁东南、翁东北饮用水源和自然保护区，为人类活动影响微弱区，物种和生态系统类型丰富，植被覆盖好，生境质量高。图 4 生物多样性综合安全格局及生态保护红线Fig.4 Comprehensive biodiversity security pattern andecological conservation redline 对比翁源县生物多样性综合安全格局与生态保护红线划定范围（图 4 及表 3），高风险低安全面积占生态保护红线划定面积的 71.32%，高风险低和较高 图 3 地形与生物多样性指数分布格局Fig.3 Topography and biodiversity index distribution patterns 1350 环境工程技术学报第 13 卷风险较低安全面积占比达到 96.69%，说明研究区生物多样性综合安全格局符合生态保护红线中南岭山地生物多样性维护-水源涵养生态保护区定位的实际情况，同时，这与原国家环境保护部 2015 年发布的中国生物多样性保护优先区域范围中指出，研究区所在的北江支流滃江流域大部分区域属于我国南岭生物多样性保护优先区域的客观事实相吻合。4讨论InVEST 模型生境质量模块基于土地利用类型的地类景观与生态威胁源的关系，通过分析生境质量（或者生物栖息地质量）间接反映生物多样性高低41，但忽视了不同土地利用强度与土地管理方式对区域生物多样性的影响，有研究表明农业有机耕作管理方式可能会增强局域物种多样性，而土地利用强度增加会减少局域物种丰富度42。因此，未来研究还需充分结合区域复杂的土地利用方式与环境生态交互影响过程，才能提升模型对区域生物生境质量的评估精度。限于县域尺度物种分布、丰富度等生物多样性相关数据的获取难度，生物多样性指数计算主要采用广泛使用的 IUCN 数据源，但有研究表明 IUCN 物种分布范围仍有较大的不确定性43，为此采用了生物多样性丰富度与重要度、人类足迹和植被总覆盖度指标修正生物多样性安全格局评价结果。此外，基于县域小尺度和数据源中能获取到的最详尽空间分辨率考量，统一将各种数据源重采样为 30 m 栅格以便于分析计算，尽管如此，更高精度的物种分布范围数据的获取用来精细刻画区域生物多样性格局仍然十分必要，野外样方调查、长期定位观测、生物多样性本底调查结合遥感观测、卫星追踪等手段获取生物地理分布数据可有效反映区域生物多样性空间分布特征。生物多样性和生境质量可分别代表物种多样性、生态系统多样性层次的区域生物多样性安全格局评价，在大中尺度景观水平上，生物多样性还包含着生态系统状况及弹性（或区域环境状况）和生态系统景观多样性等3，为此结合地形景观条件改进了区域生物多样性综合安全格局分布的评估方法，同时通过区域生态保护红线和中国生物多样性保护优先区域范围等相关资料进行验证与分析，表明改进优化后的评估结果更能较精细和客观地反映研究区生物多样性安全格局。限于篇幅仅探讨了地形景观对区域生物多样性安全格局的影响，而区域水热气候差异同样也会影响区域生物多样性富集与分布，且能更好地解释富集度的空间异质现象44，因此，在景观及更宏观的尺度上，有待进一步深入研究水热条件对区域生物多样性综合安全格局的影响。5结论（1）翁源县整体生境质量指数较高，尤其在人为扰动较弱的水系和林地最为明显，生境质量指数较低区集中在建设用地、耕地和其他农用地等人类活动影响强烈区域。（2）地形多样性综合指数与地形起伏度和地形坡度位置指数的分布格局具有较好的一致性，而与地形湿度指数相对差值分布格局相反。生物多样性指数分布格局主要与香农-威纳指数物种多样性本底分布及生物多样性丰富度、人类足迹指数、生物多样性重要度与植被总覆盖度的交互影响作用相关。（3）生物多样性综合安全格局指数与生境质量指数、生物多样性综合指数及地形多样性综合指数分布有较好的一致性，低风险/低脆弱（高安全）区主要分布在县域内城镇建设用地、交通水利和耕地等人类活动强烈区，高风险/高脆弱（低安全）区主要分布在翁北、西北、东南、东北自然保护和饮用水源等人类活动干扰微弱区。通过与区域生态保护红线等相关资料进行验证与分析，表明改进优化后的评估结果更能较精细和客观地反映研究区生物多样性安全格局。参考文献 蒋志刚,马克平.保护生物学的现状、挑战和对策J.生物多样性，2009，17（2）：107-116.JIANG Z G,MA K P.Status quo,challenges and strategy inConservation BiologyJ.Biodiversity Science，2009，17（2）：107-116.1 任海,郭兆晖.中国生物多样性保护的进展及展望J.生态科学，2021，40（3）：247-252.REN H,GUO Z H.Progress and prospect of biodiversity 2 表 3 生物多样性综合安全格局分布面积统计Table 3 Statistics on the distribution area of comprehensivebiodiversity security pattern安全格局面积/hm2面积占比/%生态保护红线内占比/%高安全 31 171.77 14.35 0.21较高安全 6 701.94 3.09 0.22中安全 27 643.77 12.73 2.89较低安全 70 503.12 32.47 25.37低安全 81 135.18 37.36 71.32总计217 155.78100100第 4 期陈朝等：考虑地形与生物多样性的 InVEST 模型及其在县域生物多样性安全格局分析中的应用 1351 conservation in ChinaJ.Ecological Science，2021，40（3）：247-252.钱云,翁丽珠,于长明,等.基于景观生态网络方法的市域生物多样性保护研究:以晋中市为例J.中国园林,2015,31(10):86-90.QIAN Y,WENG L Z,YU C M,et al.Study on city-widebiodiversity protection based on landscape ecological network:taking Jinzhong as an exampleJ.Chinese LandscapeArchitecture,2015,31(10):86-90.3 ZHAO X Q,XU X H.Research on landscape ecological securitypattern in a Eucalyptus introduced region based on biodiversityconservationJ.Russian Journal of Ecology，2015，46（1）：59-70.4 HUDSON L N,NEWBOLD T,CONTU S,et al.The database ofthe PREDICTS(Projecting Responses of Ecological Diversity InChanging Terrestrial Systems)projectJ.Ecology andEvolution，2017，7（1）：145-188.5 IUCN,Bird Life Interna