邵姝遥，曹罗丹，田鹏，赵梦琪，张柯琴，邹之晗，周静芬，周艺

(宁波大学 地理与空间信息技术系，浙江 宁波 315211)

作为国际地圈生物圈计划(IGBP)和全球变化人文计划(IHDP)两大国际合作项目共同支撑的核心研究计划，土地利用变化一直是地球科学研究的热点[1]。土地利用变化既能反映自然格局、过程与人类社会之间最直接的相互关系，更是人口、资源、环境可持续发展的基本问题[2]。目前，国内外土地利用变化研究集中在时空格局分布特征和时空开发利用强度[3]、土地利用变化驱动因素探析，以及土地利用变化对区域景观格局、生态系统稳定性的影响，包括生态风险评价、生态服务价值对土地利用的响应分析[4-6]等方面。如冯佰香等[7]基于象山港土地利用变化，分析其土地开发利用强度；刘永超等[8]对比了中国象山港和美国坦帕湾土地利用变化对生态服务价值的影响。

甬江流域坐落于我国经济发达的东南沿海区域，广布的河流水系网和充沛的河流水量孕育了河姆渡七千年文明和“宁波帮”传承精神，推快了“宁波城”走向繁荣和兴盛的步伐。目前，国内对甬江流域的研究集中于河流水污染方面[9]。随着经济飞速发展，甬江流域的土地利用正发生着急剧变化，土地利用变化破坏了流域景观格局的稳定性，对区域生物多样性和生态系统服务功能亦会产生一定的影响。在实行最严格的生态环境保护制度的背景下，了解人类活动对甬江流域土地利用类型的影响，以及从较长时间段上分析甬江流域土地利用的时空变化特征，有利于流域合理利用土地资源，促进生态环境的治理和保护。为此，特基于甬江流域土地利用变化，分析其土地利用转换方向和结构，以及土地利用开发强度，为区域合理开发和保护土地资源、维护生态环境稳定提供理论和实践指导[10]。

1 研究区概况

甬江流域位于我国东南沿海，浙江省东北部，宁绍平原东部的三江平原地区，地理坐标29°24′～30°49′N、120°49′～120°56′E。甬江流域面积约为3 788 km2，全长105 km，是宁波最主要的河流之一，由奉化江和余姚江两江在宁波市区三江口汇集而成。甬江流域处于亚热带季风气候区，地势西南高、东北低。甬江流域行政范围大部分位于宁波市的江北区、鄞州区、海曙区、镇海区、奉化区及余姚市，小部分位于北仑区、慈溪市及宁海县，区位条件优越，交通便利。甬江流域在促进宁波社会经济快速发展中起到了重要的推动作用，但在快速工业化和城市化的过程中，土地利用类型与结构变化明显。

2 数据来源与研究方法

2.1 数据来源

