赵宇豪，戎战磊，张玉凤，叶苗，蒋晖，赵传燕*

(1.兰州大学生命科学学院,甘肃 兰州730000;2.兰州大学草地农业科技学院，草地农业生态系统国家重点实验室，甘肃 兰州 730020)



近30年黑河流域草地变化分析及分布格局预测

赵宇豪1,2，戎战磊2，张玉凤1，叶苗1，蒋晖1，赵传燕2*

(1.兰州大学生命科学学院,甘肃 兰州730000;2.兰州大学草地农业科技学院，草地农业生态系统国家重点实验室，甘肃 兰州 730020)

草地是人类生存的基本资源和条件，在全球变化、全球循环和生态系统功能等方面起着非常关键的作用。本文基于GIS技术和1986、2000和2011年黑河流域草地分布数据，采用转移矩阵、草地动态变化度、空间动态度和CA-Markov模型等方法，分析了黑河流域近30年的草地变化，并对其未来分布格局做出预测。结果表明，(1) 黑河流域1986年草地面积为32997.71 km2，高覆被草地、中覆被草地和低覆被草地面积分别为4290.02、6371.80和22335.89 km2；2000年草地面积为23939.46 km2，高覆被草地、中覆被草地和低覆被草地面积分别为5508.78、5540.01和12890.67 km2；2011年草地面积为24272.50 km2，高覆被草地、中覆被草地和低覆被草地面积分别为5561.31、5715.01和12996.18 km2，近30年间黑河流域草地面积呈现先减后增的趋势。(2) 1986-2000年期间，黑河流域草地面积大幅减少，且退化严重，高覆被草地主要向中低覆被草地转换，中覆被草地主要向低覆被草地转换，低覆被草地主要向未利用地转换。(3) 2000-2011年期间，黑河流域草地面积小幅增加，各类型草地转化面积较少，保持小幅缓慢上升趋势。(4) 在现行趋势下，预测2022年黑河流域草地面积为24449.58 km2，高覆被草地、中覆被草地和低覆被草地面积分别为5560.28、5657.65和13231.65 km2，与2011年相比呈增加趋势，但下游地区草地面积和覆盖度依然较低。

黑河流域；草地；动态变化；CA-Markov模型；LUCC；GIS

草原是地球上分布最广的植被类型，总面积约为3.2×107km2，占陆地总面积的20%左右，我国是世界第二大草原国，拥有极为丰富的草地资源，面积约为 4×106km2，约占国土面积的41.7%，世界草地面积的12.5%[1-2]，主要分布在西北干旱、半干旱气候区以及青藏高原高寒气候区[3]，是我国陆地生态系统最大的天然屏障[4]。作为陆地生态系统中一个重要的子系统，草地生态系统在全球变化、全球循环和生态系统功能等方面起着非常关键的作用[6-10]，同时，草地资源是自然资源中可再生的一种，蕴藏着丰富的动植物资源，是农牧业经济的重要组成部分，不同的利用方式、管理方式对草地资源的可持续利用有重要影响[11-12]。

近年来，随着全球气候变化和人类活动干扰的增加，草地生态系统发生明显退化，现今草地退化已经成为一个全球性问题，而中国是世界上草地退化最为严重的国家之一[13]。由于长期以来人们对草地资源采取自然粗放的经营方式，目前，我国约90%的可利用天然草原呈现不同程度地退化，每年还以2×104km2的速度递增[13]，由此产生的区域生态问题越来越突出[14]。随着草地退化程度加剧，对草地退化过程进行动态监测、定量化描述与分析变得尤为必要。而传统的草地监测主要采用野外实地采样法，这种方法费时费力、效率低、成本高，还受诸多人为因素的限制，无法快速地在大范围内对草地生长状况和覆盖度变化进行动态监测。此外，这种方法常常用数量有限的监测站点代表一个区域，难以客观描述整个研究区内植被变化状况[15]。近年来，遥感(RS)和地理信息系统(GIS)等空间信息技术的迅速发展，为研究地球表面的动态变化提供了良好的手段，与传统方法相比,其具有准确、快速、覆盖面积大、价格低和周期性好等优点。实践表明，将RS和GIS等空间信息技术应用于草地退化研究是一种行之有效的方法，已经在全球不同地区草地植被动态监测和评价研究中广泛应用[16-22]。但是，目前人们多利用NDVI数据研究草地变化[23-28]，而在我国,草地的变化主要表现在草地向其他土地利用类型的转换[29]。

黑河流域内草原广布，面积约为9.5×104km2，占流域面积的60%以上，主要分为祁连山草地区，河西走廊中部草原区，河西走廊以北荒漠草场和额济纳三角洲荒漠草场。近年来随着气候变化以及人类活动的影响，黑河流域草地覆盖发生明显变化[30]。黑河上游和下游地区以牧业为主，中游地区以农业经营为主，地域经济类型呈现出山地牧业区—走廊绿洲农业区—高原牧业区组合特征，具有典型的农牧业分异和农牧交错特色，对草地资源的利用度较高，故而对黑河草地变化的研究不仅具有重要的理论意义，还具有深远的社会经济意义。因此本文采用1986、2000和2011年三期土地利用/土地覆盖变化(LUCC)数据，利用土地利用转移矩阵、草地动态变化度、空间变化度和CA-Markov模型等方法，对近30年黑河流域草地动态变化进行系统研究，主要对黑河流域草地变化特征和主要因素进行探讨，并预测其未来变化趋势，试图在更长时间尺度上展现黑河草地植被的时空变化特征，揭示黑河流域草地时空变化的特征和规律，为合理利用与保护草地资源提供重要理论依据，同时为黑河流域草地生态服务功能评价奠定基础。

1 材料与方法

图1 黑河流域土地利用及行政区划图Fig.1 Land use and administration districts of Heihe River Basin

1.1 研究区概况

