自然资源遥感, 2022, 34(3): 97-105 doi: 10.6046/zrzyyg.2021229

技术应用

基于长时间序列遥感数据的深圳景观连通性动态变化研究

唐文魁,1,2, 俞露,2, 周伟奇1, 岳隽2, 周正2

1.中国科学院生态环境研究中心,城市与区域生态国家重点实验室,北京 100085

2.深圳市城市规划设计研究院有限公司,深圳 518028

Dynamic changes in the landscape connectivity in Shenzhen City determined based on the long time series of remote sensing data

TANG Wenkui,1,2, YU Lu,2, ZHOU Weiqi1, YUE Jun2, ZHOU Zheng2

1. State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China

2. Urban Planning and Design Institute of Shenzhen, Shenzhen 518028, China

通讯作者: 俞 露(1982-),女,博士,教授级高级工程师,主要研究方向为城乡规划、海绵城市。Email:yul@upr.cn。

责任编辑: 张仙

收稿日期: 2021-07-28   修回日期: 2021-10-27  

基金资助: 国家重点研发计划课题“粤港澳大湾区城市群综合决策和协同服务研究与示范”(2019YFB2103104)

Received: 2021-07-28   Revised: 2021-10-27  

作者简介 About authors

唐文魁(1985-),男,博士后,主要研究方向为国土资源规划。Email: tangwk@upr.cn。

摘要

在过去30 a间深圳市的土地利用经历了从剧烈变化到近于饱和的过程。文章基于1988—2015年10期深圳市土地覆盖/利用数据,结合形态学空间格局分析(morphological spatial pattern analysis,MSPA)和图论模型定量分析了深圳生态用地的景观连通性。结果表明,生态用地中耕地在深圳城市化不同阶段均为城市用地扩张的主要土地来源,2005年后林地作为城市土地来源的比例显著提升。景观连通性方面,1988—2015年间深圳市生态用地等效连通面积(equivalent connected area,ECA)减少了1 175.4 km2,网络连接度(degree of network connectivity,DOC)降低了43.51%。自深圳2005年划定基本生态控制线以来,城市生境退化趋势明显放缓,但生态用地ECA依然以每年11.9 km2的速度被逐渐侵蚀。斑块重要性分析表明,阳台山和塘朗山等区域是目前维持深圳市现有景观连通性的关键枢纽节点,需加大保护力度。

关键词： 景观连通性; 图论; 形态学空间格局分析(MSPA); 土地覆盖/利用遥感影像; 深圳市

Abstract

In the past 30 years, the land use in Shenzhen City has changed dramatically until it is almost saturated now. Using the morphological spatial pattern analysis (MSPA) and the graph theory, this study quantitatively analyzed the landscape connectivity of ecological land in Shenzhen based on ten phases of remote sensing images for land cover or use from 1988 to 2015. The results show that the cultivated land was the main land source in various periods of Shenzhen’s rapid urbanization, while the proportion of forest land used for urban development had risen since 2005. For Shenzhen’s landscape connectivity from 1988 to 2015, the equivalent connected area (ECA) of the ecological land decreased by 1 175.4 km2, and the degree of network connectivity (DOC) decreased by 43.51%. Since the delineation of Shenzhen’s basic ecological control boundary in 2005, the pace of urban habitat degradation has slowed down but the ECA of the ecological land had still been gradually eroded at a rate of 11.9 km2 per year. The analysis of the importance of ecological patches shows that areas like the Yangtai Mountain and Tanglang Mountain are key ecological nodes for landscape connectivity and should be protected with greater efforts.

Keywords： landscape connectivity; graph theory; morphological spatial pattern analysis (MSPA); remote sensing image of land cover or use; Shenzhen City

