Investigation of terrain illumination effects on vegetation indices and VI-derived phenological metrics in subtropical deciduous forests

We used RapidEye and Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra data to study terrain illumination effects on 3 vegetation indices (VIs) and 11 phenological metrics over seasonal deciduous forests in southern Brazil. We applied TIMESAT for the analysis of the Enhanced Vegetation Index (EVI) and the Normalized Difference Vegetation Index (NDVI) derived from the MOD13Q1 product to calculate phenological metrics. We related the VIs with the cosine of the incidence angle i (Cos i) and inspected percentage changes in VIs before and after topographic C-correction. The results showed that the EVI was more sensitive to seasonal changes in canopy biophysical attributes than the NDVI and Red-Edge NDVI, as indicated by analysis of non-topographically corrected RapidEye images from the summer and winter. On the other hand, the EVI was more sensitive to terrain illumination, presenting higher correlation coefficients with Cos i that decreased with reduction in the canopy background L factor. After C-correction, the RapidEye Red-Edge NDVI, NDVI, and EVI decreased 2%, 1%, and 13% over sunlit surfaces and increased up to 5%, 14%, and 89% over shaded surfaces, respectively. The EVI-related phenological metrics were also much more affected by topographic effects than the NDVI-derived metrics. From the set of 11 metrics, the 2 that described the period of lower photosynthetic activity and seasonal VI amplitude presented the largest correlation coefficients with Cos i. The results showed that terrain illumination is a factor of spectral variability in the seasonal analysis of phenological metrics, especially for VIs that are not spectrally normalized.     

运用遥感指数的崩滑检测方法

以汶川地震时出现的崩滑为研究对象,以该地区震前/后TM影像为研究数据,利用变化检测的方法,比较和分析了归一化植被指数(NDVI)、归一化建筑指数(NDBI)和归一化燃烧率(NBR)应用在崩滑检测中的可行性和优劣性。结果表明,3种遥感指数均可以用于崩滑检测,但效果不同:NDVI对崩滑检测最全面但可信度最低,NBR对崩滑检测最可信但遗漏最多,NDBI处于中间水平。     

Evaluating and comparing Sentinel 2A and Landsat-8 Operational Land Imager (OLI) spectral indices for estimating fire severity in a Mediterranean pine ecosystem of Greece

The main purpose of this study is to explore the relationship between three field-based fire severity indices (Composite Burn Index-CBI, Geometrically structure CBI, weighted CBI) and spectral indices derived from Sentinel 2A and Landsat-8 OLI imagery on a recent large fire in Thasos, Greece. We employed remotely sensed indices previously used from the remote sensing fire community (Normalized Difference Vegetation Index (NDVI), Normalized Burn Ratio (NBR), differenced NDVI, differenced NBR, relative differenced NBR, Relativized Burn Ratio) and seven Sentinel 2A-specific indices considering the availability of spectral information recorded in the red-edge spectral region. The statistical correlation indicated a slightly stronger relationship between the differenced NBR and the GeoCBI for both Sentinel 2A (r = 0.872) and Landsat-8 OLI (r = 0.845) imagery. Predictive local thresholds of dNBR values showed slightly higher classification accuracy for Sentinel 2A (73.33%) than Landsat-8 OLI (71.11%), suggesting the adequacy of Sentinel 2A for forest fire severity assessment and mapping in Mediterranean pine ecosystems. The evaluation of the classification thresholds calculated in this study over other fires with similar pre-fire conditions could contribute in the operational mapping and reconstruction of the historical patterns of fire severity over the Eastern Mediterranean region.     

Developing soil indices based on brightness,darkness, and greenness to improve land surface mapping accuracy

