Investigating the impact of the temporal resolution of MODIS data on measured phenology in the prairie grasslands

ABSTRACT

The temporal resolution of vegetation indices (VIs) determines the details of seasonal variation in vegetation dynamics observed by remote sensing, but little has been known about how the temporal resolution of VIs affects the retrieval of land surface phenology (LSP) of grasslands. This study evaluated the impact of temporal resolution of MODIS NDVI, EVI, and per-pixel green chromatic coordinate (GCCpp) on the quality and accuracy of the estimated LSP metrics of prairie grasslands. The near-surface PheonoCam phenology data for grasslands centered over Lethbridge PhenoCam grassland site were used as the validation datasets due to the lack of in situ observations for grasslands in the Prairie Ecozone. MODIS Nadir Bidirectional Reflectance Distribution Function (BRDF)-Adjusted Reflectance (NBAR) data from 2001 to 2017 were used to compute the time series of daily reference and to simulate 2–32 day MODIS VIs. The daily reference and simulated multi-day time series were fitted with the double logistic model, and the LSP metrics were then retrieved from the modeled daily time series separately. Comparison within satellite-based estimates showed no significant difference in the phenological metrics derived from daily reference and multi-day VIs resampled at a time step less than 18 days. Moreover, a significant decline in the ability of multi-day VIs to predict detailed temporal dynamics of daily reference VIs was revealed as the temporal resolution increased. Besides, there were a variety of trends for the onset of phenological transitions as the temporal resolution of VIs changed from 1 to 32 days. Comparison with PhenoCam phenology data presented small and insignificant differences in the mean bias error (MBE) and the mean absolute error (MAE) of grassland phenological metrics derived from daily, 8-, 10-, 14-, and 16-day MODIS VIs. Overall, this study suggested that the MODIS VIs resampled at a time step less than 18 days are favorable for the detection of grassland phenological transitions and detailed seasonal dynamics in the Prairie Ecozone.     

Seasonal trend analysis (STA) of MODIS vegetation index time series for the mangrove canopy of the Teacapan-Agua Brava lagoon system,Mexico

Monthly time series, from 2001 to 2016, of the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) from MOD13Q1 products were analyzed with Seasonal Trend Analysis (STA), assessing seasonal and long-term changes in the mangrove canopy of the Teacapan-Agua Brava lagoon system, the largest mangrove ecosystem in the Mexican Pacific coast. Profiles from both vegetation indices described similar phenological trends, but the EVI was more sensitive in detecting intra-annual changes. We identified a seasonal cycle dominated by Laguncularia racemosa and Rhizophora mangle mixed patches, with the more closed canopy occurring in the early autumn, and the maximum opening in the dry season. Mangrove patches dominated by Avicennia germinans displayed seasonal peaks in the winter. Curves fitted for the seasonal vegetation indices were better correlated with accumulated precipitation and solar radiation among the assessed climate variables (Pearson’s correlation coefficients, estimated for most of the variables, were r ≥ 0.58 p < 0.0001), driving seasonality for tidal basins with mangroves dominated by L. racemosa and R. mangle. For tidal basins dominated by A. germinans, the maximum and minimum temperatures and monthly precipitation fit better seasonally with the vegetation indices (r ≥ 0.58, p < 0.0001). Significant mangrove canopy reductions were identified in all the analyzed tidal basins (z values for the Mann-Kendall test ≤ ?1.96), but positive change trends were recorded in four of the basins, while most of the mangrove canopy (approximately 87%) displayed only seasonal canopy changes or canopy recovery (z > ?1.96). The most resilient mangrove forests were distributed in tidal basins dominated by L. racemosa and R. mangle (Mann-Kendal Tau t ≥ 0.4, p ≤ 0.03), while basins dominated by A. germinans showed the most evidence of disturbance.     

Spectral anisotropy of subtropical deciduous forest using MISR and MODIS data acquired under large seasonal variation in solar zenith angle

