植被是生态系统的重要组成部分,在水土保持、减缓温室气体浓度上升等方面有不可代替的作用,因此,植被动态是地理学、生态学研究焦点之一。气候是山地生态系统变化的重要驱动力,植被在全球变化监测中起“指示器”作用[1]。植被是地球气候最鲜明的标志,主要植被类型表现植物对气候类型的反映[2]。气候变化直接影响植被的物候、生产力、分布格局等,区域植被特征的改变反作用于天气系统,二者相辅相成[3]。近年来,全球变暖趋势明显,山地植被对气候变化响应尤为明显[4]。研究山地植被覆盖变化可以为合理开发山区土地资源提供科学依据。

前人研究表明,遥感技术因大范围、长时间序列和经济性等特点,使其成为研究和监测植被动态的有效手段。NDVI（Normalized Difference Vegetation Index）对植被覆盖密度、植被生长状况敏感,并且与植物光合作用效率和植被初级生产密切相关[5]。因此,国内外许多研究人员从不同角度利用NDVI对植被进行研究。

大兴安岭地理位置特殊,自20世纪50年代至80年代中期,中国多次组织对该地区进行综合考察[6]。其植被动态情况不仅决定着当地的生态环境,而且对东北、内蒙古中西部及华北地区生态安全有重要影响。前人研究该区域植被时仅探讨植物区系[7,8]、植被区划[9]等问题,并只做分段讨论,鲜有对整个山脉的综合研究与分段对比分析。郑度等人也主张将大兴安岭作为整体来考察其对蒙古高原、东北平原形成和演变的作用[10]。本研究以NDVI作为植被动态指标,分析大兴安岭1982~2006年植被覆盖动态变化过程,对比不同生态地理区变化程度。分析这一区域植被动态及其与气候因子相互关系和分异规律,对生态环境监测和研究植被对气候变化响应具有重要意义。

大兴安岭位于内蒙古自治区东部、黑龙江省北部,地理位置在40°59′N~53°33′N,115°5′E~125°16′E之间,大致为东北-西南走向。北起黑龙江上游地区,与西伯利亚高原相连;南抵西拉木伦河上游山地,全长约1 500 km。东西坡向不对称,东坡山势陡峭,迅速向松辽平原下切;西侧起伏平缓,向蒙古高原逐渐过渡。大兴安岭为中国重要的气候分界线,地跨寒温、中温带两带和湿润、半湿润、半干旱3个干湿气候区。夏季海洋季风受阻于山地东坡,东坡降水多,西坡干旱,形成明显对比。此外,大兴安岭还是中国地势第二、三级阶梯北段分界线。受水热、地形等自然条件影响,植被类型丰富,不同地区植被生境差异明显,自北向南呈地带性规律分布。北部高纬地区广泛分布寒温带针叶林,随纬度降低向南逐渐过渡到中温带草原景观。

生态地理区域系统主要体现温度、水分、植被覆盖等自然要素的地域差异,可作为生态与环境特征监测的宏观框架。从生态地理分区角度研究植被覆盖变化对气候要素的响应能够更好地揭示大兴安岭区域差异特征。本文依据中国生态地理区域图[10],提取研究区范围,自北向南将其划分为4部分,如图1和表1。

表1   大兴安岭各部分生态地理特征

Table 1   Each eco-geographical region natural conditions on Greater Khingan Mountains

新窗口打开

显示原图| 下载原图ZIP| 生成PPT

图1   研究区示意图

Fig.1   Location of the study area

本文中NDVI数据来源于“中国西部环境与生态科学数据中心”网站（http://westdc.westgis.ac.cn/）,下载GIMMS（Global Inventory Modelling and Mapping Studies）NDVI数据集。该数据集与其它数据相比具有较小的误差和较高的精度,目前已经应用到了相关研究中[11,12]。时间分辨率是15 d,空间分辨率为8 km。考虑到时间序列的完整性,冬季积雪对植被NDVI的影响,并参考前人对东北地区植被物候研究成果[13],选择1982~2006年第9~20期数据进行分析。对数据进行投影变换、掩膜,得到研究区1982~2006年NDVI时间序列。气象数据来自“中国气象科学数据共享网”（http://cdc.cma.gov.cn/home.do）逐日的气温和降水数据,对少量缺值进行补插。

