随着生态文明建设不断推进, 我国生态环境保护越来越重视整体生态安全格局的建立与维护, 区域发展不仅要根据自身的需求, 同时要满足更高尺度乃至全国生态安全保障体系的要求, 由此形成了一批能够对区域生态环境起到重要生态支撑作用的区域, 该类区域以提供生态输出、维持更大尺度区域乃至全国的生态安全为区域自身发展的前提和目标, 注重人与自然间的自然性和谐^[1], 将这类区域称为“生态型区域”。党的十九大进一步强调“必须树立和践行绿水青山就是金山银山的理念”, 可见国家在生态文明建设上的决心, 在这一理念的指导下, 生态型区域如何转变发展模式, 探索一条绿色发展道路, 实现自身生态系统的完善与维护, 从而在区域上起到真正的生态支撑作用, 成为该类型区域亟待解决的问题。景观格局与过程是景观生态学研究中的核心内容^[2], 资源短缺、环境污染、生态恶化等诸多问题的出现均与景观格局演变密切相关^[3]。区域景观生态安全格局的构建不仅能够保障区域生态安全, 对于指导社会经济发展和生态环境保护均具有重要意义^[4-5], 是实现区域社会经济发展转型和生态环境建设的基础。然而, 目前通过合理设置生态用地, 有效提高城市生态服务功能和保障城市生态安全, 仍缺乏深入而又系统的研究^[3]。“源-汇”景观理论是探究不同景观类型在空间上的动态平衡对生态过程影响, 从而找到适合一个地区的景观空间格局^[2]。该理论自提出以来在非点源污染控制^[6-8]、区域空间规划^[9-12]、适宜性评价^[13-15]、城市生态安全格局构建^[16-20]等许多领域得以应用, 为建立能够保障生态服务功能得以有效实现的景观生态安全格局提供了理论支撑。

在京津冀协同发展的大背景下, 张家口市无论是在《国家主体功能区划》中, 还是在“河北省新两翼建设”以及《北京城市总体规划》中均作为生态保护的重点区域, 同时还是“国家山水林田湖生态保护修复”的试点区域, 冬奥会的申办更加剧了张家口生态环境建设的紧迫性, 如何建立合理的生态安全格局、保障其生态环境支撑和水源涵养的功能得以实现, 如何找准“绿色办奥”与本地发展的结合点, 使张家口市成为向世界展示山清水秀、蓬勃发展新中国的窗口, 这些促使张家口市成为我国生态型市域的典型区域。张家口市要实现对于区域生态环境的支撑, 走绿色发展之路, 首先必须要从总体上完成景观生态格局的合理布设, 构建区域景观生态安全格局, 并以此为基础确定社会经济发展和生态环境保护的目标和方向。其次, 张家口市目前整体社会经济发展相对滞后, 生态环境相对脆弱, 景观生态安全格局需要分批分步骤实施。本研究试图在“源-汇”理论的指导下, 构建以生态保护为主要任务的生态型市域景观生态安全格局, 识别不同等级的生态保护与生态修复区域及廊道等保护目标, 确定保障城市生态安全的重点生态区域, 并通过保护目标等级的划分, 确定区域生态保护建设的优先次序, 为实现欠发达地区分步骤建立和完善景观生态安全格局, 探索一条欠发达地区生态兴市、生态强市的路子提供范例。

1 材料与方法 1.1 区域概况

张家口市地处内蒙古高原与华北平原的过渡带, 以大马群山为界, 地貌类型分坝上、坝下两大地貌单元。截至2016年, 张家口市辖6区10县, 总面积35904.08 km², 全市总人口441.1万人。张家口不仅是承东启西、沟通南北的咽喉要地, 更是京津冀重要的水源涵养区和生态环境支撑区, 在地域文化上是晋文化、京畿文化、草原文化的连接带, 在生态系统上是内蒙古高原生态系统向华北平原生态系统的重要过渡带, 因此在区域生态安全框架中起到重要的生态屏障作用, 无论是首都经济圈还是河北新两翼建设, 都将张家口市作为重要的生态支撑区, 其生态安全格局不仅关系到区域的发展与定位, 同时对于整个国家的生态安全格局也起着至关重要的作用。张家口地区历史上曾是“风吹草低见牛羊”的生态圣地, 后因清朝末年开围放垦和建国初到20世纪80年来30多年伴随人口快速增长带来的有组织大规模开垦, 农业用地面积高达61%以上, 原本良好的生态环境遭到严重破坏, 近年来高耗水产业的兴起使原本脆弱的生态环境雪上加霜, 现状生态环境质量与支撑区域生态环境的要求有较大差距。加之该区域社会经济发展落后, 贫困人口众多, 改变现有发展模式, 分步骤建立区域景观生态安全格局显得尤为紧迫。