Soil, as one of the three basic biophysical components, has been understudied using remote sensing techniques compared to vegetation and impervious surface areas (ISA). This study characterized land surfaces based on the brightness–darkness–greenness model. These three dimensions, brightness, darkness, and greenness, were represented by the first Tasseled Cap Transformation (TC1), Normalize Difference Snow Index (NDSI), and Normalized Difference Vegetation Index (NDVI), respectively. The Ratio Index for Bright Soil (RIBS) was developed based on TC1 and NDSI, and the Product Index for Dark Soil (PIDS) was established by TC1 and NDVI. Their applications to the Landsat 8 Operational Land Imager images and 500 m 8-day composite Moderate Resolution Imaging Spectroradiometer (MODIS) in China revealed the efficiency. The two soil indices proficiently highlighted soil covers with consistently the smallest values, due to larger TC1 and smaller NDSI values in bright soil, and smaller NDVI and TC1 values in dark soil. The RIBS is capable of distinguishing bright soil from ISA without masking vegetation and water body. The spectral separability bright soil and ISA were perfect, with a Jeffries–Matusita distance of 1.916. And the PIDS was the only soil index that could discriminate dark soil from other land covers including ISA. The soil areas in China were classified using a simple threshold method based on MODIS images. An overall accuracy of 94.00% was obtained, with the kappa index of 0.8789. This study provided valuable insights into developing indices for characterizing land surfaces from different perspectives.     

Soil Moisture Estimation in a Vegetated Area Using Combination of AIRSAR and Landsat5-TM Images

High difference between dielectric constant of water (dielectric constant about 80) and dielectric constant of dried soil (dielectric constant about 2–3) makes Synthetic Aperture Radar (SAR) highly capable in soil moisture estimation. However, there are other factors which affect on radar backscattering coefficient. The most important parameters are vegetation cover, surface roughness and sensor parameters (frequency, polarization and incidence angle). In this paper, the importance of considering the effects of these parameters on SAR backscatter coefficients is shown by comparing different soil moisture estimation models. Moreover, an experimental soil moisture estimation model is developed. It is shown that this model can be used to estimate soil moisture under a variety of vegetation cover densities. The new developed model is based on combination of different indices derived from Landsat5-Thematic Mapper and AIRSAR images. The AIRSAR image is used for extraction of backscattering coefficient and incidence angle while TM image is used for calculation of Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Normalized Difference Water Index (NDWI) and Brightness Temperature. Then a soil moisture estimation model which is named as Hybrid model is developed based on integration of all of these parameters. The accuracies of this model are assessed in the NDVI ranges of 0–0.2, 0.2–0.4 and 0.4–0.7 by using SAR data in C band and L band frequencies and also in different polarizations of HH, HV, VV and TP. The results show that for instance in L band with HV polarization, R-square values of 0.728, 0.628 and 0.527 are obtained between ground measured soil moisture and estimated soil moisture values using the Hybrid model for NDVI ranges of 0–0.2, 0.2–0.4 and 0.4–0.7, respectively.     

Reconstructing satellite images to quantify spatially explicit land surface change caused by fires and succession: A demonstration in the Yukon River Basin of interior Alaska

Land surface change caused by fires and succession is confounded by many site-specific factors and requires further study. The objective of this study was to reveal the spatially explicit land surface change by minimizing the confounding factors of weather variability, seasonal offset, topography, land cover, and drainage. In a pilot study of the Yukon River Basin of interior Alaska, we retrieved Normalized Difference Vegetation Index (NDVI), albedo, and land surface temperature (LST) from a postfire Landsat image acquired on August 5th, 2004. With a Landsat reference image acquired on June 26th, 1986, we reconstructed NDVI, albedo, and LST of 1987–2004 fire scars for August 5th, 2004, assuming that these fires had not occurred. The difference between actual postfire and assuming-no-fire scenarios depicted the fires and succession impact. Our results demonstrated the following: (1) NDVI showed an immediate decrease after burning but gradually recovered to prefire levels in the following years, in which burn severity might play an important role during this process; (2) Albedo showed an immediate decrease after burning but then recovered and became higher than prefire levels; and (3) Most fires caused surface warming, but cooler surfaces did exist; time-since-fire affected the prefire and postfire LST difference but no absolute trend could be found. Our approach provided spatially explicit land surface change rather than average condition, enabling a better understanding of fires and succession impact on ecological consequences at the pixel level.     

基于植被指数模型的淡水湖泊湿地景观信息提取