Recent studies in Amazonian tropical evergreen forests using the Multi-angle Imaging SpectroRadiometer (MISR) and the Moderate Resolution Imaging Spectroradiometer (MODIS) have highlighted the importance of considering the view-illumination geometry in satellite data analysis. However, contrary to the observed for evergreen forests, bidirectional effects have not been evaluated in Brazilian subtropical deciduous forests. In this study, we used MISR data to characterize the reflectance and vegetation index anisotropies in subtropical deciduous forest from south Brazil under large seasonal solar zenith angle (SZA) variation and decreasing leaf area index (LAI) from the summer to winter. MODIS data were used to observe seasonal changes in the normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI). Topographic effects on their determination were inspected by dividing data from the summer to winter and projecting results over a digital elevation model (DEM). By using the PROSAIL, we investigated the relative contribution of LAI and SZA to vegetation indices (VI) of deciduous forest. We also simulated and compared the MISR NDVI and EVI response of subtropical deciduous and tropical evergreen forests as a function of the large seasonal SZA amplitude of 33°. Results showed that the MODIS-MISR NDVI and EVI presented higher values in the summer and lower ones in the winter with decreasing LAI and increasing SZA or greater amounts of canopy shadows viewed by the sensors. In the winter, NDVI reduced local topographic effects due to the red-near infrared (NIR) band normalization. However, the contrary was observed for the three-band EVI that enhanced local variations in shaded and sunlit surfaces due to its strong dependence on the NIR band response. The reflectance anisotropy of the MISR bands increased from the summer to winter and was stronger in the backscattering direction at large view zenith angles (VZA). EVI was much more anisotropic than NDVI and the anisotropy increased from the summer to winter. It also increased from the forward scatter to the backscattering direction with the predominance of sunlit canopy components viewed by MISR, especially at large VZA. Modeling PROSAIL results confirmed the stronger anisotropy of EVI than NDVI for the subtropical deciduous and tropical evergreen forests. PROSAIL showed that LAI and SZA are coupled factors to decrease seasonally the VIs of deciduous forest with the first one having greater importance than the latter. However, PROSAIL seasonal variations in VIs were much smaller than those observed with MODIS data probably because the effects of shadows in heterogeneous canopy structures or/and cast by emergent trees and from local topography were not modeled.     

View-illumination effects on hyperspectral vegetation indices in the Amazonian tropical forest

Because of the pointing capability of the Hyperion/Earth Observing-One (EO-1) to improve the revisit time of the scene, temporal series of narrowband vegetation indices (VIs) can be generated to study the phenology of the Amazonian tropical forests. In this study, 10 selected narrowband VIs calculated from Hyperion nadir and off-nadir data and from different view directions (forward scattering and backscattering) were analyzed for their sensitivity to view-illumination effects along the dry season on the Seasonal Semi-deciduous Forest. Data analysis was also supported by PROSAIL modeling to simulate the spectral response of this forest type in both directions. Hyperion and PROSAIL results showed that the Enhanced Vegetation Index (EVI) and Photochemical Reflectance Index (PRI) were the two more anisotropic VIs, whereas the Normalized Difference Vegetation Index (NDVI), Structure Insensitive Pigment Index (SIPI) and the Vogelmann Red Edge Index (VOG) were comparatively less sensitive to view-illumination effects. When compared to the other VIs and because of the greater dependence on the near-infrared (NIR) reflectance, EVI showed a different spectral behavior. EVI increased from forward scattering to backscattering and with decreasing solar zenith angle (SZA) towards the end of the local dry season, due to reduction in shading and enhancement of the illumination effects. On the other hand, PRI was higher with increasing shading in the forward scattering direction, as deduced from the PROSAIL simulation. Results emphasized the importance of taking into account bidirectional effects when analyzing temporal series of VIs collected over tropical forests by imaging spectrometers with pointing capability or even by multispectral sensors with large field-of-view (FOV).     

