基于GIMMS、VGT和MODIS的中国东部植被指数对比分析

GIMMS NDVI、VGT NDVI和MODIS NDVI/EVI是目前在植被变化有关研究中经常使用的植被遥感数据,它们之间的差异也得到了广泛关注。然而,在分析这些数据之间的差异时,较少有研究注意到植被本身固有的季节循环可能夸大了各数据间的相关关系。应用2000~2006年GIMMS NDVI、VGT NDVI、MODIS NDVI/EVI等不同植被遥感数据,对比了基于这些数据集的中国东部植被年际变化的差异,探讨了植被季节循环对不同遥感数据之间相关性的影响。结果表明:由不同遥感数据提取的植被年际变化特征具有明显的一致性,然而,植被本身固有的季节循环特征掩盖了不同数据集的差异。季节循环去除前,各数据集之间具有显著的相关性;季节循环去除后,各数据集的相关性明显降低,但不同数据集在北部区域依然具有较好的一致性,其差异主要出现在南部区域,差异最明显的是GIMMS与MODIS数据,二者在淮河以南的区域几乎不存在显著相关。     

北美永久湿地NDVI时序变化特征分析——以哈德逊湾西南湿地区域为例

基于2001~2009年MCD12Q1数据,提取哈德逊湾西南区连续9年都是湿地的区域作为研究区。利用1982~2006年GIMMS NDVI数据,分析此区域25 a来NDVI季节、年际变化,并与温度、降水数据进行相关分析。结果显示:研究区域NDVI季节变化呈现单峰曲线形状,夏季达到最大值。1982~2006年NDVI年际变化在春秋两季与全年平均NDVI年际变化都呈增加趋势,夏季与冬季趋势则相反。此区域NDVI季节变化与温度、降水呈极显著相关,冬季NDVI年际变化与温度以及夏季NDVI年际变化与降水的相关性不显著,而其他各时期NDVI年际变化与温度、降水也都呈极显著相关。然而,由于人类活动及其他因素的影响,此区域的NDVI年际变化驱动还存在一定的不确定性。     

1982~2015年海南岛NDVI时空变化及气候驱动力分析北大核心CSCD

为了揭示海南岛植被覆盖长期变化趋势以及进一步确定影响海南岛植被变化的主要气候驱动因子,为海南岛植被应对气候变化以及植被的良性发展提供科学指导。基于GIMMS NDVI数据,采用趋势分析法,探究1982~2015年海南岛植被的时空变化趋势;利用偏相关分析和主成分回归分析研究34 a里温度、降水和太阳辐射对海南岛植被变化的影响。研究表明:(1)空间上海南岛植被在北部和沿海地区呈明显增加趋势,而三亚及周边地区存在多处植被退化区域;(2)时间上海南岛植被整体以缓慢增长为主,增速为0.019/10 a,年际变化明显;(3)温度和太阳辐射显著主导海南岛88%地区的植被生长,其驱动因子存在如下关系:太阳辐射作用>温度作用>降水的驱动作用;(4)温度主导海南岛北部及西部地区植被的生长,太阳辐射主导驱动岛屿南部的植被生长,而降水是岛屿中部植被的主导气候驱动因子。综上所述,海南岛植被总体呈现良性发展的趋势,温度和太阳辐射是促进植被生长的主要气候因子。     

中国西北地区植被时空演变特征及其对气候变化的响应

利用GIMMS/NDVI数据分析了1982~2006年我国西北地区植被覆盖时空变化特征及其对气温和降水变化的响应。结果表明：近25 a来,中国西北地区年均植被NDVI增速为0.5%/10a,7月、8月和10月份增加趋势最显著。天山、阿尔泰山、祁连山、青海的中东部等地区植被覆盖显著增加;青海的格尔木至玉树一线、陕西的南部地区、新疆的塔里木盆地、吐鲁番、塔河、托里等地区植被退化。植被覆盖与气温、降水的年际关系都呈弱的正相关。但年内关系则都呈显著的线性关系,植被覆盖随月均温升高而增加,当月均温超过20℃时,植被NDVI呈下降趋势;月降水量在0~100 mm之间,植被NDVI随降水呈线性增长,当月降水量超过100 mm之后,不再有明显的增长趋势。     

