Normalized Difference Vegetation Index Products 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.

Spatio-temporal dynamics of vegetation over cloudy areas in Southwest China retrieved from four NDVI products

The Normalized Difference Vegetation Index (NDVI) is the most commonly used index for assessing vegetation. However, significant differences among various satellite datasets, complex terrain, and the impact of clouds on optical sensors limit vegetation change assessment based on NDVI. To address these issues, this study utilizes multi-source satellite data (GIMMS3g NDVI, CDR AVHRR NDVI, SPOT NDVI, and MODIS NDVI) to monitor vegetation dynamics at different time scales from 1990 to 2020 in Sichuan Province, China. The results indicate that over time, NDVI values from the four NDVI products in Sichuan Province have shown an upward trend. There are certain differences in the spatial distribution and spatial heterogeneity of the change rate of NDVI values among the four NDVI products at different time scales, and the differences are mainly concentrated in the Sichuan Basin (SB) and the Western Sichuan alpine plateau region (WS). Compared with the other three NDVI products, GIMMS NDVI has the highest value but the smallest increase during the study period. The SPOT NDVI value is the smallest, but the increase is relatively large. However, within the overlapping period of the four NDVI datasets, only the annual average of CDR AVHRR NDVI showed a downward trend (slope2000–2013 = −0.0001·a−1). The annual fluctuation of CDR AVHRR NDVI is the smallest, and compared to other NDVI datasets, its correlation with climate factors shows significantly weaker spatial variability. Moreover, the ability of CDR AVHRR NDVI to distinguish different vegetation land cover types is significantly poor (STD = 0.045). The findings of this study will provide a reference for further research on vegetation changes in Sichuan Province and NDVI reconstruction in cloudy areas.

Development of long‐term spatiotemporal continuous NDVI products for alpine grassland from 1982 to 2020 in the Qinghai–Tibet Plateau, China

AbstractBackgroundThe time‐series data of the Normalized Difference Vegetation Index (NDVI) is a crucial indicator for global and regional vegetation monitoring. However, the current assessment of global and regional long‐term vegetation changes is subject to large uncertainties due to the lack of spatiotemporally continuous time‐series data sets.MethodsIn this study, a long time‐series monthly NDVI data set with a spatial resolution of 250 m from 1982 to 2020 was developed by combining Moderate Resolution Imaging Spectroradiometer (MODIS) and AVHRR (Advanced Very High‐Resolution Radiometer) time‐series NDVI products using the Random Forest (RF) downscaling model.ResultsCompared to the MODIS NDVI product, the fused product shows RMSE and mean absolute error ranging from 0 to 0.075 and from 0 to 0.05, respectively, with R2 values mostly above 0.7.ConclusionsThe long time‐series NDVI products generated in this study are reliable in terms of accuracy and have great potential for long‐term dynamic monitoring of terrestrial ecosystems on the Qinghai–Tibet Plateau.

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.

Improved global 250 m 8-day NDVI and EVI products from 2000–2021 using the LSTM model

Satellite vegetation index (VI) products, such as normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI), have been widely used. However, they are severely contaminated by clouds and other factors and provide false signals of the surface vegetation conditions. In this study, the new global seamless 250 m, eight-day NDVI and EVI products from 2000–2021 were developed from Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data using a long short-term memory (LSTM) neural network method. High-quality globally representative time series VI samples were constructed to train the model using a combination of the Savitzky-Golay filter (SG), Global LAnd Surface Satellite (GLASS) leaf area index (LAI) fitting and upper envelope methods. To evaluate the proposed method and the 250 m VI products, the MODIS VI product (MOD13Q1) was used for the inter-comparisons using four widely used VI reconstruction methods. Assuming that the MODIS VI data of high quality represents the true values, the root mean square error (RMSE) for NDVI and EVI generated by the LSTM model are 0.0734 and 0.0509, respectively.

Exploiting satellite data for total direct runoff prediction using CN-based MSME model

