RESUMEN
Woody vegetation plays a vital role in regulating the water budget and energy exchange in the Earth's system. This study aimed at analyzing the spatiotemporal variability of Normalized Difference Vegetation Index (NDVI) and its response to Potential Evapotranspiration (PET), rainfall (RF), soil moisture (SM), and temperature (TEM) in the study area. The trends, correlations, and relationships between NDVI and climate variables were executed using Mann-Kendall monotonic trend (MKMT), partial correlation coefficients (PCC), and multiple linear regression (MLR) methods, respectively. Over the last 26 years, the interannual NDVI increased by 0.0065 yr-1 (R2 = 0.159, p = 0.157). The spatiotemporal MKMT and Theil-Sen slope analysis showed that interannual NDVI increased significantly in 78% of the basin's total area. Of the 78% of the basin, 31%, and 47%, of the total area showed extremely significant increasing (Zmk = 4.706, p ≤ 0.01), and significant increasing trends (Zmk = 2.378, p ≤ 0.05) respectively. The interannual variation of NDVI was well explained (R2 = 0.88, Adjusted R2 = 0.84) by the climate variables in the eastern, southeastern, and central sub-basins where agriculture, grass, sparse vegetation and barelands are the predominant land use land cover (LULC) classes. The main climatic factors that control vegetation growth and greenness during the rainy season were found to be PET, SM, and RF with 0.91, 0.99, and 0.86 PCC with NDVI respectively. The current study broadens the scientific community's understanding of the relationship between climate variables and vegetation growth in highland ecosystems. Understanding the seasonal and long-term relationship between climate and NDVI contributes to the scientific knowledge of highland ecosystems, which are extremely vulnerable to climate change.