研究表明最大值合成法MVC（Maximum Value Composites）对植被动态监测具有良好效果[14]。利用该方法将15 d的NDVI数据进行合成,获得1982~2006年生长季逐月NDVI时间序列,以确保数据受云、大气、太阳高度角等因素影响最小。在此基础上对每年生长季数据取平均值,得到各年生长季平均NDVI。

2.2.1 趋势分析方法

利用一元线性回归方程斜率分析每个像元25 a间NDVI的变化趋势。国内外许多研究中都应用了该方法[15~17],计算公式如下：

slope=n×∑i=1ni×NDVIi-∑i=1ni×∑i=1nNDVIin×∑i=1ni2-∑i=1ni2(1)

式（1）中,n为总年数,即n=25,i为n a中的第i年,其中,若slope>0时,说明NDVI是增加的,反之为减少。

2.2.2 统计分析方法

相关分析可以用来研究各地理要素之间的关系。将各气象站点月NDVI值与同期的气温和降水数据进行线性相关分析和偏相关分析。相关系数大小及正负反映了因变量与自变量之间的关系;偏相关系数则反映了多要素中,某一自变量在其它要素不变的情况其与因变量的相关关系。

2.2.3 空间自相关分析方法

空间自相关是指地物分布于不同空间位置,某一属性值之间的统计相关性,通常空间距离越近的地物相关性越大。空间自相关分为全局自相关和局部自相关两种测度指标,全局指标用于探测整个研究区域的空间模式,局部指标计算每个空间单元与相邻单元同一属性相关程度。通常采用Moran’s I进行表征,计算公式如下[18,19]：

I=n∑i=1n∑j=1nwij(yi-y̅)(yj-y̅)∑i=1n∑j=1nwij∑i=1nyi-y̅2(2)

Ii=yi-y̅S2∑jnwij(yj-y̅)(3)

式（2）、（3）分别为全局自相关和局部自相关计算公式,I为全局自相关系数,Ii为局部空间自相关系数,n为样本个数,yi和yj分别表示空间对象在i和j位置上的属性值,w为空间权重矩阵,S2为yi的离散方差。

3.1.1 NDVI年内格局变化

大兴安岭自北向南地跨13个纬度,其中温度是大兴安岭NDVI变化的重要因素之一。图2是大兴安岭月积温和与月NDVI均值随纬度变化的平均趋势。可以看出气温和NDVI都有随纬度变化的分异趋势。同一时间气温自北向南呈现升高趋势,月最低积温出现在1月,最冷月南北相差近500℃;月最高积温出现在7~8月,南北相差约290℃。NDVI年内格局与气温都有明显的季节变化格局。在2月份NDVI值最低,随温度升高,NDVI值在4月开始逐渐增加,4~7月是植被增长阶段;NDVI最大值出现在6~8月份;之后,植被进入衰败期,至10月份左右植被活动逐渐停止,进入一个相对稳定的时期。

显示原图| 下载原图ZIP| 生成PPT

图2   大兴安岭月积温（℃）（a）与月NDVI均值（b）随纬度的变化

Fig.2   The monthly accumulated temperature (a) and mean NDVI (b) varies by latitude

冬季,大兴安岭大部分地区植被NDVI值在0.2以下,仅在南段小部分地区在稍高于0.2。由图2b发现,在冬季若以北纬47°为界线,北侧NDVI值低于以南地区,夏季则相反,北侧高于南侧。冬季大兴安岭植被活动受积温影响,南部气温高于北部,植被活动稍强;在夏季,由于大兴安岭南北积温差异小,NDVI由植被类型决定,最北部是中国主要林区,中部分布针阔混交林,因此在研究区内冬季与夏季植被对气温的响应截然相反。

3.1.2 NDVI空间分布格局

为了阐明大兴安岭在1982~2006年间植被覆盖情况,统计了1982~2006年间植被生长季NDVI均值作为研究时间段内植被平均状态,如图3所示。受到天气系统、植被类型、人为干扰等因素的影响,各生态地理区差异较大。大兴安岭植被覆盖较好的地方是北段寒温带森林区和中段山地草原森林区,区内NDVI平均值分别为0.65和0.63。北部以针叶林为主要植被类型,森林茂盛;中段海拔较高的区域分布着森林植被,在低海拔区域为羊草(Leymus chinensis)和贝加尔针茅(Stipa baicalensis)草原。北段西侧是大兴安岭山地向呼伦贝尔草原过渡区域,这里曾经植被茂密,由于海拔相对较低,人类活动对自然植被干扰强烈,自20世纪30年代以后,随着牧业发展不平衡,草原存在过度放牧现象,草原退化非常明显,其NDVI平均值仅为0.57。大兴安岭南段草原区代表半干旱地区,其NDVI平均值为0.51,大兴安岭较小NDVI值集中出现在该区域中部的东坡。

