基于MODIS数据的青藏高原生长季植被NDVI的时空变化分析.pdf

2024年3月第41卷第2 期枣庄学院学报JOURNALOFZAOZHUANGUNIVERSITYMar.2024Vol.41 NO.2基于MODIS数据的青藏高原生长季植被NDVI的时空变化分析杜宜臻，崔星岩，倪杰（枣庄学院旅游与资源环境学院，山东枣庄2 7 7 16 0）【摘要】利用2 0 0 1一2 0 2 0 年分辨率1kmMODIS归一化植被指数（NDVI）产品探讨青藏高原生长季植被的时空变化特征及规律。研究发现青藏高原生长季植被NDVI在空间上具有异质性，其数值从东南向西北逐渐递减，多年均值为0.14；近2 0 a来NDVI的稳定生长区和改善区面积达到94%，青藏高原整体呈现出变绿趋势；NDVI与海拔呈负相关关系，平均海拔每升高10 0 0 m，ND V I 下降0.0 9。研究结果为保护青藏高原生态系统等方面提供一定的科学参考。【关键词青藏高原；归一化植被指数；变化规律；海拔【中图分类号TP75文献标识码】A文章编号10 0 4-7 0 7 7（2 0 2 4）0 2-0 0 49-0 70弓引言青藏高原（以下简称高原）是我国乃至世界重要的生态区域，面积辽阔、地理位置特殊，生态环境又极为脆弱与敏感1-2。当前，气候变暖给全球生态系统带来极大影响，尤其是青藏高原地区生态环境发生明显的变化，比如冻土消融、冰川退缩、湿地菱缩、湖泊扩张等现象3-5】，这会改变高原地区的土壤水热平衡，进而影响到高原生态系统的结构、稳定性等7-8，并最终影响到周边地区乃至全球的气候变化。因此，研究青藏高原地区植被的变化规律是认识青藏高原生态环境特征及其演变规律的关键，对未来生态环境保护具有重要参考价值。归一化植被指数（normalized difference vegetation Index，ND VI）作为植被空间分布密度以及植物生长状态的最佳指示因子，能够准确、连续地反映出植被的生长变化状况，是反应植被分布与植被覆盖的重要数据。近年来归一化植指数被广泛应用于植被动态变化及其影响因素的研究当中。Kawamura 等9对内蒙古锡林郭勒中部草原区进行研究，利用回归分析等方法得出MODIS-NDVI相比AVHRR-NDVI时间剖面更具有真实性。王思等10】基于MODIS-NDVI数据对广东省的植被覆盖度进行提取研究，结果表明处于中高植被覆盖度及以上的区域达到8 5%，且植被空间集聚性特征显著。董璐等11利用NDVI数据研究了新疆地区植被生长对气温的响应。青藏高原作为亚洲乃至北半球气候变化的“感应器”和“敏感区”，对中国乃至亚洲生态安全具有重要的屏障作用。当前气候变化给高原植被生长带来了一系列的影响，大量学者基于NDVI这一关键特征指标来研究青藏高原地区的植被演变特征和规律。已有研究表明，高原生长季NDVI总体变化呈上升趋势，干旱地区植收稿日期2 0 2 3-0 6-13【基金项目】山东省自然科学基金青年项目（ZR2022QD044，ZR2 0 2 2 Q D 0 16）；冰冻圈科学国家重点实验室开放基金(SKLCS-0P=2021-04)。【作者简介】杜宜臻（1992 一），女，山东临沂人，枣庄学院旅游与资源环境学院讲师，理学博士，主要从事土壤水热方面研究。通信作者：倪杰（1990 一），男，山东枣庄人，枣庄学院旅游与资源环境学院讲师，理学博士，主要从事气候变化研究。49枣庄学院学报被改善面积比较大，而且植被NDVI空间分布季节性规律明显12-13。同时，有研究发现植被NDVI受气候因子的影响有明显的地域差异性，其中降水影响高原北部的草原草甸等NDVI变化14，并对高原雅鲁藏布江流域NDVI变化起决定性作用15，此外，草地水分利用效率也会影响到草地NDVI的变化16。空间变化上研究得知，三江源区生长季植被NDVI整体表现为由西北向东南逐渐增加17，青海湖流域NDVI表现出同样的空间变化趋势，而且其季节性规律明显18。总的来说，基于MODIS产品研究青藏高原生长季植被NDVI的时空变化特征对加深青藏高原物候的认识及生态环境的保护有重要意义。1石研究区概况与数据方法1.1石研究区概况与数据来源青藏高原是世界海拔最高、最年轻的高原，平均海拔40 0 0 m以上，高原上湖泊众多，光照和地热资源充足。气候条件主要表现为气温低、日较差大、降水相对集中；并且都体现为由东南向西北逐渐减少的变化趋势。高原上冻土广布，主要的植被类型是草地、高山植被以及位于东南部的阔叶林和混交林等19。本研究所选用的数据为MODIS（M O D 13A 2）2 0 0 12 0 2 0 年青藏高原生长季NDVI与植被物候数据集2 0】，包括：2 0 0 12 0 2 0 年每年生长季平均NDVI，时间为每年59月份；生长季始期（SO S）、生长季末期（EOS），时间单位为d，空间分辨率为1km，数据均来自于国家青藏高原科学数据中心（https：/d a t a.t p d c.a c.c n)。1.2研究方法本研究主要采用趋势分析法和显著性检验法对青藏高原生长季NDVI的时空变化特征进行分析。（1）趋势分析法（Senslope）趋势分析法是一种通过对随时间变化的变量进行线性回归分析，从而来预测其变化趋势的方法。其计算方法如下：式中：Slope为变化趋势，i为第i年的平均值，n为研究的时间长度。当Slope0时，表示该像元处NDVI为增加趋势；当Slope=O时，表示该像元处NDVI基本保持不变；当Slope0表明NDVI序列呈上升趋势，Z0VVa(S)Z=30，S+1Var(S)n-1S=2,ZSgn(x,-x),(1)S=0(2)S 0杜宜臻，崔岩，倪杰基于MODIS数据的青藏高原生长季植被NDVI的时空变化分析式中：Sgn（x;）为复合函数，其表达式如下：Sgn(0)=(-1,0 02丝结果与分析2.1青藏高原植被生长季NDVI的分布特征本研究对2 0 0 1一2 0 2 0 年青藏高原植被生长季NDVI均值进行了统计，得出青藏高原多年平均NDVI数值为0.14（表1）。表1不同类型植被生长季NDVI的面积占比及NDVI值植被类型面积占比/%落叶阔叶林0.02针阔叶混交林0.25常绿针叶林0.14郁闭灌丛0.00栽培植被1.22草地82.22沼泽0.05稀疏灌丛0.98高山植被14.35结合高原植被类型分布图191分析发现，整个青藏高原植被生长状况相对较稀疏，同时具有空间异质性。受到气温和降水等方面的影响13，高原西北部生长季植被表现较差，腹地适中，东南部表现较好。统计结果表明，NDVI高值区一般为落叶阔叶林、针阔叶混交林、常绿针叶林、郁闭灌丛和栽培植被等，其中落叶阔叶林平均NDVI最大，为0.7 0，其次是常绿针叶林、针阔叶混交林和郁闭灌丛，分别为0.6 3、0.6 3和0.55。然而，这些高大林木的总面积不足青藏高原总面积的2%。低值区植被类型主要以高山植被、沼泽和草地为主，其中草地面积占整个青藏高原总面积的80%以上，而这些植被类型区的NDVI一般较低，多在0.2 左右及以下。基于前人研究结果进一步将青藏高原2 0 0 0 一2 0 2 0 年NDVI变化率进行如下分类，如表2 所示。表2 青藏高原生长季植被变化趋势面积占比NDVI变化趋势Slope-0.0040-0.0040Slope0.0007-0.0007 Slope 0.00250.0025Slope0.00580.0058Slope统计结果表明青藏高原生长季NDVI变化率为0.0 0 40.0 0 6，平均值接近于0。整体来看，青藏高原大部分地区的植被保持稳定生长，变化趋势基本不变，其面积占高原总面积的8 0%以上。同时发现，高原退化区面积占比相对较少，不到6%，而改善区（包含轻微改善）面积相对较高，占比达到14%左右。即在研究时段内，青藏高原的植被状况整体表现出变绿的趋势。510,0=0NDVI0.700.630.630.550.470.220.170.140.05程度退化轻微退化基本不变轻微改善改善（4)平均NDVI0.14面积占比/%0.295.1080.7312.651.23枣庄学院学报2.2青藏高原植被生长季NDVI的时间变化特征本研究进一步分析了青藏高原2 0 0 1一2 0 2 0 年NDVI的逐年变化趋势，结果如图1所示。2024年第2 期y=0.0008x-1.380.15R*=0.2570.10ION0.050.002001图12 0 0 12 0 2 0 年青藏高原植被生长季NDVI年际变化由图1可知，青藏高原NDVI在研究时段内总体表现为上升的趋势，其年际变化率为0.0 0 8（10 a）。M K 检验通过了信度9 9%的显著性检验，说明高原植被生长季NDVI在2 0 0 0 一2 0 2 0 年间上升趋势显著，且表现出一定的波动性。其中，NDVI的高值年份主要集中在2 0 0 9一2 0 11年以及2017一2 0 19 年间。此外，在2 0 2 0 年，植被NDVI出现了极端低值，为0.12，低于平均值的14.29%。数值上的差异体现了在相应时间段内高原气候以及局地环境的变化状况，这主要是受到高原暖湿变化的影响。也有相关研究证实青藏高原植被生长季NDV的变化总体上和气温和降水呈正相关关系12,15 植被生长季是指植被1a内从开始复苏返青（生长季始期）到植被枯黄落叶（生长季末期）的时段。为了更详细的探讨高原植被的生长变化，本研究引入植被生长季始期和末期两个参量，并基于这两个变量得到生长季长度（DOS，即生长季末期生长季始期）。基于此，分别对上述3个变量，包括SOS、EO S和DOS的时间变化趋势展开分析，结果如图2 所示。300250200P/E1501005002001Sos图2 2 0 0 1一2 0 2 0 年青藏高原植被生长季始期（SOS）、生长季末期（EOS）、生长季长度（DOS）变化由图2 可知，在研究时段内，青藏高原植被SOS提前，平均每年提前0.13d；而EOS在推迟，每年推迟0.0 5d，比SOS的变化速度略慢；进一步表明DOS在延长，其每年延长约0.14d。这种变522006SOS:y=0.13X=+388.372006一线性（SOS）2011年份2011EOS线性（EOS)年份20162020EOS:y=0.048 x=161.25DOS:y-0.144x=-161.0320162020DOS-一线性(DOS)杜宜臻，崔岩，倪杰基于MODIS数据的青藏高原生长季植被NDVI的时空变化分析化趋势与前人研究结果基本一致2 1-2 2。然而，在总的研究时段内，这些因子的变化也表现出一定的波动性。其中2 0 0 6 一2 0 0 9年变化相对较大，SOS较早，而EOS相对较晚，表明这段时间内，高原植被的DOS较其它时间段相对较长。进一步探究SOS、EO S和DOS三个因素与NDVI之间的关系，结果如图3所示。400TEOS:y=32.12+257.80R?=0.101300P/回日200DOS:y=33.63x+129.22R2=0.04211000sOS.图32 0 0 12 0 2 0 年青藏高原植被生长季始期（SOS）、生长季末期（EOS）、由图3可知：高原EOS与NDVI相关系数最高，其R为0.32；而SOS和DOS与NDVI的相关系数相对较低：其R分别为-0.15和0.2 1。统计结果表明，上述三个参数与NDVI的变化都具有一定的相关性，然而均表现为弱的相关性，其中 SOS与NDVI 的相关性表现为弱的负相关。综上所述，近2 0 a来，青藏高原植被的NDVI不仅在数值上表现为增大趋势，而且每年植被生长的持续时长也在增大。2.3青藏高原植被生长季NDVI与海拔的变化特征海拔升高导致气温降低，这会影响高原植物的生长和分布。基于此，本研究探究了青藏高原NDVI的变化与海拔之间的关系，结果如图4所示。1.00.80.6IAAN0.40.20.0研究结果表明（图4），整体上青藏高原生长季植被NDVI与海拔呈负相关，即NDVI会随着海拔的升高而降低，平均海拔每升高1km，ND V I 下降0.0 9。这主要是因为海拔越高，气温越低，53SOS:y=-13.48x+128.23R2=0.02150.20.4线性(SOS)EOS生长季长度（DOS）与NDVI的关系y=-9E-05 x+0.681R2=0.15830图4研究区海拔与NDVI变化趋势散点图0.6NDVI线性(EOS)2海拔/km0.8DOS线性（DOS)61.0枣庄学院学报NDVI与气温基本呈现正相关关系，因此NDVI也降低14,2 3。但是在不同海拔梯度下NDVI 随海拔的变化趋势表现出差异性，其中NDVI值在海拔高度0 1km与3 4km呈上升趋势，在1 3km与46km呈下降趋势。研究结果对青藏高原当前以及未来的生态保护等具有重要科学和现实意义。3结论本文基于2 0 0 1一2 0 2 0 年MODIS数据，探讨青藏高原地区生长季植被NDVI的时空变化特征，主要结论如下。（1）青藏高原植被生长季NDVI具有明显的空间异质性，NDVI高值区和低值区的植被类型存在较大差异。（2）研究时段内青藏高原生长季NDVI表现为上升趋势，除8 0%地区的植被保持稳定生长外，有接近14%的区域植被正逐年改善，高原整体表现出变绿的趋势。（3）整体上青藏高原生长季植被NDVI与海拔呈负相关，即NDVI随着海拔的升高而降低，平均海拔每升高1km，ND V I 下降0.0 9。参考文献1 LUO D,GUO D,JIN H,et al.Ecological impacts of degrading permafrost J.Frontiers in Earth Science,2022,10:967530.2 YUN H,TANG J,DIMPERIO L,et al.Warming and increased respiration have transformed an alpine steppe ecosystemon the Tibetan Plateau from a carbon dioxide sink into a source J.Journal of Geophysical Research:Biogeosciences,2022,127（1):e 2 0 2 1JC 0 0 6 40 6.3 YAO T,BOLCH T,CHEN D,et al.The imbalance of the Asian water tower J.Nature Reviews Earth&Environ-ment,2022,3（10):6 18-6 32.4罗京，牛富俊，林战举，等青藏高原多年冻土区热融滑塌发育特征及规律J冰川冻土，2 0 2 2，44（1）：96-105.5底阳平，张扬建，曾辉，等：“亚洲水塔”变化对青藏高原生态系统的影响J中国科学院院刊，2 0 19，34(11):1322-1331.6 郭正刚，牛富俊，湛虎，等青藏高原北部多年冻土退化过程中生态系统的变化特征J生态学报，2 0 0 7，27（8):3 2 9 4-3 3 0 1.7 JIN X,JIN H,IWAHANA G,et al.Impacts of climate-induced permafrost degradation on vegetation:a review JJ.Advances in Climate Change Research,2021,12(1):29-47.【8 杨兆平，欧阳华，宋明华，等青藏高原多年冻土区高寒植被物种多样性和地上生物量J生态学杂志，2010,29（4):6 17-6 2 3.9 J KAWAMURA K,AKIYAMA T,YOKOTA H,et al.Comparing MODIS vegetation indices with AVHRR NDVI for monito-ring the forage quantity and quality in Inner Mongolia grassland,China J.Grassland Science,2005,51(1):33-40.10 王思，张路路，林伟彪，等基于MODIS-NDVI的广东省植被覆盖与土地利用变化研究J生态学报，2022,42(6):92-106.11董璐，赵杰，刘雪佳，等.198 2 2 0 15年新疆地区植被生长对气温的响应J应用生态学报，2 0 19，30(7):2165-2170.12 欧阳习军，董晓华，魏榕，等青藏高原植被生长季NDVI时空变化及对气候因子的响应分析J水土保持研究，2 0 2 3，30（2)：1-10.13郝爱华，薛娴，段翰晨，等。青藏高原典型草地NDVI时空演变的季节差异及其气候驱动J生态学报，2023,43（1):3 52-3 6 3.14缪利，陆晴，刘根林，等.1999-2 0 19年青藏高原不同植被类型NDVI时空演变特征及其对气候因子的响应J水土保持研究，2 0 2 3，30（1)：97-10 5.542024年第2 期DU Yizhen,CUI Yan,NI Jie杜宜臻，崔岩，倪杰基于MODIS数据的青藏高原生长季植被NDVI的时空变化分析15 J SUN W,WANG Y,FUY,et al.Spatial heterogeneity of changes in vegetation growth and their driving forces based on sat-ellite observations of the Yarlung Zangbo River basin in the Tibetan Plateau J.Journal of Hydrology,2019,574:324-332.16 吉珍霞，裴婷婷，陈英，等.2 0 0 12 0 2 0 年青藏高原草地NDVI时空变化及驱动因子分析J草地学报，2022,30(7):1873-1881.17 徐嘉昕，房世波，张廷斌，等.2 0 0 0-2 0 16 年三江源区植被生长季NDVI变化及其对气候因子的响应J国土资源遥感，2 0 2 0，32（1）：2 37-2 46.【18 拓进宝.2 0 0 0-2 0 16 年青海湖流域NDVI时空变化及其与环境固子的关系D西宁：青海师范大学，2 0 19.【19 张慧，赵良，朱文泉，基于多源数据产品集成分类制作的青藏高原现状植被图J。北京师范大学学报（自然科学版），2 0 2 1，5 7（6）：8 16-8 2 4.【2 0 王泰华，杨大文，青藏高原生长季NDVI与植被物候数据集（2 0 0 12 0 2 0）【DB/OL国家青藏高原科学数据中心，2 0 2 2.21 J LIU Y,WANG J,DONG J,et al.Variations of vegetation phenology extracted from remote sensing data over the Tibet-an Plateau Hinterland during 2000-2014 J.Journal of Meteorological Research,2020,34(4):786797.22 CHENG M,JIN J,ZHANG J,et al.Effect of climate change on vegetation phenology of different land-cover types onthe Tibetan Plateau JJ.International Journal of Remote Sensing,2018,39(2):470-487.2 3鲍艳，魏宇晨，南素兰，等，全球2 温升背景下青藏高原植被对气候变化的响应J高原气象，2 0 2 3，42(1):49-59.【责任编辑：罗海南Analysis of Spatiotemporal Changes of NDVI in the Vegetation GrowingSeason in Qinghai-Tibet Plateau Based on MODIS Data(College of Tourism and Environmental Resources,Zaozhuang University,Zaozhuang 277160,China)Abstract:Based on the MODIS normalized difference vegetation index(NDVI)data from 2001 to 2020,the spatial-temporalvariation characteristics of vegetation and their relationship with the beginning of the growing season(SOS),the end of thegrowing season(EOS),and altitude.The results show that NDVI in the growing season of the QTP has obvious spatial hetero-geneity,and its value gradually decreases from southeast to northwest,with an annual mean of 0.14.In recent 20 years,thearea of stable growth and improvement reached totaled 94%,which means the ecological environment of the QTP showed agreen trend.What is more,there is a negative correlation between the NDVI and altitude,and NDVI decreases by O.09 for ev-ery 1 km increase in average altitude.The results aim to provide scientific references for future research on the restoration of theecological environment and sustainable development of the QTP.Key words:Qinghai Tibet Plateau;normalized difference vegetation index;variation characteristics;altitude55