基于GIS与Fragstats的安徽省生态安全格局研究

白星亮　任冲　李景宜　徐云环　马晓武　马力

摘要 从区域生态安全重要性和敏感性角度出发，结合景观生态指数，构建安徽省生态安全格局评估指标体系，绘制安徽省生态安全格局图。结果表明：安徽省生态安全格局总面积为41 083.46 km2，占全省面积的29.53%，主体部分以大别山生态屏障和黄山生态屏障为支撑，总面积为30 883.63 km2，长江干流水生态廊道为屏障链接，巢湖等其他湖泊湿地为生态格局重要组成成分，总面积为1 130.74 km2。

关键词 生态安全格局;生态重要性;生态敏感性;景观生态指数

中图分类号 X826文献标识码 A

文章编号 0517-6611（2020）13-0070-08

Abstract Combined with the landscape ecological index， the assessment index system of ecological security pattern in Anhui Province was constructed from the perspective of ecological safety importance and sensitivity. On the basis of the above research， the ecological security pattern map of Anhui Province was drawn. The results show that the total area of ecological security pattern in Anhui Province is 41 083.46 km2， accounting for 29.53% of the provinces total area. The main part is supported by the Dabie Mountain Ecological Barrier and the Huangshan Ecological Barrier， with a total area of 30 883.63 km2. The barrier link is the Yangtze River mainstream water ecological corridor， and the important components of the ecological pattern are others lakes and wetlands such as Chaohu， with the total area is 1 130.74 km2.

Key words Ecological security pattern;Ecological importance;Ecological sensitivity;Landscape ecological index

生态安全是人类赖以生存的环境安全状况。近年来，人类对自然生态环境施加了前所未有的影响，环境污染、土地退化、森林减少、生物多样性丧失等一系列严重的生态环境问题严重威胁了国家和区域生态安全。区域生态安全格局研究已成为生态安全领域前沿研究方向之一[1]。

近年来，国内外相关研究者在区域生态安全领域开展了大量系统及精细的研究。研究重点主要集中在概念及理论探讨、指标体系设计、综合模型构建、评价方法研究等领域。其中，我国学者基于生态安全重要性和敏感性理论构建生态安全格局已成为当前研究领域的最重要方法[2-7]。与此同时，随着景观生態学的发展，景观生态与生态保护研究也在逐步展开[8-9]。有学者发现城市面临的多种问题都与景观的组成分布格局有关[8]。城市化改变景观格局和生态功能导致生态系统服务（ESS）下降[9]。移动窗口法等方法用于景观分布格局的研究[10-12]。景观破碎化程度影响生物多样性程度[13]，连通性指数对生态系统服务功能的退化程度也有重要影响[14]。

为构建完整生态预警机制，保障生态环境可持续发展，安徽省依据生态环境部颁布的《生态保护红线划定技术指南》发布了全省生态保护红线[15-16]。然而，由于安徽地理位置特殊性和经济地位重要性，生态保护红线的划定与全国其他地区略有差异[17-18]。若以全国统一标准构建安徽生态红线，则会出现生态保护与生态安全格局的协调和差异问题[19-21]，生态保护红线方案的实际效益和补偿价值将有所折损[22-24]。另外，已有研究表明安徽省生态安全现状与格局分布有发生改变的迹象[2]，具体表现生态系统服务功能重要性和敏感性出现南北分布差异，生态廊道集中于长江沿岸，生态源地向皖南地区集中[25]，以西部大别山和南部黄山生态屏障为重要保障的安徽省生态安全格局逐步显现。

综上，已有研究内容多数倾向于生态安全重要性与敏感性的单独探讨，尚未能从景观生态学角度对生态安全进行全面综合的考量，使得生态安全格局与景观生态分布之间的联系略有不足。该研究将从生态安全重要性、生态安全敏感性角度出发，依托景观生态指数，构建安徽省生态安全格局，进一步丰富省域生态安全评价实证研究和分析，为区域生态文明建设和政府制定相关政策提供技术支撑和决策依据。

1 数据来源与方法

1.1 研究区概况

安徽省位于114°54′～119°37′E，29°41′～34°38′N，地处长江下游，跨长江、淮河和钱塘江三大流域，面积1.401×105 km2，约占全国总面积1.45%（图1）。秦岭-淮河地理分界线横贯全省，全省自北向南可分为淮北平原、江淮丘陵、大别山区、沿江平原和皖南山区五大自然区域。气候属亚热带向暖温带过渡型，年均气温14～17 ℃，平均日照1 800～2 500 h，无霜期200～250 d，平均降水量800～1 800 mm。全省森林覆盖率28.65%。下辖16个省辖市，8个县级市，53个县，44个市辖区，是长三角的重要组成部分，处于全国经济发展的战略要冲和国内几大经济板块的对接地带。《2016年生态文明建设年度评价结果公报》显示，安徽省绿色发展指数不足全国平均水平，为79.02，居于19位;全省人均耕地面积为全国平均水平的97.1%，人均湿地面积为全国均值的43.4%;全省轻度污染湖泊水库有29座，轻度污染河流有101条;空气优良天数比例为74.3%。

