浙江省植被物候变化及其对气候变化的响应

基于改进的Savitzky-Golay滤波算法重构了2001-2010年MODIS NDVI时间序列数据,反演了浙江近10 a来的植被生长状况,利用分段多项式拟合和动态阈值法提取了自然植被重要物候期(生长季开始时间、结束时间和生长季长度等),分析了植被物候期的年际变化趋势和空间分异特征,并结合同期气象数据,探讨了植被物候变化与气候变化的关系.主要结论如下:①浙江植被覆盖有所减少,整体下降趋势不显著,平均植被生长季为222 d,长度略有延长,其中开始时间提前趋势不显著,而结束时间推迟显著.②植被生长季为3月下旬至11月中旬,植被生长季的开始时间从北往南逐渐推迟,生长季结束时间相对集中,其中生长季延长的区域面积和缩短的区域面积相当.③气象要素与关键物候期参数的相关性分析表明冬季的热量供给是影响浙江植被生长开始的重要因子,植被生长季前期温度积累的增加有利于植被生长;生长季结束时间与当年温度呈极显著正相关,与当年降水和湿润指数的减少呈负相关,但影响不显著,从与各相关月、季度相关分析来看,秋季干湿程度对植被生长季结束时间影响相对较大.     

基于高分1号影像的森林植被信息提取

实时最新森林植被信息的提取是林业航空植保作业的必要前提。论文以安徽省蚌埠市为研究区域,探讨了基于高分1号卫星遥感数据在亚热带农林植被混合地区的森林植被信息提取。根据植被物候信息差异选择了提取森林植被信息的5个关键时期高分影像,采用分区决策树方法监测森林植被的空间分布和面积信息,并与未分区决策树法的提取结果进行比较。结果表明：采用分区决策树法和未分区决策树法对于大中尺度森林植被信息提取的总体精度均优于85%。但分区决策树森林植被提取总体精度达到90.72%,较未分区决策树法提高3.80%、4.65%,Kappa系数达到0.81,较未分区决策树法提高约0.07~0.10,结合植被物候信息的分区决策树森林植被提取法好于未分区决策树法,能够满足林业航空植保作业的精度需求。具有较高空间分辨率、宽覆盖、短重访周期的高分1号影像,对于大区域的林业航空植保当年最新森林植被信息的提取表现出较大的潜力。     

基于MODIS植被指数的浑善达克沙地植被物候变化

基于近20a MODIS的2种植被指数数据，利用Logistic曲线曲率极值法和动态阈值法，对浑善达克沙地植被物候进行了提取，分析其时空变化，并利用研究结果数据，分析植被物候对高程和气象因子的响应.结果表明，2000~2019年间研究区的植被物候呈微弱波动趋势，浑善达克沙地植被返青期（SOS）主要集中在110~140d，枯黄期（EOS）主要集中在250~280d，整体呈微弱推迟趋向（0.28d/a），生长季长度（LOS）主要集中在120~170d，整体呈微弱延长趋向（0.23d/a）.在空间分布上，占研究区51.51%的区域植被SOS呈提前趋势；占研究区67.02%的区域植被EOS呈提前趋势，并且占研究区32.98%的区域植被EOS呈推迟趋势.占研究区62.71%的区域植被LOS呈延长趋势.在海拔900~1500m区间，随着海拔升高，SOS显著推迟，EOS不显著提前，LOS显著缩短.前一年11、12月和当年1、2、4月降水量对SOS有提前作用.6、7、8、9月份气温和降水都对EOS推迟有明显的作用.总的来说，浑善达克沙地植被物候与气候因子的相应规律比较复杂，表现出季节性的差别及地形的差异性.     

基于无人机可见光影像的绿色植被提取