Local variability in the timing and intensity of tropical dry forest deciduousness is explained by differences in forest stand age

Tropical Dry Forest deciduousness is a behavioral response to climate conditions that determines ecosystem-level carbon uptake, energy flux, and habitat conditions. It is regulated by factors related to stand age, and landscape scale variability in deciduous phenology may affect ecosystem functioning in forests throughout the tropics. This study determines whether observed phenological differences are explainable by forest age in the southern Yucatán Peninsula in Mexico, where forest clearing for shifting cultivation has created a mosaic of forest stands of varying age. Matched-pair statistical tests compare neighboring forest pixels of different age class (12–22 years versus 22+ years) and detect significant differences in Moderate Resolution Imaging Spectroradiometer (MODIS) Enhanced Vegetation Index (EVI)-derived metrics related to the timing and intensity of deciduousness during three dry seasons (2008–2011). In all seasons, young forests exhibit significantly more intense deciduousness, measured as total seasonal change of EVI normalized by annual maximum EVI (p < 0.001), and larger normalized EVI change during successive dry season months relative to start-of-dry-season EVI (p < 0.001), than neighboring older forests subject to similar environmental conditions.     

Evaluating uncertainty in Landsat-derived postfire recovery metrics due to terrain,soil, and shrub type variations in southern California

ABSTRACT

Temporal trajectories of apparent vegetation abundance based on the multi-decadal Landsat image series provide valuable information on the postfire recovery of chaparral shrublands, which tend to mature within one decade. Signals of change in fractional shrub cover (FSC) extracted from time-sequential Normalized Difference Vegetation Index (NDVI) data can be systematically biased due to spatial variation in shrub type, soil substrate, or illumination differences associated with topography. We evaluate the effects of these variables in Landsat-derived metrics of FSC and postfire recovery, based upon three chaparral sites in southern California which contain shrub community ecotones, complex terrain, and soil variations. Detailed validations of prefire and postfire FSC are based on high spatial resolution ortho-imagery; cross-stratified random sampling is used for variable control. We find that differences in the composition and structure of shrubs (inferred from ortho-imagery) can substantially influence FSC-NDVI relations and impact recovery metrics. Differences in soil type have a moderate effect on the FSC-NDVI relation in one of the study sites, while no substantial effects were observed due to variation of terrain illumination among the study sites. Arithmetic difference recovery metrics – based on NDVI values that were not normalized with unburned control plots – correlate in a moderate but significant manner with a change in FSC (R ² values range 0.47–0.59 at two sites). Similar regression coefficients resulted from using Landsat visible reflectance data alone. The lowest correlations to FSC resulted from Soil-Adjusted Vegetation Index (SAVI) and are attributed to the effects of the soil-adjustment factor in sparsely vegetated areas. The Normalized Burn Ratio and Normalized Burn Ratio 2 showed a moderate correlation to FSC. This study confirms the utility of Landsat NDVI data for postfire recovery evaluation and implies a need for stratified analysis of postfire recovery in some chaparral landscapes.     

Linking satellite-based spring phenology to temperate deciduous broadleaf forest photosynthesis activity

Satellite-based remote sensed phenology has been widely used to assess global climate change. However, it is constrained by uncertain linkages with photosynthesis activity. Two dynamic threshold methods were employed to retrieve spring phenology metrics from four Moderate Resolution Imaging Spectroradiometer (MODIS) products, including fraction of Absorbed Photosynthetically Active Radiation (fAPAR), Leaf Area Index (LAI), Normalized Difference Vegetation Index (NDVI), and Enhanced Vegetation Index (EVI) for three temperate deciduous broadleaf forests in North America between 2001 and 2009. These MODIS-based spring phenology metrics were subsequently linked to the photosynthetic curves (daily gross primary productivity, GPP) measured by an eddy covariance flux tower. The 20% dynamic threshold spring onset metrics from MODIS products were closer to the photosynthesis onset metrics at the date of 2% GPP increase for NDVI and fAPAR, and closer to the date of 5% and 10% increase of GPP for EVI and LAI, respectively. The 50% dynamic threshold onset metrics were closer to the photosynthesis onset metrics at the date of 10% GPP increase for NDVI, and closer to the date of 20% GPP increase for fAPAR, LAI and EVI, respectively. These results can improve our knowledge on the photosynthesis activity status of remotely sensed spring phenology metrics.     