高邮湖湿地是江苏省重要湿地之一,对生态、环境控制、调节气候和保护生物多样性具有重要意义。采用2007年的LandsatTM影像作为遥感信息源,选择影像的光谱特征和比值植被指数(RVI)、差值植被指数(DVI)、归一化植被指数(NDVI)、归一化差异绿度指数(NDGI)、土壤调节植被指数(SAVI)和最佳土壤调节植被指数(OSAVI)6种植被指数做了光谱特征分析,从而确定出最佳指数模型,并基于决策树方法,实现研究区景观信息的遥感分类。研究结果表明,决策树分类法易于综合多种特征进行遥感影像分类,植被指数参与到决策树分类中能够提高分类的总体精度,其总体精度达到79.58%,Kappa系数为0.721 0,分类结果理想且人工参与灵活。     

Evaluation of forest fire on Madeira Island using Sentinel-2A MSI imagery

A forest fire started on August 8th, 2016 in several places on Madeira Island causing damage and casualties. As of August 10th the local media had reported the death of three people, over 200 people injured, over 950 habitants evacuated, and 50 houses damaged. This study presents the preliminary results of the assessment of several spectral indices to evaluate the burn severity of Madeira fires during August 2016. These spectral indices were calculated using the new European satellite Sentinel-2A launched in June 2015. The study confirmed the advantages of several spectral indices such as Normalized Difference Vegetation Index (NDVI), Green Normalized Difference Vegetation Index (GNDVI), Normalized Burn Ratio (NBR) and Normalized Difference Vegetation Index (NDVIreXn) using red-edge spectral bands to assess the post-fire conditions. Results showed high correlation between NDVI, GNDVI, NBR and NDVIre1n spectral indices and the analysis performed by Copernicus Emergency Management Service (EMSR175), considered as the reference truth. Regarding the red-edge spectral indices, the NDVIre1n (using band B5, 705 nm) presented better results compared with B6 (740 nm) and B7 (783 nm) bands. These preliminary results allow us to assume that Sentinel-2 will be a valuable tool for post-fire monitoring. In the future, the two twin Sentinel-2 satellites will offer global coverage of the Madeira Archipelago every five days, therefore allowing the simultaneous study of the evolution of the burnt area and reforestation information with high spatial (up to 10 m) and temporal resolution (5 days).     

Relative suitability of indices derived from Landsat ETM+ and SPOT 5 for detecting fire severity in sagebrush steppe

Remote sensing indices of burn area and fire severity have been developed and tested for forest ecosystems, but not sparsely vegetated, desert shrub-steppe in which large wildfires are a common occurrence and a major issue for land management. We compared the performance of remote sensing indices for detecting burn area and fire severity with extensive ground-based cover assessments made before and after the prescribed burning of a 3 km² shrub-steppe area. The remote sensing indices were based on either Landsat 7 ETM+ or SPOT 5 data, using either single or multiple dates of imagery. The indices delineating burned versus unburned areas had better overall, User, and Producer's accuracies than indices delineating levels of fire severity. The Soil Adjusted Vegetation Index (SAVI) calculated from SPOT had the greatest overall accuracy (100%) in delineating burned versus unburned areas. The relative differenced Normalized Burn Ratio (RdNBR) using Landsat ETM+ provided the highest accuracies (73% overall accuracy) for delineating fire severity. Though SPOT's spatial resolution likely conferred advantages for determining burn boundaries, the higher spectral resolution (particularly band 7, 2.21 μm) of Landsat ETM+ may be necessary for detecting differences in fire severity in sparsely vegetated shrub-steppe.     

Illumination effects on the differenced Normalized Burn Ratio's optimality for assessing fire severity

The influence of illumination effects on the optimality of the dNBR (differenced Normalized Burn Ratio) was evaluated for the case of the 2007 Peloponnese (Greece) wildfires using a pre/post-fire Landsat TM (Thematic Mapper) image couple. Well-illuminated pixels (south and south-east facing slopes) exhibited more optimal displacements in the bi-spectral feature space than more shaded pixels (north and north-west exposed slopes). Moreover, pixels experiencing a small image-to-image difference in illumination obtained a higher optimality than pixels with a relatively large difference in illumination. To correct for illumination effects, the c-correction method and a modified c-correction technique were applied. The resulting median dNBR optimality of uncorrected, c-corrected and modified c-correction data was respectively 0.58, 0.60 and 0.71 (differences significant for p < 0.001). The original c-correction method improved the optimality of badly illuminated pixels while deteriorating the optimality of well-illuminated pixels. In contrast, the modified c-correction technique improved the optimality of all the pixels while retaining the prime characteristic of topographic correction techniques, i.e. detrending the illumination–reflectance relationship. For a minority of the data, for shaded pixels and/or pixels with a high image-to-image difference in illumination, the original c-correction outperformed the modified c-correction technique. In this study conducted in rugged terrain and with a bi-temporal image acquisition scheme that deviated up to two months from the ideal anniversary date scheme the modified c-correction technique resulted in a more reliable change detection.     

