土地利用变化是全球变化的研究核心, 而植被是土地利用类型变化不可或缺部分^[1]。研究植被覆盖变化是研究全球气候变化对环境影响的直接条件, NDVI(Normalized Difference Vegetable Index)是指示大尺度上植被覆盖和植被生产力的良好指标^[2-3]。随着遥感监测技术的不断成熟, 利用遥感卫星对地球监测的范围也越来越广, 国内外利用遥感对植被覆盖变化的监测已经有相当多的研究^[4]。目前, 获取植被覆盖变化监测数据的卫星有许多：NOAA-AVHRR, SPOT-VGT, EOS-MODIS等, 其中MODIS-NDVI数据是从2000年开始, 时间序列相对完整, 数据分辨率也相对较高的数据集, 在植被覆盖变化监测方面有较好的应用空间^[5]。国内外对于中国NDVI长时间序列植被覆盖变化特征的研究较少, 一般是存在于片区或省市的长时间序列的研究或者是季节性的研究^[6]。对大尺度, 长时间序列对中国土地利用类型植被覆盖变化的研究存在研究的必要性, 是监测我国近年来植被响应气候或国家政策的重要依据, 也有助于相关部门对中国全局的了解和决策^[7-10]。

1 数据处理 1.1 数据来源

研究选自美国地质局(NASA)的MODIS/Terra 16天合成MOD13Q产品, 分辨率为1km, 时间刻度为2000年到2013年。数据预先采用美国NASA提供的MRT软件选取NDVI数据集和进行投影转换, 本文选用WGS84投影坐标系；然后利用Arcgis 10.0和中国矢量边界图进行裁剪, 并进行MVC最大化合成, 得到2000年—2013年年NDVI数据集^[11-17]。中国土地利用数据来源于2000年中国1:10万土地利用类型分类数据集；然后利用Arcgis10.0对土地利用类型数据集进行投影和矢量化处理, 得到土地利用类型矢量数据集。

1.2 MODIS数据处理

基于2000年至2013年年均MVC最大化数据集, 由于值在-10000到10000之间, 所以先把值缩小到-1至1取非负值得到NDVI数据集, 然后利用批量化裁剪得到逐一土地利用类型NDVI地图, 再逐一计算平均NDVI值构成3种土地利用类型(耕地, 林地, 草地)时间序列数据集, 最后对数据进行线性回归分析得到每一种土地利用类型时间序列趋势图, 趋势方程和相关系数^[18-20]。本研究数据处理均采用Arcgis10.0和python2.7处理或编程处理。

1.3 一元线性回归趋势分析

植被指数的逐年内每个栅格变化之间对应时间序列存在一定的相关性^[21-23]。随着时间的推移栅格内NDVI值有减少或增加的趋势, 为了体现这种趋势的变化, 采用一元线性回归分析方程来分析2000—2013年中国NDVI值得空间分布特征。

式中, K为斜率, NDVI_i表示第i年NDVI值, 用该方法可以计算出时间序列下植被覆盖的空间变化趋势, K>0表示植被呈改善趋势, 反之为退化趋势。这里用该方法计算中国3种土地类型的NDVI空间分布趋势图, K≤-0.001为退化的趋势, K≥0.001为改善的趋势, 其他为基本不变的区域。

2 结果与分析 2.1 3种土地利用类型NDVI变化特征分析

