Normalized Difference Vegetation Index Time Series Research Articles

Use of HANTS algorithm in Sentinel-2 images for interpolation of missing data: a case-study for soy agricultural areas in periods of growth and harvesting

Abstract. Information about the current state of an agricultural culture can be fundamental in decision making, varying from quantity of fertilizers to be used in an area to estimating the productivity of a harvest in a period to establish credits payment. Satellite images with reasonable spatial resolution like Sentinel-2 are a powerful tool to visualize agricultural areas production, especially with the application of multispectral indexes like the Normalized Difference Vegetation Index (NDVI). Yet, due to the limitation of temporal resolution of satellites and periods with large cloudy weather can make the analysis difficult. This study proposed the use of the harmonic analysis time series algorithm (HANTS) to aid the incompleteness of time series obtained by the available satellite images in a case study of soy culture in Brazil. By applying the algorithm, a harmonic interpolation is obtained to produce a full daily NDVI time series from the beginning to end of the study period, facilitating further analysis to obtain metrics of interest. In this study, time series with gaps on specified periods are created based on full time series, then interpolated and compared with the real ones based on root mean squared error (RMSE) to access its accuracy.

Modeling NPP and NDVI time series in different bioclimatic regions of Iran.

Vegetation is one of the important components of ecosystems that usually changes seasonally. An accurate parameterization of vegetation cover dynamics by developing time series models can strengthen our understanding of vegetation change. This research aims to investigate and model the temporal changes of net primary production (NPP) and normalized difference vegetation index (NDVI) across bioclimatic regions of Iran, including the Khazari, Baluchi, semi-desert, steppe, semi-steppe, and arid forests. We used Moderate Resolution Imaging Spectroradiometer (MODIS) sensor products for NPP and NDVI time series (MOD17A2 and MOD13Q1, respectively). The SARIMA (Seasonal Autoregressive Integrated Moving Average) time series model is developed for NPP and NDVI time series. The investigation of autocorrelation functions (ACF) showed a strong seasonality in NPP and NDVI at the 12-month lag time. Comparing the lag times from 1 to 24 month for different regions shows that the NPP variable has a stronger seasonality. The evaluation of error criteria which showed NPP time series models based on RMSE, R2, MRE, and CE criteria was better, while based on the ME criteria, the models perform better for NDVI time series (for example, in Khazari region for NPP and NDVI time series, respectively, ME = 3.67, 0.05, RMSE = 0.12, 0.18, R2 = 0.87, 0.63, MRE = 0.02, 0.12, and CE = 0.84, 0.12). The selected models provided a short-term forecasting of the NPP and NDVI index for study regions at 24-month time, which is useful for the planning and management to reduce vegetation degradation and preserve ecosystem and biodiversity.

An improved framework for mapping and assessment of dynamics in cropping pattern and crop calendar from NDVI time series across a heterogeneous agro-climatic region.

The absence of spatial and temporal cropping information in semi-arid regions poses a significant challenge in assessing the dynamics of agricultural systems at river basin scales. Satellite remote sensing provides qualitative and quantitative information to derive vegetation dynamics over extensive areas of inherent complexities due to limitations in the availability of field data and the diverse nature of agricultural cropping practices. Utilizing phenological information derived from MODIS Normalized Difference Vegetation Index (NDVI) time series data from 2003-2004 to 2021-2022, this study derives major crop types, and cropping calendar (sowing, maturity, and harvest dates) for each season and year at 250-m resolution. This study introduces an integrated function-fitting model based on the Fast Fourier Transform (FFT) and a Lagrangian 3-point derivative-based approach to extract phenological events from the NDVI time series automatically. The derived phenological information is used in a hybrid crop classification algorithm that combines a rule-based decision tree approach (using the phenological events/dates) and a random forest classifier (using the phenological metrics) to generate the classified crop map at a river basin scale for multiple seasons and years. The implemented approach successfully captured sowing and harvesting dates for every crop growing season over 19years, with an RMSE of 9days, as observed with the available field survey data from 2015 to 2021. The classification results from this hybrid approach demonstrate an overall 82% accuracy. The proposed method shows a substantial improvement in cropping area estimation with a 60% reduction in MAE and RMSE compared to the existing algorithm. From the long time series of derived seasonal crop data (19-year analysis period), the influence of monsoonal activities and shifts on the spatial and temporal dynamics of sowing time and cropping patterns are assessed for a large river basin, demonstrating the utility of this continuous crop calendar. Further, an extensive analysis of results highlights farmers' adaptive strategies in response to dry and wet years, providing foundational insights for future studies on assessing the impacts of climate change. Hence, the hybrid framework adopted in this study for deriving a continuous crop calendar holds immense relevance for parameterizing and utilizing such information for developing river basin scale hydrologic and crop growth models for water resources planning and assessment.

