马昊翔, 陈长成, 宋英强, 晔 沙, 胡月明,5

(1.华南农业大学 资源环境学院, 广州 510642; 2.华南农业大学地理信息工程研究所, 广州 510642; 3.国土资源部 建设用地再开发重点实验室, 广州 510642;4.广东省土地利用与整治重点实验室, 广州 510642; 5.青海大学 农牧学院, 西宁 810016)

草地是世界上分布最广泛的植被类型之一，大约占据陆地面积20%，是陆地生态系统重要的组成部分[1]。作为生物循环系统中的重要角色，草地为人类提供了重要的生态和社会服务，包括净初级物质生产、碳蓄积与碳汇、调节气候、涵养水源、保持生物多样性等功能[2]。然而，很多地区受全球气候变化和人类活动的干扰，尤其是对气候变化敏感的高寒地区的草地生态系统已经呈现沙化退化趋势，导致草地质量衰退，生物多样性降低，水土流失加剧，生产力和服务功能下降等负面影响，成为备受人类关注的环境生态问题之一[3-5]。近几十年来，全球气候持续变暖，人类活动范围不断扩大以及强度不断加剧[6]，草地生态系统在气候变化和人类活动的共同影响下发生了显著变化，其植被与两者之间的关系研究成为各国学者关注的焦点和核心[7-10]。

植被覆盖是监测草地植被变化的敏感性指标[11]，从一定程度上可以衡量草地生态环境质量的状况。利用植被指数等遥感模型获取区域的植被信息，是目前定量监测植被覆盖变化的有效手段，其中增强型植被指数(EVI)是在归一化植被指数(NDVI)的基础上开发出来的，除了继承归一化植被指数(NDVI)的优点，还改善了其高植被区饱和，大气影响校正不彻底和土壤背景等问题[12]，提高了高生物区植被的敏感度，降低了土壤背景和大气影响，对植被变化的监测具有更高的灵敏性和优越性，在草地退化监测、草地资源定量分析等研究中应用广泛[13]。目前，EVI时间序列在草地植被变化的应用研究有很多[14-16]，但主要集中在草地EVI与气候因素的响应关系，对人类活动影响的关注较少，且以定性描述分析为主，难以定量地区分气候因素和人类活动对草地植被变化的影响，从而限制了对草地植被覆盖变化的正确认识。

位于青藏高原东北部的青海省是我国高寒草地面积分布最广泛的地区之一，也是中国五大牧区之一，生态脆弱性强，是气候变化和人类活动影响的敏感区域[17]，草地沙化退化现象明显。因此，本文以青海省为研究区，提取2005—2014年增强型植被指数时间序列数据集，融合温度和降水因素，运用相关分析和残差趋势分析方法，分析青海省草地植被覆盖时空变化特征，分离量化气候因素和人类活动对草地植被覆盖变化的影响，为青海省草地植被恢复、生态建设和可持续发展提供科学依据。

1 研究区概况与数据来源

1.1 研究区概况

青海省位于中国西北地区，青藏高原东北部，地理坐标为89°35′—103°04′E，31°39′—39°19′N。全省东西长约1 200 km，南北宽约800 km，面积72.23万km2。青海地处亚欧大陆腹地，地势西高东低，平均海拔在3 000 m以上，北部有祁连山，昆仑山系和阿尼玛卿山横贯中部，南部有唐古拉山，构成青海地形地貌的基本骨架，其中分布河谷、盆地、高原和山地，形成祁连山地区、柴达木盆地和青南高原三大自然区域。境内河流众多，是黄河、长江主要水系的发源地，素有“中华水塔”和“江河源头”之称。青海省气候属于典型的大陆性高原气候，特点以高寒干旱为主，多风，日照长，雨量少，平均温度为-5～8.5℃，年均降水量为50～550 mm。

青海省草地资源丰富，是我国五大牧区之一，全省天然草地面积约4 053万hm2，占全省土地面积的56%，居全国第4，其中可利用面积约占全省天然草地面积的87%。草地类型主要分为7类，包括温性草原类、温性荒漠类、高寒草原类、高寒草甸类、高寒荒漠类、低地草甸类、山地草甸类，其中高寒草甸类面积最大、分布最广泛[18]。

