Spatiotemporal analysis of vegetation cover changes around surface water based on NDVI: a case study in Korama basin, Southern Zinder, Niger

In the modern era, vegetation dynamics is an important aspect of climate change studies. The present study examined spatiotemporal changes of (NDVI) normalized difference vegetation index in the Korama basin (Southern Zinder of Niger) from 2000 to 2018, and their correlation with climatic factors was predicted. To analyze the change of vegetation cover, geographical information system, MODIS_NDVI, remote sensing, and climate variables (e.g., temperature and precipitation) datasets were used. Further, the correlation was performed for different years of vegetation types during the growing season (June–October). Our results show an increasing trend in average maximum annual NDVI across the Korama River Basin in the years 2000 and 2018. Conversely, significantly increasing trends in most of the areas were reported. Moreover, in downstream the vegetation cover is increased in Matameye and Magaria, but with a smaller increase in the upstream rate in Mirriah. Furthermore, a decrease in the surface water was observed in the Tessaoua, Matameye, and Magaria sections of the study region in 2000 and 2018, while a rise in water surface area was observed in Matameye and Magaria in the years 2006 and 2012. During rainy and dry seasons, NDVI correlated differently with temperature and precipitation with strong seasonal variations, while the mean vegetation period of NDVI does not show any significant change. In addition, moderate increase was observed in years 2000 and 2012 (r: 0.22; P: 0.50; R2: 0.05; r: 0.31; P: 0.34, R2: 0.10, respectively), and weak decrease in 2006 and 2018 (r: 0.61; P: 0.04; R2: 0.37; r: 0.58; P: 0.06, R2:0.33, respectively). The analysis indicates that climatic parameters such as precipitation and temperature are the main limiting factors affecting the vegetation growth. Indeed, the trends calculated by the correlation analysis showed that as climate factors increased (July, August, and September), the NDVI value increased at a rate of 0.16, reflecting the best growth in vegetation and rise in water bodies, although significantly decreased during years. This study would be highly useful in choice-making for sustainable water resource management in the Korama watershed in Southern Zinder, Niger.


