近10年崇州市土地利用变化基本特征与空间格局*

潘荟交，王怀英，杨华容，彭文甫※

(1.电子科技大学成都学院，四川成都 611731； 2.四川师范大学地理与资源科学学院，成都 610101；3.西南土地资源评价与监测教育部重点实验室，四川成都 610068)

0 引言

土地系统是人类—环境交互作用复合地球系统中的重要组成部分，是人类—环境关系的纽带和桥梁，对理解和弄清人类—环境复杂关系具有重要作用[1]。土地利用/土地覆被变化(LUCC)是全球环境变化与可持续研究的重要内容，已成为国内外研究的热点问题之一[2]。土地利用变化是人类活动作用于陆地表层环境的一种重要方式和响应，建立模型对土地利用变化过程进行研究，是深入了解土地利用变化成因及过程、揭示土地利用对人类活动响应程度的重要手段和途径[3-4]。

目前已有学者利用各模型、方法从全球—国家—行政单元—流域宏观尺度上进行了土地利用变化对人类活动响应程度分析[5-7]。刘纪远等[8]采用卫星遥感信息源和相同的技术方法，对中国土地利用变化数据进行定期更新，提出并发展了土地利用动态区划的方法，基于国家尺度研究土地利用变化的空间格局与时空特征； 龙花楼[9]、马才学[10]等从省、市、区等行政单元出发讨论了土地利用变化及其驱动力。除了具有政策导向性、经济发达的“热点区域”之外，还有很多专家将研究对象放到了河海流域、喀斯特地貌山区、农牧过渡带等“脆弱区”[11]。其中山地平原过渡带，以生态脆弱、敏感著称，不仅具有复杂的地质环境，自然灾害频发，又是社会、经济、文化方面的过渡带，具有独特的地理景观[12]，是土地利用变化研究中具有代表性的选择。特别是随着土地资源研究不断深入、城镇化进程的加快及城乡统筹发展新要求的提出，作为建设过程中的重、难点区域，平原山地过渡带土地利用变化研究的重要性日益凸显。文章利用2002年、2014年两期遥感数据，选取成都崇州市作为研究区，从土地利用变化幅度、程度、方向、响应等基本特征与空间格局展开研究，为该区土地的合理利用提供数据支撑，同时为定性定量研究中小尺度下平原与山地过渡带的土地利用/土地覆盖的时空演变提供参考。

1 研究区概况

崇州市位于四川省岷江中上游成都平原西南边缘，西接邛崃山脉，北邻龙门山断裂带(图1)，地势从东南到西北逐渐升高，全境以街子镇、三郎镇、怀远镇为一线，线的西北为山地，线的东南为平原、丘陵，属于典型的平原—山地过渡带，过渡地区地质环境脆弱地震、滑坡等地质灾害频发。崇州市境内有西河、黑石河和金马河3条河流，大体形成“四山一水五分田”的土地结构[13]。崇州市属于汶川地震灾后恢复重建4个重灾县(市)的规划范围。目前现代农业初具规模，城乡建设步伐加快。

图1 研究区位置

2 数据来源及处理

2.1 数据来源及预处理

该文所涉及的数据包括遥感影像、DEM、崇州市矢量边界、野外实地采样数据、成都市和崇州市统计年鉴数据。遥感影像为2002年4月TM影像、2014年5月ETM+影像，分辨率均为30m×30m，影像质量良好，遥感数据来源于中国科学院遥感与数字地球研究所和地理空间数据云。崇州市1： 5万地形图，来源于西南土地资源评价与监测重点实验室； 土壤类型及数字高程模型(DEM)，数据来源于中国科学院资源环境科学数据中心； 野外实地采样数据为2016年7月，对崇州市进行为期1周的野外实地采样。

数据的预处理主要是大气校正、几何校正、裁剪等。基于ENVI 5.1，对遥感影像进行大气校正； 利用崇州市1： 5万地形图，对2002年TM影像进行几何校正，再应用校正后的2002年影像对2014年影像进行几何校正，误差RMS控制在0.5个像素以内。

2.2 土地利用分类及解译

土地利用分类系统是在2007年国家标准《土地利用现状分类》的规定划分基础上[14]，结合崇州市的实际情况及遥感影像可解译能力，将土地利用类型划分为建设用地、耕地、林地、水体、草地和未利用地。

影像解译采用人机交互解译方法对遥感影像进行分类，并利用Google Earth的高分辨率遥感影像，对解译过程中的不确定地物进行补判[15]，实现辅助分类； 通过野外实地采样数据，对2014年、2002年土地分类结果进行对比验证，分类结果精度分别达到85.7%、90.07%，Kappa系数均大于0.85，能够满足后期数据分析需求。

