基于SPOT VEGETATION数据的长江流域植被覆盖变化特征分析

摘要

植被覆盖动态变化是土地资源和环境变化研究的重要内容。本文以遥感技术为基础，基于1999-2009年SPOT VEGETATION的时间序列影像数据，综合应用多种数学方法，分析了长江流域植被覆盖的空间分布和动态变化特征。
　　 用平均值法得到的多年平均NDVI影像图显示长江流域植被覆盖空间分布上的平均状况，用slope分析法、Mann-kendall趋势分析法、变化矢量分析法(CVA)和Hurst指数R/S分析法对不同尺度（年、季、月）的NDVI时间序列进行数学运算，分析长江流域植被覆盖的变化特征，以及不同流域分段、水系分区、不同海拔梯度以及不同植被类型NDVI的空间分布和时间变化特征。研究结果显示：
　　 (1)长江流域植被覆盖的整体状况较好，1999-2009年的年平均NDV值＜0.1无植被区仅占流域总面积的1.8％，NDVI在0.5-0.65之间的中高植被覆盖区和NDVI＞0.65的高植被覆盖区分别占到长江流域总面积的47.4%和31.7%。季节水平上，春秋两季主要以中高植被覆盖为主，夏季以高植被覆盖为主，而冬季以中植被覆盖为主；月水平上7月中高、高植被覆盖区的面积比占到了流域总面积的88%，2月份的植被生长状况最差，年均NDVI＞0.5以上的面积比重最小，不到9%。空间分布整体上符合南多别少，东多西少的趋势。年NDVI值较少的区域主要分布在青海长江源头和中下游大型湖泊及周边大中城市。
　　 (2)年内变化上，长江流域NDVI年内变化呈单峰型，2月值最低，7月份达到最大值后开始逐渐下降。4月出现最大增幅，10月有最大减幅。在不同流域段水系区内NDVI的年内变化曲线有所差异，下游流域段的各水系NDVI年内变化呈现双峰型，尤其是下游农田的变化曲线，峰值出现在4月和8月，6月有小低谷值；中游、下游和＜200m的低平原在8月出现最大值。东部地区的生长季比西北部地区结束要晚。
　　 (3)在年际变化中，近11年来长江流域植被覆盖整体呈显著上升趋势，年均增加斜率为0.0059，11年增幅为0.0649，在1999年的基础上增长了14.25%。。其中，夏、秋两季对长江流域植被覆盖状况的改善贡献最大，季均NDVI增幅分别为0.070、0.067，上升了11.92%、13.75%。在不同植被类型中，森林、灌草地及农田NDVI均呈明显增加趋势，其中农田增幅最大为0.079，上升了17.17%。在不同海拔梯度上，500-1000米海拔带上的NDVI增幅最大为0.092，上升17.94%。从水系分区上看，仅太湖水系的植被覆盖呈退化和减少趋势。
　　 (4)变化趋势的空间分布上，近11a来，长江流域年均NDVI有92.3%像元呈上升趋势，而降低像元数仅占7.7%。表明长江流域近十来年的“长防林工程”、“天保工程”等植被恢复工程的生态效益已开始显现。表现减少趋势的区域主要分布在长江源头，金沙江流域中上游，岷沱江水系北部，局部江南丘陵地带，和部分大中城市。根据F检验，变化趋势不显著地区面积比例为39.4%，一般显著增加的区域在流域内呈零散分布，面积占18.5%，大部分地区变化呈显著增加，面积占到42.1%。春、夏、秋、冬四季呈显著增加面积分别占流域面积的54.1%、42.1%、35.5%、9.8%。12个月当中，4月份显著增加的面积最大，占到了43.5%的流域面积，12月显著增加的面积最小，只占4.5%。
　　 (5)采用变化矢量(CVA)分析法对1999-2009年、1999-2004年和2004-2009年3个时间段进行了变化矢量的模(‖Δp‖)计算，发现半数以上的像元6年或11年NDVI‖Δp‖的值位于0.2-0.4的低变化区。通过转移矩阵发现2004-2009年间的变化强度不如1999-2004年。1999-2009年的‖Δp‖空间分异：下游＞中游＞上游；乌江水系＞上游干流区＞太湖水系（前3位），金沙江水系的‖Δp‖值最小；500-1000m的低海拔山地‖Δp‖较大，海拔＞4000m的高海拔区‖Δp‖较小；4种主要植被类型上，农田＞灌丛＞森林＞草地。‖Δp‖的大小与植被造林工程、农林业生产和城市化进程等人类活动息息相关。
　　 (6)Mann-kendall的分析结果与slope分析结果很接近，两者之间对年际NDVI的变化显著与不显著上有较一致的判断，但在区分增加趋势极显著和显著水平上有10%左右的差异。根据Z值的大小分布可以较灵活的选择显著性水平是此方法的一个优势。年均NDVI序列Mann-kendall分析结果的Z值空间分异：在0.1水平上(1.64)下游NDVI增加趋势不显著；水系区内金沙江、岷沱江增加趋势不显著、太湖水系减少趋势不显著；＞4000m高海拔地区增加趋势不显著；植被类型中草地的增加趋势不显著，其余流域区段/海拔/植被类型的变化趋势显著。
　　 (7)Hurst指数分析结果显示，长江流域年均NDVI不是纯粹的随机游走而是有偏的随机游走序列，具有自相似性和时间依赖性。全流域上的Hurst指数平均值为0.8344，各流域区段/海拔梯度/植被类型像元上的平均Hurst指数都在0.75以上，表明变化趋势具有可持续性。其中下游Hurst指数比上游和中游要高；水系区内汉江水系、中游干流区和下游干流区的Hurst指数较高；海拔梯度上的Hurst指数值差别不大，以500-1000m的低海拔山地最高，＞4000的高海拔地区最低；植被类型上农田＞灌丛＞森林＞草地。说明这些地区/植被类型上的NDVI序列数据的自相似性更高，持续性相对较强，可能与人类规律性的生产活动有关。
　　 上述结论将为长江流域生态环境的保护和建设，土地资源管理和决策以及区域经济、环境、社会的可持续发展，提供一定的参考依据。