2001~2010年松木希错流域植被动态变化遥感研究

遥感在区域植被变化研究中具有十分重要的作用,能为大面积监测植被状况的演化过程提供技术支持。NDVI在高植被覆盖地区存在过饱和现象,对稀疏地区的植被变化尤其敏感。以古里雅冰帽南部的松木希错流域植被相对稀疏区域为研究区,基于MODIS NDVI数据和逐月气象观测数据,以及RS和GIS平台,对该区域2001~2010年主要植被变化趋势进行了初步研究,并对植被变化与气候驱动因子的关系进行了分析和探讨。结果表明:① 2001~2010年间该区域的植被活动有加强趋势;② NDVI表明研究区植被生长季较短(5~9月),NDVI浮动区间为0.11～0.13,低于全国水平(0.3～0.35),也低于全球稀疏灌丛的平均水平(0.2～0.4);③NDVI与年均气温整体上呈正相关,而与年降水量相关性不强。表明近年来持续升温是影响该区域植被活动加强的最主要原因。     

基于NDVI序列影像的植被覆盖变化研究

归一化植被指数NDVI是地表植被覆盖特征的重要指标之一。以新疆石河子地区2003~2006年MODIS遥感数据反演的NDVI时间序列影像为例,分析研究了植被长势的年内和年际变化,将植被长势的年内变化和年际变化分为比前一年(月)好、比前一年(月)稍好、与前一年(月)持平、比前一年(月)稍差和比前一年(月)差5个等级,得到年内和年际间植被长势的动态分布图,从植被长势分布图中NDVI的变化可以看出年际和年内植被长势的变化。并应用变化矢量分析法对2003~2006年石河子地区NDVI的变化强度进行了分析,获得了植被覆盖变化强度分布情况,研究结果表明4 a内石河子地区植被覆盖未发生大的变化,植被系统基本稳定。     

沂蒙山区植被NDVI的时空特征及其对水热条件的响应

植被是生态环境变化的综合指示器,研究其对水热条件的响应已成为当前气候变化研究中的主要内容之一。选取北方土石山区典型代表--沂蒙山区为研究对象,基于沂蒙山区1980~2010年的气温、降水和2001~2010年MODIS\|NDVI数据,结合相关分析和最小二乘法,定量分析该区植被指数的年际、年内的时空变化及其对水热条件的响应。结果表明:①近10 a沂蒙山区NDVI max的变化斜率为0.0026;②植被显著退化区和良好区分别占研究区总面积的10.52%和28.62%;③不同季节(主要是春、夏和秋季)植被状况均呈现良性发展趋势;④台站数据显示植被年际变化与年降水和年均气温的关系并不密切,而在月时间尺度上植被与气温的相关性要强于与降水的相关性。综上所述,沂蒙山区植被状况总体呈良性发展趋势,气温可能是影响该区植被生长的主导因子。     

基于MODIS和AVHRR数据源的东北地区植被NDVI变化及其与气温和降水间的相关分析 

基于逐像元一元线性回归模型,应用MODIS NDVI数据对AVHRR-GIMMS NDVI进行时间序列拓展,拓展序列通过一致性检验,基于所建立的1982~2009年植被年最大NDVI数据集,在GIS平台上进行了植被NDVI变化和NDVI与年平均气温、年降水量之间的相关分析。研究结果表明:过去28 a间,植被年最大NDVI呈3个变化阶段:1982~1992年呈小幅上升趋势,1992~2006年呈缓慢下降趋势,2006~2009年呈缓慢回升态势。由空间变异分析得出NDVI变化相对大的区域主要分布在内蒙干旱和半干旱区。21世纪初和20世纪90年代相对于80年代NDVI值升高,3个阶段平均NDVI变化幅度为±0.3。 20世纪初,赤峰地区以及松嫩平原西部地区植被NDVI呈轻度增加的面积占全区6.45%。植被年最大NDVI与年平均气温、年降水量相关性空间差异明显。偏相关系数绝对值,气温大于降水的像元数占54%;综合分析,较降水而言,气温是东北全区植被年最大NDVI的主控影响因子。对于不同植被类型年最大NDVI,受气温影响强度由大到小依次为:森林>草地>沼泽湿地>灌丛>耕地;受降水影响按草地>耕地>灌丛>沼泽湿地>森林依次减弱。     

