快速城市化深刻影响城市生态系统的功能、结构及其空间演化过程^[1-2]，尤其是大城市城区生态斑块连通性差、城郊过渡区的生态斑块破碎等问题突出。良好的生态安全格局既可以兼顾生物多样性保护又可以促进经济可持续发展，因此科学构建生态安全格局对于城市生态系统的健康和发展意义重大^[3]。

国内外针对生态安全格局的研究较为丰富，基本形成“识别源地—构建阻力面—廊道构建—安全格局构建”的主流方法^[4-5]，即“斑块—廊道—基质”的构建模式。生态源地是生物多样性较好且生境质量较高的生态用地，以往研究中源地识别直接选取自然保护区^[6]，识别地筛选标准相对单一，可采用生态系统服务功能叠加分析来识别^[7-8]，最小累积阻力模型(minimal cumulative resistance, MCR)能较好模拟生态流迁移扩散，广泛应用于阻力面构建和生态廊道模拟，以往生态阻力面构建直接采用土地覆被赋值方法^{[5, 9]}，近年来已有研究引入同类地块的空间异质性，基于土地覆被将不同的空间数据引入来修正生态阻力因子^[10-12]。北京市地处快速城镇化与工业化背景下，属超大规模城市，人口密度大^[13]，长期面临经济发展与资源环境难以协调的压力，城市发展空间回旋余地小、水资源和土地资源供需矛盾突出，区域生态格局趋向孤岛化。作为中国向世界展示绿色发展的一个窗口，北京的城市生态安全需要考虑可持续发展所面临的生态环境压力和建设发展的双重问题，因此如何有效协调城市发展和生态保护之间的矛盾，增强生态环境承载力，保护生物多样性，实现人与自然和谐共处，已成为当前北京市面临的重要挑战。本文叠加水源涵养、水土保持、防风固沙和生物多样性维持等服务功能，结合各类自然保护地范围提取生态源地，然后综合考虑北京市土地覆被、生态、交通可达性和社会经济因子，利用最小累积阻力模型生成北京市生态源地扩张阻力值和城镇用地扩展阻力差值，划分生态区、城镇建设区和生态城镇临界区，构建适合北京城市区域环境特点和社会发展水平相对完整的城市生态安全体系，以期为有效控制城市扩张发展生态安全底线^[14]和城市生物多样性保护提供科学参考。

1 研究区概况

北京市地处华北平原北端，地理坐标为39°26′—41°03′N，115°25′—117°30′E，总面积16 140.54 km²，西北部山地、东南部平原，平均海拔368 m。北京市属暖温带半湿润大陆性季风气候，年平均降水量548.86 mm，全年无霜期180~200 d，年平均气温为12.77 ℃，境内有大清河、永定河、北运河、潮白河和蓟运河五大水系，共100多条支流。全市主要包括森林、灌丛、草地、湿地、农田等5类生态系统，占全市总面积的81%。区域地带性植被为落叶阔叶林和温带针叶林。

2 研究数据与方法 2.1 研究数据

