自然资源遥感, 2022, 34(1): 198-209 doi: 10.6046/zrzyyg.2021042

技术应用

基于1990—2019年Landsat影像的干旱区绿洲土地利用变化与模拟

宋奇,1, 冯春晖1, 马自强2, 王楠3, 纪文君4, 彭杰,1

1.塔里木大学农学院,阿拉尔 843300

2.北京大学地球与空间科学学院遥感与地理信息系统研究所,北京 100871

3.浙江大学环境与资源学院,杭州 310058

4.中国农业大学土地科学与技术学院,北京 100083

Simulation of land use change in oasis of arid areas based on Landsat images from 1990 to 2019

SONG Qi,1, FENG Chunhui1, MA Ziqiang2, WANG Nan3, JI Wenjun4, PENG Jie,1

1. College of Agriculture, Tarim University, Alar 843300, China

2. Institute of Remote Sensing and Geographical Information Systems, School of Earth and Space Sciences, Peking University, Beijing 100871, China

3. College of Environmental and Resource Science, Zhejiang University, Hangzhou 310058, China

4. College of Land Science and Technology, China Agricultural University, Beijing 100083, China

通讯作者: 彭杰(1977-),男,教授,主要研究方向为干旱区资源环境遥感监测机理与模型。Email:pjzky@163.com。

责任编辑: 李瑜

收稿日期: 2021-02-7   修回日期: 2021-05-20  

基金资助: 兵团中青年创新领军人才项目“棉田土壤剖面盐渍化的卫星遥感监测研究”编号(2020CB032)
国家重点研发计划项目“土壤综合观测与智能服务平台研发与应用”共同资助编号(2018YFE0107000)

Received: 2021-02-7   Revised: 2021-05-20  

作者简介 About authors

宋奇(1996-),男,硕士研究生,主要研究方向为国土资源与遥感。Email: tarimsongqi@163.com。

摘要

探明西北干旱区典型人工绿洲阿拉尔垦区的土地利用变化及未来发展情况,为类似地区土地利用变化的调控和管理提供参考。将逐年各月份影像进行多时相合成后,运用支持向量机分类方法获得逐年土地利用分类图,从面积变化、类型转化、空间动态变化3方面进行土地利用变化分析,采用CA(cellular automaton)-Markov模型模拟2050和2080年土地利用变化情况,基于累积距平法和通径分析方法探究其突变情况和驱动因素。研究结果表明: 1990—2019年,阿拉尔垦区耕地、园地、水体和建设用地的面积呈增加趋势,其中耕地和园地面积的增加主要是由塔里木河沿岸区域之外的未利用地转换而来。至2080年,垦区东北和东南部的未利用地将被逐渐开垦,耕地、园地和建设用地的面积将大幅增加。阿拉尔垦区土地利用类型面积在2005年发生转折性变化,耕地、园地和建设用地的面积急剧增加; 垦区土地利用变化的主要驱动因素为总人口、农业生产总值和棉花价格。研究结论: 在未来土地开发利用过程中,应当制定可持续发展的耕地开发政策,严格控制建设用地面积,构建合理的土地利用结构。

关键词： 土地利用变化; 土地空间动态模拟; 长时间序列; 驱动因子

Abstract

This study aims to explore the land use change and its future development trend in the Aral reclamation area, a typical artificial oasis in the arid region in northwest China and to provide a reference for the regulation and management of land use change in similar areas. After the multi-temporal synthesis of monthly images of each year, annual land use classification maps were obtained using the support vector machine method. Then, the land use change was analyzed from the aspects of area change, type transformation, and spatial dynamic change. Finally, the cellular automaton (CA)-Markov model was used to simulate the land use change in 2050 and 2080, and the sudden changes and their driving factors were explored using the cumulative departure method and the path analysis. The results of this study are as follows. During 1990—2019, the area of arable land, garden land, water bodies, and construction land in the Aral reclamation area showed an increasing trend. Among them, the arable land and garden land increased in area mainly due to the conversion of unused land outside the areas along the Tarim River. By 2080, the unused land in the northeastern and southeastern parts of the reclamation area will be gradually reclaimed. As a result, arable land, garden land, and construction land will significantly increase. The area of various types of land use in the Aral reclamation reached a turning point in 2005, showing a sharp increase in the area of arable land, garden land, and building land. This was mainly driven by total population, gross agricultural product, and cotton prices. It can be concluded that it is necessary to develop policies on the sustainable development of arable land, to strictly control the area of construction land, and to construct a reasonable land use structure in future land development and utilization.

Keywords： land use change; dynamic simulation of land space; long time series; driving factors