3 研究方法

3.1 土地利用动态度

单一土地利用动态度是刻画不同土地利用类型在一定时间段内的变化速度和幅度的指标[16]，其表达式为：

(1)

式(1)中，K为研究区某种土地利用类型年变化率；Ua、Ub分别为研究区初期及末期某种土地类型的数量；T为研究时段长，即年数。

3.2 土地利用程度

根据刘纪远等[8]提出的土地利用程度的综合分析方法，将土地利用程度按照土地自然综合体在社会因素影响下的自然平衡状态分为若干级，并赋予分级指数，从而给出了土地利用程度综合指数及土地利用程度变化模型的定量化表达式[17]。

(2)

式(2)中，Lj为研究区土地利用程度综合指数；Ai为研究区第i级土地利用程度分级指数；Ci为研究区域内第i级土地利用程度分级面积百分比；n为土地利用程度分级数。

3.3 土地利用空间转移

转移矩阵可全面又具体地刻画区域土地利用变化的结构特征与各个土地类型的变化方向[18]，该方法来源于系统分析中对系统状态与状态转移的定量描述，大多数土地利用变化研究均采用地类间转移矩阵分析，转移矩阵的数学模型为：

(3)

式(3)中，S代表面积;n代表土地类型的数目;i、j分别表示研究区期初和期末的土地类型。转移矩阵不仅能够反映研究期间的期初与期末的土地利用结构，还能反映时段内各类土地利用类型的转移情况，即未变、转入、转出的方向。

3.4 地形因子响应

地形因子是影响土地利用方式的重要环境因子[19]，高程、地形坡度分级是土地管理、地理信息服务的重要内容[20]，该文根据2004年《土地利用更新调查技术规定》，结合相关研究成果及研究区实际情况，对高程、坡度进行了分级[21]。对研究区土地利用变化地形因子响应作研究分析。

4 结果与分析

基于TM/ETM+影像提取了崇州市2002年、2014年土地利用数据(图2)。

4.1 土地利用变化的基本特征

(1)土地利用速度变化

近10年崇州市建设用地、未利用地变化速度较快，分别为5.72%和-1.04%； 林地的变化速度最小，为-0.23%； 耕地、水体、草地增减速度较为稳定。其中崇阳街道办、大划镇等12个场镇建设用地的年增长速度较快，均达到5%以上。

(2)土地利用程度变化

图2 土地利用类型分类

图3 2002～2014年土地利用变化空间分布

表1 变化图谱单元面积

序号转移方向面积(hm2)占比(%)1耕地→建设用地68103036802林地→耕地23401812643建设用地→耕地20920511304耕地→林地1333357205未利用地→林地1245426736林地→未利用地1067675777林地→建设用地968045238草地→林地482222619水→建设用地3621619610林地→草地32328175合计17024679198未转移90753038306发生转移18508501694研究区1092615310000

近10年崇州市土地利用程度由254.25增长到260.93，其中崇阳街道办、三江镇等8个位于海拔较低的平原地区乡镇土地利用程度较大，均在310以上，土地利用处于发展时期，人类活动对土地利用的影响的程度增加； 崇阳街道办、大划镇等5个场镇土地利用程度增长量最大，除了崇阳街道办作为行政中心位于地势较低的平原地区之外，其他乡镇均位于平原向山地的过渡带。其中崇阳街道办、大划镇土地利用程度大且增长速度快，人类活动对这两镇土地利用影响剧烈； 公议乡土地利用程度变化量低于0，处于调整期。崇州市土地利用开发地域差异明显，平原地区开发程度大但发展平稳，平原山地过渡带虽然目前土地利用程度不突出，但是变化剧烈，是今后平衡土地发展的重点关注区域，土地开发方向由平原向平原山地过渡带转移。

4.2 空间转移特征

合成土地利用变化图谱表明(图3)，生成的36类图谱单元中， 6类未发生转移，面积9.075 303万hm2，占研究区的83.06%； 30类图谱单元发生变化，转换面积合计1.850 850万hm2，占研究区总面积的16.94%，计算累计百分率(表1)，排在前10位的图谱单元变化面积1.702 467万hm2，占比91.98%。

