首页 | 学术研究∶土地资源利用与保护 |
城市边缘区土地利用变化及人文驱动力机制研究* |
(南京大学 城市与资源学系,江苏 南京 210093)
区域土地利用变化已成为全球变化研究的热点之一[1]。1995 年 IGBP 和 HDP 联合提出了“土地利用与土地覆盖变化”研究计划,确定了两个主要研究方向:①通过实例研究揭示不同状态下覆被动态的特征与原因;②建立全球或区域土地利用/土地覆盖变化模型,通过模型揭示土地利用/土地覆盖与驱动力的相关关系,预测全球或区域土地利用/土地覆盖的未来趋势[2~5]。目前,国内专门分析土地利用变化与驱动力之间关系的文献尚少[6~10],但一些研究从侧面反映了社会、经济、技术及政策变化对土地利用变化的影响[11~15]。改革开放之后,中国政治、社会、经济均发生了翻天覆地的变化,这种变化又深深地烙印于土地利用格局之上。因此,运用 RS 和 GIS 技术分析土地利用动态变化的规律,不仅有助于了解引起土地利用变化的人文驱动力机制,为全球变化研究积累区域案例;亦为调整人类自身活动,优化配置资源,实现土地可持续利用服务。 无锡市马山区位于东经 120°03′12″~120°08′24″,北纬 31°21′48″~31°30′24″,地处太湖西北角,距市中心区 20km,为太湖中第二大半岛(图 1)。1987 年行政升格以后,国民经济快速、稳定、健康发展,城镇化水平大幅度提高,外商投资急剧增长,旅游事业蒸蒸日上,为分析城市化进程中土地利用变化规律及人文驱动力提供一个极佳的案例(表 1)。 图 1 研究区位置图 表 1 1988、1993、1998 年马山区社会经济变化情况 |
年份 | 人口数(人) | GDP(万元) | 三产比重 | 非农人口比重 (%) | 外商投资额(万美元) | 游客数(人次) |
---|---|---|---|---|---|---|
1988 | 9339 | 4500 | 36.0∶58.0∶6.0 | 8.60 | - | - |
1993 | 25368 | 34193 | 7.5∶64.9∶27.6 | 33.80 | 346 | 50000 |
1998 | 26872 | 106558 | 4.8∶54.9∶40.3 | 44.50 | 1124 | 1480000 |
1 土地利用变化遥感分析方法 1.1 影像解释与技术方法 国内学者已采用 RS 和 GIS 技术对土地利用/覆盖变化进行动态监测,取得一系列成果[16~21]。国外对土地利用变化的研究思路是 GIS 软件支持下,对不同时期遥感影像或土地利用图进行空间叠加运算,计算出各时期土地利用类型的转移矩阵[22],再进一步分析土地利用变化过程。本文采用国际上比较流行的土地利用变化的遥感分析、测量方法(图 2),对研究区 3 个时期的土地利用进行了定量的测量。确定 1993 年 1∶1 万地形图为基准,采用 3 次多项式及最近领域插值法对各期 LandsatTM 遥感影像进行了几何纠正,经检验配准误差小于 1 个象元,在纠正后的 TM 影像上提取了研究区有效面积为 48.21km,覆盖了马山区全域范围。再应用 SCR (Scattergram Controlled Regression)方法,以一次线性回归方程 Y=aX+b 的形式,对 3 个时期的 TM 影像进行了简易标准化处理[23]。 图 2 土地利用变化遥感监测技术路线 1.2 土地利用分类与转移矩阵计算 目前,因各种遥感影像分辨率不一,国内土地利用/土地覆盖遥感分类尚未有统一标准。本文采用 NASA 的分类体系[24],将研究区的土地利用类型划分为:建成区、开发区、农田、果园、灌草地、林地、水体、湿地和裸地,共 9 大类。依据图 2 的方法,通过与训练区得到的的结果参数相对比,计算影像各象元 ai 属于土地利用类型 xj 的概率 Pj(xj),由于可以根据最大 Pj(xj) 是否大于 0.7,将象元 ai 判定为确定象元或不确定象元。对于确定象元直接用最大似然判别法进行分类;对于不确定象元,进一步应用概率松驰法反复多次进行调整。 式中,k 为调整次数,i 为相邻的 8 个象元。直至│ Ak+1-Ak │≈0 为止(A 为不确定象元数目)。之后再用最大似然判别法进行分类。对照土地利用详查图,随机选取若干样本区,计算分类混淆矩阵和 Kappa 指数,对 3 期 TM 影像的分类结果进行检验,结果表明:Kappa 指数分别为 0.71(1988 年)、0.80(1993 年)、0.76(1998 年),均达到最低允许判别精度 0.70 的要求[25]。采用 ARC/INFO 对相邻两期土地利用分类结果进行空间叠加运算,获取土地利用变化的转移矩阵(表 2)。 表 2 马山区 1988~1998 年土地利用变化转换矩阵
注:行表示的是 k 时期的 i 种土地利用类型,列表示 k+1 时期的 j 种土地利用类型;A 表示的是 k 时期的土地利用类型转变为 k+1 时期各种土地利用类型的面积,即原始土地利用变化转移矩阵 Aij。 表示 k 时期 i 种土地利用类型转变为 k+1 时期 k 种土地类型的比例; 表示 k+1 时期的 j 种土地利用类型中由 k 时期的 i 种土地利用类型转化而来的比例。行、列的合计分别表示 k 时期和 k+1 时期各种土地利用类型的面积及其占总面积的比例。 表示 k+1 时期各种土地利用类型相对于 k 时期的变化程度。 2 土地利用变化的过程分析 2.1 土地利用变化过程的空间景观特征 根据表 2、表 3 和图 3,1988 年马山区土地利用的总体特色为农业景观,最主要类型是林地和农田,两者占全区总面积的 81.50%。1998 年该区比较典型的农业景观已被改造为城郊混合型景观,建设用地迅速扩大,成为最主要利用类型,占全区面积的 49.30%,农田和林地下降十分迅速,仅占全区面积的 41.10%。与此同时,农田、林地向果园、水体等经济效益高的利用类型转化,受人类影响低的类型向受人类影响高的类型转变。从各种景观类型的斑块数目、平均斑块面积、土地利用多样性指数等可以看出:随着经济发展,建设用地已成为最重要的类型,不仅斑块数目增加,斑块平均面积也明显增大,反映了建设用地的空间集聚不断加强;果园斑块数目和平均面积逐渐增加,总体上呈现了空间集中化特征;农田、林地和湿地的斑块表现为不断破碎化特征和分散化趋势,但湿地斑块数目却有所上升;水体、裸地和灌草地变化特征不明显。建成区和开发区的变化率最高,水体和灌草地的变化率是最小的。1993~1998 年土地利用转化速率总体上远低于 1988~1993 年的变化,反映了土地利用变化由急剧向平稳发展的过程。 表 3 个时期马山区土地利用空间景观特征参数
注:* 代表多样性指数 H, ** 代表优势度指数 D, P 是第 i 种土地利用类型占总面积的百分比,m 为研究区利用类型的总数;*** 代表破碎度指数 C, 为土地利用中所有类型斑块的总个数,A 为总面积。 图 3 不同时期遥感分类的土地利用图 2.2 建设用地的空间扩展过程 1988~1993 年,新增建设用地高达 1030.25hm2,建设用地增长了 224.47%;1993~1998 年,又增加建设用地 887.55hm2,建设用地增长了 59.60%。从图 3 中可以看出,新增建设用地除了集中分布于中心区 N-S 轴线方向的古竹公路沿线外,沿湖岸线的环半岛公路上新增了大量的旅游度假区,建设用地由中心区向外逐步推移。交通是建设用地扩展的重要依托,为了分析两者之间的关系,分别计算了交通干线的不同距离缓冲区(距交通干线一定距离的狭长地带,亦称为 buffer)内的土地利用类型的总转化率(P)及向建设用地的转化率(Q)。 P=Buffer 向土地利用类型发生变化的象元总数/Buffer 内的象元总数; Q=Buffer 向建设用地转化的象元总数/Buffer 内的象元总数。 对判读结果进行回归分析发现:P 与距交通干线距离呈线性相关关系;Q 则与距交通干线呈指数相关关系,即离交通干线越近,土地利用类型发生的转化越剧烈,向建设用地转化的可能性越大(图 4)。 图 4 土地利用转换率与交通干线的关系 3 土地利用变化的人文驱动力分析 3.1 土地利用变化的人文驱动力 选择无锡市马山区 1988~1998 年建设用地面积和农田为因变量(Y1、Y2),总人口(X1)、国民经济生产总值(X2)、固定资产投资额(X3)、外资利用额(X4)、第一、二、三产业比重(X5、X6、X7)、农副产品价格指数(X8),进行逐步回归分析,结果如下: Y1=0.688X1+0.315X3+0.457X4+0.531X7-358.46 (2) Y2=-0.902X1-0.675X7+0.333X8+3745.20 (3) 式(2)中,R2=0.975,通过 1% 的显著性检验。反映人口增长、第三产业的发展、外资的利用和固定资产投资增加是建设用地急速扩大的主要人文驱动力;式(3)中,R2=0.981,通过 1% 的显著性检验。反映人口增长、第三产业的发展、农副产品价格波动是农田变化的主要人文驱动力。马山区经济的快速增长得益于太湖国家级度假区的开发建设,中央和地方政府给予了一系列的优惠政策,固定资产投资逐年加大,外资不断地涌入,形成了“旅游为龙头”的新的产业结构。 3.2 未来土地利用变化趋势预测 随着人口的增长、旅游业发展,加大菜、果、肉等农副产品的需求。市场机制已经改变了农户的土地经营方式,追逐高利润成为农业生产的重要环节。种植业、蔬菜、果园和水产养殖的产值之比高达 1.0∶4.5∶6.5∶7.5[9],导致农业内部出现了结构调整,农田向果园和精养渔池转化,林地向果园和农田转化。这种转化不仅存在于农业内部,更多的是农业向非农建设转化,且这一趋势仍然保持着一定的增长势头。根据 1992 年土地利用详查及各年的变更调查,采用灰色系统动态模型 GM(1,1)进行模拟,取得农田(式 4)、园地(式 5)、林地(式 6)和建设用地(式 7)的预测方程如下: X1(i+1)=-131.015е-0.0538i+138.528 (4) X1(i+1)=-42.821е-0.0582i-40.162 (5) X1(i+1)=-10.186е-0.04308i+12.672 (6) X1(i+1)=-192.847е-0.0170i-189.633 (7) 可以看出,农田和林地仍然呈一定的减少趋势,但是农田减少的幅度已经相当有限;园地和建设用地呈增加趋势,但增加幅度已经远不如 1988~1998 年。 4 结论 (1)应用遥感影像,对土地利用进行分类和变化监测十分可行;以象元为基本单元,计算土地利用变化的转移矩阵,结果可以反映土地利用变化在空间上的分布格局。 (2)过去 10 多年中,马山区土地利用格局发生了根本性变化。从 1988 年比较典型的农业景观转变为 1998 年城郊混合景观类型,建设用地占用了大量的农田和林地,迅速成为区域中最为优势的土地利用形式,呈现出沿交通干线分布,且从中心区不断向周围地区扩张。农业用地内部出现了低收益向高收益转化,加剧了农田、林地的破碎化和分散化。 (3)马山区行政升格之后,人口迅速增长,外资的大量涌入,固定资产投资的逐年增加和旅游事业带动的第三产业发展是本区土地利用大规模、多样化转变的人文驱动力。同时,GM(1,1)模型模拟显示未来一段时间内土地利用变化总体趋势仍然延续着过去 10 多年的变化,但是幅度将大大减缓。 参考文献
第一作者简介 陈浮 (1974-),男,江苏射阳人,博士,从事土地利用与管理研究 ,E-mail:[email protected]。
|
| 关闭窗口 |