PDF (9825KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

宋奇, 冯春晖, 马自强, 王楠, 纪文君, 彭杰. 基于1990—2019年Landsat影像的干旱区绿洲土地利用变化与模拟[J]. 自然资源遥感, 2022, 34(1): 198-209 doi:10.6046/zrzyyg.2021042

SONG Qi, FENG Chunhui, MA Ziqiang, WANG Nan, JI Wenjun, PENG Jie. Simulation of land use change in oasis of arid areas based on Landsat images from 1990 to 2019[J]. Remote Sensing for Natural Resources, 2022, 34(1): 198-209 doi:10.6046/zrzyyg.2021042

0 引言

土地利用变化是影响生态环境的重要因素之一,已成为全球变化和可持续发展研究的重要内容[1,2]。土地利用变化改变了区域的物质循环、自然景观和各种生态过程,进而导致全球生态系统发生变化[3,4,5],使得土地利用变化与模拟成为研究区域和全球生态系统变化的主要趋势和潮流[6,7,8]。为降低区域土地利用变化对自然和社会环境的影响,有必要对连续长时序的区域土地利用变化进行分析并对未来土地利用变化情况进行模拟,以期为土地利用变化的调控和管理提供参考。

国内外研究者对土地利用变化的研究主要集中在土地利用时空变化特征[9]、土地利用变化机制或驱动力[10]、土地利用变化的环境影响[11]等方面,少有涉及连续长时间序列的土地内部结构详细转移过程与模拟。从研究区域上看,多数研究集中在东部沿海经济发达的城市区域[12,13],而着眼于西北干旱区经济较为落后的典型人工绿洲区域的研究甚少。同时,城市区域的土地利用变化研究所使用的数据源多为典型时间断面的单时相遥感影像和单一的社会统计数据[14],不仅未能有效捕捉连续性信息和解译更精细土地利用分类结果,特别是造成一些关键突变点或拐点信息的遗漏,而且也不能全面反映研究区地形地貌、路网和气象情况。因此,基于连续长时间序列的多时相遥感影像和多数据源(地形地貌、路网、气象和社会统计数据等)对人工绿洲区域的土地利用内部结构转变过程及定量模拟机制研究仍需进一步加强。

阿拉尔垦区位于我国西北干旱地区,属于典型的人工绿洲区域,垦区内荒漠区域的自然土壤表层常覆有一层盐结皮以防止水分流失,保护土壤内部结构,过度的土地开发利用破坏了表层盐结皮,加速了水分流失,影响了土壤水分和大气循环系统。同时,阿拉尔垦区紧邻塔克拉玛干沙漠,自然环境恶劣,常出现大风、扬尘和沙暴等极端天气,自然状态下的盐结皮能抵抗10级风的侵蚀,而表层盐结皮破坏后3级风就能引起扬尘,严重影响空气质量[15]。此外,植被也是研究区的天然生态屏障[16],但自然植被对生态系统的维护能力大于人工植被,而阿拉尔垦区内的植被以人工植被(农作物和经济林)为主,自然植被(林草地)为辅,且农作物多为1 a生植物,极易受到破坏。因此,有必要对垦区进行土地利用变化和模拟研究,为未来土地利用规划提供科学依据。

本研究以西北干旱区典型人工绿洲阿拉尔垦区为例,收集1990—2019年所有可用Landsat遥感影像、地形地貌、路网、气象和社会统计数据,将逐年各月份影像进行多时相合成,使用支持向量机(support vector machine,SVM)分类方法,辅以2019年高分二号遥感影像及野外勘察对多时相合成后影像进行解译分类,获得逐年土地利用分类图,从面积变化、类型转化、空间动态变化3方面进行土地利用变化分析,选用CA(cellular automaton)-Markov模型对阿拉尔垦区2010年土地利用情况进行模拟并与真实情况比较,验证模拟模型的可靠性后,模拟2050年和2080年的土地利用情景,探究长时序土地利用类型面积的突变情况和驱动因素,进一步说明了CA-Markov模拟模型的优势,为典型人工绿洲区域的生态环境综合管理和宏观的土地资源与社会经济统筹管理提供理论依据。

1 研究区概况

阿拉尔垦区位于新疆维吾尔自治区阿克苏地区,地处塔克拉玛干沙漠和天山山脉的交汇过渡地带,地理坐标为E80°30'~81°58',N40°22'~40°57'(图1),总面积4 105.92 km2。阿拉尔垦区地势相对平坦,气象多变,常有沙暴、扬尘和霜冻等灾害,全年降水少蒸发大; 生态系统属于荒漠—绿洲型,土壤表层常常附有一层盐结皮,以防土壤水分蒸发; 垦区种植结构单一,大面积种植棉花。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   研究区地理位置

Fig.1   Geographic position of study area


2 材料与方法

2.1 遥感数据源及其预处理

从美国USGS网站中(https://earthexplorer.usgs.gov/)获取1990—2019年间阿拉尔垦区(轨道号为146/32,空间分辨率为30 m)的所有Landsat遥感影像,并从中筛选出云覆盖度低于40 %的可用影像,共259景,各影像云覆盖度情况和每年影像数情况如图2所示。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   遥感数据

Fig.2   Remote sensing data


从图2可看出各影像云覆盖度情况多低于10 %,由于Landsat系列传感器不断更新,因此本文1990—1998年选用Landsat5,1999—2012年选用Landsat7,2013—2019年选用Landsat8影像,避免了多传感器合成所带来的误差,且每年影像数至少有3景。

为了进一步提高土地利用分类精确度,将所有影像逐一进行预处理,其中包括辐射定标、大气校正、配准、裁剪等操作,并提取各景影像的归一化植被指数(normalized difference vegetation index,NDVI),计算公式为:

NDVI=(NIR-R)/(NIR+R),

式中: NIR为Landsat近红外波段反射率; R为红光波段反射率。

将每年各月份影像进行NDVI最大值合成,得到逐年植被覆盖最全面的情景,从数据源方面提高分类精度,相关影像的合成流程如图3所示,对影像合成前后的局部放大对比中可以看出,NDVI最大值合成后的影像植被覆盖度更高,利于解译分类。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   多时相遥感影像合成

Fig.3   Multi-temporal remote sensing image synthesis


2.2 土地利用分类方法

参考《新疆生产建设兵团统计年鉴》(1990—2019年)中第一师各土地利用面积统计情况和已有的新疆土地利用分类体系[17,18],结合阿拉尔垦区土地利用类型特点,将研究区土地利用类型分为6类。各土地利用类型的地表特征、影像特征和实地照片情况如表1所示。

表1   阿拉尔垦区分类体系及解译标志

Tab.1  Classification system and interpretation signs of Alar reclamation area

类型地表特征影像特征影像照片
耕地棉地、小麦地、玉米地、高粱地等浅绿色,形状规则,一般为长条,均匀分布
林草地胡杨林、草场、柽柳等黄色或深黄色,形状不规则,一般分布在河流、道路两侧以及偏远地区
园地苹果地、枣地、葡萄地、香梨地等黄色,形状相对规则,多呈矩形,有明显边界线,连片分布
水体河流、沟渠、水库等白色或浅紫色,形状规则,有明显边界线
建设用地城镇、学校、广场、工业用地等灰色或深紫色,形状规则,周边常有规则的行道树,有明显边界线,连片分布
未利用地湿地、裸石、盐碱地、沙地等灰色或深紫色,形状不规则,多分布于偏远地区

新窗口打开| 下载CSV


以2019年为例,根据阿拉尔垦区分类体系及解译标志情况建立各土地利用类型的分类样本(图4),图上括号内为各类样本所占的像素个数, 共计1 024 951个分类像元数。根据各分类样本,辅以收集的相关统计数据、高分影像和实地调查,采用监督分类中的SVM分类方法进行土地利用分类。选用全连接条件随机场(fully connected conditional random field, FC-CRF)[19]这种分类后处理方法对分类后影像进行分类后处理。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   每类样本点个数和分布

Fig.4   Number and distribution of sample points for each category


由于本文采用了多源遥感影像(Landsat8,Landsat7和Landsat5)进行分类,此处以2019年的影像代表Landsat8,2010年的影像代表Landsat7,1990年的影像代表Landsat5,进一步分析不同数据源的分类结果(图5)和分类精度(表2)。从中可以看出不同数据源的总体分类精度和Kappa系数均大于80 %和0.8,并且各用地类型的制图精度和用户精度均大于75 %。整体而言,此分类方法得到的分类后结果优异,可用于后续相关研究分析。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   不同数据源的分类结果

Fig.5   Classification results of different data sources


表2   不同数据源的精度评价

Tab.2  Accuracy evaluation of different data sources

土地利用类型Landsat8
(2019年)		Landsat7
(2010年)		Landsat5
(1990年)
制图精
度/%		用户精
度/%		制图精
度/%		用户精
度/%		制图精
度/%		用户精
度/%
耕地93.1990.1491.3089.9190.5588.43
林草地87.4993.4785.8193.4185.1492.36
园地76.4694.5275.5593.5275.2292.18
水体84.1895.6582.0693.9881.3692.46
建设用地82.7682.1182.5881.0481.1780.17
未利用地98.2788.6698.1686.7497.0185.12
Kappa系数0.880.850.82
总体精度/%90.0589.8288.29

新窗口打开| 下载CSV


2.3 土地利用变化分析方法

从各土地利用类型的面积变化、类型转化和空间动态变化3方面对阿拉尔垦区土地利用变化展开分析。基于ArcGIS软件的空间分析功能,对各时间段起始年份和终止年份的分类后影像进行分析,并制作各土地利用类型变化的转移矩阵和空间动态分布图。

2.4 土地利用模拟数据的收集与处理

为提高土地利用模拟精度,参考相关文献[20,21,22]并结合阿拉尔垦区的土地利用现状,选取高程、坡度、坡向、距高速、铁路、公路和水系距离、气温、降水、人口、国内生产总值(gross domestic product,GDP)、全社会固定资产投资、第一产业、农业生产总值和棉花价格共15种影响因素进行土地利用模拟。

高程数据从中国科学院计算机网络信息中心国际科学数据镜像网站中获取,在ArcGIS软件中对高程数据进行重采样得到坡度和坡向信息,利用Near工具得到每个栅格到高速、铁路、公路和水系的距离; 气温和降水数据从中国气象科学数据共享网中获取,并进行反距离权重(inverse distance weight,IDW)空间插值得到气温和降水的空间分布情况; 人口、GDP、全社会固定资产投资、第一产业、农业生产总值和棉花价格从《新疆生产建设兵团统计年鉴》(1990—2019年)中获取,并利用多因子权重分配法将以阿拉尔垦区为基本统计单元的统计数据分布到栅格单元上。

2.5 基于CA-Markov模型的土地利用模拟

马尔科夫(Markov)模型侧重于对土地利用变化数量的预测,缺乏对各土地利用类型的空间性表达,而CA模型能表达空间系统的时空动态演化过程。将Markov模型和CA模型结合(CA-Markov模型)[23]进行土地利用类型空间动态变化的模拟和预测,其既能够对土地利用的数量变化进行表达,也能对空间动态变化进行演示。

利用IDRISI软件中的CA-Markov功能,首先利用1990—1999年10期土地利用数据和收集的15种影响因素数据,模拟阿拉尔垦区2010年的土地利用变化情况并与2010年真实土地利用情况进行比较。对模拟图像进行精度检验,以验证CA-Markov模型模拟的可靠性,使用1990—2019年30期土地利用数据和相关影响因素数据模拟出2050年和2080年的土地利用变化情况。

3 结果与分析

3.1 阿拉尔垦区土地利用变化分析

3.1.1 面积变化分析

对阿拉尔垦区1990—2019年间的最大值合成影像进行土地利用分类并统计各土地利用类型面积,得到垦区各年份土地利用变化情况及其面积变化情况(图6和图7)。分类后结果的分辨率为30 m×30 m,从图7(a)中可看出阿拉尔垦区耕地、园地、水体和建设用地面积呈增加趋势,经统计分别增加了1 147.19 km2,674.19 km2,107.58 km2和36.54 km2,林草地和未利用地面积呈减少趋势,经统计分别减少了33.67 km2和1 931.83 km2,各土地利用类型在1990—1995年间年均变化率最大(图7(b)),说明此期间垦区内土地利用变化最为剧烈。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   1990—2019年阿拉尔垦区土地利用变化情况

Fig.6   Land use change in Alar reclamation area from 1990 to 2019


图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   1990—2019年阿拉尔垦区各土地利用类型的面积变化情况

Fig.7   The area change of land use types in Alar reclamation area from 1990 to 2019


3.1.2 类型转化分析

1990—2019年阿拉尔垦区各土地利用类型的转化情况如表3所示。阿拉尔垦区耕地面积增加1 516.19 km2,减少368.99 km2,净变化量为1 147.2 km2,耕地面积的净增加主要是由未利用地转化而来; 林草地面积减少了39.63 km2,主要被转化成耕地; 园地由未利用地转化了1 160.44 km2; 水体面积增加主要由未利用地转化; 建设用地面积的增加主要是由耕地转化而来; 未利用地减少了1 934.48 km2,主要被开垦为园地和耕地。

表3   1990—2019年阿拉尔垦区各土地利用类型的变化转换矩阵

Tab.3  Transformation matrix of land use types in Alar reclamation area from 1990 to 2019(km2)

1990年
2019年耕地林草地园地水体建设用地未利用地合计
耕地—14.09832.594.940.00664.571 516.19
林草地0.38—0.360.580.004.645.97
园地349.648.81—15.690.001 160.441 534.58
水体2.866.7715.84—0.00103.36128.82
建设用地16.118.6011.590.04—1.4837.82
未利用地0.001.360.000.001.28—2.64
合计368.9939.63860.3821.251.281 934.48

新窗口打开| 下载CSV


不同时间段各土地利用类型变化特征不同(图8),1990—2000年,除未利用地外,各土地利用类型面积净变化量均增加; 2000—2010年,林草地、园地和未利用地面积净变化量均减少,分别为-14.37 km2,-272.63 km2和-536.23 km2; 2010—2019年,林草地、水体和未利用地面积净变化量为减少; 1990—2019年,耕地、园地和未利用地净变化量最明显,分别为1 147.19 km2,674.19 km2和-1 931.84 km2。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   阿拉尔垦区在不同时间段的土地利用变化特征

Fig.8   Characteristics of land use change in Alar reclamation area in different time periods


3.1.3 空间动态变化分析阿拉尔垦区不同时间段土地利用空间动态变化情况如图9所示。1990—2000年,未利用地发生了大面积开垦,主要分布在阿拉尔垦区中部; 2000—2010年,耕地、林草地和园地3种之间的相互转化主要分布在塔里木河沿岸区域; 2010—2019年,未利用地的开垦面积相对较少,主要分布在垦区东北和东南部; 1990—2019年,塔里木河沿岸区域主要发生的是耕地、林草地和园地之间的转化,沿岸区域之外的地方主要发生的是未利用地的开垦。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9   阿拉尔垦区在不同时间段的土地利用动态变化

Fig.9   Dynamic changes of land use in Alar reclamation area in different time periods


3.2 阿拉尔垦区2010年土地利用模拟和精度检验

将阿拉尔垦区2010年实际情况和模拟后结果进行对比(图10),可以看出大部分区域一致,只有一些地区存在细微差异,如中心区域的建设用地、东北和东南部的林草地。

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10   2010年阿拉尔垦区土地利用的实际及模拟情况对比

Fig.10   Comparison of real and simulated situation of land use in Alar reclamation area in 2010


通过建立混淆矩阵,进一步对2010年模拟后结果进行精度评价(表4)。表中对角线位置加黑为各类别正确分类的像素个数。通过表4可以看出,模拟精度的偏差主要是受到2005年城市规划政策新建工业园区的影响,建设用地面积增加,导致园地和建设用地精度较低。但整体看来,阿拉尔垦区的总体分类精度为89.82 %,Kappa系数为0.85,表明CA-Markov模型能够较好地模拟阿拉尔垦区土地利用变化情况,可对2050年和2080年研究区展开情景模拟。

表4   阿拉尔垦区2010年土地利用模拟的精度评价

Tab.4  The accuracy evaluation of land use simulation in Alar Reclamation Area in 2010

预测类别合计用户精度/%
真实类别耕地林草地园地水体建设用地未利用地
耕地54 61775092308672 66759 82491.30
林草地858182 0583 53610 3281 93913 453212 17285.81
园地5 2723 88065 6961 06234210 70186 95375.55
水体0594184 256040 170224 52582.06
建设用地000014 1462 98417 13082.58
未利用地08 1890208161457 741466 29998.16
合计60 747194 89670 249196 05317 4555277351067 135—
制图精度/%89.9193.4193.5293.9881.0486.74——

新窗口打开| 下载CSV


3.3 阿拉尔垦区2050年和2080年的土地利用模拟

对比研究区2050年和2080年的土地利用模拟情况(图11)和各土地利用类型的面积变化情况(表5),到2080年垦区未利用地仅占2.95%,耕地和园地大幅增加,共占垦区的75.41%,且空间分布较为均匀,到2080年建设用地面积迅速扩张到207.76 km2。这是由于阿拉尔垦区还受到总人口、农业生产总值、棉花价格、GDP、全社会固定资产投资、第一产业、年均气温和年均降水量等驱动因子的影响,导致模拟结果出现很大偏差。建议提前制定合理土地利用政策及生态修复措施,以提高各土地利用类型的资源利用率,促进土地资源均衡利用,开展生态环境治理,实现区域生态环境可持续发展。

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11   阿拉尔垦区2050年和2080年的土地利用模拟情况

Fig.11   Land use simulation situation of Alar reclamation area in 2050 and 2080


表5   阿拉尔垦区2050年和2080年各土地利用类型的面积变化情况

Tab.5  The area change of each land use types of Alar reclamation area in 2050 and 2080

年份耕地林草地园地水体建设用地未利用地
面积/km2百分
比/%		面积/km2百分
比/%		面积/km2百分
比/%		面积/km2百分
比/%		面积/km2百分
比/%		面积/km2百分
比/%
2050年1 679.3240.90388.429.461 285.5631.31345.318.41104.292.54303.027.38
2080年1 727.7742.08330.948.061 368.5033.33349.828.52207.765.06121.122.95

新窗口打开| 下载CSV


4 讨论

4.1 阿拉尔垦区土地利用类型面积变化的突变情况

探讨连续长时间序列土地利用类型面积变化的突变情况有利于区域土地信息监测、环境保护及制定可持续的土地利用政策[24]。采用累积距平法[25]分析连续长时间序列中的面积突变情况,从图12中可以看出耕地、园地、建设用地和未利用地在2005年出现明显转折,多数研究的突变情况则是由自然洪灾[26]和极端天气[27]影响所致,而本研究区的突变情况主要是由人为的政策干扰所致,这是由于阿拉尔垦区的政策导向所致,2004年推出西部大开发政策,新疆棉花产业扩大,农产品价格上涨,从2005年开始不断开荒扩地,耕地、园地和建设用地面积急剧增加,未利用地面积急剧减少。

图12

新窗口打开| 下载原图ZIP| 生成PPT

图12   阿拉尔垦区土地利用类型面积变化的突变情况分析

Fig.12   Analysis on the abrupt change of land use type area in Alar reclamation area


4.2 土地利用变化的驱动因素

对土地利用变化的驱动因素进行深入研究,有利于土地信息监测、合理利用并为土地资源可持续发展提供科学依据[28]。本研究是以西北典型的农垦区域为对象展开的土地利用变化分析,而多数研究则是以东部沿海经济较为发达的城市区域[29]为对象展开的土地利用变化分析,驱动土地利用变化的因素多为工业化[30]和城市化[31]带来的经济发展。本文利用高程、坡度等15种影响因素数据和土地利用数据,运用CA-Markov模型进行土地利用模拟,能够显示出各影响因素的驱动力大小,最终发现本文总人口、农业生产总值和棉花价格为主要驱动因素。这是由于随着人口增加,机械化作业水平提高,农业生产总值增加,棉花价格上涨,大量的土地资源被开发,导致研究区土地利用发生剧烈变化。

4.3 不同土地利用模拟模型的比较

本文以CA-Markov模型和多智能体ABM(agent-based model)模型[32]这2种最常用的土地利用模拟方法进行比较讨论,分析本文所用模型的优势。CA-Markov模型不仅能够通过简单的转化规则模拟出复杂的土地利用类型变化模型,为更高级的土地利用演化提供良好的过程分析环境,而且能与GIS和RS数据嵌合,提高GIS分析复杂自然现象和长时序动态建模的能力。ABM模型能够通过微观智能体及其与地理空间环境相互作用的方式进行模拟,从而反映土地利用变化过程中复杂的人地关系,主要呈现研究区的人为因素影响[33],缺乏对自然环境的考量。因此,CA-Markov模型是通过局部简单转换规则改变元胞状态,从宏观模拟土地利用变化结果,而ABM模型是通过智能体Agent交互作用,从微观模拟土地利用变化结果。这2种模型各具优势,但李少英等[34]分析近30 a相关土地利用与模拟的研究提出,CA-Markov模型结合了CA和Markov模型的优势,能够更精确地模拟土地利用变化情况。

5 结论

本文将每年各月份影像进行多时相合成后,运用支持向量机分类方法获得逐年土地利用分类图,从面积变化、类型转化、空间动态变化3方面展开土地利用变化分析,采用CA-Markov模型模拟2050年和2080年土地利用变化情况,探究土地利用类型面积变化的突变情况、驱动因素和不同模拟模型的优势,得出以下结论:

1)1990—2019年,阿拉尔垦区耕地、园地、水体和建设用地的面积呈增加趋势,林草地和未利用地的面积呈减少趋势; 耕地和园地面积的增加主要是由塔里木河沿岸区域之外的未利用地转换而来。

2)CA-Markov模型的总体分类精度可达89.82%,Kappa系数可达0.85,能够较好地模拟阿拉尔垦区土地利用变化情况; 模拟结果显示,至2080年垦区东北和东南部的未利用地将被逐渐开垦,耕地、园地和建设用地的面积将大幅增加。

