安徽省六安市土地利用变化及景观格局梯度效应研究.pdf

引文格式：王诗琪,周振宏,智菲,等.安徽省六安市土地利用变化及景观格局梯度效应研究 J.云南农业大学学报（社会科学）,2023,17(5):154162.DOI:10.12371/j.ynau(s).202304041.安徽省六安市土地利用变化及景观格局梯度效应研究王诗琪，周振宏*，智菲，赵铭（安徽农业大学林学与园林学院，安徽合肥230036）摘要：为揭示安徽省六安市土地利用变化特征以及不同地形分级上景观格局的演化规律，计算地形位指数、分布指数和景观指数，结合土地利用转移矩阵、动态度分析和地学信息图谱分析法，对研究区 20002020 年间 3 期遥感影像和数字高程数据进行处理。结果表明：六安市用地类型以耕地、林地和草地为主，耕地的流转面积最大；耕地和建设用地在低等级地形区呈优势分布；稳定型是景观变化中的主导类型；在 25 级地形区，CONTAG、SHEI 和 SHDI 出现较明显的随年份演变差异；2020 年，研究区景观丰度有所增加。在六安市未来规划中，应结合地形分布构建节约集约用地，协调功能用地布局；以生态为本，向西拓展绿色产业；延伸核心区功能范围，疏散人口，保护环境。关键词：土地利用；景观格局；地形梯度效应；地学信息图谱中图分类号：F301.24文献标识码：A文章编号：1004390X(2023)05015409Study on the Land Use Change and Landscape PatternGradient Effect in Lu an,Anhui ProvinceWANGShiqi，ZHOUZhenhong，ZHIFei，ZHAOMing（SchoolofForestryandLandscapeArchitecture,AnhuiAgriculturalUniversity,Hefei230036,China）Abstract:InordertorevealthecharacteristicsoflandusechangeinLuancityandtheevolutionlawoflandscapepatternindifferenttopographicgrades,thetopographicpositionindex,distributionin-dexandlandscapeindexwerecalculated,andthethreeremotesensingimagesanddigitalelevationdataofthestudyareafrom2000to2020wereprocessedbycombininglandusetransfermatrix,dy-namicattitudeanalysisandgeoscienceinformationmapanalysis.Theresultsshowedthat,themaintypesoflanduseinLuanwerecultivatedland,woodlandandgrassland,andthecirculationareaofcultivatedlandwasthelargest.Cultivatedlandandconstructionlandweredominantinlow-gradeter-rainareas;Stabletypewasthedominanttypeinlandscapechange;Ingrade25terrainarea,CON-TAG,SHEIandSHDIhadobviousdifferenceswiththeevolutionofyears;In2020,thelandscapeabundanceofthestudyareawouldincrease.InthefutureplanningofLuan,weshouldcombinethetopographicdistributiontobuildeconomicalandintensivelandandcoordinatethelayoutoffunction-收稿日期：20230418修回日期：20230509基金项目：国家自然科学基金青年科学基金项目“基于结构生态功能耦合关系的城镇绿色空间构建研究以合肥市为例”（41301650）。作者简介：王诗琪（1998），女，安徽合肥人，硕士研究生，主要从事土地利用、景观格局方面研究。*通信作者：周振宏（1979），男，安徽望江人，讲师，主要从事城乡规划、土地利用/覆盖变化和景观可持续科学研究。云南农业大学学报（社会科学），2023，17（5）：154162http:/JournalofYunnanAgriculturalUniversity（SocialScience）E-mail:alland.Basedonecology,expandgreenindustrytothewest;Extendthefunctionalscopeofthecorearea,evacuatethepopulationandprotecttheenvironment.Keywords:landuse；landscapepattern；terraingradienteffects；geoscienceinfographic土地利用及其景观格局分布与城市的发展有着紧密联系。近年来，随着我国城市化进程的加快，城市规模不断向外扩张，导致土地资源紧缺、生态环境恶化，如何科学合理地管理和规划土地资源成为未来社会发展和生态保护的中心问题14。为此，许多国内学者选取部分典型城市，针对其以往土地演变规律进行总结，抑或是通过构建模型预测其未来用地演变方向，对于土地利用景观的驱动因子研究也有了一定成果57。地形对于人类生产生活具有限制性约束作用，而人为因素作为生态环境质量的影响要素之一，对于用地景观格局的分布也具有一定的干扰效应，目前已有研究佐证了这一观点89。此前，多数学者从单一地形因子（如高程、坡度等）出发，研究其与景观格局的响应机制，也有部分学者在前人基础上选取地形位指数作为综合性地形因子的表征，探究地形对景观格局的影响作用1013，削弱了早期研究的局限性，但其结论仍囿于具体时相内的景观分异特征，而对于一个连续时段内景观格局的动态变化过程与地形的关联机制研究尚且不多。六安市位于安徽省西部，地处江淮之间，东邻省城合肥，是安徽省省会经济圈副中心城市14。六安作为安徽省的区域交通枢纽，其境内梯度差异显著，以此为研究区，加入地学信息图谱将景观格局演变路径可视化，再叠加地形位梯度，探究地形因子对用地景观分布及其变化过程的作用机理，对于制定六安市未来土地利用规划策略具有一定参考价值。一、研究区概况六安市（31013240N，1152011714E），安徽省辖地级市，下辖3 区4 县，总面积15451.2km2，地处江淮，东衔吴越，西领荆楚，北接中原，属北亚热带向暖温带转换的过渡带，四季分明，气候温和，雨量充沛。地势上西南高峻，东北低平，呈梯形分布1516。二、材料与方法（一）数据来源及预处理本研究所用数据包括地理空间数据云网站（https:/ 和 2020 年 3 期土地利用数据，并参照该系统分类标准，将研究区用地类型划分为耕地、林地、草地、水体、建设用地和未利用地6 种类型。（二）研究方法1.地形位指数和分布指数本研究中选用地形位指数作为评价地形的指标，计算公式如下17:T=lg(H/H+1)(S/S+1)（1）HS式（1）中，T 为地形位指数，H、S 表示研究区任一栅格的高程和坡度，、表示研究区任一栅格高程和坡度的平均值。高程和坡度越大，其地形位指数越大，反之越小18。通过分布指数能够展示出研究单元内各景观类型在空间上的分布差异，及其优势分布区间，计算公式如下18:P=(Mtn/Mt)/(Mn/M)（2）式（2）中，P 为分布指数，Mtn指代 t 类景观在 n 级地形所占面积，Mt为 t 类景观总面积，Mn为 n 地形区总面积，M 为研究区总面积。P1，分布表示该类景观在 n 级地形上呈劣势分布和优势分布。2.土地利用转移矩阵土地利用转移矩阵能够直观展示出研究区内始末各土地类型的转变趋势，以及不同转移方向的具体转移数量5。其公式如下5，19：Spq=S11S12S13S1mS21S22S23S2mS31S32S33S3mSm1Sm2Sm3Smm（3）第5期王诗琪，等：安徽省六安市土地利用变化及景观格局梯度效应研究155Spq式（3）中：p 与 q 分别表示研究始末的土地利用类型；m 表示土地利用类型总数；为研究期内第 p 类向第 q 类转化的总面积。3.土地利用动态度分析通过构建土地利用动态度模型，可以展示出研究时段内各种土地利用类型变化的速度和幅度2021。第一，单一土地利用动态度能够反映某个时段内研究区某用地类型在数量、空间及组合方式上的变化，计算公式为：K=UbUaUa1T100%（4）UbUa式（4）中：K 为研究时段内某一土地利用类型动态度；和分别为研究末期和初期某一土地利用类型面积,T 为研究时段，本研究中 T=10。第二，综合土地利用动态度能够反映在一定时段内研究区整体土地利用类型的变化速率，可以从全域分析土地利用的空间变化过程，计算公式为：LC=ni=1LUij2ni=1LUi1T100%（5）LUijLUi式（5）中：为第 i 类土地利用类型转为非 i 类土地利用类型面积的绝对值，为前期土地利用类型的面积，T 为 2 期间隔的时间。4.地学信息图谱地学信息图谱是一种地理时空分析方法，它能够表述区域自然过程与社会经济发展的时态演进和空间分异情况，将地学信息图谱应用于土地利用景观格局的研究中，生成景观格局信息图谱，可以对研究区用地景观的空间格局进行系统性的描述2225。三、结果与分析（一）土地利用时空变化1.土地利用类型总体特征通过 ArcGIS10.5 软件的掩膜提取功能，将安徽省土地利用数据与六安市行政边界进行叠加，得到研究区各年份土地利用类型图（图 1）。六安市土地利用类型首先以耕地、林地和草地为主，在整个研究期内，三者面积之和占研究区总面积的 90%左右；其次为建设用地和水体，在研究期20 年间，二者的平均面积占比分别为 5.96%、4.12%；未利用地占比极小，且随年份增长面积略有增长。20002020 年间六安市土地利用转移矩阵如2000年02550100 km耕地林地草地水体建设用地未利用地2010年02550100 km耕地林地草地水体建设用地未利用地2020年02550100 km耕地林地草地水体建设用地未利用地图 1 六安市各年份土地利用类型图156云南农业大学学报第17卷表 1 所示。从变化类型来看，20002010 年间耕地作为优势地类，其流转面积最大，除自转化外，主要转出类型表现为耕地转化为建设用地和林地，这两种变化类型面积之和占耕地净流出面积的 83.53%；耕地的主要转入来源包括林地和建设用地，二者转入面积之和占净流入面积的79.28%。总的来说，耕地地类的转出大于转入，即耕地面积在减小，这一现象与城镇化发展和退耕还林政策的实施密不可分。20102020 年间耕地流转趋势与前一时段保持一致。林地、草地和水体在 20002020 年间的用地面积基本维持稳定，说明其流入和转出基本保持平衡。建设用地和未利用地在研究期内用地面积均出现增长，前者以耕地流入为绝对主导，其次由林地转化而来，与城镇化建设的无序扩张有关；后者在 20002010 年间流向林地的面积大于自林地流转而来的面积，2010 年的未利用地面积主要由草地的流入得以维持，20102020 年间以林地未利用地为主，辅以耕地、水体和建设用地的少量转入，使得 2020 年的土地利用面积较之前期有所增长。2.土地利用类型动态变化将整个研究期分为 20002010 年和 20102020 年 2个等量时段，分别计算六安市 2个时段的土地利用动态度（表 2）。由表 2 可知，从综合土地利用动态度来看，后一时段近乎为前一时段的 2 倍，说明后十年土地利用类型的变化幅度翻倍增长，比之前期有较大改变。从单一土地利用动态度来看，未利用地由于用地的流入和转出达到平衡，在 20002010 年间面积无变动；20102020 年间其面积增加 2.5704km2，仅占研究区总面积的 0.02%，但因其初始基数过小，从增长幅度上来说变化巨大。此外，2 个时段内均以建设用地的单一动态度为高值，且呈现出上升趋势，说明20 年间六安市境内的建设用地表现为显著扩张趋势。其次为水体，同样表现为扩张趋势，但幅度小于建设用地。20002010 年间，耕地和林地的面积分别减少 83.82km2和 9.27km2，单一土地利用动态度分别为0.12%和0.02%，20102020 年表 1 六安市土地利用转移矩阵单位：km2时段用地类型耕地林地草地水体建设用地未利用地20002010年耕地6577.074991.305015.577235.6445168.41070.0000林地93.64235039.316051.48726.62497.83540.0027草地13.115752.74361823.40181.52280.56160.0054水体33.95165.26412.6109583.56092.27430.0000建设用地86.40900.99720.16741.1052722.36250.0000未利用地0.00000.00810.00000.00000.00000.038720102020年耕地6434.321499.078314.322644.9307211.23000.3105林地95.89054998.189658.12655.721329.97451.7316草地16.437655.14481817.67422.21041.77750.0000水体21.18156.44221.6614592.16406.84450.1647建设用地100.02601.73070.32312.9142796.07700.3735未利用地0.00000.00270.00720.00000.00000.0369表 2 六安市土地利用动态度统计单位：%时段20002010年20102020年单一土地利用动态度耕地0.120.20林地0.020.06草地0.010.01水体0.010.31建设用地1.111.60未利用地0.00549.23综合土地利用动态度0.060.11第5期王诗琪，等：安徽省六安市土地利用变化及景观格局梯度效应研究157间二者面积持续下降且降幅增大，单一动态度变为0.20%和0.06%。研究期内，草地面积先增大后减小，故其单一动态度在 2 个时段内分别为0.01%和0.01%。（二）基于地形梯度效应的景观格局时空变化1.各景观类型在地形梯度上的分布特征通过 ArcGIS10.5 的区域分析工具得到六安市地形位指数，用自然段点法分为 10 级，并界定为低（13 级）、中低（46 级）、中高（79 级）和高等级（10 级）4 种地形位梯度，得到六安市地形梯度分布如图 2。六安市以低等级地形为主，其分布范围超过了全境总面积的一半，并以2 级地形区为甚，主要分布在研究区的东北部区域；410 级地形区主要分布在西南区域，以金寨县、霍山县和舒城县西南角为主。计算六安市各种景观类型在不同地形梯度上的分布指数，其各地类分布指数值变化情况如图 3。从研究时序来看，除未利用地外，其余地类在各年份于不同梯度上的分布情况变化不大。研究期间，未利用地的分布区域由前期的 26 级地形区拓展覆盖到 1 级区域，随者年份推进，其优势分布区间也向着低等级地形位演变。由用地类型转换可知，2010 年的未利用地由林地和草地转化而来，2020 年由林地、耕地、水体和建设用地的转入，林地和草地多分布于中低级地形位，其他地类的转入拓宽了未利用地分布的梯度区间。从优势区间来看，耕地和建设用地在低等级地形区呈优势分布，表明人类的生产生活空间多集聚于此；12 级地形区为水体的优势分布区间，与六安市境内水资源分布实况相吻合；林地和草地的优势区间覆盖最广，前者包括中低和中高等级地形位，后者覆盖 510 级地形区，说明林地和草地对地形变化的适应性最强。2.景观类型变化图谱的地形梯度借助 ArcGIS10.5 软件获取研究区的景观类型变化图谱（图 4），即将研究期间六安市景观类型的演化过程可视化。将图谱变化类型分为如下低等级地形梯度中高等级地形梯度中低等级地形梯度高等级地形梯度02040 km图 2 六安市地形梯度分布图地形梯度分布指数2000年6543210109876543212010年地形梯度分布指数6543210109876543212020年地形梯度分布指数654321010987654321耕地林地草地水体建设用地未利用地图 3 六安市不同景观类型在地形梯度上的分布指数变化图158云南农业大学学报第17卷5 种：稳定型（研究区景观类型无变化）、前期变化型（研究区景观类型在 20002010 年间发生变化，20102020 年间无变化）、后期变化型（研究区景观类型在 20002010 年间无变化，20102020 年间发生变化）、波动变化型（研究区景观类型在 20002010 年和 20102020 年间均发生变化，且 2000 年和 2020 年景观类型相同）、持续变化型（研究区 3 个时期景观类型均不相同）。提取六安市景观信息图谱的各变化类型，并对其面积占比和主要变化类型情况进行统计（表 3），再通过与地形位梯度叠加获得研究区各类型图谱在地形梯度上的分布情况（图 5）。稳定型作为六安市景观变化中的主导类型，占总面积的 95.28%，其优势区间为 510 级地形区，对于地形变化的适应性更强；前期和后期变化型图谱分布分别占总面积的 0.91%和 1.49%，二者的优势区间均为低等级地形位，且主导变化类型分别为耕地建设用地建设用地和耕地耕地建设用地，这主要是由于六安市城镇化建设以及土地开发复垦政策的实施所导致。波动型在 35 级地形区呈优势分布，主要变化类型以林地耕地林地和耕地建设用地耕地为主，与前期研究得出的耕地和建设用地主要分布在低等级地形区，人类生产生活空间也多集聚于此的结论相对应。持续型图谱分布面积占比极小，其优势区间为 1 级和表 3 六安市景观类型变化图谱统计景观图谱面积/km2面积比率/%主要变化图谱类型面积/km2面积比率/%稳定型14589.705095.28耕地耕地耕地6405.2541.83林地林地林地4986.8232.57前期变化型138.60450.91耕地建设用地建设用地82.63410.54耕地水体水体19.03260.12水体耕地耕地16.89750.11后期变化型228.00941.49耕地耕地建设用地126.5830.83耕地耕地水体27.85120.18林地林地建设用地26.87580.18波动变化型351.06562.29林地耕地林地85.40550.56耕地建设用地耕地84.31350.55耕地林地耕地82.36620.54持续变化型5.84160.04耕地林地建设用地2.11810.01耕地水体建设用地1.232270.01林地耕地建设用地1.134040.0102040 km稳定型前期变化型后期变化型波动变化型持续变化型图 4 六安市景观类型变化图谱稳定型前期变化型后期变化型波动变化型持续变化型01234567891011地形梯度分布指数2.01.81.61.41.21.00.80.60.40.2图 5 各类型图谱在地形梯度上的分布指数第5期王诗琪，等：安徽省六安市土地利用变化及景观格局梯度效应研究15969 级，主要以耕地、林地和建设用地之间的相互转化为主。（三）景观格局的地形梯度效应基于景观斑块的密度、形状、多样性和聚散性等角度，选取：斑块密度 PD、最大斑块指数LPI、边缘密度 ED、景观形状指数 LSI、蔓延度CONTAG、丰富度 PR、香浓均匀度指数 SHEI 和香浓多样性指数 SHDI 共 8 个指标，对六安市不同地形位等级上景观格局指数的变化情况进行表达（图 6）。从时序上看，各景观指数随梯度的变化趋势大致相同。在 25 级地形区，CONTAG、SHEI和 SHDI 的景观指数值出现较明显的随年份演变差异，蔓延度指数呈现逐年下降趋势，香浓均匀度和香浓多样性指数逐年上升，这是由于近年来人类活动愈加频繁导致景观斑块的连通性下降，破碎化程度高，景观优势度降低，异质性增强。2020 年，丰富度指数的峰值覆盖范围由原本的地形位梯度地形位梯度分布指数分布指数分布指数分布指数分布指数分布指数分布指数分布指数PDEDCONTAGSHEILPILSIPRSHDI12010080604020025200020102020201510507006005004003002001000222018161412107570656055506.05.55.04.54.03.51.41.31.21.11.00.90.80.70.60.50.80.70.60.50.41234567891012345678910图 6 景观格局指数随梯度变化图160云南农业大学学报第17卷26 级地形区拓展到整个低中低等级地形位，表明景观丰度有所增加。35 级地形区的斑块密度明显高于其他地形位，表明此区域破碎化程度最高。最大斑块指数的显著高值分布在 1、2 级地形区，耕地作为六安市境内的优势景观多分布于此。边缘密度指数值随梯度呈波动式起伏，其低值主要分布在 9、10 级地形区上，表明高梯度地形区由于开发难度大，景观更加完整。景观形状指数在 12 级地形区内骤升，此后小起伏持续下降至最小值，说明2 级地形区的景观形状最复杂，梯度越高人类活动干扰越小、开发难度越大，景观形状越完整。四、结论与讨论本文采用六安市 2000、2010 和 2020 年 3 期土地利用数据和 DEM 高程数据，利用 ArcGIS10.5和 Fragstats4.2 软件探究研究区境内土地利用的变化规律和景观格局的空间分异特征，得出结论如下。（1）六安市用地类型以耕地、林地和草地为主。研究期间，耕地的流转面积最大，其中以耕地与林地和建设用地的相互转化为主，耕地地类的转出大于转入，即耕地面积持续减小。（2）耕地、建设用地和水体在低等级地形区呈优势分布，未利用地的优势分布区间向着低等级地形位演变，林地和草地的优势区间覆盖最广，对地形变化的适应性最强。（3）稳定型作为六安市景观变化中的主导类型，其优势区间为 510 级地形区，前期和后期变化型图谱优势区间均为低等级地形位，波动型在 35 级地形区呈优势分布，持续型优势区间为1 级和 69 级。（4）在 25 级地形区，CONTAG、SHEI和 SHDI 的景观指数值出现较明显的随年份演变差异，蔓延度指数逐年下降，香浓均匀度和香浓多样性指数持续上升，2020 年，景观丰度有所增加。由研究可知，六安市的景观格局分布明显受到地形因子驱动，因此针对研究区未来土地规划提出如下建议2629：合理优化六安市土地资源配置，应结合地形分布特征科学构建节约集约用地格局，协调布局各类功能用地，保障城市发展；以生态网络为本，向西拓展绿色产业；自城市中心向四周延展核心区功能覆盖范围，防止人口过度聚集导致区域生态环境持续恶化。参考文献1崔文兵,侯国博,陈万基,等.基于多情景模拟的伊犁河谷土地利用景观格局分析J/OL.水利水电技术（中英文）:11720230421.http:/ TM 影像的环巢湖土地利用景观格局驱动力研究J/OL.云南农业大学学报（社会科学）:1820230421.http:/ 35 年变化J.生态科学,2022,41（6）:41.DOI:10.14108/ki.1008-8873.2022.06.006.9刘彦随,璩路路.近 30 年三峡库区用地格局变化与人地系统演化J.长江流域资源与环境,2022,31（8）:1664.DOI:10.11870/cjlyzyyhj202208002.10贾磊,姚顺波,邓元杰,等.20002020 年陕西秦巴山区生境质量时空演变及其地形梯度效应J.长江流域资源与环境,2022,31（2）:398.DOI:10.11870/cjlyzyyhj202202014.11周晓艳,何依依,黄欣,等.三峡库区生境质量的地形梯度效应及对土地利用变化的响应J.农业工程学报,2021,37（11）:259.DOI:10.11975/j.issn.第5期王诗琪，等：安徽省六安市土地利用变化及景观格局梯度效应研究1611002-6819.2021.11.029.12龚熊波,杨波,刘雨先,等.湘西州山区土地利用类型的地形梯度效应与景观格局分析J.生态科学,2020,39（4）:137.DOI:10.14108/ki.1008-8873.2020.04.018.13王启名,杨昆,许泉立,等.滇池流域土地利用变化图谱的地形梯度效应J.水土保持研究,2018,25（3）:237.DOI:10.13869/ki.rswc.2018.03.033.14王瞒遇.国家园林城市六安市生活污泥处理处置现状及土地利用研究J.安徽化工,2022,48（3）:110.DOI:10.3969/j.issn.1008-553X.2022.03.031.15王庭,余佳华,袁金羽.安徽六安市提升融入淮河生态经济带发展战略与路径研究J.淮南师范学 院 学 报,2022,24（6）:47.DOI:10.3969/j.issn.1009-9530.2022.06.011.16六安市人民政府办公室关于印发六安市物流业高质量发展三年行动计划（20222024）的通知R.六安市人民政府公报,2022:46.17曾真,艾婧文,陈凌艳,等.山地城市生境质量时空演变及地形梯度效应研究:以三明市区为例J.南京师大学报（自然科学版）,2023,46（1）:100.DOI:10.3969/j.issn.1001-4616.2023.01.014.18陈铸,黄雅冰,朱志鹏,等.基于地形梯度特征的福州市景观格局演变J.应用生态学报,2018,29（12）:4135.DOI:10.13287/j.1001-9332.201812.010.19李佩君,左德鹏,徐宗学,等.基于地形梯度的雅鲁藏布江流域土地利用及景观格局分析J.山地学报,2022,40（1）:136.DOI:10.16089/ki.1008-2786.000661.20吴海中,胡刚,田晓四,等.安徽省“三生”土地利用动态变化与生态响应J/OL.贵州师范大学学报（自然科学版）:1820230421.http:/ LUCC 的生态型山区土地敏感性评价系统设计J.现代电子技术,2023,46（1）:102.DOI:10.16652/j.issn.1004-373x.2023.01.020.22陈潇菡,张磊,史丽,等.基于地学信息图谱的中缅泰老“黄金四角”地区土地利用变化分析J/OL.云南农业大学学报（自然科学）:11020230421.http:/ 2035 年远景目标纲要的通知R.六安市人民政府公报,2021:1.28孙黎.城市规划批后管理及执法问题的思考:以安徽省六安市为例J.城乡建设,2020（21）:68.29胡晓峰.县级土地利用总体规划实施保障措施研究:以六安市裕安区为例J.现代农业科技,2019（13）:257.DOI:10.3969/j.issn.1007-5739.2019.13.139.162云南农业大学学报第17卷