随着国家经济快速发展, 城市生态风险问题呈多发高发态势, 其风险形势更加复杂和难以预测, 给人民生命财产安全和生态环境质量造成了巨大损害^[1]。2020年全球新冠疫情的爆发, 从短期来看对医疗相关的生态环境保护(医疗废弃物等)提出额外的需求, 从长期来看后疫情时代生产生活和社会经济的逐渐恢复会加剧资源供需的矛盾, 引发新的城市生态风险。因此, 如何提高城市应对生态风险的能力成为实现城市可持续发展的一个重要方向。2020年11月3日, 《中共中央关于制定国民经济和社会发展第十四个五年规划和二○○三五年远景目标的建议》提出“推进以人为核心的新型城镇化。强化历史文化保护、塑造城市风貌, 加强城镇老旧小区改造和社区建设, 增强城市防洪排涝能力, 建设海绵城市、韧性城市。提高城市治理水平, 加强特大城市治理中的风险防控。”, 更表明了建设韧性城市对我国经济社会的重要作用。

韧性这一概念最早源于生态系统学, 指的是系统恢复平衡的速度, 应对危机并恢复的能力, 适应新环境的能力^[2]。城市研究的学者们后来借鉴这一概念, 对城市韧性进行了研究。尽管学界对城市韧性的理解还未形成统一的定义, 但大体都认为“城市韧性”指城市对灾害和突发性事件的反应、适应和应对能力, 以减轻损失, 并从中快速恢复^[3-5], 其内涵包含了城市生态韧性、社会韧性、经济韧性、制度韧性和基础设施韧性等多个维度^[6]。近年来, 纽约、伦敦、芝加哥、鹿特丹等城市纷纷提出最具代表性的韧性城市规划, 这些韧性城市规划主要针对不同的气候风险, 强调城市对未来气候风险的综合防护能力^[7-8]。

作为城市韧性的重要维度, 生态韧性的评估也得到了学者们的关注, 但大多数的关注点均从生态风险和环境风险的防控角度出发^[9-12]。在评估思路方面, 毕军和王华东提出要对“自然-社会-经济”复合系统进行环境风险评估^[13]；张小飞等从自然灾害、环境污染及生态退化3个方面, 构建了城市综合生态风险评价的分析框架^[14]; 王祥荣等围绕长三角城市群重点产业、城市灾害生态风险和典型生态系统, 构建集地表水、大气、土壤、生物多样性于一体, 常规监测与应急预警相结合的长三角城市群生态安全评估与风险预测预警技术体系^[15]。但是, 很少有研究从韧性的特征角度考虑城市生态韧性水平的空间变化, 也忽略了生态韧性和风险的分区管理评估, 导致其成果难以真正落实于城市规划应用。

基于以上考虑, 本研究首先从韧性的基本特征(抵抗力、适应力和恢复力)3个方面构建了城市生态韧性评估模型。其次, 以杭州为例, 评估了1995—2015年城市生态韧性时空变化格局；并利用空间自相关模型对2015年杭州生态韧性分区管理模式进行了探究；最后模拟了2035年基准和创新创业导向的两种情景下的杭州市城市扩张情况, 并评价了不同情景下的城市抵抗力、适应力、恢复力, 以及生态韧性水平。基于以上研究结果, 以期为生态韧性城市建设以及产业空间布局提出一些可供操作的政策和规划建议。

1 研究区和数据来源 1.1 研究区

杭州作为浙江省省会, 位于浙江省的北部, 属亚热带季风性气候, 位于长三角腹地是中国东南地区最发达的城市之一。杭州市西部属浙西丘陵区, 东部属于浙北平原, 全市总体地势格局为西南高、东北低, 自西南向东北倾斜, 呈阶梯下降。

改革开放以来, 杭州市经济发展迅速, 人均GDP从1995年的4614元上升到2015年的112600元, 在中国35个主要城市中排名第八。相关文献表明, 近20年来杭州市空间外延式增长迅猛^[16-18], 大量自然生态用地被占用, 引起很多环境问题, 包括城区生境质量下降、生态系统服务功能降低、城区温室效应剧增等等。此外, 杭州在谋划“十四五”发展目标时, 提出要以未雨绸缪的布局防范化解各类风险挑战, 不断增强城市发展韧性。因此, 以杭州为案例来研究其城市生态韧性时空变化具有重大意义。

1.2 数据来源

