黄土高原地处半干旱半湿润气候带,生态环境脆弱,土质松散,抗蚀性能差,加之长期乱垦滥伐,植被破坏殆尽,致使该地区水土流失严重,是我国水土保持重点区域^[1-3]。植被是联结土壤、大气和水分之间物质、能量交换的关键环节,是全球环境和气候变化的敏感指示器,具有截留降雨、减少雨滴击溅、增加土壤入渗、保土固土等功能,对减少水土流失有重要作用^[4-5]。为改善黄土高原植被覆盖状况,治理水土流失,自 1999 年开始国家实施了退耕还林(草)、封山育林以及荒山造林等大规模植被建设^{[3, 6-7]}。加强黄土高原植被覆盖时空变化研究,有助于更好地理解陆地生态系统的动态变化特征,可为黄土高原生态建设提供有用的空间信息和理论支持^[8-10]。

归一化植被指数(Normal Difference Vegetation Index,NDVI)是遥感影像近红外波段反射值(Near Infrared Reflection,NIR)和红光波段反射值(R)之差与NIR和R之和的比值,常被用来表征植被覆盖和植被生长状况,是目前科学研究中最常用的植被指数^{[5, 11-13]}。NDVI可以对植被生长动态变化进行监测,同时能够在较大时空尺度上客观反映植被信息^[12]。近年来,国内外众多学者已经对不同地区的植被NDVI时空变化趋势及其与气候变化的关系和响应特征作了大量研究^{[8, 13-15]}。例如,李本纲和陶澍利用AVHRR NDVI数据探讨了全国植被覆盖与温度和降雨的关系^[16]；Wu等利用MODIS NDVI研究了北京—天津沙源区2000—2010年植被覆盖时空变化趋势^[17]。目前对于黄土高原地区植被覆盖时空演变分析已有一些阶段性成果^{[5, 8-10]},但对整个黄土高原在大规模植被建设实施后,植被覆盖年、生长季、季节及月变化和空间分布的研究仍相对较薄弱。

黄土高原是中国气候变化的敏感地带,研究发现,20世纪80年代以来,黄土高原地区气候变化趋势明显^{[8, 18-20]}。同期,在黄土高原开展了大规模的植树造林等生态建设活动^{[2-3, 7]}。本文对该背景下的黄土高原地区植被覆盖在不同时间尺度上的变化进行分析,并在此基础上探讨气候变化对植被变化的影响。以期研究结果为该地区环境保护、生态建设和水土流失的综合治理提供理论依据。

1 数据与方法 1.1 研究区概况

黄土高原地处中国黄河中上游地区(33°43′—41°16′N,100°54′—114°33′E),平均海拔1000—1500 m,总面积约为6.4×10⁵ km²。该区属温带大陆性季风气候,四季分明,年平均温度为3.6—14.3℃。受温带大陆性气候的影响,年降水量在150—750 mm之间,全年降水主要集中在7—9月,占全年降水量的60%—80%。

1.2 数据来源及处理 1.2.1 NDVI数据

