万 方, 邓清海, 刘 莉, 李 萍, 杨 鑫

(山东科技大学 地球科学与工程学院, 山东 青岛 266590)

青藏高原被称为地球的第三极，其土地利用与土地覆被(land use and land cover, LULC)变化不仅影响高原本身的生态环境状况和区域发展，同时也会对区域和全球产生巨大的环境效应[1]。因此，青藏高原土地利用/覆被变化一直都是全球变化研究的热点之一[2-3]。耕地是河湟谷地土地利用的主要形式之一，是保障粮食安全、生态安全以及区域发展的关键。耕地变化受人为活动影响最大，其特征和动态是土地利用/覆被变化的核心内容之一，也已经成为全球变化研究的重要组成部分[4]。青藏高原受高寒和复杂地理环境的限制，耕地资源相对稀少[5]，不到整个高原土地总面积的1%[6]，且主要分布在河湟谷地、“一江两河”以及川西、藏东和滇西北地区[7]，但其对当地社会经济的发展却具有举足轻重的作用。

河湟谷地位于青藏高原东北缘青海省的东北部，是青藏高原面积最大的农区，也是青海省最重要的农区，总面积为3.30×104km2，仅占青海省面积的5.18%，却集中了全省72%的人口和73.15%的耕地[8]。可以说河湟谷地是青海省乃至整个青藏高原地区粮食安全的有力保障[9]。而且作为青藏高原地区人口最为密集、经济最发达的地区和最重要的经济增长带，近年来，随着经济和城镇化的快速发展，河湟谷地土地利用发生了显著的变化。以往对于河湟谷地耕地的研究，都是一些大尺度，基于低精度的影像的研究，并且对此地区近年来的耕地变化规律也尚不清楚。因此，本文基于最新的较高分辨率Landsat-5 TM，Landsat-7 ETM+，Landsat-8 OLI系列卫星遥感影像和数字高程模型等相关数据，并结合空间统计学等手段，来分析河湟谷地地区近30 a耕地的分布特征和时空变化规律，以期为研究区土地利用、农业发展规划、耕地可持续利用和生态环境保护等提供依据，也为全球耕地变化研究提供典型区域基础。

1 研究区概况

河湟谷地是黄河与湟水流域肥沃的三角地带，位于青藏高原东北部，青海湖以东，经纬度范围为东经100°—103°，北纬35°—38°,是青海省东部的重要农业区，其行政区划由西往东依次为湟水流域的海晏县、湟源县、湟中县、大通县、西宁市、互助县、平安县、乐都县、民和县和黄河流域的贵德县、尖扎县、化隆县、循化县，其面积约为3.30×104km2。

河湟谷地位于达坂山和拉脊山之间，是湟水河和黄河冲击形成的河谷地区，地貌类型主要以高山或极高山、洪积—冲积平原、丘陵和台地为主。相较于青藏高原其他地区，其海拔高度相对较低，平均海拔在2 000～3 000 m之间，年平均温度在5～9 ℃之间，年平均降水量在252～535 mm之间，气候温和，雨热同期，对农作物生长发育极为有利，是黄河流域人类活动最早的地区之一。河湟谷地种植的主要作物有小麦、薯类、玉米、蚕豆、青稞等作物，作物熟制为一年一熟[10]。

2 数据来源与分析方法

2.1 数据来源及处理

遥感数据来自Landsat系列陆地卫星影像，共16景，其中1990年和2010年的数据来自Landsat-5 TM，2000年的数据来自Landsat-7 ETM+，2018年的数据为Landsat-8 OLI，详细参数见表1。为保证4个时间的数据空间分辨率一致，以便于比较，我们选取多光谱波段，空间分辨率均为30 m。DEM数据为GDEMDEM 30 m分辨率数字高程模型。

表1 数据源参数

对上述Landsat遥感数据进行辐射定标、大气校正、镶嵌裁剪等预处理后，采用监督分类中的最大似然分类法对河湟谷地的土地类型进行分类[11]，分类后提取耕地。研究区耕地在遥感影像中的几何特征较为规则，地块大而且排列整齐很容易被识别，因此可以人工选择训练样本进而进行监督分类。

最后，结合Google Earth等高分辨率影像和当地的耕地实际调查信息选取了200个随机位置的控制点，对土地利用类型图进行了精度评估。

2.2 分析方法

本次研究主要分析河湟谷地耕地的时空变化以及耕地和地形的关系。其中，耕地的时间变化主要用耕地扩张动态度[12]和耕地面积相对变化率[13]来衡量；耕地的外部空间变化主要是把耕地当作一个整体来分析河湟谷地的耕地相对于整个区域的变化趋势，主要采用的指标有耕地扩张动态度、耕地面积相对变化率[14]、耕地质心转移指数等[15]；最后，结合地形数据，分析研究区耕地随着时间变化在不同海拔和坡度的分布变化。