2019土地利用数据来源于国家遥感中心(www.nrscc.gov.cn)，空间分辨率为30 m；NDVI数据来源于NASA官网(http://ladsweb.modaps.eosdis.nasa.gov)；交通路网数据来自北京大学城市与环境学院地理数据平台(http://geodata.pku.edu.cn)；气象数据来源于中国气象科学数据共享服务网(http://cdc.cma.gov.cn/satellite)；DEM数据来源于地理数据空间云平台GDEMV2数据集，空间分辨率为30 m；土壤数据来源于中国西部环境与生态科学数据中心网站(http://westdc.westgis.ac.cn)，分辨率为1 km；DMSPOLS夜间灯光影响数据来源于网站(National Geophysical Data Center, NGDC)，研究应用非辐射定标的夜间灯光影响数据集DMSP/OLSZ作为社会影响因子之一建立生态源地扩展阻力面。其中风因子和土壤湿度因子来源于网站(http://cdc.cma.gov.cn)区域内的27个国家台站的日均风速、日照时数、温度、降水等数据来计算完成。

2.2 研究方法 2.2.1 生态源地提取

研究评价了北京市水源涵养、水土保持、防风固沙、生物多样性4个生态系统服务功能，评价方法参考《生态保护红线划定技术指南(2017)》，基于GIS栅格分类功能，采用自然断点法将评价结果分为4级：一般重要区域、中等重要区域、较重要区域、极重要区域^[15]。对北京市4类生态系统服务做标准化处理，加权后得到该区域生态系统服务重要性栅格图，采用熵权法对4类生态系统服务赋权重，提取生态系统服务重要性为第4级的区域，剔除破碎小斑块(面积＜20 km²)，并矢量化提取的斑块，得到北京市生态系统服务功能极重要区域空间分布；再选取北京市现有自然保护地(自然保护区、湿地、森林公园等)，基于GIS工具箱合并功能(uion)，叠加生态功能极重要区域，最终生成北京市生态源地。

(1) 自然保护地选取。北京市共有5类79处自然保护地，研究选取北京市3类63处自然保护地和7处城市公园(表 1)。选取的保护区面积总共3 948.75 km²，其中城市公园面积最小，为17.713 km²，自然保护区面积最大，为1 402.23 km²，森林公园595.31 km²，湿地、湖泊面积184.21 km²(图 1)。

(2) 生态系统服务评价。基于生态环境部《生态保护红线划定技术指南(2017)》生态红线划定方法，选取生态服务功能(水源涵养、防风固沙功能、水土保持同能和生物多样性维持与保护功能)4个指标^[16]。

① 水源涵养。采用水量平衡方程计算水源涵养量：

$ \mathrm{TQ}=\sum\limits_{i=1}^{i}\left(P_{i}-R_{i}-\mathrm{ET}_{i}\right) \times A_{i} \times 10^{3} $

(1)

式中：TQ为总水源涵养量(m³); P为降水量(mm); R为地表径流量(mm); ET为蒸散发(mm); A为类生态系统面积(km²); i为研究区第i类生态系统类型; j为研究区生态系统型数。

② 土壤保持。采用修正通用水土流失方程(RUSLE)的水土保持服务模型开展评价：

$ A_{c}=A_{\rho}-A_{r}=R \times K \times L \times S \times(1-C) $

(2)

式中：A_c为水土保持量(t/hm²·a); A_ρ为潜在土壤侵蚀量; A_r为实际土壤侵蚀量; R为降雨侵蚀力因子〔MJ·mm/(hm²·h·a)〕; K为土壤可蚀性因子〔t·hm²·h/(hm²·MJ·mm)〕; LS为地形因子; L表示坡长因子; S表示坡度因子; C为植被覆盖因子。

③ 生物多样性。生物多样性维护服务能力指数作为评估指标，计算公式为：

$ S_{\text {bio }}=\mathrm{NPP}_{\text {mean }} \times F_{\text {pre }} \times F_{\mathrm{tem}} \times\left(1-F_{\text {alt }}\right) $

(3)

式中：S_bio为生物多样性维护服务能力指数; NPP_mean为多年植被净初级生产力平均值; F_pre为多年平均降水量; F_tem为多年平均气温; F_alt为海拔因子。

④ 防风固沙。修正风蚀方程(RWEQ)是美国农业部应用于预测农田土壤风蚀量的模型^[17]，采用RWEQ方程充分考虑气候条件、植被状况、地表土壤粗糙度、土壤可蚀性、土壤结皮的情况，计算公式为：

$ Q_{\max }=109.8\left(W \cdot E \cdot S \cdot K^{\prime} \cdot C\right) $

(4)

$ s=150.71\left(W \cdot E \cdot S \cdot K^{\prime} \cdot C\right)^{-0.3711} $

(5)

$ S_{L}=\frac{2 X}{S^{2}} Q_{\max } \mathrm{e}^{-\left(\frac{X}{S}\right)^{2}} $

(6)

式中：Q_max表示最大沙尘土通量; W表示气象因子; E表示土壤可侵蚀因子; S表示土壤结皮因子; K表示土壤糙度因子; C表示植被因子; s为关键地块长度; S_L为土壤损失量; X表示距上风向距离, 取X=50 m。通过上述公式计算得到北京市生态服务重要性空间分布如图 2所示。

2.2.2 阻力面构建

研究参考俞孔坚^[18]构建区域生态安全格局思路，分别建立北京市生态源地扩张阻力(表 2)、城镇用地扩张阻力因子^[8](表 3)，采用熵值法并结合李红波等^[19]研究方法对各评价因子赋值，相对阻力值拟定在0~100之间，并基于层析分析法(AHP)确定评价指标权重，结合北京市生态现状和城市建设扩展水平，依据不同目标选择与生态源地、城镇用地扩张过程相关程度高的影响因子，通过GIS10.2栅格计算器对评价因子采用加权求和构建阻力面模型，利用空间分析技术求取生态源地扩展阻力面和城镇用地扩张阻力面(图 3)。

采用最小累积阻力(MCR)模型来建立阻力面^[20]；

$ \mathrm{MCR}=f_{\min } \sum\limits_{j=n}^{i=m}\left(D_{i j} \times R_{i}\right) $

(7)

式中：MCR为最小累积阻力值; f为MCR与变量(D_ij×R_i)之间的正函数; D_ij为从源j扩散至空间某点穿过生态表面i的空间距离; R_i为生态表面i对源扩散方向的阻力。

$ {\rm{MCR}}_{\text {差值 }}={\rm{MCR}}_{\text {生态源地 }}-{\rm{MCR}}_{\text {城镇用地 }} $

(8)

当MCR_差值 < 0时，表示区域生态源地扩张阻力小，适合扩张生态用地，MCR_差值>0时，表示区域生态源地扩张阻力大，适合城镇用地扩张；MCR_差值=0的区域为生态源地扩张与城镇用地扩张的临界区域，即生态城镇分界区(图 4)。西北部区域生态源地扩张阻力小，该区域适合生态源地扩张，东城区、西城区、朝阳区、丰台区和海淀区域生态源地扩张阻力较大，城镇建设、交通用地密度大。生态源地扩张阻力最大值的区域位于通州、大兴部分区域、这些区域适合城市建设。

2.2.3 生态廊道构建

生态廊道是连接生态源地的线状景观，为生物迁移扩散运动提供路径^[21]。当生态阻力差值较小时，有利于生态流扩张，研究利用模型计算各生物在阻力面上迁移所需要克服的最小累积阻力值，并识别最小成本路径，构建潜在生态廊道，选取MCR_差值面上相邻两个源地之间的阻力低谷，即累积阻力差值较大(负值)的路径构建区域生态廊道，并剔除冗余性较高的重复、交叉路径，保留相邻生态源地之间的最短路径。

生态节点可以给进行迁徙的物种提供休憩地，生态作用关键^[22]，提取生态廊道的交点、生态廊道和最小阻力路径的交点作为生态节点。交通道路会增大物种在生态源地间迁移的阻力^[23]，阻碍物种在生态源地间的迁移和交流。生态廊道与交通线路交汇区域形成生态断裂点，会阻碍生物信息传播和交流。研究选取高速公路、主干道、次干道路和其他道路与生态廊道的图层叠加，识别生态断裂点。

3 结果与分析 3.1 生态源地划定

如图 5所示，北京市生态源地总面积为3 568.95 km²，总共118块，占全区总面积的21.7%，其中北京市西北区域分布较多，中心城区中零星分布的小块生态源地主要是水域、公园和城市绿地，中心城区大量土地用来城市建设，导致生境斑块破碎化，生态源地的缺失严重影响了一些分布范围较广物种的繁殖。

3.2 生态安全格局构建

当生态阻力差值较小时，有利于生态流扩张，研究基于生态源地扩张阻力模型提取生态路径，并选择累积阻力差值较大(负值)和相邻生态源地之间的最短路径作为生态廊道，并采用最小累积阻力差值，利用分位数法将生态源地扩张成本划分为高水平生态区、中水平生态区、低水平生态区、生态城镇临界区、低水平城建区、中水平城建区、高水平城建区，从而构建北京市生态安全网络(图 6)。

“高水平生态区”主要分布在北京市西北部，属于生态系统功能高值区，大部分生态源地分布在“高水平生态区”和“中水平生态区”的合围区域。“低水平生态区”集中分布在高、中水平生态区域的延伸缓冲带，大兴区和房山区南部连片分布有“低水平生态区”。“生态城镇临界区”属于生态区域和城镇区域的缓冲带，集中分布在“低水平生态区”和“低水平城建区”之间，这些区域既是生态流缓慢扩散流动的区域，也是城市发展扩张的潜在区，应合理规划该区域的生态廊道布设，避免人工建筑阻碍生物物种交流、迁徙。“中、高水平城建区”集中分布在北京市中心城区，属于城市建设高度发达区域，这些区域中零星镶嵌着“低水平生态区”斑块，这些生态源地斑块主要是城市公园、城市绿地和湿地。

生态廊道把相互独立的生态源地在空间上连接起来，是保持生态流、生态功能连通的关键生态组分构成^[24]。本研究计算得到188条生态廊道，形成了北京市生态网络空间分布(图 6)，西北部生态涵养区生态廊道数量较多，中心城区生态廊道数量较少。考虑生态廊道连续性和可落地实施要求，筛选同时服务多个生态源地与节点的生态廊道作为重要生态廊道(11条)，其他生态廊道作为一般生态廊道。生态节点是生物迁徙和交流的关键生态战略点^[25]，提取生态廊道交汇点，即关键生态节点建立“踏脚石”斑块，生态节点共有153个(图 6)，主要分布在研究区西北部，这些生态节点可以加强区域生态安全与保护。比较北京市16个区廊道网络密度(图 7)；石景山区生态网络密度最大，其次是昌平区、怀柔区和密云区，东城区生态网络密度最小。

北京城市交通发达，路网密度较大，将生态廊道与研究区路网叠加，提取生态断裂点720个。其中，一级断裂点为生态廊道与高速公路的交汇点，共计55个；二级断裂点为生态廊道与主干道的交汇点，共计224个；三级断裂点为生态廊道与次干道的交汇点，共计198个；四级断裂点为生态廊道与其他道路的交汇点，共计243个；生态断裂点集中分布中心城区(东西城区、朝阳区、海淀区、丰台区)，这些区域路网聚集，人口数密集，生态廊道断裂程度较高，其次生态断裂点较多的区域分布在昌平区、顺义和大兴区，这些区域的生态廊道也不利于生物物种迁徙扩散。针对断裂程度较高的生态廊道，不建议规划生态廊道，针对现有的珍稀濒危物种可采取就地、迁地保护措施。

4 讨论与结论 4.1 结论

本文综合考虑城市生态安全格局空间完整性和生态系统服务功能，选取北京市63处自然保护地和7处城市公园，叠加4类综合生态系统服务等级最高的区域，共提取118块生态源地和11条重要生态廊道，利用最小阻力模型差值模型(MCR_差值)将北京市划分为7大区域(高水平生态区、中水平生态区、低水平生态区、生态城镇临界区、低水平城镇区、中水平城镇区、高水平城镇区)，构建北京市生态安全格局。

(1) 生态源地集中分布在北京市西北部生态涵养区，呈包围中西城区态势，在空间上连片集中，位于“城建区”的小斑块生态源地主要是城市公园、城市绿地和水域，“城建区”生态源地破碎化程度较高，呈离散性分布。

(2) “高水平生态区”与生态源地分布区域基本叠合，总面积占研究区12.01%，属于生态核心区域；“中水平生态区”面积占研究区31.3%，主要分布在生态核心区外围，属于生态缓冲区域；“低水平生态区”面积占研究区25.63%，主要分布在“中水平生态区”外围，属于生态过渡区域；“生态城镇临界区”面积占研究区5%，分布较零碎，该区域对生态流扩散形成阻隔，因此该区域是生态化改造的关键区域。

(3) 根据生态廊道分析，由于城区密集路网和高密度人类活动影响，城区的生态廊道断裂程度较高，生态保护安全性较低，生态源地破碎程度较高，很难建立完整连通的生态廊道。西北生态涵养区域的“中、低水平建城区”呈蔓延发展趋势，占用生态绿地较多，中心城区呈组团连片发展趋势。

(4) 未来应优化生态节点空间布局，对生态环境良好且具有一定生态功能的“踏脚石”进行筛选和保护，并采用科学合理的空间格局设计，即以保护生态源地、规划生态廊道，控制生态城镇临界区建设等方式获得生态效益。

4.2 讨论

构建城市生态安全格局对于缓解城市建设和生态保护的矛盾具有积极意义，尤其是生态空间受到严重挤压的超大人口城市。本文选取北京市为研究区域，将生态资源保护与城市合理发展相协调，充分考虑人类活动和生态系统之间各组分的关联性，建立了“生态评价—生态网络构建—生态安全格局”的生态安全格局方法。国内外相关生态安全格局研究中，主要关注生态源地、缓冲区、生态廊道和生态节点的组成^[26]。超大城市的结构和功能较为复杂，城镇建设用地已经趋向饱和，生态城镇相交互的临界区应是城市生态安全关注的重点，以往的生态安全格局相关研究很少关注生态区和城镇区的过渡区域^[27-28]，本研究采用最小累积阻力模型(MCR)生成生态源地扩张阻力面与城镇用地扩张阻力面的差值面，并利用分位数法将生态源地差值成本面划分成7个区域：高水平生态区、中水平生态区、低水平生态区、生态城镇临界区、低水平城建区、中水平城建区、高水平城建区，其中生态城镇临界区分布较为零散且广泛(图 6)，是北京市生态安全格局构建需要重点关注的区域。

北京市生态格局破碎化现象明显。北京市中心城区由于建设用地较多，对生态资源造成割裂，导致生态廊道破碎化严重，集中连片的生态资源极少，生态源地破碎化程度较高，生态系统服务功能较低；西北部分地区自然条件良好，生态源地和生态节点集中分布在西北部生态涵养区，但生态涵养区生态廊道受道路基础设施阻隔，造成生态断裂点，这对于物种迁徙安全形成威胁，在建设基础设施时需要充分考虑生态廊道的保护需求，采取相关工程措施尽可能降低道路设施对生态网络的影响^[28]；部分生态廊道冗余性较高，重要生态廊道有待进一步拓宽，以保障生态要素连通交流^[29]。

近10 a来，北京市人口压力剧增，生态安全面临严重威胁，对比2008年俞孔坚^[30]对北京市生态安全格局城市增长预测，北京城市建设进一步扩张，建城区呈“摊大饼”式蔓延，中心城区绿地面积减少，中心城区与周边区域组团连片发展，城区生态空间破碎化程度较高，现亟需控制“低水平生态区”和“生态城镇临界区”城建水平，防止西北部屏障进一步割裂。研究基于生态安全格局构建范式，从维护生态系统完整性视角出发，识别生态源地、构建生态网络、分级国土空间制定生态保护措施，尤其针对城市发展空间有限的超大人口城市，划分“生态城镇临界区”，这对于阻碍城镇化无序扩张，缓解资源环境压力有重要现实意义。研究可为科学规划国土空间安全格局、实现资源优化配置提供参考，下一步研究中，应确定单一因子的城市生态安全格局构建以及边界的确定方法。