华北地区植被覆盖变化及其影响因子的相对作用分析_刘斌

第30卷第1期自然资源学报Vol.30No.1

Jan.,20152015年1月JOURNALOFNATURALRESOURCES

华北地区植被覆盖变化及其影响因子的

相对作用分析

刘斌1,2，孙艳玲1*，王中良1，赵天保3

（1.天津师范大学城市与环境科学学院，天津300387；2.南京师范大学地理科学学院虚拟地理环境教育

部重点实验室，南京210023；3.中国科学院东亚区域气候-环境重点实验室，北京100029）

摘要：利用GIMMSNDVI数据和气象数据，采用趋势分析、残差分析和相对作用分析对华北

地区1981—2006年植被覆盖时空变化特征进行了分析，并计算了气候变化和人类活动在植被

覆盖变化过程中的相对作用。结果表明，1981—2006年华北地区植被NDVI呈现显著上升趋

势，其增加速率为0.009/10a，但却存在着明显的空间差异，且植被NDVI退化区域面积大于改

善区域面积；华北地区植被覆盖变化与干燥度指数和气温有很好的相关性，说明气候变化是影

响植被覆盖变化的重要因素；此外，无论在华北地区植被改善区域还是退化区域，人类活动起

到的作用都占据了主导地位。在植被改善区，人类活动的相对作用为68.10%，大于气候变化的

相对作用（31.90%），在植被退化区，人类活动的相对作用为71.88%，也远大于气候变化的相对

作用（28.12%），且气候变化和人类活动的相对作用大小在不同空间位置上表现不同。

关键词：华北地区；植被覆盖；气候变化；人类活动；相对作用

中图分类号：Q948.1文献标志码：A文章编号：1000-3037(2015)01-0012-12

DOI：10.11849/zrzyxb.2015.01.002

植被是连结土壤、大气和水分等自然要素的纽带，能够在很大程度上代表一个地区的生态环境总体状况，并在全球变化研究中起着“指示器”的作用[1-2]。植被指数能够表征地表植被特征，有效地反映地表植被覆盖信息。其中，归一化植被指数（NDVI）是公认的植被生长状况及植被覆盖度的最佳指示因子，是监测地区或全球植被变化的有效指标之一[3-4]，被广泛应用于全球和较大区域植被状况监测、生态环境评估、作物长势检测以及自然灾害监测等方面[5-7]。目前常用的NDVI数据有NOAA-AVHRR、EOS-MODIS、SPOT-VEGETATION（VGT）和LandsatTM等。其中，GIMMSAVHRRNDVI数据从1981年开始获取，时间序列最长，覆盖范围广，常被用于监测大尺度植被的动态变化[8]。随着人类活动的加强，植被覆盖变化深刻地记录了人类活动的烙印，是气候变化和人类活动共同作用的结果[9]。近20a来，国内外学者利用NDVI数据对植被覆盖变化及其与气候变化的关系进行了深入研究，但是对人类活动对植被变化的影响关注较少，且对气候变化和人类活动对植被变化的影响研究主要停留在定性描述上[10-12]，使其影响因素分析的客观性受到影响，造成植被覆盖变化与其影响因素分析的脱节。因此，有必要对气候变化和人类活动在植被变化过程中相对作用的定量研究进行探讨分析。

收稿日期：2013-10-18；修订日期：2014-07-21。

基金项目：全球变化研究国家重大科学研究计划项目（2012CB956203）；国家科技支撑计划项目课题（2012BAC07B02）；教育部新世纪优秀人才支持计划（NCET-10-0954）；国家自然科学基金项目（41001022）。

第一作者简介：刘斌（1991-），男，硕士研究生，主要从事气候模拟与全球变化研究。E-mail:[email protected]*通信作者简介：孙艳玲（1977-），女，博士，副教授，主要从事资源环境遥感与全球变化研究。E-mail:[email protected]

1期刘斌等：华北地区植被覆盖变化及其影响因子的相对作用分析13

华北地区地处我国北方，属于沿海湿润、半湿润季风气候向内陆干旱、半干旱气候的过渡带。受气候暖干化的影响，再加上频繁的人类活动，华北地区生态环境十分脆弱[13-14]。近年来，华北地区植被覆盖变化情况备受社会公众和科学界的广泛关注。国内许多学者对华北地区及其局部区域植被变化开展了大量的研究[15-21]。然而，以往对华北地区植被变化的研究主要集中在对其植被变化特征及其与气候的关系以及人类活动对植被变化的定性探讨上，定量分析气候变化和人类活动对华北地区植被变化的影响研究还未见报道。因此，本文利用1981—2006年共26a的NDVI数据和气候数据，对华北地区植被覆盖时空变化特征进行分析，并分离出其变化过程中气候因素和人为因素的影响，在此基础上定量评价这两大影响因子在华北地区植被覆盖变化过程中的相对作用大小，为华北地区可持续发展和生态环境建设提供科学的理论依据。

