Normalized Difference Vegetation Index Time Series Data Research Articles

Reconstruction of dense time series high spatial resolution NDVI data using a spatiotemporal optimal weighted combination estimation model based on Sentinel-2 and MODIS

A regular, dense time series of the normalized difference vegetation index (NDVI) is a crucial remote sensing parameter that provides insights into the growth status of plants across different temporal and spatial scales. However, the production of high-resolution, equal-interval, and dense-time-series NDVI products currently faces technical challenges, limiting their application in agriculture and forestry-related fields. This study proposes a new method for reconstructing a dense Sentinel-2 NDVI time series based on the spatiotemporal optimal weighted combination estimation model (SOWCEM). The SOWCEM is developed using state-of-the-art spatial and temporal reconstruction algorithms. This model considers the distribution characteristics of errors in spatiotemporal dimensions, allowing for adaptive determination of the optimal combination weight at the pixel level. Then, the reconstructed regular dense-time-series NDVI data covering the region for five days/time are evaluated through a comparative analysis. The experimental results demonstrate that the regional NDVI reconstruction image based on the SOWCEM algorithm has a determination coefficient (R2) of 0.9082, a root mean square error (RMSE) of 0.0403, and a comprehensive overall error (ERGAS) of 5.7726. Compared with single spatiotemporal dimension models, the SOWCEM algorithm shows an average improvement of 5.78% in R2, while the RMSE and ERGAS decrease by 6.5% and 6.86%, respectively. Moreover, the stability and robustness of the models are also enhanced. These results indicate that the SOWCEM algorithm can effectively achieve a high-quality fusion reconstruction of regional images, significantly enhancing the accuracy of high-spatial-resolution NDVI dense temporal spectral reconstruction. Furthermore, it exhibits clear superiority and good applicability. This algorithmic framework provides a feasible approach for the dense reconstruction of high-spatial-resolution regular time-series NDVI data, enabling more accurate and efficient high-resolution NDVI remote sensing monitoring services in related field applications. The core code of SOWCEM is available at https://github.com/GISerZkun/SOWCEM.

NDVI time-series data reconstruction for spatial-temporal dynamic monitoring of Arctic vegetation structure

ABSTRACT Spatial-temporal dynamics monitoring of Arctic vegetation structure (i.e. distribution range of tundra and forest) is of great significance for evaluating global warming effect. Currently, time-series monitoring of Arctic vegetation structure relies primarily on the Normalized Difference Vegetation Index (NDVI), which is derived from optical remote sensing images. However, because of factors such as the long revisit period of satellites and the impact of climate, optical observations are severely lacking in the Arctic region. This results in NDVI time-series data highly discontinuous and difficult to reflect actual variations in Arctic vegetation structure, and the traditional time-series reconstruction method would usually fail for severe missing conditions. Therefore, this study developed a Time Series Reconstruction method considering Periodic Trend (TSR-PT), which is specifically for alleviating the severe missing observation condition in the Arctic region. It can separate the phenological change and trend change of the incomplete time series NDVI, and borrow the information from the neighboring unchanged years for compensate of the missing observations in current years, based on the learned inter-annual and intra-annual correlation. We explore its usability in monitoring vegetation structure variation in Vorkuta region (transition zone of tundra and taiga in the Arctic Circle) based on MODIS data. It is found that the proposed TSR-PT is able to reconstruct NDVI with reasonable phenological feature even the missing rate reaches over 70%, which is usually falsely constructed by traditional filtering or fitting method, and suppress them by 0.038 in terms of RMSE; besides, we find that since 21-century, the Arctic trees have continued to increase and encroach the original tundra ecosystem, which caused a largely Arctic vegetation structural change, and we believe the proposed method would largely promote the Arctic vegetation research.

Mapping and classification of Liao River Delta coastal wetland based on time series and multi-source GaoFen images using stacking ensemble model