1.2 数据来源

(1) 本文中的EVI时间序列是采用了美国国家航天局(National Aeronautics and Space Administration，NASA)提供16 d合成的MODIS(Moderate-resolution Imaging Spectroradiometer)植被指数产品：MOD13Q1(MODIS/Terra Vegetation Indices 16-Day L3 Global 250 m SIN Grid)，空间分辨率为250 m。获取了研究区内(地理范围覆盖全球正弦曲线投影系统中编号为h25v05和h26v05的两个轨道数据)2005—2014年共10 a的EVI产品数据集，并利用MRT将覆盖研究区的两轨数据进行重投影和镶嵌处理，投影转换为Albers，像素大小为250 m×250 m，校正类型是最邻近像元法。原始数据经过了一系列预处理来消除非植被的干扰。为了进一步地消除云和大气带来的噪声，本文进行了低通滤波处理，剔除噪声点，同时利用最大合成算法(maximum value composite,MVC)生成月最大EVI的时间序列数据集。根据已有的研究[7]，由于非生长季植被覆盖度差，EVI数据易受到积雪和土壤反射等影响，因此选择5—9月作为草地的生长季，从5—9月的月最大EVI序列数据集中提取每一个像元的最大EVI值，构建2005—2014年最大EVI的时间序列(记为EVImax)。

图1 研究区位置及气象站点分布

(2) 本文选取2005—2014年青海省区域内37个气象站点的月均温和月降水气象数据。该数据主要来源于中国国家气象局和中国气象数据网。气象数据与EVI植被指数同期，时间间隔单位为月。由于是在像元尺度上进行分析，因此需根据37个气象站点逐月的气象数据，对整个研究区的温度和降水进行空间插值，空间分辨率为250 m。温度的插值采用基于经纬度和DEM回归[19]的方法进行插值，得到生长季(5—9月)平均温度数据。降水的插值采用克里金插值，得到生长季(5—9月)累积降水量数据。

(3) 土地利用数据来源于中国科学院资源环境数据中心，该数据采用土地利用变化遥感信息人机交互快速提取方法，解译分辨率为100 m的青海省的Landsat 7 TM数字影像，得到青海省草地覆盖数据。

(4) 青海省气象分区数据是参考中国的主要气候分区[17]，将青海省划分4个气候区，分别为青南高原区、东部黄土高原区、环青海湖—祁连山区和柴达木盆地区。

2 研究方法

2.1 草地植被覆盖度变化分析

本文首先利用2005年和2014年土地利用数据中的草地覆盖数据作为掩膜提取青海省草地EVI像元。由于青海省植被数据存在不确定性，为了进一步消除非植被因素、裸地、稀疏植被和常绿林等影响[7]，在土地利用数据的基础上，本文通过以下4个条件：(1) 2005年和2014年土地利用数据中土地利用类型均为草地的像元；(2) 生长季(5—9月)EVI均值应大于0.1；(3) 生长季(5—9月)EVI最大值应大于0.15；(4) EVI月最大值应出现在7—9月，进一步地筛选得到青海省生长季草地EVI像元。

为了分析草地EVI逐年的空间变化规律和时空格局演变，本文利用青海省2005—2014年草地EVI数据集，采用回归分析方法，在像元的基础上以最小二乘法拟合草地EVI和年份的线性回归系数，得到不同空间位置上植被绿度变化率，来反映草地EVI的变化方向和速率，其公式[20]为：

(1)

式中：GRC是植被绿度变化率，其反映了研究区草地植被覆盖度的变化速率，n是监测的时间序列长度；n=10；i是年序号，取值范围为1～10；MEVIi是第i年草地年最大EVI值。GRC为正，则说明草地EVI在10年间呈增加趋势，反之则说明草地EVI在10年间呈降低趋势；GRC的绝对值越大，说明草地EVI的变化幅度越大，反之说明草地EVI的变化幅度越小。

2.2 相关性分析

