城市作为“社会-经济-自然”复合生态系统, 其景观具有高动态度的特征^[1-3], 具体表现为土地覆盖/土地利用在时间和空间两个维度上的高频变化^[4-5]。这种变化可能会带来环境质量的下降, 如产生扬尘、噪声等, 并降低城市生态系统的稳定性, 影响城市生态系统服务, 从而导致城市可持续发展能力下降^[6-7]。此外, 城市内部高动态地区通常是城市生态系统的敏感区, 也是城市监管的重点区域^[8]。因此, 开展城市内部景观格局动态变化的定量化研究是理解城市生态系统演变规律以及城市监管和城市规划的基础^[9]。

城市化与城市发展对景观格局的影响, 既包括城市扩张过程中对景观格局的影响, 也包括城市内部更新产生的景观格局的动态变化。基于Landsat等中等分辨率遥感影像众多的学者在城市、区域, 甚至全球尺度上研究了城市扩张引起的城市景观格局动态, 揭示城市发展过程中建设用地的增长及其带来的农田、森林、草地和湿地等其他生态系统丧失的城市化影响^[10-13]。研究中通常假设城市内部土地变化都是从透水面转化为不透水面, 且这种变化不可逆^[14-15]。但是, 城市发展过程中, 不仅存在城市空间的向外扩张, 城市建成区也存在高度的动态变化^[5], 只是这种内部动态受限于空间分辨率, 基于中等分辨率的遥感影像的方法难以揭示^{[9, 16]}。城市内部景观的高度异质性, 需要采用更高空间分辨率的数据, 在更精细的尺度, 才能准确刻画^[16-17]。街区作为城市管理的基本单元, 也是开展城市更新、内部填充等土地开发活动的基本单元^[18-19]。国际上, 如美国等国家以街区尺度开展的研究相对较多, 而中国由于城市街区尚未有明确的官方定义, 街区的划分方法还有待统一^[18], 使得在街区尺度上对城市内部动态度开展定量化的研究还十分罕见。

为揭示城市内部土地覆盖动态度变化的空间异质性, 本研究以深圳市为例, 利用多等级道路数据划分街区, 基于高分辨率卫星影像解译的土地覆盖数据对街区尺度的土地覆盖动态度及其土地利用强度进行量化, 以期回答如下3个科学问题：(1)深圳城市建成区街区尺度上土地覆盖变化的剧烈程度如何, 其在空间上如何分布？(2)城市内部街区不同土地覆盖类别动态度及其土地利用强度有何差异？(3)街区的动态度与街区大小是否存在关联？回答上述问题, 有助于深化对城市景观动态的认识, 为城市景观动态的监测、管理以及城市规划提供科学支撑。

1 研究区

深圳市地处广东省南部, 毗邻中国香港, 位于东经113°46′—114°37′, 北纬22°27′—22°52′之间, 全市面积约为1997 km², 是中国实施改革开放政策以来设立的第一个经济特区(图 1)。40余年来, 伴随着城市的飞速发展, 深圳市土地覆盖也发生着高度动态变化。1979年到2017年, 深圳市建设用地从26.49 km²增长到946.14 km², 经历了34倍的城市增长^[20]。2010年, 深圳市颁布了《深圳市城市总体规划(2010—2020)》, 控制外部增量, 挖潜内部存量。2013年深圳市的城市更新项目达72个, 占地面积5.45 km², 占深圳市建设用地面积的0.63%^[21]。深圳城市空间增长开始逐渐由“摊大饼式”的快速扩张向“存量优化”的城市更新模式转变^[22]。因此, 深圳市是研究城市内部动态度量化的理想地点。

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

