2 中国气象局气象干部培训学院, 北京100081;
3 北京师范大学, 北京100875
2 China Meteorological Administration Training Center, Beijing 100081;
3 Beijing Normal University, Beijing 100875
植被是陆地生态系统的主体,在保持水土、调节大气成分以及减缓温室气体浓度上升和维持气候稳定等方面具有不可替代的作用(Solomon and Shugart, 1993),植被变化对全球能量循环及物质的生物化学循环具有重要的影响(Liu et al., 2002; Piao et al., 2003),是全球变化研究的热点之一。遥感技术提供了长时间序列、大范围、高精度的地面植被信息,成为研究植被覆盖变化的利器,已有许多研究利用遥感数据对植被覆盖进行了监测分析和研究(Myneni et al., 1997; Menzel and Fabian, 1999; Cao et al., 2003; 张佳华等,2003; Zhou et al., 2003; 陈朝晖等,2004;丹利和谢明,2009;Gobron et al., 2010)。随着MODIS的出现,发展了一个新的植被指数:增强型植被指数EVI(Hutet et al., 1997)。EVI具备了新的特征,其优势在于减少了由土壤背景和残留气溶胶造成的污染,不易饱和,可以容忍更大的冠层密度等(Huet et al., 1997;王正兴等,2003;Zhang et al., 2004)。
中国北方(东北、华北、西北)是我国经济生活的关键区域,人口压力大,水资源匮乏,生态环境脆弱,其植被覆盖变化不仅受到气候因子的驱动,同时也受到人类活动的影响。21世纪以来气候变化更加复杂,人类活动日益剧烈。一方面国家绿化和生态建设等大规模展开,例如退耕还林是迄今为止我国投资最大、涉及面最广、群众参与度最高的一项生态建设工程,根据国家林业局统计数据(http://www.forestry.gov.cn[2012-12-01]),截至2009年工程已累计完成退耕地造林906.26×104 km2,给工程重点区域内的植被覆盖带来巨大变化;但另一方面,伴随着经济的发展,城市建设和工业建设进程不断加快,建筑用地扩增,牧地林地退化,也使一些地区植被覆盖朝着恶性方向发展(Houghton,2002; Fang et al., 2004)。北方地区植被覆盖更加深刻地记录了人类活动的烙印,在深度、广度上都 发生了巨大变化,需要及时分析中国北方地区植 被变化的阶段性、时空性。本文利用较高精度的植被覆盖数据描述了中国北方地区植被覆盖总体特征,分析了21世纪初中国北方植被的时空变化特征,研究了不同等级植被变化特征和单位面积植 被指数变化,可为生态工程建设和土地利用规划提供参考。
2 数据及分析方法 2.1 研究区概况中国北方地区包括东北、华北、西北地区15省(直辖市、自治区)(黑龙江、吉林、辽宁、内蒙古、北京、天津、河北、河南、山东、山西、陕西、甘肃、宁夏、新疆、青海),位于(33.24°N~53.28°N,75.14°E~132.58°E),人口占全国的41.39%,粮食种植面积占全国54.76%,棉花种植面积占全国73.43%,油料种植面积占全国42.82%(中国统计年鉴编委会,2009),是我国粮食生产的关键区域,在中国的经济生活中起着举足轻重的作用。中国北方地区大部分属于温带大陆性气候区,夏季盛行东南季风,冬季受西伯利亚高压控制。利用欧洲航天局(ESA)2004~2006年全球土地利用产品(Arino et al., 2007),在原分类体系基础上,考虑中国北方地区具体情况,将土地利用类型分为农田、常绿林、落叶林等9类。
图 1显示,研究区域东部以林地、农田为主,其间有一些草甸灌木分布,植被覆盖度较高。林地主要分布在内蒙古东北部、东北以及陕南等地区;农田主要分布在东北平原、华北平原、黄土高原和河谷区;山东、河北、河南、山西等大部分为农田用地。研究区中西部地区由于高海拔和气候条件限制,植被覆盖较低,内蒙古中西部、甘肃北部至新疆大部分地区为荒漠裸地区域。甘肃南部、青海东部南部和新疆南部有农田和草原灌木的存在。
![]() |
图1 中国北方地区土地利用类型 Fig. 1 Land use types of the northern part of China |
本研究所用数据来自NASA地球观测系统(EOS)提供的MODIS陆地植被指数产品MOD13A3。数据包含两种植被指数:归一化植被指数(NDVI)和增强型植被指数(EVI),时间分辨率为月,空间分辨率为1 km×1 km。MODIS-EVI经过了全面的大气校正,对残留的气溶胶做了进一步校正,纠正了土壤背景值的影响,不易饱和,可以容忍更大的冠层密度,Wang et al.(2010)利用中国北方地区两个典型生态系统通量观测数据和两种植被指数进行对比分析,发现在研究时段内MODIS-EVI比MODIS-NDVI更不易受到雪盖背景的影响,并更好地表征了实际参加光合作用的植被部分,认为MODIS-EVI更适于描述中国北方植被状况,故本文选用2000~2009年MODIS13A3-EVI数据,经过图像融合、投影转换、剪切等步骤生成中国北方EVI时序数据集。
本研究使用区域平均EVI表征当年植被平均生长状况;植被指数变化趋势EVItrend表征研究时段内当地植被生长年际变化率,正值表示植被改善,负值表示植被恶化;参考Myneni et al.(1997)关于无植被地区的定义,定义EVI小于0.05的地区为植被极稀疏地区或无植被地区,大于0.05的地区为有植被覆盖地区,在此标准下研究植被覆盖区面积变化及植被覆盖比值。
3 植被覆盖变化 3.1 植被整体分布特征多年平均中国北方地区植被覆盖的空间分布如图 2所示。北方15省面积辽阔,气候条件不同,植被覆盖空间差异显著。东部植被覆盖好、EVI数值大,西部植被覆盖差、EVI数值小,沿内蒙古—青海的东北—西南走向有一条较为明显分界线。分界线以东水热配置较好、气候适宜,植被呈现良好生长状态,内蒙古东北部、东北3省东部、山东、河南、山西一带年平均EVI达到0.3以上,夏季达到0.5以上。分界线以西的大部分地区处于内陆干旱、半干旱气候条件下,干旱少雨风大,沙漠戈壁多,不利于植被生长。内蒙古西南部、甘肃西北部、新疆大部分、青海西北部年平均EVI小于0.05,属于本文定义植被极稀疏地区或无植被地区。只有位于伊犁河谷的新疆塔城、伊宁等地水资源较为丰富,有草原草甸和针叶林的分布,以及青海湖周围高寒草原、高寒草甸的存在使当地植被指数略高。
![]() |
图2 2000~2009年中国北方地区(a)夏季平均和(b)全年平均EVI空间分布 Fig. 2 Spatial distributions of mean Enhanced Vegetation Index (EVI) over the northern part of China in (a) summer and (b) all year around from 2000 to 2009 |
多年平均中国北方地区植被覆盖的时间分布特征如图 3所示。植被随着四季轮回而生长繁茂衰弱,有明显的季节周期。中国北方地区夏季水热同期,植被蓬勃生长EVI最高,秋季温带森林落叶,作物逐渐收割,中国北方地区EVI整体下降,只有陕西南部残留一部分亚热带森林。冬季植被覆盖稀少,低EVI值多来源于城市绿化常绿植物以及由于气候变暖而提前萌发的植被。春季冰川解冻植被开始复苏,EVI又开始一轮增长。EVI逐月变化图显示中国北方地区3~7月为EVI上升期,8月至次 年2月为EVI下降期,8月为全年植被生长最旺盛的月份,5~10月为主要生长季节。
![]() |
图3 2000~2009年中国北方地区平均EVI的(a)月变化和(b)季节变化 Fig. 3 (a) Monthly variation and (b) seasonal variation of mean EVI over the northern part of China from 2000 to 2009 |
区域平均EVImean反映了一年内中国北方地区植被的平均生长状况,区域最大EVImax反映了一年内中国北方地区植被在水热条件适宜气候下最优生长状况。图 4显示了两者的年际变化特征。2000~2009年无论是区域平均EVI还是区域最大EVI都呈现上升趋势,表明中国北方地区植被活动在增强。10年间区域平均EVI呈现3个不同时期,2000~2001年为稳定期,2002~2008年为明显上升期,2009年出现反复,植被生长状况较差。10年间区域平均EVI增幅为0.009,约为5.97%。区域最大EVI受到每年温度、水分、光照等配置的影响,在不断波动中上升。2004、2005、2007年为较高时 期,2001、2006、2009年则为较低时期。10年间区域最大EVI增幅为0.035。
![]() |
图4 2000~2009年中国北方地区区域平均EVI和区域最大EVI变化特征 Fig. 4 Variation characteristics of mean EVI and max EVI over the northern part of China from 2000 to 2009 |
分析2000~2009年中国北方地区EVI的变化趋势(图 5),取5%的显著性水平,10年间植被显著改善的格点个数为1618544,占区域面积的36.19%,植被显著退化的格点个数为577013,占区域面积的12.9%,其余地区植被生长状态保持原来水平,无显著变化。可见21世纪最初10年间中国北方地区植被总体朝着改善方向发展的同时,部分地区植被产生退化或恶化。
![]() |
图5 2000~2009年中国北方地区EVI变化趋势EVItrend(白色区域植被变化未通过α=0.05显著性检验) Fig. 5 Spatial distributions of the trends in EVI over the northern part of China from 2000 to 2009 (white area: not significant at 0.05 significance level) |
中国北方地区EVI变化趋势显示了较大的空间异质性(图 5)。EVI显著增加的区域主要包括东北平原、华北平原、新疆农业生产区以及黄土高原地区,这与退耕还林还草、重建生态家园、建设高效农业等工程密切相关。内蒙古中部退化草原也得到一定程度的恢复,但部分地区植被却进一步恶化,例如阿拉善高原、乌兰察布高原。EVI显著减少的区域大多位于西部地区,但值得注意的是内蒙古东北部的常绿针叶林、黑龙江中部针阔混交林也出现了明显的退化趋势。EVI增加、减少区域的空间异质性可能源于城市化、气候影响、生态建设工程、农业生产活动等多重因子的综合作用。
3.3 季节变化特征植被生长、繁茂、衰弱有明显的季节周期,虽然每个季节植被生长条件不尽相同,但2000~2009年10年间中国北方地区植被春夏秋冬逐季节平均EVI均呈现上升趋势(图 6)。其中春季、秋季上升幅度小,夏季、冬季上升幅度大。中国北方地区夏季植被生长活动增强直接导致夏季EVI上升幅度大,而冬季EVI上升幅度大则可能与气候变暖、春季提前,植被生长开始时间提前、结束时间延后有关系。
![]() |
图6 2000~2009年中国北方地区平均EVI(a)春季、(b)夏季、(c)秋季、(d)冬季变化特征 Fig. 6 Variation characteristics of monthly EVI over the northern part of China from 2000 to 2009 in (a) spring, (b) summer, (c) autumn, and (d) winter |
为了分析各个季节EVI增长的具体原因,图 7给出中国北方地区EVI逐季节变化趋势的空间分布。图 7a显示春季EVI在华北平原地区有明显增加,这可能与气候变暖、华北平原春耕提前密切相关。东北平原则以一年一熟作物为主,春季尚未开始耕作,所以气温升高对春季农业耕作影响不大。东北地区春季EVI呈现整体下降趋势,常绿针叶林、针阔混交林退化趋势明显。
![]() |
图7 2000~2009 年中国北方地区EVI(a)春季、(b)夏季、(c)秋季、(d)冬季变化趋势EVItrend(白色区域植被变化未通过α=0.05显著性检 验) Fig. 7 Spatial distributions of the trends in monthly EVI over the northern part of China from 2000 to 2009 (white area: not significant at 0.05 significance level) in (a) spring, (b) summer, (c) autumn, and (d) winter |
中国北方地区植被生长最旺盛,植被指数最高的季节为夏季。东北平原、华北平原等主要粮食产区植被生长更加繁茂,东北大小兴安岭林区EVI数值可以达到0.7左右,陕西南部有林区存在,EVI数值可达0.5~0.7。西北部地区新疆西北部天山山麓、青海西南部及青海湖周围植被生长旺盛。10年来夏季EVI变化趋势(图 7b)显示中国北方大部分地区植被生长增强,少部分呈现退化趋势的区域为东北林区、内蒙古部分草原地区、新疆北部。
图 7c显示华北平原植被活动增强是10年来中国北方地区秋季EVI整体上升的主要原因,可能源于农作物增产、收割时间退后或种植方式改变。西部及北京、天津、河北地区植被略有下降。图 7d显示华北平原、内蒙古部分草原以及新疆北部植被活动增强是10年来中国北方地区冬季EVI整体上升的主要原因。
3.4 不同等级植被变化特征根据3.2节中的定义,利用区域年平均植被指数EVImean的年际变化和季节变化可以表征中国北方地区整体的植被发展状况,但平均EVI的使用也忽略了不同等级植被变化趋势。因为不同等级植被的变化趋势有可能与平均EVI的变化趋势不同甚至相反,所以有必要将EVI分为不同等级考察中国北方地区植被变化。
本文定义EVI小于0.05的地区为植被极稀疏地区或无植被地区(Myneni et al., 1997),则可将年平均EVI分为≤0.05、0.05~0.1、0.1~0.2、0.2~0.3、>0.3共5个等级,以区分不同植被覆盖情况。其中前两个等级为植被稀疏区域,第三和第四等级为植被一般生长区域,第五等级为植被茂盛生长区域。图 8给出了2000~2009年中国北方地区不同等级植被所占比例变化趋势,图中显示植被稀疏地区共占总面积40%左右,植被一般生长区域约占50%左右,植被茂盛生长区域占5%~10%。
![]() |
图8 2000~2009年中国北方地区不同等级植被所占比例变化趋势 Fig. 8 Variation characteristics of EVI in different grades over the northern part of China from 2000 to 2009 |
EVI小于0.05的地区为植被极稀疏地区或无植被地区,主要位于西部,约占总面积的20%左右(见图 9),EVI在0.05~0.1的地区是植被稀疏区多为荒漠草原等,约占总面积20%左右,10年来在波动中有所下降。EVI在0.1~0.3的地区为农业生产区和草甸草原等自然植被类型,EVI大于0.3的区域为林地或高效农业区,植被覆盖好,植被生长旺盛,所占总面积比例在波动中上升。
![]() |
图9 2000~2009年中国北方地区不同时段不同等级植被比例 Fig. 9 Ratio of EVI in different grades over the northern part of China from 2000 to 2009 |
定义中国北方地区总面积为单位1,使用区域年平均EVI除以植被覆盖面积百分比,得到中国北方地区单位面积EVI。图 10中2000~2009年中国北方地区单位面积EVI呈现明显的上升趋势,表明近10年间中国北方地区植被覆盖面积增加的同时,单位面积EVI增加,植被生长更旺盛。
![]() |
图10 2000~2009年中国北方地区单位面积EVI变化 Fig. 10 Variation characteristic of EVI per unit area over the northern part of China from 2000 to 2009 |
在气候变化和人类活动的双重影响下,中国北 方地区的植被覆盖在21世纪最初10年呈现上升趋势,区域年平均EVI增加0.009。东北平原、华北平原、黄土高原和新疆农业区显著增加,内蒙古东北部常绿针叶林、黑龙江中部针阔混交林退化和部分西 部地区植被退化,表明EVI变化趋势空间异质性大。
2000~2009年中国北方地区植被覆盖的动态变化呈现了季节特征,虽然逐季节平均EVI均呈现上升趋势,但幅度并不相同。春季、秋季上升幅度小,夏季、冬季上升幅度大,分析其原因,夏季植被生长增强直接导致夏季EVI上升幅度大,而冬季EVI上升幅度大则可能与全球变暖后春季提前,植被生长时间提前有关系。2000~2009年中国北方地区植被覆盖低的面积有减少趋势,植被覆盖好的面积有增加趋势,单位面积EVI增加,植被生长更旺盛。气温升高、降水变化是植被覆盖变化的重要原因,人类活动驱动下的农业生态工程同样也是植被变化的重要影响因素,21世纪最初10年来中国北方地区植被覆盖变化特征分析表明,近年来政府倡导下进行的退耕还林还草等植被建设工作带来的生态效益正在不断呈现。
[1] | Arino O, Leroy M, Ranera F, et al. 2007. GLOBCOVER-A Global Land Cover Service with MERIS [C]./ Proceedings of the ENVISAT Symposium. Montreux, Switzerland, 6pp. |
[3] | Cao M K, Prince S D, Li K R, et al. 2003. Response of terrestrial carbon uptake to climate interannual variability in China [J]. Global Change Biology, 9 (4): 536-546. |
[5] | 陈朝晖, 朱江, 徐兴奎. 2004. 利用归一化植被指数研究植被分类、面积估算和不确定性分析的进展 [J]. 气候与环境研究, 9 (4): 687-696. Chen Zhaohui, Zhu Jiang, Xu Xingkui. 2004. Recent research progress in uncertainties of the normalized difference vegetation index, land cover classification and vegetation fraction estimation [J]. Climatic and Environmental Research (in Chinese), 9 (4): 687-696. |
[7] | 丹利, 谢明. 2009. 基于MODIS资料的贵州植被叶面积指数的时空变化及其对气候的响应 [J]. 气候与环境研究, 14 (5): 455-464. Dan Li, Xie Ming. 2009. The spatiotemporal variation of leaf area index in Guizhou and its response to climate based on MODIS data [J]. Climatic and Environmental Research (in Chinese), 14 (5): 455-464. |
[9] | Fang J Y, Piao S L, He J S, et al. 2004. Increasing terrestrial vegetation activity in China, 1982-1999 [J]. Science in China (Ser. C), 47 (3): 229-240. |
[11] | Gobron N, Belward A, Pinty B, et al. 2010. Monitoring biosphere vegetation 1998-2009 [J]. Geophys. Res. Lett., 37: L15402, doi:10.1029/ 2010GL043870. |
[13] | Huet A R, Liu H Q, Batchily K, et al. 1997. A comparison of vegetation indices over a global set of TM images for EOS-MODIS [J]. Remote Sensing of Environment, 59: 440-451. |
[15] | Houghton R A. 2002. Temporal patterns of land-use change and carbon storage in China and tropical Asia [J]. Science in China (Ser. C), 45 (Suppl.): 10-17. |
[17] | Liu J Y, Liu M L, Deng X Z, et al. 2002. The land use and land cover change database and its relative studies in China [J]. Journal of Geographical Sciences, 12 (3): 275-282. |
[19] | Menzel A, Fabian P. 1999. Growing season extended in Europe [J]. Nature, 397: 659. |
[21] | Myneni R B, Keeling C D, Tucker C J, et al. 1997. Increased plant growth in the northern high latitudes from 1981 to 1991 [J]. Nature, 386: 698-702. |
[23] | Piao S L, Fang J Y, Zhou L M, et al. 2003. Interannual variations of monthly and seasonal normalized difference vegetation index (NDVI) in China from 1982 to 1999 [J]. J. Geophys. Res., 108(D14): 4401, doi:10.1029/ 2002JD002848. |
[25] | Solomon A M, Shugart H H. 1993. Vegetation dynamics and global change [M]./ Cramer W P, Leemans R. Assessing Impacts of Climate Change Vegetation Using Climate Classification System. London: Chapman and Hall, 190-217. |
[27] | Wang Z, Xiao X M, Yan X D. 2010. Modeling gross primary production of maize cropland and degraded grassland in northeastern China [J]. Agricultural and Forest Meteorology, 150: 1160-1167. |
[29] | 王正兴, 刘闯, Alfredo H. 2003. 植被指数研究进展: 从AVHRR-NDVI到MODIS-EVI [J]. 生态学报, 23: 979-987. Wang Zhengxing, Liu Chuang, Alfredo H. 2003. From AVHRR-NDVI to MODIS-EVI: Advances in vegetation index research [J]. Acta Ecologica Sinica (in Chinese), 23: 979-987. |
[31] | 张佳华, 侯英雨, 毛飞, 等. 2003. 基于遥感的华北地区作物干旱时空监测研究 [J]. 自然灾害学报, 12 (2): 255-259. Zhang Jiahua, Hou Yingyu, Mao Fei, et al. 2003. Study on spatial-temporal distribution of crop drought in North China based on remote sensing information [J]. Journal of Natural Disasters (in Chinese), 12 (2): 255-259. |
[33] | Zhang X Y, Fredl M A, Schaaf C B, et al. 2004. Climate controls on vegetation phonological patterns in northern mid and high latitudes inferred from MODIS data [J]. Global Change Biology, 10: 1133-1145. |
[35] | 中国统计年鉴编委会. 2009. 中国统计年鉴 [M]. 北京: 中国统计出版社. China Statistical Yearbook Editorial Board. 2009. China Statistical Yearbook [M] (in Chinese). Beijing: China Statistics Press. |
[37] | Zhou L M, Kaufmann R K, Tian Y, et al. 2003. Relation between interannual variations in satellite measures of northern forest greenness and climate between 1982 and 1999 [J]. J. Geophys. Res., 108(D1), 4004, doi:10. 1029/2002JD002510. |