不同传感器的模拟植被指数对水稻叶面积指数的估测精度和敏感性分析

本研究利用水稻冠层高光谱数据,模拟NOAA-AVHRR,Terra-MODIS和Landsat-TM的可见光波段反射率数据,计算各传感器的多种植被指数(NDVI,RVI,EVI,GNDVI,GRVI和Red-edge RVI),比较植被指数模型对水稻LAI的估测精度,分析不同植被指数对LAI变化的敏感性.相对于红波段植被指数,红边比值植被指数(Red-edge RVI)和绿波段指数GRVI与LAI有更好的线性相关关系,而GNDVI和LAI呈现更好的对数相关关系.MODIS的Red-edge RVI指数不仅模型拟合的精度最高,还有独立数据验证的估测精度也最高,而且它的验证精度较拟合精度下降幅度最小;其次是绿波段构建的GNDVI和GRVI植被指数的估测精度,再次是NDVI和EVI的估测精度,而RVI的估测精度最差.敏感性分析发现,13个植被指数对水稻LAI的估测能力都随着LAI的增加而下降,但归一化类植被指数和比值类植被指数对LAI变化反应的差异明显,归一化类植被指数在LAI较低时(LAI<1.5)对LAI变化的反应开始非常敏感,但迅速下降,而比值类植被指数在LAI较低时,明显小于归一化类植被指数,之后随着LAI的增大(LAI>1.5)比值类植被指数对LAI的变化敏感性,则明显高于归一化类植被指数.Red-edge RVI和绿波段指数GRVI和LAI不仅表现了很好的线性相关关系,而且在LAI大于2.9左右保持较高的敏感性.     

Evaluating the capabilities of Sentinel-2 for quantitative estimation of biophysical variables in vegetation

The red edge position (REP) in the vegetation spectral reflectance is a surrogate measure of vegetation chlorophyll content, and hence can be used to monitor the health and function of vegetation. The Multi-Spectral Instrument (MSI) aboard the future ESA Sentinel-2 (S-2) satellite will provide the opportunity for estimation of the REP at much higher spatial resolution (20 m) than has been previously possible with spaceborne sensors such as Medium Resolution Imaging Spectrometer (MERIS) aboard ENVISAT. This study aims to evaluate the potential of S-2 MSI sensor for estimation of canopy chlorophyll content, leaf area index (LAI) and leaf chlorophyll concentration (LCC) using data from multiple field campaigns. Included in the assessed field campaigns are results from SEN3Exp in Barrax, Spain composed of 35 elementary sampling units (ESUs) of LCC and LAI which have been assessed for correlation with simulated MSI data using a CASI airborne imaging spectrometer. Analysis also presents results from SicilyS2EVAL, a campaign consisting of 25 ESUs in Sicily, Italy supported by a simultaneous Specim Aisa-Eagle data acquisition. In addition, these results were compared to outputs from the PROSAIL model for similar values of biophysical variables in the ESUs. The paper in turn assessed the scope of S-2 for retrieval of biophysical variables using these combined datasets through investigating the performance of the relevant Vegetation Indices (VIs) as well as presenting the novel Inverted Red-Edge Chlorophyll Index (IRECI) and Sentinel-2 Red-Edge Position (S2REP). Results indicated significant relationships between both canopy chlorophyll content and LAI for simulated MSI data using IRECI or the Normalised Difference Vegetation Index (NDVI) while S2REP and the MERIS Terrestrial Chlorophyll Index (MTCI) were found to have the strongest correlation for retrieval of LCC.     