1.2 数据来源

基础数据主要来源于《安徽统计年鉴（2010—2016）》，高程数据来源于地理空间数据云网站GDEMV2 30 m分辨率数字高程模型，蒸散发数据来源于NASA网站（https：∥ladsweb.modaps.eosdis.nasa.gov），NPP数据来源于全国生态状况遥感调查与评估成果NPP数据集，气温和降水资料来源于中国气象科学数据共享服务网，土壤数据来源于全国生态环境调查数据库，土地利用数据来源于清华大学数据网站分辨率为10 m土地利用图（http：∥data.ess.tsinghua.edu.cn/）。

1.3 研究方法

以构建安徽省生态安全格局为目标，综合定性判定和定量评估2种分析方法，开展生态保护功能评估，以生态系统服务功能重要性、生态环境敏感性为评估指标，结合景观生态指数，合理划定识别生态安全格局重要区域。

1.3.1 区域生态安全重要性。

区域生态安全重要性评估使用NPP定量指标评估方法[15]，结合安徽省地域气候特征，不考虑防风固沙重要性评估，选择水源涵养、水土保持和生物多样性保护3项指标为生态功能重要性评估对象，分别建立评估模型，针对安徽省地域生态现状，选择不同的等级分级标准。各评估对象的模型公式如下。

式（7）含义为景观中各斑块类型所占景观面积乘以各斑块类型之间相邻的格网单元数目占总相邻的格网单元数目的比例，乘以该值的自然对数之后的各斑块类型之和，除以2倍的斑块类型总数的自然对数，其值加1后再转化为百分比的形式。其中，Pi为斑块类型占用的景观比例i;gik为基于重复计数方法的斑块类型i和k的像素之间的邻接或连接数;m为景观中存在的斑块类型数量。CONTAG值越小，表明景观中存在越多小斑块，越接近100表明景观中有连通度高的优势斑块类型存在。

2 结果与分析

2.1 生态安全重要性

运用ArcGIS软件综合计算安徽省水源涵养服务功能、水土保持服务功能和生物多样性服务功能空间分布（图2）。

安徽省水源涵养能力高值出现在六安市南部和安庆市北部的大别山区，池州市东南部和黄山市大部的黄山山脉群，以及中部北部主要江河湖地带，最大值为280 826，低值出现在北部平原区域，最低值为6 821，极差为274 005。全省平均值为45 136，整体分布态势为南多北少，西多东少，山地水源涵养能力高于平原地区，个别低值区集中出现在长江两岸滩涂地带。

安徽省水土保持能力总体良好，分布状况和水源涵养功能分布类似，高值区出现在西部大别山区和东南部黄山区，主要包括六安市、安庆市、池州市和黄山市，最高值为759，平均值为391。北部平原地带水土保持能力较好，中东部滁州市山区丘陵地带有较完整的水土保持能力集中区，总体呈南多北少，东南多西北少分布。

安徽省生物多样性分布集中于山区丘陵地带，以南部黄山市黄山区为代表，以及宣城市黄山余脉地带、安庆市中部平原山地交汇地带等，呈现南多北少的分布状况。

2.1.1 水源涵养功能重要性分析。

式（1）计算出水源涵养功能空间分布后，在ArcGIS栅格计算器（Spatial Analyst → Raster Calculator）中将安徽省水源涵养功能空间分布栅格图采用最大最小值法将数据归一化到0～1，导出栅格数据属性表，属性表记录了每一个栅格像元的水源涵养服务值，将该组服务值求和，计算累加服务值[15]。将水源涵养服务值按从高到低的顺序排列，经详细试驗后，发现将累加服务值占水源涵养服务值比例的 90%与 80%所对应的栅格值作为水源涵养服务功能评估分级的分界点，最能体现分级特征。利用ArcGIS重分类工具，将水源涵养功能重要性分为三级，即极重要、重要和一般重要[15]。提取极重要区域并剔除面积不足1 km2区域作为安徽省水源涵养能力生态安全格局分布。

水源涵养极重要区域主要集中在长江沿线、皖西大别山区、皖南黄山区、巢湖以及中北部部分河流湿地，总体呈现南多北少的分布状况（图3）。长江作为水生态安全控制点，是水源涵养极为重要的组成部分;大别山区和黄山区有充足的林木可以锁住大量水分，是全省水源涵养能力的重点区域;巢湖以及其他河湖湿地对水源涵养有着较为明显的调节作用。