Introduction
Vegetation is the key component of the earth's terrestrial ecosystems and plays a major role in soil, energy exchange, and biogeochemical cycles on earth's surface (Liu et al. 2018;Ning et al. 2015). Vegetation plays an irreplaceable function in preserving climate stability, controlling the carbon balance, and rising global greenhouse gases (GHGs), connecting soil, atmosphere, and moisture (Chuai et al. 2013;Kumar et al. 2017;Kumar and Kumar 2020). However, vegetation is prone to climate change. The dynamics of vegetation and their responses to climate change is recognized as one of the main issues of environmental change in the terrestrial ecosystem. With progress in remote sensing technologies, more information can be obtained from multispectral and multi-dated remote sensing datasets that provide advanced methods for investigating vegetation growth and inter-annual and seasonal variations (Ding et al. 2007).
The spatiotemporal variability of vegetation coverage is of great importance in the future exploring the internal relationship between vegetation and climate change, exposing the process of evolution of regional environmental conditions and forecasting future trends in growth Chauhan et al. 2020). Vegetation is the natural connection of the earth to the soil, atmosphere, and humidity. It also indicates clear signs of seasonality and inter-annual change (Chuai et al. 2013;Cui and Shi 2010) and is a touchy indicator of worldwide climate change (Weiss et al. 2004); vegetation responds simply and subtly to weather change. Nowadays, diminished soil conditions in the area affect particularly vegetation growth, in which the water system, as well as vegetation, remains vulnerable (Froend and Sommer 2010). Water resources respond to the climatological situation relatively to groundwater, providing significant protection toward climatic variability, particularly depression (Calow et al. 2010). Rouse and RH (1974) proposed the standardized difference vegetation index (SDVI), based on variations in pigment absorption characteristics in the electromagnetic spectrum's red and near-infrared regions. The complexities of vegetation have enhanced our view for vegetation seasonal growing through remotely sensed data (Meng et al. 2019). Thus, an ideal indicator of vegetation growing has called (NDVI) the normalized vegetation difference index (Piao et al. 2015). NDVI values range from − 0.1 to 1.0, that positive NDVI values suggest growing quantities of green vegetation. NDVI values close to zero and decreasing negative values suggest unvegetated characteristics such as rock and soil surfaces and water, air, ice, and clouds (Chuai et al. 2013;Schnur et al. 2010). Vegetation canopy analyzed with NDVI and remotely sensed data expose efficiently vegetation in the NIR (near-infrared) wavelengths with higher chlorophyll in the RED (red) wavelengths (Zhi-qiang and Dennis 2001). The vegetation index can largely reflect vegetation dynamics. For more details, the mathematical representation of NDVI is given as follows: NDVI = (NIR-Red)/(NIR + Red), where NIR represented reflectance radiated in near-infrared wavebands and Red in the visible reflectance wavebands from the satellite. Further, it is the primary indicator of the growth and spatial pattern of vegetation (Liu et al. 1999;Purevdorj et al. 1998), as well as plant morphology (Defries and Townshend 1994;Derrien et al. 1992). Besides, weather variables (e.g., temperature and precipitation) influence vegetation growing mechanisms (Li and Tao 2000;Li 1998;Nicholson et al. 1990;Schmidt and Gitelson 2000;Zhang et al. 2003). Land-use changes and climate factors can have several effects on vegetation, including temperature and precipitation are the main indicators used to characterize climate conditions and can have a clear impact on vegetation growth (Chuai et al. 2013). Several scientists have long used the NDVI, temperature, precipitation, and other data to study the effect of climate change on vegetation. Over the last three decades, NDVI has been widely used for analyzing climatic effects on modifications on productivity of in plant life (Chuai et al. 2013;Hu et al. 2011;Nemani 2003), whereas some recent studies examined hydrological responses to vegetation cover changes. Jiao et al. (2017) analyzed the effect of vegetation dynamics on hydrological processes in a semiarid basin using a land surface hydrology coupled model and found that vegetation dynamics induced higher ET, lower runoff, and soil moisture content and decelerated the evolving trend of discharge (Zhang et al. 2019;Wu et al. (2020).
The relationship between NDVI, precipitation, and global surface temperature was studied by Schultz and Halpert (1995), and no substantial connection was found between them. Ni et al. (2006) studied the variability of vegetation and geographic characteristic. Other research observed a great relationship between these variables and recorded changes in NDVI such as (Guo et al. 2008), which were significantly associated with both temperature and precipitation. However, some of them found a positive correlation between NDVI and precipitation and reported a strong positive correlation between NDVI and temperature at high latitude (Ichii et al. 2002). Rasmussen (1998) reported that water scarcity is the key to terrestrial ecosystem efficiency and sustainability, and changes in hydrology could have large effects on the characteristics of vegetation as well as soil. Enhanced human supply for water sources (Vorosmarty et al. 2005) may have significantly changed water supply as just a global, such as the impact of climatic change (Bartram and Cairncross 2010). The rising of CO 2 and climate variation was projected to dramatically change the global hydrological condition and has a range of effects on water source systems (Wu et al. 2012;Kumar et al. 2019). Further, the hydrological systems remain influenced by vegetation mechanisms like canopy acquisition, including water absorption from the plant (Zhao et al. 2018a). Therefore, it should be noted that the availability of water leads to good vegetation dependency, which functions are the main key by water availability (Zhao et al. 2018b).
Literature reveals that few studies have been made on water supplies and the conservation of surface water, therefore excluded from freshwater supply assessments (Calow et al. 2010). Furthermore, vegetation is the natural relation of soil, atmosphere, and humidity on Earth. It shows strong seasonal and annual shifts (Chuai et al. 2013;Cui and Shi 2010) and serves as a responsive predictor of global climate change. The change in vegetation and surface water is widely symbolic of global analysis variation in the study area where hydrological processes could better know the natural impacts and manmade disruptions on meteorological characteristics as well as predicted trends in water management (Zhao et al. 2018b). Using MODIS_NDVI satellite data, a vegetation map of Korama basin, temperature, and precipitation data, this study analyzed changes in NDVI, temperature, and precipitation and investigated the correlation between NDVI, temperature, and precipitation for different types of vegetation across 19-year wet and dry growing seasons in Korama watershed. The objectives of this study were framed to understand the spatial change of vegetation around surface water in a changing environment in the Korama watershed and to compare the correlation between NDVI, temperature, and precipitation for different vegetation types and then to clarify the reasons for the various formed relationships. Furthermore, differences between climatic variables and MODIS_NDVI satellite data throughout 2000-2018 have also been made.