本研究使用的数据来自于地理国情监测云平台(http://www.dsac.cn/)，包括浙江省的1980、1990、1995、2000、2005、2010、2015共7期1∶100 000土地利用数据。结合甬江流域行政矢量边界，叠加得到研究所需数据。根据全国土地资源分类系统、甬江流域自然生态背景及本文实际需要，把土地利用类型划分为8类：耕地、林地、草地、水域、建设用地、未利用地、滩涂，以及包括在流域边界线内的海域[11]，滩涂主要是沿海大潮高潮位与低潮位中间的潮侵地带，其他土地利用类型内涵详见文献[12]。

2.2 研究方法

2.2.1 土地利用转移矩阵

利用马尔可夫模型计算出不同时期土地开发利用状态的转移矩阵，对甬江流域每个时期的土地利用数据进行叠加，分析得到甬江流域35 a间土地利用类型的转移方向、转移面积等变化情况，具体计算公式参见文献[13]。

2.2.2 土地利用动态度

土地综合利用动态度与单一利用动态度指数可以用来描述区域和某种土地利用类型变化的强度和速度[14]。土地单一利用动态度表达式：

K=(Ub-Ua)÷(UaT)。

式中：K为研究时段内某一土地利用类型动态度，Ua、Ub分别为研究期初及研究期末某一土地利用类型的数量；T为研究时段长。土地综合利用动态度的计算公式参见文献[14]。

2.2.3 土地利用程度

参照前人的研究[15]，将土地利用程度数量化，采用土地利用程度指数来研究土地的开发利用程度，其计算公式为

式中：L为研究区的土地利用程度综合指数；Ai为研究区内第i级土地利用程度分级指数；Ci为研究区内第i级土地利用程度的面积百分率；n为土地利用程度分级数。据前人经验[15]及研究需要，划分土地利用程度分级：未利用土地级，包括的土地利用类型为未利用地、滩涂，分级指数1；林、草、水用地级，包括的土地利用类型为林地、草地水域，分级指数2；农业用地级，包括的土地利用类型为耕地、养殖用地，分级指数3；城镇聚落用地级，包括的土地利用类型为建设用地，分级指数4。

3 结果与分析

3.1 土地利用总体变化

对甬江流域1980—2015年的土地利用类型面积进行整理发现(图1)，甬江流域各土地利用类型在数量和结构上发生了较大的变化：草地面积及比例在前10 a大幅减少，此后25 a变化面积较小；耕地面积及比例在1995—2000年前有微小增加，但总体上呈持续减少态势，减少了553.51 km2；海域面积及比例在2010年及以前基本无变化，但在此后5 a内比例快速减小至消失；建设用地持续快速增加且增速明显，从1980年的3.91%增加到2015年的18.60%，共增加了556.63 km2；林地面积及比例以1995年为峰值先增加后减小，但研究初期至后期变化总体较小；水域面积持续增加，但增速不均；滩涂面积比例在1980—1990年增加，在1990—2000年保持稳定，在2000—2005年迅速减小，在2005—2010年保持稳定，于2010—2015年减小至消失；未利用地面积及比例于1980—1990年减小明显，随后5 a有较小增加，之后保持平稳。

图1 甬江流域各时期土地利用类型结构

3.2 土地利用转移分析

借助转移矩阵来分析土地利用类型在不同时期的转移方向和数量(表1)。1980—2015年，研究区内转移面积最大的为耕地转向建设用地，89.36%的耕地转移成建设用地，转移面积为521.90 km2。草地93.03%转出为林地，海域主要转出为建设用地及草地。建设用地转入来源多，转出较少，主要由耕地和林地转移而来，转入比重分别为92.61%、6.29%。林地主要转出为建设用地、耕地和草地，水域主要转出为耕地和建设用地，滩涂主要转出为水域，未利用地主要转出为建设用地。海域和滩涂在研究初期面积较小且在研究期末全部转为其他土地利用类型。甬江流域在35 a土地利用转型的过程中，各土地利用类型之间转换复杂。

表1 1980—2015年甬江流域土地利用转移矩阵

在土地利用转移空间上(图2)，甬江流域耕地转为建设用地分布较广泛，主要位于地形平坦、地理位置较好的地区，靠近人类活动强烈的城镇区域，这些区域的土地利用强度高、土地利用需求量大。在退耕还林政策背景下，耕地在部分地区转换为林地，主要发生在山地、丘陵等地形起伏相对较大且不适宜耕种的地区。耕地转为水域主要发生在东部地区，在沿海水产养殖兴起的背景下土地利用方式发生转变。在人口增加和保障粮食安全的背景下，林地转为耕地主要分布在地形相对平坦、适合耕种的地区。林地向建设用地的转换聚集在城镇的郊区，处于城市扩张的范围内，受经济发展水平和速度影响较大。草地向林地转换集中在低缓山地和丘陵地形区，处于甬江流域中部偏西地区。

3.3 土地利用速率分析

计算甬江流域35 a间各个时期各土地利用类型的综合及单一利用动态度(表2)，分析其土地利用变化速率。甬江流域35 a土地综合利用动态度为0.42%。其中，1990—1995、2000—2005年2个时间段内土地综合利用动态度较大，1980—1990、2005—2010年2个时间段土地利用综合动态度较小，反映了甬江流域土地利用变化在前期和后期较为稳定，中期比较明显。

图2 1980—2015年甬江流域土地利用转移的空间分布

表2 1980—2015年甬江流域土地利用变化速率

建设用地单一利用动态度35 a间均为正数，面积持续增长。耕地面积虽减少，但由于面积基数较大，所以动态度较小，研究期末下降了0.51%。草地、未利用地单一利用动态度起伏变化，总体上呈上升趋势。水域土地利用变化速率呈减小趋势，海域、滩涂利用动态度在2010—2015年间变化较大。林地利用动态度变化最小，虽在各个小时间段内均有变化且有个别年份变化幅度较大，但总体较为平衡。

土地综合利用动态度和单一利用动态度在不同时间段内表现出较大的差异，反映了甬江流域开发利用的不均衡性和不确定性。开垦农田与退耕还林同时进行，增大了经济成本，破坏了生态环境，耕地与建设用地的交替转换，在一定程度上降低了土地的生产力，破坏了土地的可持续利用性。

3.4 土地利用程度分析

通过计算得到甬江流域土地利用程度指数，运用ArcGIS 10.3软件对其中的1980、1995、2005、2015年土地利用程度进行空间化，再对评价结果进行分等定级，将评价结果分为弱、较弱、中、较强和强5个等级，以反映甬江流域土地利用程度的空间变化。时间上，1980—2015年间，甬江流域土地利用程度总体上呈上升趋势，反映了研究期间人类对甬江流域的土地利用程度不断上升。其中，1990—1995年，土地利用程度有所下降，而2000—2005年土地利用程度综合指数上升最快，增加了8.27。

空间上，35 a间甬江流域土地开发利用程度区域差异明显(图3)，北部高于南部，东部高于西部。整个甬江流域范围内，1980年土地利用中等以上程度合计占65.08%，较强程度的土地利用面积占41.24%，强等级土地利用程度分布在宁波市区范围的中心地带，面积占比为1.37%，弱等级土地开发利用程度的土地分布在研究区东北角。1995年土地开发利用程度为中等以上的土地面积减少，总计占甬江流域总面积的61.88%，研究区北部即余姚市东北部和宁波市区范围西北部较大面积土地的利用程度由较强级转变为较弱级。2005年土地开发利用程度为中等以上的土地面积增加，占甬江流域总面积的66.21%，原有的开发利用程度为强的余姚市中心地带、宁波市区范围的中心地带及奉化中心地带以其自身为中心得到了快速的扩张。2015年土地开发利用程度为中等以上的土地面积比例增加至最大值67.35%，土地利用程度达到最大。

图3 甬江流域土地利用程度空间分布

总体来说，研究区内土地开发利用程度普遍较高，中等以上程度面积呈增加趋势，其中开发利用程度为强的土地面积所占比例稳定持续增加且增加速度越来越快，较强开发利用程度土地面积整体上呈持续减少状态，中等程度的土地面积变化不明显，稍有减少，较弱级及弱级土地持续减少，在2015年达到最小值。由于西南部山地丘陵对土地开发利用的限制，东部及北部的土地开发利用程度越来越高。耕地中化肥农药的使用，居民点及工厂企业生活垃圾、固体废弃物的污染，使得甬江流域东部及北部的土地受到人类活动的严重影响，来自生产生活的压力越来越大。

4 小结

对甬江流域1980—2015年的土地利用时空格局变化特征进行分析，结果发现：甬江流域草地、林地、耕地、海域、滩涂、未利用地面积缩减，建设用地和水域面积持续增长。耕地主要转为建设用地、林地和水域，海域、林地和未利用地主要转为建设用地，草地主要向林地转换，土地利用转移趋向复杂化。甬江流域耕地转为建设用地主要位于地形平坦、地理位置较好的地区，靠近人类活动强烈的城镇区域。耕地转换为林地主要发生在山地、丘陵等地形起伏相对较大且不适宜耕种的地区。耕地转为水域主要发生在东部地区。林地转为耕地主要分布在地形相对平坦、适合耕种的地区。林地向建设用地转换聚集在城镇的郊区，处于城市扩张的范围内。草地向林地转换集中在低缓山地和丘陵地形区，处于甬江流域中部偏西地区。甬江流域土地综合利用动态度为0.42%。建设用地单一利用动态度35 a间均为正数，耕地单一利用动态度下降了0.51%。草地、未利用地单一利用动态度总体呈上升趋势，林地单一利用动态度变化最小。土地综合利用动态度和单一利用动态度在不同时间段内表现出较大的差异，反映了甬江流域开发利用的不均衡和不确定性。甬江流域土地利用程度不断上升，研究期末增长了57.89%，其区域分异明显，受地形条件、地理区位、经济发展等因素影响，土地利用程度空间上呈现北部高于南部、东部高于西部的特征。