Evaluation of Vegetation Indices and Phenological Metrics Using Time-Series MODIS Data for Monitoring Vegetation Change in Punjab, Pakistan

In arid and semi-arid regions, it is essential to monitor the spatiotemporal variability and dynamics of vegetation. Among other provinces of Pakistan, Punjab has produced a significant number of crops. Recently, Punjab, Pakistan, has been described as a global hotspot for extremes of climate change. In this study, the soil adjusted vegetation index (SAVI), normalized vegetation difference index (NDVI), and enhanced vegetation index (EVI) were comprehensively evaluated to monitor vegetation change in Punjab, Pakistan. The time-series MODIS (Moderate Resolution Imaging Spectroradiometer) data of different periods were used. The mean annual variability of the above vegetation indices (VIs) from 2000 to 2019 was evaluated and analyzed. For each type of vegetation, two phenological metrics (i.e., for the start of the season and end of the season) were calculated and compared. The spatio-temporal image analysis of the mean annual vegetation indices revealed similar patterns and varying vegetation conditions. In the forests and vegetation areas with sparse vegetation, the EVI showed high uncertainty. The phenological metrics of all vegetation indices were consistent for most types of vegetation. However, the NDVI result had the greatest variance between the start and end of season. The lowest annual VI variability was mainly observed in the southern part of the study area (less than 10% of the study area) based on the statistical analysis of spatial variability. The mean annual spatial variability of NDVI was <20%, SAVI was 30%, and EVI ranged between 10-20%. More than 40% of the variability was observed in the NDVI and SAVI vegetation indices.