Enhancing long-term vegetation monitoring in Australia: a new approach for harmonising the Advanced Very High Resolution Radiometer normalised-difference vegetation (NVDI) with MODIS NDVI

Abstract. Long-term, reliable datasets of satellite-based vegetation condition are essential for understanding terrestrial ecosystem responses to global environmental change, particularly in Australia, which is characterised by diverse ecosystems and strong interannual climate variability. We comprehensively evaluate several existing global Advanced Very High Resolution Radiometer (AVHRR) normalised-difference vegetation index (NDVI) products for their suitability for long-term vegetation monitoring in Australia. Comparisons with the MODIS NDVI highlight significant deficiencies, particularly over densely vegetated regions. Moreover, all the assessed products failed to adequately reproduce the interannual variability in the pre-MODIS era as indicated by Landsat NDVI anomalies. To address these limitations, we propose a new approach to calibrating and harmonising NOAA's Climate Data Record of AVHRR NDVI to the MODIS MCD43A4 NDVI for Australia using a gradient-boosting decision tree ensemble method. Two versions of the datasets are developed, one incorporating climate data in the predictors (“AusENDVI-clim”: Australian Empirical NDVI-climate) and another that is independent of climate data (“AusENDVI-noclim”). These datasets, spanning 1982–2013 at a spatial resolution of 0.05° and with a monthly time step, exhibit strong correlations (r2=0.89–0.94) and low mean errors compared with MODIS MCD43A4 NDVI (mean absolute error (MAE) = 0.014–0.028, RMSE = 0.021–0.046), accurately reproducing seasonal cycles over densely vegetated regions. Furthermore, they closely replicate the interannual variability in vegetation condition in the pre-MODIS era. A reliable method for gap-filling the AusENDVI record is also developed that leverages climate, atmospheric CO2 concentration, and woody-cover fraction predictors. The resulting synthetic NDVI dataset shows excellent agreement with the MODIS MCD43A4 NDVI and the recalibrated AVHRR NDVI time series (r2=0.82–0.95, MAE = 0.016–0.029, RMSE = 0.039–0.041). Finally, we provide a complete 41-year dataset where the gap-filled AusENDVI-clim from January 1982 to February 2000 is joined with the MODIS MCD43A4 NDVI from March 2000 to December 2022. Analysing 40-year per-pixel trends in Australia's annual maximum NDVI revealed increasing values, and shifts in the timing, of the annual peak NDVI across most of the continent, underscoring the dataset's potential to address crucial questions regarding the changing vegetation phenology and its drivers. The AusENDVI dataset can be used for studying Australia's changing vegetation dynamics and downstream impacts on the terrestrial carbon and water cycles, and it provides a reliable foundation for further research into the drivers of vegetation change. AusENDVI is open access and available at https://doi.org/10.5281/zenodo.10802703 (Burton et al., 2024).

Analysis of changes before and after forest fires with LAI, NDVI and ET time series: Focusing on major forest fires in Korea

