-
近几年,随着城镇化的推进,大量地表被人工建设用地侵蚀,导致生态斑块呈现出破碎化、孤岛化、异质化特征,降低了生态景观的连通性,正常的物种迁徙受到阻碍,不利于区域生物多样性的维护和可持续发展[1-2]。经研究论证,单纯将物种栖息地进行保护很难缓解生境破坏问题[3]。搭建生态网络可贯通被割裂的生境斑块,促进物种在残留生境斑块间的迁移、扩散和种群基因交流,维持区域内的物质动态循环和能量流动,可有效改善区域自然生态系统服务功能和保障区域生态安全[4-7]。
自20世纪70年代起,对生态网络的研究已备受国内外学者关注,在模型搭建、生态系统保护、景观连通性分析等方面已有了一定的成果[8]。生态网络的构建已形成了“源地选取—生态阻力面构建—生态廊道提取”的常见模式[9],研究手段也不断向着定量化发展,包括最小累积阻力模型、形态学空间格局分析、InVEST模型、生态系统服务价值评价、电路理论、网络分析等研究手段。其中,形态学空间格局分析(MSPA)是生态源地提取的重要手段,该方法基于数理形态原理提出,考虑斑块的整体连通性,通过土地利用二值图对栅格数据在景观结构上进行准确分类,科学识别物种重要栖息地[10-11]。最小累积阻力模型(MCR)通过构建生态阻力面来提取潜在生态廊道,能综合分析各因素与生态廊道连通性的关系,有较好的兼容性和普适性,可结合重力模型、图谱理论等对廊道相对重要成程度进行分级提取,是运用较为广泛的方法[12-14]。众多学者已通过耦合MSPA及MCR模型构建生态网络,为区域的生态系统修复规划、土地整治及城市健康发展提供了科学参考[13,15-17]。攀枝花市是长江上游的重要组成部分,区域自然资源丰富,是一座资源型城市,长期的矿业经济活动给攀枝花市生态环境造成了严重的压力,露天采矿区的增加,地质灾害风险加大,区域内水土流失面积增加[18],生境斑块破碎化加剧,面临的生态问题较为严峻[19-20]。在推进生态空间修复工程的背景下,构建攀枝花市生态网络对于维护我国西南地区区域生态系统稳定及筑牢长江上游生态屏障有重要意义。本文运用MSPA方法,结合斑块重要性( dPC) 、可能连通性( PC) 、整体连通性指数(IIC)指标科学提取市域重要生态源地;考虑人为及自然因素,依据坡度、土地利用类型、道路交通、植被覆盖等因素构建生态阻力面;结合MCR及重力模型搭建攀枝花市生态空间网络架构,多角度探究生态网络优化方向,科学提高攀枝花市域整体生态安全水平。
-
数据主要包括:攀枝花市行政边界数据、攀枝花市景区分布点数据、2021年土地利用数据、植被覆盖数据、路网矢量数据、12.5 m分辨率数字高程数据。其中,土地利用数据以2021年Sentinel-2遥感影像为基础,数据来源于ESRI官网的数据下载中心(www.arcgis.com),该数据分辨率为10 m,运用深度学习模型生成,使用超过50亿个手工标记的Sentinel-2像素进行训练,本文将其整合提取为森林、灌木、水体、建筑用地、耕地和裸露地表6种地类;植被覆盖数据基于2020年Landsat-8 OLI遥感数据(来自地理空间数据云http://www.gscloud.cn/),通过ENVI和ArcGIS软件对遥感影像进行预处理,经过FLAASH大气校正、影像镶嵌、裁剪等预处理,采用NDVI指数来表示攀枝花市植被覆盖度;路网矢量数据源自OpenStreet-Map,提取攀枝花市铁路、高速路、国道、省道4类道路信息;12.5 m分辨率数字高程数据由ALOS卫星相控阵型L波段合成孔径雷达(PALSAR)采集。
-
根据攀枝花市土地利用数据,将森林地类作为MSPA分析的前景,其他用地类型归为背景值,结合研究区大小,在满足数据精度要求的前提下,导出二值化TIFF格式文件,栅格单元大小为15 m×15 m,导入guidostoolbox软件进行MSPA景观格局分析,通过对栅格图像的斑块空间关系进行计算,得到7类景观类型数据,依次为核心区、桥接区、边缘区、孔隙、岛状斑块、支线、环道区,并对其面积占比进行统计分析。核心区即是对生态系统完整性起维护作用的重要生境斑块,一般来说,面积越大,斑块的连通性越强、生态质量越高[21],故选取面积前30的斑块作为攀枝花市潜在生态源地。
-
景观连通性是可以定量描述景观之间相互连接和延续的一种测定指标,可作为生态源地识别及评价的重要参数。采用Conefor2.6进行运算,整体连通性指数(IIC)、选用常用的可能连通性指数(PC)、斑块重要程度指数(dPC)3个景观指数,能直观看出不同斑块间的连通指数水平,公式如下:
式中:n为核心区个数;nlij指斑块i和j之间的连接数量;a i和a j为核心区i和j的面积;P ij*是物种在斑块i到斑块j之间扩散的最大可能性;AL表示研究区景观的总值;PC表示某一景观的可能连接度指数,PC值域范围在0到1,值越大,说明景观斑块的连通程度越高;dPC表示斑块的重要程度,PCremove表示删除某斑块后的可能连接度指数。
阈值的设定会影响各个指数的大小,当阈值小于不同核心区之间的距离时,会被认为连通性较差[22];因此,综合研究区的大小、核心区之间的连通情况及相关研究等,连通阈值设在2000~2500之间较为合理,经过测试,最终设置为2500,连通性概率设为0.5。最后选用dPC>4且IIC>2.5的10个核心斑块作为研究区生态源地,其余作为潜在生态源地。
-
生物在生境变化及面临生存竞争的情况下会进行迁徙,但随着生态源地间的距离扩大,在迁徙中将面临各种生态阻力。攀枝花市地势起伏较大,地形因素将迫使生态廊道沿地势平缓的区域布设。同样,交通要道也基本沿河谷延伸,为避免生态廊道的部署与主要道路重叠,在用地类型的基础上,加入道路阻力因子。因此,综合考虑自然和人为因素,分别选取土地利用类型、坡度、高程、NDVI、距铁路距离、距国道距离、距主要道路距离(包含省道及高速路)作为生态阻力因子,在GIS中通过叠加分析构建综合生态阻力面。参照研究区实际情况,将各因子划分为5级或6级,并赋予阻力值,通过层次分析确定权重(见表1)。
阻力因子
Resistance factor分级指标
Grading index阻力值
Resistance value权重
Weight阻力因子
Resistance factor分级指标
Grading index阻力值
Resistance value权重
Weight土地利用类型
Land use type森林Forests 1 0.2 国道/m
National road /m>2000 1 0.12 灌木Shrub 5 1500~2000 30 水体 Water body 20 1000~1500 50 耕地 Arable land 40 500~1000 70 裸露地表 Bare surface 60 <500 100 建筑用地 Building area 100 铁路/m
Railway /m>2000 1 0.12 坡度/°
Slope gradient /°<5 1 0.18 1500~2000 30 5~15 10 1000~1500 50 15~25 50 500~1000 70 25~35 70 <500 100 >35 100 主要干道/m
Major arterial /m>1000 1 0.1 DEM/m <1400 1 0.14 800~1000 30 1400~1700 20 500~800 50 1700~2100 50 300~500 70 2100~2600 70 <300 100 >2600 100 NDVI >0.5 1 0.14 0.43~0.5 20 0.37~0.43 50 0.27~0.37 70 <0.27 100 Table 1. Ecological resistance factor assignment
-
最小路径法指物质在流通中从起始地出发穿过各类生态阻力而到达目的地的最小消耗路径,是物种迁徙及相关生态物质交流扩散中最有效的路径。通过最小累积阻力模型可反映出物质能量及生物物种在克服各类阻力时在景观面之间流动的可能途径及变化趋向。计算公式如下:
式中:Dij是一个点j穿过景观类型到另一个点i的空间距离,Ri是生态物质经过空间i所需克服的阻力值。
基于前面分析提取的生态源地,在ArcGIS 10.8软件通过“Cost Path”工具,以生成的生态阻力面作为每个生态源地的累计成本面,构建各个源地间的最小成本路径线。结合重力模型,对生态源地斑块之间的相互作用矩阵进行分析,可定量评判各个板块之间相对重要程度,提取出攀枝花市相对重要的生态廊道。具体计算方式如下:
Gij表示斑块i和j的相互作用力大小,Ni和Nj表示两斑块的权重值,Dij是i和j斑块间潜在生态阻力的标准值,P i为斑块i的生态阻力值,S i为斑块i的面积,Lij是斑块i到j之间的累积生态阻力值,Lmax表示最大生态阻力值。
-
通过MSPA景观格局特征分析发现攀枝花市具备一定的生态基质条件,如表2、图1所示,提取出的前景面积为3693 km2,占攀枝花市总面积约49%。其中,核心区面积为2728.71 km2,占前景总面积73.88%,占全市总面积的36.79%。核心区大面积的斑块集中分布于研究区北部,南部的核心区则主要集中分布在西南方向,这些区域的山体及河流较多,受人为因素影响相对较小;中部沿着金沙江两岸是主要城区所在地,该区域建筑用地面积占比大,受人为因素影响最大,核心区分布的斑块稀少,形成了打破南北贯通的“断裂带”。其次面积占比最高的是边缘区,占前景总面积14.86%,边缘区作为核心区与外部区域的缓冲地带,起着保护核心区的关键作用,应维护该区域的稳定性。岛状斑块可直观表现生态斑块的破碎度,占前景面积比例仅为1%,说明核心区的斑块整体性较好。桥接区和环道是斑块之间的连通途径,支线具备一定的连通性能,孔隙则是核心区内部空白区域,这四个值均不大,占前景的比例依次为1.79%、2.02%、2.452%、4.00%。整体来看,研究区核心区斑块之间的独立性较强,但斑块内部存在一定量的孔隙,斑块形状不饱和,容易被外部因素干扰。
景观类型
Landscape type面积/km2
Area /km2占前景面积百分比
Percentage of foreground area占总面积百分比
Percentage of total area核心区 Core zone 2728.71 73.88% 36.79% 支线 Branch line 90.55 2.45% 1.22% 边缘区 Marginal zone 549.02 14.86% 7.40% 岛状斑块 Island-shaped patch 36.89 1.00% 0.50% 桥接区 Bridging zone 66.28 1.79% 0.89% 环道区 Circuit zone 74.45 2.02% 1.00% 孔隙 Pore zone 147.59 4.00% 1.99% Table 2. Areas of landscape types in Panzhihua city based on MSPA
对照攀枝花市的景观分布数据,提取的生态源地基本囊括了二滩国家森林公园、格萨拉生态旅游区、绿石林景区、菁河瀑布、白坡山生态保护区、啊喇自然生态旅游风景区、迤沙拉民族生态旅游区、普威镇绿野花乡景点,总面积1686.38 km2,(见图2),生态源地集中分布在北部,可见北部生境斑块的连通性最好。
-
结合生成的生态阻力面(见图3),发现攀枝花市中部、东南部、北部顶端部分的阻力值最大,阻力值较小的区域主要位于西北部及南部区域,这是由于攀枝花市是沿江流及河谷发展的城市,以水系为轴线,向外扩展城市用地、道路建设用地及一些工矿用地,人为活动的加剧影响了物种的迁徙环境;北部高山多,坡陡,加大了生态廊道建立的阻力;河流可为生物迁徙提供一定的能量保障,地形平缓、人类活动少、沿途有河流经过的区域生态阻力值小。基于最小阻力模型,提取出45条生态廊道,删除重复、冗余的生态廊道,共得到28条有效生态廊道。
生态源地斑块的相互作用力越大,意味着生态廊道的建设及维护更有价值,对物种的物质能量交换也更有利。通过重力模型构成源地间相互作用强度矩阵(见表3),对相互作用力大小进行分级,大于5的一级生态廊道7条,在2至5间的二级生态廊道10条,小于2的一般生态廊道11条,最终构建出攀枝花市生态空间网络(见图4)。可看出,源地4与源地6斑块间的相互作用力最大,为33.93,这是由于两个生境斑块的环境质量较好,两地的空间距离较近,生态累积阻力小,为生物提供了良好的迁徙环境。源地10与源地1之间的相互作用力强度最小,仅为0.27,这是由于源地1在攀枝花市北部顶端,而源地10则在最南部,相隔距离远,阻力值大。一级生态廊道主要分布在源地2,3,4,5,6之间,该区域位于西北部,呈环状、线状相连,与多个生态景点相关联,生态价值极其重要,应重点维护此区域的生态网络构造,可汇聚成更大的生态源地,作为区域生态保护的核心。其余二级生态廊道多分布在北部,串联了该区域的生态源地斑块,源地间联系紧密,需对这部分生态廊道进行重点维护,为攀枝花市北部生态屏障建设提供保障。整体来看,结合生态网络核密度分析(见图5),生态廊道核密度值在空间分布上不均匀,生态廊道分布最为密集的区域集中在西北部,但中部、南部区域的廊道分布密度较低,且生态网络未能扩散至整个研究区,生态流动性整体偏弱。
生态源地编号
Number of ecological source1 2 3 4 5 6 7 8 9 10 1 4.96 4.86 2.18 4.30 2.31 0.86 0.37 0.73 0.27 2 12.07 5.26 2.44 3.68 0.70 0.42 1.09 0.35 3 15.36 5.08 13.73 1.14 0.61 1.86 0.49 4 3.01 33.93 1.10 0.88 3.30 0.69 5 5.47 2.14 0.78 1.84 0.43 6 1.61 1.01 4.35 0.78 7 2.79 2.33 0.65 8 3.62 1.97 9 1.14 Table 3. Intensity matrix of interaction between ecological corridors
-
生态源地具备一定的生态服务范围,以提升研究区整体生态系统效益为目标,通过识别已有生态源地辐射区的较大缺口来锁定新增源地区域[23]。以10个生态源地的几何中心点为圆心,根据源地面积大小依次设置20 km、17 km、15 km半径作为源地缓冲的辐射距离,发现生态源地辐射区在3个区域有较大空缺,主要在东北部米易县、中部金沙江流域、西北部格萨拉景区附近。根据生态源地特征,将5处潜在生态源地新增优化为生态源地,通过删除重复冗余廊道,新增生态廊道16条,实现生态廊道的优化(见图6)。生态节点在生态网络中起连接邻近生态源地、为生物空间运动提供休息场所的作用。结合源地的空间分布,选取生态廊道交接点作为生态节点,并剔除位于生态源地和距离较近的生态节点,共得到26个生态节点。
-
交通道路网对生态网络的破坏较大,并形成断裂点,阻碍物种的交流。提取攀枝花市道路级别较高的铁路、国道、高速公路矢量数据,与生态廊道相交后识别出64处生态断裂点(见图6)。需注重生态断裂点的修复,如建设动物通道、天桥等设施。
-
根据攀枝花市生态空间网络生态结构特征,形成“一心一轴三片区、两带多廊道”的发展布局对策(见图7)。一心为北部连接成片的重要生态源区,应进行重点保护,划定生态保育区,严格管控人为干扰行为;一轴为跨过攀枝花市主城区的连通南北生态源地的轴线,在此轴线周边适当修建绿地公园等设施,加强绿地斑块的连通性,打造良好人居环境;三片区分别为东北部米易生态片区、西北部山地生态片区、南部迤沙拉旅游生态片区,这几个片区可作为区域生态中心,继续加强区域生态网络的流动性;两带是指安宁河流域及金沙江流域,是研究区物种迁徙的生命线,应加强水流带沿线的生态保护及治理工程等工作;多廊道指构成生态网络的各级生态廊道,维护着攀枝花市生态系统的整体稳定性,穿过山区、林区的生态廊道应加强区域生态修复及森林保育等措施,靠近或经过主干道路、地形平缓区域的生态廊道应注意生态断裂点的修补,并完善道路绿地的建设。
Research on Construction and Optimization of Ecological Spatial Network of Panzhihua City Based on MSPA-MCR
doi: 10.12172/202203300004
- Received Date: 2022-03-30
- Available Online: 2022-07-21
- Publish Date: 2022-10-26
-
Key words:
- Panzhihua city /
- Ecological network /
- Morphological spatial pattern analysis /
- Minimum cumulative resistance model /
- Gravity model
Abstract: Ecological spatial network can strengthen the connectivity between fragmented habitat patches and landscapes, which is particularly important in the ecological protection and restoration projects of mountains, forests, fields, lakes and grass. In this paper, considering the natural and human factors, Panzhihua city was selected as the study area, and the ecological network was constructed by comprehensive application of morphological spatial pattern analysis (MSPA) and minimum cumulative resistance model (MCR). The important ecological corridors were quantitatively evaluated by gravity model, and the optimization countermeasures were put forward. The results showed that the ecological network of Panzhihua city was composed of 10 ecological sources and 28 effective potential ecological corridors. Combined with gravity model, 7 first-class ecological corridors, 10 second-class ecological corridors and 11 general ecological corridors were extracted. The spatial distribution of the core density of the ecological corridor was unevenly distributed, mainly distributed in the northwest, and the ecological mobility was weak as a whole. The ecological network was optimized by adding ecological sources, identifying ecological nodes and ecological breakpoints, and the countermeasures of ecological spatial network framework layout were put forward. The research can provide scientific references for the sustainable development of ecological environment in Panzhihua city.