在近10年间，增长面积最大的图谱单元耕地、林地、水体向新增建设用地的转化达到8 140.5hm2，其中6 810.30hm2的耕地转变，占全区转换类型面积的36.796%。主要分布在城市中心的崇阳街道办等周边乡镇以及西北方向的街子镇。2009年街子镇被批复为国家4A级旅游景，受到政策影响开发利用力度增强，较大面积的耕地转为建设用地； 面积排在第2位的图谱单元类型是林地、建设用地向耕地的转换，新增耕地面积4 432.23hm2主要集中在街子镇、怀远镇、道明镇山地向平坝过度的丘陵地区，应多注意地区林地保护； 位列第3位的是由耕地、未利用地、草地转换的新增林地，面积为3 060.99hm2，未利用地向林地的转换主要集中在鸡冠山乡北部高山地区，耕地向林地的转换主要集中在公议乡，主要原因是公议乡以丘陵地区生态环境综合治理项目为契机，以发展林果为重点，大力开展退耕还林工作； 其次，新增未利用地面积1 067.67hm2，集中分布在鸡冠山乡北部高山地带。林地与未利用地的相互转换受人类的影响小，属于自然演替状态。水体与草地转出、转入面积均较小，是最稳定的土地类型。

4.3 中心城市扩展方向

缓冲分析(图4)结果表明， 2002年中心城区缓冲范围内建设用地面积由2 808hm2增长到2014年的4 928.76hm2，增加了2 120.76hm2，耕地面积由7 396.29hm2缩减到5 401.8hm2，减少了1 994.49hm2，其他土地类型较为稳定。

通过研究建设用地的扩展指数，可以发现2002～2014年城中心主要是向东、东南、东北3个方向扩展(图5)，范围包括崇阳街道办、崇平镇、羊马镇、大划镇、集贤乡等区域； 在河流东侧土地利用程度大，得到高度开发，平均增长速度是23%、18%、8%，其建设用地增长贡献率达到了71%。在北、西北、西、西南、南5个区块，河流西侧的建设用地面积均缓慢增长，可能是今后城市土地规划、管理工作关注的方向。

图4 2002年(a)和2014年(b)中心城市缓冲区

图5 各区块建设用地平均扩展速率及扩展贡献率

图6 地形因子分级

4.4 土地利用的高程、坡度响应

研究区高程、坡度分级图(图6),表明高程、坡度各等级分布均相对集中，海拔高度走向自西北向东南递减趋势，结合提取的土地利用信息，可以得出：

建设用地、耕地、水体分布在0～1 400m，其中0～500m平原地区为集中分布区域。研究期间建设用地面积大幅增加，街子镇、怀远镇在500～1 400m高程内建设用地面积增加明显； 0～500m高程的耕地面积出现了大量减少，面积为3 889.62hm2，集中在崇阳街道办及其周边乡镇； 林地面积主要分布在海拔500m以上的鸡冠山乡等5个乡镇，其中街子镇、怀远镇林地面积减少明显； 草地面积在低海拔区域比例增加明显，反映了城市建设过程绿化意识的加强； 3 000m以上的高海拔地区林地与未利用地的相互转换明显。

坡度0°～6°是人类活动最频繁区域，新增的建设用地、耕地的大量减少、草地的增加绝大部分发生在此坡度地形中； 建设用地在各坡度等级面积均有所增加，在0°～2°坡度上面积增加量最大，为4 824.27hm2； 耕地在0°～2°上面积减少了3 623.33hm2； 林地主要分布在坡度>25°的区域， 2002～2014年该坡度范围面积变化小，分别是2.528 298万hm2、2.540 79万hm2，坡度1°～6°范围面积减少明显。

表2 土地利用与地形因子相关性

相关系数2002年2014年高程坡度高程坡度建设用地-0921-0959-0937-0945耕地-0903-09300908-0934林地0716095207560952水体-0906-0967-0889-0957草地0280-09580162-0983未利用地0544-09370636-0928

由表2可见，建设用地、耕地、水体与高程、坡度之间负相关性显著， 2002年、2014年相关系数均达-0.9水平，表明地形因子对建设用地、耕地、水体的空间分布起制约作用。草地与坡度存在较高的负相关性，两期相关系数为分别为-0.958、-0.983，而与高程的相关系数较小； 未利用地与坡度存在较高的负相关性，两期相关系数为分别为-0.937、-0.928，而与高程的相关系数一般，这一结果表明坡度对草地、未利用地分布的影响程度较大且高于高程。

5 结论

(1)土地利用类型以林地、耕地、建设用地为主，建设用地变化速度最快，动态度达到5.72，崇州市整体土地利用程度由2002年的254.25增长到2014年的260.93，增加了6.68，土地直接经济价值的提升程度大，土地利用处于发展时期。以能集中反映人类活动区域差异性的乡镇为研究单元，地域差异明显，位于平原山地过渡区域的乡镇土地利用开发程度显著增强。

