在全球变化中，植被扮演着陆地生态系统中最重要的角色。植被不仅参与全球物质与能量循环、调节全球气候，同时在自然地理环境中起到“指示器”的作用^[1]。植被覆盖的变化导致了地表反照率、温室气体等地表与大气特性的变化^[2]，通过影响地气系统的水分和能量收支从而影响区域气候发生变化；气候变化则会改变植被群落的组成、结构、影响植被的生物量，进而影响到植被覆盖度的状况^[3]。目前国内外越来越关注植被对全球变化的响应与反馈机制，研究植被覆盖与气候的关系已成为全球变化研究的重要内容。归一化植被指数(NDVI)是遥感影像的近红外波段与红外波段反射率的比值参数，它准确地反映植被的覆盖程度、生长状况以及光合作用的强度^[4-5]，能在较大时空尺度上客观反映研究区植被覆盖信息^[5-6]。因此，归一化植被指数是目前研究植被变化最为常用的指标，成为当前研究生态变化的主要数据集。

近年来我国针对植被覆盖变化及其对气候响应的问题逐渐深入，不同区域进行研究得到的结论也有所差异。杨元合、朴世龙^[7]利用NOAA/AVHRRNDVI数据研究近20年来的青藏高原草地植被覆盖变化及其与气候因子的关系，结果表明高寒草地夏季NDVI的增加是夏季温度和春季降水共同作用的结果，温性草原夏季NDVI变化与气候因子并没有表现出显著的相关关系，高寒草地植被生长对气候变化的响应存在滞后效应；杜灵通、田庆久^[8]利用 SPOT-VGT NDVI数据计算宁夏近11年来的植被覆盖度并对其变化特征进行了研究，并分析其与气候因子之间的关系，结果显示植被覆盖度与气温和降水呈正相关关系；刘斌^[9]等根据GIMMS-NDVI数据对华北地区植被覆盖的时空变化特征进行分析，认为华北地区植被覆盖变化与干燥度指数和气温有很好的相关性，但人类活动的作用远大于气候变化且两者相对作用的大小在不同空间位置上表现不同；余振^[10]等探究了中国东部南北样带不同时间尺度的气候因子和植被活动变化特征，阐述12种具有代表性的植被类型对气候因子的响应方式；廖建华^[11]等以GIMMS-NDVI为基础，结合气候和人类活动数据，分析1982—2003年间长江上游植被覆盖季节变化的空间分布，发现春季、夏季降雨与气温同步增加致使植被覆盖增加，而秋冬季反之；崔林丽^[12]等探讨了近10年中国东部NDVI在全年、各个季节对气温和降水变化的旬时空响应特征，研究发现中国东部植被总体上对气温变化的响应大于降水，在空间上植被对气温变化的最大响应总体表现为北部和中部大于南部，对降水变化的最大响应表现为北部大于中部和南部。

湖北省处于中国地势第二阶梯向第三阶梯的过渡地带，位于长江中游地段，地形复杂多样，生物资源丰富，对该地区植被覆盖变化与气候关系进行单元详细研究对于林区保护、植被恢复、城市可持续发展有着重要意义。本文试图利用2001—2012年MODIS遥感数据、20个气象站点的气温和降水数据，采用线性趋势分析、相关系数等方法，分析研究区归一化植被指数的月变化和年际变化；从植被覆盖度逐年变化特征、空间分布特征、时空统计特征3个方面，多角度分析湖北省近12年来植被覆盖变化特征，并根据植被覆盖度变化趋势分析植被对气候因子的响应程度。

1 研究区域地理特征

湖北省东连安徽，南邻江西、湖南，西接重庆、陕西，河南为邻，介于29°05′—30°20′N,108°30′—116°20′E之间^[13]。地势大致为东、西、北三面环山，中间低平，略呈向南敞开的不完整盆地。全省面积185090 km²，山地约占55%,丘陵24.5%,平原仅占20%^[14]。该地区属亚热带季风性湿润气候，光能充足，热量丰富。无霜期长^[15]，年平均气温15—17℃，降水充沛，雨热同季，年降水量地域分布呈由南向北递减趋势，降水量分布有明显的季节变化，一般是夏季最多，冬季最少，年均降水量在800—1600 mm之间。植被类型主要以亚热带常绿阔叶林、亚热带常绿阔叶落叶阔叶混交林为主^[16]。