This study investigates the changes caused by forest fires and the following recovery processes, targeting four forest fire sites in Korea. The time series of two vegetation indices, leaf area index (LAI) and normalized difference vegetation index (NDVI) are investigated along with evapotranspiration (ET) time series for the study objectives. The analysis results show that LAI is most sensitive to burn severity and burned area with its change ranging from 40 % to 70 %, while the changes of NDVI and ET remain 30 % and 20 %, respectively, regardless of forest fire sites. The recovery time from forest fire also varies according to indices: the recovery time is estimated to be about 15 years when considering LAI and NDVI, while just 5 years when considering ET. Overall, LAI seems better to analyze the change in vegetation before and after forest fires. Different vegetation shift patterns after the forest fire are also noticed, mostly from evergreen needleleaf trees to deciduous broadleaf trees. However, it is also found that bad soil fertility condition and artificial afforestation help to maintain the evergreen needleleaf trees after the forest fire.

Mapping Natural Populus euphratica Forests in the Mainstream of the Tarim River Using Spaceborne Imagery and Google Earth Engine

Populus euphratica is a unique constructive tree species within riparian desert areas that is essential for maintaining oasis ecosystem stability. The Tarim River Basin contains the most densely distributed population of P. euphratica forests in the world, and obtaining accurate distribution data in the mainstream of the Tarim River would provide important support for its protection and restoration. We propose a new method for automatically extracting P. euphratica using Sentinel-1 and 2 and Landsat-8 images based on the Google Earth Engine cloud platform and the random forest algorithm. A mask of the potential distribution area of P. euphratica was created based on prior knowledge to save computational resources. The NDVI (Normalized Difference Vegetation Index) time series was then reconstructed using the preferred filtering method to obtain phenological parameter features, and the random forest model was input by combining the phenological parameter, spectral index, textural, and backscattering features. An active learning method was employed to optimize the model and obtain the best model for extracting P. euphratica. Finally, the map of natural P. euphratica forests with a resolution of 10 m in the mainstream of the Tarim River was obtained. The overall accuracy, producer’s accuracy, user’s accuracy, kappa coefficient, and F1-score of the map were 0.96, 0.98, 0.95, 0.93, and 0.96, respectively. The comparison experiments showed that simultaneously adding backscattering and textural features improved the P. euphratica extraction accuracy, while textural features alone resulted in a poor extraction effect. The method developed in this study fully considered the prior and posteriori information and determined the feature set suitable for the P. euphratica identification task, which can be used to quickly obtain accurate large-area distribution data of P. euphratica. The method can also provide a reference for identifying other typical desert vegetation.

Long-term annual mapping and spatial-temporal dynamic analysis of winter wheat in Shandong Province based on spatial-temporal data fusion (2000-2022).

Winter wheat, as one of the world's key staple crops, plays a crucial role in ensuring food security and shaping international food trade policies. However, there has been a relative scarcity of high-resolution, long time-series winter wheat maps over the past few decades. This study utilized Landsat and Sentinel-2 data to produce maps depicting winter wheat distribution in Google Earth Engine (GEE). We further analyzed the comprehensive spatial-temporal dynamics of winter wheat cultivation in Shandong Province, China. The gap filling and Savitzky-Golay filter method (GF-SG) was applied to address temporal discontinuities in the Landsat NDVI (Normalized Difference Vegetation Index) time series. Six features based on phenological characteristics were used to distinguish winter wheat from other land cover types. The resulting maps spanned from 2000 to 2022, featuring a 30-m resolution from 2000 to 2017 and an improved 10-m resolution from 2018 to 2022. The overall accuracy of these maps ranged from 80.5 to 93.3%, with Kappa coefficients ranging from 71.3 to 909% and F1 scores from 84.2 to 96.9%. Over the analyzed period, the area dedicated to winter wheat cultivation experienced a decline from 2000 to 2011. However, a notable shift occurred with an increase in winter wheat acreage observed from 2014 to 2017 and a subsequent rise from 2018 to 2022. This research highlights the viability of using satellite observation data for the long-term mapping and monitoring of winter wheat. The proposed methodology has long-term implications for extending this mapping and monitoring approach to other similar areas.

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.