The precise mapping of coastal wetlands holds great significance for monitoring carbon sequestration and storage within coastal ecosystems, particularly in light of climate change and human-induced activities. Time series and multi-source remote sensing data offer distinct advantages in spatial and temporal land use mapping, particularly in wetland systems encompassing various vegetation types. The primary aim of this study was to delineate the spatial distribution of land use within the Liao River Delta (LRD) wetland. This was achieved by employing a stacking ensemble model that integrates time-series GaoFen-1 (GF-1) optical imagery, GaoFen-3 (GF-3) synthetic aperture radar (SAR) imagery, and Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) imagery. The first step involved the application of an enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM) to fuse the GF-1 NDVI and MODIS NDVI datasets, producing time-series NDVI data. Subsequently, the parameters pertaining to vegetation phenology were obtained by employing the threshold method on the time-series NDVI data. We compiled feature datasets encompassing GF-1 spectral bands, indices, phenological parameters, and GF-3 SAR features. A Recursive Feature Elimination and Cross-Validation (RFECV) model was employed to identify and select significant features to mitigate data redundancy. Finally, a stacking ensemble model was constructed by combining five base models [K-Nearest Neighbors (KNN), Random Forest (RF), Adaptive Boosting (AdaBoost), Extreme Gradient Boosting (XGBoost), and a Light Gradient Boosting Machine (LightGBM)] to perform wetland classification. The findings were as follows: (1) ESTARFM was able to successfully fuse GF-1 NDVI and MODIS NDVI data, resulting in spatiotemporal fusion with a coefficient of determination (R2) of 0.85 and a root mean square error (RMSE) of 0.07. (2) the RFECV algorithm was employed to select relevant features in the wetland classification process. Specifically, 75 spectral band features, 89 spectral index features, 13 SAR features, and seven phenological parameters were identified as significant for this task. (3) A stacking ensemble model was constructed using the aforementioned multi-source features. This model exhibited a robust and consistent performance in wetland classification, achieving the highest overall accuracy of 94.33%. Notably, this accuracy improvement ranged from approximately 0.09% to 10.02% compared to the individual base models. Thus, the present study has the potential to be utilized for fine-scale wetland monitoring, thereby offering valuable assistance in the field of wetland environmental research.

A Robust Index Based on Phenological Features to Extract Sugarcane from Multisource Remote Sensing Data

Sugarcane is a major crop for sugar and biofuel production. Historically, mapping large sugarcane fields meticulously depended heavily on gathering comprehensive and representative training samples. This process was time-consuming and inefficient. Addressing this drawback, this study proposed a novel index, the Normalized Difference Vegetation Index (NDVI)-Based Sugarcane Index (NBSI). NBSI analyzed the temporal variation of sugarcane’s NDVI over a year. Leveraging the distinct growth phases of sugarcane (transplantation, tillering, rapid growth and maturity) four measurement methodologies, f(W1), f(W2), f(V) and f(D), were developed to characterize the features of the sugarcane growth period. Utilizing imagery from Landsat-8, Sentinel-2, and MODIS, this study employed the enhanced gap-filling (EGF) method to reconstruct NDVI time-series data for seven counties in Chongzuo, Guangxi Zhuang Autonomous Region, during 2021, subsequently testing NBSI’s ability to extract sugarcane. The results demonstrate the efficiency of NBSI with simple threshold settings: it was able to map sugarcane cultivation areas, exhibiting higher accuracy when compared to traditional classifiers like support vector machines (SVM) and random forests (RF), with an overall accuracy (OA) of 95.24% and a Kappa coefficient of 0.93, significantly surpassing RF (OA = 85.31%, Kappa = 0.84) and SVM (OA = 85.87%, Kappa = 0.86). This confirms the outstanding generalizability and robustness of the proposed method in Chongzuo. Therefore, the NBSI methodology, recognized for its flexibility and practicality, shows potential in enabling the extensive mapping of sugarcane cultivation. This heralds a new paradigm of thought in this field.

Remote sensing of evapotranspiration for irrigated crops at Yuma, Arizona, USA