PDF (6140KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

唐文魁, 俞露, 周伟奇, 岳隽, 周正. 基于长时间序列遥感数据的深圳景观连通性动态变化研究[J]. 自然资源遥感, 2022, 34(3): 97-105 doi:10.6046/zrzyyg.2021229

TANG Wenkui, YU Lu, ZHOU Weiqi, YUE Jun, ZHOU Zheng. Dynamic changes in the landscape connectivity in Shenzhen City determined based on the long time series of remote sensing data[J]. Remote Sensing for Land & Resources, 2022, 34(3): 97-105 doi:10.6046/zrzyyg.2021229

0 引言

景观连通性可以衡量景观在多大程度上促进或阻碍栖息地斑块间的生物迁移和其他生态通量,有助于更好地理解生物栖息地破碎化程度,从而避免生物多样性丧失进而严重影响区域生态系统功能和服务[1-2]。通过保护景观连通性,进行针对性的恢复和保护,可以有效改善破碎化对生物多样性的负面影响[3],从而发挥生物多样性在维持城市生态系统稳定中的重要价值。因此,景观连通性量化与评价被认为是生物多样性保护领域的重要课题[4]。

基于图论的景观网络模型自上世纪末被用于景观连通性分析,这些模型指数将复杂的景观简化为节点和链接构成的网络,到目前为止已有许多指数在各项研究中得到成功应用[5-6]。但并非所有指数都适合评估景观连通性的时间序列变化,Baranyi等[7]于2011年系统地比较了13个常用景观连通性指数,指出大部分指数可能存在信息误导传递,且对连通性变化缺乏足够的响应,因此将整体连通性指数(integral index of connectivity,IIC)和可能连通性概率指数(probability of connectivity,PC)作为评价连通性的推荐指数,其中PC指数相对更具有优越性,是唯一既包含了所有斑块空间位置信息,又能有效量化斑块对维持或改善景观功能连通性重要程度的指标[8-9]。为更准确地从空间形态上描述景观结构连通功能,Vogt等[10]结合形态学空间格局分析(morphological spatial pattern analysis,MSPA)对PC指数进行改进,提出基于景观组分的等效连通面积(equivalent connected area,ECA)和网络连接度(degree of network connectivity,DOC)2个指数。由于景观组分是通过结构网络要素进行构建,ECA指数和DOC指数一方面可以避免距离阈值等经验参数对PC指数的影响; 另一方面可以有效与景观内栖息地空间变化联系起来,更适合量化和解释由于景观要素缺失所导致的连通性变化。虽然目前PC指数已被广泛用于景观生态规划研究[11⇓-13],但ECA指数和DOC指数作为PC指数的改进,可更有效地衡量景观尺度上的生态用地可用性,特别是长时间序列下对景观连通性的评估。

经过40 a的改革开放,中国的城市化进程使建设用地大规模扩张。作为中国三大城市群之一的珠江三角洲地区城市扩张尤为剧烈,耕地、森林等土地资源快速消耗,生态负面效应逐渐凸显。深圳作为珠江三角洲地区的主要城市,自1980年成立经济特区以来,经历了快速而大规模的城市化进程。为遏制城市无序蔓延,深圳市早在2005年出台了《深圳市基本生态控制线管理规定》,在全国率先划定了基本生态控制线,并将974.5 km2(接近深圳陆域面积的50%)的土地划入城市生态保护范围,约占全市陆域面积的一半。此前已有部分学者对深圳市生态用地的景观连通性进行研究,吴健生等[14]基于1996—2008年4期深圳市土地变更数据,结合深圳市基本生态控制线,评价控制线内外生态用地景观连通性的变化; 曹翊坤等[15]基于1986—2010年5期Landsat TM遥感影像,对深圳绿色景观连通性进行评价。但在快速发展的地区,景观格局高度演替,在较粗时间尺度上生态系统功能和生态响应过程等详细信息必将大量损失,因此有必要在更长序列、更高精度尺度对深圳市景观连通性进行细致的分析和探讨,挖掘景观过程的细微波动。

本文基于深圳市1988—2015年土地覆盖/利用数据,结合MSPA和图论的景观连通性模型对深圳市生态用地进行分析,从而动态监测1988—2015年深圳市生态用地变化情况,分析其城市化进程中景观演变过程,同时基于MSPA分析和斑块重要性评价,以景观组分为结构单位,通过ECA指数和DOC指数探究生态用地结构与景观连通性的关系。

1 研究数据与方法

1.1 研究区概况

深圳位于中国广东省南部,陆域位置为E113°45'44″~114°37'21″,N22°26'59″~22°51'49″,以亚热带季风气候为主。全境东西宽、南北窄,地势东南高,西北低,大部分为低丘陵地,间以平缓台地(图1)。在1979年被撤县列市之前,深圳是一个仅有约3.3万人口的集镇。经过40 a的发展,土地利用经历了剧烈的变化和显著的城市扩张。截至2018年末,深圳市下辖9个行政区和1个新区,陆域面积为1 997.47 km2,建成区面积为927.96 km2,常住人口为1 302.66万人,城镇化率达100%,是中国第一个全部城镇化的城市。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   深圳市行政区与陆域高程分布

Fig.1   Shenzhen administrative diagram and elevation distribution


1.2 数据来源

本文采用的深圳市土地覆盖/利用数据来源于Dou等[16]。该套成果基于C4.5决策树分类器的AdaBoost算法对多时相的Landsat影像分类,平均以3 a为一个时间间隔(少量影像间隔时间为2 a或5 a),生产了深圳市1988—2015年间10期(包括1988年、1991年、1993年、1996年、2001年、2005年、2008年、2011年、2013年、2015年)的30 m地表覆盖数据。卫星传感器型号包括Landsat5 TM,Landsat7 ETM+,Landsat8 OLI,云量均为0%。所使用的遥感影像均在旱季(1月、2月、11月、12月)拍摄,这4个月份正是华南地区降雨最少、地物反射率变化最小的时期。最终将研究区划分为林地、草地、耕地、建成区、水体和裸地6种土地覆盖/利用类型,平均分类精度为90.04%,平均Kappa系数为0.890(表1)。此外,本文还收集搜集了基本生态控制线等深圳市数据开放平台资料(https://opendata.sz.gov.cn/wap/)。

表1   研究区1988—2015年土地利用分类总体精度[16]

Tab.1  Classification accuracy of land use in study area from 1988 to 2015

年份解译精度
总精度/%Kappa
1988年89.220.876
1990年89.230.898
1993年90.550.891
1996年91.310.901
2001年88.780.876
2005年89.650.893
2008年90.790.889
2011年88.270.872
2013年90.980.896
2015年91.630.912
平均值90.040.890

新窗口打开| 下载CSV


1.3 形态学空间格局分析

MSPA分析是一种基于数学形态学对栅格影像中斑块进行度量和识别的方法。该方法在二值数字图像上通过腐蚀、膨胀、开运算、闭运算进行分割,检测和定位互斥的景观形态特征类型,从而在像元层面上对具有独特功能的斑块要素进行识别[17-18]。

本文利用解译的土地覆盖/利用数据,提取出生态用地类型,包括耕地、林地、草地和水域。将生态用地像元集设置为前景,其他用地类型像元集作为背景,使用GuidosToolbox数字图像工具[10],基于八邻域(点和边)分析识别前景像元集与结构要素之间的空间拓扑关系,最终将前景像元集分割为核心类(core)、孤岛类(islet)、孔隙类(perforation)、边缘类(edge)、桥接类(bridge)、环道类(loop)和支线类(branch)7种互不重叠的斑块类型。在MSPA分析中,不同边缘宽度的设置会对不同斑块类型的空间形态产生较大影响,且斑块边缘区域由于受到相邻斑块和周围环境的影响,通常具有显著的边缘效应,本文结合生态廊道宽度设置要求[19],将边缘宽度设置为60 m。最后根据斑块类型的特征及其含义,识别出对研究区生态用地具有重要连通意义的核心类和桥接类2种斑块类型[20],进一步研究其在景观连通性中的作用。

1.4 景观连通性分析

在MSPA分析所得的7种斑块类型中,只有核心类和桥接类斑块与图论理论相关。核心类和桥接类可以基于图论构建景观网络,其中核心类要素对应于网络节点,桥接类要素对应于这些节点之间的链接。组分是指结构上相互连接的节点和连接节点的链接组成的整体,不同组分之间空间上缺少连接,彼此孤立。

将以组分为景观结构单位,通过ECA指数和DOC指数探究生态用地结构与景观连通性的关系,其中ECA指数具有面积单位,可以很好地与研究区栖息地面积进行比较。ECA和DOC这2个指数均是基于组分面积和数量进行计算,具体公式为[10]:

ECA=∑i=1NCCi2
DOC=ECA∑i=1NCCi×100%

式中: Ci为生态用地中第i个景观组分的面积大小,与研究区具有相同的面积单位量纲; NC为组分数量。景观网络要素之间的联系越来越紧密,往往呈现出更少的组分,当所有要素相互连接成单个组分时,DOC为100%。

1.5 斑块重要性评价

PC指数是基于景观网络的栖息地可用性度量指数,PC值的大小直接反映景观中节点相互连接概率的高低。通过移除或添加景观中某一斑块所引起的PC值变化(dPC)作为相对度量,以衡量单个斑块对整体景观连通性的贡献量[21]。具体公式为:

PC=∑i=1n∑j=1naiajpij*/AL2
dPCk=PC-PCremove,kPC×100%

式中: n为景观中斑块总数; ai和aj分别为斑块i和斑块j的面积;pij*为物种在斑块i和斑块j之间所有可能路径中扩散的最大可能性; AL为整个研究区面积; PCremove,k指除去第k个斑块后PC值。

2 结果与分析

2.1 土地利用变化分析

图2为1988—2015年深圳市土地利用变化。总体变化来看,1988—2015年期间生态用地逐渐趋向萎缩,其中耕地面积从1988年的624.91 km2减少到2015年的151.00 km2,平均每年减少0.88%; 林地减少197.55 km2,年均减少率为0.35%; 水体面积减少有一定程度的波动,从1988年8.7%减少到2015年4.56%; 草地占比极小,占比在0.05%~2.28%之间。与生态用地萎缩相对应,深圳市建成区面积占比从1988年的5.63%快速增加到2015年的41.77%,面积增加721.1 km2,增长6.4倍。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   1988—2015年深圳市土地利用变化

Fig.2   Land use change in Shenzhen City from 1988 to 2015


表2显示了建成区不同土地来源的占比情况。在来源占比分析中,考虑裸地作为一种过渡性城市用地类型,一般由其他土地利用类型转换而来,也能直接反映建设规模和城市扩张方向,因此将裸地和建成区进行归并处理。总体来看,建成区主要是由生态用地中的耕地转化而来,耕地来源的占比达到38.22%~74.08%,在2001—2005年间达到峰值,说明城市建设用地迅速扩张侵占了大量耕地。2005年以后,随着可用耕地资源大大减少,林地作为补充,来源占比显著提升,为城市扩张提供土地资源。除了陆地扩张,自1991年开始深圳市一直通过填海造陆占据水域,导致生态用地中水体面积逐渐减少。

表2   1988—2015年深圳市建成区不同来源占比

Tab.2  Area percentages of other land-use types transforming into built-up areas in each period(%)

阶段来自林地来自草地来自耕地来自水体
1988—1991年35.700.2464.070.00
1991—1993年12.272.3370.8614.53
1993—1996年38.000.0954.467.46
1996—2001年23.542.3864.289.80
2001—2005年12.535.0574.088.33
2005—2008年37.718.1838.2215.89
2008—2011年41.379.6842.386.57
2011—2013年40.663.6542.5113.18
2013—2015年29.9110.6749.859.57

新窗口打开| 下载CSV


根据深圳市建成区面积增长速率的不同,将深圳市的城市扩张大致分为3个阶段(图2(k)): 早期城市化阶段(1988—1993年)、快速城市化阶段(1993—2005年)和集约城市化阶段(2005—2015年)。早期城市化阶段(1988—1993年),建成区范围限于南山区南部、福田区东部以及罗湖区西南部等深圳市南部沿海地区; 快速城市化阶段(1993—2005年),建设用地增长主要集中在深圳市中西部区域,以西部的宝安区和北部的龙岗区、龙华区最为显著。集约城市化阶段(2005—2015年),深圳市城镇化率达到较高水平,特别2013年之后城市用地增长趋于饱和,深圳市逐渐进入存量发展时期。

2.2 基于MSPA的景观格局分析

基于MSPA分析将研究区生态用地识别为7种不同的斑块类型。一般来说,不同斑块类型具有不同的生态学意义[22]。核心类斑块是区域内较大的栖息地斑块,具有生态源地的功能; 孔隙类和边缘类斑块是核心类斑块与非生态区的过渡区域,具有边缘效应; 桥接类斑块是指衔接不同核心类斑块的连通区域,可作为构建生态安全格局的廊道; 环道类斑块与桥接类斑块类似,是在同一核心区内部起联系作用。孤岛类、支线类是指面积远小于核心区且连通性较低的细小、破碎斑块。由图3可见,核心类斑块作为主要生境类型,面积占比最高,而支线类、孤岛类斑块由于破碎化程度高,斑块数量最多,这与其斑块类型的生态学意义相对应。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   1988—2015年深圳市MSPA斑块类型面积和数量的时间序列变化

Fig.3   Time series changes of the area and number of MSPA types in Shenzhen City from 1988 to 2015


对不同斑块类型的面积占比与斑块数量进行Spearman相关性分析,在0.05的置信水平下,核心类斑块面积占比和斑块数具有显著负相关(图3(a)),表明随着核心类斑块面积减少,斑块数量逐渐增加,破碎化程度越高。相反,除边缘类和环道类外,孤岛类、孔隙类、桥接类和支线类斑块面积和斑块数均在时间序列上有相同的趋势性,相关关系表现为显著正相关(图3(b)、(c)、(e)、(f)),即随着斑块数量增加,这些斑块类型的面积也逐渐增加。

从时间序列的变化趋势上,不同斑块类型表现出不同的变化趋势。核心类斑块面积呈逐渐下降的趋势,从80.56%降至37.48%。与之相反,孤岛类和支线类斑块表现出持续增长的趋势,面积占比分别从0.24%增长至2.17%,0.78%增长至2.01%。孔隙类、边缘类、桥接类、环道类等其余类型斑块变化轨迹较复杂,呈先增长后降低的趋势,均在1996年或2001年达到增长峰值。

结合相关分析,说明在早期和快速城市化阶段,核心类斑块被城市用地侵占的同时,部分核心类斑块转化为孔隙类等过渡型斑块或桥接类等狭长型斑块,还有部分转化孤岛类等破碎型斑块,导致孤岛类、孔隙类、桥接类等斑块类型在这一阶段面积占比和斑块数量都呈增长趋势。而在2001年后,政府为保障城市生态安全,开始对仅存的核心类斑块进行保护,但这一阶段随着深圳市人口规模增长,人地矛盾压力增大,缺乏保护的孔隙类和桥接类等斑块类型更容易受到城市化扰动,进一步孤岛化,导致这一阶段孔隙类和桥接类斑块成为孤岛类和支线类斑块的主要来源,斑块面积占比和数量同步降低。

2.3 基于MSPA的连通性分析

以景观组分为结构单位,通过ECA指数和DOC指数探究生态用地结构与景观连通性的关系。景观组分是基于生态网络中能相互链接的节点构建的组团结构,景观连通性反映了景观中不同组分在空间格局上的联系。图4和表3显示了1988—2015年深圳市生态用地的组分分布和连通性变化情况,其中不同组分在图4中以不同颜色交替显示,颜色本身无实际具体含义,因此不同年份间的相同颜色无可比性,颜色越多表明生态用地缺乏网络要素联系,被分割成更多的组分。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   深圳市生态用地网络组分分布

Fig.4   Distribution of ecological land components and the basic ecological controlling line in Shenzhen City


表3   1988—2015年深圳市生态用地及各类用地连通性变化

Tab.3  Dynamic of ecological land connectivity in Shenzhen City from 1988 to 2015

年份生态用地林地耕地水体
组分数/个ECA/km2DOC/%ECA/km2DOC/%ECA/km2DOC/%ECA/km2DOC/%
1988年931 589.6198.90395.7562.06150.5795.7863.7156.88
1991年1391 457.7898.78404.1460.4253.3160.7476.9066.56
1993年2011 300.8297.78397.0360.4224.3939.3165.8062.63
1996年2051 018.2387.34374.2659.0619.0536.5643.0249.90
2001年238897.2688.17266.3947.609.3332.2133.7446.86
2005年390533.5258.33348.3862.882.1124.9123.6343.85
2008年412545.1561.30341.0260.751.4222.2615.3032.50
2011年358463.6552.60260.7547.012.7832.1816.4133.95
2013年585348.5749.58240.8348.771.1424.9313.8134.84
2015年449414.2155.39248.7149.771.0720.9814.9938.02

新窗口打开| 下载CSV


由图4和表3可知,1988—2015年期间深圳市斑块数量从457逐渐增加至1 524,表明生态用地随时间序列逐渐被分割隔离。通过式(1)—(2)计算ECA和DOC指数对连通性进行评估。1988—2015年期间随着组分增加,生境斑块之间的连通可能性减小,生态用地ECA值相应减少1 175.4 km2,DOC值降低43.51%。生态用地连通性与耕地连通性呈现较好的相关性(R=0.794,p=0.006),与林地、水体相关性不显著,表明耕地面积减少是导致生态用地连通性降低的主要原因。

对应深圳市城市化不同阶段,在早期城市化阶段(1988—1993年),景观之间的联系紧密,呈现出较少的组分,研究区生态组分数量在93~201个之间,ECA值减少288.79 km2,每年减少约57.76 km2,DOC值基本在97%以上。这一阶段虽然区域内ECA值有所减少,但并未造成DOC指数大幅降低。快速城市化阶段(1993—2005年),深圳市中西部区域、东北部地区的生态组分分布趋于分散,组分数量增加至390个,保存较为完整的组分主要分布在坡度大于25%的丘陵台地区,以大鹏新区、盐田区为主。相应的,这段时期ECA值每年约减少63.94 km2,且DOC指数下降明显,特别在2001年与2005年之间存在明显的拐点,DOC值从88.17%突降至58.33%,下降了约29百分点。2005年之后,深圳市区域生态组分分布逐渐趋于稳定,ECA值降低119.31 km2,DOC值维持在49.58%~61.30%范围。

研究表明,为了避免地区生物多样性被毁灭性破坏,40%的栖息地丧失被认为是允许的最大损失阈值[23]。本文研究表明,2005年以前深圳生态用地网络连接度DOC指数下降明显,基本已触发40%的损失阈值; 2005年基本生态控制线划定后,景观连通性下降得到遏制,DOC指数维持在50%上下波动,现存的生态组分边界基本与深圳市划定的基本生态控制线保护范围相一致(图4(k))。

2.4 斑块重要性分析

景观图论网络允许评估每个单独的网络要素对于维持生态用地连通性的重要程度。基于MSPA分析提取出核心类斑块,通过式(3)—(4)可计算斑块重要性指数dPC。对计算出的dPC进行归一化处理,然后对重要斑块进行1—5自然断点分级,值越大表示斑块的重要性越高。部分关键年份的结果如图5所示。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   深圳市生态用地斑块连通重要性评价

Fig.5   Evaluation of the connectivity importance of ecological patches in Shenzhen City


由图5可知,研究区高重要性等级(4—5级)斑块在时间序列上呈现出逐渐被等级低的区域蚕食的趋势,以西部宝安区、光明区和中部龙岗区、龙华区斑块重要值下降最明显。1988年,研究区范围内4—5级重要斑块分布较广,除宝安区、南山区、福田区、罗湖区部分沿海地区外,其余地区几乎全被4—5级重要斑块覆盖。1996年,4—5级重要斑块主要分布于深圳东部大鹏新区、盐田区、光明区、坪山区全区以及龙岗区北部区域,宝安区和南山区北部区域也有部分分布; 2001年,4—5级重要斑块大幅下降,仅维持在深圳市中西部小部分范围,东部区域大范围转变为3级斑块。2008年斑块等级分布与2001年大致类似,深圳市中部、西部区域大部分被等级1的区域取代。至2015年,3级斑块仅在盐田区、大鹏新区和坪山区北部区域分布; 4—5级重要斑块仅维持在龙华区与宝安区交界的阳台山、龙华区与福田区交界的塘朗山等区域。

3 讨论

深圳市是目前世界上人口密度、建设强度最高的城市之一,土地资源面临难以为继的难题。生态用地遭到侵蚀,景观破碎度增加,ECA指数减少是近30 a深圳生态景观连通性的主要变化特征。

2005年深圳市首先划定基本生态控制线,对城市自然生态空间起到一定的优化和锚固作用。但本文研究表明2005年以后深圳市生态用地ECA指数降低119.31 km2,仍以每年11.9 km2的速度被逐渐侵蚀。说明按照基本生态控制线现有的管控要求,依然无法阻止过度城市化所造成的生态破坏。虽然后续政府部门多次修订《深圳市人民政府关于进一步规范基本生态控制线管理的实施意见》加强管控,但控制线作为一条动态的底线频频突破,暴露了目前基本生态控制线存在刚性管控和精细化管控不足的问题。比如,实施意见中允许生态环境保护相适宜的重大道路交通设施、市政公用设施等项目经审批后,占用基本生态控制线范围进行建设的。然而这些项目类别的内涵却没有细化说明与界定,导致许多开发项目因定义模糊而不断在生态控制线内新增建设布局。至2015年,线内新增建设用地面积46.1 km2,年均新增4.6 km2,至2015年深圳市基本生态控制线内生态用地减少至841.8 km2,仅占86.4%[24]。说明目前基本生态控制线对深圳市城市化进程只起到抑制和延缓作用,现有政策约束很难从根本上提高深圳市整体生态质量。这一结果也与其他学者分析结果相一致[25-26]。

深圳市生态用地面积和连通度下降会带来一系列生态风险,因此有必要基于连通性评估,重新梳理深圳市现有的生态网络结构,预控现有生态廊道内关键衔接点空间。本文分析表明,至2015年阳台山、塘朗山等区域是目前维持深圳市现有生态网络格局的关键枢纽节点。这些节点一旦遭到破坏会使深圳市生态网络进一步恶化,因此需要加大保护力度。

4 结论

1988—2015年深圳市土地覆盖/利用变化表明,27 a间深圳市建成区面积增加了6.4倍。耕地在深圳市城市化不同阶段均为城市用地扩张的主要土地来源,2005年后林地作为城市土地来源的比例显著提升。1988—2015年间深圳市生态用地等效连通面积降低了73.9%,网络连接度降低43.5%。2005年以来,深圳市城市生境退化趋势明显放缓,但生态用地等效连通面积依然以每年11.9 km2的速度在逐渐减少。至2015年重要“踏脚石”斑块主要位于龙华区与宝安区交界的阳台山、龙华区与福田区交界的塘朗山等区域。

本文将MSPA形态学和景观图论模型有效结合,以景观组分为结构单位探究生态用地结构与景观连通性的关系,为区域景观连通性评价提供了一种快速而合理的方法。评价结果可以进一步对廊道修复区域进行优先级识别,为后续景观生态保护提供切实可行的方案。

参考文献

Cotter M, Häuser I, Harich F K, et al.

Biodiversity and ecosystem services:A case study for the assessment of multiple species and functional diversity levels in a cultural landscape

[J]. Ecological Indicators, 2017, 75:111-117.

DOI:10.1016/j.ecolind.2016.11.038      URL     [本文引用: 1]

Liu W, Hughes A C, Bai Y, et al.

Using landscape connectivity tools to identify conservation priorities in forested areas and potential restoration priorities in rubber plantation in Xishuangbanna,Southwest China

[J]. Landscape Ecology, 2020, 35(2):389-402.

DOI:10.1007/s10980-019-00952-2      URL     [本文引用: 1]

Acevedo M A, Sefair J A, Smith J C, et al.

Conservation under uncertainty:Optimal network protection strategies for worst-case disturbance events

[J]. Journal of Applied Ecology, 2015, 52(6):1588-1597.

DOI:10.1111/1365-2664.12532      URL     [本文引用: 1]

齐珂, 樊正球.

基于图论的景观连接度量化方法应用研究——以福建省闽清县自然森林为例

[J]. 生态学报, 2016, 36(23):7580-7593.

[本文引用: 1]

Qi K, Fan Z Q.

Evaluation method for landscape connectivity based on graph theory:A case study of natural forests in Minqing County,Fujian Province

[J]. Acta Ecologica Sinica, 2016, 36(23):7580-7593.

[本文引用: 1]

Wang L, Young S S, Wang W, et al.

Conservation priorities of forest ecosystems with evaluations of connectivity and future threats:Implications in the Eastern Himalaya of China

[J]. Biological Conservation, 2016, 195:128-135.

DOI:10.1016/j.biocon.2015.12.044      URL     [本文引用: 1]

Shanthala D B S, Murthy M S R, Debnath B, et al.

Forest patch connectivity diagnostics and prioritization using graph theory

[J]. Ecological Modelling, 2013, 251:279-287.

DOI:10.1016/j.ecolmodel.2012.12.022      URL     [本文引用: 1]

Baranyi G, Saura S, Podani J, et al.

Contribution of habitat patches to network connectivity:Redundancy and uniqueness of topological indices

[J]. Ecological Indicators, 2011, 11(5):1301-1310.

DOI:10.1016/j.ecolind.2011.02.003      URL     [本文引用: 1]

Saura S, Pascual-Hortal L.

A new habitat availability index to integrate connectivity in landscape conservation planning:Comparison with existing indices and application to a case study

[J]. Landscape and Urban Planning, 2007, 83(2):91-103.

DOI:10.1016/j.landurbplan.2007.03.005      URL     [本文引用: 1]

Saura S, Torné J.

Conefor Sensinode 2.2:A software package for quantifying the importance of habitat patches for landscape connectivity

[J]. Environmental Modelling and Software, 2009, 24(1):135-139.

DOI:10.1016/j.envsoft.2008.05.005      URL     [本文引用: 1]

Vogt P, Riitters K.

GuidosToolbox:Universal digital image object analysis

[J]. European Journal of Remote Sensing, 2017, 50(1):352-361.

DOI:10.1080/22797254.2017.1330650      URL     [本文引用: 3]

陈昕, 彭建, 刘焱序, 等.

基于“重要性—敏感性—连通性”框架的云浮市生态安全格局构建

[J]. 地理研究, 2017, 36(3):471-484.

DOI:10.11821/dlyj201703006      [本文引用: 1]

生态安全格局构建是保障城市生态安全,实现城市可持续性发展的重要途径,也是景观生态学研究的热点领域之一。生态源地识别与空间阻力面构建一直是生态安全格局构建中的技术难点。以广东省云浮市为例,从生态系统服务重要性、生态敏感性与景观连通性三个方面识别生态源地,利用DMSP/OLS夜间灯光数据修正基本生态阻力面,并运用最小累积阻力模型判定生态廊道,从而综合构建云浮市生态安全格局。研究表明：云浮市生态源地占全市总面积的36.47%,主要分布在西部与南部的山林地。云浮市生态廊道总长度为508.87 km,其中景观廊道总长度为315.58 km,组团廊道总长度为193.29 km,呈环状辐射分布于植被覆盖相对较好的山区地带。全市16个自然保护区基本都位于生态源地范围内。构建的“重要性—敏感性—连通性”框架可以为区域生态安全格局构建提供新思路,从而有效指引相关空间规划。

Chen X, Peng J, Liu Y X, et al.

Constructing ecological security patterns in Yunfu City based on the framework of importance-sensitivity-connectivity

[J]. Geographical Research, 2017, 36(3):471-484.

[本文引用: 1]

倪庆琳, 侯湖平, 丁忠义, 等.

基于生态安全格局识别的国土空间生态修复分区——以徐州市贾汪区为例

[J]. 自然资源学报, 2020, 35(1):204-216.

[本文引用: 1]

Ni Q L, Hou H P, Ding Z Y, et al.

Ecological remediation zoning of territory based on the ecological security pattern recognition:Taking Jiawang District of Xuzhou City as an example

[J]. Journal of Natural Resources, 2020, 35(1):204-216.

DOI:10.31497/zrzyxb.20200117      URL     [本文引用: 1]

王博娅, 刘志成.

北京市海淀区绿地结构功能性连接分析与构建策略研究

[J]. 景观设计学, 2019, 7(1):34-51.

[本文引用: 1]

Wang B Y, Liu Z C.

Functional connectivity analyses and construction strategies of the structure of green space network in Haidian District of Beijing

[J]. Landscape Architecture Frontiers, 2019, 7(1):34-51.

DOI:10.15302/J-LAF-20190104      URL     [本文引用: 1]

吴健生, 刘洪萌, 黄秀兰, 等.

深圳市生态用地景观连通性动态评价

[J]. 应用生态学报, 2012, 23(9):2543-2549.

[本文引用: 1]

Wu J S, Liu H M, Huang X L, et al.

Dynamic evaluation on landscape connectivity of ecological land:A case study of Shenzhen,Guangdong Province of South China

[J]. Chinese Journal of Applied Ecology, 2012, 23(9):2543-2549.

[本文引用: 1]

曹翊坤, 付梅臣, 谢苗苗, 等.

基于LSMM与MSPA的深圳市绿色景观连通性研究

[J]. 生态学报, 2015, 35(2):526-536.

[本文引用: 1]

Cao Y K, Fu M C, Xie M M, et al.

Landscape connectivity dynamics of urban green landscape based on morphological spatial pattern analysis (MSPA) and linear spectral mixture model (LSMM) in Shenzhen

[J]. Acta Ecologica Sinica, 2015, 35(2):526-536.

[本文引用: 1]

Dou P, Chen Y.

Dynamic monitoring of land-use/land-cover change and urban expansion in Shenzhen using Landsat imagery from 1988 to 2015

[J]. International Journal of Remote Sensing, 2017, 38(19):5388-5407.

DOI:10.1080/01431161.2017.1339926      URL     [本文引用: 2]

Velázquez J, Gutiérrez J, Hernando A, et al.

Evaluating landscape connectivity in fragmented habitats:Cantabrian capercaillie (tetrao urogallus cantabricus) in Northern Spain

[J]. Forest Ecology and Management, 2017, 389:59-67.

DOI:10.1016/j.foreco.2016.12.006      URL     [本文引用: 1]

Soille P, Vogt P.

Morphological segmentation of binary patterns

[J]. Pattern Recognition Letters, 2009, 30(4):456-459.

DOI:10.1016/j.patrec.2008.10.015      URL     [本文引用: 1]

朱强, 俞孔坚, 李迪华.

景观规划中的生态廊道宽度

[J]. 生态学报, 2005, 25(9):2406-2412.

[本文引用: 1]

Zhu Q, Yu K J, Li D H.

The width of ecological corridor in landscape planning

[J]. Acta Ecologica Sinica, 2005, 25(9):2406-2412.

[本文引用: 1]

史芳宁, 刘世梁, 安毅, 等.

城市化背景下景观破碎化及连接度动态变化研究——以昆明市为例

[J]. 生态学报, 2020, 40(10):3303-3314.

[本文引用: 1]

Shi F N, Liu S L, An Y, et al.

Changes of landscape fragmentation and connectivity with urbanization:A case study of Kunming City

[J]. Acta Ecologica Sinica, 2020, 40(10):3303-3314.

[本文引用: 1]

Saura S, Rubio L.

A common currency for the different ways in which patches and links can contribute to habitat availability and connectivity in the landscape

[J]. Ecography, 2010, 33(3):523-537.

[本文引用: 1]

于亚平, 尹海伟, 孔繁花, 等.

南京市绿色基础设施网络格局与连通性分析的尺度效应

[J]. 应用生态学报, 2016, 27(7):2119-2127.

DOI:10.13287/j.1001-9332.201607.006      [本文引用: 1]

以南京市为研究区,基于ArcGIS、Erdas、GuidosToolbox和Conefor等软件平台,采用形态学空间格局分析(MSPA)和景观连通性分析方法,通过在MSPA中设置不同的粒度、边缘宽度和在景观连通性分析中设置不同的扩散距离阈值,对2013年南京市绿色基础设施网络格局变化的尺度效应、边缘效应与距离效应进行评价.结果表明: 基于MSPA获取的景观类型构成存在明显的尺度效应和边缘效应,且边缘效应对MSPA景观类型的影响较尺度效应更为明显.不同扩散距离对景观连通性的影响很大,对于南京市来说,2 km或2.5 km是关键的扩散距离阈值.当输入数据选择粒度30 m、边缘宽度30 m时,可以得到南京市城市绿色基础设施(UGI)网络更为详尽的景观信息.基于MSPA与景观连通性方法,分析尺度效应、边缘效应、距离效应对研究区UGI网络景观类型的影响,有助于选择合适的粒度、边缘宽度及扩散距离,并更好地理解UGI网络的空间格局和与生态过程相关的尺度效应和距离效应,从而使得UGI网络时空格局变化分析时的参数设置更为科学合理.研究结果可为中尺度范围内UGI景观网络时空格局分析时的参数设置提供重要的参考和依据,对其他地区UGI景观网络的分析也具有一定的借鉴意义.

Yu Y P, Yin H W, Kong F H, et al.

Scale effect of Nanjing urban green infrastructure network pattern and connectivity analysis

[J]. Chinese Journal of Applied Ecology, 2016, 27(7):2119-2127.

[本文引用: 1]

Yin D, Leroux S J, He F.

Methods and models for identifying thresholds of habitat loss

[J]. Ecography, 2017, 40(1):131-143.

DOI:10.1111/ecog.02557      URL     [本文引用: 1]

邓金杰, 黄伟坚.

城市生态用地管控:问题及规划管理建议——以深圳市为例

[J]. 城市观察, 2018(3):131-138.

[本文引用: 1]

Deng J J, Huang W J.

Problems and suggestions on urban ecological land use control:A case study of Shenzhen

[J]. Urban Insight, 2018(3):131-138.

[本文引用: 1]

张甜, 王仰麟, 刘焱序, 等.

1987—2015年深圳市主城区景观演变过程多时相识别

[J]. 地理学报, 2016, 71(12):2170-2184.

DOI:10.11821/dlxb201612008      [本文引用: 1]

城市扩张具有典型的阶段性特征,城市化初期的不透水面快速蔓延和中后期的生态环境约束往往呈现非线性关系。基于多时相数据对典型城市发展中关键指标的变化拐点进行识别,将有助于认识城市景观演变的时间节点,理解城市化的生态响应过程。本文选用深圳市西部1987-2015年长时序Landsat影像共27期,逐年提取新构建的归一化裸露指数（MNDBI*）与归一化植被指数（NDVI*）,从而在时间上寻找城市增长的转折点,在空间上识别不同空间位置的景观演变特征。研究结果表明,1987-2015年深圳市城市增长十分明显,并以2003年为拐点经历了“快速增长”至“平稳约束”的变化;与此对应,植被指数的关键拐点与城市化拐点基本重合,印证了城市建设用地扩张和生态系统响应的时空关联特征。此外,从空间分异来看,深圳市南部各区的城市化约束期出现相对更早,拐点一般在1995-1998年间;而北部各区的快速城市化时期持续更长,一般在2003-2006年后趋于平稳。在空间上,以市图书馆为中心刻画西、北、西北、东北剖面线方向的指标动态,发现采样中心附近的景观类型变化较小,而更远辐射半径经历了更大幅度的城市化;这说明深圳市南部城市发展较早达到饱和,且近30年的城市扩张以由南向北的放射式蔓延为主。面临城市化带来的生态威胁,相关环保措施的有力实施仅能延缓城市化导致的生境退化步伐,城市扩张所带来的生态破坏依然不容小觑,合理而有力的政策颁布、实施与监管在未来的城市发展中极为必要。

Zhang T, Wang Y L, Liu Y X, et al.

Multi-temporal detection of landscape evolution in western Shenzhen City during 1987—2015

[J]. Acta Geographica Sinica, 2016, 71(12):2170-2184.

[本文引用: 1]

吴健生, 黄乔, 曹祺文.

深圳市基本生态控制线划定对生态系统服务价值的影响

[J]. 生态学报, 2018, 38(11):3756-3765.

[本文引用: 1]

Wu J S, Huang Q, Cao Q W.

Effects of basic ecological control line policy on ecosystem services value in Shenzhen

[J]. Acta Ecologica Sinica, 2018, 38(11):3756-3765.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-2
版权所有 © 2015 《自然资源遥感》编辑部
地址:北京学院路31号中国国土资源航空物探遥感中心 邮编：100083
电话:010-62060291/62060292 E-mail:zrzyyg@163.com
本系统由北京玛格泰克科技发展有限公司设计开发