Field validation of NDVI to identify crop phenological signatures

Purpose and MethodsCrop identification using remotely sensed imagery provides useful information to make management decisions about land use and crop health. This research used phonecams to acquire the Normalized Difference Vegetation Index (NDVI) of various crops for three crop seasons. NDVI time series from Sentinel (L121-L192) images was also acquired using Google Earth Engine (GEE) for the same period. The resolution of satellite data is low therefore gap filling and smoothening filters were applied to the time series data. The comparison of data from satellite images and phenocam provides useful insight into crop phenology. The results show that NDVI is generally underestimated when compared to phenocam data. The Savitzky-Golay (SG) and some other gap filling and smoothening methods are applied to NDVI time series based on satellite images. The smoothened NDVI curves are statistically compared with daily NDVI series based on phenocam images as a reference.ResultsThe SG method has performed better than other methods like moving average. Furthermore, polynomial order has been found to be the most sensitive parameter in applying SG filter in GEE. Sentinel (L121-L192) image was used to identify wheat during the year 2022–2023 in Sargodha district where experimental fields were located. The Random Forest Machine Leaning algorithm was used in GEE as a classifier.ConclusionThe classification accuracy has been found 97% using this algorithm which suggests its usefulness in applying to other areas with similar agro-climatic characteristics.

Enhancement of Comparative Assessment Approaches for Synthetic Aperture Radar (SAR) Vegetation Indices for Crop Monitoring and Identification—Khabarovsk Territory (Russia) Case Study

Crop identification at the field level using remote sensing data is a very important task. However, the use of multispectral data for the construction of vegetation indices is sometimes impossible or limited. For such situations, solutions based on the use of time series of synthetic aperture radar (SAR) indices are promising, eliminating the problems associated with cloudiness and providing an assessment of crop development characteristics during the growing season. We evaluated the use of time series of synthetic aperture radar (SAR) indices to characterize crop development during the growing season. The use of SAR imagery for crop identification addresses issues related to cloudiness. Therefore, it is important to choose the SAR index that is the most stable and has the lowest spatial variability throughout the growing season while being comparable to the normalized difference vegetation index (NDVI). The presented work is devoted to the study of these issues. In this study, the spatial variabilities of different SAR indices time series were compared for a single region for the first time to identify the most stable index for use in precision agriculture, including the in-field heterogeneity of crop sites, crop rotation control, mapping, and other tasks in various agricultural areas. Seventeen Sentinel-1B images of the southern part of the Khabarovsk Territory in the Russian Far East at a spatial resolution of 20 m and temporal resolution of 12 days for the period between 14 April 2021 and 1 November 2021 were obtained and processed to generate vertical–horizontal/vertical–vertical polarization (VH/VV), radar vegetation index (RVI), and dual polarimetric radar vegetation index (DpRVI) time series. NDVI time series were constructed from multispectral Sentinel-2 images using a cloud cover mask. The characteristics of time series maximums were calculated for different types of crops: soybean, oat, buckwheat, and timothy grass. The DpRVI index exhibited the highest stability, with coefficients of variation of the time series that were significantly lower than those for RVI and VH/VV. The main characteristics of the SAR and NDVI time series—the maximum values, the dates of the maximum values, and the variability of these indices—were compared. The variabilities of the maximum values and dates of maximum values for DpRVI were lower than for RVI and VH/VV, whereas the variabilities of the maximum values and the dates of maximum values were comparable for DpRVI and NDVI. On the basis of the DpRVI index, classifications were carried out using seven machine learning methods (fine tree, quadratic discriminant, Gaussian naïve Bayes, fine k nearest neighbors or KNN, random under-sampling boosting or RUSBoost, random forest, and support vector machine) for experimental sites covering a total area of 1009.8 ha. The quadratic discriminant method yielded the best results, with a pixel classification accuracy of approximately 82% and a kappa value of 0.67. Overall, 90% of soybean, 74.1% of oat, 68.9% of buckwheat, and 57.6% of timothy grass pixels were correctly classified. At the field level, 94% of the fields included in the test dataset were correctly classified. The paper results show that the DpRVI can be used in cases where the NDVI is limited, allowing for the monitoring of phenological development and crop mapping. The research results can be used in the south of Khabarovsk Territory and in neighboring territories.