基于GLASS数据的青藏高原2001—2018年蒸散发时空变化分析

准确认知青藏高原蒸散发时空变化特征,为当地可持续农业的水资源规划及理解高原气候变化具有重要现实意义。研究基于GLASS陆表潜热通量产品,采用Mann-Kendall趋势分析方法,结合青藏高原生态地理分区方案,分析了2001—2018年青藏高原蒸散发的时空变化特征及其与气温、降水和植被的关系。结果表明：①GLASS ET产品可以较好地表征青藏高原蒸散发的时空分布特征;②青藏高原多年平均蒸散发为296.52 mm,整体上呈现出东南高西北低的空间格局,其中东喜马拉雅南翼最高（690.94 mm）,柴达木盆地最低（163.47 mm）;③近18 a来,青藏高原蒸散发年际变化呈波动性上升,只有东喜马拉雅南翼在下降;④研究期间,青藏高原蒸散发以显著性增长趋势为主,占47.44%,主要位于高原东部边缘和中西部腹地,呈显著性减小趋势的地区占3.82%,主要集中于东喜马拉雅南翼;⑤蒸散发的空间分布在干旱区与气温呈负相关,在湿润区呈正相关,与降水空间格局总体呈正相关;⑥蒸散发与NDVI的空间分布呈较好的正相关,与NDVI的变化趋势相关性较为复杂,大部分呈正相关,小部分呈负相关。     

基于MODIS数据的玉米植被参数估算方法的对比分析

基于实测数据建立了FPAR、LAI的植被指数估算模型（NDVI、RVI、NDWI）,并将其应用于MODIS BRDF数据对德惠地区玉米FPAR、LAI进行估算,然后将MODIS 15A2 FPAR/LAI产品值分别与BRDF估算值、地面实测值进行对比分析。主要得出以下结论：植被指数NDVI、RVI都能较好地用于实测数据和MODIS BRDF数据的FPAR、LAI估算;NDWI虽然在实测数据中估算玉米FPAR、LAI的效果优于NDVI、RVI,但其应用于MODIS BRDF数据估算FPAR、LAI时,效果却较差。BRDF数据估算FPAR与MODIS 15A2 FPAR值的关系因生长时期不同而异,在玉米生长前期,前者高于后者,而生长后期两者却较相近;BRDF估算LAI值一直都高于MODIS 15A2 LAI产品值。生长季前期,MOD15A2 FPAR、LAI值接近实测值,而在后期却高于实测值。通过分析也表明,玉米苗期MODIS 15A2 FPAR数值变化范围较小,产品算法对实际FPAR变化尚不够敏感,这可能是影响MODIS FPAR产品精度的一个原因。     

Analysis of the Spatio-temporal Variation in FPAR of the Tibetan Plateau from 1982 to 2015

