基于MODIS NDVI和多方法的青藏高原植被物候时空特征分析

doi:10.11873/j.issn.1004-0323.2018.3.0486

遥感技术与应用

2018, Vol. 33

Issue (3): 486-498 DOI: 10.11873/j.issn.1004-0323.2018.3.0486

遥感应用

基于MODIS NDVI和多方法的青藏高原植被物候时空特征分析

周玉科¹，刘建文²

(1.中国科学院地理科学与资源研究所生态系统网络观测与模拟院重点实验室，北京 100101；
2.福州大学空间信息工程研究中心数据挖掘与信息共享教育部重点实验室，福建福州 350003)

Spatio-temporal Analysis of Vegetation Phenology with Multiple Methods over the Tibetan Plateau based on MODIS NDVI Data

Zhou Yuke¹，Liu Jianwen^2

(1.Ecology Observing Network and Modeling Laboratory，Institute of Geographic and Nature Resources Research，Chinese Academy of Sciences，Beijing 100101，China；2.Key Laboratory of Spatial Data Mining and Information Sharing of Ministry of Education，Fuzhou University，Spatial Information Research Center of Fujian Province，Fuzhou 350002，China)

全文: PDF(12408 KB)

摘要：

基于遥感的植被物候分析具有连续观测、覆盖面广的优点，已经成为全球气候变化背景下植被生理物理变化研究的重要手段。基于青藏高原96个气象站点和对应的MODIS NDVI时序数据，探索分析了2000~2014年间站点、典型植被覆盖和高原整体区域上的植被物候变化趋势及地理环境影响。首先采用三次样条(Spline)函数法、双逻辑斯蒂(D\|L)函数法和奇异谱分析法(SSA)进行NDVI时序数据重构，其次利用导数法(Derivative method，Der)及阈值法(Threshold method，Trs)提取植被物候关键参数(生长季长度LOS、生长季开始点SOS、生长季结束点EOS)，进而对比分析6种方法提取结果的差异和适用条件；再利用M\|K test趋势分析法计算了各站点和区域的物候发展趋势，并探索其与海拔、降水、温度间的相关性；最后将本文遥感LOS与气温阈值法得到的生长期长度指标(GSL)进行了对比验证。结果表明：①阈值法有较强的稳定性及适应性，导数法结果差异较大，整体来看阈值法优于导数法。草原及森林覆被类型中3种物候参数使用SSA、Spline、D\|L与阈值法的组合提取结果较优；②多种方法提取的多种物候参数趋势空间差异性较大，小尺度上趋势比较一致。青藏高原东南部湿润半湿润的灌木草原区及西北荒漠草原区，呈现SOS和EOS推迟、但LOS仍然延长的趋势；青藏高原西南部湿润区呈现SOS延后、EOS提前、LOS缩短的趋势；青藏高原分布广泛的草原区，各物候参数并未呈现显著变化的趋势特征；③温度与物候参数发展趋势呈相关关系，温度升高的情况下呈现SOS提前、EOS延后的现象。由于高原地貌、气候的复杂性，多数站点的物候发展趋势与海拔、降水的相关性并不显著，GSL与LOS的相关性也呈现相似的特征。本文证实在植被物候参数遥感提取中，多种方法差异较大，很少有在全部区域均适用的方法；站点尺度上青藏高原在近15年来植被物候并未呈现明显春季提前或秋季延后趋势。

关键词： 青藏高原; NDVI时序重构; 植被物候参数; 趋势分析; 阈值法; 导数法

Abstract: In the context of global climate change，vegetation phenology analysis based on remote sensing has become an critical method for studying the characteristics of physical and physiological changes of vegetation.This paper uses the MODIS NDVI time\|series data of 96 meteorological stations over the Tibetan Plateau during 2000\|2014 to explore the development trend of vegetation phenological and geographical environment factors of each meteorological station，typical vegetation coverage and the whole plateau region.Firstly，using three cubic spline function method (Spline)，double logistic function method(D\|L)and singular spectrum analysis (SSA)，NDVI time\|series data is reconstructed，then using the derivative method (Der)and threshold method (Trs)，the key parameters of phenological information is extracted，after that differences and application conditions between the six methods are analyzed and compared.Secondly，using M\|K test trend analysis method，the phenological development trend of each site and area were calculated，the relationship between phenological development trend and altitude，precipitation，temperature is studied.Finally，by the Growing season length(GSL)obtained by temperature threshold method，LOS is compared and verified.in grassland and forest land cover types，SSA，Spline，D\|L combined with threshold method to get the Start of Season(SOS)，end of season(EOS)，Length of season (SOS)respectively is a good combination strategy.(2)The spatial differences of various phenological parameters extracted by different methods are large，and the trend is relatively consistent at small scales.Southeast humid and semi\|humid shrub steppe region and northwestern desert steppe in the Tibetan Plateau，SOS and EOS delayed，but LOS prolonged；southwestern humid region，SOS and EOS delayed，LOS shortened；widely distributed grassland，the phenological parameters did not show significant tendency.(3)Temperature is related to the development trend of phenological parameters.With temperature increasing，the phenomena of SOS advance，EOS lag are presented.Because of the complexity of the plateau landform and climate，there was no significant relationship between phenological development trend for most of the site with the altitude and precipitation，only a few sites have strong correlation，the correlation between GSL and LOS also showed similar characteristics.For remote sensing based phonological analyses，our study identify there is no method existing here that is a adaptive across all the Tibetan Plateau.in addition，at point scale the phenological parameters do not represent a significant earlier or later trend.