3 结果与分析

3.1 精度检验结果

结合Google Earth等高分辨率影像和当地的耕地实际调查信息对土地利用类型图进行了精度评估，共选取了200个随机位置的控制点，4期图像的精度评价见表2。

用户精度是从分类结果中任选一个样本，其与地面实际类型相一致的概率。生产者精度是从地面上任选一个样本，其与分类结果相同的概率。本次分类的结果生产者精度和用户精度都在70%以上，基本符合精度要求。kappa系数是用来评价图像的分类精度，主要用于精确性评价和图像的一致性判断。如果两幅图像差异很大，则其kappa系数小，当它大于0.8时，认为分类质量为极好。由表2可看出，4期影像的kappa系数均在0.7～0.81之间，其中1990年最低，2000年最高，总体来说kappa系数比较高，分类质量较好。总体分类精度是正确分类的像元数与总像元的比值，该值越接近100%，分类结果越精确。总体精度最低的1990年的影像，为81.04，其次是2018年和2010年的影像，这也说明了1990年的Landsat-5影像精度最低。4期影像的总计精度都大于80%说明分类结果比较准确。4期影像的总体精度和kappa系数都符合分类的精度要求，表明分类结果具有一定的可信度。

表2 1990-2018年河湟谷地影像精度评价结果 %

3.2 耕地面积序列变化

通过对近30 a来的Landsat遥感数据采用监督分类的最大似然分类法处理[16]，然后对分类结果进行聚类处理，将土地类型分为耕地、建设用地、水域、其他类型4类，由于本文主要是对耕地进行提取，故没有对草地和裸地进行细分。分别得到1990，2000，2010，2018年的河湟谷地土地利用类型图，并从土地利用类型中提取耕地类别得到耕地变化图(图1)。

图1 1990—2018年河湟谷地耕地变化特征

利用土地利用类型分类后的数据提取耕地，再计算耕地扩张动态度指数，可得到河湟谷地耕地面积时序变化[17](表3)。从表3可知，在1990—2010年，河湟谷地的耕地面积是处于一直在增长的趋势，而且在1990-2000年是低速增长的状态，而2000年以后处于一个快速扩张的状态，但是从2010年以后，耕地面积又出现了减少的势态，结合土地类型转移矩阵可知是由于伴随着人口的增加，有12%的农业用地转变为了建设用地，此外，青海省政府积极响应退耕还林还草政策，有17%的耕地转换为了林草地和裸地。

表3 1990-2018年河湟谷地耕地面积时序变化

根据耕地相对变化率公式可得出1990—2018年河湟谷地各县区耕地面积的相对变化率，进而利用ArcGIS 10.3得到1990—2018年河湟谷地耕地相对变化率图区位分析图[18](图2)。

(1)

式中：R为耕地面积的相对变化率；Ka、Kb代表整个研究区中某指定区域两个时间段的耕地面积；Ca和Cb分别代表整个研究区两个时间段的耕地面积。

从图2可看出，1990—2018年年来河湟谷地总体上耕地面积的变化幅度比较大，尤其是湟源县、贵德县和循化县3县相对变化率最大，而乐都县近30 a来变化幅度最小。在1990—2000年之间，河湟谷地大部分地区相对变化率都大于1，说明各个县区的耕地面积变化的幅度都大于整个研究区域的变化幅度，说明在1990—2000年之间，耕地变化还是比较大的；在2000—2010年之间，除贵德县相对变化率系数较大，耕地面积变化较大以外，其余11个县区耕地面积相对变化率较小；而在2010—2018年之间，大部分地区耕地面积相对变化系数都大于1，说明各地区耕地面积变化较大[19]。在这30 a的变化中，部分县的相对变化率出现了负数，这是因为这些县相对于整个区域的变化趋势相比，出现了负增长的态势。

图2 1990-2018年河湟谷地耕地相对变化率

3.3 河湟谷地耕地空间格局变化

3.3.1 耕地转换特征 根据解译得到的土地利用数据，可计算出耕地类型转移矩阵(表4)。

表4 1990-2018年河湟谷地耕地转移矩阵 %

从表4中可以看出，2000年耕地相较1990年转入了57.83%，转出了39.79%有10.42%的耕地转为建设用地；2010年耕地相较于2000年转入了53.94%，转出了31.74%，其中有12.06%的耕地转化为建设用地；2018年耕地相较于2010年转入了46.31%，转出了46.9%，有21.94%的耕地转化为建设用地，表明在2000—2018年之间河湟谷地地区城镇化进度比较快。近几年经济水平的提高造成了建设用地的不断增加，由1990年的3 480 km2增加至2018年的6 250 km2。