2.1.2 水土保持功能重要性分析。

式（2）计算出水土保持功能分布后，数据归一化导出服务值属性表并按从大到小顺序排列，经试验分析确定将总值占水土保持服务总值比例的60%与40%所对应的服务值作为生态系统服务功能评估分级的分界点，并将其重分类为三级重要性。提取极重要区域并剔除面积不足1 km2区域作为安徽省水土保持能力生态安全格局分布。

水土保持极重要区域集中在皖西大别山区、皖东南黄山区和皖东滁州市中部山区（图4），该区域林木较多，尤其适宜乔木灌丛等植被生长，植被对水土的依附性又可以保证水土不被流失。

2.1.3 生物多样性保护重要性分析。

计算安徽省生物多样性保护功能服务值，经试验重要性分界点选择80%和95%，重分类后得出生物多样性保护重要性分级。提取极重要区域并剔除斑块面积不足1 km2区域，作为安徽省生物多样性保护生态安全格局分布。

生物多样性表现为动植物种类数量和集聚程度，安徽省从北向南依次为平原、丘陵、山地分级，平原地区生物多样性明显弱于南部山地，黄山区为全省生物多样性高值区，总体呈现南多北少分布态势，西部大别山南麓出现斑块状生物多样性高值区（图5）。山区林木茂盛，自然环境丰富多样，适宜更多生物生长繁殖。

2.2 水土流失敏感性分析

根据通用水土流失方程的基本原理，依据评估模型计算得水土流失敏感性指数，并对其归一化处理，结合水土流失强度分级标准和专家知识[33-34]，将安徽省水土流失敏感性指数分级为一般敏感、敏感、极敏感，得出安徽省水土流失敏感性分级。提取极敏感区域并剔除斑块面积不足1 km2区域作为安徽省生态敏感性格局。

综合全省水土流失敏感性指数来看，极敏感区域主要集中在皖南黄山区，皖西大别山部分区域，以及大别山、黄山余脉延伸区，整体呈现东南多西北少的趋势（图6）。黄山和大别山作为重要旅游资源地，在旅游开发过程中难免存在生态与发展间的矛盾，水土流失问题虽不甚于我国西北地区，但从全省来看，两大山区应注重保持水土等相关举措，严防水土流失敏感性问题扩大化。

2.3 景观生态指数分析

2.3.1 景观破碎度与景观完整性分析。

景观分裂指数（DIVISION）与斑块密度指数（PD）类似，所得均可显示区域内景观种类的复杂程度。通过对1 km、500 m、250 m的窗口调试分析，250 m窗口大小较能反映全省景观的分布情况。在Fragstats软件中得出全省DIVISION和PD的空间分布格局（图7）。

全省范围内分裂指数和密度指数分布相似。高值区主要集中于淮河支流区域，西部大别山北麓和南麓，长江沿线，以及马鞍山市、芜湖市、宣城市交汇地带（图7）。对照安徽省土地利用图，发现北部高值区成因在于河流周围均为耕地而凸显;长江沿线和大别山南北麓高值区内土地利用类型复杂多样，河流、农田、城镇以及其他土地利用类型均有分布，三市交汇地区土地利用更为复杂。剔除分布格局中土地利用类型为建设用地和农田的区域，按照从小到大顺序排列其余指数值，计算累加指数值。由于景观分裂指數与斑块密度指数表达意义相反：分裂指数越高，区域内景观分布越离散，破碎度越高;密度指数越高，区域内景观种类越多，景观完整性越低。通过对指数值占累加值比例的40%、60%、80%作为分级关键点试验分析后，发现60%～100%区域内指数值最能作为景观破碎度高值区和景观完整性低值区（图8）。

剔除农田和建设用地后，景观破碎度高值区和景观完整性低值区出现在大别山南麓，长江沿线，马鞍山、芜湖、宣城三市交汇地带（图8），表明当前区域内林地河湖等土地利用类型分布较为分散，景观特征表现为区域内同种土地利用类型独立分布，其中长江沿线极值区成因与长江沿线斑块状分布的沿江防护林有关。省内其他地区的极值区与单独出现的河湖湿地有关。

2.3.2 蔓延度分析。使用Fragstats计算出安徽省蔓延度，高值说明景观中某种优势斑块类型形成良好的连接性，反之则表明景观是具有多种要素的分布格局。剔除城镇农田所在区域，经试验将蔓延度值低于60%的区域选为低值区。

相比于DIVISION和PD指标，蔓延度低值区有更广的分布范围，除大别山和黄山地区外，东部滁州市境内山区也存在低蔓延度地区（图9），表明该地区景观类型较为丰富，各斑块间连接性较低。