(2)近10年来空间发生转换的土地类型面积合计1.850 850万hm2，占研究区总面积的16.939%，建设用地、耕地、林地在空间转移中是最活跃的土地类型。中心城区建设用地面积大幅增长，2002～2014年中心城区缓冲范围内建设用地面积增加了接近1倍，建设用地空间扩展规律呈现自河流向东、东南、东北3个方向扩展，这3个方向需要更加重视土地集约化管理。

(3)地形条件对土地资源的空间布局起着制约作用，建设用地、耕地、水体与高程、坡度之间负相关性显著。海拔低、坡度较缓的平原地区依然是人类活动最频繁的区域，建设用地增加、耕地大量减少绝大部分发生在该地区； 山地平原过渡带是今后土地利用发展的重点地区。该区域耕地、建设用地的增加主要来自林地的减少； 高海拔地区林地与未利用地的相互转换属于自然演替状态。

掌握区域土地利用结构变化规律，深化对城镇发展建设历史的科学认识，对该区域土地利用调控和管理具有一定的现实意义，对城镇建设意见、经济社会政策、城市建设布局规划的制定具有重要参考意义，为土地利用结构更加合理提供依据。

[1] 吴琳娜，杨胜天，刘晓燕，等.1976年以来北洛河流域土地利用变化对人类活动程度的响应.地理学报， 2014, 69(1): 54～63

[2] 张新荣， 刘林萍，方石，等.土地利用、覆被变化(LUCC)与环境变化关系研究进展.生态环境, 2014，(12)： 2013～2021

[3] 徐勇， 孙晓一，汤青.陆地表层人类活动强度：概念、方法及应用.地理学报， 2015, 70(7): 1068～1079

[4] Morse N B,Wollheim W M.Climate variability masks the impacts of land use change on nutrient export in a suburbanizing watershed.Biogeochemistry， 2014, 121(1): 45～59

[5] 曹小敏， 李爱农，雷光斌，等.尼泊尔土地覆被遥感制图及其空间格局分析.地球信息科学学报， 2016, 18(10): 1384～1398

[6] 马才学， 孟芬，赵利利.1990～2005年武汉市土地利用时空变化及其政策驱动因素分析.水土保持研究， 2015，(2)： 117～122

[7] 樊敏， 李富程，郭亚琳，等.退耕还林对岷江上游高山聚落区生态服务价值变化的影响.山地学报， 2016, 34(3): 356～365

[8] 刘纪远， 张增祥，徐新良，等.21世纪初中国土地利用变化的空间格局与驱动力分析.地理学报， 2009, 64(12): 1411～1420

[9] 龙花楼， 王文杰，翟刚，等.安徽省土地利用变化及其驱动力分析.长江流域资源与环境， 2002, 11(6): 526～530

[10]马才学， 孟芬，赵利利.1990～2005年武汉市土地利用时空变化及其政策驱动因素分析.水土保持研究， 2015, 22(2): 117～122

[11]徐嘉兴， 王继尧，梁涛，等.近18年太湖流域土地利用变化及区域差异分析.地理空间信息， 2009, 7(4): 48～51

[12]钟兆站， 李克煌.山地平原交界带与区域可持续发展.地理科学进展， 1998, 17(2): 23～31

[13]成都市土地利用总体规划(2006-2020).国土资源部信息中心， 2012

[14]焦琨， 杨子生.我国土地利用分类新旧标准的对比分析.现代农业科技， 2008,(23)： 296～297

[15]郑泽忠， 范东明，李玉霞，等.多种对地观测技术及Google Earth技术应用于土地利用/覆被变化建库.地球物理学进展， 2008, 23(4): 1326～1331

[16]葛静， 王丹，杨萍果.临汾市土地利用变化的格局分析及其对生态系统服务价值的影响.中国农业资源与区划， 2016, 37(10): 77～82

[17]杨华容， 王怀英，彭文甫，等.区域土地利用/覆被时空动态变化研究——以四川省金堂县为例.中国农业资源与区划， 2016, 37(8): 37～46

[18]张琪， 罗格平，李龙辉，等.基于土地利用/覆被变化表征的现代绿洲演变过程——以天山北坡三工河流域为例.地理学报， 2016, 71(7): 1157～1171

[19]尚正永， 张小林，周晓钟.基于RS/GIS的城市空间扩展与外部形态演变研究——以江苏省淮安市为例.经济地理， 2012, 32(8): 64～70

[20]刘艳艳， 吴大放，董玉祥，等.珠海市土地利用空间格局与地形的相关性研究.地理与地理信息科学， 2010, 26(3): 68～72

[21]马群， 赵庚星，李玉环.小尺度地形因子对农地土壤质量的影响研究.地理与地理信息科学， 2009, 25(3): 77～80