Detection of subtle deforestation due to logging using satellite remote sensing in wet and dry savanna woodlands of Southern Africa

Tropical deforestation through logging activities poses a direct threat to biodiversity. However, the detection of logging has remained a challenge. Based on study sites in Zimbabwe and Zambia, we tested whether the Normalized Difference Vegetation Index (NDVI) and the Coefficient of Variation in NDVI (CVNDVI) derived from high and medium spatial resolution satellite data could be used to detect logging in dry and wet miombo woodlands. Separately, we integrated NDVI and CVNDVI in logistic regression to test whether each can be used to successfully predict logging in the study sites. We tested whether the spatial resolution of satellite data has an effect in detection of logging using NDVI and CVNDVI derived from Landsat 8 and Worldview-2. Based on the ROC curves, we concluded that remotely sensed data could provide an effective predictive tool for detecting logging. However, in wet miombo woodlands the predictive power of remotely sensed data is weak.     

Fuzzy Based Approach for Moist Deciduous Forest Identification Using MODIS Temporal Data

Fuzzy based soft classification have been used immensely for handling the mixed pixel and hence to extract the single class of interest. The present research attempts to extract the moist deciduous forest from MODIS temporal data using the Possibilistic c-Means (PCM) soft classification approach. Temporal MODIS (7 dates) data were used to identify moist deciduous forest and temporal AWiFS (7 dates) data were used as reference data for testing. The Simple Ratio (SR), Normalized Difference Vegetation Index (NDVI), Soil Adjusted Vegetation Index (SAVI), and Transformed Normalized Difference Vegetation Index (TNDVI) were used to generate the temporal vegetation indices for both the MODIS and the AWiFS datasets. It was observed from the research that the MODIS temporal NDVI data set1, which contain the minimum number of images and avoids the temporal images corresponding to the highest frequency stages of onset of greenness (OG) and end of senescence (ES) activity of moist deciduous forest have been found most suitable data set for identification of moist deciduous forest with the maximum fuzzy overall accuracy of 96.731 %.     

Soil salinity and vegetation cover change detection from multi-temporal remotely sensed imagery in Al Hassa Oasis in Saudi Arabia

Detecting soil salinity changes and its impact on vegetation cover are necessary to understand the relationships between these changes in vegetation cover. This study aims to determine the changes in soil salinity and vegetation cover in Al Hassa Oasis over the past 28 years and investigates whether the salinity change causing the change in vegetation cover. Landsat time series data of years 1985, 2000 and 2013 were used to generate Normalized Difference Vegetation Index (NDVI) and Soil Salinity Index (SI) images, which were then used in image differencing to identify vegetation and salinity change/no-change for two periods. Soil salinity during 2000–2013 exhibits much higher increase compared to 1985–2000, while the vegetation cover declined to 6.31% for the same period. Additionally, highly significant (p < 0.0001) negative relationships found between the NDVI and SI differencing images, confirmed the potential long-term linkage between the changes in soil salinity and vegetation cover.     

Mapping rubber trees based on phenological analysis of Landsat time series data-sets

This study proposes a strategy for accurate mapping of rubber trees through the analysis of Landsat time series datasets. The phenological dynamics of rubber trees were derived from the Normalized Difference Vegetation Index (NDVI) to verify the three important phenological metrics of rubber trees; defoliation, foliation and their growing stages. A decade (2006–2015) ago, Landsat time series NDVIs were used to study the strength of relationship between rubber trees, evergreen trees and oil palm trees. Two important results that could discriminate these three types of vegetation were found; firstly, a weak relationship of NDVIs between rubber trees and evergreen trees during the defoliation period (r² = 0.1358) and secondly between rubber trees and oil palm trees during the growing period (r² = 0.2029). This analysis was verified using Support Vector Machine to map the distribution of the three types of vegetation. An accurate mapping strategy of rubber trees was successfully formulated.     