黑河流域是我国第二大内陆河流域，横跨高山湿冷及半干旱区、温带和暖温带干旱区3个自然地理带，介于37°45′-42°40′ N，96°42′-102°4′ E之间(图1)，南北跨度大，流域内自然环境条件差异显著。其发源于南部祁连山区，全长821 km，流域面积约14.29万km2，主要覆盖三省区的11个县市，上游包括青海省祁连县和甘肃省肃南裕固族自治县，中游为甘肃嘉峪关、山丹、民乐、张掖、临泽、高台、肃南、金塔和肃州等市县，下游含甘肃肃北蒙古族自治县和内蒙古自治区额济纳旗。流域受大陆性气候和青藏高原气候的综合控制，降水少、蒸发量大，气候干燥，太阳辐射强烈，昼夜温差大，气候具有明显的东西差异和南北差异。上游多年平均降水量大于350 mm，降水量由东向西递减，年平均蒸发量为800～1200 mm，年平均气温-5～4 ℃。中游降水由东部的250 mm降低至西部的50 mm，蒸发量则从东部的2000 mm增加至西部的4000 mm以上，年平均气温2.8～7.6 ℃。下游属于干旱荒漠气候，年降水量小于50 mm，年蒸发量大于3500 mm，年平均气温8～10 ℃。巨大的自然条件差异，导致自东南向西北形成了链状山地—斑块绿洲—广域荒漠的生态系统，由此导致植被的地带性和非地带性的分布并存，也为草地类型的多样性奠定了环境基础，由北向南，草原类型依次为荒漠、草原化荒漠、荒漠化草原、山地森林草原、山地高寒草甸草原和高寒草甸。

1.2 土地利用/土地覆盖变化(LUCC)数据

土地利用/土地覆盖是生态环境中的一个敏感因子，受自然因素和人类活动的双重影响，在生态环境监测和可持续发展研究中占有非常重要的地位[31]。土地利用/土地覆盖变化(LUCC)不仅客观地记录了人类改变地球表面特征的空间格局，再现了地球表面景观的时空动态变化过程，而且与区域和全球环境变化、生物多样性减少，生态安全等问题息息相关[32-36]。目前，在国际地圈与生物圈计划(IGBP) 和国际全球环境变化的人文因素计划(IHDP)的推动下，国内外众多学者对其进行了大量研究[37-40]。