图 1可以看出3种土地类型：林地, 耕地和草地14年来年平均NDVI变化趋势图, 图中可以看出3种土地利用类型的植被覆盖指数总的趋势都是明显增加的, 中国近十几年来植被总体呈现出增长的趋势, 而且增长的显著性较好。其中, 林地NDVI的增长趋势为0.002, 相关系数为0.8211, 耕地NDVI的增长趋势为0.0022, 相关系数为0.685；草地NDVI的增长趋势为0.0024, 相关系数为0.6403。然而林地有4年NDVI处于波峰然后有明显下降, 分别是2002年, 2004年, 2007年和2010年；耕地有五年NDVI处于波峰然后有明显下降, 分别是2002年, 2004年, 2008年, 2010年和2012年；草地有三年NDVI处于波峰然后有明显下降, 分别是2003年, 2008年和2012年, 可以看出各种土地类型植被的生长趋势都与各年份的环境因子有关, 才会在同一时间段有增长和下降的趋势。3种土地利用类型中, 均在2009年NDVI值最低, 而在2008年中国大陆全国范围内发生冰灾, 显而易见这种极端气候对3种土地利用类型的影响非常明显。由于国家政策的支持, 如六大林业重点工程、对荒漠地区的开发与利用、以及城市绿化、植树造林等, 使得我国各种土地类型的植被覆盖指数明显呈现出了增长的趋势, 其中林地的增长趋势最为稳定, 说明退耕还林、对天然林地的保护等政策对林地植被覆盖植树的改善非常明显, 而草地的退牧还草及荒漠地草原化等也取得了不小的成就。

2.2 3种土地利用类型再分类分析

耕地分为水田和旱地；林地分为有林地, 灌木林, 疏林地和其他林地；草地分为高覆盖草地, 中覆盖草地和低覆盖草地。从图 2和表 1中可以看出, 水田有四个波峰分别是2002年、2005年、2007年和2010年, 相应波谷在2003年、2006年、2009年和2013年, 总趋势线有增长的趋势, 相关系数为0.3114, 说明水田NDVI变化与时间关系性不大, 可能与人为因素和气候因素影响较多；旱地有三个波峰在2005年、2008年和2012年, 而波谷在2006年和2009年, 相关系数为0.7632, 增长的趋势较明显；有林地的相关系数为0.7457, 增长趋势较显著, 有三个波峰在2002年、2005年和2010年和两个明显波谷在2006年和2009年；灌木林, 疏林地和其他林地的相关系数分别为0.8278和0.8268, 增长趋势显著, 而且灌木林和疏林地峰值分布相同(波峰在2002年、2007年和2010年；波谷在2003年、2006年和2009年), 趋势线相似分别为0.0022、0.0023和0.0022；高覆盖草地有三个波峰在2003年、2008年和2012年, 两个波谷在2007年和2009年, 相关系数为0.3836, 增长趋势显著性较低, 总趋势有增长的趋势；中覆盖草地和低覆盖草地的趋势k值相近分别为0.0028和0.0027, 相关系数分别为0.7115和0.7109, 增长趋势较明显, 两种草地的波峰分布相同在2003年、2008年、2010年和2012年, 两种草地的波谷分布相同在2006年、2009年和2011年。可以看出各种植被类型在2009年均呈现波谷阶段, 除高覆盖草地在2007年出现波谷外均在2006年呈现波谷阶段。由3种土地类型再分类分析可以看出, 近14年来由于国家对环境保护的重视, 退耕还林、退牧还林以及对荒漠地的开发与利用使得各土地类型植被覆盖指数均有增长的趋势, 而林地的增加使得水田面积的减少, 水田的增长趋势波动性最大, 同时高覆盖草地的波动性最大是由于高覆盖草地对环境如降水、温度等的响应非常敏感。

2.3 3种土地利用类型NDVI时间序列空间变化特征分析