Study area
The study area or basin of the Korama having a catchment area of 120,000 km 2 corresponds to the extreme west of the great basin of Lake Chad and is located on the borders of Niger and Nigeria. The basin is located from east to west and limited to north and east by outcrops of the Precambrian basement of Damagaram Mounio, to the west through sedimentary deposits of the basin of Iullemenden, and to the south by way of the border of Niger-Nigeria (Fig.1). With a The climate in the region is the sort of Sahel with a short rainy season (June through September) and a long dry season (October through May). The annual precipitation amounts increase (500-600 mm) from north to south, with an average over the period 1930-2011, from 500 mm to the synoptic station of Zinder airport, north of the basin, and over 600 mm in Magaria climate station, south of the basin. Korama River and its tributaries are part of the properly established hydrographic community within the endorheic basin. Such streams are only represented in areas, even in the rainy season, by strings of ponds (Sandao 2010). The geological map of Niger State basement complex and sedimentary basins has been in the literature (Amadi et al. 2012).

Data description
The current study used remote sensing and satellite data to analyze vegetation and water body changes based on NDVI via spatial approach over the period from 2000 to 2018. A (DEM 90 m of resolution) was obtained from NASA satellite, (https ://earth explo rer.usgs.gov/). MODIS/TERRA MOD13Q1vegetation indices 16 days L3 Global 250 m, with high-resolution satellite data for the years 2000, 2006, 2012, and 2018, respectively, were obtained from NASA's earth observation system (https ://earth explo rer.usgs.gov/). The climate variable (e.g., temperature, precipitation) from 1900 to 2017 was obtained from CRU web satellite data, (https ://www.cru.uea.ac.uk/data). These datasets have been designed for a well appreciative of our research area analysis by using datasets from the years 2000, 2006, 2012, and 2018, respectively. Hydrometeorological data collected for this study are given in Table 1.

Approach for the spatial and temporal analysis of NDVI
The method to improve the reflectance radiated in the NIR near-infrared bands and the reflectance radiated in the red visible bands was used to synthesize monthly NDVI data based on MODIS data from the preprocessed time and obtain the annual NDVI dataset. The method of trend analysis was used to analyze patterns in the annual, interannual, and spatial NDVI variability. This approach employs a minimum square method to match the NDVI variance trend slope over a certain period of time, and it can be used to completely represent the spatiotemporal characteristics of changes in vegetation. Remote sensing and the various computational software and/or statistical tools, such as Excel, word, Arc-GIS, Origin Lab, were used for the satellite data from NASA Earth Observation System (EOS) for February 21, 2019, for extraction of hydrological features in the GIS system (Arc Map 10.3 software). The area of water was extracted from the attribute tables by calculating geometry from 2000 to 2018. Further, NDVI could be analyzed from imagery data and recognized as valid plants covering indicators (Olusegun and Adeyewa 2013; Hess et al. 1996). Buffer analysis within the GIS environment was accomplished to know the variability in the NDVI response and further estimate the temporal changes around the surface water. Thus, vegetation patterns were widely investigated due to variability in the water cycle and climatic parameters. The amount of NDVI from − 1 and 1, where after classifying a greater value, indicates greater vegetation, whereas low value generally means less vegetation cover (Stow al. 2003). The watershed was delineated from DEM and further created a separate basin (Fig. 2).
To be consistent in this study, MODIS/TERRA MOD13Q1vegetation indices 16 days L3 Global 250 m, of the year 2000, 2006, 2012, and 2018 would be utilized to analyze the correlation between climatic parameters and vegetation growing. The correlations coefficient and band x mean value of NDVI in the year, y i monthly climatological parameter, − y mean value of climatological parameters in the year, n = 12, r xy mean correlation coefficient among climatic variables and MODIS_NDVI. NDVI helps to examine growth and greenness of vegetation in a particular region. It is the ratio of the difference in reflectance near-infrared bands and reflectance in the visible bands to their sum (Rouse and RH 1974;Weiss et al. 2004;Kumar et al. 2017), and the mathematical representation is given in Eq. 2.
where NIR is the reflectance in the near-infrared band and RED the reflectance in the visible band. Its value ranges + 1 to − 1, positive values represent healthy vegetation cover, and negative values represent poor vegetation cover ). Vegetation cover is calculated through the NDVI value mapped in the studied catchment. The NDVI usually varies around 0.1 and 0.7 in heavy vegetation fields, as different to green vegetation. Moreover, water has more visible reflectance, that NIR arising in a negative value of NDVI, whereas some areas such as bare soil and rock provide the equal value of NDVI in both bands and close to zero (Rouse and RH 1974). In Equation 3, the interaction among NDVI and FVC scaled was shown (Carlson and Repley 1997): where FCV (fractional vegetation cover) represents a reasonable approximation of vegetation cover area pixel percentage and generally represents a good indicator for vegetation cover estimation. Thus, fractional vegetation cover varies from zero to one.
Generally, MODIS/TERRA MOD13Q1vegetation indices 16 days L3 Global 250 m provides an interesting view of variation across ecosystems with higher spatial frequency, thereby providing a strong benchmark for future change assessment.