1数据与研究方法

1.1研究区概况

华北地区行政上包括北京市、天津

市、河北省、山西省和内蒙古自治区，

如图1所示。该地区面积广阔、地形复

杂，植被类型多样，受温带季风气候影

响显著，夏季高温多雨，冬季寒冷干

燥，春秋季较短，从沿海到内陆降水逐

渐减少，四季分明，光照充足，适宜农

作物及植被生长。此外，华北地区是我

国的政治文化中心，人口众多，自然资

源丰富，经济发展迅速，但该地区干旱

化趋势显著，生态环境较为脆弱且受人

类活动的影响较大。图1华北地区范围和气象站点分布1.2NDVI数据来源与处理Fig.1LocationofNorthChinaandmeteorologicalstationsNDVI数据集来源于美国国家航天航

空局（NationalAeronauticsandSpaceAdministration，NASA）全球观测与制图研究组（GlobalInventoryModelingandMappingStudies，GIMMS）正式发布的基于NOAA气象卫星数据，从1981年7月至2006年12月的每半月合成的数据集，空间分辨率为8km×8km。采用最大值合成法（MaximumValueComposites，MVC）获取年最大NDVI值，用以代表当年植被生长的最好状况，最终得到1981—2006年共26a的华北地区年NDVI时间序列数据。MVC法可以进一步减少来自云、气溶胶、视角以及太阳高度角等带来的影响[22]。其中，1981年的NDVI数据从7月开始，但由于华北地区植被生长状况最好的时期是7、8月，因而缺失的数据不会对当年最大NDVI值产生影响。

1.3气象数据来源与处理

气象数据来源于中国气象科学数据共享服务网提供的1980—2006年间华北地区及其周边的223个气象站点逐月的气温和降水数据（图1）。首先逐站点计算1980—2006年间每年的年均温、年降水量、干燥度指数、大于0℃

月均温之和、生物热量指数和生物干

14自然资源学报30卷湿度指数6种气候指标。然后利用克里格插值获得华北地区这6种气候指标在1980—2006年间共计162幅空间栅格图像，其投影方式和空间分辨率与NDVI时间序列数据一致。其中，干燥度指数和生物热量指数、生物干湿度指数的计算详见参考文献[23-24]。

1.4研究方法

1.4.1趋势分析

趋势分析采用的是一元线性回归趋势分析，该方法可以模拟出每个栅格的变化趋势，进而综合反映整个区域的时空变化规律[25-27]，其计算公式如下：

Slope=nöænöænn×(i×MNDVI,i)-ç∑i÷ç∑MNDVI,i÷

i=1øèi=1øèi=1

nnöæn×∑i2-ç∑i÷

i=1èi=1ø

式中：Slope为像元NDVI回归方程的斜率，变量i为1~26的年序号，n为年跨度，MNDVI,i（1）代表第i年最大NDVI值。若Slope>0，表示随时间变化NDVI值呈上升趋势，且Slope值越大，表示上升趋势越明显；反之，若Slope

1.4.2残差分析

残差分析法是由Evans和Geerken提出[28-29]，目前在研究中得到了广泛的应用[9,30-34]。该方法是对每个栅格像元的NDVI与气候指标做回归分析，从而得到每个像元NDVI的预测值，该预测值可以视为气候因子对NDVI的影响，然后利用遥感观测的NDVI真实值减去NDVI预测值，就得到了人类活动对NDVI的影响，从而剥离开气候变化和人类活动分别对植被覆盖变化的影响。其表达式如下：

ε=NDVI真实值-NDVI预测值（2）

式中：ε为NDVI中人类活动所贡献的部分，即为NDVI残差值，当ε＞0时，表示人类

活动产生正面影响；当ε＜0，表示

人类活动产生负面影响；当ε≈0，

表示人类活动影响微弱。

其中，气候指标选取的是与

NDVI相关性最好的当年的干燥度指

数和前一年的年平均气温，并对研

究区内每个栅格像元逐个建立NDVI

与这两个气候指标的二元一次线性

回归方程。

1.4.3相对作用分析

植被覆盖变化过程中气候变化

和人类活动的相对作用的计算主要

采用许端阳等提出的计算相对作用