1.2 “源-汇”理论

“源-汇”理论最初源于大气科学, 后被引入到景观生态学领域, 该理论将空间格局和生态过程有机结合在一起, 为探索格局与过程的关系提供较好的理论基础^[21]。“源”景观是指在格局与过程研究中, 那些能促进生态过程发展的景观类型；“汇”景观是那些能阻止延缓生态过程发展的景观类型^[2]。不同的景观类型对整个区域所起到的作用不同, 景观类型按照对整个区域的物种流、物质流和能量流起到的作用分为“源”和“汇”, 同时一些景观类型起到传输的作用, 如林地、草地、湿地、水域等起到“源”的作用；另外一些景观类型如生物生产用地、建设用地等起到“汇”的作用；各种河流、绿色廊道等起到传输的作用。本研究中采用面积较大的景观组分和良好格局信息的林地核心斑块作为景观“源景观”^[22], 其余作为“汇”。

1.3 最小累计阻力模型

最小累计阻力模型是以“源-汇”理论为基础的模型, 最小累积阻力指扩展主体从源经过不同阻力的景观阻抗到达目的地(汇)所消耗的费用或克服阻力所作的功, 反映了一种潜在可达性^[9]。最小累计阻力模型是把不同景观类型相互转化和景观格局的改变看作景观要素通过克服阻力来实现, 它反映了景观要素空间运动的趋势^[23]。最小累积阻力模型(minimum cumulative resistance, 简称MCR)可以用来建立阻力面, 该模型考虑三个方面的因素, 源、距离和景观介质特征^[24], 基本公式如下^{[5, 25]}:

(1)

式中, D为距离, R为阻力系数, 最关键的过程在于如何确定扩展源和阻力系数^[9]。

扩展源：在景观中有一些能够促进生态过程发展的景观类型称之为“源景观”^[26], 生态源景观是由具有高生态服务功能的景观类型组成的景观组分, 常由林地、水域和草地构成^[27], 可将土地利用方式典型且相对稳定的景观单元作为“源”^[17]。张家口市为典型的生态型区域, 选取作为生态系统重要支撑的具有一定规模的草地、林地和水域作为扩展源, 借鉴以往研究成果^{[26, 28]}, 选取面积 > 20 hm²的空间上具有连续性的核心林地和草地以及面积大于5 hm²的水域作为源景观。

阻力面：根据以往研究以及张家口市自然地理特征, 本研究分别选取草地、林地和水域作为“源”, 耕地、园地、建筑用地、交通用地和未利用地作为“汇”, 根据土地利用的特征设定阻力值。本文设定1为最小阻力系数, 表示某一景观单位可以自由移动, 在区域为最适宜空间, 100为最大阻力系数, 景观单元无法移动, 为最不适宜空间^[21]。根据张家口市生态功能支撑区的定位, 将促进生态系统支撑的景观类型设定为低阻力值, 对于生态系统源类型扩张起到阻碍作用的景观类型设定为阻力高值, 具体而言, 对于生态系统类型的保护建筑用地和地交通用地阻力最大, 取值100、未利用地次之, 取值70, 耕地和园地具有不同生态系统支撑功能, 其阻力分别为60和50, 该区域草地生态系统面临较为严重的退化风险, 因此对于生态系统的支撑具有一定的阻力, 取值10, 林地和水域阻力最小, 取值1^[21]。根据公式(1)可以通过ArcGIS的cost-distance功能实现阻力面的构建。

1.4 数据来源

数据来源于张家口市国土部门提供的《2015年张家口市土地利用分类》栅格数据, 精度为1 km×1 km, 涉及土地利用类型包括耕地、园地、林地、草地、交通用地、水域、未利用地、建筑用地。生态保护地源于《河北省主体功能区规划》中的禁止开发区, 主要包括保护区、水源地、森林公园、湿地公园、重要湿地、风景名胜区等。

2 结果与分析 2.1 源景观分布特征