This paper explores the potential to enhance the functionality of the modified Sahu-Mishra-Eldho model (MSME-CN) using indirect soil moisture measurements derived from satellite data. The current version of the MSME-CN model is not applicable in ungauged watersheds due to the necessity of calibrating the crucial parameter α, which reflects soil saturation, based on measured rainfall-runoff events. We hypothesize that the Normalized Difference Vegetation Index (NDVI) can serve as an indirect indicator of soil moisture to assess the soil saturation parameter α in the MSME model. This hypothesis was tested across five different watersheds, three located in the southeastern USA and two in southern Poland. The NDVI product, developed from data obtained from the Advanced Very High-Resolution Radiometer (AVHRR), was utilized in this study. Results indicate that NDVI is a robust indicator of soil moisture for representing the α parameter in the MSME model. The correlation coefficient between α and NDVI a day prior to a rainfall event was around 0.80 for the WS80 and Kamienica watersheds and nearly 0.60 for the other watersheds. The analysis corroborates the hypothesis that NDVI can serve as an indirect parameter of soil moisture to assess the soil saturation parameter α in the MSME-CN model. Based on Nash-Sutcliffe Efficiency (NSE) statistics, the total direct runoff predicted by the MSME-CN model, with the α parameter updated using NDVI, was rated ‘very good’ for the WS80 and AC11 watersheds, ‘good’ for the Kamienica watershed, ‘satisfactory’ for Stobnica, and ‘unsatisfactory’ for the high forest density WS14 watershed, potentially highlighting the model's limitation in such watersheds.

Spatiotemporally consistent global dataset of the GIMMS Normalized Difference Vegetation Index (PKU GIMMS NDVI) from 1982 to 2022

Abstract. Global products of remote sensing Normalized Difference Vegetation Index (NDVI) are critical to assessing the vegetation dynamic and its impacts and feedbacks on climate change from local to global scales. The previous versions of the Global Inventory Modeling and Mapping Studies (GIMMS) NDVI product derived from the Advanced Very High Resolution Radiometer (AVHRR) provide global biweekly NDVI data starting from the 1980s, being a reliable long-term NDVI time series that has been widely applied in Earth and environmental sciences. However, the GIMMS NDVI products have several limitations (e.g., orbital drift and sensor degradation) and cannot provide continuous data for the future. In this study, we presented a machine learning model that employed massive high-quality global Landsat NDVI samples and a data consolidation method to generate a new version of the GIMMS NDVI product, i.e., PKU GIMMS NDVI (1982–2022), based on AVHRR and Moderate-Resolution Imaging Spectroradiometer (MODIS) data. A total of 3.6 million Landsat NDVI samples that were well spread across the globe were extracted for vegetation biomes in all seasons. The PKU GIMMS NDVI exhibits higher accuracy than its predecessor (GIMMS NDVI3g) in terms of R2 (0.97 over 0.94), root mean squared error (RMSE: 0.05 over 0.09), mean absolute error (MAE: 0.03 over 0.07), and mean absolute percentage error (MAPE: 9 % over 20 %). Notably, PKU GIMMS NDVI effectively eliminates the evident orbital drift and sensor degradation effects in tropical areas. The consolidated PKU GIMMS NDVI has a high consistency with MODIS NDVI in terms of pixel value (R2 = 0.956, RMSE = 0.048, MAE = 0.034, and MAPE = 6.0 %) and global vegetation trend (0.9×10-3 yr−1). The PKU GIMMS NDVI product can potentially provide a more solid data basis for global change studies. The theoretical framework that employs Landsat data samples can facilitate the generation of remote sensing products for other land surface parameters. The PKU GIMMS NDVI product is open access and available under a Creative Commons Attribution 4.0 License at https://doi.org/10.5281/zenodo.8253971 (Li et al., 2023).

Trend of Changes in Phenological Components of Iran’s Vegetation Using Satellite Observations