在分析健康绿色植被光谱特性及无人机可见光影像典型地物各波段像元值差异的基础上,提出一种综合利用红、绿、蓝3个可见光波段信息的新型绿色植被指数——差异增强植被指数（DEVI）.利用该指数及其他8种常见可见光植被指数结合不同阈值方法提取研究区域绿色植被信息,并采用地表真实感兴趣区和基于SVM的监督分类方法进行精度量化评价.结果表明：由DEVI计算的植被指数灰度影像直方图具有良好双峰形态,可利用双峰直方图阈值法快速确定阈值,且阈值一般位于0.9~1之间;同时,DEVI提取精度明显优于其余8种植被指数,且采用双峰直方图阈值法时,总体精度为98.98%,Kappa系数为0.9791,相对误差为1/83.为验证DEVI是否具有良好的可适用性及可靠性,选取3种典型植被覆盖区域进行可行性验证,结果表明：利用DEVI可高精度提取建筑密集区域和植被零散分布区域的绿色植被信息,总体精度分别为98.42%和98.56%,Kappa系数分别为0.9610和0.9635,相对误差分别为1/125和1/91;而植被集中分布区域提取精度略低于上述2种典型区域,总体精度为97.40%,Kappa系数为0.9371,相对误差为1/53.因此,提出的差异增强植被指数——DEVI可以有效、高精度、低成本提取不同植被覆盖典型区域无人机可见光影像中的绿色植被信息,为陆地生态系统中的绿色植被监测研究提供一种可行性方法.     

京津冀城市扩张对植被和地表城市热岛的影响

基于MODIS-LST、EVI以及土地利用/覆盖数据,分析了2000~2015年京津冀13个城市生长季及不同季节城市扩张对植被（ΔEVI）和地表城市热岛影响的时间变化趋势以及二者之间的关系,结果表明,京津冀13个城市主城区生长季及不同季节的城市扩张对植被均存在消极影响（ΔEVI<0）,对夏季植被的影响最大（ΔEVI=-0.131）;白天和夜间的地表城市热岛强度（SUHII）最大值均出现在夏季,分别为2.594℃和1.859℃;不同季节的ΔEVI均值呈现显著下降趋势,其速率分别为-0.0038,-0.0033,-0.0045和-0.0018/a（P<0.01）;生长季、春季和夏季白天SUHII均呈显著增加趋势,速率分别为0.076,0.093和0.106℃/a（P<0.01）,不同季节夜间SUHII均呈显著增加的趋势,其增加速率分别为0.05,0.055,0.049,0.054和0.046℃/a;京津冀13个城市主城区的ΔEVI均值和SUHII均值呈显著负相关（P<0.01）,其相关系数介于-0.959~-0.592之间;白天SUHII与ΔEVI的相关性最高为生长季（r=-0.959）,夜间SUHII与ΔEVI的相关性最高为春季（r=-0.936）.     

基于RS和GIS的贺兰山高山林线初步研究

论文基于RS和GIS方法,运用Landsat 8-OLI影像和DEM数据获取贺兰山高山林线分布信息,并对林线植被NDVI特征和林线分布特征进行研究。结果表明：1）生长季初期山地针叶林和高山灌丛草甸NDVI值差异明显,两种植被NDVI栅格像元的分布及其数量在山地垂直空间的变化指示出贺兰山2 900~3 200 m海拔范围出现高山林线现象。2）贺兰山高山林线主要分布在中部顶峰周边区域,林线上边界平均高度3 106 m,下边界平均高度3 022 m。主山脊两侧区域的东坡和北坡林线的平均高度较高,西坡和南坡低。3）山地针叶林与高山灌丛草甸NDVI栅格数据可以作为影像分类数据的补充,对影像分类方式确定的高山林线分布区域进行修正,能够提高林线提取的准确性。     

20世纪初以来沂蒙山区森林植被动态及其对气候变化的响应——基于遥感和树轮的研究