Estimation of change detection in NDVI (2000-2018)
The NDVI was used to determine green vegetation production and for detecting changes in vegetation. It is also widely used to calculate the interaction among atmospheric variation and the level of vegetation growing. After classifying high and low vegetation, a pixel with high NDVI values indicates high vegetation, whereas low NDVI value generally means less vegetation (Fig. 3). Furthermore, negative NDVI values are the best factor for water bodies. Hence, the annual variation in the growing season NDVI, precipitation, and temperature for different types of vegetation is illustrated in Fig 3. The mean NDVI vegetation period has not substantially changed. But a moderate increase was observed in the years 2000 and 2012 (r: 0.226; P: 0.503; R 2 : 0.051; r: 0.317; P: 0.341, R 2 : 0.100, respectively) and substantially decrease in 2006 and 2018 (r: 0.610; P: 0.046; R 2 : 0.372; r: 0.581; P: 0.060, R 2 :0.337, respectively). Further, Fig. 3 shows changes in the growing season (NDVI from 2000 to 2018 in Korama basin, Zinder, Niger). The investigations show that most of the water sources have reduced in their aerial extent and some of them have even disappeared. Significantly, surface water areas have been reducing from general to particular in the study area. It is due to the intensification and expansion of the agricultural area, urbanization, and agricultural land use of chemicals (pesticides), which causes additional vegetation reduction.   Fig. 4. NDVI tends to be greatly affected by temperature in these areas at the start of the growing season. Areas of negative correlations with NDVI temperature correlated well with areas with positive NDVI-precipitation correlations, often in semiarid regions. In these regions, inter-annual variability of NDVI was determined by temperature as well as precipitation. As the correlation of NDVI-precipitation is higher than that of NDVI-temperature, precipitation is considered to have the dominant effect on the growth of vegetation in the semiarid area. Figure 4 shows that, for all vegetation types, NDVI correlated positively with temperature in the Magaria and Tessaoua parts of the study area. NDVI is also positively correlated with precipitation for all parts of vegetation cover except for the Southern part of the study area (i.e., Mirriah). However, these results suggest that the warmer environment and increased precipitation in Korama basin are critical for most vegetation cover areas to grow in the wet season. Furthermore, higher temperatures and low precipitation are advantageous for the growth of vegetation in the wet season, as they are situated under conditions of ample precipitation, and more precipitation will lower temperature, as described above. However, it is more complex to measure the correlation with surface water and climatic parameters variable. So in this work, from each NDVI values, the change analysis is an important way to define the variations detected from 2000 to 2018. This research suggests how varying patterns and procedures may be detected at more than one scale the usage of alternate in NDVI values.

Interannual changes of vegetation (2000-2015)
Water deficit is positively correlated NDVI as water-deficit stress occurs, the soil becomes green (Fig. 4), and NDVI is decreasing, followed by vegetation undergoes a regression. Furthermore, interannual trends during the period of 2000-2018 were analyzed based on the climate correlation maps obtained in this study. Results show the rise in NDVI characterized by a positive NDVI temperature dependence. Thus, temperature is one of the most significant factors influencing NDVI. In the years of 2006 and 2018, the decline in precipitation is correlated with a decrease in NDVI, as precipitation plays a vital role in plant production in these semiarid regions. Although there was a strong coefficient of correlation between NDVI and precipitation in the area because it is known that the analysis of detection changes of vegetation cover is due to the availability of water in the soil. However, some variation of NDVI values in the most study area is widely shown during 2000-2018.
In all years from 2000 to 2018, there was a respectable vegetation cover around a water body in Tessaoua, Matameye, and Magaria parts of the study area, while a decrease in vegetation was observed in Mirriah part of the years 2000 and 2006. This might be due to a good tendency of precipitation with a decrease in temperature in the parts with good vegetation coverage and a decrease in precipitation with higher temperature during the time of vegetation decrease.
The correlation was explored between the interannual differences throughout the NDVI as well as the climate variables. It has been shown that the interannual temperature variability, especially in the increasing season, has a significant effect on ecosystem activity in the high latitude area. The growth of vegetation in semiarid regions is mainly restricted through precipitation, although water resource, as the primary cause during plant growth, is often not directly linked to the arid and semiarid areas. But for two reasons in a regular climate with abundant wet-season water followed by a noticeable dry season. First, throughout the rainy season, the groundwater is sufficiently charged (Adimalla et al. 2020). Secondly, related to higher rainfall and resource for higher production, an improved performance can restore and spread carbon through the process, enabling plants to have greater implications for deeper resource access (Liu et al 1999;Rana et al. 2020), although the good reaction of vegetation to a climatic parameter such as precipitation suggests that ecosystem growth and development are linked by precipitation.