3.3.2 耕地质心变化 由耕地质心转移指数得出河湟谷地质心扩展变化转移指数(表5)，并将其质心转移状况在遥感影像中标出，得到河湟谷地耕地质心转移图(图3)。

表5 1990-2018年河湟谷地耕地质心扩展变化转移指数

图3 1990-2018年河湟谷地耕地重心转移特征

由表5，图3可知，在1990—2000年，河湟谷地耕地的质心往东南方向移动，移动距离为16 650 m，以正东方向为0度，逆时针为正方向，则角度为-30.74°。2000—2010年，质心移动的距离较小，为8 686 m，说明2000—2010年之间耕地变化不大，方向是向东北方向移动[20]； 2010—2018年，耕地质心偏移量较大，为16 287 m，方向接近于向西移动，这与湟水的流经方向基本一致，这也可以看出新增的耕地大多在西部，说明耕地扩张有向西移动的趋势[21]。

3.3.3 耕地空间分布叠加分析 将解译后的土地利用图利用ArcGIS建立色彩映射表，并进行类别提取，从而将耕地类别提取出来再将每10 a的数据进行叠加分析，得到不同时期耕地空间分布叠加分析结果(图4)。

从图4可以看出，2000年相对于1990年，新增耕地主要在研究区西北地区，即海晏县，其他地区耕地分布基本保持不变；2010年相对于2000年新增耕地主要在研究区的北部，并且变化的趋势是向着原有耕地边界周围的扩展，西南地区的耕地也有一定的增加，相较于2010年耕地的分布更加密集，而且此10 a的耕地增加的范围比较多，说明此10 a的耕地增加比较迅速；2018年相较于2010年的耕地主要向西北和东南发展，西北部的耕地范围进一步扩大，东南地区也有了新的耕地地块的出现。

图4 河湟谷地不同时期耕地空间分布叠加分析结果

3.4 耕地变化与地形的关系

(1) 不同海拔的耕地变化。耕地在不同海拔高度下的分布反映了高程对耕地具有很大影响[22]，一般来说耕地分布在海拔较低的平原地区，但河湟谷地地处黄河与湟水流域肥沃的三角地带，所以即使海拔较高处，农业也十分发达，是中国有名的高原农业区。将从影像中提取到的1990—2018年河湟谷地的耕地数据与DEM高程数据做叠加分析，即可统计出不同高程下的耕地分布[23]，结果见图5。

由图5可以看出，耕地随高程变化基本符合高斯正态分布，耕地大多都分布在海拔在2 500～3 000 m的范围内，由1990—2018年3个10 a的耕地高程分布情况可知，在1990—2010年耕地均值都是集中在海拔2 700 m的高度，但到了2018年耕地随高程变化的曲线有所转移，耕地的分布更向着高海拔地区转移，这是由于2010年后青海地区经济发展较为迅速，城镇建设占用了海拔较低的耕地，但耕地是保障人类生存的根本，所以需要开垦新的耕地来保障粮食需求，再加上全球变暖，近30 a来河湟谷地各区域呈现温度上升趋势(图6)，耕地质量提高，作物的种植海拔上限上升，耕地质量提高，在高海拔地区也有了适宜作物种植的温度环境，所以耕地向着高海拔地区开垦。由耕地面积随高程分布变化曲线可知，1990—2000年耕地随高程变化分布基本不变，曲线基本一致。但在2000—2010年的10 a里，耕地开垦速度加快，峰值增大，但是曲线形态不变，即只是耕地面积的增加，但耕地开垦的海拔不变，还是集中在海拔2 700 m的地区，到了2010年后，耕地才向着高海拔地区转移。

图5 1990-2018年河湟谷地耕地面积随高程变化分布

图6 1990-2018年河湟谷地不同海拔气温变化趋势

(2) 不同坡度的耕地变化分析。坡度对耕地分布有着重大影响，一般耕地都分布在坡度较缓的区域内[24]，从DEM高程数据中提取坡度信息，重分类后与耕地分布图相叠加，得到不同坡度的耕地变化情况(图7)。

图7 1990—2018年不同坡度耕地面积统计

由图7中可以看出，绝大部分耕地都分布在坡度在0°～16°的范围内，1990，2000，2010，2018年在此范围内耕地所占比例为分别为73%，64%，65%，63%。而且在0°～8°的坡度下耕地相对更多，而大于32°耕地面积很小；从时间上来看，在1990—2010年，更多新增耕地在8°～23°的范围内，坡度0°～8°的耕地基本没有变化或略有减少。一方面随着人口的增加需要开垦更多的耕地，低坡度地区已经不能满足新增耕地的需求，同时随着科技水平的提高，新的种植模式开发，向着坡度更陡的地区去开垦新的耕地。在2010—2018年可以看出在低坡度0°～8°的区域耕地继续增加，高坡度地区耕地很明显有减少的趋势，这是因为高坡度地区不宜大规模耕作，造成水土流失。在国家政策的影响下，当地政府采取退耕还林还草措施，以更好地保持水土，所以高坡度的耕地有减少的趋势。