PROSAIL模型模拟下的植被指数土壤调节能力比较与适用环境分析

为了消除土壤背景信息对植被指数的影响,近几十年发展了土壤调节植被指数系列(SAVI family)。在不同环境条件下,不同指数抗土壤影响的能力不同。在总结了以消除土壤影响为目的的植被指数建立过程的基础上,利用PROSAIL辐射传输模型模拟的两组数据集,比较分析了NDVI、SAVI、TSAVI、MSAVI、OSAVI和GESAVI在不同叶面积指数(LAI)对应不同土壤背景的情况下抗土壤干扰、表达植被信息的能力,指出了不同植被指数应用的最适环境条件。结合植被指数—信噪比图,将这6种植被指数分成3类:在中低LAI值下,若植被覆盖度均匀,OSAVI和TSAVI有较强的消除土壤影响、表达植被信息的能力;当区域LAI分布不均、植被类型混杂时,MSAVI在表达植被信息时具有较好的稳定性。根据每类植被指数的特征,利用MODIS-VI和MODIS-LAI产品初步验证了上述结论的有效性。     

Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops

Leaf area index (LAI) and biomass are important indicators of crop development and the availability of this information during the growing season can support farmer decision making processes. This study demonstrates the applicability of RapidEye multi-spectral data for estimation of LAI and biomass of two crop types (corn and soybean) with different canopy structure, leaf structure and photosynthetic pathways. The advantages of Rapid Eye in terms of increased temporal resolution (∼daily), high spatial resolution (∼5 m) and enhanced spectral information (includes red-edge band) are explored as an individual sensor and as part of a multi-sensor constellation. Seven vegetation indices based on combinations of reflectance in green, red, red-edge and near infrared bands were derived from RapidEye imagery between 2011 and 2013. LAI and biomass data were collected during the same period for calibration and validation of the relationships between vegetation indices and LAI and dry above-ground biomass. Most indices showed sensitivity to LAI from emergence to 8 m²/m². The normalized difference vegetation index (NDVI), the red-edge NDVI and the green NDVI were insensitive to crop type and had coefficients of variations (CV) ranging between 19 and 27%; and coefficients of determination ranging between 86 and 88%. The NDVI performed best for the estimation of dry leaf biomass (CV = 27% and r² = 090) and was also insensitive to crop type. The red-edge indices did not show any significant improvement in LAI and biomass estimation over traditional multispectral indices. Cumulative vegetation indices showed strong performance for estimation of total dry above-ground biomass, especially for corn (CV ≤ 20%). This study demonstrated that continuous crop LAI monitoring over time and space at the field level can be achieved using a combination of RapidEye, Landsat and SPOT data and sensor-dependant best-fit functions. This approach eliminates/reduces the need for reflectance resampling, VIs inter-calibration and spatial resampling.     

Exploring the use of vegetation indices to sense canopy nitrogen to phosphorous ratio in grasses