大小的思路和方法[35]，如图2所示。

图2植被覆盖变化过程中气候变化和人类活动

相对作用的评价流程

Fig.2Flowchartofassessingtherelativeroleofclimatechangeand

humanactivitiesintheprocessofvegetationcoverchange首先通过残差分析，在每个栅格像元上分别得到了气候变化和人类活动造成的NDVI

值，其中，气候变

1期刘斌等：华北地区植被覆盖变化及其影响因子的相对作用分析15化造成的影响用NDVI预测值来衡量，人类活动造成的影响用NDVI残差值来衡量。然后对NDVI预测值和NDVI残差值分别作趋势分析，再结合NDVI真实值在1981—2006年间的变化趋势，分别建立了3种不同可能情景下计算气候变化和人类活动相对作用大小的评价方法，如表1所示。

表1各种情景下气候变化和人类活动在植被覆盖变化过程中的相对作用评价方法

vegetationcoverchangefordifferentscenarios

Table1Methodsforassessingtherelativeroleofclimatechangeandhumanactivitiesintheprocessof

注：ΔNDVIC为t+i时刻气候变化造成的NDVI预测值与t时刻NDVI预测值的差值；ΔNDVIH为t+i时刻人类活动造成的NDVI残差值与t时刻NDVI残差值的差值。

2结果与分析

2.1华北地区植被覆盖时空变化特征

由图3可以看出，1981—2006年间华北地区植被覆盖状况整体上呈明显上升趋势，NDVI线性倾向率为0.009/10a，通过了

P

来植被覆盖情况总体有所改善。从不同时

间段来看，华北地区植被NDVI在20世纪

80年代期间呈波动性上升趋势，90年代期

间总体表现为先升后降，21世纪以来则显

著上升，植被恢复持续好转。

尽管华北地区植被NDVI整体呈现增

长趋势，但却存在着明显的空间差异。通

过对华北地区逐个栅格像元进行一元线性

趋势分析，并对其结果进行标准差重分图31981—2006年间华北地区植被NDVI变化

Fig.3NDVIchangeinNorthChinafrom1981to2006

16自然资源学报30卷

类，得到华北地区植被NDVI变化趋势

图。由图4和表2可看出，1981—2006年

间华北地区植被覆盖状况得到改善区域占

总面积的25.14%，退化区域占27.77%。

且植被改善区域主要分布在内蒙古赤峰和

通辽南部地区、兴安盟东部、河套平原、

锡林郭勒盟、山西北部以及河北西北部和

东南部地区。植被退化区域主要分布在大

兴安岭北部、呼伦贝尔草原、满洲里西北

部、呼和浩特和包头以北、阿拉善左旗、

阿拉善右旗、山西中部、河北唐山市附近

地区、北京市、天津市等。

2.2植被改善区气候变化和人类活动的相

对作用大小

本文中提到的植被改善区包括表2中from1981to2006

的植被明显改善区域、中度改善区域和表21981—2006年华北地区植被变化趋势结果统计

Table2StatisticalresultoftrendofvegetationchangeinNorth轻微改善区域。根据植被改善区域气候

变化和人类活动相对作用的3种可能情Chinafrom1981to

2006

景，计算得到1981—2006年间气候变化

和人类活动在华北地区植被改善区域的相对作用。从整体上来看，在植被改善区

域，气候变化的相对作用为31.90%，人

类活动的相对作用为68.10%，可见26a

间人类活动是促使华北地区植被改善的主要因素。同时，气候变化和人类活动的相对作用的大小在不同空间位置上的

表现不同（图5），以人类活动为主（相对作用在50%以上）导致植被改善的区域面积占总改善面积的74.33%，主要分布于兴安盟东部、锡林郭勒盟、赤峰西部、河套平原、呼和浩特和包头以南地区以及山西西北部；而以气候变化为主导致植被改善的区域面积只占总改善面积的25.67%，主要分布于河北南部以及内蒙古科尔沁右翼中旗附近地区。此外，在通辽和鄂尔多斯附近等其他地区，气候变化和人类活动对植被恢复的作用则呈现出交互错杂的现象。从各类变化区域来看，在植被明显改善区域、中度改善区域和轻微改善区域，以人类活动为主导致植被改善的面积占各自区域总面积的比例都已超过70%，尤其在明显改善区域，达到82.45%。

2.3植被退化区气候变化和人类活动的相对作用大小