本研究利用1995、2005与2015年的遥感影像数据通过监督分类结合目视解译获得分类结果。每幅遥感图像都经过图像裁剪、大气校正、几何校正的预处理。首选通过监督分类得到林地、耕地、水域与湿地、建设用地这4个类别, 精度检验结果表示分类的Kappa系数均达到85%以上。此外, 杭州市重点排污单位名录来自杭州市生态环境局, 其坐标信息来自其企业名称与POI兴趣点的匹配。

2 核算方法与模型构建 2.1 城市生态韧性评估模型构建

国际区域可持续发展协会(ICLEI, 2002)定义韧性城市为对于危害能够及时抵御、吸收、快速适应并做出有效反应的城市, 重点关注如何增强城市对抗冲击的能力, 提升自身韧性强度, 提高自身恢复能力, 由此可以总结韧性城市包括抵抗力、恢复力、适应力三大特征。相对于风险的事前预警特性, 韧性可用于事前评估或事后回溯分析, 因此韧性城市若干风险防控管理可以达到“事半功倍”的功效。因此, 本研究从抵抗力、适应力和恢复力三个方面构建城市生态韧性评估模型。

抵抗力表示城市生态系统抵抗外力干扰的能力。很多研究指出, 生态系统抵抗力与其生态系统服务功能密切相关^[19-21]。生态系统服务功能是指生态系统与生态过程所形成的维持人类赖以生存的自然环境条件与功能^[22-24]。Costanza等将全球生态系统类型、服务功能分成了17个类别^[25]。本研究从水、土、气等多介质出发选择气候调节、气体调节、水源涵养、土壤的形成与保护、废物处理等调节型生态系统价值构成了生态系统的抵抗能力(P)。

本研究利用生态系统服务价值的当量系数, 进而获得生态系统服务价值系数。生态系统服务价值当量因子是生态系统潜在服务价值的相对贡献率, 等于每年每公顷粮食价值的1/7 ^[26]。杭州的生态系统服务价值的计算公式如下：

(1)

式中, ESV为生态系统服务价值；A_k是土地利用类型k的面积；VC_fk是土地利用类型k的第f项生态系统价值系数, 其系数来自于长三角地区土地利用类型的部分生态系统服务价值系数^[27]。一个越稳定的生态系统, 其适应力(A)越高^[28]。因此, 本研究利用生态系统景观结构稳定性相关指标来表述其适应力。生态系统的景观结构稳定性取决于与空间异质性和景观连通性有关的景观格局组织稳定性^[29-31]。本研究利用景观指数从景观异质性和景观连通性来评价景观稳定性, 景观指数可以反映景观生态空间格局, 已经广泛应用于景观生态区划、生态规划、生态监测、生态预测、生态影响评价等, 可以对景观生态空间做出有目的和适时的调控^[32-33]。景观异质性和景观连通性分别描述生态系统景观结构的两个方面, 彼此不能替代, 因此其权重可假设为相等^[34], 具体指标如表 1。

生态系统恢复力也可以称为生态弹性, 体现生态系统遭受危害而恢复原样的能力与潜力^[35]。当遇到外部压力(例如气候或地质灾害)时, 接近自然生态系统的土地利用类型更容易恢复, 而人为主导的建设用地类型在面对外部压力时恢复能力较低, 因此遭受的破坏更大。其测算公式参考根据Peng等人提出的生态弹性模型与系数^[34], 具体如下:

(2)

式中, R是生态系统恢复力, A_k是土地利用类型, 土地利用类型k的生态恢复力系数, 系数参考Peng等人的研究^[34]。

此外, 抵抗力、适应力和恢复力的计算具有不同的单位, 因此在使用多个指标进行评估时需要数据标准化。本研究将这三者标准化到[0, 1]这个区间, 计算韧性水平的具体公式如下：

(3)

式中, A为适应力, P为抵抗力, R为恢复力。

2.2 空间自相关模型

空间自相关分析是用于衡量空间变量的分布是否具有集聚性, 是分析空间格局的重要方法和有效手段。生态韧性与自然地理要素分布、社会经济发展状况相关, 这些因素在空间上均具有随机性和结构性, 具有地学特征, 存在着一定的空间关联性, 可以运用空间自相关进行分析。此外, 空间自相关模型在生态系统分区管理方面也已经进行了大量应用^[36-39]。空间自相关分析可分为全局空间自相关分析和局部空间自相关分析, 一般利用Moran′s I指数来表示。本研究采用采用局部空间自相关的LISA来表现韧性水平在局部空间集聚规律^[40]。

2.3 城市土地利用变化模型