Key words: Tibetan Plateau；NDVI time\ series reconstruction；Vegetation phenological trend analysis；Threshold method；Derivative method

收稿日期: 2017-05-14 出版日期: 2018-07-04

TP 79

基金资助: 国家自然科学基金项目(41601478)，国家重点研发计划项目(2016YFC0500103)，中国科学院STS项目(KFJ\|SW\|STS\|167)，资源与环境信息系统国家重点实验室开放基金项目(LREIS，2016)。

作者简介: 周玉科(1984-)，男，山东济宁人，博士，助理研究员，主要从事生态遥感研究。Email：zhouyk@igsnrr.ac.cn。

	服务
	把本文推荐给朋友
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	周玉科
	刘建文

引用本文:

周玉科，刘建文. 基于MODIS NDVI和多方法的青藏高原植被物候时空特征分析[J]. 遥感技术与应用, 2018, 33(3): 486-498.

Zhou Yuke，Liu Jianwen. Spatio-temporal Analysis of Vegetation Phenology with Multiple Methods over the Tibetan Plateau based on MODIS NDVI Data. Remote Sensing Technology and Application, 2018, 33(3): 486-498.

链接本文:

http://www.rsta.ac.cn/CN/10.11873/j.issn.1004-0323.2018.3.0486 或 http://www.rsta.ac.cn/CN/Y2018/V33/I3/486

[1]	孟梦，牛铮. 近30 a内蒙古NDVI演变特征及其对气候的响应[J]. 遥感技术与应用, 2018, 33(4): 676-685.
[2]	陈思宇，巩垠熙，梁天刚. 星载激光雷达在青藏高原湖泊变迁中的应用研究[J]. 遥感技术与应用, 2018, 33(2): 351-359.
[3]	马敏娜，袁文平. 青藏高原总初级生产力估算的模型差异[J]. 遥感技术与应用, 2017, 32(3): 406-418.
[4]	王丽娟,郭铌,王玮，芦亚玲，沙莎. 基于TESEBS模型估算高原地区地表蒸散发[J]. 遥感技术与应用, 2017, 32(3): 507-513.
[5]	黄田进，梁丁丁，贾立，张静潇,卢静，周杰. 青藏高原地区湖泊面积插补迭代自动提取#br#[J]. 遥感技术与应用, 2017, 32(2): 289-298.
[6]	杨志刚,达娃,除多. 近15 a青藏高原积雪覆盖时空变化分析[J]. 遥感技术与应用, 2017, 32(1): 27-36.
[7]	冯莉,李柳华,郭松,卢荻. HJ-1A NDVI与MODIS NDVI时间序列提取植被物候特征对比研究[J]. 遥感技术与应用, 2016, 31(6): 1158-1166.
[8]	安培浚,高峰,王立伟. 青藏高原冰川、积雪与地质灾害空间观测研究态势分析[J]. 遥感技术与应用, 2016, 31(6): 1223-1230.
[9]	封敏,李再明,邱炳文,王重洋,罗钰涵. 21世纪初中国农作物与自然植被物候时空差异[J]. 遥感技术与应用, 2016, 31(5): 1003-1012.
[10]	王玉丹,南卓铜,陈浩,吴小波. 基于K 最近邻模型的青藏高原CMORPH日降水数据的订正研究[J]. 遥感技术与应用, 2016, 31(3): 607-616.
[11]	杨成松,车涛,欧阳斌. 青藏高原地表温度时空变化分析[J]. 遥感技术与应用, 2016, 31(1): 95-101.
[12]	肖林,车涛. 青藏高原积雪对气候反馈的初步研究[J]. 遥感技术与应用, 2015, 30(6): 1066-1075.
[13]	唐志光,王建,梁继,李朝奎,王欣. 基于MODIS的青藏高原雪线高度遥感监测[J]. 遥感技术与应用, 2015, 30(4): 767-774.
[14]	杨麦,王雪梅. 基于GIS的图书馆文献计量分析[J]. 遥感技术与应用, 2015, 30(4): 819-824.
[15]	赵志强,李爱农,边金虎,郭文静,刘倩楠,赵伟,黄成全. 山区可见光—近红外遥感影像浓密植被暗像元自动识别方法研究[J]. 遥感技术与应用, 2015, 30(1): 58-67.

Viewed

Full text

Abstract

Cited

Shared

Discussed