A satellite-based vegetation index model that tracks daily crop growth and evapotranspiration (ETc) is developed, tested, and validated over irrigated farms in Yuma irrigation districts of Arizona and California. Model inputs are remotely sensed normalized difference vegetation index (NDVI) images, crop type maps, and local weather. The utility and novelty of the model is a more accurate assessment of ETc than currently provided by the US Bureau of Reclamation’s evapotranspiration modeling system. The model analyzes NDVI time series data from the European Space Agency’s Sentinel-2 satellites using the Google Earth Engine, constructs FAO-56 style crop growth stages from NDVI, and then estimates daily ETc using pre-defined crop coefficients (Kc) and grass reference evapotranspiration (ETos). Four crops were selected to test and evaluate model performance: short-season broccoli, mid-season cotton and wheat, and perennial alfalfa. Comparison of model results showed that Reclamation reports overestimate alfalfa and wheat ETc by 21–25%, cotton ETc by 6%, and underestimate broccoli ETc by 21%. Variability resolved by the model ranged 6–18% of median total ETc. Comparison of model results with those obtained from 13 eddy covariance sites showed validation discrepancies ranging 1–14%: average total actual ETc differences were 12, − 14, 78, and 87 mm/season, respectively, for alfalfa, broccoli, cotton, and wheat. The wide availability of Sentinel-2 data, collected every 5 days or less, and the rapid processing via Google Earth Engine make the vegetation index model implementation fast and practical. Its accuracy and ability to resolve ETc for every field would benefit the Reclamation water accounting system and provide valuable consumptive water use data for any Colorado River stakeholder.

Appraisal of long-term responsiveness of normalized difference vegetation index to climatic factors using multiscale time–frequency decomposition in an arid environment

An arid climate is a unique condition that has a significant impact on the growth of crops and natural vegetation. The normalized difference vegetation index (NDVI) is a crucial remotely sensed measurement of greenness due to its strong correlation with crop and vegetation growth and productivity. In the present study, the spatiotemporal dynamics of NDVI were analyzed from 2000 to 2021 in the segment of the arid western plain zone of Rajasthan, India. NDVI time-series data, as well as data related to climatic factors, viz., precipitation, soil moisture, evapotranspiration, and 2-m air temperature, were collected from Giovanni, the Goddard Earth Science dataset. The Mann–Kendall (MK) trend test and Sen’s slope depicted the long-term continuous time–frequency trend, while Karl Pearson’s correlation analysis depicted the significant relationship between all the factors except 2-m air temperature. The seasonal and mean monthly results of all the factors except 2-m air temperature showed considerable coherence with NDVI. The multiscale time–frequency decomposition or wavelet analysis depicted the fifth to the seventh month and the ninth to the 15th month of the cycle, showing the significance of the cropping pattern and the natural vegetation growth cycle. The cross-wavelet analysis further depicted important coherence, leading, and lagging phases among climatic factors and NDVI. Our research provided significant insights into the long-term variability and coherence of various climatic factors with NDVI that are applicable on regional and global scales.

A spatio-temporal fusion method based on Landsat and MODIS normalized vegetation index data

Abstract. Normalized difference vegetation index (NDVI) can simply and effectively reflect the growth status of plants, and has a linear relationship with vegetation coverage, thus it is an important indicator to identify vegetation growth status and coverage. However, the commonly used MODIS NDVI or Landsat NDVI cannot simultaneously achieve the high temporal and spatial resolution. To this end, we proposed a generic and automatic method to fuse MODIS and Landsat satellite vegetation index (NDVI) into daily high-resolution products which takes time series NDVI data from MODIS/Landsat as input, filters noise pixels, and generates high-resolution and high-frequency products through a fusion process of temporal-linearly interpolation to increase the temporal resolution and spatial filtering to alleviate the blocky artifacts from MODIS pixel border. The fused NDVI maps are visually compared to evaluate the performance of this method, and more quantitative assessment will be continued in the next step. This study can provide important support for large-scale long-term surface vegetation monitoring.

Trends of greening and browning in terrestrial vegetation in China from 2000 to 2020