The absorbed and utilized Fraction of Photosynthetically Active Radiation(FPAR) reflects the capacity of carbon fixation and oxygen release by vegetation, which may vary over space and time in large scale. Analysis of spatial-temporal variation in FPAR is an important topic of plant ecology. Based on GIMMS NDVI3g (1982~2015) and MODIS NDVI (2001~2015) data in the Tibetan plateau, here we used the nonlinear, semi-theoretical and semi-empirical models to inverse and analyze the spatial and temporal variation in FPAR. The results showed that (1) The spatial distributions of FPAR derived from GIMMS NDVI3g and MODIS NDVI were highly consistent, with the correlation coefficient being 0.82 (P<0.01). The area in which the trends of inter-annual change in the two inversion data were consistent for at least 6 years made up 80% of the studying area. (2) FPAR in Tibetan Plateau was greatly affected by slope and altitude. Changes in FPAR were fastest at slopes of 15~35 degrees and highest at altitude of 700~2 100 m. The effect of slope direction on FPAR was limited. There was little difference in FPAR among different slope directions except for the south where the FPAR was relatively lower. (3)The FPAR data from 1982 to 2015 demonstrated seasonal variation. The inter-annual variation in FPAR was most significant in winter, in which FPAR in about 78.5% of the area increased. FPAR declined most significantly in the fall. (4) FPAR derived from both the MODIS NDVI and GIMMS NDVI data demonstrated a small, temporary decline in 2012. The trend of inter-annual variation in FPAR was largely different among different vegetation types. In conclusion, the FPAR data from 1982 to 2015 in the Tibetan plateau demonstrated both spatial and seasonal variation, which may have important implications for further studies concerning climate and environmental changes in the region.     

Evaluation of earth observation based long term vegetation trends — Intercomparing NDVI time series trend analysis consistency of Sahel from AVHRR GIMMS, Terra MODIS and SPOT VGT data

AVHRR (Advanced Very High Resolution Radiometer) GIMMS (Global Inventory Modelling and Mapping Studies) NDVI (Normalized Difference vegetation Index) data is available from 1981 to present time. The global coverage 8 km resolution 15-day composite data set has been used for numerous local to global scale vegetation time series studies during recent years. Several aspects however potentially introduce noise in the NDVI data set due to the AVHRR sensor design and data processing. More recent NDVI data sets from both Terra MODIS and SPOT VGT data are considered an improvement over AVHRR and these products in theory provide a possibility to evaluate the accuracy of GIMMS NDVI time series trend analysis for the overlapping period of available data. In this study the accuracy of the GIMMS NDVI time series trend analysis is evaluated by comparison with the 1 km resolution Terra MODIS (MOD13A2) 16-day composite NDVI data, the SPOT Vegetation (VGT) 10-day composite (S10) NDVI data and in situ measurements of a test site in Dahra, Senegal. Linear least squares regression trend analysis on eight years of GIMMS annual average NDVI (2000-2007) has been compared to Terra MODIS (1 km and 8 km resampled) and SPOT VGT NDVI data 1 km (2000-2007). The three data products do not exhibit identical patterns of NDVI trends. SPOT VGT NDVI data are characterised by higher positive regression slopes over the 8-year period as compared to Terra MODIS and AVHRR GIMMS NDVI data, possibly caused by a change in channels 1 and 2 spectral response functions from SPOT VGT1 to SPOT VGT2 in 2003. Trend analysis of AVHRR GIMMS NDVI exhibits a regression slope range in better agreement with Terra MODIS NDVI for semi-arid areas. However, GIMMS NDVI shows a tendency towards higher positive regression slope values than Terra MODIS in more humid areas. Validation of the different NDVI data products against continuous in situ NDVI measurements for the period 2002-2007 in the semi-arid Senegal revealed a good agreement between in situ measurements and all satellite based NDVI products. Using Terra MODIS NDVI as a reference, it is concluded that AVHRR GIMMS coarse resolution NDVI data set is well-suited for long term vegetation studies of the Sahel-Sudanian areas receiving < 1000 mm rainfall, whereas interpretation of GIMMS NDVI trends in more humid areas of the Sudanian-Guinean zones should be done with certain reservations.     

Validation of the collection 5 MODIS FPAR product in a heterogeneous agricultural landscape in arid Uzbekistan using multitemporal RapidEye imagery