2 数据准备与研究方法 2.1 数据来源与预处理

遥感数据采用2001年1月至2012年12月Terra-MODIS13Q1 16d合成的产品数据，空间分辨率250 m，数据来源美国国家航空航天局(https://ladsweb.nascom.nasa.gov/)。数字高程数据(DEM)分辨率为30 m(图 1)，数据来源地理空间数据云平台(http://www.gscloud.cn/)。使用ENVI软件将下载的数据批处理进行格式和投影转换，同时完成图像的拼接和裁切。

气象数据为2001—2012年湖北省地区20个气象站点(图 1)的气温和降水数据，该数据来源于中国气象科学数据共享服务平台(http://www.escience.gov.cn/)。利用ArcGIS 10.0软件对研究区20个气象站点的年均气温和年降水量数据分别进行Kriging插值，得到湖北省地区2001—2012年间年均气温和年降水量。

2.2 研究方法 2.2.1 NDVI均值法

采用最大值合成法将各旬数据合成为月数据，并对12个月NDVI数据求平均值。得到的年均NDVI值可以较好的反映该年份地表植被覆盖状况，而且消除了不同地区季节变化对地表植被覆盖变化的影响^[17]。

2.2.2 NDVI趋势分析

为了定量研究NDVI的变化趋势，本文采用一元线性回归分析的方法，该方法可以模拟每个栅格的变化趋势，以单个像元时间变化特征反映整个空间变化规律，综合反映区域时空格局演变^[18]。对时间自变量和NDVI因变量数据，采用最小二乘法，计算数据集中所有像元的NDVI与时间的回归斜率^[19]。

Slope是像元NDVI回归方程的斜率，若Slope＞0，表示随时间变化植被指数是增加的趋势，并且数值越大植被覆盖度增加越明显；反之，若Slope＜0则代表植被指数呈现下降的趋势。

式中,i为年序号；n为年跨度，即12年；NDVI_i为第i年最大化NDVI值。

2.2.3 相关分析

采用最小二乘法拟合NDVI和气候要素随时间的变化速率，即NDVI和气候要素y随时间t的线性回归系数^[20]。

计算两个变量相关系数公式^[21]如下：

式中,n为统计年数(12a)；i为年序号；NDVI为对应像元12a平均NDVI值；x为气象要素(气温或降水)。

2.2.4 统计分析

为探究NDVI与气候因子之间的相互关系，利用统计分析软件SPSS采用Pearson相关系数和偏相关系数对NDVI与各气候因子进行相应的分析。

3 结果分析 3.1 湖北省NDVI时间变化 3.1.1 NDVI月变化

根据2001—2012年月平均MODIS/NDVI资料统计出NDVI平均值，图 2为湖北省近12年NDVI月平均值，明显看出，随月份变化NDVI值也随之发生改变，8月NDVI值最高为0.80，1月NDVI值最低为0.43。2—8月全省NDVI呈明显上升趋势，主要是因为2月份开始鄂中部平原地带种植棉花、水稻秧苗处于生长的阶段、鄂北岗地的冬小麦和山区的油菜正处于生长时期^[22-23]、鄂东、鄂西山地的树木也抽枝发芽；至8月份，研究区内阳光充足、雨水丰沛植被处于生长旺盛阶段^[24]，各类植被均处于生长旺盛阶段NDVI达到一年中最高值；9月份开始农作物逐渐成熟，大部分植被落叶，植被覆盖减少，NDVI下降，至次年1月降到全年最低值。

3.1.2 NDVI年际变化

图 3是NDVI年均值曲线及其线性趋势，2001—2012年湖北省地区植被覆盖整体呈增加趋势，增速为1%每10年，变化趋势不显著。2001年NDVI值为0.588，达到最低点，这主要是由于三峡工程的实施导致长江上、中游地区植被遭到破坏，生态移民对植被造成了一定的影响；而后NDVI持续上升，2004年NDVI值上升到极点0.619，2005—2006年经历小幅波动后直线上升且在2007年NDVI值达到最高值0.633；2007—2012年间NDVI值呈现先下降减小后上升增加的趋势，其中2008—2010年间有明显的下滑趋势，主要是该时期全国天气异常，遭遇极端天气，南方遭受冰冻、寒潮、冰冻等灾害性天气导致植被受灾，NDVI不断下降^[25-27]。研究区这种波动性的变化与所实施的城市建设、水利林业工程有关，后期一系列的退耕还林、造林绿化的工程实施取得了明显的成效。

3.2 湖北省NDVI空间变化

图 4为湖北省2001—2012年平均NDVI空间分布图。由图可知，湖北省NDVI空间分布呈现自东向西逐渐递增的趋势。NDVI高值区主要分布在鄂西北、鄂西南部地区，其中最大值高达0.821；NDVI低值区分布于鄂中部及东部大范围区域。出现这种分布差异现象的主要原因是由于鄂中部、东部地势较为平坦，河网水系密布，利于居民居住和城市发展，主要城市如武汉市、黄冈市、黄石市、鄂州市、孝感市都聚集分布于此，故中、东部地区植被覆盖度比较低。而西部地势较高，从西南到西北分布着大巴山、武当山，植被类型丰富多样，人类活动相对而言较少所以该地区NDVI值高。

图 5为湖北省2001—2012年年均NDVI变化趋势图，图中发现，2001—2012年湖北省地区植被覆盖变化有明显的区域差异。其中植被覆盖度减小的区域主要集中在武汉、黄石、鄂州、黄冈、孝感、荆州等城市；植被覆盖度增加的地区广泛分布在恩施西南部、十堰西北部、随州东北部等地区。造成这种明显分化的主要原因是自2007年以来“武汉城市圈”的兴起发展，城市的建设用地侵占了耕地、林地资源，导致植被覆盖度严重降低。

根据已有数据的分布规律，采用等间隔等分法将变化斜率分为7个等级(表 1，图 6)。结果显示，严重退化的区域主要是武汉市及其周边地区，退化面积为10.67676km²，占总面积的0.0057%；中度退化的区域零散分布于东部地区，退化的面积为348.8389km²，占总面积的0.19%；轻度退化的区域范围较广，几乎分布于整个湖北省地区，但集中连片的区域主要分布于中东部地区，以武汉市、黄冈市、赤壁市、仙桃市为主，退化为面积11308.63km²，占总面积的6.08%；基本不变的区域的面积为172670.4km²，占总面积92.88%；轻度改善面积1426.35km²，占总面积0.77%；中度改善面积126.5681km²，占总面积的0.6%；植被覆盖明显改善的区域仅占全省陆域面积的0.0046%。NDVI有所改善的区域零星分布于鄂西北角、鄂北部顶端以及鄂东北部小部分地区。据上述分析可以得出，湖北省在近12年中植被覆盖状况基本保持不变，这主要是因为该省东部、北部、南部、植被以常绿阔叶林为主，且西部为神农架原始生态区，神农架地区是保持完好的亚热带森林生态系统，物种多样。

3.3 气候因子变化特征

图 7为湖北省2001—2012年均气温、年均降水量变化及其线性趋势。图中显示，近12年湖北省年均气温变化幅度小，多年年均气温约为16.64℃；2001—2003年呈下降趋势，2004年反弹回峰值，2005年又下降至谷底，2006、2007年维持在峰值阶段，2008—2012年持续下降。从整体上来看年均气温呈下降趋势，大致呈W型，下降的速率大致为0.63℃/10a，通过0.05的显著水平检验。在全球气候变暖的大背景下，区域气候呈现变冷趋势的原因需要进一步探讨。

由图可知，近12年湖北省年均降水量波动较大，2001、2003、2006、2009、2011这5年的降水量相对于其他年份的降水量而言偏低，分别为736、749、873、1035、921mm；而2002年降水量最高，为1285mm；2004—2005年、2007—2008年降水量相对比较平稳，基本维持在同一水平；2010年降水量则达到一个小高峰为1238mm。从整体而言湖北省近12年的年均降水量是呈现上升的趋势，平均增长幅度为11.95mm/a，12年来降水量变化并不显著，未通过0.05的显著水平检验。

3.4 植被覆盖与气候因子的响应 3.4.1 年变化尺度上NDVI与气候因子的相关性分析

图 8为湖北省NDVI与气温和降水的相关关系图，分析图 8可知，NDVI与降水正负相关并存，整体上呈正相关关系。湖北省NDVI与降水量呈正相关的区域面积占湖北省总面积83%，但置信度超过0.05水平的区域不多，仅占总面积的9.47%。这些区域零星分布于湖北的东北部、鄂西北、鄂西南以及鄂东部地区。

图 8显示了NDVI与气温和降水的空间相关性系数，从图看出，气温和NDVI之间是正负相关并存，呈负相关的区域面积较正相关的区域稍大。湖北省内NDVI与气温的空间相关没有达到0.05显著性水平。由此可见湖北省NDVI受温度影响小。

相关系数的平方即为拟合优度，可以相互进行比较而且能直观判断拟合优劣。因此可以利用相关系数平方之差判断NDVI变化主要受气温还是降水量变化的影响^[28]。NDVI与降水量的相关系数平方减其与气温相关系数平方的空间分布结果(图 8)。其中值为正表明NDVI与降水量的相关性大于其与气温的相关性，值越大表示降水量的影响力越强；反之，当值为负时，则表明气温的相关性大于降水量的相关性，负值的绝对值越大表明气温对NDVI的控制越强。分析图可知，NDVI受气温控制的区域面积占研究区总面积的0.46%远远小于降水量对NDVI控制的区域面积。由此可见，气温对湖北省NDVI变化影响小，在少数区域降水对NDVI变化有显著影响。

3.4.2 生长季各月NDVI年际变化对气候因子的响应及其滞后性

滞后性是指植被不能对影响其生长的因子的动态变化做出及时的反应。气温、降水不仅可以影响植被当前生长状况，同时还存在一定的滞后性，可用滞后期表示气象因子对植被变化的影响^[29]。为了研究湖北省植被NDVI对气象因子的滞后响应，选用2001—2012年湖北省地区植被生长季(4—9月)的NDVI与同期气候因子做偏相关分析，同时以2个月滑动平均(上月与本月)和3个月滑动平均(前2个月与本月)与该月NDVI作偏相关分析。

表 2为湖北省近12年来生长季各月NDVI与气温和降水的相关关系。由表 2可知，4月份降水丰富，此时NDVI与降水的相关性好，与当月降水偏相关系数为0.543，而与前2个月和前3个月的降水偏相关系数分别高达0.807、0.633，有着显著相关(P＜0.05)，与前2个月的气温呈显著性负相关(P＜0.05)，与前3个月的气温呈中度负相关，这表明前期的水分累积对植被生长有较大的影响。5月份植被生长达到较为稳定的状态，受气候因子的影响不大。6月份气温升高、植被蒸腾作用变强，植被NDVI与当月的气候因子均无相关性，与前2个月的气候因子有显著的相关性(P＜0.05)，与前3个月的降水有中度相关性。7月是研究区的盛夏季节，温度不是植被生长的限制因子，植被NDVI与当月的降水偏相关系数为0.7，与前3个月的降水偏相关系数为0.73，均达到显著性水平(P＜0.05)。8、9月份植被NDVI与降水的相关性较好。研究表明，植被NDVI变化与气候因子的相关性有着明显的季节性，同时植被NDVI对气候因子的响应存在滞后性，不同月份对气候因子的滞后时间存在差异。

4 结论

本文利用RS和GIS技术，采用多种方法揭示了湖北省地区植被覆盖变化状况以及对气候因子的响应程度。该研究区雨热同期，并且水热条件相对较好，因此气温和降水都是植被生长的主要驱动因素。气温上升和降水量的增加在总体上对多数植被生长有积极作用，但是对于略小部分地区的植被却带来了不良影响导致植被覆盖度与气温、降水呈负相关状态。

湖北省植被覆盖整体较好，2001—2012年间该省植被覆盖整体上呈现增加的趋势，其增长速度为1%/10a，期间有小幅度的下降趋势但下降区域范围不大；从植被覆盖月变化来看，冬季植被呈增加趋势，主要由于冬季该地区农业生产活动加强所致；空间上，植被覆盖格局呈现“西高，东低”的分布特征；根据逐像元的变化斜率分析发现，湖北省植被覆盖在近12年里呈现出总体保持稳定，局部轻微改善的趋势。植被覆盖稳定区分布广泛，其面积占全省面积86.42%；轻度和明显改善的面积共占全省面积10%，主要分布在湖北省西部以及西北部、东北部地区；退化区主要集中在湖北省主城区，长江沿岸地带。

由年际尺度分析可知，湖北省12年年均气温、年均降水与年均NDVI的相关系数较低，但NDVI与降水量的关系较气温更为密切，降水量的适量增加有利于植被的生长；月变化分析表明植被NDVI变化与气候因子的相关性有着明显的季节性也有滞后性。