气候变化和人类活动不仅使华北地区部分区域的植被覆盖有所增加，同时也使部分区域的植被覆盖有所减少，生态环境恶化。本文中的植被退化区包括表2中的植被严重退化区域、中度退化区域和轻微退化区域。根据植被退化区域气候变化和人类活动相对Fig.4TrendsofNDVIchangeinNorthChina图41981—2006年华北地区NDVI变化趋势作用的3种可能情景，计算得到1981—2006年间气候变化和人类活动在华北地区植被退化区域的相对作用。从整体上来看，在植被退化区域，气候变化的相对作用为28.12%

，

1期刘斌等：华北地区植被覆盖变化及其影响因子的相对作用分析17

图51981—2006年华北地区气候变化和人类活动在植被改善区域的相对作用

coverinNorthChinafrom1981to

2006Fig.5Therelativeroleofclimatechangeandhumanactivitiesintheincreasedareaofvegetation

人类活动的相对作用为71.88%。由此可见，在1981—2006年间人类活动是导致该地区植被退化的主要因素。而且，二者相对作用的大小还存在着明显的空间差异。从图6中可以看出，以人类活动为主（相对作用在50%以上）而导致植被退化的区域面积占华北地区植被退化总面积的73.38%，主要分布在河北中部和东北部、内蒙古包头和呼和浩特以北地区、呼伦贝尔中南部、阿拉善左旗和乌海市附近地区，分布广泛；以气候变化为主而导致植被退化的区域面积占退化总面积的26.62%，主要分布在内蒙古阿拉善右旗和满洲里附近等地区，分布范围小且较分散。而其他植被退化区域，如北京、天津、山西中部、呼伦贝尔北部，气候变化和人类活动的相对作用则比较复杂。从各类变化区域来看，在植被严重退化区域、中度退化区域和轻微退化区域中，以人类活动为主导致植被退化的面积占各自区域总面积的比例都已超过70%

，由此也可以看出人类活动也是导致

图61981—2006年华北地区气候变化和人类活动在植被退化区域的相对作用

coverinNorthChinafrom1981to2006Fig.6Therelativeroleofclimatechangeandhumanactivitiesinthedecreasedareaofvegetation

18自然资源学报30卷华北地区植被退化的主要因素。

3讨论

3.1植被覆盖变化与气候因子的相关性

本文通过残差法剥离植被覆盖变化过程中气候变化的影响。根据残差法的原理可知，只有NDVI与气候指标显著相关时，才能较好地剥离气候变化对NDVI的影响。目前对NDVI与气候指标的相关分析研究主要集中在气温和降水这两个气候指标，而对具有显著生态学意义的气候指标与NDVI的关系研究则较少[36-38]。高大伟等指出NDVI与生态气候指标（生物热量指数和生物干湿度指数）具有较好的相关性[39]。因此，本文的气候指标不仅考虑了气温和降水，还考虑了干燥度指数、生物热量指数、生物干湿度指数和>0℃月均温之和4个综合气候指标来进行NDVI与气候指标的相关分析。结果发现，1981—2006年间干燥度指数和当年NDVI呈显著负相关，相关系数为-0.682，通过P

图71981—2006年华北地区NDVI与气候指标相关系数的空间分布

NorthChinafrom1981to2006Fig.7SpatialdistributionofthecorrelationcoefficientsbetweenNDVIandclimateindicatorsin

3.2人类活动对植被覆盖变化的影响

除气候变化因素外，人类活动也是植被覆盖变化过程中不可忽视的重要驱动因素。通过对华北地区NDVI残差序列进行趋势分析，得到NDVI残差年际变化趋势的空间分布，如图8所示。1981年以来，华北地区人类活动对植被覆盖变化的影响逐渐增强。一

方面表现为人类活动对植被覆盖的增加起到了促进作用，另一方面表现为人类活动对植

1期刘斌等：华北地区植被覆盖变化及其影响因子的相对作用分析19被覆盖起到了破坏作用。由图8统计可

知，NDVI残差变化趋势为正的区域占总

面积的55.16%，主要分布在内蒙古兴安盟

东部、赤峰南部、呼和浩特南部、河套平

原、河北东南部和山西北部地区，说明在

这些地区人类活动对植被生长起到了积极

作用。如在内蒙古河套平原，农药化肥的

大量使用以及农田灌溉水利设施的建设使

得NDVI的变化强烈依赖于引黄灌溉，对

降水的依赖性较小[12]；河北环抱京津，其

绿化工程是三北防护林体系的重点工程，

三北防护林工程实施以来，河北省项目区

森林覆盖率逐渐增大，“泽被当地，护卫图81981—2006年华北地区NDVI残差变化趋势