根据源景观识别方法, 三类生态源景观分布如图 1所示。草地源景观总面积10684 km², 占全市草地面积的79.20%, 其中天然牧草占比4%、人工草地占比16%、其他草地占比80%, 天然牧草主要分布在坝上地区, 包括行政区域主要涉及尚义、张北、康保和沽源县, 面积占比低, 草原退化严重, 破碎度较高, 坝下丘陵区主要以其他草地类型为主。林地源景观总面积7066 km², 占全市林地面积的64.07%, 其中有林地面积占比32%、灌木林占比29%、其他林地占比39%, 主要分布在坝下山地, 行政区域涉及崇礼、赤城、怀来、逐鹿县, 近年来张家口市对林地建设较为重视, 特别是申奥成功以来, 提出了明确的森林覆盖率要求, 但目前对于森林质量的建设还有待增强。水域源景观总面积657.79 km², 占全市水域面积的85.49%, 主要分布在张家口的几大流域, 包括永定河流域、潮白河流域、大清河流域和内陆河流域, 张家口市近些年的水资源环境状况堪忧, 水污染问题特别是面源污染问题严峻, 许多流域出现河道断流, 流域内水土流失严重, 区域水源涵养能力不足。

2.2 阻力面分布特征

通过MCR分析不同源景观的阻力分布情况(图 2), 草地阻力值分布范围为0—42562, 阻力高值区主要集中在人类活动干扰较高的市域周边以及蔚县, 坝上草原为天然牧草的主要集中分布区, 而坝上区域的草地源扩张阻力中等偏上区域分布广泛, 康保县、张北以及尚义县的北部阻力值较高, 使草地源扩张受到严重阻碍, 区域呈现高度破碎化。林地阻力值分布范围为0—37992, 阻力高值和低值区分布较为集中, 阻力高值区主要分布在坝下丘陵地区。水域阻力值分布范围为0—99180, 阻力高值和低值区分布最为集中, 呈现由向北向东南逐渐升高的趋势。将不同源景观阻力在ArcGIS中进行叠加, 叠加后阻力值分布范围为0—34648, 综合阻力高值区主要分布在东南区域, 大清河流域、潮白河流域和潮白河与永定河分水岭(图 2)。

2.3 区域生态功能分区

阻力值越大表明该区域生态要素特征越不稳定, 要进行生态修复, 反之则为生态特征较为稳定, 要保护与稳定^[21], 因此, 根据综合阻力值突变点将张家口市生态功能分成5个等级, 每个功能区具有不同的生态功能和保护重点, 从而形成张家口市生态安全的总体格局, 其分级结果如图 3所示。生态阻力低值区为生态保护区域, 高值区域为生态修复区域。生态保护核心区面积2057.821 km², 占区域面积的6%, 生态保护核心区在全市分布较为均匀, 构成了区域保护的核心框架, 其生态阻力值低, 具有能够支撑保护源扩展的足够大的“源景观”类型, 是建立保护区域的优先区域。生态保护重要区面积19661.82 km², 占区域面积的54%, 形成区域生态安全保障的重要支撑区域, 起到连接生态保护核心区以及作为生态保护核心区缓冲区的作用。生态保护过渡区面积8601.364 km², 占区域面积的23%, 是生态保护区域向生态修复区域过渡的重要区域, 对于生态型市域, 在社会经济发展允许的基础上, 应逐步建设为生态保护重要区, 增大森林、草地和湿地等生态源景观的面积, 降低区域破碎度, 提高生态系统的数量和品质；生态修复重要区面积5444.50 km², 占区域面积的15%, 该区域受到人类干扰频繁, 区域生态支撑能力不足, 甚至遭到破坏, 是进行生态修复的主要区域, 重点进行水土保持、水源涵养等生态修复。生态修复极重要区648.052 km², 占区域面积的2%, 生态阻力最高的区域, 同时也是优先进行生态修复的区域。

2.4 重点生态修复区

为了确定生态修复区的修复优先次序和修复重点, 根据阻力值的分布确定突变点, 将草地阻力高于8000的区域作为草地阻力生态修复的重点区域, 草地阻力修复重点区域面积1538.64 km², 占全市面积的4.18%；将林地阻力高于10000的区域作为林地阻力生态修复的重点区域, 林地阻力修复重点区域面积2121.17 km², 占全市面积的5.76%；将水域阻力高于20000的区域作为水域阻力生态修复的重点区域, 水域阻力修复重点区域面积4247.89 km², 占全市面积的11.54%, 各区域分布见图 4。