The fraction of photosynthetically active radiation (FPAR) absorbed by a vegetation canopy is an important variable for global vegetation modelling and is operationally available from data of the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensor starting from the year 2000. Product validation is ongoing and important for constant product improvement, but few studies have investigated the specific accuracy of MODIS FPAR using in situ measurements and none have focused on agricultural areas. This study therefore presents a validation of the collection 5 MODIS FPAR product in a heterogeneous agricultural landscape in western Uzbekistan. High-resolution FPAR maps were compiled via linear regression between in situ FPAR measurements and the RapidEye normalized difference vegetation index (NDVI) for the 2009 season. The data were aggregated to the MODIS scale for comparison. Data on the percentage cover of agricultural crops per MODIS pixel allowed investigation of the impact of spatial heterogeneity on MODIS FPAR accuracy. Overall, the collection 5 MODIS FPAR overestimated RapidEye FPAR between approximately 6% and 15%. MODIS quality flags, the underlying biome classification and spatial heterogeneity were investigated as potential sources of error. MODIS data quality was very good in all cases. A comparison of the MODIS land-cover product with high-resolution land-use classification revealed a significant misclassification by MODIS. Yet, we found that the overestimation of MODIS FPAR is independent of classification accuracy. The results indicate that the amount of background information, present even in the most homogeneous pixels (～70% crop cover), is most likely the reason for the overestimation. The behaviour of pure pixels could not be investigated due to a lack of appropriate pixels.     

Downscaling real-time vegetation dynamics by fusing multi-temporal MODIS and Landsat NDVI in topographically complex terrain

Canopy phenology is an important factor driving seasonal patterns of water and carbon exchange between land surface and atmosphere. Recent developments of real-time global satellite products (e.g., MODIS) provide the potential to assimilate dynamic canopy measurements with spatially distributed process-based ecohydrological models. However, global satellite products usually are provided with relatively coarse spatial resolutions, averaging out important spatial heterogeneity of both terrain and vegetation. Therefore, bias can result from lumped representation of ecological and hydrological processes especially in topographically complex terrain. Successful downscaling of canopy phenology to high spatial resolution would be indispensable for catchment-scale distributed ecohydrological modeling, aiming at understanding complex patterns of water, carbon and nutrient cycling in mountainous watersheds. Two downscaling approaches are developed in this study to overcome this issue by fusing multi-temporal MODIS and Landsat TM data in conjunction with topographic information to estimate high spatio-temporal resolution biophysical parameters over complex terrain. MODIS FPAR (fraction of absorbed photosynthetically active radiation) is used to provide medium spatial resolution phenology, while the variability of vegetation within a MODIS pixel is characterized by Landsat NDVI. The algorithms depend on the scale-invariant linear relationship between FPAR and NDVI, which is verified in this study. Downscaled vegetation dynamics are successfully validated both temporally and spatially with ground-based continuous FPAR and leaf area index measurements. Topographic correction during the downscaling process has a limited effect on downscaled FPAR products except for the period around the winter solstice in the study area.     

Atmospheric conditions for monitoring the long-term vegetation dynamics in the Amazon using normalized difference vegetation index

This study examined the effect of biomass-burning aerosols and clouds on the temporal dynamics of the normalized difference vegetation index (NDVI) exhibited by two widely used, time-series NDVI data products: the Pathfinder AVHRR land (PAL) dataset and the NASA Global Inventory Monitoring and Modeling Studies (GIMMS) dataset. The PAL data are 10-day maximum-value NDVI composites from 1982 to 1999 with corrections for Rayleigh scattering and ozone absorption. The GIMMS data are 15-day maximum-value NDVI composites from 1982 to 1999. In our analysis, monthly maximum-value NDVI was extracted from these datasets. The effects were quantified by comparing time-series of NDVI from PAL and GIMMS with observations from the SPOT/VEGETATION sensor and aerosol index data from the Total Ozone Mapping Spectrometer (TOMS), and results from radiative transfer simulation. Our analysis suggests that the substantial large-scale NDVI seasonality observed in the south and east Amazon forest region with PAL and GIMMS is primarily caused by variations in atmospheric conditions associated with biomass-burning aerosols and cloudiness. Reliable NDVI data can be typically acquired from April to July when such effects are relatively low, whereas there is a few effective NDVI data from September to December. In the central Amazon forest region, where aerosol loads are relatively low throughout the year, large-scale NDVI seasonality results primarily from seasonal variations in cloud cover. Careful treatment of these aerosol and cloud effects is required when using NDVI from PAL and GIMMS (or other source) to determine large-scale seasonal and interannual dynamics of vegetation greenness and ecosystem-atmosphere CO₂ exchange in the Amazon region.     