3)土地利用类型面积在2005年发生转折性变化,耕地、园地和建设用地的面积急剧增加,未利用地面积急剧减少,这是由于阿拉尔垦区的政策导向所致; 垦区土地利用变化的主要驱动因素为总人口、农业生产总值和棉花价格。

(责任编辑: 李瑜)

参考文献

董祚继.

从土地利用规划到国土空间规划——科学理性规划的视角

[J]. 中国土地科学, 2020, 34(5):1-7.

[本文引用: 1]

Dong Q J.

From land-use planning to territorial space planning:A perspective of scientific and rational planning

[J]. China Land Science, 2020, 34(5):1-7.

[本文引用: 1]

胡苏李扬, 李辉, 顾延生, 等.

基于高分辨率遥感影像的神农架大九湖湿地土地利用类型变化及其驱动力分析——来自长时间尺度多源遥感信息的约束

[J]. 自然资源遥感, 2021, 33(1):221-230.doi: 10.6046/gtzyyg.2020100.

DOI:10.6046/gtzyyg.2020100      [本文引用: 1]

Hu S L Y, Li H, Gu Y S, et al.

An analysis of land use changes and driving forces of Dajiuhu wetland in Shennongjia based on high reso-lution remote sensing images:Constraints from the multi-source and long-term remote sensing information