地理系统是一个多要素的复杂系统，一个要素的变化必然引起其他要素的变化，进行地理要素之间的相关分析就是为了揭示地理系统中各要素之间相互关系的密切程度。在多要素的地理系统中，当单独研究两个要素之间的相互关系的密切程度时，可利用偏相关系数来测定，显著性检验采用t检验；当研究几个要素同时与某一要素的相关关系的密切程度时，可利用复相关系数来测定，显著性检验采用F检验[21]。草地EVI与平均温度或累积降水量相关系数计算公式如下：

(2)

草地EVI与平均温度或累积降水量偏相关系数计算公式如下：

(3)

(4)

T检验计算公式为：

(5)

式中：ret,rer,rtr分别为草地EVI与温度、草地EVI与降水、温度与降水的单相关系数；ret·r为草地EVI与温度的偏相关系数；rer·t为草地EVI与降水的偏相关系数；r为偏相关系数；n为样本数；m为自变量个数。

草地EVI与平均温度和累积降水量复相关系数计算公式如下：

(6)

F检验计算公式为：

(7)

式中：re·tr为草地EVI与平均温度和累积降水量复相关系数；ret为草地EVI与温度的单相关系数；rer·t为草地EVI与降水的偏相关系数；r为复相关系数；n为样本数；m为自变量个数。

2.3 基于EVI的残差趋势分析

残差趋势法是基于生态系统中净初级生产力和气候因素的相关关系进行分析的。在长时间序列的NDVI数据中，假设气候变化是唯一的驱动因素，采用基于像元的相关分析方法，检测NDVI的残差趋势，从而剔除气候因素的影响来分析人类活动导致的植被覆盖变化[22]。本文利用监测植被变化能力更优的EVI数据代替NDVI数据来进行残差趋势分析，主要分为以下两个步骤：

(1) 基于局部空间区域构建长时间序列EVI数据和气候因素的回归分析模型。本文采用平均气温和累积降水量两个气候因素进行回归分析，即对每个像元建立EVI与气温和累积降水量的响应关系，来模拟仅受气候变化影响下每个像元的年EVI。但由于草地植被覆盖数据存在空间非平稳性，而且难以从历史数据中获取受人类活动影响较小的样点作为基准，本文假设在一定局部空间范围内，植被、气候、地形等自然条件是相似的，植被覆盖度与气温和降水具有相同的响应关系，该空间内与气温和降水相关系数较高且生物量最大的像元可以看作是受人类活动影响较小的样点。因此，本文以每一个像元为中心设立5 km×5 km邻域空间，即该邻域空间内植被、气候、地形等自然条件是相似的，搜索该邻域内相关系数大于0.7且EVImax最大的像元作为训练样点，构建2005—2014年气候因素和年EVImax的回归分析模型，分析气温和降水对EVI的综合影响，模拟该区域仅受气候变化影响下每个像元的年EVI，其公式如下：

(8)

(9)

(10)

将每个像元的残差项按年序排列，以线性回归拟合其年际变化趋势，进行人类活动对草地植被覆盖变化影响的趋势分析，其公式如下：

ε=at+c

(11)

式中：ε为每个像元的EVI残差项；a为拟合趋势的斜率，如果a>0，则表明像元在去除气候变化的影响下，人类活动对草地植被覆盖变化有正面的影响，如果a<0，则表明像元在去除气候变化的影响下，人类活动对草地植被覆盖变化有负面的影响，如果a=0，则表明像元在去除气候变化的影响下，人类活动对草地植被覆盖变化的影响较小。

3 结果与分析

3.1 青海省生长季草地EVI时空变化特征

3.1.1 青海省生长季草地EVI时间变化特征 为了研究青海省生长季草地植被覆盖随时间变化的特征，取2005—2014年EVI数据中草地EVI平均值，代表当年的植被覆盖状况，得到10年间青海省不同区域生长季草地EVI时序变化图(图2)。从图2中可以看出，青海省草地EVI值为0.35～0.45，呈波动增加趋势，年增长速率为0.002 3，增长趋势不明显(R2=0.3181，p>0.05)。