将不同源景观生态修复区域进行叠加分级, 结果如图 4所示, 修复区总面积8932.05 km², 占全区面积的24.27%。修复顺序按照林草水修复区7.53 km²、林草修复区519.09 km²、林水修复区99.18 km²、水草修复区398.55 km²、林地修复区2121.17 km²、草地修复区1538.64 km²和水域修复区4247.89 km²。结合生态功能分区生态修复重点区域的分布范围, 其重点修复顺序为水域修复区>林草水修复区>林草修复区>林水修复区>水草修复区>林地修复区>草地修复区。

重点生态修复区主要集中在永定河流域, 永定河流域最终汇入官厅水库, 是北京的重要水源涵养地, 特别是永定河重要分支的洋河流域, 其流经张家口市城市发展的中心区域, 面临的生态修复压力巨大, 该区主要为林地、草地和林草修复区；另一个重要的生态修复片区分布在永定河的另一个支流桑干河流域, 该区域同样为林地、草地和林草修复区；第三大生态修复区位于大清河流域, 其修复重点为水域、草地和草水修复区；此外还有处于永定河和潮白河分水岭的水土保持重要修复区、位于潮白河流域的水土保持和草地重要修复区以及零散分布于内陆和流域的草地生态修复区。

2.5 重点生态保护核心区识别

为了确定区域生态保护的优先次序, 根据生态保护核心区的分布范围, 分析现有生态保护地的分布范围, 如图 5所示。将张家口市现有保护地作为中心, 其10 km范围内作为生态保护辐射区域^[29], 与生态功能分区中的生态保护核心区进行叠加, 未在生态保护辐射区域内的核心区即为重点生态保护缺失区域, 如图 6所示张家口市生态保护区缺失严重, 重点生态保护缺失区域总面积达1555.77 km², 占生态保护核心区面积的75.60%。重点生态保护缺失区域大体分为三个片区, 坝上区域主要涉及到尚义县、张北县、康保县和沽源, 以及崇礼和赤城的坝上区域, 第二片区域为市中心周边区域, 第三片区主要分布在阳原县和蔚县。坝上区域主要为湿地链保护的不完善和草原生态系统保护的严重缺失, 市域周边片区主要为水土保持与森林保护体系的缺失, 坝下丘陵区主要为湿地保护的缺失。

2.6 保护廊道构建

生态廊道是生态流在源间运行的高速通道, 有利于生态系统的物质和能量交换, 对于维持生态系统的稳定性、保障生态安全至关重要, 与生态源景观形成环状的空间格局为最佳^[30]。根据综合阻力面的分布大小, 识别现有保护地与生态保护核心区之间的最小加权路径, 结果如图 7所示。根据所得廊道分布特点, 识别主要廊道建设缺失区域, 建立完整廊道体系, 与生态保护核心区共同构成“四横两纵”保护网络, 如图 7所示。廊道建设重点在于保护现有廊道, 补充重点缺失廊道, 逐渐形成复杂的廊道联通网络体系, 根据张家口市现有廊道体系的分布特征, 重点识别出如图 7所示的3个廊道缺失区域, 1区分布在坝上区域, 是草原生态系统保护体系的重要连接通道, 结合前面分析得知该区域的草原生态系统保护面临着破碎化的严峻威胁, 建议在该区域建设国家草地公园, 使其成为草原廊道稳定的连接点；2区位于坝下山地区域, 重点以森林生态系统的联通建设为基础, 并水土保持和水源涵养功能保护相结合, 构建网络化的廊道体系, 其中2区廊道关系到冬奥会会场所在地崇礼的生态连接度, 可与奥运会场生态建设与主要道路绿化相结合, 进行优先建设；3区依托永定河流域水系联通体系, 建立城市通向生态保护区域的复合廊道体系。

3 讨论

本研究基于“源-汇”理论, 就张家口市景观生态安全格局构建进行了探讨, 建立了以五级生态功能区域为基础, 生态保护核心区和重点生态修复区为重点, “四横两纵”生态廊道为连接的“点-线-面”三位一体的立体保护格局。五级生态功能区划明确了全域的保护重点, 其中生态保护核心和重点生态修复区分别为区域重点“源景观”和“汇景观”, 共占整个区域面积的30%, 是支撑整个张家口市生态安全格局的基本骨架, 对重点“源景观”的保护和对重点“汇景观”的修复是构建区域生态安全格局的重中之重。