遥感能够直接反映植被生长信息,但目前观测时间较短;树轮资料可以弥补遥感手段的不足,但往往缺乏与植被生长状况（如植被生产力）的直接关联。已有研究往往计算不同植被指数指标和树轮指数的相关性,并基于此来重建长时间尺度的植被动态信息。但基于相关性所选取的指标在不同地区存在较大差异,不利于空间对比研究。以沂蒙山区为研究区,选取能够有效表征植被年初级生产力的生长季NDVI累积值指标,利用Bootstrap法建立了其与树轮宽度的关系,重建了20世纪初期以来植被动态时间序列,并分析了变化特征及其与气候变化的关系。结果表明：沂山地区生长季NDVI累积值的多年平均值为7.36,低于蒙山和塔山地区;蒙山地区植被状况呈现好转趋势,特别是20世纪90年代之后植被好转趋势更加明显,而塔山和沂山地区植被无显著变化。小波分析结果表明研究区植被动态存在较显著的2 a、4 a或8 a尺度的周期变化,与生长季帕尔默干旱指数及平均温度的周期变化相一致,但干旱指数与植被动态具有更高的相关性。本研究综合运用树轮与遥感技术发展了长时间尺度植被动态时间序列重建方法,有助于更好地理解植被动态及其对气候变化的响应。     

2000年以来内蒙古生长季旱情变化遥感监测及其影响因素分析

为探讨内蒙古旱情状况及其影响因素,利用MODIS 16 d合成的植被指数产品数据MOD13A2和8 d合成的地表温度产品数据MOD11A2构建Ts-NDVI特征空间,计算温度植被干旱指数（TVDI）,基于内蒙古2000—2017年生长季每16 d的TVDI分析了近18年来内蒙古生长季旱情时空分布特征,结合气温和降水资料初步探讨了旱情变化的影响因素。结果表明：（1）2000—2017年内蒙古生长季TVDI平均值为0.6,重旱和中旱所占面积最大,其中2007年、2010年为旱情最为严重的年份。内蒙古干旱空间分异明显,西南部以轻旱为主,中部地区以中旱为主,大兴安岭以西的呼伦贝尔草原等地旱情严重。（2）近18年内蒙古生长季干旱程度呈现轻微加重趋势,年际变化值θ_slope介于-0.07~0.7之间,阿鲁科尔沁旗东北部至霍林河一带旱情加重趋势最为严重,阿荣旗和扎兰屯等农业生产地区旱情有轻微加重趋势。（3）2017年内蒙古生长季以区域性和局域性干旱为主,6月和9月干旱最为严重,呼伦贝尔草原和鄂尔多斯高原西部干旱发生频率高且程度重。（4）内蒙古干旱影响因子分析结果表明,TVDI值与气温呈正相关、与降水和坡度呈负相关、与小于1300 m的高程呈正相关、与大于1300 m的高程呈负相关关系。内蒙古生长季TVDI与气温和降水偏相关分析结果表明,锡林郭勒盟苏尼特左旗北部、呼伦贝尔鄂伦春自治旗和呼伦贝尔草原等地旱情与气温正相关关系较为显著（P<0.01）,锡林郭勒盟东北部干旱情况与降水负相关关系较为显著（P<0.01）,其中,气温对旱情的影响强于降水。     

基于MODIS数据估算海河流域植被生态用水方法探讨

基于中分辨率成像光谱仪(MODIS)数据和地面气象数据,建立了区域植被生态用水模型．通过地表温度和植被覆盖度的“现状”三角关系,应用该模型计算了地表温度．植被覆盖度指数(ITVC),结合植被系数法计算了地表植被生态用水,并用该模型的计算结果分析了海河流域的生态用水．结果显示,基于“现状”三角法和植被系数法计算的植被现状环境用水与地表大型蒸渗仪实测结果较一致．海河流域植被生态用水具有明显的区域分布特征,西部山区和中东部黑龙港及运东平原和大清河淀东平原区植被生态用水量比较低,而在太行山山前平原、东南部平原以及北部燕山一带植被生态用水量比较高．     

荒漠草原区油气管道建设对地表植被格局的影响