从不同气候分区来看，生长季草地EVI差异明显，其中东部黄土高原区草地EVI平均值明显高于其他地区，生长季草地EVI平均值为0.441 2，说明东部黄土高原区草地植被覆盖比其他地区要好；青南高原区和环湖及祁连山区生长季草地EVI平均值比较接近，分别为0.409 3，0.381 4；柴达木盆地生长季草地EVI平均值为0.246 5，是所有地区中最低的，说明总体上柴达木盆地是青海省草地植被覆盖最低的地区。2005—2014年青南高原区、东部黄土高原区和柴达木盆地生长季草地EVI均呈现增长趋势，环湖和祁连山区草地EVI呈现退化趋势。东部黄土高原区年变化速率为0.005 3，达到了0.01的显著性水平，呈现显著增加趋势；青南高原区、柴达木盆地年变化速率均为0.002 2，增加趋势不明显；环湖和祁连山区年变化速率为-0.001 4，退化趋势不明显。

图22005-2014年青海省生长季草地EVI变化趋势

3.1.2 青海省生长季草地EVI空间变化特征 尽管在年际变化上青海省生长季草地EVI的变化趋势不明显，但却存在明显的空间差异。如图3所示，青海省生长季草地EVI的空间分布呈现东南部高，西北部低的分布特征。由东南至西北，平均EVI逐渐降低。其中青南高原区南部和东南部、环湖和祁连山区温度较低，具有良好的降水条件，适合高寒草甸类、高寒草原类的植被生长，具有较高的EVI值；东部黄土高原区多河谷，多盆地，气候适宜，但由于城镇密集，人口密度大，对草地植被的破坏较严重，因此EVI值较低；青南高原区北部和西部海拔较高，温度较低，柴达木盆地夏季降水较少，蒸发量大，不利于草地植被的生长，因此，这两个区域的EVI值最低。

利用像元统计的方法对10年间青海省EVI平均值进行分级统计分析，结果表明：草地低植被覆盖区(平均EVI值介于0.15～0.35)占草地总面积的40.38%；草地中植被覆盖区(平均EVI值介于0.35～0.55)占草地总面积的42.61%，其中0.35～0.45的区域占20.93%，0.45～0.55的区域占21.68%；草地高植被覆盖区(平均EVI值大于0.55)占17.01%，其中0.55～0.65的区域占13.46%，大于0.65的区域占3.55%。