[J]. Remote Sensing for Land and Resources, 2021, 33(1):221-230.doi: 10.6046/gtzyyg.2020100.

DOI:10.6046/gtzyyg.2020100      [本文引用: 1]

Yu H, Liu X M, Kong B, et al.

Landscape ecology development supported by geospatial technologies:A review

[J]. Ecological Informatics, 2019, 51:185-192.

DOI:10.1016/j.ecoinf.2019.03.006      URL     [本文引用: 1]

高文龙, 苏腾飞, 张圣微, 等.

矿区地物分类及土地利用/覆盖变化动态监测——以海流兔流域为例

[J]. 自然资源遥感, 2020, 32(3):232-239.doi: 10.6046/gtzyyg.2020.03.30.

DOI:10.6046/gtzyyg.2020.03.30      [本文引用: 1]

Gao W L, Su T F, Zhang S W, et al.

Classification of objects and LUCC dynamic monitoring in mining area:A case study of Hailiutu watershed

[J]. Remote Sensing for Land and Resources, 2020, 32(3):232-239.doi: 10.6046/gtzyyg.2020.03.30.

DOI:10.6046/gtzyyg.2020.03.30      [本文引用: 1]

Guo Y, Fang G, Xu Y P, et al.

Identifying how future climate and land use/cover changes impact streamflow in Xinanjiang Basin,East China