以油气管道经过的乌鲁木齐荒漠草原地区为研究对象,选择施工前、施工期以及恢复期3个不同施工阶段的遥感影像数据,解译地表植被类型,采用象元二分法提取植被覆盖度,通过比较时间上植被覆盖格局和土地利用格局的变化,分析油气管道建设对生态系统影响的范围和特点.研究发现:不同建设时期管道两边荒漠草原区土地利用类型差异较大,草地面积持续降低,共计减少32.53%,灌木林地面积增加到24.61%,建设用地面积在建设期增加,植被恢复期降低.研究区复杂程度提高,原有生态结构发生改变;建设期管道施工作业带及两侧10m范围内植被覆盖度显著低于其他区域,而管道两侧50~100以及100~200m范围内的植被覆盖度无明显差异,但与200~300m范围相比差异显著.管道建设主要影响范围为管道上方施工作业带以及两侧10m的范围内,间接影响范围为施工作业带两侧200m内;恢复期人为措施对研究区植被覆盖度的主要影响范围为施工作业带;管道建设改变了荒漠草原生态系统的组成结构,是管道周边土地利用类型和荒漠草地系统变化的主要推动力;人为恢复措施无法使受影响区域的植被恢复到建设前的水平.     

退耕还林（草）工程实施前后黄土高原植被覆盖时空变化分析

基于GIMMS NDVI3g（the third generation of Global Inventory Modeling and Mapping Studies Normalized Difference Vegetation Index）数据,辅以趋势分析、Mann-Kendall检验、Hurst指数等方法,识别了1982—2013年及1982—1999、2000—2013年黄土高原植被覆盖时空演变特征,并探讨其驱动因素。研究发现：1）1982—2013年及1982—1999、2000—2013年期间黄土高原生长季NDVI分别以0.019/10 a（P<0.01）、0.016/10 a（P<0.05）和0.057/10 a（P<0.001）的速率增加;2）除1999年以前林地外,所有植被类型的生长季NDVI均呈现显著的增加趋势,2000—2013年尤为明显;3）黄土高原生长季NDVI呈现由东南向西北递减的趋势,1982—2013年及1982—1999、2000—2013年NDVI显著上升的面积分别占74.94%、24.26%和53.34%,主要集中在黄土高原的北部和中部地区;4）研究区未来生长季NDVI呈持续性和反持续的比重分别为33.32%和66.68%,其中持续改善和由改善变为退化的面积分别占31.08%和61.88%;5）2000年以后降水增多与生长季NDVI上升相对应,大规模的生态工程建设对2000—2013年生长季NDVI增加有重要影响。     

柴达木盆地植被生长时空变化特征及其对气候要素的响应

以2001-2010 年MODIS NDVI 植被数据为基础,并结合该区同期气温和降水量数据,采用线性趋势分析、偏相关分析、Hurst 指数等数理分析方法,研究了近10 a 来柴达木盆地植被时空变化特征、影响因素及未来可能的变化趋势。结果表明：①2001-2010 年柴达木盆地植被生长季NDVI平均值（NDVI）呈显著增加趋势,线性增长率为0.041/10 a。柴达木盆地主要植被类型灌丛、草原、荒漠NDVI 的线性增长率分别为0.043/10 a、0.034/10 a、0.028/10 a;②柴达木盆地植被变化具有阶段性特征,整体呈“S”型增长,具有两次明显的植被快速增长期;③柴达木盆地植被改善面积占研究区总面积的12.43%,并主要分布在布尔汉布达山、祁漫塔格山、鄂拉山、柴达木山、宗务隆山的高山区,冬给措纳湖周边和绿洲核心区。植被退化面积仅占研究区总面积的0.35%,并主要分布在绿洲边缘区,特别是柴达木河绿洲和那陵格勒河绿洲的中下游地区;④植被NDVI 与1-5 月平均气温以及5-8 月可利用降水量呈显著正相关关系,植被生长与温度呈显著正相关的区域主要分布在绿洲核心区及高山区,占植被区面积的8.36%,植被生长与可利用降水量呈显著正相关的区域主要分布在低山区及山地河谷地带,占植被区面积的30.95%;⑤植被生长季延长和生长加速是柴达木盆地植被NDVI 增加的主要原因,气候的暖湿化是促使柴达木盆地植被改善的主要驱动力;⑥柴达木盆地植被改善具有强持续性,未来大部分区域植被将持续改善。     

塔里木河下游区域植被时空变化