Reduced availability of plant nutrients such as nitrogen (N) and phosphorous (P) has detrimental effects on plant growth. Plant N:P ratio, calculated as the quotient of N and P concentrations, is an ecological indicator of relative N and P limitation. Remote sensing has already been widely used to detect plant traits in foliage, particularly canopy N and P concentrations and could be used to detect canopy N:P faster and at lower cost than traditional destructive methods. Despite the potential opportunity of applying remote sensing techniques to detect canopy N:P, studies investigating canopy N:P remote detection are scarce. In this study, we examined if vegetation indices developed for canopy N or P detection can also be used for canopy N:P detection. Using in situ spectrometry, we measured the reflectance of a common grass species, Yorkshire fog (Holcus lanatus L.), grown under different nutrient ratios and levels. We calculated 60 VIs found in literature and compared them to optimized VIs developed specifically for this study. The VIs were calculated using both the original narrow band spectra and the spectra resampled to the band properties of six satellite sensors (MSI – Sentinel 2, OLCI – Sentinel 3, MODIS – Terra/Aqua, OLI – Landsat 8, WorldView 4 and RapidEye) to investigate the influence of bandwidths and band positions. The results showed that canopy N:P was significantly related to both existing VIs (r² = 0.16 - 0.48) and optimized VIs (r² = 0.59 – 0.72) with correlations similar to what was observed for canopy N or canopy P. Existing VIs calculated with MSI and OLI sensors bands showed higher correlation with canopy N:P compared to the other sensors while the correlation with optimized VIs was not affected by the differences in sensors’ bands. This study might lead to future practical applications using in situ reflectance measurements to sense canopy N:P in grasslands.     

Detecting winter wheat phenology with SPOT-VEGETATION data in the North China Plain

Monitoring phenological change in agricultural land improves our understanding of the adaptation of crops to a warmer climate. Winter wheat–maize and winter wheat–cotton double-cropping are practised in most agricultural areas in the North China Plain. A curve-fitting method is presented to derive winter wheat phenology from SPOT-VEGETATION S10 normalized difference vegetation index (NDVI) data products. The method uses a double-Gaussian model to extract two phenological metrics, the start of season (SOS) and the time of maximum NDVI (MAXT). The results are compared with phenological records at local agrometeorological stations. The SOS and MAXT have close agreement with in situ observations of the jointing date and milk-in-kernel date respectively. The phenological metrics detected show spatial variations that are consistent with known phenological characteristics. This study indicates that time-series analysis with satellite data could be an effective tool for monitoring the phenology of crops and its spatial distribution in a large agricultural region.     

Retrieval of EVI from Oceansat 2 Data and Comparison with MODIS Derived EVI

The vegetation index is derived using many remote sensing sensors. Vegetation Index is extensively used and remote sensing has become the primary data source. Number of vegetation indices (VIs) have been developed during the past decades in order to assess the state of vegetation qualitatively and quantitatively. Analysis of vegetation indices has been carried out by many investigators scaling from regional level to global level using the remote sensing data of varying spatial, temporal and radiometric resolutions. There are as many as 14 VIs in use. Globally operational algorithms for generation of NDVI have utilized digital counts, at sensor radiances, ‘normalized’ reflectance (top of the atmosphere), and more recently, partially atmospheric corrected (ozone absorption and molecular scattering) reflectance. Presently NDVI and EVI are standard MODIS data products which are widely used by the scientific community for environmental studies. The OCM sensor in Oceansat 2 is designed for ocean colour studies. The OCM sensor has been used for studying ocean phytoplankton, suspended sediments and aerosol optical depth by many investigators. In addition to its capability of studying the ocean surface, OCM sensor has also the potential to study the land surface features. In a past EVI has been retrieved using OCM sensor of Oceansat 1. However, there is slight change in the band width of Oceansat 2—OCM sensor compared with OCM of Oceansat 1 sensor. In the present paper an attempt has been made to derive EVI using Oceansat 2 OCM sensor and the results have been compared with MODIS data. The enhanced vegetation index (EVI) is calculated using the reflectance values obtained after removing molecular scattering and ozone absorption component from the total radiance detected by the sensor. The band-2, Band-3, band-6 and band-8 corresponding to Blue, Red and Infrared part of the visible spectrum have been used to determine EVI. The result shows that Oceansat 2 derived EVI and MODIS derived EVI are well correlated.     

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.     

A hyperspectral experiment over tropical forests based on the EO-1 orbit change and PROSAIL simulation

ABSTRACT