[J]. Science of the Total Environment, 2020, 710:136275.

DOI:10.1016/j.scitotenv.2019.136275      URL     [本文引用: 1]

Chen L D, Liu Y, Lyu Y H, et al.

Pattern analysis in landscape ecolo-gy:Progress,challenges and outlook

[J]. Acta Ecologica Sinica, 2008, 28(11):5521-5531.

DOI:10.1016/S1872-2032(09)60011-1      URL     [本文引用: 1]

杜金龙, 朱记伟, 解建仓, 等.

基于GIS的城市土地利用研究进

[J]. 自然资源遥感, 2018, 30(3):9-17.doi: 10.6046/gtzyyg.2018.03.02.

DOI:10.6046/gtzyyg.2018.03.02      [本文引用: 1]

Du J L, Zhu J W, Xie J C, et al.

An overview of urban land use research based on GIS

[J]. Remote Sensing for Land and Resources, 2018, 30(3):9-17.doi: 10.6046/gtzyyg.2018.03.02.

DOI:10.6046/gtzyyg.2018.03.02      [本文引用: 1]

张瑜, 赵晓丽, 左丽君, 等.

黄土高原土地利用变化对生态系统服务价值的影响

[J]. 自然资源遥感, 2019, 31(3):132-139.doi: 10.6046/gtzyyg.2019.03.17.