Relationship between NDVI and climatic variables
There is a need to note a similarity in the seasonal of the NDVI index in the study area during the period 2000-2018. The correlation between annual average NDVI and climatic variables is shown in Fig. 3; the figure shows time-analyzed and distribution of the relationship between NDVI and climatic factors. A positive correlation between NDVI and temperature was found in Southern Niger during the wet period, while in the dry season negative correlations were detected. Figure 4 shows that seasonal NDVI is positively correlated with temperature in the Korama basin (Magaria and Tessaoua sections). During the wet season, the correlation was significantly positive, slowly weakening in the drying season. The correlation between NDVI and temperature during the dry season was slightly negative, especially in the summer months (Fig. 5). Similarly in the rainy season (August-October), the correlation was greater than in the summer. In summary, areas with positive correlations across the entire Korama basin were greater than with negative correlations. On the other hand, the NDVI-precipitation correlation is distinct from the NDVI-temperature correlation (Fig. 6). In southern Korama, the positive correlation is pronounced, which is distinctly different from the rest of the cold climate zone. Besides the systematic examination of the correlation, coefficient values, we noticed that the coefficient spatial correlation between NDVI and precipitation was very different for NDVI depending on the period of the growing season, the time interval, and the period during which precipitation was accumulated, and the form of land cover. Therefore, climate change can be attributed as a cause of changes in climatic parameters, especially temperature and precipitation. These parameters are vital factors for plant growth .
The NDVI levels start to rise in 0.11 and hit their peak values to 0.16 and then decline thereafter, representing the best vegetation growth in the area (Fig. 6), decrease for the NDVI values indicating less vegetation growth. Furthermore, Table 2 shows the difference of correlation coefficient between the Pearson coefficient and R-square, which were calculated, however, that the slope is substantially different to 0, while NDVI is already at 0.05 rate, which indicates the coefficient of Pearson correlation among climatic parameters and satellite data throughout the rising time of the year from 2000 to 2018. Besides, trough analysis by using the technique of change detection based on its characteristics such as vegetation in the area, hence the detection of vegetation cover variation with MODIS_NDVI values threshold of 0.73, 0.72, 0.65, and 0.69, respectively. The most precipitation, in the three months (July, August, and September), of cumulative rainfall has decreased, while NDVI values have significantly decreased during analysis (Fig. 6). Further, the correlation between climate variability and NDVI time's variable with 500-600 mm of mean precipitation in the area was good as shown in Fig. 7. Overall, the relationship with a certain precipitation threshold is roughly linear, but NDVI values change moderately as do annual precipitation. So the absolute values of precipitation were greater than the temperature in the wet season, which means, during growing season, NDVI is more susceptive with precipitation. Figure  3 shows the temporal variation of the NDVI trend from 2000 to 2018. (Fig. 8) The impacts of the climatic variable on plant growth, during 2000-2018 were analyzed. In this regard, the interaction between climate factors and NDVI was assessed. The coefficient of correlation was carried out for climate lead time for precipitation and temperature. Since a positive NDVI-climate variable correlation was observed as shown in Fig. 6, in lead time duration, the correlation was conducted throughout the climate variable and NDVI time series. Analysis results found that the variability in vegetation growth is linked to climate variables (e.g., precipitation and temperature) and human activities (e.g., fertilization and irrigation). However, Fig. 7 shows the change of meteoroidal graphs from 1901 to 2017 (annual temperature ( °C ) and annual precipitation total (mm)), which was used to analyze the climate variability during a period 1901-2017 by using observed and CRU data.
The interaction between NDVI and its cumulative effects on the future and/or regional climate has been analyzed by applying band collection statistics and basic statistical parameters (e.g., minimum, maximum, mean, and standard deviation) for every layer. The results of the correlation between NDVI are given in Table 3. Table 3 shows an analysis of a set of NDVI raster calculations from 2000 to 2018 using the correlation matrix option. Further to proceed only the basic statistical measurements (mean, minimum, maximum, as well as standard deviation), we will be calculated for each layer. This is a default value, other than the standard, which calculates the dataset of NDVI raster bands from 2000 to 2018. However, the correlation matrix between each layer (NDVI) will also be determined. Furthermore, the correlation matrix results for bands collection statistics between NDVI 2000NDVI , 2006NDVI , 2012NDVI , and 2018, respectively, showed as layers: layers 1 and 2 (0.528), layers 1 and 3 (0.507), layers 1 and 4 (0.047), layers 2 and 3 (0.543), layers 2 and 4 (0.306), and layers 3 and 4 (0.421), respectively. NDVI option allows greater results for vegetation that vary in density and also for vegetation dispersed from spatial images of different time series, whereas the study suggests that variability of vegetation growth pattern is linked to climate variable (e.g., precipitation and temperature).