We designed a unique hyperspectral experiment from the Earth Observing One (EO-1) orbit change to evaluate solar illumination effects over tropical forests in Brazil. Ten nadir-viewing Hyperion images collected over a fixed site and period of the year (July to August) were selected for analysis. We evaluated variations in reflectance and in 16 narrowband vegetation indices (VIs) with increasing solar zenith angle (SZA) from the pre-drift (2004–2008) to the EO-1 drift period (2011–2016). To detect changes in reflectance and shadows, we applied spectral mixture analysis (SMA) and principal component analysis (PCA) and calculated the similarity spectral angle (θ) between the vegetation spectra measured with variable SZA. The magnitude of the illumination effects was also evaluated from change-point analysis and nonparametric Mann-Whitney U tests applied over the time series. Finally, we complemented our experiment using the PROSAIL model to simulate the VIs variation with increasing SZA resultant from satellite drift. The results showed significant changes in Hyperion reflectance and VIs, especially when the EO-1 crossed the study area at earlier times and larger SZA in 2015 (9:05 a.m.; SZA = 59°) and 2016 (8:30 a.m.; SZA = 67°). Compared to the pre-drift period (10:30 a.m.; SZA = 45°), the SZA differences of 14° (2015) and 22° (2016) increased the shade fractions and decreased the vegetation brightness. PCA separated the pre-drift and drift reflectance datasets, showing shifts in scores due to changes in brightness. θ increased with SZA, indicating changes in the shape of the vegetation spectra with drift. For most VIs, the change-point analysis indicated 2015 (SZA = 59°) as the predominant year of detected changes. Compared to the EO-1 original orbit, the Plant Senescence Reflectance Index (PSRI), Anthocyanin Reflectance Index (ARI) and Structure Insensitive Pigment Index (SIPI) presented the largest positive changes during drift, while the Photochemical Reflectance Index (PRI), Visible Atmospherically Resistant Index (VARI) and Enhanced Vegetation Index (EVI) had the largest negative changes. The effect size of the illumination geometry on these VIs was large, as indicated by increasing values of the Cohen’s r metric toward 2016. The anisotropy of the Hyperion VIs was generally consistent with that from PROSAIL in the simulated pre-drift and drift periods. Focusing on structural indices, it affected the relationships between VIs and simulated leaf area index (LAI) at large SZA.     

MODIS增强型植被指数EVI与NDVI初步比较

利用东亚地区典型地带性植被和MODIS数据，对广泛使用的植被指数NDVI和新开发的增强型植被指数EVI进行了对比分析。由MODIS开发的NDVI和EVI对干旱-半湿润环境下低覆盖植被的描述能力相似，但对湿润环境下高密度植被的描述有明显差别：NDVI年时间过程的季节性不明显，表现为全年高平的曲线；而EVI仍然有季节性，表现为钟形曲线，与月平均温度关系更密切。EVI的这一特征为研究高覆盖植被的季节性变化提供了新的思路。     

Temporal optimisation of image acquisition for land cover classification with Random Forest and MODIS time-series

The analysis and classification of land cover is one of the principal applications in terrestrial remote sensing. Due to the seasonal variability of different vegetation types and land surface characteristics, the ability to discriminate land cover types changes over time. Multi-temporal classification can help to improve the classification accuracies, but different constraints, such as financial restrictions or atmospheric conditions, may impede their application. The optimisation of image acquisition timing and frequencies can help to increase the effectiveness of the classification process. For this purpose, the Feature Importance (FI) measure of the state-of-the art machine learning method Random Forest was used to determine the optimal image acquisition periods for a general (Grassland, Forest, Water, Settlement, Peatland) and Grassland specific (Improved Grassland, Semi-Improved Grassland) land cover classification in central Ireland based on a 9-year time-series of MODIS Terra 16 day composite data (MOD13Q1). Feature Importances for each acquisition period of the Enhanced Vegetation Index (EVI) and Normalised Difference Vegetation Index (NDVI) were calculated for both classification scenarios. In the general land cover classification, the months December and January showed the highest, and July and August the lowest separability for both VIs over the entire nine-year period. This temporal separability was reflected in the classification accuracies, where the optimal choice of image dates outperformed the worst image date by 13% using NDVI and 5% using EVI on a mono-temporal analysis. With the addition of the next best image periods to the data input the classification accuracies converged quickly to their limit at around 8–10 images. The binary classification schemes, using two classes only, showed a stronger seasonal dependency with a higher intra-annual, but lower inter-annual variation. Nonetheless anomalous weather conditions, such as the cold winter of 2009/2010 can alter the temporal separability pattern significantly. Due to the extensive use of the NDVI for land cover discrimination, the findings of this study should be transferrable to data from other optical sensors with a higher spatial resolution. However, the high impact of outliers from the general climatic pattern highlights the limitation of spatial transferability to locations with different climatic and land cover conditions. The use of high-temporal, moderate resolution data such as MODIS in conjunction with machine-learning techniques proved to be a good base for the prediction of image acquisition timing for optimal land cover classification results.     