Unpacking satellite pixels: UAVs reveal fine-scale drivers of land surface phenology in a winter rainfall shrubland

Land surface phenology (LSP) can reveal important connections between vegetation dynamics and climate but remains poorly understood in evergreen winter-rainfall shrublands globally. Field-based studies have indicated diverse plant functional strategies in shrublands, but further work is required to link LSP to vegetation functional composition in these regions. We analysed time-series of the normalised difference vegetation index (NDVI) in fynbos shrublands of South Africa using multi-spectral imagery from satellites and unmanned aerial vehicles (UAVs). We investigated the climate drivers of seasonal vegetative phenology and long-term NDVI trends at multiple spatial scales ranging from the landscape to individual species. At coarse spatial resolutions, NDVI time-series indicated rainfall-driven vegetation dynamics in fynbos, both at inter and intra-annual time scales. However, high-resolution time-series from UAVs exposed an underlying divergence in vegetative phenology and long-term NDVI trends between shallow and deep-rooted growth forms. Phenophases and NDVI trends of isolated, deep-rooted, overstory shrubs were decoupled from rainfall relative to dense overstory patches and shallow-rooted understory growth forms. Variations in growth form phenology were not detected at coarse spatial scales due to scaling and competitive effects based on the functional composition of the vegetation.

Reconstruction of NDVI based on Larix gmelinii tree-rings during June–September 1759–2021

Investigating the long-term dynamics in the canopy proves to be difficult due to the short observational records of the normalized difference vegetation index (NDVI). To explore the linkage between tree growth, NDVI dynamics and large-scale atmospheric circulation in the Greater Khingan Mountain, Northeast China, we established a chronology of Larix gmelinii tree ring width at three elevations (870–920 m, 1,100–1,150 m and 1,270–1,320 m) in the northern foothills of the mountain range. We then calculated the correlations between the tree ring chronologies and NDVI and climate factors, and reconstructed the NDVI time series from June to September 1759–2021 in the region based on the middle-elevation tree ring chronology. The results identify the positive effect of temperature (r = 0.56, p < 0.01) and the negative effect of precipitation (r = −0.44, p < 0.01) in the growing season as the main influencing factors of NDVI for the study period (1981–2019). The 11-year moving average of the reconstructed NDVI series reveals two periods of high canopy vigor (1898–1926 and 2009–2013) and three periods of low canopy vigor (1860–1962, 1882–1888 and 1968–1977) in the last 263 years. These periods correspond to drought events recorded in the historical literature. Wavelet analysis shows that the reconstructed sequences exhibited 11–13, 23–25, and 39–42 years period variations. Integrating this with spatial correlation analysis reveals that tree growth in the Mangui region was impacted by the combined effect of the North Pacific Decadal Oscillation and the North Atlantic Multidecadal Oscillation. The results of this paper provide a reference for the study of vegetation change patterns in the northern foothills of the Greater Khingan Mountains.

Using Geostationary Satellite Observations to Improve the Monitoring of Vegetation Phenology

