1. 前言
随着气候变化、人口增长、环境污染、能源短缺等诸多全球性问题的日益突兀 [1],相较于发达国家,我国人均土地可利用资源稀少的现实与当前全面城市化趋势的冲突也越来越不可避免。在当前背景之下,人们开始逐渐把城市空间拓展和土地利用结构调整视为解决这个问题的答案。城镇用地扩展是土地利用和土地覆被变化最重要的类型之一。作为城市空间扩展最有力的土地表征类型,城镇用地变化的形态、结构等特征,在学者研究城市扩展这一方面有着极其重要的指示作用 [2] [3] [4] [5]。同时,由此引发的快速城市化不仅显著影响着城市原有景观格局,甚至导致如生物多样性减少和栖息地破坏 [6]、热岛效应增强 [7]、环境污染、水土流失 [8] 等一系列环境问题。景观格局(Landscape pattern)是诠释并表现空间结构特征的专业术语,是景观组成单元的类型,数目与空间配置。是景观异质性在空间上的综合表现 [9]。其中,景观格局是指不同大小和形状的斑块在空间上的排列状况,是人为因素和自然因素共同作用的结果,具有显著的时空特征 [10] [11]。目前,前人大多在大尺度、宏观方面对LUCC的景观格局进行研究,采用多样化的方法和模式:通过GIS技术获取影像数据中土地利用数据,通过矩阵转化反映土地景观格局的变化,应用生态学景观指数方法定量分析LUCC的规律,从而为区域尺度上LUCC的研究积累了丰富的成果和可以借鉴的经验。例如,马晓勇等 [12]、韩锦涛等 [13]、哈孜亚•包浪提将等 [14]、阳文锐等 [15]、何丹等 [16] 分别对太原市、乌鲁木齐市、北京市、京津冀都市圈的城市景观格局变化进行了研究。同时,本文还采用ARCGIS与城镇用地定量分析相结合的方法探寻青岛市城镇用地扩张的变化。李晓文等人以海为目标,针对性的采用城市用地扩展的强度以及空间分层的特征来分析,对上海的城镇用地变化情况的拓展区域和模式做出了总结 [17];Mundia CN等人利用遥感影像开始研究并最终得出了结论,城镇用地扩展的代价是部分绿地的减少,并且和道路网有着直接相关的线性关系,并以此为依据提出一些提高未来城市空间使用效率的建议 [18];Tan KC等人同样使用类似的方法研究了马来西亚的槟榔屿,不同的是他们采用了陆地表面温度和植被指数这两个指数来分析工地和居住用地对城镇用地扩展的影响 [19];Trusilova K等人以气候变化对城市扩展的影响这一角度切入,对欧洲的相关领域做了初步的探索 [20]。本文基于青岛市重要的社会、经济发展地位以及其城市景观格局变化的研究还属于空白的原因,采用GIS技术与量化分析指数相结合的方法分析青岛市一定时期内二者的相关变化,为未来城镇用地扩展方向提供指导与建议。
2. 研究区概况与研究方法
2.1. 研究区概况
本文所讨论的城市青岛位于中国山东省境内,地理位置在山东半岛东南部,胶东半岛的东部,青岛市濒临黄海,与朝鲜韩国隔海相望,在陆地上与烟台市,潍坊市和日照市接壤。青岛市的总面积为11282平方公里,市一下有七个市,同时管理着三个县级市。青岛地理位置优越,是我国沿海及其重要的中心城市港口城市,同时也是我国发达的沿海旅游度假城市,是我国极其重要的海洋发展先进去,更是亚洲东北地区的航运纽带,因为其独特的优越性,被我国“一带一路”政策选定为新亚洲欧洲经济走廊的主要关节城市,也是海上战略合作重要指点。如今我国“一带一路”战略政策实施良好,国家“一三五”计划不断推行,青岛市走在了经济发展的前沿,这使得青岛市的城市化进程不断提高,因此青岛市成为了山东省划定的蓝色经济区的火车头,旅游业也逐渐发展成为亚洲旅游胜地。截止到2018年,青岛市总人口达到了939万,完成了12002亿元的GDP,在旅游业的开发商也成为了国家的模范。
2.2. 数据源处理与方法
本文所用的数据是青岛市历年土地利用数据,来源于中国科学院提供的1973、1980、1990、2000、2010和2018年等6期青岛市矢量土地利用数据。
依据影像光谱特征,结合野外勘察的资料,同时参照相关地理资料,与Google Earth对照,在专家参与下,从地物的几何形状,颜色特征、纹理特征和空间分布情况进行分析,并在综合各位专家意见后,建立遥感影像解译标志。通过解译最终获得青岛市土地利用数据,将景观类型分为耕地、林地、草地、水域、未利用地和城镇用地共6类,获得不同时期景观类型分类图后,在ArcGIS软件支持下,将矢量文件转换为栅格数据文件,为保持数据处理的精度将像元大小设置为60 m。在Fragstats4.0软件的支持下,将栅格数据文件导入,选取景观格局指数并计算,将处理结果导入XCEL中进行进一步的分析与处理(如图1 )。
2.3. 对应指标选取
2.3.1. 城市景观对应指标选取
对城市景观的研究方法具有多样性,经过文献阅读综述和对城市景观指标总结筛选之后,在景观类型变化方面选取结构和变化幅度模型、土地利用转移模型,在景观格局方面选取景观指数的方法,具体建立并选取以下模型和指标来进行分析。
变化程度模型
景观类型结构模型(1)是反应不同的景观类型占总体的比重。其数学表达式为:
Ua = Sm/S * 100%
式中Sm代表某景观类型的面积,S代表某时段区域总面积,Ua为某景观结构比例。
景观类型变化幅度模型(2)同景观类型面积变化的程度。其数学表达式为:
![](//html.hanspub.org/file/7-1040093x10_hanspub.png)
Figure 1. Land use map of Qingdao from 1973 to 2018
图1. 1973~2018年青岛市土地利用图
式中Sm、Sn分别代表某景观类型的面积,Ud为某一景观类型的变化幅度。
矩阵转移模型
土地利用转移矩阵反映了某一区域某一时段初期和末期各景观类型面积之间相互转化的动态过程信息 [21] [22] [23] [24] [25]。
2.3.2. 城镇用地指标变化选取
在城镇用地的变化上,选取了以下几种指标对应分析
扩展强度指数。扩展强度指数的定义是平均每年城镇用地的扩展面积占研究区域总面积的比重 [26] [27] [28]。主要表现城镇用地扩展状况相对于研究区的强弱程度,公式为:
式中UII是城市扩展强度指数;Ua、Ub为起止时相 城镇用地数量;TLA为研究单元土地总面积;T为研究时段的时间跨度。•需要注意的是,由于数据间隔不同,T取值也不同。
紧凑度指数。
公式为
BCI为城市外部形态的紧凑度,A为建成区面积,P为建成区外围轮廓的周长。城市紧凑度指数值越大,说明其形状越有紧凑性,反之,说明形状的紧凑性越差
分维指数
公式为
D为城市建设用地斑块的分维指数,P、A分别表示斑块的面积和周长。D值越大表明形态越复杂、边界越不规则。
3. 结果分析
3.1. 景观类型结构组成及变化分析
结合青岛市1973~2018时段6期的景观类型分布图(图1)、结构模型和变化幅度模型得到1973~2018青岛市各景观类型构成及变化(表1)及1973~2018年青岛市各景观类型结构。
![](Images/Table_Tmp.jpg)
Table 1. Composition and changes of various landscape types in Qingdao from 1973 to 2018 Unit: Km2
表1. 1973~2018青岛市各景观类型构成及变化 单位:Km2
从各景观类型构成来看,1973年到2018年间共计6年的主要景观类型都是耕地、建设用地和草地,所占比例不同,每一年的类型分别为67.81%、12.66%、8.89%;67.78%、12.75%、8.85%;67.69%、13.25%、8.78%;66.76%、14.30%、8.76%;67.74%、18.65%、2.97%;65.69%、21.45%、2.92%攻击六份数据。
总体来看1973~2018年各时段的主要景观类型均为耕地、建设用地和草地,但在各时段所占比例不同,耕地面积在1973~1980年变化幅度是0.06%,呈现增长状态,在1980~1990年变化幅度是减少0.08%,在1990~2000年变化幅度是减少0.94%,在2000~2010年变化幅度是增长1.04%,在2010~2018年变化幅度是减少1.79%,呈现下降状态,在1973~2018年变化幅度是减少1.72%,总体表现为下降趋势。
林地面积在1973~1980年变化幅度为0,表现为稳定,在1980~1990年变化幅度是减少0.02%,呈现下降趋势,在1990~2000年变化幅度为0,表现为稳定,在2000~2010年变化幅度是减少1.10%,呈现为下降趋势,在2010~2018年变化幅度是增加0.03%,呈现为增长趋势,在1973~2018年变化幅度为较少1.09%,表现为下降趋势。
草地面积在1973~1980年变化幅度为−0.03%,呈现下降趋势,在1980~1990年变化幅度为−0.08%,呈现下降趋势,在1990~2000年变化幅度为−0.02%,呈现下降趋势,在2000~2010年变化幅度为−5.78%,呈现下降趋势,在2010~2018年变化幅度为−0.04%,呈现下降趋势,在1973~2018年变化幅度为−5.95%,总体表现为下降趋势。
水域面积在1973~1980年变化幅度为0.00%,表现为稳定,在1980~1990年变化幅度为−0.29%,呈现下降趋势,在1990~2000年变化幅度为−0.09%,呈现下降趋势,在2000~2010年变化幅度为1.37%,呈现增长趋势,在2010~2018年变化趋势为−0.69%,呈现下降趋势,在1973~2018年变化幅度为0.67%,总体表现为增长趋势。
建设用地面积在1973~1980年变化幅度为0.10%,呈现增长趋势,在1980~1990年变化幅度为0.51%,呈现增长趋势,在1990~2000年变化幅度为1.05%,呈现增长趋势,在2000~2010年变化幅度为4.36%,呈现增长趋势,在2010~2018年变化趋势为2.88%,呈现增长趋势,其中2000~2010年增长幅度最大,约为其他时段的2~4倍,在1973~2018年变化幅度为8.92%,总体表现为增长趋势。
未利用地面积在1973~1980年变化幅度为0.00%,表现为稳定,在1980~1990年变化幅度为−0.05%,呈现下降趋势,在1990~2000年变化幅度为0.00%,表现为稳定,在2000~2010年变化幅度为−0.17%,呈现下降趋势,在2010~2018年变化趋势为−0.01%,呈现下降趋势,在1973~2018年变化幅度为−0.22%,总体表现为下降趋势。
景观类型变化转移矩阵
根据获取的1973、1980、1990、2000、2010、2018年6期的影像数据,通过ArcGIS10.3软件对1973~1980、1980~1990、1990~2000、2000~2010、2010~2018共5个时段的景观类型变化进行转移矩阵分析,得到1973~2018年青岛市各时段景观类型变化转移矩阵(表2),各时段主要景观类型转移变化(图2)。
![](Images/Table_Tmp.jpg)
Table 2. Transfer matrix of landscape type change in Qingdao city from 1973 to 2018
表2. 1973~2018年青岛市各时段景观类型变化转移矩阵
![](//html.hanspub.org/file/7-1040093x14_hanspub.png)
Figure 2. Structure of various landscape types in Qingdao from 1973 to 2018
图2. 1973~2018年青岛市各景观类型结构
由表2可以得到,在1973~1980年景观类型转化主要为耕地转化为建设用地,其他少量建设用地转化为耕地和水域。
在1980~1990年景观类型转化主要为耕地转化为建设用地,其他少量草地转化为耕地、水域转化为耕地、林地转化为草地。
在1990~2000年景观类型转化主要为耕地转化为建设用地,其他为少量耕地转化为水域。
在2000~2010年景观类型转化主要为建设用地转化为草地、草地转化为耕地、建设用地转化为耕地、林地转化为耕地,未利用地得到利用向其他景观类型转化。
在2010~2018年景观类型转化主要为建设用地转化为耕地、耕地转化为水域,建设用地转化为水域,其他为少量林地转化为耕地、耕地转化为草地。
以上各时段主要景观类型转移变化具体如图2所示。
从总体看来,1973~2018年青岛市主要景观类型变化为耕地向建设用地的转移,在1973~1980、2000~2010两个时段内转移面积最大,说明耕地是建设用地扩张的主要来源。其中耕地的来源主要为林地、草地,以及少部分的水域。此外,从耕地转林地和林地转耕地、耕地转草地和草地转耕地、耕地转水域和水域转耕地的对比中可以发现,主要转化类型为耕地向林地和草地的转移,耕地向林地、草地、水域的转化也出现一定幅度的增加。综合以上,可以得出1973~2018年青岛市各景观类型结构如图3所示。
3.2. 城镇用地变化分析
3.2.1. 城镇用地扩展指数
通过对实验数据分析得出表3、图4。
![](Images/Table_Tmp.jpg)
Table 3. Percentage growth and expansion index of urban and rural land in different periods
表3. 城乡用地不同时段增长百分比和扩展指数
![](//html.hanspub.org/file/7-1040093x22_hanspub.png)
Figure 4. Change in percentage of urban land use growth and expansion intensity index in different time periods
图4. 不同时间段城镇用地增长百分比与扩展强度指数变化
青岛市城镇用地在1973~1980年增长很低,甚至出现负增长现象,负增长百分比为−33.04%,扩展强度指数为−0.16,城镇用地面积出新减少,属于负扩展,这在研究的所有时间段内都是极其少见的状况。1980~1990年增长较多,增长了46.44%,扩展指数强度为0.11,属于中速扩展,在1990~2000年增长速度增长稍弱,为38.84%,扩展指数强度为0.13,为缓慢增长,在2000年至2010年内增长面积最多,并且增长速度最快,由519498924.00 增长到了1112193431.00平方米,增长了114.08%,扩展强度指数为0.53,为高速扩展,在2010~2019年经过城镇用地高速发展之后,增长速度放缓,无论是增长百分比还是扩展强度指数都到了瓶颈,增长百分比为2.46%,扩展强度指数为0.03。由城镇用地的变化可以看出来,城镇用地的增长,明显与政府的宏观调控与政策相挂钩。尤其是2010~2018年这段时期内,因为国家重视生态城市的建设,并且顺应经济发展的需要,城镇用地的增长较为保守。
3.2.2. 城镇用地紧凑度与分维度指数
从空间上来看,青岛城镇用地扩张的紧凑度在因为地形以及海岸线等因素,导致明显偏低,由表4紧凑度来看,2000年的青岛紧凑度最高,1973年城镇发展用地的数值最低,为0.0087。除去最高的2000年的城镇发展紧凑度与最低值1973年的城镇用地紧凑度来看,其余时期内的城镇用地发展的紧凑度呈现波浪形起伏但是稳定在0.09左右,这表明了青岛市的城镇用地紧凑度在其余时期内有着稳定的发展规律。从分维指数来看,青岛市城镇用地的分维指数一直低于1.6,说明城镇用地建设的增长具有比较强的紧凑性。图5、图6,从整体来说,青岛的城镇用地建设因为地形以及海岸线等客观因素呈现较为分散的现象,但随着时间的增长分散程度趋于稳定,并且,其扩展形态和边界一直较为规则。这说明,青岛市的此次城镇用地发展尽管受地形所限,但一直都有所规划,并非是盲目的发展。
![](Images/Table_Tmp.jpg)
Table 4. Index of changes in urban land compactness and sub-dimension in different periods
表4. 不同时期下城镇用地紧凑度与分维度指数变化
![](//html.hanspub.org/file/7-1040093x23_hanspub.png)
Figure 5. Changes in urban land use compactness index in different time periods
图5. 不同时间段城镇用地紧凑度指数变化
3.3. 城镇用地发展影响因素分析
3.3.1. 自然地理因素
人类居住空间的组织形式开始向邻里模式转变,人控制环境的范围和认知能力会限定居住空间规模 [29]。青岛是海滨丘陵城市,地势总体为东高西低,南北两侧隆起,中间显示低凹的状态。其中,山地约为青岛市总面积的15.5%,丘陵占2.1%,平原占37.7%,洼地占21.7%。青岛市海岸的种类包括岬湾相间的山基岩岸、山地港湾泥质粉砂岸及基岩砂砾质海岸等3种基本类型。这基本上决定了青岛市城镇建设的位置以及大体的走向。
![](//html.hanspub.org/file/7-1040093x24_hanspub.png)
Figure 6. The change of urban land use index by dimension in different time periods
图6. 不同时间段城镇用地分维度指数变化
3.3.2. 经济发展条件
经济的发展是城镇用地发展的基础,其增长快慢决定着城区的扩展速度以及空间格局的演变,可以说,经济发展是推动城区发展的主要推动力。由上述实验数据观察到,在2000~2010年间城镇建设用地发展最快,选取此时段青岛市GDP发展指标与全国同一时期的GDP发展指标先对比(表5、图7),发现此时青岛市的经济发展曲线高于全国,说明了经济因素对于城镇用地的发展有着强力的推动作用。
![](Images/Table_Tmp.jpg)
Table 5. Comparison of economic development from 2001 to 2010
表5. 2001~2010年经济发展对比
4. 结论
在青岛市土地利用变化中,本文对青岛市1973~2018年各时段的土地利用类型从结构比例、转移类型选取指标进行了分析,得到了在1973~2018年各时段土地利用覆被变化、转移、景观格局和城镇利用土地的时空变化趋势,得到以下结论。
![](//html.hanspub.org/file/7-1040093x25_hanspub.png)
Figure 7. Changes in Qingdao’s economic development indicators from 2000 to 2010
图7. 2000~2010青岛市经济发展指标变化
4.1. 青岛市的景观格局变化
青岛市的景观格局的变化过程是动态的。在不同的时期内,各类景观用地的变化也不相同。由整体时间期限观察,可以发现一些基本的规律。自1978年~2018年,耕地、林地、草地以及未利用地面积总体呈现下降状态,水域以及城镇建设用地的面积一直呈现上升趋势,这说明了青岛市的自1973年到2018年为止的土地利用总体趋势主要是将耕地、林地、草地以及未利用地转化为建设用地与水域用地。在实现生态化城市建设的同时,政府需要审慎地对于各类用地的转化数量加以控制,加强关于耕地、林地以及草地的数量,同时保留出一定的未利用地,这能够更加有效率地完成生态城市的建设。
4.2. 青岛市的城镇土地利用转移
就城市景观格局变化,青岛市城镇土地利用转移趋于稳定。在城市化背景下,随着城市建设对建设用地需求的增加,其他土地利用类型向建设用地转移,但随着城镇化建设的巅峰时期过去,人们逐渐拾起了对于生态城市建设的关注,这在2010年~2018年其余土地领用转型为建设用地为历史上最低,同时期内,其余土地转化数量水域用地最高,其次为草地。可以预见的是,土地转化类型数量在生态用地的建设中可能会有所变化,但是城镇建设用地增长数量将大体保持不变。
4.3. 城镇用地紧凑度稳定
从整体景观格局变化角度分析,城镇用地紧凑度维持在一个稳定的数据上下浮动,这说明青岛市的城镇用地发展规划始终规划,尽管受限于地形因素,但是其规划由数据来看已经摸索出一些因地制宜的发展经验,城市形态趋向紧凑,城市结构向着标准化发展。而且,城镇用地紧凑度不能单独分析,而是要与外部大环境相结合。例如,城镇用地紧凑度在2000年~2010年呈现下降趋势,但是与此同时也要注意到,该时期内的青岛城镇用地处于扩张的巅峰时期,该数值的显降表明了这段时期内,由于城镇用地的大发展,改短时期内的规划显然需要注意一定的问题。同时,城乡用地分维数的降低,建设用地斑块复杂性降低,说明政府开始有意识地提高城市发展规划的科学性。政府的主导因素是促使景观类型分布区域简单的主要原因。
4.4. 自然地理条件决会对城镇用地的发展起到双向作用
经济因素推动着城镇用地的发展,而政策因素则对于经济因素有着直接的影响,还对城镇用地的发展起着导向作用。在城镇用地发展这方面观察,城镇用地在10年以后发展保守,这一方面是生态城市建设,另一方面,由城市建设的紧凑度指数可以看出来,城镇用地扩张收到了地形的限制作用。从经济方向来看,自1973~2000年,城镇用地的发展还呈现波浪形的发展趋势,但到了下一时期内,城镇用地发展呈现快速上升趋势。自2000年至2010年,根据图表可以明显看出此时期内青岛市GDP增长率曲线上升趋势超过全国GDP增长率曲线,此时的青岛市城镇用地的发展也呈快速发展状态上升状态。在政策方面来说,2010年是“十一五”计划的最后一年,青岛市土地利用结构在“十一五”期间进行了调整,新的土地利用模式初现成果,在城镇用地扩展方面,扩展速度变缓以适应新模式,其余景观用地也相应地在局部增加,同时也说明两者共同响应于政策的号召。