The Assessment of Climate Change on Rainfall-Runoff Erosivity in the Chirchik–Akhangaran Basin, Uzbekistan

Changes in the frequency or intensity of rainfall due to climate always affect the conservation of soil resources, which leads to land degradation. The importance of assessing past and future climate differences plays an important role in future planning in relation to climate change. The spatiotemporal variability of erosivity depending on precipitation using the rainfall erosivity (R) of Universal Soil Loss Equation under the global circulation model (GCM) scenarios in the Chirchik–Akhangaran Basin (CHAB), which is in the northeastern part of the Republic of Uzbekistan, was statistically downscaled by using the delta method in Representative Concentration Pathways (RCPs) 4.5 and 8.5 during the periods of the 2030s, 2050s and 2070s. The (R) was used to determine the erosivity of precipitation, and the Revised Universal Soil Loss Equation (RUSLE) itself determined the effects of changes in erosivity. Ten weather station observational data points for the period from 1990 to 2016 were used to validate the global circulation models (GCMs) and erosion model. The assessment results showed an increase in precipitation from the baseline by an average of 11.8%, 14.1% and 16.3% for all models by 2030, 2050 and 2070, respectively, while at the same time, soil loss increased in parallel with precipitation by 17.1%, 20.5 % and 23.3%, respectively, in certain scenarios. The highest rainfall was observed for the models ACCESS1–3 and CanESM2 on both RCPs and periods, while more intense rainfall was the main reason for the increase in the spatial and temporal erosion activity of the rainfall-runoff. This study is a useful reference for improving soil conservation, preventing water erosion and ensuring the future sustainability of agricultural products, as well as improving the operational management and planning of agriculture.


