1. 引言
1.1. 研究目的
近年来,土地合理利用日益受到关注,土地利用也是景观标志的突出体现。随时代的进步与经济发展,城市化与工业化日益推进,为土地带来诸如土地生产力下降、河流干枯、林地减少等环境、生态问题,也带来很多负面影响,土地利用面临转型,这些种种变化必定引起景观格局变化。对不同时期的景观格局进行分析,导致其变化的因素及变化的结果,最终得出有参考价值的数据分析成果。土地利用空间具有生产、生态、生活三类功能分区,基于此分类体系,在时间方面上,分析十年间的变化,可看出区域内生产用地、生活用地以及生态用地的转型情况,即土地利用的景观格局变化,发现问题;在空间方面上,分析景观内部的异质性与差异性 [1],景观生态的脆弱性,探讨三生空间之间的转换与联系与对景观结构的影响 [2]。分析景观格局在时空上的演变,对其进行指数分析 [3],可以为区域内土地利用的转型方向提供建议,并为区域内部环境与生态的改善提供参考数据和方案依据。
1.2. 研究区概况
栖霞市隶属于山东省烟台市,位于北纬37˚05'N~37˚32'N,120˚33'E~121˚15'E之间,东西63千米,南北51千米。地处山地丘陵地区,山峦起伏,平均海拔178.71米,四季分明。年平均气温11.4℃,无霜期约211 d,年平均降水量648~846 mm。受地形影响,坡短湍急,汛期容易形成山洪,也容易形成旱灾。
1.3. 研究方法
在地理空间数据云中下载2009年10月与2019年4月的Landsat TM遥感图像,云量均小于10。采用ENVI对遥感影像进行大气校正、几何校正进行图像预处理;将预处理后的图像导入eCognition,进行样本点选取和土地利用类型划分,将该区域分为7类土地利用类型,分别为耕地、林地、园地、建设用地、居民点、水域、工矿用地,对图像进行修正后导出栅格文件。之后利用Arcgis10.2进行图像页面处理,并将两个年份的图像叠加融合,导入Excel表格制作转移矩阵,得出变化数据。再以土地利用类型图为底图,进行栖霞市“三生空间”分类,分析动态变化。最后,将土地利用类型图像导入Fragstats4.2,进行各类景观格局指标的叠加计算,得出数据,进行分析。
2. 国内外研究概况
景观生态学是从上世纪80年代开始在国内兴起,研究方法也有越来越大的进展。由于景观格局的演变因素是自然因素与人为因素的共同结果,通常,都是通过景观指数分析法或者建立空间数字模型 [4] 来对景观的动态变化进行分析研究。近年来,国内外对景观格局分析最常用到的软件为Frastats、Pach以及Apack这三类软件对景观指数进行计算分析,如今,能够反映不同优势度、破碎度、多样性的景观指数多达200种,在进行不同区域的景观格局分析时,需要酌情选取。
在我国,生态愈加脆弱,在该研究领域涉及越来越多,如赵培强、陈明霞等人对祁连山景观多样性的指数分析 [5];包慧娟集成地理信息系统与遥感,获得景观要素补丁,获取补丁的数目、周长等自身属性特征,运用景观生态学的原理与方法,计算研究所需的景观指数,对结果进行分析并揭示景观现状和动态 [6];樊彦丽、田淑芳等人对天津市滨海新区的景观格局时空演变研究等 [7]。多种多样的研究方式为景观生态学的研究贡献力量。
国外在景观格局的研究中,也有诸多贡献。Baskent等在2001年便论证了景观格局变化的驱动力,并分为自然和社会两方面进行论证 [8];Joseph Tzanopoulos等探讨了海中小岛景观的变迁与模式,主要研究了物理和人为之间的相互作用 [9];Alphan Hakan等选择研究区内的三个点并对其图像分类,用分类结果的交叉表对土地变化做定量分析,通过对早期和后期面积和形状指标的计算值对景观变化进行评估 [10]。当今年代,有更多年轻学者对景观动态变化有了更精准、更优质的分析,在研究过程与方法中总结出更多经验。
3. 土地利用变化研究
人与土地关系一直以来是长期土地发展的关键问题,分析土地利用变化是非常重要的一方面,通过土地利用覆被变化情况,可以得到利用类型变化趋势以及覆被变化趋势。从而对各类用地合理利用以及生产、生活、生态空间进行合理整治,促进土地利用的可持续发展。
3.1. 土地利用类型变化分析
利用Arcgis10.2对图像进行空间叠加分析,导入Excel表格得出土地利用类型转移矩阵,计算土地利用综合指数,对栖霞市土地利用的类型变换做出分析。
土地利用综合指数(La):以定量的方式表示土地利用根据地形分类表现出的差异变化规律显示出的综合利用程度。在此,将各类土地利用类型分为4级,并为7类土地利用类型赋值。园地与耕地为4,林地与水域为3,建设用地、居民点、工矿用地为2,未利用地为1。
其公式为:
[8]
其中,Ai表示分级指数,Ci表示利用面积百分比。
通过计算得出,栖霞市内占地面积比重最大的土地利用类型为园地,占地65%以上,另外以耕地与林地居多,其优势土地利用类型显著。在利用程度指数运算后,得出各级土地利用程度处于较优水平,但不同类型有一定差异。根据不同地形地貌,其分异的程度与特征各不相同,如:建设用地周边土地利用程度较高,而滩涂等地利用程度则较低。
图1和图2为栖霞市2008年与2018年土地利用类型状况;表1为转移矩阵所反映的地类转变情况,从中得出以下分析。
由于区域内多山地丘陵,栖霞市土地利用类型以园地和林地为主,园地面积占比约65.38林地面积占比约12。2009年至2019年十年之间,栖霞市园地减少约378.44平方千米,呈现下降趋势,下降约4.4%,而由于栖霞市颁布《水土保持工作条例》,与“国家级生态示范区建设”重点工程的实施,积极贯彻落实退耕还林措施,林地面积增加426.57平方千米,呈上涨态势,涨幅为30.45。相比之下,耕地面积基本持平,仅有微量减少;而水域面积有大幅度缩小,具体缩减78.95平方千米,下降36.56。在建设用地方面,根据“栖霞市国民经济和社会发展第十个五年计划”等。栖霞市近年注重制造业基地建设,建设用地面积有上涨趋势,用地面积增长141.07平方千米,增幅为41,居民点有小幅缩减,但总体呈平稳态势。
![](//html.hanspub.org/file/9-1150715x12_hanspub.png)
Figure 1. Land use status of Qixia city in 2009
图1. 2009年栖霞市土地利用类型
![](//html.hanspub.org/file/9-1150715x13_hanspub.png)
Figure 2. Land use status of Qixia city in 2019
图2. 2019年栖霞市土地利用类型图
![](Images/Table_Tmp.jpg)
Table 1. Land use type transfer matrix of Qixia in 2009 and 2019
表1. 2009、2019年栖霞市土地利用类型转移矩阵
在土地利用指数方面,由于近年来一些政策出台,尤其是2017年又颁布了退耕还林补助政策,更加推动了林区复种的推行,土地利用程度指数有下降趋势,其园地面积、耕地面积等都有减少。因此,通过单一的土地利用变化,我们可以看到栖霞市政策在过去十年中对土地利用变化的影响。在土地利用结构中,没有明显的变化,园地仍然是优势地类,但优势度有缩减趋势。退耕还林将稳定该地区的生态系统。
3.2. 生产、生态、生活分区利用变化分析
以土地利用类型为底图,根据其功能的不同,及土地利用的侧重点不同,从而形成生产、生活、生态相互关联与影响的整体,并利用Arcgis10.2将不同土地利用类型划分为“三生空间”分类体系,并制成出图。分析各级面积与结构变化。根据下表2可以看出,将其分为两级分类体系,第一级为生产、生活、生态用地;第二级为农业生产空间、城市生产空间、城乡生活用地、重点生态用地、生态容纳地。而下图3~6中可直观看出各级分区的空间面积变化。
通过图像显示,栖霞市十年间生产、生活、生态用地面积有显著变化,生产空间,尤其是农业生产空间有大幅度缩减,城市生产空间也有所扩张;与之相反,生态空间,尤其是重点生态用地有明显提升,生态容纳地有微小幅度缩减。其原因是烟台市生态保护红线的划出以及生态空间布局优化措施实施,栖霞市生态保护红线保护面积约278.87平方公里,占总面积的13.8%,对于生态空间优化方面,保护重点生态区域,维护生态安全与可持续发展,对于水域,重点保护清洋河、白洋河、漩河、黄水河的河流网体系,优化空间布局,调整生态景观结构,在发展经济与推进城市化的同时,与生态空间的保护双管齐下,通过各类别管控措施,减缓了生境方面的破碎化,促进烟台市内土地空间向“双核、一带、多轴”的景观格局发展。
![](Images/Table_Tmp.jpg)
Table 2. Classification table of production-life-ecology of qixia city
表2. 栖霞市三生空间分类表格
![](//html.hanspub.org/file/9-1150715x14_hanspub.png)
Figure 3. Land division map of Qixia city in 2009
图3. 2009年栖霞市土地一级分区图
4. 景观格局信息提取与分析
景观格局指数的提取与分析,是分析景观格局动态变化的方法之一,也是必要环节。根据不同地域选择不同景观格局指数,根据时空变化显示的不同特征,能分析出区域的优势和劣势,从而对未来的趋势进行预判,并采取措施修补景观劣势,突出景观优势。
4.1. 景观格局指数提取
景观格局指数可以清晰体现某区域内景观的结构与景观的差异性,可提取出高浓度景观信息,因此,
![](//html.hanspub.org/file/9-1150715x15_hanspub.png)
Figure 4. Land division map of Qixia city in 2019
图4. 2019年栖霞市土地一级分区图
![](//html.hanspub.org/file/9-1150715x16_hanspub.png)
Figure 5. The second-level zoning map of Qixia city in 2009
图5. 2009年栖霞市土地二级分区图
![](//html.hanspub.org/file/9-1150715x17_hanspub.png)
Figure 6. The second-level zoning map of Qixia city in 2019
图6. 2019年栖霞市土地二级分区图
在此我们选择利用景观格局指数分析法来做出分析。本文选取了斑块数量(NP)、类型斑块面积(CA)、斑块所占面积比例(PLAND)、散布与并列指数(IJI)、香农多样性指数(SHDI)、香农均匀度指数(SHEI)、蔓延度指数(CONTAG)这7种景观格局指数指标。
4.2. 景观格局指数的计算与分析
分别将栖霞市2009年、2019年土地利用类型图,以tif格式导入Fragstats4.2,选择对应的指标进行计算,从而得出以下图表3。
(一) 景观结构变化分析
![](Images/Table_Tmp.jpg)
Table 3. Landscape structure index of Qixia city in 2009 and 2019
表3. 栖霞市2009年、2019年景观结构指数
注:表内CA表示类型斑块面积,PLAND表示斑块面积所占比例。
2009年与2019年栖霞市景观结构指数,如表2所示,10年之间各类景观用地都发生了变化。首先,总体斑块数量由3935转变为3661,斑块数量有下降趋势,说明该地区景观结构有紧密性趋势,破碎化程度降低。分开来看,耕地、园地斑块面积有所下降,由于栖霞市地处山地丘陵,种植果树面积较大,因此耕地与园地是该区域内面积比中最大的部分,10年之间该地类斑块所占比例有降低趋势,这与生态保护政策密切相关;由此,园地面积比例下降3个百分点,而林地斑块面积比例上涨了3.17%,面积上涨幅度较大。水域方面,水域面积占地比例有下降趋势,根据图表可以看出,10年之间有河流干枯情况,干枯河流成为滩涂或荒地,使原本的水域面积减小,这也是生态破坏的表现,如何将干涸水域再利用也是生态保护的重点所在。另外,建设用地板块面积上涨幅度较大,这与经济飞速发展不可分割;城市高层的扩建,使居民点占地面积缩小,土地利用程度变大,因而居民点斑块面积有微幅缩小。工矿用地与未利用地变化面积不大,保持平稳。
(二) 景观类型水平分析
首先是耕地类型水平分析,由表2可得出,耕地斑块面积由23,645.61 km2下降为23,360.85 km2,下降幅度非常小,斑块面积比例下降0.2%左右,总体未发生太大变化。由于耕地红线保护,耕地不会急剧减少。景观内耕地分布破碎程度加大,说明人为因素影响下,即使耕地面积无明显变化,耕地的密集程度降低。
第二为园地类型水平分析,园地斑块总面积由134,753.7 km2减少为129,078.3 km2,相对于耕地,园地的变化幅度非常大,下降5675 km2,变化速度较快,为所有地类中变化幅度最大的地类,但园地始终为该地区的主导地类。总体看来,农业生产用地方面,耕地与园地有缩减趋势,但不会有大幅度减少,且会有回暖,根据2017年栖霞市规划完善下达永久基本农田任务为72.5730万亩,截止到2014年,土地整理下达的补充耕地面积仅完成不到20%,必须坚守住基本农田红线不动摇,并加快耕地补充。
另外是建设用地和工矿用地,此类城市生产用地方面,有扩大趋势,面积比例总体上涨1.1%,根据栖霞市2018年《政府工作报告》,栖霞市将加快总投资,建设重点项目,发展实体经济,坚持制造业和服务业双轮驱动,促进城市化率,因此,今后栖霞市建设用地规模与工矿用地规模将不断扩大。
城乡生活用地,居民点方面,有缩减态势,面积比例减少0.9%,栖霞市的土地节约集约利用程度加强,重建高层居民点,节约土地利用面积,该类斑块面积减少是必然。
生态用地方面,根据表格得出,林地十年增长6397.26 km2,为增速最快的土地利用类型,斑块面积比例增大,说明其优势型越发明显,在退耕还林措施的实施情况下,该种地类将有持续上升趋势。水域方面,虽然有减少趋势,但是变化不大,维护水域生态面积是维护生物多样性的重要一环。
(三) 景观多样性变化分析
![](Images/Table_Tmp.jpg)
Table 4. Landscape diversity index of Qixia in 2009 and 2019
表4. 栖霞市2009年、2019年景观多样性指数
注:表内IJI表示散布与并列指数,SHDI表示香农多样性指数,SHEI表示香农均匀度指数,CONTAG表示蔓延度指数。
在2009年与2019年景观多样性指数表格中(表4),可以得出:首先,IJI指数处于适中数值,说明该地有明显异质性,这源于该地位于山地丘陵,确实受垂直地带性影响;香农多样性指数减小幅度约0.04,表明景观斑块有向大型斑块发展趋势,破碎程度降低,异质性亦有减弱态势;香农均匀度指数下降约0.02,变化很小,维持稳定状态。香农多样性指数与香农均匀度指数整体变化幅度较小,原因是退耕还林战略的实施,林地面积增加迅速,优势景观优势度增强,优势景观也加强其控制能力。蔓延度指标显示的是区域内部景观的聚集程度或是连接度,该地蔓延度指数上升了约2%,可看出该地各景观之间的连通性确实有加强,优势景观园地、林地起较大的控制作用。
5. 景观格局变化驱动力分析
分析栖霞市景观格局变化的驱动力可以更好把握区域内的动态变化因素,从而分析若以此种方向发展是否有利于人类生产生活以及生态环境保护和可持续发展。若不利于,则可更换思路,转变发展方向,为栖霞市的土地规划提供有意义的参考数据。
每个区域的发展内因都出于两个大方面,一为自然,二为社会。自然因子包括土壤、坡度、气候、地形等因素;社会因子包括一、二、三产业产值,GDP、城镇化水平、人口密度等因素。较为稳定的因素为自然因素,对景观造成的为长期且持续影响;而人为因素则较为活跃,发挥的为短期且剧烈的作用。
5.1. 自然因素
本文选取地形、气温、降水三个方面进行分析。地形方面,栖霞市为山区丘陵地形,素有“六山一水三分田”之说,由于地形高低起伏不平,高程与坡度对于景观的分布与土地利用造成很大制约。而在气温方面,据山东省栖霞市气象局数据显示,10年来,年平均气温上涨0.4℃,≥10℃积温日期由226 d缩减为224 d,霜冻日期加长,也影响了园地内农作物产量。降水方面,年降水量的变率有所增大,年降水量有下降趋势,年降水量小于400 mm的年份增加,因此栖霞市发生干旱的概率增加,影响农业生产用地产出与生态用地保护。
5.2. 社会因素
经济驱动因子:经济发展永远是催化土地利用类型改变与景观格局演变的根本推动力,它对于栖霞市景观格局的动态变化的作用最为广泛、深远、长久。
![](Images/Table_Tmp.jpg)
Table 5. GDP and population of Qixia in 2009 and 2019
表5. 2009、2019年栖霞市GDP与人口数值表
根据表5可以看出,2009~2019期间,栖霞市人口增长速度较快,GDP在十年内翻倍,经济发展速度迅速,国民生产总值增长必然是城市化推进的成果,体现在景观格局方面,土地利用中,建设用地出现飞速增长态势,这与园地的退化速率呈正相关,经济的增长加速了农业生产用地退化,因而也引起生态用地退化。
人口驱动因子:人口也是影响景观格局动态变化的重要因素之一,人口增长对住房、交通、工作等方面有了更多要求,因此会占据更多土地开发为公共基础设施或社会服务用地,这必然导致生产、生态用地减少,建设用地、生活用地增加,人口猛增其实对生态用地存在很大威胁。
政策驱动因子:根据栖霞市《土地利用总体规划》等一系列重要文件显示,栖霞市将会更加重视耕地保护与生态环境保护,政策颁布与强制性实施有利于景观格局在人为控制下,向可持续方向发展。
6. 结论
本文在ENVI、Arcgis10.2、eCognition、Frastats4.2的支持下,绘制了栖霞市土地利用类型图、土地利用类型分区图,并对其景观格局指数进行分析和驱动力分析,得出以下结论:
2009~2019年,土地利用类型在农业生产用地方面,转变较大,尤其是主导类型园地方面,减少迅速;其余有小幅变化,城市生产用地增加,生态用地减少,但总体趋势维稳。在景观格局指数方面,十年间斑块总体数量减少,斑块面积增大,说明破碎化程度降低,整体性加强;从多样性方面来看,香农多样性指数降低,蔓延度指数升高,异质性加强,连通性加强。
在自然因素与社会因素的共同作用下,栖霞市的景观格局发生较大变化,自然因素是推动变化的基础,而社会因素则是推动景观演变的催化剂,是推动农业生产用地、生态用地退化的主要驱动力,经济增长速度越快,越容易导致区域内生态用地畸形,因此,需要国家政策来强制调控。对于该市景观的演变过程与结果,有关部门应及时落实,调整措施进行景观格局的调整。
参考文献