齐齐哈尔市是东北地区工业发展的重要枢纽,也是黑龙江省西部经济、政治及文化中心。地处黑龙江省西南部松嫩平原,东经123°—126°,北纬45°—48°,总体地势平坦,北部地势较高,东部及南部地势较低; 市域内包含嫩江等主要河流廊道,本土植被主要为草甸草原。齐齐哈尔市具有典型东北地区的温带大陆性季风气候,四季分明,春季气候干旱多风,夏秋两季短暂而多雨,冬季漫长而寒冷。近年来,哈大齐工业走廊相关建设和产业发展对城市群起到了极大的经济带动作用,齐齐哈尔市以非金属矿为工业发展特色,快速工业发展使得主城区建设用地急剧扩张^[20],侵占城乡生态空间,湿地、林地、草地斑块布局零散,主城区东南部与大庆市接壤的大面积湿地斑块逐渐呈现破碎化,水土流失较为严重,亟待修补与资源整合。

随着城镇化的推进,建设用地不断扩张,研究区生态斑块总面积减少,区域生境呈现破碎化的趋势。根据分析结果(表2),2003—2017年,研究区域建设用地面积不断增加,未利用地减少,各生态斑块也处于不断变化中。2003年耕地面积为1 553.11 km²,所占比重最大,湿地所占的面积比仅次于耕地,建设用地面积为354.28 km²,占总面积的8.59%; 到2017年,建设用地占研究区域总面积的10.66%,林地、草地、湿地面积均减少,区域生态环境亟待保护与修整。

表2 研究区域2003-2017年各类土地利用变化

以2017年土地利用类型栅格数据中的湿地、林地、草地作为MSPA的前景要素,耕地、建设用地、未利用地作为背景要素,将栅格转化为tiff格式的二值图,导入Guidos识别出7种景观类型,分别为核心区、桥接区、岛状斑块、环道区、边缘区、支线、孔隙,由表3看出,核心区面积为850.43 km²,所占比例最大,为54.16%,其次为边缘区,占比15.91%,而岛状斑块、环道区与孔隙所占比例较小,表明齐齐哈尔主城区景观构成以较大型斑块为主,且斑块具有边缘复杂、形态破碎等特征,支线的总面积为122.19 km²,占比为7.78%,大型斑块的外围分布较多支线,表明斑块较易与外界进行物质能量交流。

表3 不同景观要素的生态学含义及面积占比

结合前人研究^[9,11,23],设置斑块连通距离阈值设为500 m,连通概率为0.5,通过Conefor 2.6计算,根据结果进行景观分级,将核心区分为四级:生态源地dPC>0.1,一级核心区0.1>dPC>0.05,二级核心区0.05>dPC>0.01,三级核心区0.01>dPC。提取面积大于2.5 hm²的桥接区并分为3个等级:一级桥接区dPC>0.8,二级桥接区0.8>dPC>0.2,三级桥接区0.2>dPC。提取面积大于2.5 hm²的岛状斑块作为潜力节点,综合评价见图1。为有效区分,核心区dPC值取小数点后三位,结合文献^[13,16,23]选取dPC>0.1的16个斑块作为生态源地(表4)。

图1 基于dPC指数的核心区、桥接区及岛状斑块分级评价

表4 核心区(生态源地)景观连通性重要程度排序

基于最小累积阻力模型(MCR),根据16个生态源地与前期分析得出的研究区综合阻力面,在ArcGIS中计算每两斑块间的成本路径,16个生态源地间存在共计120条生态廊道。通过重力模型计算16个生态源地之间的相互作用矩阵(表5),参照前人研究^[16,37],将相互作用力大于5的廊道提取出来,作为重要生态廊道,共计14条,其余作为一般廊道,共计106条(图2)。表5表明,斑块10与斑块11间相互作用力最大,为326.32,说明此两斑块间距离较短,物种迁徙扩散阻力相对较小,斑块关联度较高,在未来规划中应加强斑块10与斑块11间生态廊道的围合保护。

图2表明,生态廊道集中于研究区东南部,东南部生境斑块面积较大,景观连通性较好,而西北部廊道体系连通性较差,结合

参考文献^[16],在研究区西北部重新选取具有重要生态意义的斑块,构建新廊道来完善研究区生态网络体系。根据前期MSPA分析结果,核心区具较高的生态效益,选取研究区西北部除生态源地外,剩余核心区中面积较大且连通性高的斑块作为补充生态源地。对于研究区边界处的核心斑块应重点考虑,因其位于边界,其与内部源地建设生态廊道后,廊道体系将全面覆盖研究区,此生态网络体系更为完善^[16]。最终选取7处补充生态源地(表6),与原有源地进行补充计算,得出新增133条规划廊道(图3)。