Introduction
With the advent of agriculture, the value of soil has increased dramatically in people's lives, since land resources are one of the main bases and non-renewable resources for human development [1]. gradually decreases. As we mentioned above, the study area covers the Chatkal (Repu blic of Kyrgyzstan) and Ugam (Republic of Kazakhstan) Rivers, since the formation of water resources begins at the Chatkal and Pskem ranges (Figure 1b). Based on this, the total area of the study area is 22120.9 thousand km 2 . The soils in the CHAB physical-geographical region are diverse, depend on the soil-forming rocks, and are displayed based on geomorphologic zoning [29]. Soil-geomorphological regions are made up of each belt of soil types and subtypes. In the territory of the basin, light brown meadowsteppe soils, mountain brown soils, and gray earth types of soils are identified [30]. Intense and heavy rainfall in the investigated area creates appropriate environmental conditions for the luxurious growth of vegetation [31]. Chirchik-Akhangaran valley vegetation including various herbs: feather grassy, fescue, mallow. It also includes shrubs and trees: almonds, hawthorn, barberry, juniper, maple, wild cherry, nuts, wild apple trees, cherry plums, poplar, willow, birch, spruce, etc. Where light and typical gray soils are developed, ephemerals and ephemerids predominately grow, e.g., bulbous, couch-grass, mountain plum bitter almonds. The main branches of agriculture are cotton growing, horticulture, grape gardening, silkworm breeding and non-rain farming (mainly wheat and barley)-in the foothills of the mountains [32].
The climate in the Chirchik-Akhangaran valley is heterogeneous and is typically semiarid and arid with a strong increase in both precipitation and temperature from the mountains to the plains from the north to south. Considerable impacts are induced by cyclones, humid western air masses, and cold northern air masses, thus causing a lower air temperature and a decline in precipitation. The amount of precipitation in the CHAB (Figure 1) is distributed unevenly, with more rain in the northeast and in the region closer to the mountains [33]; even more precipitation falls in the The soils in the CHAB physical-geographical region are diverse, depend on the soil-forming rocks, and are displayed based on geomorphologic zoning [29]. Soil-geomorphological regions are made up of each belt of soil types and subtypes. In the territory of the basin, light brown meadow-steppe soils, mountain brown soils, and gray earth types of soils are identified [30]. Intense and heavy rainfall in the investigated area creates appropriate environmental conditions for the luxurious growth of vegetation [31]. Chirchik-Akhangaran valley vegetation including various herbs: feather grassy, fescue, mallow. It also includes shrubs and trees: almonds, hawthorn, barberry, juniper, maple, wild cherry, nuts, wild apple trees, cherry plums, poplar, willow, birch, spruce, etc. Where light and typical gray soils are developed, ephemerals and ephemerids predominately grow, e.g., bulbous, couch-grass, mountain plum bitter almonds. The main branches of agriculture are cotton growing, horticulture, grape gardening, silkworm breeding and non-rain farming (mainly wheat and barley)-in the foothills of the mountains [32].
The climate in the Chirchik-Akhangaran valley is heterogeneous and is typically semiarid and arid with a strong increase in both precipitation and temperature from the mountains to the plains from the north to south. Considerable impacts are induced by cyclones, humid western air masses, and cold northern air masses, thus causing a lower air temperature and a decline in precipitation. The amount of precipitation in the CHAB (Figure 1) is distributed unevenly, with more rain in the northeast and in the region closer to the mountains [33]; even more precipitation falls in the mountains. The minimum amount of precipitation falls in the southwest of the district and amounts to 250-300 mM per year; in Sustainability 2020, 12, 3369 4 of 21 the foothill northeastern part, the precipitation reaches 550 mM. In the western part of the mountains, where the air masses experience a forced rise and as a result lead to cooling and reaching a saturation state, as a result of the orographic increase in precipitation [34], the precipitation intensity reaches up to 3000 mM. The average January temperature is from −1.3 • C, to −1.8 • C, the lowest temperature is from −34 • C (the plane,) to −38 • C (mountain peaks), and the average July temperature on the plain is +26.8 • C. In the mountainous part, it is cooler, more humid and the temperature is +20 • C; in contrast, the highest air temperature in summer in the plains reaches +47 • C, and it is much cooler in the mountains [35]. During intense rains, the streams expand their channels and merge in the valley. Figure 2 shows the average values of multiple measured-based climate data of mean monthly precipitation and maximum and minimum temperature in various elevation zones in the basin under consideration during the 1990-2016 time frame. mountains. The minimum amount of precipitation falls in the southwest of the district and amounts to 250-300 mm per year; in the foothill northeastern part, the precipitation reaches 550 mm. In the western part of the mountains, where the air masses experience a forced rise and as a result lead to cooling and reaching a saturation state, as a result of the orographic increase in precipitation [34], the precipitation intensity reaches up to 3000 mm. The average January temperature is from −1.3 °C, to −1.8 °C, the lowest temperature is from −34 °C (the plane,) to −38 °C (mountain peaks), and the average July temperature on the plain is +26.8 °C. In the mountainous part, it is cooler, more humid and the temperature is +20 °С; in contrast, the highest air temperature in summer in the plains reaches +47 °C, and it is much cooler in the mountains [35]. During intense rains, the streams expand their channels and merge in the valley. Figure 2 shows the average values of multiple measured-based climate data of mean monthly precipitation and maximum and minimum temperature in various elevation zones in the basin un der consideration during the 1990-2016 time frame. Monthly average observed precipitation and temperature changes in different zones in the study area based on the data were collected from the Center for Hydrometeorological Service of the Republic of Uzbekistan (UZHYMET), where the basin is divided into two parts: upstream and downstream, which is explained further using the second picture ( Figure 2).
The analysis expressed in Figure 2, clarifies that in the upper stream river when the air masses cool down and reach saturation, in February-April, the orographic increase in precipitation reaches up to 3000 mm, (as occurred during 1990-2016). At the same time, in the downstream river, minimum precipitation reaches 60 mm, (as occurred during 1990-2016).
For all the climatic stations located above and below the river, the average monthly maximum temperatures during the period 1990-2016 were from +27 °C to +36 °C, however, in the summer season in the lower reaches of the river, the maximum temperature can reach +47 °C. The minimum temperature are mainly observed in the winter season and can reach up to −8 °C .