FLUS模型是基于GeoSOS软件开发, 主要结合人工神经网络(ANN)算法、自适应惯性与竞争机制对传统CA模型进行了改进, 更加适用于多地类情景模拟^[41]。FLUS模型主要包括两个方面：人工神经网络(ANN)和自适应惯性竞争元胞自动机。人工神经网络(ANN)是通过训练样本来估计每个单元网格上的每种土地利用类型的可能性, 而自适应惯性和竞争机制是用来处理土地利用类型之间的相互作用。FLUS模型已经在城市土地利用扩张模拟方面进行了广泛应用^[42-45]。

本研究设置二种情景来模拟2035年杭州市城市扩张情况, 这二种情景分别为基准情景以及创新创业导向的“创新天堂”情景, 并利用GeoSOS-FLUS软件实现FLUS模型的土地利用格局模拟。

基准情景根据2005年、2015年两期的土地利用变化情况, 采取原始土地利用转移矩阵预测得到的土地利用规模, 综合其他研究对于土地利用变化驱动因子的分析及探究^[46-47], 以及规划实践过程中实际影响杭州城市发展的因素, 最终选择从社会经济、自然条件、交通条件、规划因素等方面选择驱动因子(表 2)。

“创新天堂”情景这个概念来自于2017年《杭州市总体城市设计修编(草案)》所提出的“风雅钱塘, 诗画江南, 创新天堂”的杭州未来城市定位。“创新天堂”情景是一种创新产业发展视角且耕地保护比较严格, 与基准情景相比, 该情景设置耕地转为建设用地的比例下降15%。此外, “创新天堂”情景增加了创新驱动因素(表 2), 根据《杭州市智能制造产业发展“十三五”规划》与《杭州城东智造大走廊发展规划纲要》提取相关的规划制造业创新中心与主要产业平台, 计算每个栅格到创新中心与主要产业平台的距离并赋值, 得到杭州创新中心与主要产业平台驱动因子。最后, 本研究将以上创新因素驱动图层作为为创新产业发展视角下额外因子, 输入至人工神经网络, 计算适宜性概率图层, 进行后续创新产业发展视角下的用地空间格局的预测模拟。

3 结果与讨论 3.1 生态韧性时空变化

图 1表示了1995—2015年期间, 每间隔10年的杭州城市生态韧性水平的时空变化。如图 1所示, 从时间上来看, 城市生态韧性水平呈现下降趋势, 从网格水平上统计其均值下降了32.86%, 最高值下降了20.27%, 标准差也下降了34.33%。

空间分布来看, 1995年杭州生态韧性水平区域差异不大, 大部分区域处于中等和高值生态韧性水平, 低值生态韧性水平区域只要分布在城市东北区域和中心区域偏西北、东北部分。随着城市扩张, 2005年生态韧性低值区域向东北和东南方向扩张, 东北和东南区域原本连续分布的高值区域变为零散分布。由于缺乏相对适宜的整体规划和统筹安排, 位于东南和东北部分的萧山区乡镇企业发展迅速, 这些工厂多呈小规模散落分布, 占用了大量的自然资源, 生态服务功能降低, 景观破碎度提高, 导致其韧性水平大幅度下降。此外, 城市东北部靠近杭州湾, 原本存在大量优质湿地资源, 其东北部围垦过度大大破坏了自然生态景观。2015年杭州城市总体生态韧性水平不高, 大部分区域处于低值水平, 高值区域零星分布, 而中等区域仅在城市东南分布, 西北部和西部原本存在的中等韧性区域也转为低值水平。快速土地城镇化占用了大量的耕地和林地, 大量农村务工人员涌入城市, 致使住房用地需求增加, 部分耕地被占, 转变为建设用地, 又进一步加剧了耕地资源不足的严峻挑战。

杭州1995—2015年生态韧性水平时空变化与大连^[49]、京津冀^[50]等地的相关研究呈现出一定相似：随着城镇化发展, 生态韧性水平总体分布规律为城市中心城区周边的逐渐降低, 韧性水平低值区域从面状零散分布变化至面状连片分布。尽管杭州与大连、京津冀等地交通、经济等区位相差较大, 但是我国城市早期“摊大饼”式的城市扩张模式都趋同。以工业园区建设为主导的粗放的城市扩张模式导致城市周边大量优质耕地被占用, 而后期的“耕地占补平衡”政策虽然在表面上保证了18亿亩耕地红线, 但是在地方实际操作中往往出现“占优补劣”等现象, 导致林地、湿地等大量生态资源被侵占作为重要耕地补充来源, 进一步导致中心城区周边生态用地景观破碎, 降低了其生态韧性水平。