Terrestrial vegetation condition is altering generally as a result of climate change and anthropogenic activity during the past few decades. To reveal the impact of long-term climate factors and artificial protection on multiple vegetation types, it is crucial to understand the spatial distribution of vegetation greening and browning and the effect of national ecological restoration programs. In this study, we established a persistent vegetation change index (P-value) to characterize greening (restoration) and browning (degradation) in China in 2000–2020. Firstly, we generated annual time-series normalized difference vegetation index (NDVI) data from MODIS product by averaging the monthly maximum NDVI values for each year. Secondly, we calculated the P-value to investigate the continuous change in vegetation state by incremental time interval. Finally, patterns and trends of greening and browning in forests, shrublands, and grasslands were quantified and mapped at pixel and sample point levels. The findings of our study revealed that Chinese wild vegetated lands greened up by ∼3.4 × 104 km2 (25%) and turned brown in ∼1.6 × 104 km2 (11%) between 2000 and 2020. Net greening was detected in all biomes, most conspicuously in several ecological program regions in northern China. The NDVI time-series data in 31% of field plots showed a consistent result, 11% of field plots showed a browning trend, and 58% of field plots showed a stable state. These results indicated a synergistic effect on forests, shrublands, and grasslands, but with regional variations attributed to differences in precipitation abundance, the implementation of positive ecological programs by the government and negative human activities. Additionally, these findings provide valuable insight into large-scale terrestrial vegetation transitions and have practical applications for decision-making and policy development in the assessment and restoration of ecosystems aimed at reducing carbon emissions, mitigating climate change, preserving biodiversity, and conserving water resources in China.

Changes in Vegetation Resistance and Resilience under Different Drought Disturbances Based on NDVI and SPEI Time Series Data in Jilin Province, China

Extreme drought is increasing in frequency and intensity in many regions globally. Understanding the changes in vegetation resistance and resilience under aggravated drought is essential for maintaining regional ecosystem stability. In this study, a drought event–vegetation response framework was developed to explore vegetation resistance and resilience changes. The normalized difference vegetation index (NDVI) was correlated with the standardized precipitation evapotranspiration index (SPEI) at multiple timescales to screen out the vegetation response time to drought. Then, the SPEI for the response time was detected using run theory to identify drought events during the period 2000–2017. Finally, drought-induced NDVI anomaly changes were identified using a sliding window to explore the changes in resistance and resilience to drought. This study focuses on Jilin province, China, which contains a famous environmentally vulnerable area. The results illustrate that the response time of vegetation to drought is 3 months. The northwest of Jilin province is considered to be drought-vulnerable because it has suffered the most drought events, i.e., 19–21 times, with severities in the range of 2.6–3.2 and durations in the range of 3.6–4.1 months. Grassland shows the weakest resistance and the strongest resilience, and tree cover shows the strongest resistance and the weakest resilience under severe drought disturbance among all vegetation. As the severity and duration of drought increase, the resistance decreases, and the resilience increases. During the growing season, the drought from May to July significantly impacts the vegetation resistance. Drought occurring from June to July has much less impact on resilience. Drought in August to September has less impact on resistance and a more significant impact on resilience. The results of this study may increase our knowledge regarding the response of vegetation to drought and guide ecosystem stability restoration.

Using Enhanced Gap-Filling and Whittaker Smoothing to Reconstruct High Spatiotemporal Resolution NDVI Time Series Based on Landsat 8, Sentinel-2, and MODIS Imagery

Normalized difference vegetation index (NDVI) time series data, derived from optical images, play a crucial role for crop mapping and growth monitoring. Nevertheless, optical images frequently exhibit spatial and temporal discontinuities due to cloudy and rainy weather conditions. Existing algorithms for reconstructing NDVI time series using multi-source remote sensing data still face several challenges. In this study, we proposed a novel method, an enhanced gap-filling and Whittaker smoothing (EGF-WS), to reconstruct NDVI time series (EGF-NDVI) using Google Earth Engine. In EGF-WS, NDVI calculated from MODIS, Landsat-8, and Sentinel-2 satellites were combined to generate high-resolution and continuous NDVI time series data. The MODIS NDVI was employed as reference data to fill missing pixels in the Sentinel–Landsat NDVI (SL-NDVI) using the gap-filling method. Subsequently, the filled NDVI was smoothed using a Whittaker smoothing filter to reduce residual noise in the SL-NDVI time series. With reference to the all-round performance assessment (APA) metrics, the performance of EGF-WS was compared with the conventional gap-filling and Savitzky–Golay filter approach (GF-SG) in Fusui County of Guangxi Zhuang Autonomous Region. The experimental results have demonstrated that the EGF-WS can capture more accurate spatial details compared with GF-SG. Moreover, EGF-NDVI of Fusui County exhibited a low root mean square error (RMSE) and a high coefficient of determination (R2). In conclusion, EGF-WS holds significant promise in providing NDVI time series images with a spatial resolution of 10 m and a temporal resolution of 8 days, thereby benefiting crop mapping, land use change monitoring, and various ecosystems, among other applications.