Data Set
All of the stations produced daily observed data during 1990-2016 that were collected from the Center for Hydrometeorological Service of the Republic of Uzbekistan (UZHYMET), with consideration of the basin. We divided the basin into two parts: upstream and downstream (  Monthly average observed precipitation and temperature changes in different zones in the study area based on the data were collected from the Center for Hydrometeorological Service of the Republic of Uzbekistan (UZHYMET), where the basin is divided into two parts: upstream and downstream, which is explained further using the second picture ( Figure 2).
The analysis expressed in Figure 2, clarifies that in the upper stream river when the air masses cool down and reach saturation, in February-April, the orographic increase in precipitation reaches up to 3000 mM, (as occurred during 1990-2016). At the same time, in the downstream river, minimum precipitation reaches 60 mM, (as occurred during 1990-2016).
For all the climatic stations located above and below the river, the average monthly maximum temperatures during the period 1990-2016 were from +27 • C to +36 • C, however, in the summer season in the lower reaches of the river, the maximum temperature can reach +47 • C. The minimum temperature are mainly observed in the winter season and can reach up to −8 • C.

Data Set
All of the stations produced daily observed data during 1990-2016 that were collected from the Center for Hydrometeorological Service of the Republic of Uzbekistan (UZHYMET), with consideration of the basin. We divided the basin into two parts: upstream and downstream (Figure 2), to clarify where more precipitation falls. Pr Upstream and Pr Downstream indicate how the precipitation upstream and downstream of the river basin changed during the period of 1990-2016, respectively. Tmax Upstream, Tmax Downstream, Tmin Upstream, and Tmin Downstream correspond to the variability of the maximum and minimum temperatures upstream and downstream, respectively, of the investigation area during the 1990-2016period.
As sources of future meteorological data, we obtained daily precipitation outputs from 5 GCMs in the CMIP 5 archive. The study had three series of experiments for the future periods of the 2030s (2020-2039), 2050s (2040-2069) and 2070s (2060-2099). A GCM dataset was classified by using the CMIP5 model, and as the scenarios were employed, -RCP 4.5 and RCP 8.5 were statistically downscaled using the delta method (See Table 1). For assessment of the GCMs performance in the Chirchik River basin, we used historical data from 1975-2005 [36] and then interpolated, and the resulting data sets were used as input data. While RCP 4.5 involves the use of a number of technologies and strategies to reduce greenhouse gas emissions and stabilize radiation exposure until 2100 [37], the RCP 8.5 scenario characterizes an increase in greenhouse gas emissions over time due to high population growth and increasing energy demand [38]. According to Luo [39], in the GCM, various runs have a similar capacity to conduct climate projections focusing on Central Asia. All data sets were statistically reduced and resampled into 0.5 × 0.5 grids to unify their spatial resolutions [9,37]. Table 1. General information on selected GCMs, obtained from the Earth System Grid Federation (ESGF, https:// esgf-node.llnl.gov/ projects/ esgf-llnl/ ) portal.

Model
Institution Country Abbreviation

Determination of Soil Erosion Rate
Annual soil loss due to precipitation is a common process in soils [40]. However, when compiling and updating soil maps, in order to better understand the consequences of changes in precipitation intensity and to quickly track, these changes, we selected and used a revised model of the Universal Soil Loss Equation (RUSLE) based on several criteria, as described in [41]: where: A-is the average annual potential soil erosion (tha −1 y −1 ) of the grid located at (i, j); R-is the factor of average rainfall erosivity, (MJ mM ha −1 h −1 yr −1 ); K-is the soil erosion compliance factor (Mg h −1 MJ mM −1 ); and L and S are topographic factors that take into account the length and steepness of the slope, respectively [42]. The LS, C and P values are dimensionless.

Rainfall-Runoff Erosivity Factor (R)
When calculating the precipitation coefficient, the unit of the index is expressed as an indicator of the erosion force of a certain amount of precipitation. In this process, determining as a function of the volume of erosion over any period, a storm or a series of storms is calculated for the intensity and duration of precipitation [43]. However, due to the lack of a measuring station or the lack of data in remote regions, the regression Equations (2) and (3) between the coefficient R and the annual precipitation requires additional measurements of the parameters. Therefore, simplified methods were chosen that ensured the exceptional simplicity of studying the spatial and temporal variability of Sustainability 2020, 12, 3369 6 of 21 precipitation erosion. The rainfall activity data used to obtain the coefficient R are gridded GCM data precipitation from CMIP 5 archive [44].
where: R is the erosion index of rain; and P is the annual rainfall in mM. In two cases, it is recommended to use the following equation: if the annual precipitation < 850 mM, Equation (2) is used, otherwise, if the annual precipitation ≥ 850 mM, then Equation (3)

Soil Erodibility Factor (K)
The soil erodibility coefficient (K), which reflects the rate of soil loss depending on the erosion flow (R factor), and calculated on the basis of soil textures [11], is an empirical measure of soil erosion [45,46] and represents the susceptibility of the soil to erosion [47,48]. The structure and permeability of the soil profile and organic matter are the main soil properties affecting K [46], and the value of K is characterized by the soil texture and permeability of organic compounds depending on the soil type [48] and is modeled with the aid of an equation. We deduced the value of K from the equation of EPIC as follows Equation (4) [48].
The determination of the K-factor for soil, is given in Equation (4), and is based on its texture [11], representing the shares (%) of CLA (clay), SAN (sand) and SIL (silt), respectively, and this result explains the reason for the decrease in the K content in soils with a high content of organic carbon. Here, the value of K was converted from conventional US units to SI units by multiplying it by 0.1317 [48].

Slope Length and Steepness Factor LS
Calculating the slope length according to the digital elevation model (DEM) is a problematic process when modeling erosion. The need to calculate the length and slope of the site provides the most important value for the erosion estimator and depends on the parameters of the factors L and S in the equation, which are usually combined as LS [49], Equation (5) since the physical properties of the Earth have the greatest influence on the occurrence of erosion. To estimate the factor LS, Equations (6), and (7) [50] were applied using the Shuttle Radar Topography Mission (SRTM) with a resolution of 30 m provided by NASA.
2.3.5. Cover Management Factor (C Factor) Among all the risk factors for soil erosion, the C factor is an important indicator that helps preserve soil composition and reduce erosion [41,51]. The factor C, which is considered in the RUSLE model, was used as a factor that determines the vegetation, efficiency and productivity of soil cultures [52]. The calculation of factor C in this study and normalized difference vegetation index (NDVI) data provided by NASA [53] were used according to Equation (8).
where α and β are unitless parameters that determine the curve's shape compared to the NDVI and the C-factor.

Conservation Practice (P Factor)
Factor P is an indicator of the surface conditions that affect the flow and the hydraulic system [54], and reflects the effectiveness of the practice of soil conservation, which promotes a decrease in the volume and speed of streaming water and helps reduce the rate of soil erosion [55]. In the study area, the P-factor was defined based on the studies of Wischmeier and Smith [41], and the values of P given in Table 2 were used on the basis of an analysis of the results in some studies in mountainous regions.

Model Validation
The validity of the derivation of R based on the data used was obtained from the mean monthly rainfall and average annual precipitation (1990-2016) from the daily data collected from the Center for Hydrometeorological Service of the Republic of Uzbekistan.
The validity of the model was evaluated statistically, through comparison using the baseline climate data  and observation data where the Taylor diagrams were used for comparison with the observational data for all 10 stations, when the correlation coefficient (CC) was approximately 0.98% ( Figure 3). The treatment of the rainfall-runoff erosivity model was identified by comparing the R factors from the four factors of measures [56].

Future Precipitation Change
In different periods, the rainfall frequency change for GCMs in the future was determined by baseline data, and all stations showed monthly average changes that were almost identical for RCP 4.5 and RCP 8.5 ( Figure 4). In all the models, the precipitation mainly fell in winter and autumn; in general, this trend was stronger in RCP 4.5 than in RCP 8.5, except for at Tuyabuguz station ( Figure  4i). Despite the decrease in rainfall activity in the summer period, July showed rain activity at all stations for both scenarios. In March, according to the RCP 4.5 scenario, the amount of precipitation at the Bekabad (  (Figure 4j ), stations, recessions began in April. Despite the increase in rainfall intensity in winter at almost all stations in February, a significant decline in the average monthly rainfall was observed in the RCP 8.5 scenario, with a greater proportion of decline at the Bekabad station. When the validity of the derivation of the R factor was evaluated with mean monthly data, no significant differences were observed between the GCMs and theUzhydromet data.

Future Precipitation Change
In different periods, the rainfall frequency change for GCMs in the future was determined by baseline data, and all stations showed monthly average changes that were almost identical for RCP 4.5 and RCP 8.5 (Figure 4). In all the models, the precipitation mainly fell in winter and autumn; in general, this trend was stronger in RCP 4.5 than in RCP 8.5, except for at Tuyabuguz station (Figure 4i). Despite the decrease in rainfall activity in the summer period, July showed rain activity at all stations for both scenarios. In March, according to the RCP 4.5 scenario, the amount of precipitation at the Bekabad  (Figure 4j), stations, recessions began in April. Despite the increase in rainfall intensity in winter at almost all stations in February, a significant decline in the average monthly rainfall was observed in the RCP 8.5 scenario, with a greater proportion of decline at the Bekabad station. When the validity of the derivation of the R factor was evaluated with mean monthly data, no significant differences were observed between the GCMs and theUzhydromet data.  Yangiyul. Figure 5 shows scatter-plots the provide a quantitative comparison of mean monthly precipitation. The scatterplots indicate the corresponding coefficient of determination (R 2 ) and root mean square error (RSME,); for the ACCESS1-3, bcc-csm1-1-m, CanESM2, CSIRO-Mk3-6-0, and GISS-E2-P models, the RSME values were 12.05 mm/month, 6.71 mm/month, 5.33 mm/month, 9.75 mm/month and 3.2 mm/month ( Figure 5), respectively, and the R 2 values were 0.64, 0.92, 0.92, 0.76 ( Figure 5) and 0.88, respectively.  The scatterplots indicate the corresponding coefficient of determination (R 2 ) and root mean square error (RSME,); for the ACCESS1-3, bcc-csm1-1-m, CanESM2, CSIRO-Mk3-6-0, and GISS-E2-P models, the RSME values were 12.05 mM/month, 6.71 mM/month, 5.33 mM/month, 9.75 mM/month and 3.2 mM/month ( Figure 5), respectively, and the R 2 values were 0.64, 0.92, 0.92, 0.76 ( Figure 5) and 0.88, respectively.

Rainfall Erosivity Change under Climate Change
The simulation showed that depending on the rainfall, the spatial distribution of the R values changed. Table 3 shows the effects of precipitation on past and future erosion and erosion intensities

Rainfall Erosivity Change under Climate Change
The simulation showed that depending on the rainfall, the spatial distribution of the R values changed. Table 3 shows the effects of precipitation on past and future erosion and erosion intensities on the CHAB; in all GCM ensembles, precipitation and rainfall-runoff erosion significantly increased. All the scenarios showed that the average value increased from the baseline level in all periods. The average values were 471.6 mM, 489.2 mM and 504.8 mM in the 2030s, 2050s and 2070s, respectively. As a rule, in the 2070s, there was high growth in all RCP scenarios, and the highest growth was forecasted with the ACCESS1-3 and CanESM2 models. Additionally, the lowest value of the increase in the coefficient of radioactivity was observed in the bcc-csm1-1-m model, while the average score corresponded with the GISS-E2-R model. The design results showed that the average rainfall erosion increased to 1081.5 Mj·mm·ha −1 h −1 yr −1 in the 2030s, 1165.2 Mj·mm·ha −1 h −1 yr −1 in the 2050s and 1231.6 Mj·mm·ha −1 h −1 yr −1 in the 2070s, with the above periods and scenarios. Different GCMs had different precipitation activities and rainfall erosivity. All of the models showed no continuous increase or decrease in precipitation activities or rainfall erosivity; however, the average erosion value in all models showed a steady upward trend [57]. Some models indicated a decrease in precipitation and, possibly, a parallel decrease in erosion (Table 3), while some models excluded the possibility of erosion reduction with reduced precipitation. This study showed that an increase in erosivity density and erosion was associated not only with an increase in the amount of precipitation but also with an increase in the frequency of precipitation, since the frequency of precipitation is very important for the erosion process ( Figure 6, Figures 8 and 9 and Table 3). on the CHAB; in all GCM ensembles, precipitation and rainfall-runoff erosion significantly increased.
All the scenarios showed that the average value increased from the baseline level in all periods. The average values were 471.6 mm, 489.2 mm and 504.8 mm in the 2030s, 2050s and 2070s, respectively. As a rule, in the 2070s, there was high growth in all RCP scenarios, and the highest growth was forecasted with the ACCESS1-3 and CanESM2 models. Additionally, the lowest value of the increase in the coefficient of radioactivity was observed in the bcc-csm1-1-m model, while the average score corresponded with the GISS-E2-R model. The design results showed that the average rainfall erosion increased to 1081.5 Mj·mm·ha −1 h −1 yr −1 in the 2030s, 1165.2 Mj·mm·ha −1 h −1 yr −1 in the 2050s and 1231.6 Mj·mm·ha −1 h −1 yr −1 in the 2070s, with the above periods and scenarios. Different GCMs had different precipitation activities and rainfall erosivity. All of the models showed no continuous increase or decrease in precipitation activities or rainfall erosivity; however, the average erosion value in all models showed a steady upward trend [57]. Some models indicated a decrease in precipitation and, possibly, a parallel decrease in erosion (Table 3), while some models excluded the possibility of erosion reduction with reduced precipitation. This study showed that an increase in erosivity density and erosion was associated not only with an increase in the amount of precipitation but also with an increase in the frequency of precipitation, since the frequency of precipitation is very important for the erosion process ( Figures 6, 8, 9 and Table 3).

Factors Affecting the Erosion Process
According to observations from 1990-2016 period, the highest value of the average annual rainfall was 3085 mm (Figure 7a), which corresponded to the north-western and middle parts of the basin. As a rule, the lowest amount of precipitation was observed downstream of the basin in the south and southwest regions of the basin. In the study area, the solubility coefficient of the soil (K) (Figure 7b) ranged from 0 to 0.021; as a rule, 0 means water, and water bodies. Soils in the basin had

Factors Affecting the Erosion Process
According to observations from 1990-2016 period, the highest value of the average annual rainfall was 3085 mM (Figure 7a), which corresponded to the north-western and middle parts of the basin. As a rule, the lowest amount of precipitation was observed downstream of the basin in the south and southwest regions of the basin. In the study area, the solubility coefficient of the soil (K) (Figure 7b) ranged from 0 to 0.021; as a rule, 0 means water, and water bodies. Soils in the basin had sand (58.9%), silt (29.9%) and clay (36.8%). Potential soil degradation and risk factor (%) LS (Figure 7c) values of slope categories (< 5, 5-20, 21-40, 41-70, > 70) [58] were tied to the mountain tops that are up to 4299 m high. Land use/land cover in the study area was classified into 10 types, with the largest proportion (45.62%) being in grassland. It should be noted that the study area is located in a mountainous area. Cropland (22.46%), bare-rock gravel (12.76%), build-up area (7.67%), forest (6.31%), bare (3.53%), permanent glacier snow (0.80%), waterbody (0.48%), wetland (0.28%) and shrubland (0.05%) were the remaining types. According to some studies, the density of factor C (Figure 7d) gradually increased from lower to upper areas, and the highest density was observed in dense forest on mountains. The C-factor in the study object ranged from 0 to 75.4. Analyzing the soil-erosion processes, it can be concluded that the problem of soil conservation is becoming more urgent [59].  [58] were tied to the mountain tops that are up to 4299 m high. Land use/land cover in the study area was classified into 10 types, with the largest proportion (45.62%) being in grassland. It should be noted that the study area is located in a mountainous area. Cropland (22.46%), bare-rock gravel (12.76%), build-up area (7.67%), forest (6.31%), bare (3.53%), permanent glacier snow (0.80%), waterbody (0.48%), wetland (0.28%) and shrubland (0.05%) were the remaining types. According to some studies, the density of factor C ( Figure 7d) gradually increased from lower to upper areas, and the highest density was observed in dense forest on mountains. The C-factor in the study object ranged from 0 to 75.4. Analyzing the soilerosion processes, it can be concluded that the problem of soil conservation is becoming more urgent [59].  Table 3 shows the changes in past and future erosion caused by precipitation activity. When the erosion tendency is observed in all periods compared with the baseline climate, CanESM2 shows that the largest increase is greater than 60, while at the same time it reflects a decrease with a slight difference and the largest decrease falls involving bcc-csm1-1-m were by RCP 8.5, -9.1% by the 2070s. The equivalent values from CSIRO-Mk3-6-0 were RCP 8.5, -3.4% by the 2030s; those from bcc-csm1-1-m were RCP 4.5, -5.7%, those from CSIRO-Mk3-6-0 were RCP 4.5, and -1.9% by the 2050s and those from CSIRO-Mk3-6-0 were RCP 4.5, and -2.6% by the 2070s, which showed a decline, despite high erosion (Table 3). However, the erosivity density ( Figure 8 and Table 3) is not entirely suitable for erosion, and the density is not only associated with the rainfall-runoff activities of rains, but also with other factors affecting the density. According to several studies, an erosion density higher than 1 causes greater erosion from precipitation than does precipitation activity [60,61]. Compared to the base climate, the highest density level occurred in the GISS-E2-R model, with a maximum of 33.3%  Table 3 shows the changes in past and future erosion caused by precipitation activity. When the erosion tendency is observed in all periods compared with the baseline climate, CanESM2 shows that the largest increase is greater than 60, while at the same time it reflects a decrease with a slight difference and the largest decrease falls involving bcc-csm1-1-m were by RCP 8.5, -9.1% by the 2070s. The equivalent values from CSIRO-Mk3-6-0 were RCP 8.5, -3.4% by the 2030s; those from bcc-csm1-1-m were RCP 4.5, -5.7%, those from CSIRO-Mk3-6-0 were RCP 4.5, and -1.9% by the 2050s and those from CSIRO-Mk3-6-0 were RCP 4.5, and -2.6% by the 2070s, which showed a decline, despite high erosion (Table 3). However, the erosivity density ( Figure 8 and Table 3) is not entirely suitable for erosion, and the density is not only associated with the rainfall-runoff activities of rains, but also with other factors affecting the density. According to several studies, an erosion density higher than 1 causes greater erosion from precipitation than does precipitation activity [60,61]. Compared to the base climate, the highest density level occurred in the GISS-E2-R model, with a maximum of 33.3% for GISS-E2-R and -8.5 in the 2050s, and the largest decrease with a slight difference in the case of the base climate was with CSIRO-Mk3-6-0 at -8.5, which marked a decrease by 0.7% to 2030. Despite the erosivity increases shown with the CSIRO-Mk3-6-0 -4.5, by the values for the 2030s and the 2050s show erosion density as slightly decreasing. Soil erosion caused by rainfall intensity can be an increase or decrease [62,63] due to global climate change, which is true since climate variability is determined by what is happening in the atmosphere.

Discussion
The change in the intensity and amount of precipitation occurs mainly due to the changing processes for water vapor in the atmosphere, which circulates from the oceans through the earth [37]. According to several studies [9,64], changes in sea surface temperature have a significant impact on the direction of water vapor from the oceans to land [65]. The formation of climate in the territory under discussion is largely determined by atmospheric processes that are characteristic of the Central Asian region as a whole [35,66]. The duration of precipitation in the study area and the features of their distribution are characteristic of the Central Asian region, with a small amount of falls in the lower part of the region annually. As a rule, more torrential rainfall falls on the mountain areas compared to the downstream in the investigated area. Since mountain slopes help to streamline and enhance convective processes, air masses on the windward slopes experience a forced rise, resulting in the cooling of the moist air masses, thereby increasing precipitation [34,35]. Significant impacts are also caused by cyclones, wet western air masses, and cold northern air masses, leading to lower air temperatures and increased rainfall and, ultimately, can cause erosion with high values which explains the spatiotemporal distribution of rainfall-runoff activity. Erosion also varies due to land use/land cover changes, as land use/land cover plays a very important role in soil erosion [52]. Based

Discussion
The change in the intensity and amount of precipitation occurs mainly due to the changing processes for water vapor in the atmosphere, which circulates from the oceans through the earth [37]. According to several studies [9,64], changes in sea surface temperature have a significant impact on the direction of water vapor from the oceans to land [65]. The formation of climate in the territory under discussion is largely determined by atmospheric processes that are characteristic of the Central Asian region as a whole [35,66]. The duration of precipitation in the study area and the features of their distribution are characteristic of the Central Asian region, with a small amount of falls in the lower part of the region annually. As a rule, more torrential rainfall falls on the mountain areas compared to the downstream in the investigated area. Since mountain slopes help to streamline and enhance convective processes, air masses on the windward slopes experience a forced rise, resulting in the cooling of the moist air masses, thereby increasing precipitation [34,35]. Significant impacts are also caused by cyclones, wet western air masses, and cold northern air masses, leading to lower air temperatures and increased rainfall and, ultimately, can cause erosion with high values which explains the spatiotemporal distribution of rainfall-runoff activity. Erosion also varies due to land use/land cover changes, as land use/land cover plays a very important role in soil erosion [52]. Based on the above considerations, the spatial distribution of precipitation in the Chirchik River basin, in turn, can explain why the erosion process is not the same in the entire research area, and this explains the change in the intensity and amount of precipitation, which in turn shows the impact on erosion.
Studies around the world, predict a significant increase or decrease [67] in erosion under future climate conditions and show that the erosion process is affected by changes in precipitation intensity and quantity as a result of climate change. For instance, estimates in accordance with GCMs [68] suggest that more than 70 percent of the territory of the European continent, will experience a 19% increase in precipitation erosion by 2050. According to [69,70], erosion varied significantly in different scenarios without any particular direction or magnitude for Eastern Africa. Given that not all GCM models show a different direction or amount of erosion, they attribute this to the uncertainty inherent in GCM scenarios. Meanwhile, studies in the United States [71] under in all GCM scenarios, the value of rainfall changes increases over time. Thus, with the inability to capture, precipitation variability equally due to the inability of GCMs, this study reflects the spatial variability of erosivity in the Chirchik-Akhangaran basin based on different GCMs (Table 3, Figure 9).  Figure 9). The Chirchik-Akhangaran district is exposed to outbreaks of various air masses [35]. Climate change caused by global warming as a result of urbanization and industrialization in the middle and southern parts of the basin may have a direct impact on airflow, resulting in simultaneous in increased rainfall and increased erosion. In addition, soil erosion depends not only on climate and total rainfall but is also closely related to water resources and soil cover, its composition and biological effects, physical and geographical features, plant dynamics and soil composition. Moreover, this has been confirmed in other prediction studies (e.g., [9,35,70,72]). External influences are also important in the process of erosion. Deforestation, mechanical measures, land reclamation, land development and, ultimately, improper use of water resources or inappropriate distribution of water resources also play an important role in soil erosivity. The Chirchik-Akhangaran district is exposed to outbreaks of various air masses [35]. Climate change caused by global warming as a result of urbanization and industrialization in the middle and southern parts of the basin may have a direct impact on airflow, resulting in simultaneous in increased rainfall and increased erosion. In addition, soil erosion depends not only on climate and total rainfall but is also closely related to water resources and soil cover, its composition and biological effects, physical and geographical features, plant dynamics and soil composition. Moreover, this has been confirmed in other prediction studies (e.g., [9,35,70,72]). External influences are also important in the process of erosion. Deforestation, mechanical measures, land reclamation, land development and, ultimately, improper use of water resources or inappropriate distribution of water resources also play an important role in soil erosivity.
In the study area, slopes also play an important role, which indicates a high rate of rainfall erosion in mountainous regions in the northwest as well as in the eastern and middle parts of the basin (Figure 9). The average annual R coefficient showed differences in different GCMs, while in the ACCESS1-3 and CanESM2 (both RCPs and periods) models the highest values were observed in the northeast and east of the basin when the average values were observed in the north-western and eastern parts of the basin, and the lowest values were observed in the low part, according to the GISS-E2-R (Table 3) model. Hence, the observation data gradually increased in the middle and upstream areas of the basin and decreased in the downstream.
Studies around the world show that rainfall increases over time and space [67]. All this may be due to changes in the frequency and intensity of precipitation, an increase in temperature and changes in land use, which indicates that future rainfall will significantly affect soil erosion. The erosion process in the study area is exclusively associated with precipitation. The results show that the GCMs selected for the study area do not have the same growth; however, wind changes affect the air mass, resulting in a change in the spatial distribution of rainfall intensity meaning changes in the expected erosion in the region can occur.
The impacts of physical and mechanical factors and the lack of a system to prevent erosion can further accelerate rainfall erosion in the future, and this can lead to a strong deposition on the channel of the Chirchik River, which is the Charvak reservoir; the consequences of this trend will critically affect the long-term management of agriculture and water resources in the region. Evaluation of the inverse relationship between soil cover and precipitation intensity allows for predicting future trends in soil erosion for the study area, which allows us to quickly make decisions on the operational management of sustainable agriculture and water resources management in the region.

Conclusions
This study demonstrated climate models combined with factors in the RUSLE model, which is a practical and relevant method for assessing the spatiotemporal variability of soil erosion and for determining long-term measures in rural management. The use of climate models and appropriate scenarios is crucial for research that addresses the uncertainties of the future and makes quick management decisions.
According to an estimate, the most affected areas by the influence of rainfall erosivity shown in the models ACCESS1-3 and CanESM2. As a result, growing trends in annual rainfall erosivity from the base climate to GCMs, show a positive change in average annual rainfall erosivity by 11.8%, 14.1% and 16.3% in the 2030s, 2050s, and 2070s, respectively compared to the baseline . This was found despite a decline in rainfall-runoff erosivity on the GISS-E2-R model of about 31% and 32% from the baseline climate for both RCPs, respectively, while there is an increase in average erosion from the baseline maximum of 18.8%. Overall, there is an increase in the average erosion for most of the model ensembles from the 2030s through the 2070s period. Despite the increase in the average value in all GCMs, when compared with the baseline climate in some models, average annual soil loss, erosivity density and rainfall activity decreased by small values that were incompatible with each other in the 2030s, 2050s and 2070s, respectively, which once again confirms the importance of soil erosion studies with the rainfall intensity. The cumulative, average annual rainfall intensity and erosion in all climate models and under all scenarios showed steady growth compared to the base climate.
Assessing future soil loss using the GCMs within the empirical RUSLE model gives an idea of the influence of rainfall intensity changes in the Chirchik River basin and allows for operational management of long-term agriculture, thereby improving the ecological situation in the region and ensuring the sustainability of the local economy. The result of this study may be typical but may also be useful in predicting the occurrence of a negative effect caused by climate change in a region. It is recommended to use other GCMs and scenarios in the future for assessing erosion processes in the basin, taking into account any factors affecting erosion, such as wind erosion. Hydrotechnical and agrotechnical measures, the restoration of vegetation cover, creation of a green shield, as well as organizational and economic, hydro-technical and forest reclamation should all be directed at preventing rainfall-runoff erosivity.