本研究利用ArcGIS 10.2对杭州市重点排污单位进行了核密度分析, 并利用Geoda分别对韧性水平和韧性水平—排污企业密度进行了单变量和双变量空间自相关分析, 其结果如图 2所示。图 2显示, 重点排污企业主要集中在城市北部和东北部, 分别靠近余杭临平工业园和萧山大江东产业集聚区, 这两个区域地势平坦、产业优惠政策丰富且交通便利, 承接了杭州主城区高污染产业的转移：杭州市从2002年起陆续确定杭州中策、杭氧集团、油脂化工、杭州卷烟厂、杭州电缆有限公司、杭州汽车发动机厂、杭华油墨、杭州重型机械有限公司等数家企业的整体搬迁计划, 这些企业大多将迁入配套齐全和产业集聚的江东、临江工业园区和余杭、钱江经济开发区。此外, 韧性水平的空间自相关分析结果(图 2)表明生态韧性低值区域在城市中部和西北部出现集聚, 而高值区域在城市东南、东北以及中心区域偏西南部分出现集聚, 近年来, 政府投入大量资金以提高靠近杭州湾的城市东北部分和城市中心区域的西湖风景区的湿地质量, 湿地生态资源可以通过发挥其生态屏障作用来阻断不利因素和灾害的影响, 对于提高生态韧、防范化解生态风险具有重大作用。城市生态韧性水平—排污企业密度的双变量空间自相关结果表明(图 2), 城市的东北区域呈现低韧性水平—高排污企业密度的集聚分布, 这与杭州整体向东发展的城市策略密切相关。杭州东北部地势平坦, 且紧邻杭州湾, 与宁波、上海等长三角核心城市交通便利, 形成大江东产业集聚区, 同时承接了周边城市的以印染化工为主的传统制造业的转移。产业集聚区的发展带动东北区域建设用地的扩张, 侵占周边生态用地资源, 使得生态用地斑块破碎化, 导致该区域生态韧性水平降低。同时, 该区域生态韧性水平比较低, 而高排污企业密度高, 说明其环境污染风险比较大, 但是其防范化解能力比较低。因此需要加大该区域环境监测和生态治理的资金投入, 严防重大环境污染事件的发生。

3.2 城市扩张背景下的生态韧性情景模拟

本研究基于2005、2015年两期杭州城市土地利用数据模拟2035年其空间格局, 这10年期间土地利用转移情况如表 3所示。2005—2015年期间林地、建设用地规模比较稳定, 很难转化为其他用地, 而耕地、水域比较容易向其他地类进行转换。耕地是建设用地的主要来源, 占全部向建设用地转移的70.49%。此外, 水域向耕地转移占全部向耕地转移的68.52%, 这与杭州东北区域湿地复垦密切相关。

图 3是基于FLUS模型模拟的2035年杭州市城市土地利用空间格局结果。图 3显示, 基准情景的2035年建设用地更多地以边缘式发生在主城区的近郊, 且并未呈现明显集聚, 是一种类似摊大饼城市扩张模式。“创新天堂”情景是一种紧凑增长模式, 建设用地主要增长在城市西北和东南部。这些部分是未来科技城、云栖小镇、大江东产业集聚区、杭州经济开发区、萧山经济开发区、余杭经济开发区(钱江经济开发区)、萧山临空经济示范区、杭州仁和先进制造业基地所在的区域, 说明“创新天堂”情景更好的契合了创新产业发展规划的产业发展需求。

图 4和图 5分别是创新创业导向的“创新天堂”和基准情景的杭州生态韧性空间分布水平。相比基准情景, “创新天堂”城市发展情景下杭州整体生态韧性水平更高。从抵抗力、适应力和恢复力3个方面来看, 两者的恢复力水平差异不是很大, 而“创新天堂”情景的抵抗力水平和适应力水平明显高于基准情景。特别是适应力水平, 基准情景下区域分异尤为明显, 展现由四周向中心梯次递减的趋势, 两极化现象较为明显, 在城市中心区域出现明显低值集聚现象, 而“创新天堂”情景的适应力水平的空间差异不大, 没有出现明显的低值集聚现象。这与两种情景的城市扩张方式相关, “创新天堂”情景是一个更为紧凑填充式的城市增长模式, 而基准情景的城市增加以边缘式发生在建成区的近郊。虽然边缘式增加可以提高景观多样性, 提高景观异质性, 但是会提高其破碎程度并降低景观连通性。本研究的适应力指标体系设计中(表 2), 景观破碎度权重(0.5)大于景观多样性(0.25), 这是因为考虑到快速城镇化的后期, 各种景观类型相互渗入以及建设用地小斑块, 逐渐变成以建设用地为主导的景观格局, 其景观多样性会降低, 而这种稳定状态对适应力提高并没有正面作用。