张家口市生态安全格局除占区域面积6%的生态保护核心区外, 生态保护重要区占区域面积的54%, 其生态系统类型主要为林地、湿地和草原, 同样需要不断加大保护力度, 防止区域破碎化程度加强, 有关研究表明, 区域内至少有50%以上的土地应作为保护区才有利于物种的空间运动不受景观破碎化的影响^{[5, 31]}, 因此维护张家口市生态型地区的支撑地位, 保护生态重要区的生态系统安全同样是核心。应不断加大退耕还林、还草、还湿力度, 削减区内基本农田面积, 为扩大森林、草原和湿地生态系统留足空间, 从而起到保障水源涵养功能得以有效发挥。研究不仅考虑了区域生态系统现状, 同时考虑了现有保护体系, 通过与现有保护体系进行对比, 识别了现有保护区域的缺失部分, 同时对于生态修复重点区域进行了等级划分, 这些有助于区域实行分类分批建立和完善区域生态保护体系。

通过本研究确定景观生态安全格局的构建流程如图 8所示, 首先要明确区域的发展定位, 生态型城市要为区域乃至国家的发展提供生态支撑, 因此必须优先保障其特定的生态功能得以实现, 从而能够真正起到对区域的生态支撑作用, 其发展定位要优先考虑生态功能得以保障, 生态支撑作用得以发挥。就张家口市而言, 作为重要的首都水源涵养区和生态环境支撑区域, 必须保障足够的生态“源景观”-即生态保护区域, 同时对于受到人类干扰较高的区域, 要重点考虑建立生态修复区, 从而保障生态功能得以正常发挥, 此外还要结合重点保护区和关键生态廊道建设等生态保护关键区, 形成“点-线-面”三位一体的立体生态安全保护格局体系。

“源-汇”理论为区域景观生态安全格局的构建提供了新的途径, 其对生态过程的关注使得区域生态安全格局的构建更能体现格局和过程的关系, 本研究基于该理论提出了张家口市景观生态安全格局的分步建设对策, 对于经济发展相对落后, 生态系统功能要求较高的生态型市域而言具有较高的现实性意义。本研究选取不同的景观类型作为源汇景观过程中, 其阻力等级的确定仅参考了现有研究, 对于不同阻力值的划定是否对生态安全格局产生影响尚有待进一步研究, 对于生态安全格局的有效性可以进一步在时间和空间上进行验证。

4 结论

张家口市是重要的生态型区域, 对于保护区域和全国的生态安全具有重要意义, 冬奥会的成功申办不仅成为张家口社会经济与生态保护的新契机, 更对张家口的生态环境提出了更高的要求。本研究针对张家口生态系统特征, 构建了张家口市景观生态安全格局, 形成了以五级生态功能分区为基础, 生态保护核心区和重点生态修复区为重心, “四横两纵”生态廊道为连接的“点-线-面”三位一体生态安全格局, 通过建立完整的生态保护体系、生态修复体系和廊道体系, 实现张家口市整体生态功能提升。

张家口生态保护体系建设重点以完善现有生态保护地为主, 坝上地区以建立草原保护区和湿地保护链为主, 坝下山地区域以森林保护区建设为主, 坝下丘陵地区主要为永定河上游的湿地保护体系建设。生态修复体系应重点加强永定河流域综合整治和生态修复, 加强上游水源涵养林建设, 注重加强幼龄林抚育管理, 推进低质低效林改造, 提升森林质量, 加大水土保持整治力度, 结合流域生态功能定位, 对于承担重要水生态服务功能的水源涵养区、河流廊道、重要湿地、河岸防护林等区域, 科学划定河湖管理保护范围以及生态保护红线, 保障并维持其自然结构特征和生态功能。充分利用河流水系和公路沿线建立生态廊道, 重点完善生态廊道缺失点的廊道体系建设, 结合生态保护体系, 逐步形成“四横两纵”的保护网络体系。

本研究结合张家口现有的景观格局体系与生态保护体系, 指出了张家口市生态保护的重点与缺口, 为城市未来发展格局调整提供了理论基础, 同时本研究提出了张家口市分步骤建设的重点, 即优先建设重点生态修复区、重点生态保护区和生态廊道, 逐渐扩大保护体系建设的范围, 这有助于欠发达地区逐步改善社会经济发展体系和筹集资金。此外本研究表明, 根据生态型区域的发展定位, 张家口市为了完成支撑区域生态安全的重任, 其生态缺口巨大, 在今后的生态环境建设中应着重建立生态补偿机制。