Spatial–temporal changes in vegetation coverage in the global coastal zone based on GIMMS NDVI3g data

ABSTRACT

In this paper, we used the Global Inventory Modelling and Mapping Studies (GIMMS) third-generation Normalized Difference Vegetation Index (NDVI) (GIMMS NDVI3g) dataset. Based on GIMMS NDVI3g data over the global coastal zone from 1982 to 2014, the spatial–temporal characteristics of vegetation coverage were analysed by plotting the spatial pattern and monthly calendar of NDVI; furthermore, historical trends and future evolutions of vegetation coverage change at the pixel scale were studied by performing the Mann-Kendall trend test and calculating the trend slope (β) and Hurst index (H) of NDVI. The main findings are as follows: 1) Vegetation density exhibits dramatic differences in the global coastal zone. Specifically, desert belts mostly have perennial non-vegetation or low vegetation coverage, and tundra belts principally have moderate or high vegetation coverage; additionally, forest belts mainly have dense vegetation coverage. 2) In the global coastal zone, intra-annual variations in vegetation coverage show a ‘∩’-shaped curve with an obvious peak from June to September (maximum in July or August), while inter-annual variations show a fluctuating but generally slowly increasing trend over the entire study period; accordingly, variations in different subregions show significant differences. 3) At monthly, seasonal and annual scales, the overall vegetation coverage increases in the global coastal zone, while there are relatively few areas with decreasing vegetation coverage; furthermore, change trends of vegetation coverage in most areas will demonstrate relatively strong positive persistence in the future. 4) The increasing trend in high-latitude coastal tundra is extremely significant in the growing season because vegetation in the tundra belts is highly sensitive to climate change. 5) Areas with a decreasing trend of vegetation coverage exhibit spatial patterns of aggregation in the ‘circum urban agglomeration’ and ‘nearby desert belt’ regions, that is, the decreasing trend of vegetation coverage is relatively high in coastal urban agglomeration areas and desert belt peripheries. This paper is expected to provide knowledge to support vegetation conservation, ecosystem management, integrated coastal zone management and climate change adaptation in coastal areas.     

Downscaling MODIS LST in the East African mountains using elevation gradient and land-cover information

Land-surface temperature (LST) is strongly affected by altitude and surface albedo. In mountain regions where steep slopes and heterogeneous land cover are predominant, LST can vary significantly within short distances. Although remote sensing currently provides opportunities for monitoring LST in inaccessible regions, the coarse resolution of some sensors may result in large uncertainties at sub-pixel scales. This study aimed to develop a simple methodology for downscaling 1 km Moderate Resolution Spectroradiometer (MODIS) LST pixels, by accounting for sub-pixel LST variation associated with altitude and land-cover spatial changes. The approach was tested in Mount Kilimanjaro, Tanzania, where changes in altitude and vegetation can take place over short distances. Daytime and night-time MODIS LST estimates were considered separately. A digital elevation model (DEM) and normalized difference vegetation index (NDVI), both at 250 m spatial resolution, were used to assess altitude and land-cover changes, respectively. Simple linear regressions and multivariate regressions were used to quantify the relationship between LST and the independent variables, altitude and NDVI. The results show that, in Kilimanjaro, altitude variation within the area covered by a 1 km MODIS LST pixel can be up to ±300 m. These altitude changes can cause sub-pixel variation of up to ±2.13°C for night-time and ±2.88°C for daytime LST. NDVI variation within 1 km pixels ranged between –0.2 and 0.2. For night-time measurements, altitude explained up to 97% of LST variation, while daytime LST was strongly affected by land cover. Using multivariate regressions, the combination of altitude and NDVI explained up to 94% of daytime LST variation in Kilimanjaro. Finally, the downscaling approach proposed in this study allowed an improved representation of the influence of landscape features on local-scale LST patterns.