Linear and Non-Linear Vegetation Trend Analysis throughout Iran Using Two Decades of MODIS NDVI Imagery

Vegetation is the main component of the terrestrial Earth, and it plays an imperative role in carbon cycle regulation and surface water/energy exchange/balance. The coupled effects of climate change and anthropogenic forcing have undoubtfully impacted the vegetation cover in linear/non-linear manners. Considering the essential benefits of vegetation to the environment, it is vital to investigate the vegetation dynamics through spatially and temporally consistent workflows. In this regard, remote sensing, especially Normalized Difference Vegetation Index (NDVI), has offered a reliable data source for vegetation monitoring and trend analysis. In this paper, two decades (2000 to 2020) of Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI datasets (MOD13Q1) were used for vegetation trend analysis throughout Iran. First, the per-pixel annual NDVI dataset was prepared using the Google Earth Engine (GEE) by averaging all available NDVI values within the growing season and was then fed into the PolyTrend algorithm for linear/non-linear trend identification. In total, nearly 14 million pixels (44% of Iran) were subjected to trend analysis, and the results indicated a higher rate of greening than browning across the country. Regarding the trend types, linear was the dominant trend type with 14%, followed by concealed (11%), cubic (8%), and quadratic (2%), while 9% of the vegetation area remained stable (no trend). Both positive and negative directions were observed in all trend types, with the slope magnitudes ranging between −0.048 and 0.047 (NDVI units) per year. Later, precipitation and land cover datasets were employed to further investigate the vegetation dynamics. The correlation coefficient between precipitation and vegetation (NDVI) was 0.54 based on all corresponding observations (n = 1785). The comparison between vegetation and precipitation trends revealed matched trend directions in 60% of cases, suggesting the potential impact of precipitation dynamics on vegetation covers. Further incorporation of land cover data showed that grassland areas experienced significant dynamics with the highest proportion compared to other vegetation land cover types. Moreover, forest and cropland had the highest positive and negative trend direction proportions. Finally, independent (from trend analysis) sources were used to examine the vegetation dynamics (greening/browning) from other perspectives, confirming Iran’s greening process and agreeing with the trend analysis results. It is believed that the results could support achieving Sustainable Development Goals (SDGs) by serving as an initial stage study for establishing conservation and restoration practice