4 讨论与结论

4.1 讨 论

(1) 本研究得到的结果与研究区前人相关成果较为一致。具体而言，青海省国土资源厅发布的数据表明：由于青海省实施的国家水土保持重点工程——坡耕地水土流失综合治理项目，2010年在同仁县、互助县、大通县等县兴修水平梯田28.2 km2，这与本文结果2010年后河湟谷地地区耕地向着高海拔地区扩展是一致的。张镱锂等[1,6]通过对1992—2015年已有LUCC数据与成果的再分析，得到了河湟谷地地区耕地扩张明显，82%的耕地分布在海拔2 400～4 400 m之间的河谷地区，这也与我们的结果基本一致。曾永年等[25]研究表明青海高原东部的耕地主要分布在坡度0°～15°之间，占该地区总耕地面积的77.08%，而坡度大于25°的耕地面积为仅占总耕地面积的4.50%。这与本文的结果(2010年坡度在0°～16°之间的耕地面积占总耕地面积的65%)基本一致。此外，与上述研究相比，本文时间序列较长，详细分析了1990—2018年河湟谷地耕地的时空演变特征，应用价值较高。

(2) 耕地变化受到多个因素的影响。就河湟谷地而言，耕地变化主要受到耕地保护政策、生态退耕工程、土地开发整理工程、重大建设工程、人口、城镇化进程等因素的影响。中国和青海省出台了一系列的耕地保护政策，这对河湟谷地耕地保护起到了促进的作用，也是河湟谷地耕地面积总体增加的主要原因。土地开发整理是河湟谷地耕地面积增加的主要来源，尤其是重大的土地开发整理工程，如青海省2010—2017年实施的“青海东部黄河谷地百万亩土地开发整理重大项目”，不但新增耕地136 km2，而且有效改善了黄河谷地的农业生产条件[26]。生态退耕工程、重大建设工程及城镇化进程是河湟谷地地区耕地减少的主要原因：青海省从1999年开始开展生态退耕工作，在2000—2005年期间，生态退耕面积1.47×105hm2，年均退耕面积2.45×104hm2，是青海省耕地减少的主要原因[25]。研究区内退耕还林/草区域主要分布于湟水及支流谷地周围坡度较大的区域以及拉脊山、达坂山山麓。2005年以后，生态退耕工程逐渐减少。由于“西部大开发”的实施，研究区内有很多重大工程建设项目，这些项目占用了大量耕地。此外，城镇化也是造成耕地减少的一个重要因素，研究区在青藏高原城镇空间格局演变中一直处于核心地位[27]，每年城镇建设都会占用不少耕地。

(3) 河湟谷地是青藏高原最大的农区，随着经济的不断发展，越来越多的耕地转化为了建设用地，建设用地和耕地之间的竞争更加激烈，很多优质耕地被建设占用，影响了这里的粮食安全。分析此区域耕地的空间分布特征和时空变化规律，可以得出不同时期的耕地开垦状态和耕地发展趋势，对未来利用和保护耕地具有重要的参考价值。但本研究在从遥感影像提取耕地的过程中，没有区分不同的耕地类型，也没有对耕地质量变化进行分析，但这些因素对于耕地资源稀少的青藏高原地区耕地的可持续发展无疑是重要的。而且，由于没有足够的数据，本文对于耕地变化的驱动因素，也没有做定量的相关性分析。今后应进一步探讨研究区耕地类别、耕地质量及耕地驱动因素对耕地变化的影响。

4.2 结 论

(1) 1990—2018年，随着城镇化水平的提高，河湟谷地部分耕地转化为了建设用地，尤其在2000—2010年这10 a里城镇化速率最高。耕地也是处于不断扩张中，也是在2000—2010年扩张速度最快，在这十年里部分耕地转化为了建设用地，新增耕地向着高海拔、高坡度地区开垦，2010年后耕地出现了减少趋势。

(2) 耕地的区位变化情况为，在2000年前耕地重心是向东移动的，而2000年后耕地重心偏移量向西北移动，这说明仅靠东边的耕地已经不能满足粮食需求，需要开垦更多的耕地，只能向气候环境热量条件稍差的西北部开垦。由于全气候变暖，降水增加，使得热量、水分条件变好，耕地质量提高，在高海拔地区也有了适宜作物种植的温度环境，所以近30 a来耕地向着高海拔地区开垦，近10 a在海拔3 000 m高度左右也出现了新增耕地。此外在一些坡度较陡的地区也出现了新增耕地，这也是因为气候变暖和人类扩张耕地的双重因素所造成的。在近10 a内坡耕地的增加有减缓的趋势，这也是受国家水土保持和退耕还林政策的影响。