通过廊道网络闭合指数(α指数)、网络连接度指数(β指数)、网络连通率指数(γ指数)对优化前后的生态网络体系进行对比。结果如下:优化前的生态源地数(网络节点)为16,廊道总数为120,α,β,γ指数分别为3.89,7.50,2.86; 优化后的生态源地数(网络节点)为23,廊道总数为253,α,β,γ指数分别为5.63,11.00,4.02,生态网络连通性有较大提升,表明生态网络结构得到优化,规划的潜在生态廊道提高了区域生态源地连接水平,增强了区域生态系统稳定性。

表5 基于重力模型构建的生态源地相互作用矩阵

表6 补充生态源地生态属性及选择依据

图2 研究区域生态廊道识别与分级

图3 优化后研究区域生态网络体系

将优化后的“源地—廊道”生态网络体系进行分级缓冲区分析,生态网络体系50 m缓冲区内为极敏感区,面积为584.91 km²; 50～200 m缓冲带为高敏感区,面积为455.33 km²; 200～500 m缓冲带为中敏感区,面积为631.60 km²; 500～1 000 m缓冲带为低敏感区,面积为779.61 km²; 1 000 m缓冲带外的其他地区面积共1 652.72 km²,为不敏感区(图4)。极敏感区占研究区域总面积的14.18%,其比重大于高敏感区的11.04%,在未来生态规划中,潜在生态网络较近范围内的缓冲区应格外注重相关法律法规的制定,禁止随意开发建设。在低敏感区与不敏感区可根据实地情况,在避让生境斑块的前提下进行适度的开发建设。

将从Open Street Map网站和QGIS软件获取的矢量路网在ArcGIS中与生态廊道叠加,得出断裂点共计352个。其中,一级断裂点为生态廊道与铁路交汇点,共计116个; 二级断裂点为生态廊道与高速路交汇点,共计108个; 三级断裂点为生态廊道与国道交汇点,共计128个(图5)。生态断裂点在平齐铁路沿线上较为集中,在未来生态保护工作中需注重在相关生态廊道断裂点处修建桥上、桥下通道,为物种扩散提供通道; 可结合实地情况建设立体绿道活动空间,调节区域小气候,满足居民日常活动需要。

图4 生态网络体系缓冲区评价

图5 研究区域重要道路与生态断裂点分布

构建并优化生态网络能够增强区域生态系统的物质能量流动,缓解工业城市在发展建设中对生态环境产生的破坏,对城乡生态可持续发展具有重要意义。本文通过探究齐齐哈尔市主城区2003—2017年的土地利用变化,得出现阶段研究区域存在的生境斑块减少、斑块破碎化等生态问题,为生态网络的构建与优化提供建设方向,在东北重工业城市因地制宜的生态规划策略方面进行了一定的创新。基于形态学空间格局分析(MSPA)方法与景观连通性识别出生态源地,并结合区域景观类型与地形要素构建阻力面,基于最小累积阻力模型(MCR)识别潜在生态廊道。通过重力模型计算对廊道进行重要性分级,通过相关网络指数科学评价生态网络体系优化结果,并构建生态网络体系缓冲区作为研究区域生态敏感性分区,为后续相关部门的生态规划工作提供建设参考。研究结果表明:齐齐哈尔市主城区在城镇化进程中,生态斑块逐渐破碎化,生境面积减小,边缘趋于规则。核心区生态源地集中于研究区东南部的自然湿地以及穿过主城区的嫩江支流,原提取出生态源地共16处。重要生态廊道共计14条,一般生态廊道共计106条,主要集中于研究区东南部,在东北部亦有少量分布,而西北部的生态廊道存在连通水平不均、缺乏闭合环路等问题,生态网络体系有待完善。通过计算补充了7个生态源地与133条规划生态廊道,根据网络指数计算得出,优化后的生态网络结构连通性有了较大提升。优化后的生态网络体系共计23处生态源地、253条生态廊道,旨在为区域生态保护区划定与防护林廊道建设提供参考。根据网络体系缓冲区的生态敏感性分区得出未来规划中的重点保护区域与适度开发区域,根据352个生态断裂点布局得出立体绿道拟修建位置,为齐齐哈尔市相关生态规划设计策略的制定提供选址与方法借鉴。

结合前期分析,研究区域西北部耕地斑块面积较大,作为人为活动较强的区域,可能对斑块内小型林草地产生干扰破坏,应加大对此类小型斑块的保护力度,通过科学补植修复破碎林草地生境,提高区域生态系统的物质能量流动与生物多样性。本研究对东北老工业城市的生态网络构建与优化仍处于尝试探索阶段,整体研究过程中存在以下不足:首先,源地选择没有综合考虑到研究区内不同物种的生存规律差异,源地的扩充还需进一步斟酌; 其次,阻力面的构建缺乏对城市人口迁移等影响因子的考虑。希望在未来研究中能够不断完善,为齐齐哈尔等工业城市的绿色建设提供更为科学的生态网络构建与优化策略。