川西高原植被状况及其与气候变化关系遥感研究.docVIP

川西高原植被状况及其与气候变化关系遥感研究 　　摘要：利用近年来川西高原遥感影像和气象数据，分别以年和季为时间分辨率分析川西高原的植被长势及其与气候变化关系。研究结果显示：近年来川西高原植被长势呈显著增加趋势；除夏季外，各季植被长势增长显著；川西高原夏季，反应植被长势的NDVI值与季均温、季平均积温、季平均降雨量表现出明显的相关性；春、冬季NDVI值与温度、积温的相关性优于降水；春、夏两季平均活动积温和平均降雨量对川西高原NDVI值的综合影响比季均温和季平均降雨量的综合影响要强，冬季则相反。 　　关键词：川西高原；植被遥感；归一化差值植被指数；气候变化 　　中图分类号：Q948 　　文献标识码：A文章编号2016 　　1引言 　　川西高原位于四川省西部，属青藏高原的东南缘（97°E～104°E、27°N～34°N），由阿坝藏族羌族自治州和甘孜藏族自治州2个行政州组成。其地处长江中上游生态屏障区，是全国五大牧区之一的川西北牧区的重要组成部分，也是我国天然林的主要分布区域。该区生态环境特殊多样，植被资源丰富，其84.14%的土地面积分布着高寒草甸、高寒灌丛草甸、乔木林地和灌木林地[1]。 　　植被是维持川西高原生态环境和畜牧业生产的重要资源，其生长发育、空间分布都与气候因子关系密切。在气候变化的背景下研究川西高原植被变化情况以及其与气候变化关系有助于指导高原畜物业生产和植被资源保护。在以往的研究中专门针对川西高原开展的相关研究较少，虽然已有学者对川西高原的植被覆盖与气候因子的关系进行过研究，但研究以年为时间单位且研究的时间范围距今较久[2]。本文使用1999年1月～2013年2月的相关遥感资料和气象数据，主要以季节为时间尺度，开展相关研究。 　　2植被状况的表征 　　由Rouse等提出的归一化差值植被指数NDVI值，是目前已有的40多种植被指数中应用最广的一种，该指数对绿色植被反应敏感且可比性强，常用来进行区域和全球的植被状态研究[3～7]。可反应川西高原的植被长势状况。NDVI的取值范围在-1.0～1.0之间，大于0.1为有植被覆盖，值越大表示植被长势越好。其计算公式为： 　　NDVI=A2-A1A2+A1 　　式中，A1为红光波段反射率， 　　A2为近红外波段反射率。 　　3遥感资料来源及处理 　　使用SPOT卫星VGT传感器NDVI10日合成遥感数据，数据空间分辨率为1 km，下载地址为http：//www.vito-eodata.be。由于SPOT影像NDVI数据被转换成0？250的NDV值，在ENVI平台上采用公式NDVI=NDV×0.004-0.1进行波段运算，将数据还原，再与川西高原边界做掩膜处理，得到川西高原NDVI旬数据。为保证陆面信息的完整性，采用最大值合成法计算NDVI月最大值。NDVI月最大值等权平均后得到季均、年均NDVI值。为保证研究时间段的完整性，使用1999年1月～2012年12月的数据进行年际平均分析（表1），选用1999年3月～2013年2月的数据进行纵横向的季平均分析。 　　4川西高原植被长势变化情况 　　4.1以年为时间尺度的植被长势变化情况 　　在ENVI平台上统计1999～2012年川西高原区域年均NDVI值，并对其做年际变化趋势分析，发现川西高原区域NDVI均值随年份呈增长趋势，增长趋势系数为0.0027，达到 0.01显著性水平。研究时间段内，1999年植被长势最差，NDVI均值仅为0.399，最高值出现在2011年，NDVI均值为0.442。对研究时段内川西高原NDVI年均值图像进行单像元分析，发现研究区域内以单像元为单位的最大NDVI值呈现出更加明显的年际增长趋势，相关系数为0.841。单像元最大NDVI值出现在2011年，为0.801，见图1。 　　ENVI平台上计算植被长势最好年（2011年）与最差年（1999年）的NDVI差值，运算结果如图2所示。经过统计发现，植被覆盖状况最好年与最差年之间的差值在-0.287～0.230之间。2011年川西高原仅零星区域植被状况差于1999年。经统计，94.6%的区域植被长势优于1999年，其中，NDVI差值在0.03～0.06之间的区域最多，约占整个高原面积的43.1%；其次是差值 　　5川西高原植被状况与气候变化的关系 　　5.1气候变化分析 　　选择与植被生长密切联系的温度、积温、降水3个气候因子，分析气候变化情况。统计1999年春季～2012年冬季川西高原31个气象站点各季实测日均温度、季活动积温、季降雨量，再进行等权平均得到各因子区域均值。由于川西高原海拔较高，植被以冷季禾本科牧草为主，计算活动积温的生物学下限温度取零摄氏度。 　　分别对各季气候