Investigating vegetation changes, especially plant phenology, can yield valuable information about global warming and climate change. Time series satellite observations and remote sensing methods offer a great source of information on distinctions and changing aspects of vegetation. The current study aimed to determine the trend and rate of changes in some phenological components of Iran’s vegetation. In this regard, the current study employed the daily NDVI (Normalized Difference Vegetation Index) product of the AVHRR sensor with a spatial resolution of 0.05° × 0.05°, named AVH13C1. Then, using the HANTS algorithm, images of amplitude zero, annual amplitude, and annual phase were prepared annually from 1982 to 2019. Using TIMESAT software, the starting, end, and length of time of growing season were calculated for each pixel time series to prepare annual maps. The Mann–Kendall statistical test was used to investigate the significance of changes during the study period. On average in the entire area of Iran, the annual phase was declining with a trend of −0.6° per year, and the time for the start and end of the season was declining by −0.3 and −0.65 days per year, respectively. Major changes were noticed in the northeast, west, and northwest regions of Iran, where the annual phase declined with a trend of −0.9° per year. Since the annual growth cycle of the plant (equivalent to 356 days) was in the form of a sinusoidal signal, and the angular changes in the sine wave were between zero and 360°, each degree of change was equivalent to 1.01 days per year. Therefore, the reduction in the annual phase by −0.9 degrees almost means a change in the time (due to the earlier negative start phase) of the start of the annual growth signal by −0.9 days per year. The time of the start and end of the growing season declined by −0.6 and −1.33 days per year, respectively. The reduction in annual phase and differences in time of the starting season from 1982 to 2019 indicate the acceleration and earlier initiation of various phenological processes in the area.

Polar-facing slopes showed stronger greening trend than equatorial-facing slopes in Tibetan plateau grasslands

The orientation of slopes in alpine zones creates microclimates, e.g. equatorial-facing slopes (EFSs) are generally drier and warmer than are polar-facing slopes (PFSs). The vegetation growing in these microhabitats responds divergently to climatic warming depending on the slope orientation. We proposed a spatial metric, the greenness asymmetric index (GAI), defined as the ratio between the average normalized difference vegetation index (NDVI) on PFSs and EFSs within a given spatial window, to quantify the asymmetry of greenness across aspects. We calculated GAI for each non-overlapping 3 × 3 km2 (100 × 100 Landsat pixels) grid, and seamlessly mapped it on Tibetan Plateau (TP) grassland using NDVI time series from the Landsat-5, -7 and -8 satellites. PFSs were greener than EFSs (GAI > 1) in warm and dry areas, and EFSs were greener than PFSs (GAI < 1) in cold and wet areas. We also detected a stronger greening trend (0.0040 vs 0.0034 y−1) and a higher sensitivity of NDVI to temperature (0.031 vs 0.026 °C−1) on PFSs than EFSs, leading to a significant positive trend in GAI (0.00065 y−1, P < 0.01) in the TP from 1991 to 2020. Our results suggest that global warming exacerbated the greenness asymmetry associated with the slope orientation: PFSs are more sensitive to warming and have been greening at a faster rate than EFSs. The gradient of EFSs and PFSs provided a “natural laboratory” to study interaction of water and temperature limitations on vegetation growth. Our study is the first to detect the effect of aspect on the greening trend in the TP. Future research needs to clarify the full biotic and abiotic determinants for this spatial and temporal asymmetry of greenness across aspects with the support of extensive field measurements and refined high-resolution NDVI products.

Temporal and spatial analysis of vegetation cover change in the Yellow River Delta based on Landsat and MODIS time series data.

Based on the Landsat normalized difference vegetation index (NDVI) and the NDVI product of MODIS, this study synthesized two kinds of time-series images. The features were selected according to the characteristics of the time series, and the random forest algorithm was used for classification. Based on the classification results and GIS spatial analysis, the temporal and spatial changes in vegetation cover in the Yellow River Delta from 2000 to 2020 were studied. The results showed that from 2000 to 2020, the vegetation first increased and then decreased, and the dynamic degree of land cover change was generally low. The monthly average minimum NDVI values during the vegetation growth period mostly occurred before 2010, and the maximum values occurred after 2010. From the spatial perspective, the average vegetation area of the Yellow River Delta accounted for 31.54% of the total study area; specifically, the spatial pattern of vegetation distribution was relatively fixed, and the fixed vegetation area accounted for 63.90% of the total vegetation area. The spatial distribution had significant differences, and the vegetation was distributed radially from the center of the Yellow River to the periphery, with significant fragmentation found outside the watershed. The Yellow River had a strong interference with vegetation growth, and the stable vegetation distribution areas were concentrated near the Yellow River. The correlation coefficient between vegetation distribution and the location of the Yellow River was - 0.9964.