Geostationary satellite data enable frequent observations of the Earth’s surface, facilitating the rapid monitoring of land covers and changes. However, optical signals over vegetation, represented by the vegetation index (VI), exhibit an anisotropic effect due to the diurnal variation in the solar angle during data acquisition by geostationary satellites. This effect, typically characterized by the bi-directional reflectance distribution function (BRDF), can introduce uncertainties in vegetation monitoring and the estimation of phenological transition dates (PTDs). To address this, we investigated the diurnal variation in the normalized difference vegetation index (NDVI) with solar angles obtained from geostationary satellites since the image had fixed observation angles. By establishing a temporal conversion relationship between instantaneous NDVI and daily NDVI at the local solar noon (LSNVI), we successfully converted NDVIs obtained at any time during the day to LSNVI, increasing cloud-free observations of NDVI by 34%. Using different statistics of the time series vegetation index, including LSNVI, daily averaged NDVI (DAVI), and angular corrected NDVI (ACVI), we extracted PTD at five typical sites in China. The results showed a difference of up to 41.5 days in PTD estimation, with the highest accuracy achieved using LSNVI. The use of the proposed conversion approach, utilizing time series LSNVI, reduced the root mean square error (RMSE) of PTD estimation by 9 days compared with the use of actual LSNVI. In conclusion, this study highlights the importance of eliminating BRDF effects in geostationary satellite observations and demonstrates that the proposed angular normalization method can enhance the accuracy of time series NDVI in vegetation monitoring.

NDVI estimation using Sentinel-1 data over wheat fields in a semiarid Mediterranean region

ABSTRACT Annual crop monitoring is a key parameter for managing agricultural strategies. Several studies have relied on remote sensing products such as the normalized difference vegetation index (NDVI) as a vegetation dynamic metric. However, the dependence of optical data on weather conditions limits its availability. In this study, we reconstruct the NDVI time series of wheat fields using the moving averages of the Sentinel-1 normalized VH/VV cross-polarization ratio (IN) and the interferometric coherence in VV polarization over wheat selected fields in a semiarid site in Tunisia during two seasons, from 2018 to 2020. The crop cycle is divided into two periods: before and after the heading phase, which occurs in approximately the middle of March. Due to the volume-scattering impact, the second phase is divided into the ripening and maturation phase (NDVI ≥0.4) and senescence phase (NDVI <0.4). To estimate the NDVI values, different methods are used: curve-fitting equations and machine learning regressors such as the random forest (RF) and the support vector regressor (SVR). Low root mean square error (RMSE) values characterize NDVI estimation during the first period. In the second period, the RMSE values reach 0.06 when the NDVI is lower than 0.4. When the NDVI values exceed 0.4 in the second period, lower accuracy marks the NDVI estimation using the curve-fitting equations as a function of IN or coherence. Relative low accuracy characterizes the regression algorithms’ estimations when NDVI ≥ 0.4 compared to their performance during the aforementioned periods. The proposed approach was tested on different wheat fields. The NDVI estimations are characterized by RMSE values varying between 0.12 and 0.19. The use of RF and SVR outperformed the curve-fitting methods with an RMSE equal to 0.12. The present findings revealed the high accuracy of the proposed approach to estimate the missing values of wheat fields NDVI values during the vegetation development period until heading and the senescence phase. The presence of the mutual effect of the vegetation water content and its volume complicated the NDVI estimation using the C-band data.

Harmonized NDVI time-series from Landsat and Sentinel-2 reveal phenological patterns of diverse, small-scale cropping systems in East Africa