本研究所采用的NDVI数据来源于比利时佛莱芒技术研究所(Flemish Institute for Technological Research,VITO)发布的10日最大化合成SPOT-VEGETATION(SPOT-VGT)S10 NDVI产品,空间分辨率为1km,时间尺度为1998年4月至2014年5月。该数据集已经过严格处理,通过辐射校正、几何校正、大气校正等处理之后,可保证NDVI数据的质量。本文使用的数据为1999年1月至2013年12月覆盖东亚的SPOT-VGT S10 NDVI数据集(可在http://free.vgt.vito.be免费下载)。

下载后的数据采用VITO提供的VGT Extract软件来提取,设置输出范围和输出格式(ENVI,unsigned Integer)后,通过批处理的方式提取SPOT-VGT S10 NDVI数据集。在ArcGIS 10.0中,利用研究区边界截取黄土高原地区NDVI时序数据集,根据公式NDVI=DN×0.004-0.1^[21-22]将旬DN值转化为对应的旬标准NDVI真实值,并将数据投影系统转为Albers投影,即可得研究区的逐旬NDVI时序数据集。为进一步降低噪声并提高数据质量,本文对研究区逐旬NDVI时序数据集,通过IDL7.0提供的Savitzky-Goaly滤波方法^[12]重建高质量NDVI时间序列数据。

1.2.2 气象数据

1998年1月—2013年12月年黄土高原地区51个气象站点的日降水和温度数据来源于中国气象科学数据共享服务网(http://cdc.cma.gov.cn),本研究中选取的气象站点空间分布见图 1。最后,利用ArcGIS10.0软件,对51个气象站点的温度和降雨数据分别进行克里格(Kriging)插值,获得整个黄土高原地区的平均温度和平均降雨量。

1.3 分析方法

根据研究需要,将研究时段划分为年际、生长季、季节、月份四种时间尺度。季节划分采用气象学上的标准：春季是3—5月,夏季是6—8月,秋季是9—11月,冬季是12月和次年1—2月,生长季是4—9月^[21]。采用最大合成法(Maximum Value Composites,MVC)获得各研究时段内最大NDVI 值作为该研究时段的 NDVI 值。

NDVI随时间的变化趋势利用slope分析,即对每个栅格的NDVI数值进行一元线性回归分析,计算公式为：

式中,slope为变化速率,即斜率；i表示年序号,n为研究序列的长度,NDVI_i为第i年的NDVI值。若slope＞0,则该像元15年间的植被覆盖度变化趋势是增加的,且数值愈大植被覆盖度增加趋势愈明显；反之则表明植被变化趋势是减少的；slope=0表示此像元的植被覆盖度无变化。

此外,本研究中的偏相关分析与线性拟合分析分别借助SPSS17.0软件和Origin8.0 软件完成。

2 结果与分析 2.1 黄土高原地区NDVI变化特征 2.1.1 黄土高原地区NDVI时空分布特征

黄土高原地区月NDVI的最小值出现在1月份和2月份(均为0.16),随着降雨的增加以及气温的升高,NDVI也逐渐增大。其中NDVI值较高的月份是7—9月,8月份研究区NDVI值达到最高(0.50),之后NDVI开始下降。图 2为研究区不同季节,生长季和年最大NDVI平均值(1999—2013年)的空间分布图。从不同季节1999—2013年平均NDVI可以看出,冬季黄土高原植被覆盖度最低(平均NDVI为0.19),绝大部份地区NDVI处于0.1—0.2之间,占黄土高原面积的59.06%。随着温度回升,黄土高原春季植被覆盖逐步增加,研究区NDVI平均值为0.35。在夏季,黄土高原地区降水量增加,植被茂密,NDVI明显增大并达到最大值(平均NDVI为0.51),大部分地区NDVI主要分布在0.6—0.7之间,占整个研究区的17.34%。秋季黄土高原大多数植被开始落叶、凋零,农作物开始收割,地表植被减少,NDVI与夏季相比也有减小趋势,整个研究区平均NDVI为0.46。生长季NDVI和全年NDVI的空间分布范围较一致,研究区NDVI值主要分布范围为0.5—0.6,分别占整个研究区的17.81%和17.80%,另外整个研究区平均NDVI分别是0.52和0.31。

从图 2黄土高原地区NDVI的空间分布图中可以看出,在各研究时段内,黄土高原东南部地区NDVI值较高,其次是东部,而西北部干旱半干旱地区(荒漠、沙地以及裸岩)NDVI值较低。黄土高原地区东南部植被覆盖明显高于西北部的原因,可能是东南部水热条件较好,植被类型主要为针阔叶林、寒温性灌丛和草丛^[9]；而西北部由于东南季风带来的暖湿气流不能深入内陆,主要的植被类型是草原、草甸和高山稀疏植被,导致NDVI较低^[9]。

2.1.2 黄土高原地区NDVI变化趋势

图 3为黄土高原地区NDVI从1999年到2013年的变化趋势。从图中可以看出,不同季节、生长季和全年的NDVI整体上呈增加的趋势。除冬季(P=0.020)外,其他季节NDVI与生长季和全年NDVI均呈极显著(P＜0.01)的增加趋势。据统计,研究区四季NDVI平均增加速度由快到慢依次为夏季(0.0087)＞秋季(0.0081)＞春季(0.0064)＞冬季(0.0018)。可见,黄土高原地区植被覆盖增加主要发生在夏秋两季,研究结果与张宝庆等^[5]的研究结果相同。另外,生长季NDVI平均增长速度(0.0087)与夏季NDVI增长速度无差异。年NDVI的变化速度最大,为0.0099。此外,夏季和生长季NDVI与全年NDVI的变化较为一致,1999—2001年黄土高原地区NDVI处于相对较低的水平,自2001年开始,NDVI迅速上升,在2005年有短暂的下降之后,进一步增长。整体上来看,从1999年开始黄土高原地区NDVI呈增加趋势,说明退耕还林(还草)政策的实施,使得黄土高原地区的生态环境逐年改善,植被状况总体呈好转的趋势。

为进一步研究黄土高原地区植被覆盖变化的幅度和空间分布,对各像元NDVI与时间进行一元线性回归分析。结果显示(图 4),黄土高原地区1999—2013年年最大NDVI变化斜率最大的地区是黄河中游区域,尤其是水土流失最为严重的丘陵沟壑区,其植被恢复态势最为明显。而中西部部分区域NDVI的变化斜率较小。从6个不同研究时段来看,黄土高原地区NDVI都表现出不同程度的增加趋势,大部分区域NDVI变化斜率主要在0—0.1之间,分别占整个黄土高原的68.13%(春季),60.32%(夏季),62.87%(秋季),65.87%(冬季),61.04%(生长季)和51.79%(全年)。此外,本研究将NDVI变化斜率小于零的区域定义为植被退化区。1999—2013年黄土高原植被退化面积为19.1×10³km²,仅占研究区总面积的2.99%。

2.2 黄土高原地区NDVI对水热条件的响应 2.2.1 黄土高原地区气候的时间变化特征

从表 1月平降雨量和温度可以看出,两者的变化趋势基本一致。降雨量和温度的最小值均出现在12月份到次年2月份,其中降雨量较集中和温度较高的月份均在6—9月,在7月份达到峰值。由月NDVI的变化趋势和表 1可以看出,月NDVI、月降雨量和月温度的增速相似,且存在NDVI对于降雨和温度的滞后效应。

图 5为不同季节、生长季的降雨量和温度变化趋势。从图中可以看出,冬季降雨量和温度最低,而夏季降雨量和温度最大,其次为生长季。从降雨量的变化趋势可以看出,除春季降雨量有减小趋势外,其他季节降雨量均呈增加趋势,但变化趋势均不显著(P＞0.05)。冬季降雨量的年际波动较小,而夏季、秋季和生长季的降雨量年际波动较大。从温度的变化趋势可以看出,不同研究时段的温度,均有降低的趋势,但不显著(P＞0.05)。

对黄土高原地区1998—2013年51个气象观测站点的降雨量和温度实测资料分析发现,黄土高原气候变化呈现“冷湿化”趋势。降雨量呈增加的趋势,变化斜率为3.023；温度呈现降低的趋势,变化斜率为-0.046(图 6)。但是,由于降雨量和温度年际波动较大,研究区降雨量增加(P=0.293)温度降低(P=0.061)的趋势不显著。

2.2.2 研究区NDVI对水热条件的响应

从NDVI与降雨量、温度的偏相关分析结果来看,在春季、夏季、秋季、生长季和全年,NDVI与降雨量均呈正相关,而冬季降水与NDVI呈负相关(表 2)。在植被变化对温度变化的敏感性方面,四季、生长季和全年的平均NDVI和平均温度均呈负相关。原因可能是：(1)温度是植被覆盖增加的限制因子。温度上升可以加速地表蒸散发过程,由此引起土壤干化从而抑制植被生长,降低NDVI；(2)与其他因素相比,人类活动对NDVI的影响占绝对主导地位,如最近十几年研究区退耕还林草措施实施和农业生产水平的提高促使植被NDVI显著增加,而最近十几年的年气温呈现不显著的微弱下降趋势(图 8),即气温对NDVI的负效应是一种伪效应^[21]。具体原因,还需要进一步深入定量的分析研究。然而,研究区NDVI与降雨量和温度的相关性均不具有统计的显著性(P＞0.05),说明研究区植被的季节和年际变化分别与季节和年降雨量和气温的关系并不密切。由于降雨和温度变化具有时间尺度上的差异,仅仅通过降雨和温度与NDVI在季节和年际变化尺度上的关系,并不能充分说明植被变化对水热条件的响应^[12]。并且,考虑到水热条件对植被生长影响的滞后性,有必要从月时间尺度上分析植被变化对水热条件的响应。

为了进一步研究NDVI对降雨和温度的响应,本研究分别计算了1999—2013年180个月份NDVI与降雨及温度的偏相关系数。气候因子分别选取当月、上月和上上月3个特征时间的值与当月NDVI计算偏相关系数,结果见表 3。从表中可以看出,研究区当月NDVI与当月、上月及上上月的降雨和温度存在显著的相关关系(P＜0.05)。NDVI与当月降雨的相关性高于与当月温度的相关性。但是,NDVI与上月及上上月温度的相关性要强于与降雨的相关性。另外,NDVI与降雨和温度的最大的偏相关系数分别出现在当月(R_{NDVI 降水·气温}=0.606)和上月(R_{NDVI 气温·降水}=0.783)。可见,植被生长对温度的响应存在一个月的滞后期,而降雨无滞后效应。

3 结论

(1) 退耕还林(草)工程实施以后,黄土高原地区植被覆盖有所改善,1999—2013年研究区年最大NDVI的平均值为0.31,研究区大部分地区年NDVI的平均值在0.5—0.6之间。NDVI在空间分布的差异较明显,其中黄土高原东南部水热条件有利于植被生长,NDVI较高,而西北部植被覆盖度较低。

(2) 自1999年以来,黄土高原地区年最大NDVI呈极显著(P＜0.01)的增加趋势。四季与生长季NDVI也呈上升趋势,四季NDVI增加速度由快到慢依次为夏季(0.0087)＞秋季(0.0081)＞春季(0.0064)＞冬季(0.0018)。生长季NDVI的增长速度与夏季一样,为0.0087。

(3) 1998—2013年,黄土高原地区气候呈现降雨增加温度降低的“冷湿化”趋势。但由于降雨和温度的年际波动较大,降雨和温度的变化趋势不明显。

(4) 在季节、生长季和全年时间尺度上,NDVI与降雨和温度的相关性不显著。在月时间尺度上,NDVI与当月降雨量的相关性要强于与当月温度的相关性。另外,植被生长对温度的响应存在一个月的滞后期,而降雨无滞后效应。