黑河流域叶面积指数的遥感估算

研究利用Landsat7ETM⁺遥感数据获取黑河流域植被叶面积指数(LAI)空间分布的可行性。该研究是基于黑河流域分布式水文模型的一个重要输入项———LAI空间分布数据的需要而产生的。文章在详尽的野外观测数据基础上,分别探究实测LAI与同时相ETM⁺3、4、5、7波段反射率及相关植被指数(SR、NDVI、ARVI、RSR、SAV I、PVI、GESAVI)的相关关系,率定最佳的LAI遥感反演及其空间分布方案。研究发现,针对特定的自然条件,将研究区分为植被覆盖度小的稀疏立地和覆盖度大的密集立地,分别采用土壤调节植被指数(SAVI)和大气阻抗植被指数(ARVI)进行2种林地的LAI估算最为可靠,在此基础上,提出黑河地区LAI估算及其空间分布的遥感制图方案。     

Comparison of the Landsat Surface Reflectance Climate Data Record (CDR) and manually atmospherically corrected data in a semi-arid European study area

This study contributes to the quality assessment of atmospherically corrected Landsat surface reflectance data that are routinely generated by the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS). This dataset, named Landsat Surface Reflectance Climate Data Record (Landsat CDR), is available at global scale and offers unprecedented opportunities to land monitoring and management services that require atmospherically corrected Earth observation (EO) data. Our assessment is based on the comparison of the Landsat CDR data against a set of Landsat and DEIMOS-1 images processed to a high degree of accuracy using an industry-standard atmospheric correction algorithm (ATCOR-2). The software package has been used for many years and its correction procedures can be considered consolidated and well-established. The dataset of Landsat and DEIMOS-1 images was acquired over a semi-arid agricultural area located in Lower Austria and was independently corrected by using a manual fine-tuning of ATCOR-2 parameters to reach the highest possible accuracy. Results show a very good correspondence of the surface reflectance in each of the six reflective spectral channels as well as for the NDVI (Normalized Difference Vegetation Index). An additional comparison against a NDVI time series from MODIS revealed also a good correspondence. Coefficients of determination (R²) between the two multi-year and multi-seasonal Landsat/DEIMOS datasets range between 0.91 (blue band) and 0.98 (nIR, SWIR-1 and SWIR-2). The results obtained for our semi-arid test site in Austria confirm previous findings and suggest that automatic atmospheric procedures, such as the one implemented by LEDAPS are accurate enough to be used in land monitoring services that require consistent multi-temporal surface reflectance data.     

Wetland Monitoring,Serving as an Index of Land Use Change-A Study in Samaspur Wetlands,Uttar Pradesh,India

Wetlands are among the most productive ecosystems in the world and any alterations might lead to changes in their bio-physical, socio-economic and climatic conditions. Wetland dynamics as an index of land use change were studied. Satellite remote sensing was utilized to understand the periodic and seasonal dynamics of Samaspur wetlands using Landsat and RESOURCESAT-1 temporal data. Index-based (i.e., Normalized Difference Water Index (NDWI) and Normalized Difference Vegetation Index (NDVI)) classification resulted in meaningful discrimination of wetland classes. Results indicate (i) effective water spread areas have increased to optimum capacity at 1990 due to the influence of Sharda canal, (ii) expansion of the agricultural area has led to reduction of the wetland buffer area, and (iii) increase in vegetation biomass due to pesticide-fertilizer runoff and sedimentation load. We also reiterate (i) free availability of the Landsat satellite data in public domain facilitating such monitoring studies and (ii) availability and utility of SWIR band information in wetland classification exercise. The study concludes that policy-driven measures have both long and short term impacts on land use and its natural wetland ecosystems; and the characterizing the later serves as indictor of the former and perhaps vice versa.     

Vegetation recovery monitoring over the waste dump in Haizhou opencast coalmine area,China, during 1975–2000 using NDVI and VF index