研究基于OpenStreetMap(https://www.openstreetmap.org)下载的2017年5等级道路矢量作为划分深圳市街区的边界, 不同等级道路分别为：高速公路(motorway)、主干道(trunk)、一级道路(primary)、二级道路(secondary)、和三级道路(tertiary)(表 1)。为了刻画深圳市高时空动态度特征, 选取2017、2018年SPOT 6和GF-1号(全色波段分辨率为1.5 m, 多光谱分辨率为6 m)高分辨率遥感影像解译的土地覆盖结果(图 2), 进行街区动态度的量化。遥感影像土地覆盖解译结果的类别共有9类, 分别为林地、草地、水体、建筑、在建用地、道路、其他不透水地表、裸土、农田。参考Qian等^[16]精度验证的方法, 获得两期影像分类总体精度分别为89.67%和90.67%。

2.2 研究方法

本研究在街区尺度上量化深圳市动态度, 主要包括以下三个方面：(1)基于深圳市多等级道路矢量划分深圳市街区；(2)利用2017、2018年土地覆盖分类结果量化街区的动态度及其土地利用强度；(3)统计分析不同街区动态度与街区大小以及土地覆盖类型的关系。

2.2.1 城市街区的划分

研究中通常将街区定义为具有正方形或长方形形状, 建筑填充其中, 以闭合街道为边界的密集空间^[23]。作为城市居民社会活动的基本单元, 街区已成为城市生态、城市规划、城市管理等领域的重要分析尺度。近些年, 国内已有学者在考虑空间、功能维度等方面的基础上提出街区的划分方法^{[18, 24-25]}, 但不同方法仍存在复杂、不易应用, 以及与实际不符等情况。本研究借鉴Xia等的方法^[24], 利用5个等级的道路矢量来划分城市街区。具体方法如下：在ArcGIS10.1软件中，将获取的不同等级的道路矢量数据进行融合(Merge)处理, 手动修改不同等级道路间的拓扑错误, 利用线转面(Feature To Polygon) 工具将道路线状矢量转化为街区面状矢量。由于不同等级道路间相互交错, 产生许多较小面积的地块, 这些地块并不具有城市生态学和城市管理学意义, 对初步获得的街区进行如下处理：(1)将面积小于1hm²的地块合并到相邻街区地块中；(2)2005年起深圳市划定了生态控制线, 管制绿色空间^[22], 部分由道路直接产生的地块内部全部为林地, 不符合街区定义, 故将其删掉, 最终得到深圳市街区分析单元。

2.2.2 城市街区动态度及其土地利用强度量化

为反映街区内土地覆盖变化的频率与幅度, 研究使用街区动态度指数对其进行量化。城市内部土地覆盖变化包括地类转出和地类转入两个过程, 只关注地类的净变化无法反映土地覆盖真实的动态度^{[16, 26]}。动态度的量化需要同时关注这两个过程。本研究借鉴Wang等量化城市绿地动态度的方法, 利用土地利用动态度指数对深圳市街区尺度上的土地覆盖动态度进行量化^[9]。首先, 基于深圳市2017年和2018年的土地覆盖分类结果, 进行两年间土地覆盖的变化检测分析。然后, 在街区尺度上对深圳市9种地类间的转化结果进行动态度指数(Dynamic Index, DI)的计算：

(1)

其中，DI为街区内土地覆盖动态变化的比例；s为遥感影像解译的期数，本研究中s为2；i为土地覆盖类别代码；j为土地覆盖类别数量，本研究中j为9；P(i→not i)为土地覆盖类别i转变为非i类别的比例。

土地利用强度指数不仅可以反映土地利用程度的高低，也可以反映土地利用开发的环境效应，被广泛应用于土地覆盖动态变化研究中^[27-28]。研究使用刘纪远等提出的土地利用强度综合指数^[29]，其计算公式为：

(2)

土地利用强度变化可定量反映土地利用的综合水平及其变化趋势，计算公式如下：

(3)

其中，I为土地利用强度综合指数，A_i为第i级的土地利用程度分级指数; C_i为第i级土地利用程度分级面积百分比；ΔI为土地利用强度综合指数的变化，该值为正数时，表明土地利用程度增加，而该值为负数时，表明土地利用程度降低。C_ia和C_ib分别为第i级土地利用程度在时间a、b时的分级面积百分比。

根据庆大方等对土地利用分级的描述^[29]，以及深圳市土地利用的实际情况，对深圳市土地利用强度等级进行以下分级：

2.2.3 统计分析

热点分析(hotspot analysis)是一种识别感兴趣变量聚集分布区域的空间分析方法，广泛应用于土地覆盖动态及生态系统服务研究^[30-31]。本研究利用局部Gi^*统计量(Getis-Ord Gi^*)对街区尺度动态度进行热点分析，识别其热点和冷点的空间分布，分析其空间聚集特征，可为城市内部高动态度区域的监管提供数据支撑。为探究街区动态度在不同街区大小上的差异，研究基于街区动态度与街区面积的散点图及二者数据分布的箱式图，进行回归分析，从而为城市街区动态监测提供参考。土地覆盖转变过程中，不同地类对动态度的贡献度往往不同，将每种类别转移占总变化的比例定义为动态度贡献度(公式4)，统计每个街区中不同类别转移对于街区动态度贡献度，分析街区不同土地覆盖转移类别的差异。

(4)

其中，P(i→j)为街区内部类别i转化为j的比例。

3 结果 3.1 街区动态度空间分布与热点分析

深圳市6个等级道路共划分出5044个地块, 去除掉位于生态控制线内的地块, 共得到4796个街区(图 3)。街区面积共1061.43 km², 占深圳市国土面积的53.15%。计算街区动态度指数, 发现2017—2018年深圳市共有716个街区发生土地覆盖变化, 占街区总数的14.93%, 表明深圳市街区尺度上的土地覆盖变化仍处于高度动态中。其中, 动态度较高的街区主要分布在盐田区、大鹏新区、南山区、宝安区, 街区平均动态度分别为18%、17.6%、17.5%和17.1%, 而龙岗区、福田区、坪山新区等区的土地覆盖变化动态度则相对较低, 街区平均动态度分别为10%、12.2%和12.7%(图 4), 说明深圳市西部沿海地区总体上仍处于较为剧烈的开发中。

2017—2018年深圳市在不同置信度下动态度的热点分布如图 5。本文以P < 0.05(95%可信度)作为判断显著性的标准, 动态度热点主要分布在宝安区西北部、宝安区与南山区交界的前海新区, 以及坪山新区的中部和东部。这些区域主要是深圳市近些年原特区外重点开发的区域。而作为原特区的南山区、罗湖区、福田区, 由于开发程度较高, 则为深圳市动态度冷点区域。

3.2 街区动态度与土地覆盖变化类型的关系

分析街区内部土地覆盖类别对街区动态贡献度可以发现：不同类别土地覆盖间的转移对街区动态的贡献度具有较大差异(图 6)。其中, 街区贡献度最高的为林地到裸土的转化, 为12.59%；其次为在建用地向建筑、林地向建筑用地, 以及其他不透水地表向裸土的转移, 其动态度贡献度分别为7.21%, 7.00%, 6.20%；而草地和农田与其他类别间的转化对街区动态的贡献度相对较小, 分析原因主要为城市内街区草地与农田的比例很少。值得注意的是, 其他不透水地表向裸土的变化的动态度贡献度较高, 说明城市更新引起的城市内部动态占据重要地位。

对比街区的土地利用强度变化可以发现, 不同街区的土地利用强度变化差异较大, 土地利用强度变化值的范围为-400—337.2。在深圳市4796个街区中, 土地利用强度发生变化的有735个。其中, 土地利用强度变化为负值的有390个, 占变化街区总数的8.13%, 主要分布在城市内部, 如宝安区东部、龙岗区南部、坪山新区、罗湖区、福田区等。而土地利用强度变化为正值的有345个, 占变化街区总数的7.19%, 主要分布在城市外围, 如宝安区西部沿海、龙华新区、龙岗区与坪山新区相邻的区域、南山区北部等。以上街区土地利用强度变化表明, 深圳市城市内部总体上土地开发强度降低, 而城市外围的土地开发强度增加, 可能的原因为城市更新更注重城市生态环境的改善, 增加了街区的生态用地, 使街区的土地开发强度有所降低。

3.3 街区动态度与街区大小的关系

分析街区大小与街区动态度数据可以发现：深圳市街区大小主要集中在1—85.03hm²之间, 街区平均面积为42.22 hm²(图 8), 表明深圳市不同街区的大小存在较大差异, 且街区大小在均值上下分布范围较大。街区动态度则主要分布在0.1%—43.95%, 街区的平均动态度为14.63%(图 8), 说明不同街区的动态度大小也存在较大差异, 需要对不同街区进行的差异化的监管。从街区动态度和街区大小的散点图可以得到：随着街区面积的增加, 街区动态度先急剧下降, 之后下降程度逐渐放缓。利用对数回归进行拟合, 得到街区动态度随街区大小变化的回归方程, 并通过统计上的显著性检验(P < 0.01)。街区动态度与街区大小间的非线性关系表明动态高的街区主要为面积相对较小的街区, 这一结果也为街区尺度上的城市监管提供了参考。

4 讨论 4.1 街区动态度加深了城市动态度变化的理解

研究通过对深圳市城市街区动态度的量化, 在街区尺度上揭示了深圳市内部土地覆盖动态度特征, 加深了对城市生态系统动态变化的理解, 丰富了城市发展过程中的土地覆盖时空动态的变化规律。对比前人研究结果^{[20, 32]}, 深圳市在城市化初期与后期, 城市土地覆盖均具有高度的动态度, 但两个阶段动态变化特征存在明显差异。从动态度空间分布上看, 城市化初期, 深圳市大规模进行围海造陆等城市扩张活动, 土地覆盖变化主要发生在西部沿海和北部的林地^[20], 城市动态度热点主要分布在城市外围。而进入城市化后期, 深圳市城市发展从“增量扩张”向“存量优化”转变, 城市高动态度热点由城市外部向城市内部转移。南山区、盐田区作为原经济特区, 面临着以城市更新为主要方式的“存量优化”改造。《深圳市实施东进战略行动方案(2016—2020)》中提出加快东部城区城市化进程, 打造东部发展新轴带, 实现东西均衡发展^[33]。大鹏新区东部沿海地区作为深圳市着力打造旅游名片, 盐田区作为东进的重要通道都具有较高的动态度。

从变化类别来看, 以往的研究表明深圳市快速的城市化主要通过侵蚀林地、农田、湿地来增加建设用地^[20], 因此向外部蔓延的城市扩张成为深圳市城市化初期动态度的主要贡献者。而现阶段的城市动态则主要来自于城市内部的存量式开发, 如城市更新活动等, 表现为其他不透水地表转化为裸地的变化也具有较高的动态度贡献, 并且城市内部的土地利用强度变化相对城市外围较低, 一定程度上改善了城市人居环境。因此, “存量优化”阶段下的城市动态度监测在兼顾城市绿色基础设施(如城市公园、广场植被等)和蓝色基础设施(如河道、公园水体等)的同时, 也应该关注城市内部建筑以及城市棕地的动态^[34], 以改善城市人居环境, 实现城市的可持续发展。

4.2 城市动态度与城市生态监管

及时准确获得城市动态度信息是城市生态监管和城市规划的基础与前提^[35]。街区作为城市居民活动的主要场所, 已成为城市生态与城市规划研究中的基本基元^{[25, 36]}。街区大小是街区的基本属性之一, 其直接影响城市居民的社会生活, 如出行方式等, 对于城市的社会经济生活具有重要意义^{[18, 24, 37]}。本研究结果显示, 不同大小的街区, 其动态度存在较大差异, 面积较小的街区, 城市动态度较高, 而面积较大的街区具有较低的动态度, 二者存在非线性关系。分析可能的原因为面积较小的街区其内部景观组成较为单一, 对其进行土地开发, 动态度的影响是全局的；而面积较大的街区, 其景观组成相对复杂, 内部往往有较大面积的城市公园等, 土地开发活动频度相对较低, 其动态度影响通常是局部的。同时, 这一结果也为城市动态监测提供了参考。在实际监测过程中, 监管人员可将高动态度的小街区合并到临近街区, 形成一定尺度的“街区组”；而对于动态度低的大街区可进行更为细致的“子街区”分割, 进而在“街区组”和“子街区”水平上进行监管。此外, 监管中应对面积较大且动态度较高的街区进行高频率的监测, 而对于面积较小且动态度较低的街区开展低频率的监测, 从而实现街区的差异化监管, 提高城市监管的效率。

5 结论

城市发展过程中不仅往外扩张, 城市建成区内部也发生变化, 量化城市内部的这种变化, 是理解城市发展的生态环境影响, 有效开展城市生态监管与城市规划的前提和基础。本研究以深圳市为例, 基于街区量化相邻两年城市的土地覆盖动态度, 得到如下结论：(1)深圳市城市内部具有较高的动态度, 2017—2018年间14.93%的街区发生了不同程度的变化。动态度的热点区域主要分布在南山区、宝安区、盐田区等, 冷点则主要分布在龙岗区、坪山新区、福田区等。(2)相较于城市化初期, 现阶段深圳市动态度热点由西部向东部转移, 城市更新引起的土地覆盖变化(如不透水地表转变为裸土等)的动态度贡献较高, 是导致城市内部景观动态的重要因素。与此同时, 深圳市街区的土地利用强度也存在较大差异, 且表现为土地利用强度增加的街区多分布在城市外围, 而街区土地利用强度降低多分布在城市内部。(3)深圳市街区尺度上空间动态度随着街区面积的增大, 先急剧下降, 后趋于平缓, 二者关系呈现非线性关系。研究结果揭示了深圳市内部街区尺度上的动态度特征, 加深了对城市景观动态发展的认识, 为城市规划和生态空间监管提供了重要支撑。