Impact on vegetation growth of surface water quality parameters in Korama Basin
Since 1990s, vegetation variability response to water soil has been an actively studied area of research at a regional and/or global scale in the context of changing climate and anthropogenic activities. The area under study climate conditions and topography is mainly affected by the dynamic changes in the vegetation cover. Moreover, climate change directly affects the growth of vegetation, while landscape conditions influence the distribution of water and heat at the surface. Vegetation and ecosystem services regulation in the watershed development and evolution is linked to the hydrological situation system  (Zipper et al. 2018). Sandao et al. (2010) conducted the study to see the impact of surface water quality parameters on vegetation growth based on physicochemical data (conductivity, salinity, residual sodium (SAR) and residual carbonates) and stated that water is poorly mineralized in the area and suitable for irrigation. Besides, they calculated the sodium absorption ratio (SAR) and completed a Table of Classification of water for irrigation, according to the SAR values and conductivity. Furthermore, they applied the method of residual sodium carbonate (RSC) and categorized the suitability of water for irrigation into three classes: (i) RSC below or equal to 1.25 indicates a good quality water for irrigation; (ii) RSC between 1.25 and 2.5 indicates average quality water for irrigation, and (iii) RSC above or equal to 2.5 indicates salty water and unsuitable for irrigation. The values of groundwater physical parameters, measured in situ (conductivity, pH and temperature), show that the values of groundwater conductivity in the Korama basin vary from 51 to 825 µs/ cm, with an average of 221 µs/cm. These values show that more than 95% of the water has conductivity below or equal to 750 µs/cm. This quality of water is coming within the good and excellent classes for irrigation. Groundwater in the area is generally acidic or rarely neutral (pH below or equal to 7 pH units). This parameter has not been specifically studied by the use of such water for irrigation. Further, the values of groundwater temperature, between 20 °C and 27 °C, are close to annual average values of the atmospheric temperature in the area. These temperatures do not seem to have side effects on the suitability of water for irrigation. In addition to water quality, climate variables and anthropogenic disturbance, NDVI may also be affected by soil degradation, air pollution, slope, and    and key factors affecting NDVI will be the main direction of our upcoming research.

Conclusion
In this study, the trend analysis method was used, based on temporal MODIS_NDVI data, to analyze the spatiotemporal variability in the vegetation cover throughout the Korama basin from 2000 to 2018. Based on meteorological and topographical evidence, the principal factors influencing the vegetation cover variations were investigated. Results reveal that the overall growth rate of NDVI in the study area across 19 years has been 0.16 in the downstream areas of the basin, while vegetation cover is found larger compared to the upstream areas. Similarly in the upstream areas, NDVI changes are small, with little fluctuation, whereas downstream changes are larger. The areas with increasing trends are distributed primarily in the downstream region, while the areas with declining trends are concentrated in the upstream portion of the basin. Areas with substantial increases in the spatial distribution of vegetation cover are primarily located near the section of Magaria, and major decreases in cover distribution are found near the central section and Mirriah. The NDVI shows a significant positive correlation with a climatic variable in the area under study. The associations between the NDVI and temperature and precipitation are variable in different regions. The precipitation in the upstream of the river greatly affects the forest areas, but meanwhile, the temperature in the upper basin reaches greatly affects the lower region. All the climate factors have little influence on the plain's agricultural areas, while they have a greater combined impact on the basin's central region with smaller individual impacts. Therefore, this study helps to understand the relationship between the NDVI, environment, and human factors and may offer suggestions for sustainable management in the Korama watershed. Further, the results of this study provide a reference for sustainable environmental management and response to vegetation change in the context of climate change.