显示原图| 下载原图ZIP| 生成PPT

图3   1982~2006年生长季平均NDVI空间分布

Fig.3   Spatial distribution of average NDVI of growth season of Greater Khingan Mountains in 1982-2006

3.1.3 大兴安岭NDVI变化趋势

线性趋势能够排除短期气候对植被覆盖的影响,并能反映环境变化对植被动态的综合影响[20]。图4是大兴安岭1982~2006年植被覆盖变化线性趋势空间分布。大兴安岭植被整体变化幅度较大,斜率的变化范围在-0.002 87~0.004 32之间,平均值是0.000 74,总体方差是0.000 89。整体上来看,研究时段内大兴安岭地表植被覆盖整体改善的区域比退化的区域面积大,NDVI值呈现增加趋势和减少趋势的像元分别约占总数的80%和20%。

显示原图| 下载原图ZIP| 生成PPT

图4   1982~2006年生长季平均NDVI变化趋势

Fig.4   The variation rate of average NDVI of growth season of Greater Khingan Mountains in 1982-2006

进一步分析不同生态地理区NDVI变化趋势,在表2中统计各生态地理区变化趋势的栅格比例。其中,北段森林区植被覆盖下降趋势最明显,有26.79%的面积呈现下降趋势。另外,中段草原森林区也达到了21.57%。以上两区呈下降趋势的栅格集中在两区交界处的根河市、大兴安岭地区以及鄂伦春自治旗3个地区的交界处。有研究表明,大兴安岭北部和中部人为干扰对景观影响较大,主要表现为农业开垦、道路修建、森林砍伐、居民点建设和采矿[21]。北段西侧作为研究区内生态较为敏感的区域,近年来随着“退耕还林还草”、“三北防护林工程”、“围栏封育”等一些措施的实施,生态保护工作已经初显成效[22],有84.54%面积的NDVI呈增加趋势。南部草原区植被覆盖有明显转好趋势,呈增加趋势栅格数占该区总数的85.81%,其中最东南区域的增幅较大。近年来在这里先后建立了塞罕坝、大海陀等国家级保护区,使植被得到了很好的保护与恢复。

表2   大兴安岭各生态地理区NDVI变化趋势

Table 2   Variation of NDVI in different eco-geographical region on Greater Khingan Mountains

新窗口打开

3.2.1 典型气象站点NDVI变化

在各生态地理区中选取典型气象站点（尽可能在各区中心位置）,分别是图里河、额尔古纳右旗、博克图和林西,分析各典型站点1982~2006年间植被生长季NDVI变化情况。如图5所示,各个气象站点NDVI值均具有波动性,图里河、额尔古纳右旗、博克图3个气象站点NDVI年际变化波动幅度较大,生长季平均NDVI最高与最低值之差在0.1以上,图里河与博克图站点附近NDVI最高。林西气象站点相对平稳,在0.354~0.430之间波动,NDVI值最低。林西气象站点的NDVI值有明显的波动上升趋势,博克图、额尔古纳右旗气象站点在波动中稍有升高,图里河有轻微下降。各气象站点的NDVI趋势与其所在生态地理区植被动态情况基本一致。

显示原图| 下载原图ZIP| 生成PPT

图5   典型气象站点生长季平均NDVI变化

Fig.5   NDVI changing trend of growth season in typical region meteorology stations

3.2.2 NDVI与气象因子的相关性

为了揭示大兴安岭NDVI与气象因子之间关系,将气象站点植被生长季各月NDVI值与相应时间月平均气温、月累计降水量降水进行线性和偏相关分析,如表3所示。大兴安岭NDVI与气温和降水都有显著的正相关性,全部通过PⅢ（0.721）>Ⅰ（0.501）。整体性相对较好的大兴安岭,NDVI空间自相关性自北向南逐渐升高,NDVI分布随机性强,植被分布相对破碎。

3.3.2 NDVI局部空间自相关及动态变化