Reconstruction of a Monthly 1 km NDVI Time Series Product in China Using Random Forest Methodology

The normalized difference vegetation index (NDVI) is one of the most common metrics used to describe vegetation dynamics. Unfortunately, low-quality pixels resulting from contamination (by features including clouds, snow, aerosols, and mixed factors) have impeded NDVI products’ widespread application. Researchers have thought of several ways to improve NDVI quality when contamination occurs. However, most of these algorithms are based on the noise-negative deviation principle, which aligns low-value NDVI products to an upper line but ignores cases where absolute surface values are low. Consequently, to fill in these research gaps, in this article, we use the random forest model to produce a set of high-quality NDVI products to represent actual surface characteristics more accurately and naturally. Climate and geographical products are used as model inputs to describe environmental factors. They represent the random forest (RF) model that establishes relationships between MODIS NDVI products and meteorological products in high-quality areas. In addition, auxiliary data and empirical knowledge are employed to meet filling requirements. Notably, the random forest (RF) algorithm exhibits a mean absolute error (MAE) of 0.024 and a root mean squared error (RMSE) of 0.034, in addition to a coefficient of determination (R2) value of 0.974. Furthermore, the MAE and RMSE of the RF-based method decreased by 0.014 and 0.019, respectively, when compared to those of the STSG (spatial–temporal Savitzky–Golay) plan and by 0.013 and 0.015, respectively, when compared to the LSTM (long short-term memory) method. R2 increased by 0.039 and 0.027, respectively, compared to the STSG and LSTM methods. We introduced a novel series of NDVI products that demonstrated consistent spatial and temporal connectivity. The novel product exhibits enhanced adaptability to intricate environmental conditions and promises the potential for utilization in investigating vegetation dynamics within the Chinese region.

Monitoring vegetation dynamics (2010–2020) in Shengnongjia Forestry District with cloud-removed MODIS NDVI series by a spatio-temporal reconstruction method

Shengnongjia Forestry District is the national natural reserve with abundant biological resources in China. Monitoring its variations of vegetation for the intimate connection with eco-environmental changes is of great significance. In this paper, the 16-day composite MODIS normalized difference vegetation index (NDVI) products (MOD13A1) with 500 m resolution from 2011 to 2020 were selected to investigate the vegetation dynamics in Shengnongjia Forestry District. To alleviate the cloud contamination in NDVI products, a spatio-temporal prefill method with harmonic analysis of time series (ST-HANTS) is proposed. ST-HANTS first prefills the raw NDVI time series using spatio-temporal information to reduce data gaps, which effectively improves the reconstruction performance of the subsequent HANTS algorithm. In the simulation experiments, ST-HANTS shows the highest average correlation coefficient and the lowest average root mean square error compared with other mainstream methods, including HANTS, Savitzky–Golay filter, wavelet transform, and data assimilation. The reconstruction curves are close to the upper envelope of the NDVI time series, which is more consistent with the vegetation phenology and can effectively capture the key points in the growth cycle. By analyzing the cloud-free NDVI time series reconstructed with ST-HANTS, results reveal the overall trend of Shennongjia vegetation coverage is high in the middle while low at the edge. Except for the population centers and Hongping Airport, the NDVI of most areas is above 0.7 and shows a remarkable increasing tendency. Moreover, the fluctuation degree of NDVI in the whole study area is low, indicating that the ecological environment of Shennongjia is relatively stable. The vegetation variations are influenced by land surface temperature and precipitation, and the vegetation growth response to precipitation exhibits an apparent hysteresis.

The response and sensitivity of global vegetation to water stress: A comparison of different satellite-based NDVI products

