Study on Relationship of Land Cover Changes and Ecohydrological Processes of the Tuul River Basin

: The Tuul River Basin is the most important socioeconomic and political base area of Mongolia. Therefore, studying the interrelationships between changes in the ecohydrological processes of this basin and its land cover is of great importance for maintaining sustainability and the environment. This study investigated the annual average air temperature, total annual precipitation, and river discharge variability, and land cover changes at selected stations of the basin by using the hydrometeorological analysis, satellite analysis, and land cover determination statistical analysis. During the study period, the average annual air temperature rose from − 1.5 ◦ C to +0.3 ◦ C (1.8 ◦ C). The average annual precipitation exhibits relatively low change during this period. River discharge varied during the study period. A signiﬁcant decreasing trend in river discharge was observed at the Terelj ( ϕ = − 2.72) and Ulaanbaatar ( ϕ = − 5.63) stations, whereas the other stations, Altanbulag, Lun, and Orkhontuul, showed a signiﬁcant increasing trend. During the study period, changes in land cover were directly related to main hydrometeorological parameters. Between 2000 and 2020, the amount of grassland decreased by 319.67 km 2 , while the area of water bodies increased by 28.36 km 2 . In the study area, mainly water bodies and sensitive areas of the land cover types were changed due to changes in precipitation. Studies in the arid and semiarid regions of Central Asia show that changes of ecohydrological processes have a signiﬁcant impact on land cover changes.