Vegetation condition monitoring has been done from 1975 to 2000 in the waste dump of Haizhou opencast coalmine area, China, using remote sensing techniques with the objective of improving our understanding of the temporal and spatial variation of vegetation recovery in the mining dump. Four historical vegetation indexes (NDVI, VF, soil brightness and vegetation greenness) from two Landsat 2 MSS images and two Landsat 5 TM images are extracted and analyzed. For the purpose of comparison and analysis two improved techniques such as normalization grading of change slope and image segmentation were used in this study. Based on the results obtained through the above analysis two conclusions are derived: (1) vegetation recovery in the study area is in an improved condition, (2) two remote sensing based vegetation indexes such as VF and NDVI are the optimal parameters to monitor vegetation condition, which could be used as the indicators of land reclamation progress in the mining area.     

Characterizing spring phenology of temperate broadleaf forests using Landsat and Sentinel-2 time series

Vegetation phenology has a great impact on land-atmosphere interactions like carbon cycling, albedo, and water and energy exchanges. To understand and predict these critical land-atmosphere feedbacks, it is crucial to measure and quantify phenological responses to climate variability, and ultimately climate change. Coarse-resolution sensors such as MODIS and AVHRR have been useful to study vegetation phenology from regional to global scales. These sensors are, however, not capable of discerning phenological variation at moderate spatial scales. By offering increased observation density and higher spatial resolution, the combination of Landsat and Sentinel-2 time series might provide the opportunity to overcome this limitation.In this study, we analyzed the potential of combined Sentinel-2 and Landsat time series for estimating start of season (SOS) of broadleaf forests across Germany for the year 2018. We tested two common statistical modeling approaches (logistic and generalized additive models using thin plate splines) and the two most commonly used vegetation indices, the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI).We found strong agreement between SOS estimates from logistic and spline models (r_EVI = 0.86; r_NDVI = 0.65), whereas agreement was higher for EVI than for NDVI (RMSD_EVI = 3.07, RMSD_NDVI = 5.26 days). The choice of vegetation index thus had a higher impact on the results than the fitting method. The EVI-based SOS also showed higher correlation with ground observations compared to NDVI (r_EVI = 0.51, r_NDVI = 0.42). Data density played an important role in estimating land surface phenology. Models combining Sentinel-2A/B, with an average cloud-free observation frequency of 12 days, were largely consistent with the combined Landsat and Sentinel-2 models, suggesting that Sentinel-2A/B may be sufficient to capture SOS for most areas in Germany in 2018. However, in non-overlapping swath areas and mountain areas, observation frequency was significantly lower, underlining the need to combine Landsat and Sentinel-2 for consistent SOS estimates over large areas. Our study demonstrates that estimating SOS of temperate broadleaf forests at medium spatial resolution has become feasible with combined Landsat and Sentinel-2 time series.     

Imaging and visualizing the spectral signatures from Landsat TM and ‘τ’ value-based surface soil microzonation mapping at and around Agartala (India)

This work aims to assess the soil microzonation of Agartala city and its surrounding areas based on spectral geophysical signatures. Different spectral resolutions of Landsat TM have been used for assessing the Normalized Difference Vegetative Index, spatial thermal emission representation and plant water moisture representation. Normalized Difference Vegetative Index (NDVI) was measured from band 4 (near-infrared (NIR)) and band 3 (photosynthetically active radiation (PAR)). The Digital Number (DN) values of thermal infrared band (TIR) were used for measuring spatial variation of thermal representation in the city area. A very simple model was developed for measuring thermal emission representative index from NDVI and automated classified TIR band. Overlaid NDVI and classified TIR shows the spatial distribution of thermal emission representative values. Classified mid-wave infrared band (MWIR) was used for measuring the surface geotherm units (τ) which are related with different types of soil. On the basis of spatial distribution of τ value which is clearly visible in a thermal emission representative map overlaid by classified MWIR, the soil microzonation map of the study area was prepared. This soil microzonation map shows that Agartala and its surrounding areas are characterized by four types of soil which are related to different geomorphic and geological units. The soil of this area is classified as dry sandy soil and sandy clay soil of the highland areas and wet sandy alluvium and clayey alluvium of the flood plain area.