Quantitatively Computing the Influence of Vegetation Changes on Surface Discharge in the Middle-Upper Reaches of the Huaihe River, China

Changes in meteorology, hydrology, and vegetation will have significant impacts on the ecological environment of a basin, and the middle-upper reach of Huaihe River (MUHR) is one of the key regions for vegetation restoration in China. However, less studies have quantitatively accounted for the contribution of vegetation changes to land surface discharge in the MUHR. To quantitatively evaluate the influence of vegetation changes on land surface discharge in the MUHR, the Bernaola–Galavan (B–G) segmentation algorithm was utilized to recognize the mutation year of the Normalized Difference Vegetation Index (NDVI) time sequence data. Next, the functional relationship between the underlying surface parameter and the NDVI was quantitatively analyzed, and an adjusted Budyko formula was constructed. Finally, the effects of vegetation changes, climate factors, and mankind activities on the surface discharge in the MUHR were computed using the adjusted Budyko formula and elastic coefficient method. The results showed the following: (1) the surface runoff and precipitation from 1982 to 2015 in the MUHR presented a falling trend, yet the NDVI and potential evaporation presented an upward trend; (2) 2004 was the mutation year of the NDVI time series data, and the underlying surface parameter showed a significant linear regression relationship with the NDVI (p < 0.05); (3) the vegetation variation played a major role in the runoff variation during the changing period (2005–2015) in the MUHR. Precipitation, potential evaporation, and human activities accounted for −0.32%, −15.11%, and 18.24% of the surface runoff variation, respectively.

An L1-regularized variational approach for NDVI time-series reconstruction considering inter-annual seasonal similarity

Long-term normalized difference vegetation index (NDVI) data are extensively applied in environmental and ecological researches. However, cloud-induced interference and contamination seriously affects the quality of the current NDVI products, bringing huge uncertainty for subsequent applications. Although plenty of temporal methods have been developed to reconstruct NDVI time series, they still struggle to solve the challenge of temporally continuous gaps, due to the insufficient utilization of the prior characteristics in the time series. This paper develops a variational-based method to reconstruct multi-year NDVI time series by jointly regularizing local smoothness and inter-annual Seasonal similarity using the L1-norm (termed the SeasonL1 method). The proposed method innovatively introduces the information from adjacent years combined with assistance from temporal neighbors, and the L1-norm is imposed to establish the two corresponding regularization terms, to better characterize their statistical distributions. Simulated and real-data experiments were conducted in the Yangtze River Economic Belt of China using Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI data, to validate the performance of the SeasonL1 method by comparing with five classic methods. The results demonstrate that the proposed SeasonL1 method can achieve a satisfactory performance with an acceptable time cost in terms of both the quantitative indicators and spatio-temporal visual effect. In particular, the SeasonL1 method shown obvious advantages in recovering temporally continuous missing values and preventing over-smoothing in the inflection points of vegetation growth. We expect that the SeasonL1 method will become a useful filtering approach for obtaining high-quality NDVI time-series data in large-scale applications.

A novel method using explainable artificial intelligence (XAI)-based Shapley Additive Explanations for spatial landslide prediction using Time-Series SAR dataset

As artificial intelligence (AI) techniques are becoming more popular in landslide modeling, it is important to understand how decisions are made. Fairness, and transparency becomes ever more vital due to ethical concerns and trust. Despite the popularity of machine learning (ML) algorithms in landslide modeling, the explainability of these methods are often considered as black box. This paper aims to propose an explainable artificial intelligence (XAI) for landslide prediction using synthetic-aperture radar (SAR) time-series data, NDVI (normalized difference vegetation index) time-series data and other geo-environmental factors such as DEM (digital elevation model) derivatives. We employed a Shapley Additive Explanations (SHAP) approach to understand how and what decisions ML-based models are making. 37 features were extracted from various sources such as ALOS-PALSAR (ALOS Phased Array type l-band Synthetic Aperture Radar), ALOS-2 (SAR), Landsat-8, topographic maps, and DEM for landslide susceptibility mapping in a landslide prone area in Chukha, Bhutan as a test site. The result was then compared using two standard ML methods: random forest (RF) and support vector machine (SVM). As per results, the RF model outperformed (0.914) the SVM. Moreover, the higher reliability of the RF model was proved by the area under the curve (AUC) of 0.941. XAI results revealed, features like altitude, aspect, NDVI-2014, NDVI-2017, and NDVI-2018 were the most effective features for landslide prediction by both models. Interestingly, among those features, NDVI-2014, aspect, and NDVI-2017 negatively correlated with the landslide prediction; whereas positively correlated when SVM was utilized. This interpretation ability indicates the advantages of XAI over the conventional methods as it measures the impact, interaction and correlation of conditioning factors within a model. The current research finding can provide more transparency and explainability when working with MLs in landslide studies. This could help to build trust among the geoscientists and decision-makers while making geohazard prediction.