为了对比生态输水工程实施前后塔里木河下游植被变化,分析植被对生态输水的时空响应,基于1986-2017年研究区上方无云的22期Landsat NDVI数据,辅以波动分析、趋势分析,探索研究区植被的时空变化特征,并对其驱动因素进行分析。结果表明：（1）空间上,植被的分布与水资源的分布密切相关,主要沿河道和积水区域分布。（2）1986-2000年,NDVI呈现缓慢下降趋势,2000-2017年NDVI总体呈现波动增加趋势;植被对生态输水有一定的滞后性,且与输水量密切相关。（3）研究区NDVI时间序列的方差介于0.0077~0.1957间,存在明显的空间差异,其波动的大小与植被种类密切相关。（4）1986-2000年,植被的NDVI空间趋势以减小为主,并主要发生在英苏和老英苏附近;2000-2017年,植被的NDVI空间趋势以增加为主,且主要发生在沿河道区域和喀尔达依东侧。     

神东矿区AVHRR/NDVI的时空、开采强度和气候效应

研究采矿扰动区长时序、多维度NDVI变化趋势,有助于认识非自然生态区植被演替规律,对高强度开采条件下西北脆弱生态区地表植被自然修复和人工修复工作具有指导意义.研究通过设立直接影响区、间接影响区和自然生态校验区,利用长时序GIMMS AVHRR/NDVI（1981～2006）植被指数数据集,结合同期温度、降水和煤炭产量信息,从时间、空间、气候和开采强度方面开展对比分析,并以长时序MODIS NPP/NDVI（2000～2010）数据进行结果验证.时序分析表明,全球变化背景下,神东矿区植被生长季二度延长;空间分析表明,神东矿区NDVI增量低于缓冲区,缓冲区NDVI增量低于自然生态区;开采强度分析表明,随着神东矿区开采强度增强,NDVI增速放缓,增长速率低于自然生态区;气候变化分析表明,神东矿区NDVI受温度升高、降水减少共同作用,与温度相关性较高,与降水相关性较低;同期MODIS NPP/NDVI验证结果表现出与AVHRR/NDVI一致的规律性.     

2000~2020年黄河流域植被时空演化驱动机制

以归一化植被指数（NDVI）作为植被覆盖及生长状况指标,基于2000～2020年MODIS NDVI数据及同时期气象数据,采用Theil-Sen斜率估算、Mann-Kendall检验、相关性分析和残差分析等方法研究了2000～2020黄河流域植被时空演化驱动机制.结果表明,2000～2020年黄河流域生长季NDVI均值以0.005 a^-1的速率波动上升,植被明显改善的区域主要分布于流域中游的秦岭山系、陕北高原和吕梁山系;黄河流域生长季NDVI与降水和气温的偏相关系数均值分别为0.57和0.49,降水对植被的影响高于气温;人类活动对植被生长起明显改善的区域主要分布在流域中部的陕北高原、吕梁山系和宁夏南部等区域,对植被生长起抑制作用的区域主要分布在银川、包头、西安、洛阳、郑州和太原等人类活动强烈的城市区域;人类活动和气候变化分别对黄河流域植被变化贡献了72%和28%,在人类活动和气候变化的驱动下,黄河流域植被生长得到改善的面积占流域面积的96.4%,其中人类活动贡献率大于80%的区域面积占34.3%,主要分布在流域中部和东南部.气候变化贡献率大于80%的区域面积占4....     

2000-2013年塔里木河流域生长季NDVI时空变化特征及其影响因素分析

受塔里木河流域综合治理工程实施和近期气候变化的影响,流域植被覆被时空分布产生一定的变化,厘清植被覆被与流域气候变化及人类活动的关系可以为塔里木河流域生态维护与治理提供科学参考依据。为此,论文以NDVI为指示因子,运用趋势分析、R/S分析、偏相关分析以及残差分析等方法,分析了2000-2013年综合治理工程期间NDVI的时空变化特征,并探讨及区分降水、气温气候因子以及人类活动对植被覆被变化的影响范围和程度,结果表明：1）2000-2013年,塔里木河流域植被生长季NDVI总体呈现增加趋势,增加速率为0.8%/10 a,平原区增速明显高于山区;且开都河-孔雀河流域山区、塔里木河干流的上、中游部分地区呈现比较明显的退化趋势。与此同时,塔里木河干流下游生长季NDVI持续改善。2）山区植被覆被变化主要受气候变化的影响,其中温度是高山区植被生长的主要限制因子,温度的增加促进植被的生长;中低山区以及出山口平原地区植被生长季NDVI变化是降水和温度共同作用的结果,且主要受降水的影响。降水与植被生长季NDVI变化呈正相关,温度与植被生长季NDVI变化呈负相关。3）平原绿洲区植被生长季NDVI增加主要是绿洲灌区面积不断增加以及塔里木河流域生态治理工程对植被恢复的结果,人类活动是该区域植被生长的主要驱动力。4）塔里木河干流生态闸工程在恢复下游植被的同时,也在一定程度上影响了上、中游地区的用水,尤其是导致中游植被出现退化趋势,退化速率约为0.1%/10 a。相关部门应进一步加强水资源的合理配置,充分发挥生态闸工程的水资源调度调控作用。     