为了更直观的观察大兴安岭NDVI的空间分布情况,计算多年平均（1982~2006年）、前5 a（1982~1986年）和后5 a（2002~2006年）平均NDVI的局部空间自相关系数,分别代表研究区的平均状况、初始状态和最后的状态,生成LISA聚集图（图6）。该结果反映出大兴安岭NDVI的空间依赖性和异质性,同时表明植被动态的格局和演变过程。从空间分布来看,大兴安岭植被形成4个主要聚集区,NDVI值高-高关联分布在大兴安岭北段、中段和山脉的最南端;低-低关联分布在北段西侧和南段;高-低、低-高关联在研究区内极少的零星分布,体现了植被覆盖的局部离群区域。与之前的分析结果具有较高的一致性。

显示原图| 下载原图ZIP| 生成PPT

图6   大兴安岭NDVI局部空间自相关格局及动态变化

Fig.6   Local spatial autocorrelation pattern and change of NDVI in Greater Khingan Mountains

纵向比较25 a间大兴安岭NDVI局部自相关变化特征,Moran’s I系数有一定程度的下降。位于高-高象限的像元数目减少约10%,位于低-低象限的像元数略有增加。高-高像元减少区域分布在研究区最北端的中部;升高区域在山脉最南端,植被覆盖转好。在北段西侧区域低-低自相关类型有扩张的趋势。说明研究区内NDVI差异在减小,各生态地理区的植被状况趋于相似。

本文基于GIMMS NDVI数据,采用趋势分析法、统计分析、空间自相关分析研究了大兴安岭自1982~2006年的NDVI变化情况,对大兴安岭整体及各生态地理区得出如下结论：

1） 从年内格局、空间分布格局以及动态趋势3个方面对大兴安岭NDVI时空格局进行分析。表明研究区冬季NDVI自北向南升高;在夏季自北向南降低。从空间分布来看,该区北段和中段NDVI值较高,植被覆盖较好。整体上植被有明显的增加趋势,北段和中段退化趋势的区域较大。

2） 大兴安岭NDVI与气温和降水相关关系明显,都通过了0.01的显著性检验。NDVI与气温的相关性更强,在偏相关分析中表现突出,因此气温是影响大兴安岭植被的主要气候因子。通过对不同生态地理区比较发现,由于植被类型的差异,南段草原区与气温的相关性最低。

3） 大兴安岭NDVI整体上具有较高的聚集性,Moran’s I系数达到0.889,相似的NDVI值形成了4个主要的聚集区域。1982~2006年间高-高聚集的区域呈现动态减少,减少的区域大兴安岭北段,增加区域分布在山脉的最南端,低-低自相关变化不大。

The authors have declared that no competing interests exist.

中国植被覆盖变化及其与气候因子关系——基于NOAA时间序列数据

[本文引用: 1]     

植被—气候关系遥感分析研究进展

[本文引用: 1]     

1982~2003年内蒙古植被带和植被覆盖度的时空变化

[本文引用: 1]     

Evidence of response of vegetation to environmental change on high-elevation sites in the Swiss Alps

[本文引用: 1]     

植被归一化指数(NDVI)及气候因子相关起伏型时间序列变化分析

[本文引用: 1]     

[本文引用: 1]     

大兴安岭中段东西两侧植被的分异性

[本文引用: 1]     

大兴安岭南部山地植物区系多样性研究

[本文引用: 1]     

东北植被区划及其分布格局

[本文引用: 1]     

中国生态地理区域系统研究

[本文引用: 2]     

Analysis of Sahelian vegetation dynamics using NOAA-AVHRR NDVI data from 1981-2003

[本文引用: 1]     

中国陆地生态系统近30年NPP时空变化研究

[本文引用: 1]     

东北地区植被物候期遥感模拟与变化规律

[本文引用: 1]     

Multiresolution spatial characterization of land degradation phenomena in southern Italy from 1985 to 1999 using NOAA-AVHRR NDVI data

[本文引用: 1]     

Remote sensing of vegetation and land-cover change in Arctic Tundra Ecosystems

[本文引用: 1]     

1998~2008年中国东部沿海植被覆盖变化特征

长白山区植被生长季NDVI时空变化及其对气候因子敏感性

[本文引用: 1]     

Notes on continuous stochastic phenomena

[本文引用: 1]     

Local indicators of spatial association—LISA

[本文引用: 1]     

Spatio-temporal changes of NDVI and its relation with climatic variables in the source regions of the Yangtze and Yellow rivers

[本文引用: 1]     

大兴安岭北部林区景观格局变化及其影响分析

[本文引用: 1]     

呼伦贝尔森林-草原交错区景观格局时空动态

[本文引用: 1]     

1982~2003年中国东北地区不同类型植被NDVI与气候因子的关系研究

[本文引用: 1]