Vegetation changes in Yellow River Delta wetlands from 2018 to 2020 using PIE-Engine and short time series Sentinel-2 images

Vegetation is the functional subject in the wetland ecosystem and plays an irreplaceable role in biodiversity conservation. It is of great significance to monitor wetland vegetation for scientific assessment of the impact of vegetation on ecological environment and biodiversity. In this paper, a method for extracting wetland vegetation based on short time series Normalized Difference Vegetation Index (NDVI) data set was constructed. First, time series NDVI data were constructed using Sentinel-2 images. Then, the Support Vector Machine (SVM) classifier was used to classify the wetland vegetation types. The distributions of the main wetland vegetation in the study area in 2018 and 2020 were got. Finally, the land cover transfer matrix was calculated to analyze the spatial pattern and change of wetland vegetation emphatically from 2018 to 2020. Based on 46 Sentinel-2 images acquired in 2018 and 2020, the spatial pattern and change of vegetation in the Yellow River Delta wetlands were extracted and analyzed in this paper. The results show that: (1) The method for extracting wetland vegetation in estuary delta based on PIE-Engine platform and short time series NDVI data constructed in this paper can effectively extract the wetland vegetation information. The overall accuracy of the classification results reached 90.47% in 2018 and 80.30% in 2020. The Kappa coefficient of the classification results are 0.874 in 2018 and 0.739 in 2020 respectively. Compared with the results from the random forest classification method and the maximum likelihood classification method, the accuracy is improved by 6.40% and 13.04%, and the Kappa coefficient is improved by 0.055 and 0.069. (2) There were significant changes in vegetation coverage in the Yellow River Delta wetlands from 2018 to 2020. The Spartina alterniflora increased by 3.74km2. The Suaeda salsa degraded seriously, and the total area decreased by 20.38km2. In addition, the increase of Spartina alterniflora effectively guaranteed the stability of the coastline in the study area. This study can provide a theoretical basis for wetlands vegetation classificaton, and the classificaton results can provide scientific reference for protecting the ecological environment of wetlands and maintaining ecological stability.

Remote sensing based characterisation of community level phenological variations in a regional forest landscape of Western Ghats, India

The use of remote sensing for examining phenological variation in tropical forests is scarce. The major objectives of the study were to characterize the intra-annual variability of phenological cycle of the Biligiri Ranganathaswamy Temple Tiger Reserve (BRT) and the potentiality of these phenological metrices in defining species assemblages by classifying the forest. Sentinel-2 derived temporal Normalized Difference Vegetation Index (NDVI) data of 2019 was used to extract the vegetation trends and to derive phenological metrics using CropPhenology R package. Seasonal trends revealed that the highest greenness was associated with high NDVI values in September and October. We identified seven vegetation classes in the region and used Random Forest classifier to prepare a community level classification map with an overall classification accuracy of 68.9%. Our results revealed that incorporating the field sampling data and NDVI data can be effectively used for identifying, mapping and monitoring phenology of the BRT landscape..

Reconstructing High-Spatiotemporal-Resolution (30 m and 8-Days) NDVI Time-Series Data for the Qinghai–Tibetan Plateau from 2000–2020