4 结论与建议

本研究从抵抗力、适应力和恢复力3个方面构建了城市生态韧性评估模型, 并以杭州为例, 评估了1995—2015年城市生态韧性时空变化格局, 再者利用空间自相关模型对2015年杭州生态韧性分区治理模式进行了探究。最后, 基于FLUS模型模拟了2035年基准和创新创业导向的两种情景下的城市土地利用空间格局, 并评估了不同城市发展情景下的生态韧性水平空间分布。基于以上研究, 本文得到以下结论：

(1) 总体来看, 杭州1995—2015年城市生态韧性水平呈现下降趋势, 从网格水平上统计其均值下降了32.86%, 最高值下降了20.27%, 标准差也下降了34.33%, 说明更多网格呈现低值集聚。从空间分布来看, 1995—2015年期间, 全域城市生态韧性水平低值区域呈零散分布变化至面状连片分布：1995—2005年, 生态韧性低值区域向东北和东南方向扩张, 这与东北部湿地围垦过度以及被建设用地侵占密切相关。2005—2015年, 杭州西北部和西部原本存在的生态韧性中等水平区域也转为低等水平, 该区域大量优质林地生态资源的被侵占是其变化的重要原因。

(2) 从空间分布来看, 2015年生态韧性低值区域在城市中部和西北部出现集聚, 而高值区域集聚在城市东南、东北以及中心区域偏西南部分。同时, 城市生态韧性水平—排污企业密度的双变量空间自相关结果表明杭州的东北区域显著呈现低韧性水平—高排污企业密度的集聚分布, 这说明此区域环境污染风险比较大, 依靠生态系统自身韧性来防范化解风险能力比较低。

(3) 相比基准情景, 紧凑填充式的城市增长的“创新天堂”城市发展情景下2035年杭州整体生态韧性水平更高。尽管两者的恢复力水平差异不是很大, 但是“创新天堂”情景的抵抗力水平和适应力水平明显高于基准情景。

根据以上研究结果, 本文为杭州生态韧性城市建设提出以下三点建议：

(1) 促进市域均衡开发, 严格控制湿地开发, 加强生态用地质量与数量建设, 提升城市整体韧性：杭州应当在全域范围内统筹发展资源, 避免在中心城区形成严重的景观压力和城市蔓延问题, 降低过度集中导致的生态韧性降低, 同时也使发展机会得以分享。此外, 杭州紧接杭州湾, 拥有大量湿地资源, 一方面要尽量避免城市东北部的侵占湿地的围垦活动, 加强湿地景观的重塑；另一方面尽量避免城市建设用地对城市中心城区部湿地资源的进一步侵占和切割, 保护生态用地, 重点在于保持西溪湿地、西湖、以及东北部杭州湾湿地的完整性, 发挥其生态屏障作用, 阻断不利因素和灾害的影响。

(2) 实施城市环境生态治理分区管理模式：生态韧性低值集聚的中部和西北区域, 应对环境生态风险的能力较低, 未来国土空间规划中要通过加强该区域生态空间保护与建设来消除差异化城市生态韧性空间格局；东北区域呈现显著低韧性水平—高排污企业密度集聚分布的部分亟需尽快搬移高污染排放企业, 加速当地产业升级, 同时加强该区域的环境生态质量监测水平, 以此降低局部生态风险, 防止局部风险扩散引起全域连锁风险效应。

(3) 优化新区发展模式, 鼓励紧凑的城市发展：城市东部是杭州未来重要的发展空间, 要防止形成“摊大饼”式的开发局面, 并大力建设绿色景观。此外, 建议杭州政府采用产业升级和转型策略来缓解中心城区的人口和土地压力, 并投入更多的资金用于常态下环境治理与生态保护, 加大技术与人才流通, 努力消除未来城市增长带来的城市内部的生态韧性差异化格局。在未来城市向东部发展的过程中, 建议杭州政府逐步探索城市内部的区县政府之间以生态风险转移为主体的生态补偿制度, 缓和城市内部区域生态韧性差异化格局, 提高城市整体防范化解生态风险的能力。