Normalized Difference Vegetation Index (NDVI) is an effective proxy for vegetation activity and terrestrial productivity, and the widely used global NDVI products include the third generation Global Inventory Modelling and Mapping Studies Advanced Very High-Resolution Radiometer NDVI (NDVI3g), Satellite Pour l’Observation de la Terre VEGETATION NDVI (NDVIsp) and Terra Moderate Resolution Imaging Spectroradiometer NDVI (NDVItr). Here, for the first time, we investigated the differences between the three NDVIs in response and sensitivity to water stress represented by root zone soil moisture, vapor pressure deficit, self-calibrating Palmer Drought Severity Index and Standardized Precipitation Evapotranspiration Index across climate zones, humidity gradients, vegetation types and tree cover gradients. Over the temperate and arid climates region, the three NDVIs show higher consistent response fractions (CRFs) to water stress in the water-limited regions, with their sensitivity being in the order of NDVI3g < NDVIsp < NDVItr. The three NDVIs tend to exhibit inconsistent responses to water stress in the areas with higher humidity and tree cover, and their sensitivities to water stress fall as humidity and tree cover increase. In addition, the three NDVIs display higher CRFs and sensitivities to water stress over the non-forested areas. Our results also show that, in the cold climates region, the three NDVIs respond in the opposite way to water stress as that in water-limited areas, and NDVI3g has higher sensitivity to water stress than NDVIsp and NDVItr. This study highlights the differences between three NDVIs products for qualifying and quantifying their response and sensitivity to water stress.

AVHRR NDVI Compositing Method Comparison and Generation of Multi-Decadal Time Series—A TIMELINE Thematic Processor

Remote sensing image composites are crucial for a wide range of remote sensing applications, such as multi-decadal time series analysis. The Advanced Very High Resolution Radiometer (AVHRR) instrument has provided daily data since the early 1980s at a spatial resolution of 1 km, allowing analyses of climate change-related environmental processes. For monitoring vegetation conditions, the Normalized Difference Vegetation Index (NDVI) is the most widely used metric. However, to actually enable such analyses, a consistent NDVI time series over the AVHRR time-span needs to be created. In this context, the aim of this study is to thoroughly assess the effect of different compositing procedures on AVHRR NDVI composites, as no standard procedure has been established. Thirteen different compositing methods have been implemented; daily, decadal, and monthly composites over Europe and Northern Africa have been calculated for the year 2007, and the resulting data sets have been thoroughly evaluated according to six criteria. The median approach was selected as the best-performing compositing algorithm considering all the investigated aspects. However, the combination of the NDVI value and viewing and illumination angles as the criteria for the best-pixel selection proved to be a promising approach, too. The generated NDVI time series, currently ranging from 1981–2018, shows a consistent behavior and close agreement to the standard MODIS NDVI product. The conducted analyses demonstrate the strong influence of compositing procedures on the resulting AVHRR NDVI composites.

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.

Mapping of Land Degradation Vulnerability in the Semi-Arid Watershed of Rajasthan, India

Global soils are under extreme pressure from various threats due to population expansion, economic development, and climate change. Mapping of land degradation vulnerability (LDV) using geospatial techniques play a significant role and has great importance, especially in semi-arid climates for the management of natural resources in a sustainable manner. The present study was conducted to assess the spatial distribution of land degradation hotspots based on some important parameters such as land use/land cover (LULC), Normalized Difference Vegetation Index (NDVI), terrain characteristics (Topographic Wetness Index and Multi-Resolution Index of Valley Bottom Flatness), climatic parameters (land surface temperature and mean annual rainfall), and pedological attributes (soil texture and soil organic carbon) by using Analytical Hierarchical Process (AHP) and GIS techniques in the semi-arid region of the Bundi district, Rajasthan, India. Land surface temperature (LST) and NDVI products were derived from time-series Moderate-Resolution Imaging Spectroradiometer (MODIS) datasets, rainfall data products from Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), terrain characteristics from Shuttle Radar Topography Mission (SRTM), LULC from Landsat 9, and pedological variables from legacy soil datasets. Weights derived for thematic layers from the AHP in the studied area were as follows: LULC (0.38) &gt; NDVI (0.23) &gt; ST (0.15) &gt; LST (0.08) &gt; TWI (0.06) &gt; MAR (0.05) &gt; SOC (0.03) &gt; MRVBF (0.02). The consistency ratio (CR) for all studied parameters was &lt;0.10, indicating the high accuracy of the AHP. The results show that about 20.52% and 23.54% of study area was under moderate and high to very high vulnerability of land degradation, respectively. Validation of LDV zones with the help of ultra-high-resolution Google Earth imageries indicates good agreement with the model outputs. The research aids in a better understanding of the influence of land degradation on long-term land management and development at the watershed level.

A Global 250-m Downscaled NDVI Product from 1982 to 2018