As the largest and highest alpine ecoregion in the world, the Qinghai–Tibetan Plateau (QTP) is extremely sensitive to climate change and has experienced extraordinary warming during the past several decades; this has greatly affected various ecosystem processes in this region such as vegetation production and phenological change. Therefore, numerous studies have investigated changes in vegetation dynamics on the QTP using the satellite-derived normalized-difference vegetation index (NDVI) time-series data provided by the Moderate-Resolution Imaging Spectroradiometer (MODIS). However, the highest spatial resolution of only 250 m for the MODIS NDVI product cannot meet the requirement of vegetation monitoring in heterogeneous topographic areas. In this study, therefore, we generated an 8-day and 30 m resolution NDVI dataset from 2000 to 2020 for the QTP through the fusion of 30 m Landsat and 250 m MODIS NDVI time-series data. This dataset, referred to as QTP-NDVI30, was reconstructed by employing all available Landsat 5/7/8 images (>100,000 scenes) and using our recently developed gap-filling and Savitzky–Golay filtering (GF-SG) method. We improved the original GF-SG approach by incorporating a module to process snow contamination when applied to the QTP. QTP-NDVI30 was carefully evaluated in both quantitative assessments and visual inspections. Compared with reference Landsat images during the growing season in 100 randomly selected subregions across the QTP, the reconstructed 30 m NDVI images have an average mean absolute error (MAE) of 0.022 and a spatial structure similarity (SSIM) above 0.094. We compared QTP-NDVI30 with upscaled cloud-free PlanetScope images in some topographic areas and observed consistent spatial variations in NDVI between them (averaged SSIM = 0.874). We further examined an application of QTP-NDVI30 to detect vegetation green-up dates (GUDs) and found that QTP-NDVI30-derived GUD data show general agreement in spatial patterns with the 250 m MODIS GUD data, but provide richer spatial details (e.g., GUD variations at the subpixel scale). QTP-NDVI30 provides an opportunity to monitor vegetation and investigate land-surface processes in the QTP region at fine spatiotemporal scales.

LSTM-Based Prediction of Mediterranean Vegetation Dynamics Using NDVI Time-Series Data

Vegetation index time-series analysis of multitemporal satellite data is widely used to study vegetation dynamics in the present climate change era. This paper proposes a systematic methodology to predict the Normalized Difference Vegetation Index (NDVI) using time-series data extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS). The key idea is to obtain accurate NDVI predictions by combining the merits of two effective computational intelligence techniques; namely, fuzzy clustering and long short-term memory (LSTM) neural networks under the framework of dynamic time warping (DTW) similarity measure. The study area is the Lesvos Island, located in the Aegean Sea, Greece, which is an insular environment in the Mediterranean coastal region. The algorithmic steps and the main contributions of the current work are described as follows. (1) A data reduction mechanism was applied to obtain a set of representative time series. (2) Since DTW is a similarity measure and not a distance, a multidimensional scaling approach was applied to transform the representative time series into points in a low-dimensional space, thus enabling the use of the Euclidean distance. (3) An efficient optimal fuzzy clustering scheme was implemented to obtain the optimal number of clusters that better described the underline distribution of the low-dimensional points. (4) The center of each cluster was mapped into time series, which were the mean of all representative time series that corresponded to the points belonging to that cluster. (5) Finally, the time series obtained in the last step were further processed in terms of LSTM neural networks. In particular, development and evaluation of the LSTM models was carried out considering a one-year period, i.e., 12 monthly time steps. The results indicate that the method identified unique time-series patterns of NDVI among different CORINE land-use/land-cover (LULC) types. The LSTM networks predicted the NDVI with root mean squared error (RMSE) ranging from 0.017 to 0.079. For the validation year of 2020, the difference between forecasted and actual NDVI was less than 0.1 in most of the study area. This study indicates that the synergy of the optimal fuzzy clustering based on DTW similarity of NDVI time-series data and the use of LSTM networks with clustered data can provide useful results for monitoring vegetation dynamics in fragmented Mediterranean ecosystems.

Development of a new index for automated mapping of ratoon rice areas using time-series normalized difference vegetation index imagery