Isolating type-specific phenologies through spectral unmixing of satellite time series

Vegetation phenology is commonly studied using time series of multi-spectral vegetation indices derived from satellite imagery. Differences in reflectance among land-cover and/or plant functional types are obscured by sub-pixel mixing, and so phenological analyses have typically sought to maximize the compositional purity of input satellite data by increasing spatial resolution. We present an alternative method to mitigate this ‘mixed-pixel problem’ and extract the phenological behavior of individual land-cover types inferentially, by inverting the linear mixture model traditionally used for sub-pixel land-cover mapping. Parameterized using genetic algorithms, the method takes advantage of the discriminating capacity of calibrated surface reflectance measurements in red, near infrared, and short-wave infrared wavelengths, as well as the Normalized Difference Vegetation Index (NDVI) and the Normalized Difference Water Index. In simulation, the unmixing procedure reproduced the reflectances and phenological signals of grass, crop, and deciduous forests with high fidelity (RMSE?相似文献   

Spatially detailed retrievals of spring phenology from single-season high-resolution image time series

Vegetation indices derived from satellite image time series have been extensively used to estimate the timing of phenological events like season onset. Medium spatial resolution (≥250 m) satellite sensors with daily revisit capability are typically employed for this purpose. In recent years, phenology is being retrieved at higher resolution (≤30 m) in response to increasing availability of high-resolution satellite data. To overcome the reduced acquisition frequency of such data, previous attempts involved fusion between high- and medium-resolution data, or combinations of multi-year acquisitions in a single phenological reconstruction. The objectives of this study are to demonstrate that phenological parameters can now be retrieved from single-season high-resolution time series, and to compare these retrievals against those derived from multi-year high-resolution and single-season medium-resolution satellite data. The study focuses on the island of Schiermonnikoog, the Netherlands, which comprises a highly-dynamic saltmarsh, dune vegetation, and agricultural land. Combining NDVI series derived from atmospherically-corrected images from RapidEye (5 m-resolution) and the SPOT5 Take5 experiment (10m-resolution) acquired between March and August 2015, phenological parameters were estimated using a function fitting approach. We then compared results with phenology retrieved from four years of 30 m Landsat 8 OLI data, and single-year 100 m Proba-V and 250 m MODIS temporal composites of the same period. Retrieved phenological parameters from combined RapidEye/SPOT5 displayed spatially consistent results and a large spatial variability, providing complementary information to existing vegetation community maps. Retrievals that combined four years of Landsat observations into a single synthetic year were affected by the inclusion of years with warmer spring temperatures, whereas adjustment of the average phenology to 2015 observations was only feasible for a few pixels due to cloud cover around phenological transition dates. The Proba-V and MODIS phenology retrievals scaled poorly relative to their high-resolution equivalents, indicating that medium-resolution phenology retrievals need to be interpreted with care, particularly in landscapes with fine-scale land cover variability.