图中可以看出耕地, 林地, 草地近14年来的生长趋势均以改善为主。图 3耕地中改善面积占总耕地的64.21%, 退化的区域占18.50%, 改善比较明显的区域位于甘肃以南、宁夏、陕西以北、东北三省以西和新疆少数地区, 而退化比较明显的区域位于江苏、浙江、上海交汇地带的繁华区以及大部分省会城市附近。2000年来, 耕地改善区域主要位于西北地区, 而东部繁华区以及人口聚集的城市中心区域, 耕地植被覆盖指数明显减少, 说明耕地由东北地区向西北地区移动, 东部建设用地、生态用地对耕地的占用以及退耕还林等政策对耕地的影响非常明显；图 4林地改善的区域占总林地的54.21%, 退化的区域占20.13%, 其中改善比较明显的区域主要在甘肃以东、陕西以北、山西、吉林和辽宁以北、南部少部分地区, 退化明显的区域位于黑龙江部分地区、内蒙古以东、湖南东部、四川中部、新疆以北、西藏以东、台湾中南部地区。东北大小兴安岭、长白山, 西南地区以及东南沿海地区、四川盆地及东南丘陵植被覆盖指数的增强与退耕还林有关, 而华北地区、黄土高原等地植被覆盖指数增强与草地林地转换有关, 而林地减少区域主要集中在东北传统林区以及东部沿海的人口聚集区, 这是由于人类对林地的砍伐以及建设用地等对林地的占用造成的, 但总体以改善为主, 说明退耕还林, 建立天然林地保护区等政策对林地的保护起到了至关重要的作用；图 5草地改善的区域占55.53%, 退化的区域占18.23%, 改善明显的区域主要位于甘肃东部、宁夏、陕西以北、山西以西、内蒙古以东部分地区, 而退化的区域主要位于内蒙古东南部地区、新疆北部地区和西藏中部地区。草原区域主要集中在北部, 由于国家退牧还草工程等实施, 草地植被覆盖指数有明显的增强的趋势, 而部分地区由于降水量的减少, 加上其承载力比较脆弱, 以及部分放牧的影响, 草地植被覆盖指数有减少的情况。

3 结论

(1) 通过3种土地利用类型平均趋势图, 分析3种土地类型的趋势变化。3种土地类型植被指数在2000—2013年来均有增长的趋势, 其中草地增长趋势最大, 从显著性来看, 林地增长更为显著, 说明近14年来林地增长速率稳定, 而草地增长最快但波动性也比较明显。六大林业重点工程、对荒漠地区的开发与利用、以及城市绿化、植树造林等, 使得我国各种土地类型的植被覆盖指数明显呈现出了增长的趋势, 其中林地的增长趋势最为稳定, 说明退耕还林、对天然林地的保护等政策对林地植被覆盖植树的改善非常明显, 而草地的退牧还草及荒漠地草原化等也取得了不小的成就。

(2) 通过3种土地类型再分类分析, 分析3种主要土地类型的再分类类型的变化特征。耕地中, 水田NDVI的趋势显著性最低, 增长趋势较低, 旱地NDVI显著且增长速率大；林地中, 4种林地的NDVI趋势显著性较高, 增长稳定, 其中有林地NDVI增长速率最低；草地中, 高覆盖草地NDVI趋势显著性最低, 说明波动性很大, 而增长速率也比较低, 中, 低覆盖草地NDVI显著并增长速率大。近14年来由于国家对环境保护的重视, 退耕还林、退牧还林以及对荒漠地的开发与利用使得各土地类型植被覆盖指数均有增长的趋势, 而林地的增加使得水田面积的减少, 水田的增长趋势波动性最大, 同时高覆盖草地的波动性最大是由于高覆盖草地对环境(如降水、温度等)的响应非常敏感。

(3) 由3种土地利用类型植被覆盖空间图, 分析3种土地利用类型2000—2013年来空间变化特征。耕地改善区域主要位于西北地区, 而东部繁华区以及人口聚集的城市中心区域, 耕地植被覆盖指数明显减少；林地减少区域主要集中在东北传统林区以及东部沿海的人口聚集区, 这是由于人类对林地的砍伐以及建设用地等对林地的占用造成的, 但总体以改善为主, 说明退耕还林, 建立天然林地保护区等政策对林地的保护起到了至关重要的作用；由于国家退牧还草工程等政策的实施以及草原化荒漠建设, 使得草地植被覆盖指数有明显的增强的趋势。3种主要土地类型NDVI均以改善为主, 其中耕地植被覆盖指数改善区域最大占总耕地面积的64.21%；3种土地类型均有所改善且改善明显的区域主要集中在甘肃以南, 陕西以北和东北部分地区。