Over recent decades, the global demand for food has continued to grow, owing to population growth and the loss of arable land. Rice ratooning offers new opportunities for increasing rice production and has received renewed interest because of the minimal additional labor input required for its adoption. Regular, regional-scale monitoring of the spatial patterns of both traditional and ratoon rice cropping systems provides essential information for agricultural resource management and food security studies. However, the similar phenological characteristics of traditional double rice and ratoon rice cropping systems make it challenging to accurately classify these cropping practices based on satellite observations alone. In this study, we first proposed an improved phenology-based rice cropping area detection algorithm using moderate resolution imaging spectroradiometer (MODIS) normalized difference vegetation index (NDVI) imagery. A new index, ratoon rice index, was then developed to automatically delineate ratoon rice cropping areas with the aid of a base map of rice in Hubei Province, China. The accuracy assessment using ground truth data showed that our approach could map both traditional and ratoon rice cropping areas with high user accuracy (91.25% and 91.43%, respectively). The MODIS-retrieved rice cropping areas were validated using annual agricultural census data, and coefficient of determination ( R 2 ) values of 0.60 and 0.41 were recorded for traditional and ratoon rice cropping systems, respectively. The total area of ratoon rice was estimated to be 1 283.6 km 2 , 5.0% of the total rice cropping area, in Hubei Province in 2016. These demonstrated the feasibility of extracting the spatial patterns of both traditional and ratoon rice cropping systems solely from time-series NDVI and field survey data and strides made in facilitating the timely and routine monitoring of traditional and ratoon rice distribution at subnational level. Given sufficient historical satellite and phenology records, the proposed algorithm had the potential to enhance rice cropping area mapping efforts across a broad temporal scale ( e.g ., from the 1980s to the present).

Enhanced Spatial–Temporal Savitzky–Golay Method for Reconstructing High-Quality NDVI Time Series: Reduced Sensitivity to Quality Flags and Improved Computational Efficiency

The Spatial–Temporal Savitzky–Golay (STSG) method for noise reduction can address the problem of temporally continuous Normalized Difference Vegetation Index (NDVI) gaps and effectively increase local low NDVI values without overcorrection. However, STSG largely depends on the quality flags of the NDVI time-series data, and inaccurate quality flags yield misleading final results. STSG also requires extensive computing time when used in large-scale applications. This study proposes an enhanced method, called cuSTSG, to address the aforementioned limitations of STSG. First, cosine similarities between the annual NDVI time series were used to identify and exclude the NDVI values with inaccurate quality flags from the NDVI seasonal growth trajectory. Second, computational performance was improved by reducing redundant computations and parallelizing computationally intensive procedures using the Compute Unified Device Architecture (CUDA) on graphics processing units (GPUs). Experiments on four MODIS NDVI time-series datasets of various sizes and regions showed that compared with the original STSG, cuSTSG reduced the mean absolute errors of the final products by 4.90%, 7.77%, 11.76% and 2.06%, respectively. The results also showed that cuSTSG on a GPU achieved 75+ speed-up compared with the Interactive Data Language-implemented STSG, and 30+ speed-up compared with the C++-implemented STSG. cuSTSG can effectively mitigate the impacts of inaccurate quality flags on final products and generate high-quality NDVI time series at large scales with high accuracy and performance. The source code of cuSTSG is available at https://github.com/HPSCIL/cuSTSG.

Genetic dissection of seasonal vegetation index dynamics in maize through aerial based high-throughput phenotyping.

Plant phenotyping under field conditions plays an important role in agricultural research. Efficient and accurate high-throughput phenotyping strategies enable a better connection between genotype and phenotype. Unmanned aerial vehicle-based high-throughput phenotyping platforms (UAV-HTPPs) provide novel opportunities for large-scale proximal measurement of plant traits with high efficiency, high resolution, and low cost. The objective of this study was to use time series normalized difference vegetation index (NDVI) extracted from UAV-based multispectral imagery to characterize its pattern across development and conduct genetic dissection of NDVI in a large maize population. The time series NDVI data from the multispectral sensor were obtained at five time points across the growing season for 1,752 diverse maize accessions with a UAV-HTPP. Cluster analysis of the acquired measurements classified 1,752 maize accessions into two groups with distinct NDVI developmental trends. To capture the dynamics underlying these static observations, penalized-splines (P-splines) model was used to obtain genotype-specific curve parameters. Genome-wide association study (GWAS) using static NDVI values and curve parameters as phenotypic traits detected signals significantly associated with the traits. Additionally, GWAS using the projected NDVI values from the P-splines models revealed the dynamic change of genetic effects, indicating the role of gene-environment interplay in controlling NDVI across the growing season. Our results demonstrated the utility of ultra-high spatial resolution multispectral imagery, as that acquired using a UAV-based remote sensing, for genetic dissection of NDVI.