北部湾沿海地区植被覆盖对气温和降水的旬响应特征

论文分析北部湾沿海地区植被覆盖对气候变化的响应,为该区域植被恢复和植被生产力研究提供依据.基于研究区2000--2011年423景SPOT-VEGETATION逐旬NDVI数据及逐日气温和降水数据,利用像元二分模型,相关分析,偏相关分析和时滞偏相关分析等数理统计方法,对研究区植被覆盖时空变化特征及与旬降水和旬均温的相关性及滞后性进行分析.结果表明:1)近12 a来,北部湾沿海地区植被覆盖度平均值呈增长趋势,由2000年的65.23%增加到2011年的72.02%,增加了6.79%;2)研究区植被生长季旬NDVI均值介于0.21~0.67之间,在不同时期变化是不同的,其值呈现出"降低--增长--降低"3种变化过程;3)各种植被类型与温度具有显著的相关关系,全部通过了0.01的显著性水平,且NDVI与温度的显著性水平高于NDVI与降水的显著性水平,说明北部湾沿海地区植被覆盖NDVI对气象因子中的温度更为敏感;4)NDVI与气温的时滞偏相关系数显著高于NDVI与降水的时滞偏相关系数,旬NDVI与旬降水的滞后时间多集中于6~9旬之间,而旬NDVI与旬气温的滞后时间多以0~5旬为主;5)不同类型植被的生长对气温和降水的响应时间不一致,但与水热条件时滞偏相关系数越高的植被类型响应时间越短.上述研究结果表明,近12 a来北部湾沿海地区植被处于恢复状态,且植被对降水和气温具有明显的阈值和滞后性.     

1982-2013年潞安矿区NDVI3g变化趋势及气候响应

为了揭示全球变化背景下矿业开发活动对矿区生态环境的影响,采用长时序GIMMS AVHRR NDVI3g（1982-2013年）全球植被指数数据集和气候信息（年均降水量和气温数据）,运用IDL编程实现数据合成运算、线性回归和趋势拟合,从时间、空间、气候三方面对矿区、缓冲区（10 km、20 km）、校验区的NDVI平均值和总值进行比较研究,并推算出生长期变化趋势.时序分析表明,32 a来矿区NDVI总量随开采年限延长呈先增后减/波动中下降的趋势,下降速率为0.18/（10 a）.矿区植被返青期滞后3 d,枯黄期提前30 d,生长期缩短33 d,缩减速率为10.3 d/（10 a）.空间分析表明,除潞安矿区外其他三区生长季均有所延长,研究区平均NDVI年增长率依次为矿区（1.09%） < 10 km缓冲区（2.16%） < 20 km缓冲区（8.86%） < 校验区（9.87%）,矿区NDVI总量自1995年后开始减少,非开采区NDVI总量呈增加趋势.气候变化分析表明,校验区NDVI对温度敏感性高于降水量,矿区NDVI对降水量敏感性高于温度,其中温度对两区植被生长有明显滞后性.研究显示,矿业开发活动抑制了矿区及周边区域NDVI的增长,NDVI年增长率远低于校验区,受温度升高、降水量减少共同作用,自然生态下校验区NDVI呈增加趋势,生长季延长;而受开采扰动影响下的矿区植被活动呈减弱趋势,生长期也有所缩短.