Fig.8ResidualtrendofNDVIinNorthChinafrom京津”的生态效果已经显现[17]；在山西北

1981to2006部和西南部，三北防护林生态工程、天然

林保护工程和雁门关生态恢复项目的建设以及长期的退耕还林有利于水土保持，植被恢复情况良好，绿化荒山、封山育林、植树造林等政策也起到积极作用[41]。NDVI残差变化趋势为负的区域占总面积的44.84%，主要分布在大兴安岭中北部、呼伦贝尔大部分地区、科尔沁沙地、阿拉善左旗、阿拉善右旗、克什克腾旗附近地区、呼和浩特北部、鄂尔多斯西部、巴彦淖尔高原、山西中南部，此外，在河北中东部、北京、天津也有零星分布，这些地区人类活动对植被生长起到了破坏作用。如在科尔沁沙地，近年来人口的急剧增加和过度垦荒、超载放牧，使得草地退化严重，而且还加剧了土壤风蚀[42]；阿拉善右旗南部由于20世纪50年代以来从甘肃民勤县涌入大批农民，在沙质土地上挖野菜、药材，使本就脆弱的生态环境和植被覆盖遭到严重破坏[12]；山西中部地区由于近年来工业化发展迅速，造成植被状况有所恶化[41]；京津唐地区由于快速的城市化进程，导致城市规模扩大，人口激增，绿地和原始土地被大量用于住宅建设和开发区建设，植被覆盖有所减少[43-44]。当然，本文是在大尺度上区分了华北地区植被覆盖变化过程中人类活动的信号，要进行具体的植被建设恢复工作，还需要在更小的尺度上进一步实地调查验证。4结论

本文利用1981—2006年的GIMMSNDVI时间序列数据，结合华北地区及其周边223个气象站点的气候数据，对该地区的植被覆盖时空变化进行了分析。在此基础上，分别对植被改善区和植被退化区构建了不同情景下气候变化和人类活动两大影响因素相对作用的评价方法，对气候变化和人类活动在华北地区植被变化过程中的相对作用进行了研究，主要得出以下结论：①在气候变化和人类活动的共同作用下，1981—2006年间华北地区植被NDVI呈明显上升趋势，植被覆盖状况整体上有所改善，但却存在明显的空间差异，且植被改善区面积小于植被退化区面积。②1981—2006年间，华北地区植被ND-VI与干燥度指数和年平均气温具有很好的相关性，且与干燥度指数呈显著负相关，与年平均气温呈显著正相关，说明气候变化是影响华北地区植被覆盖变化的重要因素。③

除

20自然资源学报30卷去气候因子对植被覆盖变化的影响外，人类活动对华北地区植被覆盖变化也起到了重要作用。1981—2006年间，在华北地区植被改善区域，人类活动的相对作用为68.10%，远大于气候变化对植被增加的作用（31.90%），而且以人类活动为主（相对作用在50%以上）导致植被改善的区域面积占改善区总面积的74.33%。在植被退化区域，人类活动的相对作用为71.88%，也是大于气候变化的相对作用（28.12%），同样以人类活动为主（相对作用在50%以上）导致植被退化的区域面积占退化区总面积的73.38%。由此可见，人类活动无论在华北地区植被覆盖改善过程中还是退化过程中都已经成为主导因素。参考文献(References)：

[1]孙红雨,王长耀,牛铮,等.中国地表植被覆盖变化及其与气候因子关系——基于NOAA时间序列数据分析[J].遥

感学报,1998,2(3):204-210.[SUNHong-yu,WANGChang-yao,NIUZheng,etal.Analysisofthevegetationcov-erchangeandtherelationshipbetweenNDVIandenvironmentalfactorsbyusingNOAAtimeseriesdata.JournalofRemoteSensing,1998,2(3):204-210.]

[2]信忠保,许炯心,郑伟.气候变化和人类活动对黄土高原植被覆盖变化的影响[J].中国科学D辑:地球科学,2007,

37(11):1504-1514.[XINZhong-bao,XUJiong-xin,ZHENGWei.ResponseofvegetationcoverchangetoclimatechangeandhumanactivitiesinLoessPlateau.ScienceinChinaSeriesD:EarthSciences,2007,37(11):1504-1514.]

[3]孙智辉,刘志超,雷延鹏,等.延安北部丘陵沟壑区植被指数变化及其与气候的关系[J].生态学报,2010,30(2):

533-540.[SUNZhi-hui,LIUZhi-chao,LEIYan-peng,etal.ThevariationsofNDVIandtherelationwithclimateinhillyandgullyregionofNorthernYan