Capturing the phenological responses of crops in diverse, small-scale agricultural systems is one of the most notable remote sensing research gaps in Eastern Africa. Because of the heterogeneous and complex agricultural environments, remote sensing methodologies successfully implemented elsewhere yield poor results in Eastern Africa. The reliance on local rainfall and dry spells can lead to alterations in planting patterns, impacting growing season timing and yield, underscoring the importance of high temporal resolution for capturing crop phenology. We generated a high-resolution land use/land cover (LULC) classification map of 10 m using Sentinel-1 and Sentinel-2 imagery together with a large set of reference data acquired in situ in 2022. The LULC classification yielded an overall accuracy of 87.9% and Cohen's kappa value of 0.83. Based on the LULC map, we investigated the phenological metrics of the most important crops in the Kamuli district of Eastern Uganda by combining a multi-sensor harmonization approach, thresholding method, and cloud computing. Previously unavailable Landsat 9 band adjustment coefficients were derived. By harmonizing Landsat (LS) 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 and 9 Operational Land Imager (OLI), with Sentinel-2 Multispectral Instrument (MSI), the phenological patterns during the March to May growing season were captured using a dense temporal Normalized Difference Vegetation Index (NDVI) time series. The seasonal onset for most crop types was in the end of March to middle of April while the end of the growing season occurred from the end of June to middle of July. Investigation of the 13 crop types using a high-resolution LULC map and dense, harmonized NDVI time series signifies the first field-level phenological metrics estimation of smallholder farming systems in Eastern Africa. Future research should focus on expanding the methodology in space and time to provide insights into past and current growing season development.

Extracting the winter wheat using the decision tree based on time series dual-polarization SAR feature and NDVI.

Winter wheat is one of the most important crops in the world. It is great significance to obtain the planting area of winter wheat timely and accurately for formulating agricultural policies. Due to the limited resolution of single SAR data and the susceptibility of single optical data to weather conditions, it is difficult to accurately obtain the planting area of winter wheat using only SAR or optical data. To solve the problem of low accuracy of winter wheat extraction only using optical or SAR images, a decision tree classification method combining time series SAR backscattering feature and NDVI (Normalized Difference Vegetation Index) was constructed in this paper. By synergy using of SAR and optical data can compensate for their respective shortcomings. First, winter wheat was distinguished from other vegetation by NDVI at the maturity stage, and then it was extracted by SAR backscattering feature. This approach facilitates the semi-automated extraction of winter wheat. Taking Yucheng City of Shandong Province as study area, 9 Sentinel-1 images and one Sentinel-2 image were taken as the data sources, and the spatial distribution of winter wheat in 2022 was obtained. The results indicate that the overall accuracy (OA) and kappa coefficient (Kappa) of the proposed method are 96.10% and 0.94, respectively. Compared with the supervised classification of multi-temporal composite pseudocolor image and single Sentinel-2 image using Support Vector Machine (SVM) classifier, the OA are improved by 10.69% and 5.66%, respectively. Compared with using only SAR feature for decision tree classification, the producer accuracy (PA) and user accuracy (UA) for extracting the winter wheat are improved by 3.08% and 8.25%, respectively. The method proposed in this paper is rapid and accurate, and provide a new technical method for extracting winter wheat.

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.

Modelling vegetation phenology at six field stations within the U.S. Great Plains: constructing a 38-year timeseries of GCC, VCI, NDVI, and EVI2 using PhenoCam imagery and DAYMET meteorological records

Vegetation phenology models still rely on temperature as the primary limiting factor to growth. They generally do not recognize the importance of photoperiod and water availability, which can cause them to under-perform. Moreover, few models have used machine learning algorithms to find relationships in the data. In this paper, four Vegetation Indexes (VIs), namely the green chromatic coordinate (GCC), the vegetation contrast index (VCI), the normalized difference vegetation index (NDVI) and the two-band enhanced vegetation index (EVI2), are predicted for the North American Great Plains. This is possible by using six PhenoCams, Daily Surface Weather and Climatological Summaries (DAYMET), processing them with the machine learning algorithm XGBoost (XGB) and comparing them with seven phenophase stages throughout a growth cycle. Examining the results, GCC was the best fitting model with an R2 of 0.946, while EVI2 was the poorest with an R2 of 0.895. Also, the results indicate that changing temperature and precipitation patterns are driving a significant change in phenology of the grasslands. We developed a model capable of explaining 90 to 93% of the variability in four VIs across six grassland PhenoCam sites over the growing season using the XGB regression. Our model demonstrates the importance of including photoperiod, temperature, and precipitation information when modeling vegetation phenology. Finally, we were able to construct a 38-year phenology record at each PhenoCam location.