DOI:10.6046/gtzyyg.2019.03.17      [本文引用: 1]

Zhang Y, Zhao X L, Zuo L J, et al.

The impact of land use change on ecosystem services value in Loess Plateau

[J]. Remote Sensing for Land and Resources, 2019, 31(3):132-139.doi: 10.6046/gtzyyg.2019.03.17.

DOI:10.6046/gtzyyg.2019.03.17      [本文引用: 1]

Wang Y H, Dai E F, Yin L, et al.

Land use/land cover change and the effects on ecosystem services in the Hengduan Mountain region,China

[J]. Ecosystem Services, 2018, 34(a):55-67.

DOI:10.1016/j.ecoser.2018.09.008      URL     [本文引用: 1]

Wang J, Chen Y Q, Shao X M, et al.

Land-use changes and policy dimension driving forces in China:Present,trend and future

[J]. Land Use Policy, 2012, 29(4):737-749.

DOI:10.1016/j.landusepol.2011.11.010      URL     [本文引用: 1]

Wu J G.

Landscape ecology

[J]. Encyclopedia of Ecology(Second Edition), 2019, 4:527-531.

[本文引用: 1]

Antrop M.

The language of landscape ecologists and planners:A comparative content analysis of concepts used in landscape ecology

[J]. Landscape and Urban Planning, 2001, 55(3):163-173.

DOI:10.1016/S0169-2046(01)00151-7      URL     [本文引用: 1]

Liu Y L, Luo T, Liu Z Q, et al.

A comparative analysis of urban and rural construction land use change and driving forces:Implications for urbanerural coordination development in Wuhan

[J]. Central China.Habitat International, 2015, 47:113-125.

[本文引用: 1]

Kabisch N, Selsam P, Kirsten T, et al.

A multi-sensor and multi-temporal remote sensing approach to detect land cover change dynamics in heterogeneous urban landscapes

[J]. Ecological Indicators, 2019, 99:273-282.