Satellite-based normalized difference vegetation index (NDVI) time series data are useful for monitoring the changes in vegetation ecosystems in the context of global climate change. However, most of the current NDVI products cannot effectively reconcile high spatial resolution and continuous observations in time. Here, to produce a global-scale, long-term, and high-resolution NDVI database, we developed a simple and new data downscaling approach. The downscaling algorithm considers the pixel-wise ratios of the coefficient of variation (CV) between the coarse- and fine-resolution NDVI data and relative changes in the NDVI against a baseline period. The algorithm successfully created a worldwide monthly NDVI database with 250 m resolution from 1982 to 2018 by translating the fine spatial information from MODIS (Moderate-resolution Imaging Spectroradiometer) data and the long-term temporal information from AVHRR (Advanced Very High Resolution Radiometer) data. We employed the evaluation indices of root mean square error (RMSE), mean absolute error (MAE), and Pearson’s correlation coefficient (Pearson’s R) to assess the accuracy of the downscaled data against the MODIS NDVI. Both the RMSE and MAE values at the regional and global scales are typically between 0 and 0.2, whereas the Pearson’s R values are mostly above 0.7, which implies that the downscaled NDVI product is similar to the MODIS NDVI product. We then used the downscaled data to monitor the NDVI changes in different plant types and places with significant vegetation heterogeneity, as well as to investigate global vegetation trends over the last four decades. The Google Earth Engine platform was used for all the data downscaling processes, and here we provide a code for users to easily acquire data corresponding to any part of the world. The downscaled global-scale NDVI time series has high potential for the monitoring of the long-term temporal and spatial dynamics of terrestrial ecosystems under changing environments.

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 (&gt;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.

New Investigation and Challenge for Spatiotemporal Drought Monitoring Using Bottom-Up Precipitation Dataset (SM2RAIN-ASCAT) and NDVI in Moroccan Arid and Semi-Arid Rangelands

Abstract Remotely sensed soil moisture products showed sensitivity to vegetation cover density and soil typology at regional dryland level. In these regions, drought monitoring is significantly performed using soil moisture index and rainfall data. Recently, rainfall and soil moisture observations have increasingly become available. This has hampered scientific progress as regards characterization of land surface processes not just in meteorology. The purpose of this study was to investigate the relationship between a newly developed precipitation dataset, SM2RAIN (Advanced SCATterometer (SM2RAIN-ASCAT), and NDVI (eMODIS-TERRA) in monitoring drought events over diverse rangeland regions of Morocco. Results indicated that the highest polynomial correlation coefficient and the lowest root mean square error (RMSE) between SM2RAIN-ASCAT and NDVI were found in a 10-year period from 2007 to 2017 in all rangelands (R = 0.81; RMSE = 0.05). This relationship was strong for degraded rangeland, where there were strong positive correlation coefficients for NDVI and SM2RAIN (R = 0.99). High correlations were found for sparse and moderate correlations for shrub rangeland (R = 0.82 and 0.61, respectively). The anomalies maps showed a very good similarity between SM2RAIN and Normalized Difference Vegetation Index (NDVI) data. The results revealed that the SM2RAIN-ASCAT and NDVI product could accurately predict drought events in arid and semi-arid rangelands.

Short-term agricultural drought prediction based on D-vine copula quantile regression in snow-free unfrozen surface area, China

This study aims to propose an agricultural drought prediction model based on D-vine copula quantile regression considering several factors related to agricultural drought occurrence mechanisms in China. The proposed model is applied for short-term agricultural drought prediction (represented by 1-month time scale of the standardized soil moisture index (SSI-1)) considered antecedent soil moisture (represented by SSI-3), real-time rainfall (represented by 1-month time scale of the Standardized Precipitation Index (SPI-1)) and vegetation cover (represented by the Normalized Difference Vegetation Index (NDVI)) based on monthly soil moisture data from the Climate Change Initiative (CCI) program of European Space Agency (ESA), precipitation data from the Tropical Rainfall Measuring Mission (TRMM) and NDVI products from MODIS. The proposed model was employed and evaluated in China and results showed it performed well in snow-free unfrozen surface area such as south-east China. The outcome of this study can contribute to early warning for agricultural drought.