由于草地EVI平均值只能反映草地植被覆盖度的分布状况，因此，根据公式(1)计算每一个像元10年来的植被绿度变化率，进一步地分析青海省生长季草地EVI的空间变化。根据草地植被变化趋势的不同并参考自然断裂法[23]，将研究区内草地植被变化情况划分为5种变化类型：明显改善(Ⅰ级，GRC≥0.005)；轻微改善(Ⅱ级，0.002

图32005-2014年青海省生长季草地EVI均值空间分布

图4 2005-2014年青海省生长季草地EVI绿度变化程度表1 2005-2014年青海省草地EVI均值变化类型

3.2 青海省生长季草地EVI对气候变化的响应关系

3.2.1 生长季草地EVI与气候因子的偏相关性分析 温度和降水是青海省草地生长的主要限制因子，对草地植被覆盖的空间分布起着决定性的作用。通过计算2005—2014年青海省生长季草地EVI与同期的平均温度和累积降水量的复相关系数可知，其相关系数为0.738，通过p<0.05的显著性检验，表明青海省生长季草地EVI与温度和降水密切相关。

为了进一步分析温度和降水对生长季草地EVI影响的空间分布，本文采用偏相关分析方法逐像元的计算平均温度和累积降水量与草地EVI的偏相关系数，得到结果如图5所示，可以看出青海省生长季草地EVI与平均温度和累积降水量的关系是正相关和负相关并存。利用ArcGIS 10.2软件对区域内的像元进行统计，得到生长季草地EVI与平均温度的偏相关系数为-0.881 3～0.891 7，生长季EVI与累积降水量的偏相关系数为-0.882 4～0.883 0，空间差异明显。同时计算得到生长季草地EVI与平均温度、累积降水量的平均偏相关系数分别为0.271 0，0.137 3，表明总体上生长季草地EVI与平均温度和累积降水量都呈正相关关系。

研究区内生长季草地EVI与平均温度呈正相关关系的占全区的72.02%，呈负相关关系的占27.98%；与温度呈正相关关系的主要分布在青南高原区，环湖和祁连山区，其中31.34%的像元通过显著性检验(p<0.05)，集中分布在青南高原西部、扎陵湖和鄂陵湖环湖区域、阿尼玛卿山东西侧、茶卡—共和盆地西部和祁连山区北部，说明这些地区的草地EVI随着温度的升高而增大。这些地区属于青海省典型的高寒区和山地区，海拔较高，平均海拔在4 500 m以上，温度较低，平均温度在-2℃以下，降水较为充沛，主要的草地植被类型为高寒草甸类和高寒草原类，对温度的敏感性较高。当温度升高，有利于冻土层的融化，促进草地植被的光合作用，有利于草地的生长。因此，在这些地区，热量是草地生长的主要限制因子，温度与生长季草地EVI呈明显的正相关关系。而与温度呈负相关关系的主要分布在柴达木盆地和东部黄土高原区，其中7.15%的像元通过显著性检验(p<0.05)，集中分布在柴达木盆地东部、茶卡—盐湖环湖区域、沙珠玉河流域和河湟谷地，说明这些地区的草地生长随着温度的增加而退化。这是由于这些地区海拔较低，平均温度较高，降水相对稀少，主要的草地类型为温性草原类，当温度升高会使地表蒸发量进一步增加，加剧该地区水分的缺乏，从而抑制草地的生长。

研究区内生长季草地EVI与累积降水量呈正相关关系的占全区的60.38%，呈负相关关系的占39.62%。与降水呈正相关关系的主要分布在青南高原区、柴达木盆地和青海湖环湖区域，其中24.67%的像元通过显著性检验(p<0.05)，集中分布在青南高原西部和北部，柴达木盆地东北部、青海湖西侧和沙珠玉河流域，说明这些地区降水与草地EVI呈明显的正相关关系，降水是草地生长主要限制性因素，原因是这些地区降水量相对高山区和北部高寒区较少，平均温度较高，蒸发量比较大，对降水的敏感性较高。与降水呈负相关关系的主要分布在环湖和祁连山区、青南高原区东北部，其中11.55%的像元通过显著性检验(p<0.05)，集中分布在阿尼玛卿山区和祁连山区，说明这些地区降水与草地EVI呈明显的负相关关系，降水越多反而导致草地的退化。这可能是因为这些地区属于降水较为充沛的高寒山地区，降水量增多，云量增加，会导致太阳辐射减少，草地植被的光合作用减弱，导致草地生长受到制约。

图5 2005-2014年青海省生长季草地EVI与温度、降水偏相关系数

3.2.2 不同变化类型生长季草地EVI与气候因子的驱动关系 由于青海省生长季草地EVI与气候因素的相关性存在明显的空间差异，因此本文进一步地探讨不同草地植被变化类型下草地EVI与气候因素的驱动关系，结果如表2所示。

柴达木盆地研究结果表明，基本不变区域(Ⅲ级)和轻微退化区域(Ⅳ级)生长季草地EVI与降水的偏相关和与温度、降水的复相关显著，与温度的偏相关不显著，主要受温度和降水共同驱动影响，降水的影响比较大；轻微改善区域(Ⅱ级)生长季草地EVI与温度、降水的复相关显著，主要受温度和降水共同驱动影响；明显退化区域(Ⅴ级)和明显改善区域(Ⅰ级)与温度、降水的偏相关和复相关均不显著，说明气候因素并不是主要的驱动因素。

东部黄土高原区研究结果表明，基本不变区域(Ⅲ级)和轻微退化区域(Ⅳ级)生长季草地EVI与降水偏相关显著，与温度的偏相关和与温度、降水的复相关不明显，主要受降水驱动；明显改善(Ⅰ级)、明显退化(Ⅴ级)和轻微退化区域(Ⅳ级)生长季草地EVI与温度、降水的偏相关和复相关均不明显，说明这些区域主要受非气候因子驱动。

环湖和祁连山区研究结果表明，明显改善(Ⅰ级)和轻微退化区域(Ⅳ级)生长季草地EVI与温度、降水的复相关显著，主要受温度、降水驱动，其中明显改善区域(Ⅰ级)和轻微退化区域(Ⅳ级)与降水的偏相关明显，受降水的影响较大；基本不变区域(Ⅲ级)生长季草地EVI与温度偏相关显著，与降水偏相关和与温度、降水的复相关不明显，主要受温度驱动；轻微改善区域(Ⅱ级)和明显退化区域(Ⅴ级)生长季草地EVI与温度、降水的偏相关和复相关均不明显，说明该区域主要受非气候因子驱动。

青南高原区研究结果表明，基本不变区域(Ⅲ级)和轻微改善区域(Ⅱ级)生长季草地EVI与温度、降水的复相关显著，说明这些区域主要受温度、降水驱动；轻微退化区域(Ⅳ级)生长季草地EVI与温度偏相关显著，与降水和与温度、降水的复相关不明显，说明该区域主要受气温驱动；明显改善区域(Ⅰ级)和明显退化区域(Ⅴ级)生长季草地EVI与温度、降水的偏相关和复相关均不明显，说明这些区域主要受非气候因子驱动。

值得注意的是，明显改善(Ⅰ级)和明显退化区域(Ⅴ级)生长季草地EVI均值基本上与温度、降水的相关性均不显著，可能的解释是这些区域主要是人类活动比较大的城市区和畜牧区，除了气候因素对草地生长的影响之外，人类活动对草地生态系统的扰动也比较剧烈。因此为了进一步探讨人类活动对草地植被生长的影响，利用残差趋势法分析去除气候变化的人类活动对草地植被生长影响的空间分布。

3.3 青海省生长季草地EVI与人类活动的关系

本文基于局部回归的残差趋势法计算反映人类活动的残差趋势序列，得到2005—2014年青海省生长季草地EVI的残差趋势图(图6)。残差趋势绝对值越大表明该区域草地生长受人类活动的干扰影响越大，其中人类活动对草地生长的影响主要有两种表现形式：一是残差趋势值为正且越大，表明人类活动对草地生长具有促进作用(如生态工程等)，二是残差趋势值为负且越小，表明人类活动对草地生长具有破坏作用(如放牧、城市扩张和采矿等)；残差趋势绝对值越小表明该区域草地生长受人类活动的干扰影响较小。

2005—2014年青海省生长季草地EVI年平均残差趋势为0.005 7,表明人类活动对青海省草地植被生长呈正影响。从图6中可以看出残差趋势绝对值大的区域集中在人口活动强度比较大的城镇区和放牧区。残差趋势呈显著增加的集中分布在青南高原区西南部、扎陵湖—鄂陵湖区域、海南洲和东部黄土高原区的河湟谷地区域，表明该区域草地植被生长不能仅仅用温度和降水的变化解释，很大程度上还受人类活动的影响，且人类活动对草地植被生长起了促进作用。残差趋势呈显著减少的集中分布在杂多县和治多县交界处、囊谦县、称多县、玛多县南部、阿尼玛卿山两侧区域、班码县、久治县和青海湖区域，表明这些区域的草地EVI持续减少，草地正面临退化的风险，且主要是由人类活动引起的。

图7反映了不同变化类型下人类活动对草地EVI变化的影响差异，图中可知随着草地退化程度的降低或草地改善程度的增加，草地EVI残差趋势值呈显著降低或显著增加的趋势。明显退化区域草地EVI残差趋势值基本小于0且差异明显，说明在明显退化区域(Ⅴ级)，人类活动对草地生长起着主要的破坏作用；同样的，明显改善区域草地EVI残差趋势值基本在0以上且差异明显，说明在明显改善区域(Ⅰ级)，人类活动对草地生长起着主要的促进作用。此外通过比较图4和图6可以发现，青海省生长季草地EVI的变化趋势和人类活动的残差趋势在空间分布上十分相似，再通过表2的结果表明，明显退化和明显改善区域的生长季草地EVI与温度、降水的相关性均不明显，说明草地生长变化除了响应气候变化之外，很大程度上是受人类活动的影响，尤其是在草地植被明显变化的区域，人类活动是影响草地植被覆盖变化的主要限制因素。

表2 不同变化类型生长季草地EVI均值与气候因子的偏相关系数及复相关系数

注：*表示通过0.05的置信度检验；**表示通过0.01的置信度检验。

图6 2005-2014年青海省生长季草地EVI残差变化趋势

图7 草地植被不同变化类型的残差变化趋势

4 讨论与结论

(1) 2005—2014年青海省草地植被整体呈波动增加趋势，生长季草地EVI的年增长速率为0.002 3，其中东部黄土高原区整体呈显著增加趋势，青南高原区、柴达木盆地整体呈增加趋势，环湖和祁连山区整体呈退化趋势。青海省草地EVI呈东南部高，西北部低的分布特征，草地植被改善区域集中分布在青南高原区西部、海南洲和东部黄土高原区的河湟谷地、扎陵湖和鄂陵湖环湖区域，祁连山区北部，草地植被退化区域集中分布在阿尼玛卿山两侧、青海湖盆地区域、玛多县和称多县。

(2) 青海省生长季草地EVI与温度和降水密切相关，总体上与平均温度和累积降水量均呈正偏相关关系，31.34%的像元与温度呈显著正相关，集中分布在青南高原西部、扎陵湖和鄂陵湖环湖区域、阿尼玛卿山东西侧、茶卡—共和盆地西部和祁连山区北部；24.67%的像元与降水呈显著正相关，集中分布在青南高原西北部、柴达木盆地东北部、青海湖西侧和沙珠玉河流域。因此，温度和降水在一定程度上是草地植被生长的主要限制因素，且对于高寒地区的青海省来说，温度与EVI的关系更为密切。

(3) 青海省草地植被变化主要是气候变化与人类活动共同作用的结果。柴达木盆地(Ⅱ级、Ⅲ级和Ⅳ级区域)。环湖和祁连山区(Ⅰ级和Ⅳ级区域)、青南高原区(Ⅱ级和Ⅲ级区域)主要受气温和降水驱动；东部黄土高原区(Ⅲ级和Ⅳ级区域)主要受降水驱动；环湖和祁连山区(Ⅲ级区域)、青南高原区(Ⅳ级区域)主要受气温驱动；其他区域主要受非气候因子驱动，人类扰动对这些区域的草地植被覆盖变化影响较大。人类活动对青海省草地植被的生长总体呈正影响，其中起着促进作用的区域集中分布在青南高原区西南部、扎陵湖—鄂陵湖区域、海南洲和东部黄土高原区的河湟谷地区域，起着破坏作用的集中分布在杂多县和治多县交界处、囊谦县、称多县、玛多县南部、阿尼玛卿山两侧区域、班码县、久治县和青海湖区域。在草地植被明显变化的区域，人类活动对草地生态系统的建设和破坏起着主导作用。

青海省草地植被覆盖变化是气候和人类活动等因素共同耦合作用的结果。本文虽量化了人类活动对草地植被变化的影响，但仅把人类活动作为一个影响因子来考虑。而事实上，人类活动也是一个多因素共同影响的整体，在多因素影响下，其对草地植被变化存在建设和破坏两种不同的效应。一方面是建设效应[24]，如在青南高原区西部、扎陵湖—鄂陵湖区域，自2005年以来，国家批准启动三江源自然生态保护和建设，建立了索加—曲麻河保护区和扎陵湖—鄂陵湖保护区，实施了退耕退牧还草、生态移民等多项重大工程，有效地改善草地的植被覆盖状况。另一方面是破坏效应[25-26]，如在称多县、玛多县、阿尼玛卿山区域，海拔比较高，地形较为复杂，不利用生态工程建设的实施，人类的活动方式主要以放牧为主，超载放牧导致草地生态系统持续退化；在青海湖区域，近年来环湖旅游资源和矿产资源得到有效的开发，但同时也给青海湖环湖区域的草地生态环境带来了很大压力。因此，在更小尺度上进行实地调查研究，进一步细化反映人类活动的影响因素，量化各人为因素对青海省草地植被变化的相对影响仍有待深入地分析和探讨。