DOI:10.1016/j.ecolind.2018.12.033      [本文引用: 1]

With global changes such as climate change and urbanization, land cover is prone to changing rapidly in cities around the globe. Urban management and planning is challenged with development pressure to house increasing numbers of people. Most up-to date continuous land use and land cover change data are needed to make informed decisions on where to develop new residential areas while ensuring sufficient open and green spaces for a sustainable urban development. Optical remote sensing data provide important information to detect changes in heterogeneous urban landscapes over long time periods in contrast to conventional approaches such as cadastral and construction data. However, data from individual sensors may fail to provide useful images in the required temporal density, which is particularly the case in mid-latitudes due to relatively abundant cloud coverage. Furthermore, the data of a single sensor may be unavailable for an extended period of time or to the public at no cost. In this paper, we present an integrated, standardized approach that aims at combining remote sensing data in a high resolution that are provided by different sensors, are publicly available for a long-term period of more than ten years (2005-2017) and provide a high temporal resolution if combined. This multi-sensor and multi-temporal approach detects urban land cover changes within the highly dynamic city of Leipzig, Germany as a case. Landsat, Sentinel and RapidEye data are combined in a robust and normalized procedure to offset the variation and disturbances of different sensor characteristics. To apply the approach for detecting land cover changes, the Normalized Difference Vegetation Index (NDVI) is calculated and transferred into a classified NDVI (Classified Vegetation Cover-CVC). Small scale vegetation development in heterogeneous complex areas of a European compact city are highlighted. Results of this procedure show successfully that the presented approach is applicable with divers sensors' combinations for a longer time period and thus, provides an option for urban planning to update their land use and land cover information timely and on a small scale when using publicly available no cost data.

Zhao L, Huang W, Chen J, et al.

Land use/cover changes in the Oriental migratory locust area of China:Implications for ecological control and monitoring of locust area

[J]. Agriculture,Ecosystems and Environment, 2020, 303:107110.

DOI:10.1016/j.agee.2020.107110      URL     [本文引用: 1]

Shen S G, Yue P, Fan C J.

Quantitative assessment of land use dynamic variation using remote sensing data and landscape pattern in the Yangtze River delta,China

[J]. Sustainable Computing:Informatics and Systems, 2019, 23:111-119.

DOI:10.1016/j.suscom.2019.07.006      URL     [本文引用: 1]

Wang X C, Dong X B, Liu H M, et al.

Linking land use change,ecosystem services and human well-being:A case study of the Manas River basin of Xinjiang,China

[J]. Ecosystem Services, 2017, 27(a):113-123.

DOI:10.1016/j.ecoser.2017.08.013      URL     [本文引用: 1]

Jia B Q, Zhang Z Q, Ci L G, et al.

Oasis land-use dynamics and its influence on the oasis environment in Xinjiang,China

[J]. Journal of Arid Environments, 2004, 56(1):11-26.

DOI:10.1016/S0140-1963(03)00002-8      URL     [本文引用: 1]

Tian G J, Liu J Y, Xie Y C, et al.

Analysis of spatio-temporal dynamic pattern and driving forces of urban land in China in 1990s using TM images and GIS

[J]. Cities, 2005, 22(6):400-410.

DOI:10.1016/j.cities.2005.05.009      URL     [本文引用: 1]

岳东霞, 杨超, 江宝骅, 等.

基于CA-Markov模型的石羊河流域生态承载力时空格局预测

[J]. 生态学报, 2019, 39(6):1993-2003.

[本文引用: 1]

Yue D X, Yang C, Jiang B H, et al.

Spatio-temporal pattern prediction of the biocapacity in the Shiyang River basin on the basis of the CA-Markov model

[J]. Acta Ecologica Sinica, 2019, 39(6):1993-2003.

[本文引用: 1]

郭瑞琦, 陆波, 陈恺霖.

基于CLUMondo模型的多情景土地利用变化动态模拟——以广西沿海城市为例

[J]. 自然资源遥感, 2020, 32(1):176-183.doi: 10.6046/gtzyyg.2020.01.24.

DOI:10.6046/gtzyyg.2020.01.24      [本文引用: 1]

Guo R Q, Lu B, Chen K L.

Dynamic simulation of multi-scenario land use change based on CLUMondo model:A case study of coastal cities in Guangxi

[J]. Remote Sensing for Land and Resources, 2020, 32(1):176-183.doi: 10.6046/gtzyyg.2020.01.24.

DOI:10.6046/gtzyyg.2020.01.24      [本文引用: 1]

季翔, 黄炎和, 林金石, 等.

基于CA-Markov模型与ANUDEM内插法的崩岗侵蚀量预估

[J]. 农业工程学报, 2018, 34(21):128-136.

[本文引用: 1]

Ji X, Huang Y H, Lin J S, et al.

Estimation of erosion amount in collapsed gully based on CA-Markov model and ANUDEM interpolation

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(21):128-136.

[本文引用: 1]

成超男, 胡杨, 冯尧, 等.

基于CA-Markov模型的城市生态分区构建研究——以晋中主城区为例

[J]. 生态学报, 2020, 40(4):1455-1462.

[本文引用: 1]

Chen C N, Hu Y, Feng Y, et al.

Construction of urban ecological zones based on CA-Markov model:A case study of the main urban of Jinzhong

[J]. Acta Ecologica Sinica, 2020, 40(4):1455-1462.

[本文引用: 1]

Matsa M, Mupepi O, Musasa T, et al.

A GIS and remote sensing aided assessment of land use/cover changes in resettlement areas:A case of ward 32 of Mazowe district,Zimbabwe