Materials
Hydrometeorology data such as Tuul River Basin air temperature, precipitation, and river discharge were obtained from the Information and Research Institute of Meteorology, Hydrology and Environment. In addition, some data were obtained from NOAA's National Centers for Environmental Information (NCEI) (https://ngdc.noaa.gov/) to verify the data of some stations.
The soil map, topographic map, ecological landscape potential map, vegetation map, and forest map were selected for accuracy testing validation. Landsat Enhanced and Thematic Mapper (ETM+) and Landsat Thematic Mapper (TM) satellite data with a resolution of 30 m obtained for years 2000 to 2020 were employed for land cover classification. A global land cover dataset (GlobeLand30) product provides satellite data for land cover in the study area, which is from 2014 and described by the National Geomatics Center of China. All Landsat TM and ETM + images derived from the USGS (http://landsat.usgs.gov/) were employed to interpret the spatial land cover situation during the previous 10 years.

Methods
Different trends in hydrometeorology variables  and land cover changes (2000, 2010, and 2020) were determined and relationships assessed using the innovative trend analysis method, Mann-Kendall test, Sen's slope estimator method, satellite imagery, and correlation calculations. To evaluate the reliability of the innovative trend analysis method (ITAM), the results were compared with the Mann-Kendall test (MK) and Sen's slope estimator test (β). ArcGIS software was used to identify and process the land cover type (LC type). Transformation matrix (TM) of land covers were used to identify differences in land cover changes and to detect changes in water-related components. The relationship between hydrometeorology variables and land cover was determined by correlation and regression analysis based on various measurement data and statistical calculations (SC) (Figure 2).

Hydrometeorology Methods
The study of long-term trends in water and climate change allows us to identify and improve its causes and effects [19,20]. The following popular trend analysis methods, innovative trend analysis method, Mann-Kendall test, and Sen's slope estimator method, were used to detect the various trends of climate and river discharge in the periods from 1979 to 2019. In the study, significance levels at 10%, 5%, and 1% were taken to assess the hydrometeorology times series data. In this study, the trend of changes in annual average in air temperature, precipitation, and river discharge data were analyzed separately by the following methods.
The innovative trend analysis method divides a time series into two equal parts and it sorts both subseries in ascending order [21,22]. Following that, the two halves are placed on a coordinate system : 1,2,3, … /2 on x-axis and : /2 1, /2 2, … on y-axis. If the time series data on a scattered plot are collected on a 1:1 (45°) straight line, it indicates no trend. On the other hand, the trend is increasing when data points accumulate above the 1:1 straight line and decreasing when data points accumulate below the 1:1 straight line. Mark the first half of the time series sorted along the ( ) axis and the second half of the time series along the ( ) axis (Figure 3).

Hydrometeorology Methods
The study of long-term trends in water and climate change allows us to identify and improve its causes and effects [19,20]. The following popular trend analysis methods, innovative trend analysis method, Mann-Kendall test, and Sen's slope estimator method, were used to detect the various trends of climate and river discharge in the periods from 1979 to 2019. In the study, significance levels at 10%, 5%, and 1% were taken to assess the hydrometeorology times series data. In this study, the trend of changes in annual average in air temperature, precipitation, and river discharge data were analyzed separately by the following methods.
The innovative trend analysis method divides a time series into two equal parts and it sorts both subseries in ascending order [21,22]. Following that, the two halves are placed on a coordinate system (x i : i = 1, 2, 3, . . . n/2) on x-axis and x j : j = n/2 + 1, n/2 + 2, . . . n on y-axis. If the time series data on a scattered plot are collected on a 1:1 (45 • ) straight line, it indicates no trend. On the other hand, the trend is increasing when data points accumulate above the 1:1 straight line and decreasing when data points accumulate below the 1:1 straight line. Mark the first half of the time series sorted along the (x j ) axis and the second half of the time series along the (y j ) axis ( Figure 3). The mean value difference among and give the trend magnitude of the data series. The ITAM indicator is expressed as: 1 10 The mean value difference among x i and x j give the trend magnitude of the data series. The ITAM indicator is expressed as: where ϕ is the trend indicator, n is the number of each subseries, and µ is the mean of the first data series. The value of ϕ is positive if it tends to increase. If the value is negative, it tends to decrease, and if i is between 1:1, it does not have a significant trend [23,24].
To evaluate the reliability of the innovative trend analysis method, the results were compared with the Mann-Kendall test and Sen's slope estimator test.
Mann-Kendall (MK) rank-based nonparametric statistical test applied in hydrometeorology studies is used for detecting trends in time series data [25]. It detects trends in a series without specifying whether the trend is linear or nonlinear. The variables S to be tested are calculated as: The trend test is applied to x i data values (i = 1, 2, . . . n − 1) and x j (j = i + 1, 2, . . . n). The value of each x i is used as a reference point to compare with the value of x j , which is given as: where x j and x i are the values in period j and i. When the number of data series is greater than or equal to ten (n ≥ 10), the MK test is then characterized by a normal distribution with the mean E(S) = 0 and variance Var(S) equated as [26,27]: where m is the number of the tied groups in the time series and t k is the number of ties in the kth tied group. The test statistic Z is: When Z is greater than zero, it indicates an increasing trend and when Z is less than zero, it is a decreasing trend [28].
The trend magnitude is calculated by Sen's slope estimator methods [21,23]. The slope Q i between two data points is given by the equation: where x j and x k are data points at time j and (j > k), respectively. When there is only a single datum in each time, then N = n(n−1)

2
; n is the number of time periods. However, if the number of data in each year is many, then N < n(n−1)

2
; n is the total number of observations. The N values of the slope estimator are arranged from smallest to largest. Then, the median of slope (β) is computed as: when N is odd Q[(N/2) + Q(N + 2)/(2)/(2)] when N is even (8)  The sign of β shows whether the trend is increasing or decreasing.

Land Cover Methods
Spatial-temporal changes and impact analysis were performed by duplicating maps to determine changes in land cover over time and to establish correlations between land cover patterns. The magnitude of the land cover change is determined by comparing the satellite image at the beginning (2000) and end (2020) of the study period [5,29]: where TA is total area, CA is changed area, CE is the extent of change, t 1 is beginning time, and t 2 is ending time. Identified changes in land cover were calculated, along with the area size of the types of land cover transferred. To evaluate its accuracy, the modified area was confirmed using the kappa coefficient. In other words, the kappa coefficient is a statistical technique that was applied in the present study for assessing the accuracy of land cover changes. It confirms the agreement between predefined producer ratings and user-assigned ratings [30]: where K is the kappa coefficient, P(A) is the number of times the K raters agree, P(E) is the number of times the K rates are expected to agree only by chance, A and D are unchanged categories, A1 and B1 are subject's categories, and N is the change of results.

Analysis and Estimation of Hydrometeorology and Land Cover Variables Influence in the Tuul River Basin
Determining the relationship between hydrometeorological variables and land cover in the study area is important for understanding the ecohydrological process. To do this, it is important to establish the relationship between the various parameters of the water for statistical calculations and to identify the interrelationship between the interactions. From the hydrometeorological variables, the interaction between precipitation, air temperature, river discharge, and water bodies of land cover is calculated by the correlation and regression relationship [31,32].

Analysis of Hydrometeorology
The mean annual air temperature, mean annual precipitation, and average annual river discharge were calculated for each station on the study site by determining the long-term conditions of hydrometeorology change. Additionally, the correlations were calculated using the mean value of hydrological variables.
The result shows that the ITAM curve annual air temperature reveals a statistically sharp increasing trend in Ulaanbaatar from 1979 to 2019 (ϕ = 14.11), a statistically abrupt increasing trend in Terelj from 1979 to 2019 (ϕ = 5.08), a statistically upward trend in Altanbulag from 1979 to 2019 (ϕ = 2.82), a statistically sharp upward trend in Ondorshireet from 1979 to 2019 (ϕ = 4.65), a statistically sharp upward trend in Lun from 1979 to 2019 (ϕ = 3.29), a statistically increasing trend in Bayannuur from 1979 to 2019 (ϕ = 4.25); in Orkhontuul, a statistically increasing trend was observed with (ϕ = 4.75) from 1979 to 2019. An overall statistically significant upward trend was observed in the average (five stations) (ϕ = 7.02) ( Figure 4). In the study period, the average annual air temperature rose from −1.5 • C to +0.3 • C (1.8 • C). This rise could have a direct and indirect significant effect on the ecohydrological processes of the arid and semiarid regions of Central Asia. (a) Terelj  The trend analysis of annual air temperature for all stations is presented in Table 1. The ITAM curve for mean annual precipitation reveals a statistically sharp downward trend in Terelj from 1979 to 2019 (ϕ = −2.44); in Ulaanbaatar from 1979 to 2019 a downward trend (ϕ = −0.49); a statistically increasing trend observed in Altanbulag (ϕ = 1.7) and Ondorshireet from 1979 to 2019 (ϕ = 1.17); a decreasing trend in Lun (ϕ = −0.27); and increasing trend was observed in Bayannuur from 1979 to 2019 (ϕ = 1.38); in Orkhontuul, an upward trend was obtained (ϕ = 0.96). Overall, a slightly increasing trend was observed at the average condition (ϕ = 0.60) ( Figure 5). In the study period, the mean annual precipitation was 238 mm and the change was not significant. The annual trend analysis of precipitation at all stations is presented in Table 2. From 1989 to 2019, the ITAM curve mean annual river streamflow shows a statistically sharp downward trend in Terelj (ϕ = −2.72), decreasing trend was observed in Ulaanbaatar (ϕ = −5.63), and increasing trend was obtained at the Altanbulag, Orkhontuul, and Lun stations: (ϕ = 5.34), (ϕ = 5.07), and (ϕ = 3.42), respectively. Overall, the average value for discharge shows a slightly increasing trend (ϕ = −0.92) ( Figure 6).  The annual trend analysis of river discharge for all stations is presented in Table 3.

Analysis of Land Cover
Determining land cover changes using satellite data is important because it allows us to determine land cover changes, shifts, and causes. Study areas were classified using two Sustainability 2021, 13, 1153 9 of 16 commonly used supervised classification methods. The maximum-likelihood classification of the Tuul River Basin was used to analyze land cover maps for three different years (2000, 2010, and 2020). To confirm this, there results were compared with the soil, topographic, ecological landscape potential, vegetation, and forest maps. With appreciable accuracy, the Tuul River Basin landscape was classified according to eight land classes (Figure 7). * Trends at 0.1 significance level. ** Trends at 0.05 significance level. *** Trends at 0.01 significance level.

Analysis of Land Cover
Determining land cover changes using satellite data is important because it allows us to determine land cover changes, shifts, and causes. Study areas were classified using two commonly used supervised classification methods. The maximum-likelihood classification of the Tuul River Basin was used to analyze land cover maps for three different years (2000, 2010, and 2020). To confirm this, there results were compared with the soil, topographic, ecological landscape potential, vegetation, and forest maps. With appreciable accuracy, the Tuul River Basin landscape was classified according to eight land classes (Figure 7). The number of classified pixels varied greatly in several types of land cover. It is noteworthy that most of the maximum changes occurred in the grassland and wetland area from 2000 to 2010. During those 10 years, the grassland area and wetland increased by 59.00 km 2 and 11.89 km 2 . In contrast, the water bodies and cultivated lands decreased by 45.89 km 2 and 28.00 km 2 . However, the change in land cover was even greater between 2010 and 2020, which was almost the opposite of the previous decade. The amount of grassland significantly decreased by −378.67 km 2 during the second decade, while the area of water bodies increased by 74.25 km 2 (Table 4).  The number of classified pixels varied greatly in several types of land cover. It is noteworthy that most of the maximum changes occurred in the grassland and wetland area from 2000 to 2010. During those 10 years, the grassland area and wetland increased by 59.00 km 2 and 11.89 km 2 . In contrast, the water bodies and cultivated lands decreased by 45.89 km 2 and 28.00 km 2 . However, the change in land cover was even greater between 2010 and 2020, which was almost the opposite of the previous decade. The amount of grassland significantly decreased by −378.67 km 2 during the second decade, while the area of water bodies increased by 74.25 km 2 (Table 4). There was also a sharp increase in the size of the bareland, which was not evident in the previous decade. In this case, the high variability of the land cover indicates that there are some factors that strongly influence it.
During 2000 to 2010, the forms of land cover that are most sensitive to water changed significantly. Therefore, the study of land cover transitions and changes may be important in detecting the effects of change (Table 5). From 2000 to 2010, about 69.95 km 2 area of cultivated land was converted into grassland and 0.22 km 2 area of cultivated land was converted into shrubland. The conversion of forest land to grassland, wetland, and water bodies was 1.67 km 2 , 3.70 km 2 , and 0.48 km 2 , respectively. The conversion of grassland to cultivated land, forest, wetland, water bodies, artificial surfaces, and shrubland was 41.02 km 2 , 5.21 km 2 , 12.00 km 2 , 4.14 km 2 , 2.78 km 2 , and 1.73 km 2 , respectively. The conversion of wetland to grassland, water bodies, and shrubland was 8.21 km 2 , 0.43 km 2 , and 0.32 km 2 , respectively. The conversion of water bodies to forest land, grassland, wetland, and shrubland was 2.04 km 2 , 44.89 km 2 , 3.75 km 2 , and 1.03 km 2 respectively. The conversion of artificial surfaces to grassland, wetland, and shrubland was 0.37 km 2 , 0.07 km 2 , and 0.04 km 2 , respectively. The conversion of shrubland to cultivated land, grassland, wetland, water bodies, and artificial surfaces was 0.93 km 2 , 3.56 km 2 , 1.25 km 2 , 0.93 km 2 , and 1.79 km 2 , respectively. The most sensitive types of land cover are grassland and water bodies. The highest increase in land cover type was 59 km 2 for grassland and the highest decrease was 45.89 km 2 for water bodies. From 2000 to 2020, the forms of land cover that are most sensitive to temperature and water changed significantly (Table 6). From 2000 to 2020, about 5.13 km 2 area of artificial surfaces land was converted into grassland. The conversion of cultivated land to artificial surfaces land, grassland, and shrubland was 2.77 km 2 , 292.95 km 2 , and 3.57 km 2 , respectively. The conversion of forest land to bareland, grassland, water bodies, and wetland was 2.19 km 2 , 438.02 km 2 , 4.02 km 2 , and 3.71 km 2 , respectively. The conversion of grassland to artificial surfaces land, bareland, cultivated land, forest land, shrubland, water bodies, and wetland was 96.63 km 2 , 31.93 km 2 , 79.02 km 2 , 557.73 km 2 , 2531.30 km 2 , 62.13 km 2 , and 14.03 km 2 , respectively. The conversion of shrubland to artificial surfaces land, bareland, cultivated land, grassland, water bodies, and wetland was 8.70 km 2 , 70.48 km 2 , 1.19 km 2 , 2177.62 km 2 , 13.91 km 2 , and 1.11 km 2 , respectively. The conversion of water bodies to forest land, grassland, shrubland, and wetland was 4.08 km 2 , 36.55 km 2 , 9.29 km 2 , and 2.26 km 2 respectively. The conversion of wetland to grassland, shrubland, and water bodies was 3.59 km 2 , 4.09 km 2 , and 1.55 km 2 , respectively. Numerous studies to estimate land cover displacement confirm the confidence value by the kappa coefficient [5,[33][34][35]. If the accuracy value of the kappa coefficient classification is 0.61-0.80, it substantially ensures the accuracy of the land cover transition [36][37][38]. In this study, land cover classification for the years 2000 and 2020 was produced and the accuracy checked for the transformation matrix; the kappa coefficient was substantial, 0.78. This means that the estimates of land cover changes in the Tuul River Basin are accurate and interrelationships can be further analyzed.

Analysis of Interrelationships of Ecohydrological Processes
In the Tuul River Basin, air temperature statistically increased sharply at all stations, with Terelj (ϕ = 5.08) and Ulaanbaatar (ϕ = 14.11) the highest (Table 1). This greatly increases evaporation and has a direct effect on the reduction of water and moisture in the surface area. In the semiarid and arid regions of the Central Asian Plateau, precipitation is the main source of surface water, and surface water levels may increase as precipitation increases. It is important to examine whether there is a statistical relationship between hydrological processes and land cover components in the basin. This is especially important in understanding the ecohydrological process by establishing the relationship between river discharge and precipitation, and the components associated with the surface area of water bodies. In the study area, it was observed with a sharp increase in the average annual air temperature and a moderately negative correlation with a decrease in the water body type of land cover (r = −n.54). It was observed that the average annual precipitation has a strong positive, statistically significant correlation with river water discharge (r = 0.70). The interrelationships between the other sections were weak (Figure 8).
The ecohydrological processes in Mongolia are directly related to the amount of precipitation [39]. As a result of direct and indirect human and natural factors, the hydrological processes of Mongolia's major rivers are changing and river discharge is declining [40]. In particular, global warming increases evaporation and affects the land surface, decreasing the water surface [41]. Therefore, this may have had a significant impact on reducing the surface area of water bodies (45.89 km 2 ) in the Tuul River Basin from 2000 to 2010. However, between 2010 and 2020, a slight increase in precipitation in the basin increased the area with water bodies (74.25 km 2 ). Between 2000 and 2020, the amount of grassland decreased by 319.67 km 2 , while the area of water bodies increased by 28.36 km 2 . In addition, due to temperature and evaporation, the amount of grassland area decreased and the amount of bareland increased. Therefore, it is important to establish an interrelationship by comparing the main hydrometeorological parameters, which have a major impact on the ecohydrological processes, with water bodies and sensitive areas of the land cover types (Figure 9).

Analysis of Interrelationships of Ecohydrological Processes
In the Tuul River Basin, air temperature statistically increased sharply at all stations, with Terelj (φ = 5.08) and Ulaanbaatar (φ = 14.11) the highest (Table 1). This greatly increases evaporation and has a direct effect on the reduction of water and moisture in the surface area. In the semiarid and arid regions of the Central Asian Plateau, precipitation is the main source of surface water, and surface water levels may increase as precipitation increases. It is important to examine whether there is a statistical relationship between hydrological processes and land cover components in the basin. This is especially important in understanding the ecohydrological process by establishing the relationship between river discharge and precipitation, and the components associated with the surface area of water bodies. In the study area, it was observed with a sharp increase in the average annual air temperature and a moderately negative correlation with a decrease in the water body type of land cover (r = −n.54). It was observed that the average annual precipitation has a strong positive, statistically significant correlation with river water discharge (r = 0.70). The interrelationships between the other sections were weak (Figure 8). Figure 8. Regression between hydrometeorology and water variable interrelationships, such as (a) air temperature and water bodies type, (b) temperature and river discharge, (c) precipitation and river discharge, and (d) precipitation and water bodies type.
The ecohydrological processes in Mongolia are directly related to the amount of precipitation [39]. As a result of direct and indirect human and natural factors, the hydrological processes of Mongolia's major rivers are changing and river discharge is declining [40]. In particular, global warming increases evaporation and affects the land surface, decreasing the water surface [41]. Therefore, this may have had a significant impact on reducing the surface area of water bodies (45.89 km 2 ) in the Tuul River Basin from 2000 to 2010. However, between 2010 and 2020, a slight increase in precipitation in the basin increased the area with water bodies (74.25 km 2 ). Between 2000 and 2020, the amount of grassland decreased by 319.67 km 2 , while the area of water bodies increased by 28.36 km 2 . In addition, due to temperature and evaporation, the amount of grassland area decreased and the amount of bareland increased. Therefore, it is important to establish an interrelationship by comparing the main hydrometeorological parameters, which have a major impact on the ecohydrological processes, with water bodies and sensitive areas of the land cover types (Figure 9). Studies in the arid and semiarid regions of Central Asia show that climate change has a significant impact on land cover changes [5,42]. When linking changes and shifts in this land cover to a hydrological system, it may be important to consider the spatial and temporal interrelationships of the key factors that affect it. In addition, changes in land cover have occurred in basin areas as a result of human and natural impacts ( Figure 10). The areas that changed the most were within the river zone. In particular, the size of grasslands along the river basin changed significantly and shifted to other forms. There is also a not less of a change in cultivated land and artificial surfaces that are changed by human activities. This is a manifestation of changes in the types of land cover due to human activities [43]. On the other hand, precipitation and water bodies of the basin are clearly declining due to global warming [44]. This is a natural factor affecting the land cover [45].
This baseline study is to determine the impact of climate change on ecohydrological processes, land cover changes, ecosystems, and surface aquifers in a semiarid region. It is consistent with the results of other studies on the hydrological processes and land cover of the Central Asian Plateau with a semiarid climate.

Conclusions
In the Tuul River Basin, hydrometeorological variables, satellite data analysis and determination of land cover, and statistical analysis were used to calculate the values of the interrelationships of ecohydrological processes.
During the study period, the average annual air temperature rose from −u.5 °C to +0.3 °C (1.8 °C). This increase could have a significant effect on the ecohydrological processes and ecosystems of the arid and semiarid regions of Central Asia. The average annual air temperature showed a significantly increasing trend at all stations. In particular, Terelj (φ = 5.08) and Ulaanbaatar (φ = 14.11) had the highest increases. During the study period, the average annual precipitation was 238 mm. The average annual precipitation change during this period was relatively small. A significant decreasing trend was observed at the Terelj (φ = −2.44) and Ulaanbaatar (φ = −0.49) stations, whereas the other stations, Altanbulag, Ondorshireet, Bayannuur, and Orkhontuul stations, demonstrated little increasing trend. The areas that changed the most were within the river zone. In particular, the size of grasslands along the river basin changed significantly and shifted to other forms. There is also a not less of a change in cultivated land and artificial surfaces that are changed by human activities. This is a manifestation of changes in the types of land cover due to human activities [43]. On the other hand, precipitation and water bodies of the basin are clearly declining due to global warming [44]. This is a natural factor affecting the land cover [45].
This baseline study is to determine the impact of climate change on ecohydrological processes, land cover changes, ecosystems, and surface aquifers in a semiarid region. It is consistent with the results of other studies on the hydrological processes and land cover of the Central Asian Plateau with a semiarid climate.

Conclusions
In the Tuul River Basin, hydrometeorological variables, satellite data analysis and determination of land cover, and statistical analysis were used to calculate the values of the interrelationships of ecohydrological processes.
During the study period, the average annual air temperature rose from −1.5 • C to +0.3 • C (1.8 • C 361 • C). This increase could have a significant effect on the ecohydrological processes and ecosystems of the arid and semiarid regions of Central Asia. The average annual air temperature showed a significantly increasing trend at all stations. In particular, Terelj (ϕ = 5.08) and Ulaanbaatar (ϕ = 14.11) had the highest increases. During the study period, the average annual precipitation was 238 mm. The average annual precipitation change during this period was relatively small. A significant decreasing trend was observed at the Terelj (ϕ = −2.44) and Ulaanbaatar (ϕ = −0.49) stations, whereas the other stations, Altanbulag, Ondorshireet, Bayannuur, and Orkhontuul stations, demonstrated little increasing trend.
Changes in land cover varied in relation to hydrometeorological changes. It is noteworthy that most of the maximum changes occurred in the grassland and wetland area from 2000 to 2010. During those 10 years, the grassland area and wetland increased by 59.00 km 2 and 11.89 km 2 . In contrast, water bodies and cultivated lands decreased by 45.89 km 2 and 28.00 km 2 . However, the change in land cover was even greater between 2010 and 2020, which was almost the opposite of the previous decade. In other words, the amount of grassland decreased significantly by −378.67 km 2 during the study period, while the area of water bodies increased by 74.25 km 2 . Between 2000 and 2020, the amount of grassland decreased by 319.67 km 2 , while the area of water bodies increased by 28.36 km 2 . During the study period, changes in land cover were directly related to hydrometeorological main parameters. Water bodies and sensitive areas of the land cover types changed, mainly due to changes in precipitation in the study area. Studies in the arid and semiarid regions of Central Asia show that climate change has a significant impact on land cover changes.
This baseline study is to determine the impact of climate change on ecohydrological processes, land cover changes, ecosystems, and surface aquifers in a semiarid region.