Monitoring Grassland Variation in a Typical Area of the Qinghai Lake Basin Using 30 m Annual Maximum NDVI Data

The normalized difference vegetation index (NDVI) can depict the status of vegetation growth and coverage in grasslands, whereas coarse spatial resolution, cloud cover, and vegetation phenology limit its applicability in fine-scale research, especially in areas covering various vegetation or in fragmented landscapes. In this study, a methodology was developed for obtaining the 30 m annual maximum NDVI to overcome these shortcomings. First, the Landsat NDVI was simulated by fusing Landsat and MODIS NDVI by using the enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM), and then a single-peaked symmetric logistic model was employed to fit the Landsat NDVI data and derive the maximum NDVI in a year. The annual maximum NDVI was then used as a season-independent substitute to monitor grassland variation from 2001 to 2022 in a typical area covering the major vegetation types in the Qinghai Lake Basin. The major conclusions are as follows: (1) Our method for reconstructing the NDVI time series yielded higher accuracy than the existing dataset. The root mean square error (RMSE) for 91.8% of the pixels was less than 0.1. (2) The annual maximum NDVI from 2001 to 2022 exhibited spatial distribution characteristics, with higher values in the northern and southern regions and lower values in the central area. In addition, the earlier vegetation growth maximum dates were related to the vegetation type and accompanied by higher NDVI maxima in the study area. (3) The overall interannual variation showed a slight increasing trend from 2001 to 2022, and the degraded area was characterized as patches and was dominated by Alpine kobresia spp., Forb Meadow, whose change resulted from a combination of permafrost degradation, overgrazing, and rodent infestation and should be given more attention in the Qinghai Lake Basin.

The Retrieval of Ground NDVI (Normalized Difference Vegetation Index) Data Consistent with Remote-Sensing Observations

The Normalized Difference Vegetation Index (NDVI) is widely used for monitoring vegetation status, as accurate and reliable NDVI time series are crucial for understanding the relationship between environmental conditions, vegetation health, and productivity. Ground digital cameras have been recognized as important potential data sources for validating remote-sensing NDVI products. However, differences in the spectral characteristics and imaging methods between sensors onboard satellites and ground digital cameras hinder direct consistency analyses, thereby limiting the quantitative application of camera-based observations. To address this limitation and meet the needs of vegetation monitoring research and remote-sensing NDVI validation, this study implements a novel NDVI camera. The proposed camera incorporates narrowband dual-pass filters designed to precisely separate red and near-infrared (NIR) spectral bands, which are aligned with the configuration of sensors onboard satellites. Through software-controlled imaging parameters, the camera captures the real radiance of vegetation reflection, ensuring the acquisition of accurate NDVI values while preserving the evolving trends of the vegetation status. The performance of this NDVI camera was evaluated using a hyperspectral spectrometer in the Hulunbuir Grassland over a period of 93 days. The results demonstrate distinct seasonal characteristics in the camera-derived NDVI time series using the Green Chromatic Coordinate (GCC) index. Moreover, in comparison to the GCC index, the camera’s NDVI values exhibit greater consistency with those obtained from the hyperspectral spectrometer, with a mean deviation of 0.04, and a relative root mean square error of 9.68%. This indicates that the narrowband NDVI, compared to traditional color indices like the GCC index, has a stronger ability to accurately capture vegetation changes. Cross-validation using the NDVI results from the camera and the PlanetScope satellite further confirms the potential of the camera-derived NDVI data for consistency analyses with remote sensing-based NDVI products, thus highlighting the potential of camera observations for quantitative applications The research findings emphasize that the novel NDVI camera, based on a narrowband spectral design, not only enables the acquisition of real vegetation index (VI) values but also facilitates the direct validation of vegetation remote-sensing NDVI products.