[J]. Journal of Environmental Management, 2020, 276:111312.

DOI:S0301-4797(20)31235-4      PMID:32927188      [本文引用: 1]

The research assessed land use/cover changes and their socio-economic impacts on local resident farmers in ward 32 of Mazowe District. A mixed method research design that combines both qualitative and quantitative approaches was adopted. Questionnaires, interviews and observations were employed to gather data from the field. GIS and Remote sensing techniques were used for land use/cover change detection and calculations. Data collected from the field was analyzed using Statistical Package for Social Sciences (SPSS) version 25.0 and Microsoft excel. Arc map 10.5 was used to analyze Remote sensing data. Four major land use/covers namely forest, grassland, water and cultivated were identified as baseline land covers of Mazowe ward 32 in 2000. The research revealed that there has been significant loss in forest, grassland and wetland/water cover and gain in cultivated and bare land cover between 2000 and 2018 in ward 32 of Mazowe District. Scarcity of high energy tree species, fruit trees, construction poles and extinction of wild animals are some of the observed impacts of forest cover loss. The study recommends that Mazowe EMA, RDC and Forestry Commission should cooperate in monitoring natural resources in Mazowe District to ensure sustainable utilization.Copyright © 2020 Elsevier Ltd. All rights reserved.

夏芳.

钱塘江流域气候变化及其对水文径流的影响

[D]. 杭州:浙江大学, 2016.

[本文引用: 1]

Xia F.

Climate change in the Qiantang River basin and its influence on local runoff

[D]. Hangzhou:Zhejiang University, 2016.

[本文引用: 1]

Zhuang Q W, Wu S X, Yan Y Y, et al.

Monitoring land surface thermal environments under the background of landscape patterns in arid regions:A case study in Aksu River basin

[J]. Science of the Total Environment, 2020, 710(25):136336.

DOI:10.1016/j.scitotenv.2019.136336      URL     [本文引用: 1]

陶帅, 邝婷婷, 彭文甫, 等.

2000—2015年长江上游NDVI时空变化及驱动力——以宜宾市为例

[J]. 生态学报, 2020, 40(14):5029-5043.

[本文引用: 1]

Tao S, Kuang T T, Peng W F, et al.

Analyzing the spatio-temporal variation and drivers of NDVI in upper reaches of the Yangtze River from 2000 to 2015:A case study of Yibin City

[J]. Acta Ecologica Sinica, 2020, 40(14):5029-5043.

[本文引用: 1]

Fan Q D, Ding S Y.

Landscape pattern changes at a county scale:A case study in Fengqiu,Henan Province,China from 1990 to 2013

[J]. Catena, 2016, 137:152-160.

DOI:10.1016/j.catena.2015.09.012      URL     [本文引用: 1]

袁静文, 武辰, 杜博, 等.

高分五号高光谱遥感影像的城市土地利用景观格局分析

[J]. 遥感学报, 2020, 24(4):465-478.

[本文引用: 1]

Yuan J W, Wu C, Du B, et al.

Analysis of landscape pattern on urban land use based on GF-5 hyperspectral data

[J]. Journal of Remote Sensing, 2020, 24(4):465-478.

[本文引用: 1]

朱永森, 曾永年, 张猛.

基于HJ卫星数据与面向对象分类的土地利用/覆盖信息提取

[J]. 农业工程学报, 2017, 33(14):258-265.

[本文引用: 1]

Zhu Y S, Zeng Y N, Zhang M.

Extract of land use/cover information based on HJ satellites data and object-oriented classification

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(14):258-265.

[本文引用: 1]

朱金峰, 周艺, 王世新, 等.

1975—2018年白洋淀湿地变化分析

[J]. 遥感学报, 2019, 23(5):971-986.

[本文引用: 1]

Zhu J F, Zhou Y, Wang S X, et al.

Analysis of changes of Baiyangdian wetland from 1975 to 2018 based on remote sensing

[J]. Journal of Remote Sensing, 2019, 23(5):971-986.

[本文引用: 1]

刘敬杰, 夏敏, 刘友兆, 等.

基于多智能体与CA结合模型分析的农村土地利用变化驱动机制

[J]. 农业工程学报, 2018, 34(6):242-252.

[本文引用: 1]

Liu J J, Xia M, Liu Y Z, et al.

Driving mechanism of rural land use change based on multi-agent system and cellular automata

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(6):242-252.

[本文引用: 1]

王兆林, 牙升业, 蒲海霞, 等.

基于改进CA-Markov模型的山地城市边缘区土地利用变化模拟

[J]. 农业工程学报, 2020, 36(16):239-248.

[本文引用: 1]

Wang Z L, Ya S Y, Pu H X, et al.

Simulation of spatiotemporal variation of land use in mountainous-urban fringes based on improved CA-Markov model

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2020, 36(16):239-248.

[本文引用: 1]

李少英, 刘小平, 黎夏, 等.

土地利用变化模拟模型及应用研究进展

[J]. 遥感学报, 2017, 21(3):329-340.

[本文引用: 1]

Li S Y, Liu X P, Li X, et al.

Simulation model of land use dynamics and application:Progress and prospects

[J]. Journal of Remote Sensing, 2017, 21(3):329-340.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-2
版权所有 © 2015 《自然资源遥感》编辑部
地址:北京学院路31号中国国土资源航空物探遥感中心 邮编：100083
电话:010-62060291/62060292 E-mail:zrzyyg@163.com
本系统由北京玛格泰克科技发展有限公司设计开发