黑河流域LUCC数据[41-43]来源于黑河计划数据管理中心(http://westdc.westgis.ac.cn),所有数据采用统一的投影和坐标系，分辨率为1∶10万。采用国家土地分类系统，包括6个一级分类(耕地、林地、草地、水域、城镇和农村居民及工矿用地和未利用土地)，25个二级分类。其中草地分类体系为：高覆盖草地，指覆盖度>50%的天然草地、改良草地和割草地，此类草地一般水分条件较好，草被生长茂密，一般分布在山区、丘陵及河间滩地、沙地丘间等；中覆盖草地指覆盖度20%～50%的天然草地和改良草地，此类草地一般水分不足，草被较稀疏主要分布在较干燥地方；低覆盖草地指覆盖度 5%～20%的天然草地，此类草地水分缺乏，草被稀疏，牧业利用条件差 主要生长在干燥地方。草地向未利用地转换，以及高覆被草地向中低覆被草地转换，中覆被草地向低覆被草地转换意味着草地退化。

1.3 土地利用转移矩阵

土地利用转移矩阵是国内外常用的一种揭示土地利用/覆被类型间变化方向的方法[44-49]，其将土地利用变化的类型转移面积按矩阵的形式加以列出，不但包括静态的一定区域某时间点的各地类面积数据，而且能够细致地反映某一区域某一时段起止各土地类型面积之间相互转化的动态过程信息，具有丰富的统计学意义。在土地利用变化分析和地类模拟方面具有重要作用，具体形式如下。

A代表面积；n,m代表转移前后的土地利用类型数；ij(i,j=1,2,…,n)分别代表转移前与转移后的土地利用类型；Aij表示转移前的i地类转换成以后的j地类的面积。矩阵中的每一行元素代表转移前的i地类向转移后的各地类的流向信息，矩阵中的每一列元素代表转换以后的j地类面积从转移前的各地类的来源信息。为方便研究，通常转移前后采用相同的分类体系和分类精度，这样Aij的行数列数相同，为一个n阶方阵。本文运用ArcGIS对黑河流域3期土地利用数据进行统计和叠加分析，建立黑河流域1986-2000年，2000-2011年土地利用转移矩阵，分析近30年间草地转化情况。

1.4 草地动态变化度

单一土地利用类型动态度指某研究区一定时间范围内某种土地利用类型的数量变化[50]，本研究根据其定义，采用草地动态变化度来衡量草地变化的速率，计算公式如下。

LT=ΔP×Pa-1×T-1×100%

式中,LT为研究时段内草地动态变化度，ΔP为研究期内草地面积变化,Pa为研究期初期草地面积，T为研究时期长。动态变化度越大说明在研究期间内草地面积相对于自身的面积变化幅度越大。

1.5 草地空间动态度

空间动态度表示土地利用动态变化的空间过程和强烈程度[51],参考其定义,本研究使用草地空间动态度来衡量草地变化的强烈程度，计算公式可表示为：

式中，Li为研究时段内草地空间动态度；C1为研究期内草地转换为其他土地利用类型的面积，C2为研究期内非草地土地利用类型转换为草地的面积，C0为研究初期的草地面积，(C1+C2)/C0为研究期内草地的总变化量与研究初期草地面积的比值。

1.6 CA-Markov模型

马尔可夫链是用于土地利用变化建模的传统方法，系统的将来状态仅依赖于当前状态,通过创建从时期1到时期2的土地利用变化转换概率矩阵来实现未来预测,其表达式为。

S(t+1)=Pij×S(t)

式中：S(t+1)和S(t)分别表示t+1和t时刻土地利用系统的状态,Pij为状态转移矩阵。

元胞自动机(CA)是时间、空间和状态都离散的动力学模型。具有强大的空间运算能力，可以有效地模拟复杂系统的空间变化。CA 模型可用下式表示[52]。

S(t+1)=f(S(t),N)

式中：S为元胞有限、离散的状态集合；t、t+1为不同时刻；N为元胞的邻域；f为局部空间的元胞转化规则。

CA 和 Markov 模型两者均为时间、状态离散型的动力学模型，都被广泛应用于空间模拟中。 CA模型与Markov模型各具特色，但是又都存在局限，Markov模型侧重于数量方面的预测能力，而其空间参数比较弱，不能对空间上的变化进行表达，CA模型则具有较强的空间概念，可以模拟复杂的空间演变能力，但其局限于局部元胞间的相互作用，整体估算能力较差。将二者结合，既提高了土地利用类型转化的预测精度,又可以有效地模拟土地利用格局的空间变化，被广泛探讨和应用[53-56]。本文在软件IDRISI 17.0中使用CA-Markov模型,以2000和2011年土地利用类型为基准，对2022年草地分布进行预测分析。

1.7 统计年鉴

各省统计年鉴是一部按年度连续出版的大型统计资料书，其中含有大量的统计数据，可以全面反映省内的社会、经济和科技发展变化情况。本研究查阅了1987、2001和2012年青海省、甘肃省和内蒙古自治区统计年鉴，挑出包含在黑河流域中各区县的人口与牲畜存栏量情况，用以分析草地变化的驱动因素。

2 结果与分析

2.1 草地动态变化度和空间动态度

近30年间，黑河流域草地面积呈现先减后增的趋势，1986、2000和2011年的草地面积如图2所示。1986年,流域内草地面积为32997.71 km2，占流域面积的23.14%，其中高覆被草地、中覆被草地和低覆被草地面积分别为4290.02、6371.80和22335.89 km2；2000年流域内草地面积为23939.46 km2，占流域面积的16.79%，其中高覆被草地、中覆被草地和低覆被草地面积分别为5508.78、5540.01和12890.67 km2，2011年流域内草地面积为24272.50 km2，占流域面积的17.02%，其中高覆被草地、中覆被草地和低覆被草地面积分别为5561.31、5715.01和12996.18 km2。

1986-2011年，流域内草地面积减少8725.21 km2，动态变化度为-1.06%，空间动态度3.51%。其中高覆被草地面积增加1271.29 km2，动态变化度为1.19%，空间动态度为6.36%，表明其年均面积改变不大，但年均空间转换活动较为剧烈。中覆被草地面积减少656.79 km2，动态变化度为-0.41%，空间动态度为5.86%，年均面积改变程度和空间转换活动与高覆被草地类似。低覆被草地面积变化最大，减少了9339.71 km2，降幅达到了41.82%，动态变化度和空间动态度分别为-1.67%和4.86%，表明其年均面积改变和空间转换活动程度一般。1986-2000年，流域内草地面积减少9058.25 km2，动态变化度为-1.96%，空间动态度6.28%。高覆被草地面积增加1218.76 km2，动态变化度为2.03%，空间动态度为11.22%，表明其年均面积变化和空间转换活动均较为剧烈，中覆被草地面积减少831.80 km2，动态变化度为-0.93%，空间动态度为10.27%，表明其年均面积改变不大，但其年均空间转换活动较为剧烈，低覆被草地面积变化最大，减少了9445.26 km2，降幅高达57.71%，动态变化度为-3.02%，空间动态度8.69%，表明其年均面积变化较为剧烈，空间转换活动程度一般；2000-2011年间，所有草地类型面积均比2000年有所上升，但均未回升到1986年的水平，所有草地类型的空间动态度和动态变化度值均较小，表明近10年来，草地面积变化和空间转换较为平缓(图3)。

2.2 各区县草地变化

由表1可知，近30年期间，黑河上游各区县高覆被和低覆被草地面积均有增长，中覆被草地面积均有下降，其中肃南县的低覆被草地面积增加较多；在中游地区，玉门市和嘉峪关市境内无高覆被草地，临泽县和甘州区高覆被草地呈大幅下降趋势，其余区县呈增长趋势，其中民乐县增长较多。中覆被草地面积在中游各区县中有增有减，在甘州区、肃州区、高台县和金塔县中面积增加，在其余区县中面积减少。低覆被草地面积除在金塔县和山丹县有较大幅度的增加外，在其余区县面积减少；下游地区，肃北县境内依然无高覆被草地，而额济纳旗境内开始出现高覆被草地。中覆被草地在额济纳旗中面积增加，在肃北县中减少。下游各区县中的低覆被草地均有较大程度的减少。

1986-2000年期间，黑河上游除祁连县中覆被草地有所减少外，草地面积呈现增加趋势，而在中游，除山丹县、民乐县和金塔县草地面积有小幅增加以外，其余区县草地面积呈现下降趋势；黑河下游额济纳旗、肃北蒙古族自治县草地面积大幅缩减，其中额济纳旗草地面积由11766.71 km2减小到1840.44 km2，降幅达到84.36%，肃北蒙古族自治县草地面积由1324.37 km2降至70.23 km2，降幅高达94.70%。上述草地面积大幅下降的主要原因是因为其境内低覆被草地的剧烈减少。2000-2011年期间，黑河上游祁连县草地面积小幅下降，基本与2000年持平，其中高覆被草地与低覆被草地面积小幅降低，中覆被草地面积小幅升高，肃南县各类草地面积均有增加；中游区县中，甘州区、金塔县和玉门市草地面积有小幅上升，其余各区县草地面积均呈下降趋势，其中肃州区下降趋势比较明显；下游地区，各区县草地总面积均有所升高，且各类草地面积均有增加。

图2 1986,2000,2011年黑河流域草地分布Fig.2 1986,2000,2011 Heihe River Basin grassland distribution

图3 近30年各类草地面积空间变化分布Fig.3 Grassland types change spatial variation distribution over past 30 years

表1 1986,2000,2011年黑河流域主要区县草地面积Table 1 1986,2000,2011 grassland area of primary districts in Heihe River Basin km2

黑河上游地区因缺少祁连县数据，以肃南县数据为代表分析上游情况，由表2可知，1986-2011年期间，全流域内各区县人口和牲畜存栏量均有大幅增加。1986-2000年期间，上游地区人口增加，牲畜存栏量减少，中游地区，各区县人口和牲畜存栏量均呈增加趋势，下游地区各区县人口均呈增加趋势，除肃北蒙古族自治县牲畜存栏量略有下降之外，其余各县市牲畜存栏量增加，其中额济纳旗牲畜存栏量增幅巨大，2000年牲畜存栏量是1986年的2.9倍。2000-2011年期间，上游人口和牲畜存栏量呈增加趋势，中游地区，山丹县、民乐县、高台县、临泽县和玉门市人口减少，但其牲畜存栏量都有较大程度的增加，甘州区和肃州区人口和牲畜存栏量均有上升，且牲畜存栏量增幅较大。下游地区各县市人口增加牲畜存栏量降低，其中额济纳旗牲畜存栏量减少较多，减少了一半左右。

2.3 草地转化

由表3可知，1986-2011年期间，高覆被草地主要向中低覆被草地转化，剩余转出面积由大到小依次为未利用地、水域、林地、耕地和居民地，转入面积中，最多的为林地，面积为1933.04 km2。其次依次为中覆被草地，未利用地、低覆被草地、水域、耕地和居民地；中覆被草地向低覆被草地和未利用地的转化较多，转出面积分别为2403.32和1175.09 km2，剩余转换面积依次为高覆被草地、水域、耕地、林地和居民地，转入面积中，最多的为低覆被草地，面积为1396 km2；其次依次为未利用地、林地、高覆被草地、水域、耕地和居民地；低覆被草地绝大部分面积向未利用地转换，草地变为未利用地，是指其向沙地、戈壁、沼泽地、裸土地、裸岩等转换，表明低覆被草地退化明显，剩余转换面积为中覆被草地、耕地、水域、高覆被草地、林地和居民地，转入面积中最多的未利用地，其次依次为中覆被草地、林地、高覆被草地、耕地、水域和居民地。

由表4可知，1986-2000年期间，高覆被草地主要向中低覆被草地转化，剩余转换面积由大到小依次为未利用地、水域、林地、耕地和居民地，转入面积中，最多的为林地，面积为1903.39 km2。其次依次为中覆被草地，未利用地、低覆被草地、水域、耕地和居民地；中覆被草主要向低覆被草地转化，转化面积高达2347.42 km2，同时未变化和向未利用地转换的也较多，剩余转换面积依次为高覆被草地、水域、耕地、林地和居民地，转入面积中，最多的为低覆被草地，面积为1364.48 km2；其次依次为林地，未利用地、高覆被草地、水域、耕地和居民地；低覆被草地绝大部分面积向未利用地转换，剩余转换面积为中覆被草地、水域、耕地、高覆被草地、林地和居民地，转入面积中最多的未利用地，其次依次为中覆被草地、林地、高覆被草地、耕地、水域和居民地。

表2 1986-2011黑河主要区县人口和牲畜存栏量Table 2 1986-2011 population and grazing capacity of primary districts in Heihe River Basin

注：祁连县无1986年统计数据，用1987年数据代替。牲畜存栏量为当年年末大牲畜和羊的存栏量，大牲畜包括牛，马，驴，骡和骆驼等。
Note:The 1986 statistical data lacked in Qilian county,which was replaced by the 1987 statistical data.Cattle stocks refer to the stocks of sheep and large livestock which includes cows,horses,donkeys,mules and camels,etc.at the end of year.

表 3 1986-2011年黑河流域土地利用转移矩阵Table 3 Land use transfer matrix of Heihe River Basin in 1986-2011 years km2

表4 1986-2000年黑河流域土地利用转移矩阵Table 4 Land use transfer matrix of Heihe River Basin in 1986-2000 years km2

由表5可知，2000-2011年期间，绝大部分高覆被草地没有变化，剩余主要向中低覆被草地转换，向水域、林地、耕地和居民地转换较少，转入面积中，主要为中低覆被草地和未利用地、其余土地类型转入较少。中覆被草地绝大部分没有变化，其主要向低覆被草地和高覆被草地转换，剩余转换面积由大到小依次为未利用地、耕地、林地、水域和居民地，转入面积中未利用地最多，其次主要为低覆被草地和高覆被草地，其余土地利用类型转入中覆被草地较少；低覆被草地绝大部分未变化，转换面积由大到小依次为未利用地、耕地、中覆被草地、高覆被草地、林地、水域和居民地。

表5 2000-2011年黑河流域土地利用转移矩阵Table 5 Land use transfer matrix of Heihe River Basin in 2000-2011 years km2

2.4 草地分布格局

运用CA-Markov模型，预测黑河流域2022年草地面积，结果如图4所示，草地总面积为24449.58 km2，高覆被草地、中覆被草地面积和低覆被草地面积分别为5560.28、5657.65和13231.65 km2，其草地总面积大于2000和2011年面积，但不及1986年草地总面积。2022年高覆被草地面积，中覆被草地面积低于1986和2011年的面积，较2000年草地面积略有升高，低覆被草地面积较2000和2011年有所升高，但是远不及1986年的水平。表明在此趋势下，2022年黑河流域草地面积呈现总体好转，局部下降的趋势。

图4 黑河流域2022年草地分布Fig.4 Grassland distribution of Heihe River Basin in 2022 years

3 讨论

1986-2000年期间，黑河流域草地总面积减少，但其高覆被草地面积增加，这主要可能因为气候、政策和人为因素的驱动。高覆被草地主要位于黑河上游地区，1987年以来西北气候开始由暖干向暖湿转型，黑河上游位于气候显著转型和轻度转型区，气温升高，降水量增加[57-58]，在高寒草甸地区，温度的升高有利于草地的生长[59]，且1988年祁连山国家级自然保护区成立，黑河上游大部分地区位于保护区内，其实施了一系列生态保护政策的实施，如封山育林、退耕还林还草工程等，使人类对黑河上游草地干扰程度大幅降低。期间虽然上游区县人口有所增加，但牲畜存栏量有显著降低，所以这些因素共同造成了1986-2000年间高覆被草地面积增加的情况。中覆被草地面积降低比例不大，主要向低覆被草地转换，减少区域与周伟等[30]的研究结果相似，这可能是因为上游地区非保护区域内过度放牧所致[60-62]。低覆被草地面积减少较多，这主要是因为额济纳旗境内低覆被草地的大量减少，这与邹亚荣等[29]的研究结果相似，主要的驱动因素可能是由于人为的干扰。黑河流域中上游拦坝蓄水，中游地区工农业用水激增[63-64]，导致河流下泄量日趋减少，使下游河床、湖泊干涸[65]，地下水位下降，植被衰退[66]。其次下游地区人口和牲畜存栏量大量增加，额济纳旗尤为明显，其牲畜存栏量增加了将近3倍，地表原始植被在牲畜的过度践踏和啃食下，丧失了再生能力,逐渐退化并且消失。同时，下游地区草原有丰富的野生食用和药用植物资源，因其较高的经济价值使人们过度挖采[67]，使本来很脆弱的草原生态系统又遭到严重破坏。

2000-2011年期间，黑河流域草地总面积和各类草地面积均小幅上升，这主要是由政策因素驱动导致。2000年以来国家在上游地区实行的一系列生态治理政策[39]，中游地区草地面积呈现下降趋势，据转移矩阵显示期间低覆被草地转向耕地的占比有所增加，这可能是因为中游地区人口稠密，是我国西北主要的商品粮基地，较大的农垦压力占用了一部分草地，同时据年鉴统计，2000-2011年期间中游地区牲畜存栏量大幅增加，过度放牧也使得草地面积有所减小。下游地区草地面积有所上升，这可能是因为2000年以来，国家相继实施了退耕还林还草工程、退牧还草工程以及草原生态保护补助奖励政策，期间额济纳旗的牲畜存栏量下降了近一半，这极大地减轻了草地的负担，同时国家在黑河中游启动了水量统一调度方案[39]，使得下游地区的水资源短缺问题得到了一定程度的缓解，因此，草原面积有一定程度的回升。

黑河流域以农牧业为主,草地面积与盖度的减少会对其经济带来制约[48]，根据目前的趋势发展，黑河流域的草地面积会逐渐得到恢复，但其恢复程度还不够，尤其是下游地区，长时间的植被退化，使得下游沙漠化严重，植被恢复困难，因此应该在保证中上游正常生产工作的同时，适当削减水坝数量和用水量，严禁过度建坝蓄水，尽量加大下游地区的供水量，使下游地区年过水量增大，令自然植被得以恢复，同时加大草原生态保护补助奖励政策，尽量降低放牧对草原的破坏。

4 结论

(1) 近30年间黑河流域草地呈现先减后增的趋势，但远未到达1986年的草地面积，这主要是由低覆盖草地大量减少导致。

(2) 1986-2000年期间，黑河流域草地面积大幅减少，且退化严重，高覆被草地主要向中低覆被草地转换，中覆被草地主要向低覆被草地转换，低覆被草地主要向未利用地转换。

(3) 2000-2011年期间，黑河流域草地面积小幅增加，各类型草地转出面积较少，保持小幅缓慢上升趋势。

(4) 在现行趋势下，未来黑河流域草地面积呈增加趋势，但下游地区草地面积和覆盖度依然较低，国家应制定相关政策加大下游地区草地的恢复和保护力度。

References：

[1] Fan J,Zhong H,Harris W,etal.Carbon storage in the grasslands of China based on field measurements of above and below-ground biomass.Climatic Change,2008,86(3):375-396.

[2] Huang Y,Sun W J,Zhang W,etal.Study on carbon budget in grassland of china:progress and perspectives.Quaternary Sciences,2010,30(3):456-465.黄耀,孙文娟,张稳,等.中国草地碳收支研究与展望.第四纪研究,2010,30(3):456-465.

[3] Zhou W,Gang C C,Li J L,etal.Spatial-temporal dynamics of grassland coverage and its response to climate change in China during 1982-2010.Acta Geographica Sinica,2014,69(1):15-30.周伟,刚成诚,李建龙,等.1982-2010年中国草地覆盖度的时空动态及其对气候变化的响应.地理学报,2014,69(1):15-30.

[4] Ren J Z,Liang T G,Lin H L,etal.Study on grassland’ s responses to global climate change and its carbon sequestration potentials.Acta Prataculturae Sinica,2011,20(2):1-22.任继周,梁天刚,林慧龙,等.草地对全球气候变化的响应及其碳汇潜势研究.草业学报,2011,20(2):1-22.

[5] Liu X Y,Long R J,Shang Z H,etal.Evaluation method of ecological services function and their value for grassland ecosystems.Acta Prataculturae Sinica,2011,20(1):167-174.刘兴元,龙瑞军,尚占环.草地生态系统服务功能及其价值评估方法研究.草业学报,2011,20(1):167-174.

[6] Liu X N,Sun J L,Zhang D G,etal.A study on the community structure and plant diversity of alpine meadow under different degrees of degradation in the Eastern Qilian Mountains.Acta Prataculturae Sinica,2008,17(4):1-11.柳小妮,孙九林,张德罡,等.东祁连山不同退化阶段高寒草甸群落结构与植物多样性特征研究.草业学报,2008,17(4):1-11.

[7] Liang S,Peng S S,Lin X,etal.NDVI-based spatial-temporal change in grassland growth of China from 1982 to 2010.Acta Scientiarum Naturalium Universitatis Pekinensis,2013,49(2):311-320.梁爽,彭书时,林鑫,等.1982-2010年全国草地生长时空变化.北京大学学报(自然科学版),2013,49(2):311-320.

[8] Scurlock J,Hall D O.The global carbon sink:a grassland perspective.Global Change Biology,1998,4(2):229-233.

[9] Scurlock J,Johnson K,Olson R J.Estimating net primary productivity from grassland biomass dynamics measurements.Global Change Biology,2002,8(8):736-753.

[10] Yan Z Q,Qi Y C,Dong Y S,etal.Nitrogen cycling in grassland ecosystems in response to climate change and human activities.Acta Prataculturae Sinica,2014,23(6):279-292.闫钟清,齐玉春,董云社,等.草地生态系统氮循环关键过程对全球变化及人类活动的响应与机制.草业学报,2014,23(6):279-292.

[11] A Y G L.Dynamic Changes Analysis and Value of Ecosystem Services Assessment of Grassland of the Western Jilin[D].Changchun:Jilin University,2009.阿依古丽.吉林西部草地动态变化分析与生态系统服务功能价值评价[D].长春:吉林大学,2009.

[12] Dai S P,Zhang B,Wang H J,etal.Analysis on the spatio-temporal variation of grassland cover using SPOT NDVI in Qilian Mountains.Progress in Geography,2010,29(9):1075-1080.戴声佩,张勃,王海军,等.基于SPOT NDVI的祁连山草地植被覆盖时空变化趋势分析.地理科学进展,2010,29(9):1075-1080.

[13] Xiao T,Wang C Z,Feng M,etal.Dynamic characteristic of vegetation coverage in the Three-River Source Region from 2000 to 2011.Acta Agrestia Sincia,2014,22(1):39-45.肖桐,王昌佐,冯敏,等.2000-2011年青海三江源地区草地覆盖度的动态变化特征.草地学报,2014,22(1):39-45.

[14] Li Y F.The Zoning Evaluation of Grassland Change in China and its Effect on Ecological Security[D].Wuhan:Hubei University,2014.李燕飞.中国草地变化及其生态安全保障功能的分区评价[D].武汉:湖北大学,2014.

[15] Yu H.Dynamics of Grassland Growth and Its Response to Climate Change on Tibetan Plateau[D].Lanzhou:Lanzhou University,2013.于惠.青藏高原草地变化及其对气候的响应[D].兰州:兰州大学,2013.

[16] Ma L Y,Cui X,Feng Q S,etal.Dynamic changes of grassland vegetation coverage from 2001 to 2011 in Gannan Prefecture.Acta Prataculturae Sinica,2014,23(4):1-9.马琳雅,崔霞,冯琦胜,等.2001-2011年甘南草地植被覆盖度动态变化分析.草业学报,2014,23(4):1-9.

[17] Yang H F,Gang C C,Mu S J,etal.Analysis of the spatio-temporal variation in net primary productivity of grassland during the past 10 years in Xinjiang.Acta Prataculturae Sinica,2014,23(3):39-50.杨红飞,刚成诚,穆少杰,等.近10年新疆草地生态系统净初级生产力及其时空格局变化研究.草业学报,2014,23(3):39-50.

[18] Liu J Y,Xu X L,Shao Q Q.The spatial and temporal characteristics of grassland degradation in the Three- River Headwaters Region in Qinghai Province.Acta Geographica Sinica,2008,63(4):364-376.刘纪远,徐新良,邵全琴.近30年来青海三江源地区草地退化的时空特征.地理学报,2008,63(4):364-376.

[19] Wang Y,Xia W T,Liang T G,etal.Spatial and temporal dynamic changes of net primary product based on MODIS vegetation index in Gannan grassland.Acta Prataculturae Sinica,2010,19(1):201-210.王莺,夏文韬,梁天刚,等.基于MODIS植被指数的甘南草地净初级生产力时空变化研究.草业学报,2010,19(1):201-210.

[20] Luo C F,Xu C J,You H Y,etal.Analysis on grassland degradation in Qinghai Lake Basin during 2000—2010.Acta Ecologica Sinica,2013,33(14):4450-4459.骆成凤,许长军,游浩妍,等.2000-2010年青海湖流域草地退化状况时空分析.生态学报,2013,33(14):4450-4459.

[21] Kirsten M B,Braden C O,Julian J P.Disturbance analyses of forests and grasslands with MODIS and Landsat in New Zealand.International Journal of Applied Earth Observation,2016,45:42-54.

[22] Xie H L,Li X B,Zhang Y T,etal.Spatial difference of regional grass changes based on ESDA at county level in Beijing-Tianjin-Hebei Area.Journal of Natural Resources,2012,27(7):1224-1232.谢花林,李秀彬,张燕婷,等.基于ESDA的京津冀地区草地变化空间分异.自然资源学报,2012,27(7):1224-1232.

[23] Ding M J,Zhang Y L,Liu L S,etal.Temporal and spatial distribution of grassland coverage change in Tibetan Plateau since 1982.Journal of Natural Resources,2010,25(12):2114-2122.丁明军,张镱锂,刘林山,等.1982-2009年青藏高原草地覆盖度时空变化特征.自然资源学报,2010,25(12):2114-2122.

[24] Zhang Q Y,Wu S H,Zhao D S,etal.Responses of growing season vegetation changes to climatic factors in Inner Mongolia grassland.Journal of Natural Resources,2013,28(5):754-764.张清雨,吴绍洪,赵东升,等.内蒙古草地生长季植被变化对气候因子的响应.自然资源学报,2013,28(5):754-764.

[25] Guo G M,Xie G D,Zhen L.The relationship between NDVI change and precipitation in Guyuan City.Resource Science,2007,29(2):178-182.郭广猛,谢高地,甄霖.泾河上游固原地区的NDVI变化与降水的相关性研究.资源科学,2007,29(2):178-182.

[26] Xin X P,Zhang B H,Li G,etal.Variation in spatial pattern of grassland biomass in China from 1982 to 2003.Journal of Natural Resources,2009,24(9):1582-1592.辛晓平,张保辉,李刚,等.1982-2003年中国草地生物量时空格局变化研究.自然资源学报,2009,24(9):1582-1592.

[27] Jia W X,Liu M,Yang Y H,etal.Estimation and uncertainty analyses of grassland biomass in Northern China:Comparison of multiple remote sensing data sources and modeling approaches.Ecological Indicators,2016,60:1031-1040.

[28] Zhang G L,Xu X L,Zhou C P,etal.Responses of vegetation changes to climatic variations in Hulun Buir grassland in past 30 years.Acta Geographica Sinica,2011,66(1):47-58.张戈丽,徐兴良,周才平,等.近30年来呼伦贝尔地区草地植被变化对气候变化的响应.地理学报,2011,66(1):47-58.

[29] Zou Y R,Zhang Z X,Zhou Q B,etal.Spatial pattern and Its Analysis of China’ s grassland change in recent ten years using remote sensing and GIS.Journal pf Remote Sensing,2003,7(5):428-432.邹亚荣,张增祥,周全斌,等.遥感与GIS支持下近十年中国草地变化空间格局分析.遥感学报,2003,7(5):428-432.

[30] Zhou W,Wang Q,Zhang C B,etal.Spatiotemporal variation of grassland vegetation NDVI in the middle and upper reaches of the Hei River and its response to climatic factors.Acta Prataculturae Sinica,2013,22(1):138-147.周伟,王倩,章超斌,等.黑河中上游草地NDVI时空变化规律及其对气候因子的响应分析.草业学报,2013,22(1):138-147.

[31] Zheng L Y.Study on Dynamic Change of NPP and Grassland Change in Northern Tibet Based on Remote Sensing and Biological Processes Model BEPS[D].Beijing:Chinese Academy of Meteorological Sciences,2006.郑凌云.基于卫星遥感与BEPS生态模式的藏北草地变化及NPP动态研究[D].北京:中国气象科学研究院,2006.

[32] Shi P J,Song C Q,Jing G F.Strengthening the study of land use/cover change and its impact on eco-environmental security——The trend of the study of the dynamics of Human-nature system based on “global change open science conference 2001” in Amsterdam,Netherlands.Advance in Earth Sciences,2002,17(2):161-168.史培军,宋长青,景贵飞.加强我国土地利用/覆被变化及其对生态环境安全影响的研究——从荷兰“全球变化开放科学会议” 看人地系统动力学研究的发展趋势.地理科学进展,2002,17(2):161-168.

[33] Hou L S,Cai Y L.An essential analysis and review on land use/cover change research.Progress in Geography,2004,23(6):96-104.后立胜,蔡运龙.土地利用/覆被变化研究的实质分析与进展评述.地理科学进展,2004,23(6):96-104.

[34] Du L T.Study on land use/ cover change in mountain areas of southern Ningxia based on RS technology.Progress in Geography,2006,25(6):94-101.杜灵通.基于遥感技术的宁夏南部山区LUCC 研究.地理科学进展,2006,25(6):94-101.

[35] Chen B M,Liu X W,Yang H.Review of most recent progresses of study on land use and land cover change.Progress in Geography,2003,22(1):22-29.陈百明,刘新卫,杨红.LUCC 研究的最新进展评述.地理科学进展,2003,22(1):22-29.

[36] Liu J H,Gao J X,Geng B,etal.Study on the dynamic change of land use and landscape Pattern in the farming-pastoral region of Northern China.Research of Environmental Sciences,2007,20(5):148-154.刘军会,高吉喜,耿斌,等.北方农牧交错带土地利用及景观格局变化特征.环境科学研究,2007,20(5):148-154.

[37] Zhang K,Yu Z,Li X,etal.Land use change and land degradation in China from 1991 to 2001.Land Degradation and Development,2007,18(2):209-219.

[38] Thielen D R,Jose J J,Montes R A,etal.Assessment of land use changes on woody cover and landscape fragmentation in the Orinoco savannas using fractal distributions.Ecological Indicators,2008,8(2008):224-238.

[39] Liu J W,Jin T T,Liu G H,etal.Changes in land use and soil and water conservation of the upper and middle reaches of Heihe river basin during 2000—2010.Acta Ecologica Sinica,2014,34(23):7013-7025.刘金巍,靳甜甜,刘国华,等.黑河中上游地区2000-2010年土地利用变化及水土保持服务功能.生态学报,2014,34(23):7013-7025.

[40] Wang F,Gao J E,Shao H,etal.Response of ecosystem service values to land use change based on GIS and ecological compensation in Loess Plateau.Science of Soil and Water Conservation,2013,11(1):25-31.王飞,高建恩,邵辉,等.基于GIS 的黄土高原生态系统服务价值对土地利用变化的响应及生态补偿.中国水土保持科学,2013,11(1):25-31.

[41] Wang Y M,Yan C Z,Wang J H.Landuse/Landcover data of the Heihe River Basin in 1986[C].Heihe Plan Science Data Center,2011.doi:10.3972/heihe.021.2013.db.王一谋,颜长珍,王建华.黑河流域1980年代末土地利用/土地覆盖数据集[C].黑河计划数据管理中心,2011.doi:10.3972/heihe.021.2013.db.

[42] Wang J H,Zhao J,Hu X L,etal.Landuse/Landcover data of the Heihe River Basin in 2000[C].Heihe Plan Science Data Center,2014.doi:10.3972/heihe.039.2014.db.王建华,赵军,胡晓利,等.2012年新编黑河流域2000年土地利用/土地覆盖数据集[C].黑河计划数据管理中心,2014.doi:10.3972/heihe.039.2014.db.

[43] Wang J H,Ge C M,Zhao J,etal.Landuse/Landcover data of the Heihe River Basin in 2011[C].Heihe Plan Science Data Center,2014.doi:10.3972/heihe.093.2014.db.王建华,盖春梅,赵军,等.黑河流域2011年土地利用/土地覆被数据集[C].黑河计划数据管理中心,2014.doi:10.3972/heihe.093.2014.db.

[44] Liu R,Zhu D L.Methods for detecting land use changes based on the land use transition matrix.Resources Science,2010,32(8):1544-1550.刘瑞,朱道林.基于转移矩阵的土地利用变化信息挖掘方法探讨.资源科学,2010,32(8):1544-1550.

[45] Qiao W F,Sheng Y H,Fang B,etal.Land use change information mining in highly urbanized area based on transfer matrix:A case study of Suzhou,Jiangsu Province.Geographical Research,2013,32(8):1497-1507.乔伟峰,盛业华,方斌,等.基于转移矩阵的高度城市化区域土地利用演变信息挖掘——以江苏省苏州市为例.地理研究,2013,32(8):1497-1507.

[46] Yue D X,Du J,Liu J Y,etal.Spatio-temporal analysis of ecological carrying capacity in Jinghe Watershed based on Remote Sensing and Transfer Matrix.Acta Ecologica Sinica,2011,31(9):2550-2558.岳东霞,杜军,刘俊艳,等.基于RS 和转移矩阵的泾河流域生态承载力时空动态评价.生态学报,2011,31(9):2550-2558.

[47] Hao H M,Hao Y L,Ren Z Y.Analysis on dynamic and pattern of land use/cover change in Guanzhong Area.Scientia Agricultura Sinica,2011,44(21):4525-4536.郝慧梅,郝永利,任志远.近20年关中地区土地利用/覆盖变化动态与格局.中国农业科学,2011,44(21):4525-4536.

[48] Xiang B,Zhang Z X,Buhorsor.A study on arable land dynamic changes and management in ecotone of farming,forestry and animal husbandry by using remotely sensed data -taking Horqin Right Wing Front Banner in Inner Mongolia as an example.Resources Science,1999,21(6):9-12.香宝,张增祥,布和敖斯尔.基于遥感的农林牧交错带耕地动态变化及管理研究——以内蒙古科右前旗为例.资源科学,1999,21(6):9-12.

[49] Bai G C,Xia J G,Wang C Q,etal.Land use transition in dongpo district,meishan City in China based on a tendentious model for analyzing spatial transition of land use types.Resources Science,2009,31(10):1793-1799.白根川,夏建国,王昌全,等.基于地类空间转化趋势模型的眉山市东坡区土地利用转化分析.资源科学,2009,31(10):1793-1799.

[50] Wnag X L,Bao Y H.Study on the methods of land use dynamic change research.Progress in Geography,1999,18(1):81-87.王秀兰,包玉海.土地利用动态变化研究方法探讨.地理科学进展,1999,18(1):81-87.

[51] Liu S H,He S J.A spatial analysis model for measuring the rate of land use change.Journal of Natural Resources,2002,17(5):533-540.刘盛和,何书金.土地利用动态变化的空间分析测算模型.自然资源学报,2002,17(5):533-540.

[52] He D,Zhou J,Gao W,etal.An integrated CA-Markov model for dynamic simulation of land use change in Lake Dianchi Watershed.Acta Scientiarum Naturalium Universitatis Pekinensis,2014,50(6):1095-1105.何丹,周璟,高伟,等.基于 CA-Markov 模型的滇池流域土地利用变化动态模拟研究.北京大学学报(自然科学版),2014,50(6):1095-1105.

[53] Cheng G,Zhang Z L,Lv J S.Landscape pattern analysis and dynamic prediction of Sanchuan basin in East China based on CA-Markov model.Chinese Journal of Ecology,2013,32(4):999-1005.程刚,张祖陆,吕建树.基于 CA-Markov 模型的三川流域景观格局分析及动态预测.生态学杂志,2013,32(4):999-1005.

[54] Guan D J,Li H F,Inohae T,etal.Modeling urban land use change by the integration of cellular automaton and Markov model.Ecological Modelling,2011,222(s 20-22):3761-3772.

[55] Tong S,Sun Y,Thushara R,etal.Predicting plausible impacts of sets of climate and land use change scenarios on water resources.Applied Geography,2012,32(2):477-489.

[56] Wang Y S,Yu X X,He K N,etal.Dynamic simulation of land use change in Jihe watershed based on CA-Markov model.Transactions of the CSAE,2011,27(12):330-336.王友生,余新晓,贺康宁,等.基于CA-Markov 模型的藉河流域土地利用变化动态模拟.农业工程学报,2011,27(12):330-336.

[57] Shi Y,Shen Y,Kang E,etal.Recent and future climate change in northwest China.Climatic Change,2007,80(3):379-393.

[58] Cheng Y,Xu D X,Guo N.Analysis on the vegetation change in the Qilian Mountains since recent 22 years.Arid Zone Research,2008,25(6):772-777.程瑛,徐殿祥,郭铌.近20年来祁连山区植被变化特征分析.干旱区研究,2008,25(6):772-777.

[59] Guo N,Zhu Y J,Wang J M,etal.The relationship between NDVI and climate elements for 22 years in different vegetation areas of northwest China.Journal of Plant Ecology (Chinese Version),2008,32(2):319-327.郭铌,朱燕君,王介民,等.近22年来西北不同类型植被NDVI变化与气候因子的关系.植物生态学报,2008,32(2):319-327.

[60] Zhang H,Shen W T,Wang Y S,etal.Study on grassland capacity in the Heihe River Basin.Journal of Natural Resources,2005,20(4):514-521.张慧,沈谓涛,王延松,等.黑河流域草地承载力研究.自然资源学报,2005,20(4):514-521.

[61] Lan Y C,Sun B M,Ding Y J,etal.Studies on ecological environment changes of Heihe River Basin and its influence factors.Journal of Arid Land Resources and Environment,2004,18(2):32-39.蓝永超,孙保沐,丁永建,等.黑河流域生态环境变化及其影响因素分析.干旱区资源与环境,2004,18(2):32-39.

[62] Chang J,Wang G X,Wang Y B.Driving factors of land use change in the Heihe River:case study of Zhangye prefecture.Journal of Glaciology and Geocryology,2005,27(1):117-123.常娟,王根绪,王一博.黑河流域土地利用变化的影响因素-以张掖地区为例.冰川冻土,2005,27(1):117-123.

[63] Hao X W,Zhang Y F,Kang J F.Study on transition and desertification of natural geographical landscape pattern in Heihe River Basin.Journal of Anhui Agri,2015,43(7):265-268.郝兴旺,张永福,康建锋.黑河流域自然地理景观格局变迁与荒漠化研究.安徽农业科学,2015,43(7):265-268.

[64] Wang G X,Cheng G D.Land Desertification status and developing trend in the Heihe River Basin.Journal of Desert Research,1999,12(4):368-374.王根绪,程国栋.黑河流域土地荒漠化及其变化趋势.中国沙漠,1999,12(4):368-374.

[65] Xiao S C,Xiao H L,Lan Y C,etal.Water issues and integrated water resource management in Heihe River Basin in recent 50 years.Journal of Desert Research,2001,31(2):529-535.肖生春,肖洪浪,蓝永超,等 近50 a 来黑河流域水资源问题与流域集成管理.中国沙漠,2001,31(2):529-535.

[66] Li S,Li F,Sun W,etal.Modern desertification process in Ejina Oasis and its dynamic mechanism.Scientia Geographica Sinica,2004,24(1):61-67.李森,李凡,孙武,等.黑河下游额济纳绿洲现代荒漠化过程及其驱动机制.地理科学,2004,24(1):61-67.

[67] Qi S Z,Wang T.Current status and causes of land desertification in middle and lower reaches of Heihe River Basin.Journal of Soil and Water Conservation,2003,17(4):98-101,109.齐善忠,王涛.黑河流域中下游地区土地沙漠化现状及其原因分析.水土保持学报,2003,17(4):98-101,109.

Analysis of change in grassland area in the Heihe River basin over the past 30 years and prediction of future trends

ZHAO Yu-Hao1,2,RONG Zhan-Lei2,ZHANG Yu-Feng1,YE Miao1,JIANG Hui1,ZHAO Chuan-Yan2*

1.SchoolofLifeSciences,LanzhouUniversity,Lanzhou730000,China;2.StateKeyLaboratoryofGrasslandAgro-ecosystems,CollegeofPastoralAgricultureScienceandTechnology,LanzhouUniversity,Lanzhou730020,China

Grassland,is one of the most widespread vegetation types in the world,and accounts for nearly 20% of the global land surface.China has abundant grassland resources,and the second largest grassland area of any country.However,in recent years,China’s grassland ecosystem has suffered dramatic impairment through human activities and climate change.Grassland in the Heihe River Basin provides necessary resources for the livelihood of local communities.It is therefore important to monitor grassland change.In this study,we assessed historical grassland cover in Heihe River Basin using GIS technology and evaluated land use and land cover change (LUCC) since 1986.We also predicted the future grassland distribution pattern using a CA-Markov Model and land-use transfer matrix,reflecting change in spatial distribution and metabolic indices of grassland.In 1986,the total area of grassland in the Heihe River Basin was 32998 km2,of which 4290 km2were classed as having high coverage,6372 km2were classed as having medium coverage and 22336 km2were classed as having low coverage.In 2000,the grassland area was 23939 km2,and areas with high,medium and low coverage were,respectively,5509,5540 and 12891 km2.In 2011,the grassland area was 24273 km2,and areas with high medium and low coverage were,respectively,5561,5715 and 12996 km2.During the period 1986-2000,the area of grassland in the Heihe River Basin was considerably reduced and severe degeneration occurred.During this period,there was substantial transformation of grassland from high to medium coverage,from medium to low coverage,and from low coverage to unutilized land.During the period 2000-2011,the grassland area in Heihe River Basin increased slightly and transition between the three coverage classes was small.Finally,in 2022 the predicted grassland area of the Heihe River Basin is 24450 km2with,respectively,5560,5658 and 13232 km2with high,medium,and low coverage,respectively.

Heihe River Basin;grassland;dynamic change;CA-Markov model;LUCC;GIS

10.11686/cyxb2016294 http://cyxb.lzu.edu.cn

赵宇豪,戎战磊,张玉凤,叶苗,蒋晖,赵传燕.近30年黑河流域草地变化分析及分布格局预测.草业学报,2017,26(6):1-15.

ZHAO Yu-Hao,RONG Zhan-Lei,ZHANG Yu-Feng,YE Miao,JIANG Hui,ZHAO Chuan-Yan.Analysis of change in grassland area in the Heihe River basin over the past 30 years and prediction of future trends.Acta Prataculturae Sinica,2017,26(6):1-15.

2016-07-19；改回日期：2016-12-23

干旱气象科学研究基金项目(IAM201608)资助。

赵宇豪(1991-)，男，山东淄博人，在读硕士。E-mail:zyhcn@hotmail.com

*通信作者Corresponding author.E-mail:nanzhr@lzb.ac.cn