High greenhouse gas emission scenarios increase soil erosion risk: A case study of the Min-Tuo River Basin, China

High greenhouse gas emission scenarios increase soil erosion risk: A case study of the Min-Tuo River Basin, China | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article High greenhouse gas emission scenarios increase soil erosion risk: A case study of the Min-Tuo River Basin, China Nan Jiang, Fuquan Ni, Yu Deng, Mingyan Wu, Mengyu Zhu, Yuxuan Wang, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4952297/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract To quantify the current status of soil erosion in the Min-Tuo River Basin and predict future soil erosion conditions, the RUSLE model and GIS technology were employed to assess soil erosion in the Min-Tuo River Basin from 2000 to 2020. Based on the CMIP6 model, precipitation in the Min-Tuo River from 2021 to 2050 was predicted, and the CA-Markov model was used to forecast land use types from 2030 to 2050, thereby predicting the spatial distribution trends of soil erosion in 2030, 2040, and 2050. The results indicate that: (1) The overall soil erosion intensity in the study area increased from 2000 to 2020; (2) Future precipitation in the basin, predicted using the CMIP6 model, showed a fluctuating upward trend compared to historical levels. The SSP5-8.5 scenario had the largest fluctuation; (3) In terms of future land use types, an increase in the area of converted cropland may increase the risk of soil erosion; (4) The area affected by soil erosion in the future will increase, with the highest annual average soil erosion modulus and total amount occurring under the SSP5-8.5 high emission scenario; (5) The dynamic evolution of future soil erosion levels indicates that the stability of soil erosion levels under the SSP2-4.5 scenario is higher than under the other scenarios. In the future high CO2 emission scenario, the risk of soil erosion in the Min-Tuo River Basin will increase, necessitating the active implementation of various prevention and control measures to prevent further exacerbation. Soil erosion CMIP6 RUSLE model CA-Markov model Min-Tuo River Basin Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 1. Introduction Soil erosion is a global environmental threat (Keno & Suryabhagavan, 2015). It is the main cause of soil degradation and is related to soil compaction, soil structure and organic matter loss, and soil acidity (Hessel et al., 2011 ). Almost no form of land on our planet is free from soil erosion (Bekele, 2019 ). Approximately 30% of the global land area and 80% of ag-ricultural land experience soil erosion, with around 3 billion people living on degraded land (Mirzabaev et al., 2016 ). Soil erosion is particularly severe in developing regions with high population density and ecological vulnerability (Renard et al., 2011 ). The characteristics of the Min-Tuo River Basin align with these conditions. Additionally, severe forest resource de-struction, the 2008 Wenchuan earthquake (Guo et al., 2012 ), and extreme rainfall events (Piao et al., 2019 ) have further worsened the region's ability to retain soil and resulted in local soil desertification and lithification, leading to serious soil and water loss in the area. Therefore, quantifying soil loss and predicting future soil erosion is of great significance as it can provide scientific basis for watershed soil and water loss management and macro-scale decision-making for soil and water conservation, thus promoting ecological environment protection and restoration in the watershed. The issue of soil erosion has become a subject of study for numerous scholars. Among these, the use of soil erosion models has proven to be the most widely effective method for researching soil erosion. Several scholars have successively developed models such as the USLE (Wischmeier & Smith, 1958 ), WEPP (Woodward, 1999 ), and LISEM (De Roo, 1996 ), with the USLE model being the most widely utilized. Subsequently, scholars have made supplementary revisions to the USLE model. In 1997, the United States Department of Agriculture proposed the RUSLE model (Renard & Ferreira, 1993 ), which comprehensively reflects the influencing factors of soil erosion and has been widely used (Rymszewicz et al., 2015 ; Ma et al., 2023 ; Wang & Su, 2020 ). With the development of GIS and RS technologies, their combined use allows for the comprehensive exploration of spatial information and patterns related to soil erosion, thereby providing a technological possibility for studying the spatiotemporal variations in soil erosion. Over the past century, global climate has been trending significantly warmer (Suresh, 2023 ; Hairsine & Rose, 1992 ), leading to increased regional differences in precipitation. Rainfall is one of the important factors influencing soil erosion, and global climate change can cause changes in regional precipitation, thereby affecting soil erosion conditions.Soil erosion is a crucial aspect of sustainable development. Analyzing soil erosion factors and processes using soil erosion models can provide a theoretical basis for formulating regional soil and water conservation measures (Nearing et al., 1999 ). Therefore, this study is based on the RUSLE model to investigate the current status of soil erosion in the Min-Tuo River Basin and predict future soil erosion conditions based on the CMIP6 model. In the past, many related studies have studied the current situation of soil erosion in the Min-Tuo River Basin. For instance, Zhong et al. ( 2022 ) used the RUSLE model to assess soil erosion in the Tuojiang River from 2000 to 2018 and explored the spatiotemporal dynamic evolution patterns. Deng & Shi ( 2022 ) quantitatively evaluated soil erosion and the spatial distribution of soil and water conservation ecosystem services in the Min-Tuo River Basin using the USLE model. Zhang et al. ( 2020 ) generated a soil erosion modulus map for the Min-Tuo River Basin using RUSLE based on supervised classification of the current land use map, indicating that the basin is mainly characterized by mild and moderate erosion areas. Liu et al. ( 2019 ) investigated the erosion status and characteristics of the upper reaches of the Min-Tuo River Basin and explored erosion risk in the region using a GIS weighted overlay analysis method. However, there is currently very little research on the overall soil and water loss situation in the entire Min-Tuo River Basin after 2010. Given the background mentioned above, this study focuses on the Min-Tuo River Basin as the main area. Through the use of remote sensing images, precipitation records, DEM, soil types, vegetation types, land use, and other data, the study employed the RUSLE model to analyze soil erosion in the Min-Tuo River Basin from 2000 to 2020. Furthermore, based on the CMIP6 model, six models and three different emission scenarios were selected. After bias correction and downscaling using the multi-model ensemble (MME) approach, the study predicted the precipitation in the basin from 2021 to 2050. Additionally, with the assistance of the CA-Markov model, the study forecasted land use types for the years 2030, 2040, and 2050, and subsequently predicted the soil erosion situation in the Min-Tuo River Basin from 2030 to 2050. This research can provide theoretical references for soil erosion control and soil conservation measures in the study area under future climate change scenarios. It can also scientifically predict the soil erosion in large river basins in China at a spatial scale, aiming to provide a scientific and rational research basis for the ecological environment and comprehensive management projects in the Min-Tuo River Basin. 2. Materials and Methods 2.1 Study area The The Min-Tuo River Basin (Fig. 1 ) is mainly composed of four parts: the Minjiang, Dadu River, Qingyijiang, and Tuojiang (Hu et al., 2020 ). Particularly, the upper reaches of the Minjiang and the arid valleys of the Dadu River have extremely fragile ecosystems. Under the dual pressure of population and economic development, the soil erosion in the Min-Tuo River Basin is becoming more and more serious. With the acceleration of urbanization, industrialization and agricultural modernization, the expansion of economic scale in the future will have an important impact on soil erosion in the Min-Tuo River Basin. At the same time, with the implementation of national strategies such as the Yangtze River Economic Belt and the Chengdu-Chongqing Economic Circle, higher requirements have been put forward for soil and water conservation in river basins. 2.2 RULSE model The In recent years, with the maturity of "3S" technology, the difficulty of obtaining geographical information by traditional research methods has been overcome, the efficiency of soil erosion assessment has been improved, and the RUSLE model has been better applied. The RUSLE model expression is as follows: \(\:A=R\times\:K\times\:L\times\:S\times\:C\times\:\) P (1) Where \(\:\text{A}\) is soil and water conservation quantity (t/hm2·a); \(\:\:\text{R}\) is rainfall erosion factor (MJ·mm/hm2·h·a); \(\:\text{K}\) is soil erodibility factor [t·h/(MJ·mm)]; \(\:\:\text{L}\text{S}\) is topographic factor, dimensionless; \(\:\:\text{C}\) is vegetation cover factor, dimensionless; P is soil and water conservation measure factor, dimensionless. The rainfall data partly comes from the National Science and Technology Infrastructure of China, using the 1km resolution annual precipitation data from 2001 to 2020 in China. This data is derived by synthesizing the yearly cumulative sum based on the monthly precipitation data set at a 1km resolution by Professor Peng Shouzhang from 1901 to 2020, as referenced in Peng et al. ( 2017 ) and Su et al. ( 2022 ). Another part is sourced from the Resource and Environment Data Center of the Chinese Academy of Sciences, which is the spatially interpolated data set of annual mean temperature and annual precipitation since 1980 in China. This data set is based on daily observational data from over 2400 meteorological stations in China, processed through collection, statistics, and spatial interpolation calculations. The NDVI data is sourced from the National Ecological Science Data Center of China, as referenced in Yang et al. ( 2019 ). The land use/cover data is sourced from the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences. The soil data is sourced from the World Soil Database, where the HWSD (Harmonized World Soil Database) can be queried for soil texture data, including sand, silt, clay, and organic carbon content. The terrain data uses a 30m resolution Digital Elevation Model (DEM) data sourced from the CGIAR-CSI SRTM Elevation Database. For the calculation of soil and water conservation, it is necessary to ensure consistency in the resolution and projection coordinate system of each factor when using them collectively. Please refer to Table 1 for the sources and usages of the relevant data. Table 1 Data sources. Data classification Processing Resolution Data sources Rainfall data Estimated R-value 1km Data Center for Resources and Environment, CAS/National Earth System Science Data Center Soil data Estimated K value 1km World Soil Database Terrain data Estimated LS value 30m CGIAR-CSI SRTM Elevation Database NDVI data Estimated C value 30m National Ecological Science Data Center Land use/cover data Estimated P value 30m Institute of Geographical Sciences and Natural Resources, CAS 2.2.1. Rainfall erosivity factor( \(\:R\) ) Rainfall erosivity factor is one of the potential forces causing soil erosion, which can be judged by annual rainfall data.Guidelines for Calculation of Soil Loss Amount of Production and Construction Projects (SL773-2018).The formula for calculating \(\:\text{R}\) is as follows: $$\:R=0.067{P}_{d}^{1.627}$$ 2 Where \(\:R\) is the erosive factor of annual average rainfall, and the unit is MJ·mm/(hm²·h); \(\:{P}_{d}\) is annual average rainfall, unit is mm. 2.2.2. Soil erodibility factor( \(\:K\) ) Soil erodibility factor is one of the essential parameters of RUSLE model, which is used to reflect soil sensitivity to erosion. \(\:\text{K}\) can be determined by soil texture data by Rao et al.(2020)、Williams ( 1990 ) and Fan et al.(2012). An algorithm for estimating soil erodibility factor in EPIC model is proposed to calculate \(\:\text{K}\) . The specific formula is as follows: Where \(\:\text{K}\) is soil erodibility factor, the unit is t·h/(MJ·mm); \(\:\:\text{S}\text{A}\text{N}\:\) is sand content (%); \(\:\text{S}\text{I}\text{L}\) is silt content (%); \(\:\text{C}\text{L}\text{A}\) is clay content (%); \(\:\text{C}\) is organic carbon content (%). 2.2.3. Topographic factors ( \(\:LS\) ) Topographic factor is also one of the main causes of soil erosion. Slope length factor ( \(\:\text{L}\) ) and slope gradient factor ( \(\:\text{S}\) ) are selected in RUSLE model. \(\:\text{L}\text{S}\) reflects the influence of topography fluctuation on soil erosion. The greater the \(\:\text{L}\text{S}\) , the faster the soil erosion rate Although it is difficult to measure \(\:\text{L}\text{S}\) in large-scale research, the influence of topographic factors on soil erosion can still be expressed by topographic fluctuation rate, that is, the maximum height change in a certain height area of the ground (Wischmeier & Smith, 1978 ). The specific formula is as follows: $$\:H={H}_{max}-{H}_{min}$$ 5 Where \(\:H\) topographic relief, \(\:{H}_{max}\) is the highest value of elevation, \(\:{H}_{min}\) is the lowest value of elevation. 2.2.4. Vegetation cover factor( \(\:C\) ) Vegetation cover factor plays an important role in soil and water conservation, \(\:\text{C}\:\) ranges from Between 0 and 1(Yang & Huang, 2021 ). When \(\:\text{C}\) =1, the ground is completely bare without vegetation; when \(\:\text{C}\) =0 the ground vegetation coverage is good. It is a common way to obtain vegetation coverage from remote sensing images and calculate vegetation coverage factors according to vegetation coverage. Cai et al.(2000)The proposed algorithm is used to estimate the vegetation cover factor through artificial and natural rainfall experiments. The specific formula is as follows: fvc \(\:=\frac{NDVI-{NDVI}_{min}}{{NDVI}_{max}-{NDVI}_{min}}\) (6) $$\:C=\left\{\begin{array}{c}1,\:fvc\le\:0.095\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\\\:0.6508-0.3436Inc,0.095\:\:<\:fvc<0.783\:\:\\\:0,\:fvc\ge\:0.783\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\:\end{array}\right.\:\:$$ 7 Where \(\:\text{N}\text{D}\text{V}\text{I}\) is normalized vegetation index(Wei et al., 2022 ); \(\:{\text{N}\text{D}\text{V}\text{I}}_{\text{m}\text{i}\text{n}}\) and \(\:{\text{N}\text{D}\text{V}\text{I}}_{\text{m}\text{a}\text{x}}\) are the minimum and maximum values of \(\:\text{N}\text{D}\text{V}\text{I}\) respectively; \(\:\text{f}\text{v}\text{c}\) is vegetation cover index (%); \(\:\:\text{C}\) is vegetation cover factor. 2.2.5. Soil and water conservation measures( \(\:P\) ) Soil and water conservation measures can regulate the amount of soil erosion in a region. The factor of soil and water conservation measures is difficult to calculate directly. Usually, the value of \(\:\text{P}\) factor of a region is determined according to different land use types. The range of \(\:\text{P}\) is between 0 and 1. The larger \(\:\text{P}\) is, the smaller the amount of soil erosion is. When \(\:\text{P}\) =0, it means that no erosion occurs. When \(\:\text{P}\) =1, it means that no soil and water conservation measures are taken (Humphrey et al., 2021). According to the research conclusions of relevant scholars (Peng et al., 2007 ). The following values were assigned to the study areas with different land use properties (Table 2 ). Using ArcGIS software to reclassify the assigned data, you can obtain \(\:\text{P}\) values and their distribution maps for different land use types within the basin. Table 2 Factor values of soil and water conservation measures in the study area. land use types cultivated land forest land grass waters construction land unused land soil and water conservation factor 0.25 1 1 0 0 0 2.3. CMIP6 climate model data According to the applicability of CMIP6 in the study area (Lu et al., 2022 ; Wu et al., 2023 ), six atmospheric circulation models with better performance in CMIP6 (Table 3 ) were selected. Due to the differences among the selected six climate models, the bilinear interpolation method was used to unify the spatial resolution of all models. After the interpolation, Delta downscaling (Chen et al., 2011 ) was performed on the interpolated model data at each grid. Finally, the daily meteorological data for the Min-Tuo River Basin under three shared socioeconomic pathways (Table 4 ) were obtained for the future period. Table 3 Information on the six CMIP6 patterns used in this study. Number Schema name Countries Spatial resolution (°) 1 ACCESS-ESM1-5(ACC) Australia 1.2°× 1.8° 2 BCC-CSM2-MR(BCC) China 2.8°× 2.8° 3 CanESM5(CAN) Canadian 2.8°× 2.8° 4 IPSL-CM6A-LR(IPSL France 1.3°× 2.5° 5 MRI-ESM2-0(MRI) Japan 1.1°× 1.1° 6 MIROC6(MIR) Japan 1.4°× 1.4° Table 4 Introduction to SSPs Number Path name Forced category Shared socio-economic pathways 2100 year radiative forcing/(W/m²) 1 SSP 1-2.6 low forcing sustainable development path 2.6 2 SSP 2-4.5 low forcing intermediate path 4.5 3 SSP 5-8.5 high forcing Traditional fossil fuel-based pathways 8.5 2.2.4. CMIP6 data processing The Delta method (Chen et al, 2011 ) is used to downscale and bias correct the CMIP6 data of the selected 6 models. This method compares the differences between GCMs output and observational data, and downscales the climate projection from large scale to smaller spatial scales. It has the advantages of strong interpretability, high accuracy, and high computational efficiency. When describing the changes in climate elements, precipitation is represented using the rate of change, while temperature is represented using the amount of change. The calculation formulas for future precipitation and temperature scenarios at meteorological stations are as follows: $$\:{P}_{f}={P}_{o}\bullet\:\frac{{P}_{Gf}}{{P}_{Go}}$$ 8 $$\:{T}_{f}={T}_{o}+\left({T}_{Gf}-{T}_{Go}\right)$$ 9 In the equations, \(\:{\text{P}}_{\text{f}}\) 、 \(\:{\text{T}}_{\text{f}}\) represent the reconstructed future precipitation and temperature sequences; \(\:{\text{P}}_{\text{G}\text{f}}\) 、 \(\:{\text{T}}_{\text{G}\text{f}}\) represent the future precipitation and temperature sequences estimated by the climate model; \(\:{\text{P}}_{\text{G}\text{o}}\) 、 \(\:{\text{T}}_{\text{G}\text{o}}\) represent the multi-year average precipitation and temperature simulated by the climate model during the reference period; \(\:{\text{P}}_{\text{o}}\) 、 \(\:{\text{T}}_{\text{o}}\) represent the multi-year average precipitation and temperature of the observed field during the reference period. The standardized Taylor diagram (Taylor, 2001) is used to assess the model performance of different climate models. It effectively and intuitively illustrates the differences in performance among multiple models and the magnitude of the errors between simulated and actual values, and has been widely applied in climate model evaluation research. Single-model to multi-model ensemble prediction is one of the effective ways to improve model accuracy, and the calculation formula for using MME is as follows: $$\:{F}_{MME}=\frac{1}{n}\sum\:_{i=1}^{n}{F}_{i}$$ 10 In the equation, \(\:{\text{F}}_{\text{M}\text{M}\text{E}\:}\) represents the average result of the multi-model ensemble, and \(\:\:{\text{F}}_{\text{i}\:}\:\) represents the simulation results of each individual model. Quantitative evaluation was conducted on the future precipitation results of the Min-Tuo River Basin estimated by the six climate models and the multi-model ensemble average, which intuitively shows the relationship between the predicted precipitation and observed precipitation under the six models and the ensemble average(Fig. 2 ) Therefore, overall, the future precipitation simulation values after the multi-model ensemble average are closer to the observed precipitation values than those of a single model, and the reliability of the simulation results is higher. Therefore, in this study, the multi-model ensemble average was used as the future precipitation data, and the future soil erosion was predicted based on this. 2.2.5. CA-Markov Model In order to estimate the soil erosion modulus of the Min-Tuo River Basin under future climate conditions, it is necessary to predict the land use types in order to obtain the P factor values. Because the Cellular Automata (CA) model and Markov model complement each other well, the CA-Markov model is widely used for future land prediction (Abijith & Saravanan, 2021 ; Aqil et al., 2022 ; Wang & Wang, 2022 ). Therefore, the CA-Markov model is used to predict future land use in the study area. The \(\:\text{K}\text{a}\text{p}\text{p}\text{a}\) coefficient can be used to test the consistency between the simulation results and the actual data, and is commonly used to assess the accuracy of predicting land use changes. The formula for calculating the \(\:\text{K}\text{a}\text{p}\text{p}\text{a}\) coefficient is as follows: $$\:Kappa=\frac{{p}_{a}-{p}_{c}}{{p}_{p}-{p}_{c}}$$ 11 The formula where \(\:{\text{p}}_{\text{a}}\) represents the probability of correct simulation, \(\:{\text{p}}_{\text{c}}\) represents the probability of simulated prediction, and \(\:{\text{p}}_{\text{p}}\) represents the probability of ideal simulation. The \(\:\text{K}\text{a}\text{p}\text{p}\text{a}\) coefficient ranges between [-1,1], and when 0.6≤ \(\:\text{K}\text{a}\text{p}\text{p}\text{a}\) ≤0.8, it indicates significant consistency; when 0.8< \(\:\text{K}\text{a}\text{p}\text{p}\text{a}\) ≤1, it indicates very good consistency. 3. Results and analysis 3.1 Soil Erosion Historical Analysis Based on the calculated soil erosion impact factors, the soil erosion modulus for the years 2000, 2010, and 2020 was obtained, and the intensity of soil erosion was analyzed. According to the "Classification and Grading Standards for Soil Erosion in China" (SL190-2007) (2008), the soil erosion intensity in the study area was divided into six erosion levels: slight erosion, mild erosion, moderate erosion, severe erosion, very severe erosion, and extreme erosion (Table 5 ). Table 5 Classification of soil erosion modulus. level micrometric mild moderate strong extremely strong violently Erosion modulus [t/(km²·a)] 15000 As seen from Fig. 3 , the erosion in Ya'an, Liangshan Prefecture, Ganzi Prefecture, and the southern part of Aba Prefecture is showing an increasing trend. According to Table 6 , it is known that the soil erosion modulus reached its highest value in 2020, at 34939.8t/(km²·a). The overall trend over the 20 years from 2000 to 2020 first decreased and then increased. The increase from 2010 to 2020 was greater than the decrease from 2000 to 2010, at 7.1% and 2.6% respectively. Table 6 Statistical table of soil erosion modulus in the study area from 2000 to 2020. time Minimum t/(km²·a) Maximum value [t/(km²·a)] Range of change (%) 2000 0 33493.6 2010 0 32621.4 -2.6 2020 0 34939.8 7.1 From Table 7 , it can be observed that: (1) In 2000, slight erosion and moderate erosion were the main types, covering areas of 112688.41km² and 20007.15km², accounting for 69.72% and 12.38% of the total erosion area, respectively; (2) In 2010, slight erosion and moderate erosion were predominant, covering areas of 97840.53km² and 25076.86km², accounting for 60.54% and 15.52% of the total erosion area, respectively; (3) In 2020, slight erosion and moderate erosion were predominant, covering areas of 54167.59km² and 48477.77km², accounting for 33.52% and 29.99% of the total erosion area, respectively; (4) From 2000 to 2020, the area of slight erosion showed a decreasing trend, reaching its minimum in 2020 at 54167.59km², accounting for 33.52% of the total erosion area. The areas of mild, moderate, and severe erosion showed an increasing trend over the 20 years, reaching their maximum values in 2020 at 28669.73km², 48477.77km², and 24358.86km², accounting for 17.74%, 29.99%, and 15.07% of the total erosion area, respectively. The areas of very severe and extremely severe erosion showed an increasing trend over the 20 years, but with a smaller magnitude, reaching their maximum values in 2020 at 5862.79km² and 84.56km², accounting for 3.63% and 0.05% of the total erosion area, respectively. Table 7 Statistical Table of Soil Erosion Area in Study Area from 2000 to 2020. Year 2000 2010 2020 Area (proportion) Area/km² Ratio/% Area/km² Ratio/% Area/km² Ratio/% type Micro erosion 112688.41 69.72 97840.53 60.54 54167.59 33.52 Mild erosion 17808.95 11.02 17830.72 11.03 28669.73 17.74 Moderate erosion 20007.15 12.38 25076.86 15.52 48477.77 29.99 Intense erosion 8677.3 5.37 16063.14 9.94 24358.86 15.07 Extreme erosion 2403.79 1.49 4743.96 2.94 5862.79 3.63 Severe erosion 35.7 0.02 66.1 0.04 84.56 0.05 The Sankey diagram (Fig. 4 ) provides a more intuitive view of the erosion grade transitions. 3.2. Soil Erosion Future Forecast and Analysis 3.2.1. Precipitation Forecast The precipitation forecast based on the multi-model ensemble average is shown in Fig. 5 . In spatial terms, Wang et al. ( 2021 ) used China as the study area, and the results indicated that the spatial estimates of precipitation and extreme precipitation indices from CMIP6 were higher than the actual values. In this study, based on the MME simulations, the future high precipitation area in the Min-Tuo River basin is located in the southwest of the study area. This may be due to the enhanced disturbance of climate in the high-altitude areas of the Min-Tuo River basin under the selected future climate change scenarios, resulting in abnormally high precipitation in these high-altitude regions. The low precipitation area in the future is located in the northern part of the study area, which aligns with the actual spatial distribution of annual precipitation in the study area at present. 3.2.2. Land Use Forecast With the help of IDRISI software, the Kappa index for simulated land use and actual land use in the study area in 2010 and 2020 was calculated to be 0.87 and 0.86, respectively, indicating good consistency. Considering that the forecast period of this study is from 2030 to 2050, we selected representative years 2030, 2040, and 2050 with a 10 year interval to visually display the evolution of future land use types in the Min-Tuo River basin (Fig. 6 ). Under future climate change scenarios, the area of cropland conversion increases. However, if agricultural activities are not managed properly, this could potentially increase the risk of soil erosion in the Min-Tuo River basin. 3.2.3. Soil Erosion Time Distribution Simulation Using the MME method, the future soil erosion values for the Min-Tuo River basin were calculated for three scenarios. Considering the length of the time series, representative years were selected at 10-year intervals to calculate the average soil erosion modulus and total erosion for the Min-Tuo River basin (Table 8 ). Over time, there are differences in soil erosion in the Min-Tuo River basin under different scenarios. Comparatively, in the SSP5-8.5 scenario, the multi-year average soil erosion modulus and total amount in the Minjiang River basin are ranked as follows: SSP5-8.5 > SSP1-2.6 > SSP2-4.5. Table 8 Soil erosion statistics of Min-Tuo River in the future. SSPs Year Average soil erosion modulus [t/(km²·a)] Total erosion(10⁸t) SSP1-2.6 2030 442.88 7.22 2040 442.86 7.22 2050 445.46 7.26 SSP2-4.5 2030 437.39 7.13 2040 437.99 7.14 2050 440.64 7.18 SSP5-8.5 2030 446.58 7.28 2040 443.36 7.23 2050 447.74 7.3 In the SSP5-8.5 scenario, the future changes in soil erosion in the Min-Tuo River basin exhibit a "V" shape, meaning that the total soil erosion decreases before 2030, but gradually increases after 2030. In terms of the year when the maximum erosion occurs, in all three scenarios, the maximum soil erosion occurs in 2050. This may be due to the fact that the maximum rainfall also occurs in the 2040s in all three scenarios. Over the future time scale, the increase in heavy rainfall or precipitation variability leads to soil erosion and substantial soil loss. Figure 7 presents the average soil erosion intensity level in the Min-Tuo River basin during the forecast period under the three different scenarios. Therefore, it is necessary to establish an early warning system for soil erosion forecasting to address the future changes in soil erosion. 3.2.4. Prediction of spatial evolution of soil erosion Furthermore, the spatial distribution of soil erosion in Min-Tuo River in the future is discussed in three scenarios on the spatial scale (Fig. 8 ). Under SSP1-2.6 scenario, the annual average soil erosion intensity grade and intensity change of Min-Tuo River in the future are shown in Fig. 8 (a). Under SSP1-2.6 scenario, the area of soil erosion in Min-Tuo River reaches 5.35×10⁴km² (above mild erosion), accounting for 32.81% of the total area of the basin. The areas with increased soil erosion intensity are scattered in Ya'an city, Liangsha, Ganzi and the south of Aba in the southwest of the basin. Therefore, under the low emission forcing scenario of SSP1-2.6 in the future, it is necessary to focus on preventing further increase of soil erosion intensity in the corresponding areas mentioned above. As far as the soil erosion grade change under SSP1-2.6 scenario is concerned, the soil erosion grade of 91.76% area of Min-Tuo River will remain stable in the future, 7.13% area will be aggravated, and 1.11% area will be alleviated. Under SSP2-4.5 scenario, the spatial characteristics of soil erosion intensity in Min-Tuo River are shown in Fig. 8 (b). SSP2-4.5 medium emission scenario, the soil erosion grade distribution tends to be consistent with SSP1-2.6 scenario, but in the future, the soil erosion intensity grade is different from SSP1-2.6 scenario, specifically reflected in: 91.89% of the Min-Tuo River soil erosion grade will be stable, 7.00% of the area soil erosion degree will be aggravated, erosion degree will be slightly reduced. Similarly, under SSP5-8.5 scenario, the spatial characteristics of soil erosion intensity of Min-Tuo River are similar to those of SSP1-2.6 and SSP2-4.5 scenarios, as shown in Fig. 8 (c), except that the change of soil erosion grade is different from the other two models. The differences are reflected in: In future SSP5-8.5 high emission scenario, the area of soil erosion grade increase in Min-Tuo River will further increase on the basis of SSP2-4.5 scenario, reaching 12.50%, in which the area of soil erosion intensity increase sharply is the same as that in SSP2-4.5 scenario, and the soil erosion grade in 91.39% area of Min-Tuo River will remain stable in the future, and the soil erosion degree will increase by 0.52% compared with SSP2-4.5 scenario. Therefore, it is necessary to identify potential areas with abrupt increase in soil erosion intensity rate under future high emission scenarios, and to prevent slight soil erosion areas from transforming into severe soil erosion areas. 4. Discussion This study demonstrates the dynamic changes in soil erosion patterns in the Min-Tuo River basin over a 10-year period and forecasts future soil erosion using the recently developed CMIP6 model for approximately 30 years. The selected CMIP6 model is part of the World Climate Research Programme, and the approved climate change simulation and estimation data ensure the scientific credibility of the predictions. 4.1. Soil erosion aggravation in Min-Tuo River in the future The severity of soil erosion in the Min-Tuo River basin is expected to intensify in the future, leading to increasingly serious issues of soil and water loss. The relationship between future climate and soil erosion is closely intertwined, with climate factors directly or indirectly influencing soil erosion. Among these factors, precipitation has the most significant impact on soil erosion and forms its foundation. This study considers future precipitation conditions, which are projected to increase under three scenarios compared to historical periods, especially with intensified and prolonged precipitation events, including heavy rainfall, which will exacerbate soil and water loss (Yao et al., 2023 ). Additionally, as future temperatures are expected to rise continuously (Song & Yan, 2022 ), leading to high-temperature droughts and extreme heatwaves, soil cracking will occur, resulting in serious consequences such as debris flows when rainfall infiltrates. Addressing the increasingly severe soil erosion caused by future high-temperature heavy rainfall is a matter of significant concern and requires attention and resolution. 4.2. Strengthening soil and water conservation measures in farmland Under future climate scenarios, the increased magnitude of precipitation in the basin will exacerbate soil erosion, and changes in land use may further deteriorate soil erosion in certain areas. Considering the compounded effects of land use change, a significant conversion of grassland to cropland may potentially worsen soil erosion. Therefore, in the process of soil and water erosion control in the future Min-Tuo River basin, attention should be paid to maintaining forest and grassland coverage in the watershed and ensuring stable growth of vegetation cover to mitigate or alleviate the current soil and water loss situation. In the context of increased rainfall in the future, it is advisable to increase the proportion of grassland appropriately, leveraging the role of maintaining surface cover to mitigate soil and water loss while reducing ecological water consumption in the region to safeguard water resources. Furthermore, to prevent the worsening of soil and water loss, the conversion of land from forests, grasslands, and wetlands to cropland should be avoided (Fang, 2021 ). It is recommended to consolidate the achievements of returning farmland to forest and grassland projects, continue to maintain stable growth of vegetation cover to mitigate future soil erosion caused by surface cover factors, and strengthen comprehensive management of small watersheds with slope farmland treatment and accompanying measures for slope surface water systems to enhance watershed water conservation capacity in the upstream areas of the basin. 4.3. Uncertainty of simulation results The simulation results exhibit a certain degree of uncertainty, stemming from both the uncertainties in climate simulation and soil erosion simulation. Firstly, the uncertainty in climate simulation arises from the inability to access high-resolution climate simulation data, leading to reliance on downscaled climate simulation data for the study area. Enhancing the resolution of climate simulation in the future could improve the accuracy of simulations. Secondly, the uncertainty in soil erosion simulation is intrinsic to the soil erosion models themselves, affecting the precision of future soil erosion simulations. In this study, the classic RUSLE model was employed to calculate the soil erosion modulus in the Min-Tuo River basin. Compared to other methods for assessing soil loss, the RUSLE model tends to overestimate erosion due to its omission of the sedimentation process in soil erosion calculations, focusing solely on soil erosion quantity. Therefore, it is necessary to consider how to adjust the model to better suit the soil erosion conditions in the Min-Tuo River basin. In estimating the soil erosion modulus, it is essential to select more reasonable and applicable methods for calculating the factors of the RUSLE model based on actual conditions in the study area. Specifically, the factors are determined based on land use type, but the allocation values are somewhat subjective due to the lack of consideration for different soil and water conservation measures for different land use categories. Implementing appropriate land management interventions to reverse the trends in land use/land cover change and soil erosion in the study watershed is crucial. Additionally, precision in the calibration and calculation of soil erosion factors needs further improvement. For instance, in calculating the rainfall erosivity factor, using daily or even sub-daily rainfall data when available can enhance accuracy. Moreover, for the calculation of soil erodibility values, field sampling to determine basic soil properties and subsequent verification can enhance precision. Furthermore, while this study only considered land use type as a human factor, factors such as population density and GDP could also be included to further investigate the impact of socio-economic development on soil erosion in the Min-Tuo River basin, enriching the research content. Regarding the prediction of future soil erosion in the Min-Tuo River basin, this study only calculated the soil erosion modulus and severity at the annual temporal scale, without elucidating the seasonal or monthly variations in future soil erosion conditions. Thus, the seasonal and monthly variations in soil erosion severity under future scenarios remain unclear. Future research could explore the spatiotemporal distribution characteristics of soil erosion at seasonal and monthly scales. 5. Conclusions Based on the prediction of future precipitation and land use types of Min-Tuo River, this study predicts the soil erosion modulus in the future period (2021 ~ 2050) by using the CMIP6 and scenario data published by the World Climate Research Program, and expounds its temporal and spatial evolution characteristics as follows: The analysis of soil erosion in the study area shows that the soil erosion types from 2000 to 2020 are mainly slight erosion, light erosion and moderate erosion, and the area of light erosion decreases especially from 69.72% in 2000 to 33.52% in 2020. The area of soil erosion reduction is less than that of soil erosion intensification in the study area from 2000 to 2010 and 2010 to 2020. It indicates that soil erosion will be aggravated gradually from 2000 to 2020; Six CMIP6 models and three different emission scenarios were selected to predict the future precipitation of Min-Tuo River. The prediction results showed that the future precipitation of Min-Tuo River fluctuated and increased under three scenarios, and the fluctuation amplitude and increment of future precipitation of Min-Tuo River were the largest under SSP5-8.5 scenario; As far as the future land use types are concerned, the top three categories are: cultivated land to forest land, forest land to cultivated land and grassland to cultivated land, with the conversion areas of 5339km², 40874km² and 27393km² respectively. Under the future climate change scenario, the conversion area of cultivated land will increase, but if the agricultural management activities are improper, it is likely to increase the risk of soil erosion in Min-Tuo River; Compared with the present stage, the area of soil erosion in Min-Tuo River will increase under the future three climate scenarios. Comparatively speaking, the annual average soil erosion modulus and total amount of Min-Tuo River will be the largest under SSP5-8.5 low emission scenario, with the variation between 446.58 ~ 447.74 t/(km²·a) and 7.28 ~ 2.30×10⁸t respectively, followed by SSP1-2.6 scenario, and the soil erosion degree of Min-Tuo River will be the lightest under SSP2-4.5 scenari. Spatially, high-intensity soil erosion will occur in Ya'an, Liangshan, Ganzi and the south of Aba in the southwest of the basin under three different scenarios in the future The stability of soil erosion grade in SSP2-4.5 scenario is higher than that in SSP1-2.6 scenario and SSP5-8.5 scenario; Comparing the changes of soil erosion grade under three scenarios, we can see that the risk of soil erosion degree will become larger and larger with the emission grade of greenhouse gases in the future. Declarations Credit authorship contribution statement Nan Jiang: Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing–original draft. Fuquan Ni: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing–review & editing. Yu Deng: Data curation, Formal analysis, Investigation, Methodology, Resources. Mingyan Wu: Formal analysis, Investigation, Methodology, Resources, Software. Mengyu Zhu: Formal analysis, Investigation, Methodology, Resources, Software. Yuxuan Wang: Data curation, Formal analysis, Investigation, Methodology, Software. Huazhun Ren: Formal analysis, Investigation, Methodology, Resources, Software. Ziying Yue: Data curation, Formal analysis, Investigation, Methodology, Software. Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Data availability The data used in this study is as follows: the parameters for the RULSE equation are sourced from Table 1, where the rainfall data is obtained from National Geographic Resource Science SubCenter, National Earth System Science Data Center, National Science & Technology Infrastructure of China, it can be found by the DOI:10.12041/geodata.269728669112188.ver1.db of the data used; Soil data is sourced from the World Soil Database, with the website address http://www.fao.org/soils-portal/so, after entering the URL, select Soil Maps and Databases in the Date Hub to download the relevant data; Terrain data is derived from the global 30 m resolution DEM data set, which is derived from the CGIAR-CSI SRTM Elevation Database(http://srtm.csi.cgiar.org/download); The NDVI utilizes NASA's MODIS satellite data, visit NASA's MODIS data download website (https://modis.gsfc.nasa.gov/data/dataprod/mod13.php) and select the NDVI data required; Land use / cover data can be accessed through http://www.globeland30.org, and data download requires registration and application submission, this data is prepared for estimating the P factor, and the value of the P factor is ultimately determined by empirical methods based on relevant literature of the study area(Peng et al. 2017). The data for CMIP6 is obtained from https://esgf-node.llnl.gov/search/cmip6/, the models and scenarios used are listed in Table 2 and Table 3, respectively. Funding The authors declare that no funds, grants, or other support were received during the preparation of this manuscript. References Abijith, D., & Saravanan, S. (2021). Assessment of land use and land cover change detection and prediction using remote sensing and CA Markov in the northern coastal districts of Tamil Nadu, India. Environmental Science and Pollution Research, 29(57) , 86055-86067. https://doi.org/10.1007/s11356-021-15782-6 Aqil, T., Jianguo, Y., & Faisal, M. (2022). Land change modeler and CA-Markov chain analysis for land use land cover change using satellite data of Peshawar, Pakistan. Physics and Chemistry of the Earth 128, 103286. https://doi.org/10.1016/j.pce.2022.103286 Bekele, T. (2019). Effect of Land Use and Land Cover Changes on Soil Erosion in Ethiopia. Peertechz Publications Private Limited, 20(1), 79-91. https://doi.org/10.1080/15715124.2020.1767636 Cai, C. F., Ding, S. W., Shi, Z. H., Huang, L., & Zhang, G. Y. (2000). Study of applying USLE and geographical information system IDRISI to predict soil erosion in small watershed. Journal of Soil and Water Conservation, 14, 19–24. doi: 10.13870/j.cnki.stbcxb.2000.02.005 Chen, J., Brissette, F. O. P., & Leconte, R. (2011). Uncertainty of downscaling method in quantifying the impact of climate change on hydrology. Journal of Hydrology 401(3-4), 190-202, https://doi.org/10.1016/j.jhydrol.2011.02.020 Deng, B.,& Shi, L. (2022). Remote Sensing Quantitative Research on Soil Erosion in the Upper Reaches of the Minjiang River. Frontiers in Earth Science 10 , 930535, https://doi.org/10.3389/feart.2022.930535 Woodward, D. E. (1999). Method to predict cropland ephemeral gully erosion - ScienceDirect. Catena 37(3-4), 393-399. https://doi.org/10.1016/S0341-8162(99)00028-4 Fan, L., Lu, C., & Chen Z. (2012). A Review of EPIC Model and Its Applications. Progress in Geography, 31(05), 584-592. doi: 10.11820/dlkxjz.2012.05.006 Suresh, B. S. (2023). Adaptation in the Scientific Method: An Outline for Mutual Learning and Knowledge Co-production in Climate Science. Science, Technology and Society 28(4), 621-638. https://doi.org/10.1177/09717218231197448 Guo, B., Tao, H., Liu, B., & Jiang, L. (2012). Characteristics and analysis of soil erosion in Li country after Wenchuan earthquake based on GIS and USLE. Transactions of the Chinese Society of Agricultural Engineering , 28(14), 118-126. doi: 10.3969/j.issn.1002-6819.2012.14.01 Hairsine, P. B.,& Rose, C. W. (1992). Modeling water erosion due to overland flow using physical principles: 2. Rill flow. Water Resources Research 28(1), 245–250. https://doi.org/10.1029/91WR02381 Fang, H. (2021). Changes in Cultivated Land Area and Associated Soil and SOC Losses in Northeastern China: The Role of Land Use Policies. International Journal of Environmental Research and Public Health 18(21) , 11314. https://doi.org/10.3390/ijerph182111314 Hessel, R., Jetten, V. G., Liu, B., & Qiu, Y. (2011). Evaluating effects of soil and water management and land use change on the Loess Plateau of China using LISEM. Handbook of Erosion Modelling. Chichester, West Sussex: Wiley-Blackwell, 223-248. https://doi.org/10.1002/9781444328455.ch12 Humphrey, O. S., Osano, O., Aura, C. M., Marriott, A. L., Dowell, S. M., Blake, W. H., et al. (2022). Evaluating spatio-temporal soil erosion dynamics in the Winam Gulf catchment, Kenya for enhanced decision making in the land-lake interface. Science of the Total Environment, 815, 151975. https://doi.org/10.1016/j.scitotenv.2021.151975 Hu, J., Ma, J., Nie, C., Xue, L., Zhang, Y., Ni, F., et al. (2020). Attribution Analysis of Runoff change in Min-tuo River Basin based on SWAT model simulations, china. Scientific reports, 10(1), 2900. doi: 10.1038/s41598-020-59659-z Zhang, J., Li, H., Zhang, X., Yu, H., Liang, F., Qin, Z., et al. (2020). Sensitivity evaluation of soil erosion based on land use types: A case study of Minjiang River Basin. Journal of Intelligent & Fuzzy Systems, 38(5), 5697-5705. https://doi.org/10.3233/JIFS-179658 Keno, B., & Suryabhagavan, K. V. (2014). Multitemporal remote sensing of landscape dynamics and pattern change in Dire district, southern Ethiopia. Journal of Geomatics, 8(2), 189-194. doi: 10.4172/2157-7617.1000226 Liu, J., Ni, F., Deng, Y., Li, L., Wang, Z., Zhou, D., et al. (2019). Risk assessment of soil erosion in the upper reaches of Minjiang Rive. South-to-North Water Transfers and Water Science & Technology, 17(01), 105-112. doi: 10.13476/j.cnki.nsbdqk.2019.0015. Lu, K., Arshad, M., Ma, X., Ullah, I., Wang, J., & Shao, W. (2022). Evaluating observed and future spatiotemporal changes in precipitation and temperature across China based on CMIP6‐GCMs. International Journal of Climatology, 42(15), 7703-7729. https://doi.org/10.1002/joc.7673 Ministry of Water Resources of the People's Republic of China. (2008). SL190-2007 , classification and grading standards for soil erosion. Beijing: China Water & Power Press. Ministry of Water Resources of the People's Republic of China. (2018). SL773-2018 , guidance on estimation of soil loss in production and construction projects. Beijing: China Water & Power Press. Mirzabaev, A., Goedecke, J., Dubovyk, O., Djanibekov, U., Le, Q. B., & Aw-Hassan, A. (2016). Economics of land degradation in Central Asia. Economics of land degradation and improvement–A global assessment for sustainable development, 261-290. doi: 10.1007/978-3-319-19168-3_10 Nearing, M. A., Simanton, J. R., Norton, L. D., Bulygin, S. J., & Stone, J. (1999). Soil erosion by surface water flow on a stony, semiarid hillslope. Earth Surface Processes and Landforms: The Journal of the British Geomorphological Research Group, 24(8), 677-686. https://doi.org/10.1002/(SICI)1096-9837(199908)24:83.0.CO;2-1 Peng, J., Li, D. D., & Zhang, Y. Q. (2007). Analysis of spatial characteristics of soil erosion in mountain areas of northwestern Yunnan based on GIS and RUSLE. Journal of Mountain Science, 25(5), 548-556. ISSN:1008-2186/CN:51-1516 Peng, S., Ding, Y., Wen, Z., Chen, Y., Cao, Y., & Ren, J. (2017). Spatiotemporal change and trend analysis of potential evapotranspiration over the Loess Plateau of China during 2011–2100. Agricultural and forest meteorology, 233, 183-194. https://doi.org/10.1016/j.agrformet.2016.11.129 Wang, Q., & Wang, H. (2022). An integrated approach of logistic-MCE-CA-Markov to predict the land use structure and their micro-spatial characteristics analysis in Wuhan metropolitan area, Central China. Environmental Science and Pollution Research, 1-24. doi: 10.1007/s11356-021-17750-6 Rao, L., Xu, Y., Hu, J., Wang, Z., Yao, W., & Shen, Z. (2020). Study on Soil Erodibility Factoe K on Soil Cover Area of Piaha Sandstone Region. Journal of Basic Science and Engineering, 28(04), 763-773. doi: 10.16058/j.issn.1005-0930.2020.04.002 Renard, K. G., & Ferreira, V. A. (1993). RUSLE model description and database sensitivity. Journal of environmental quality, 22(3), 458-466. https://doi.org/10.2134/jeq1993.00472425002200030009x Renard, K. G., Yoder, D. C., Lightle, D. T., & Dabney, S. M. (2011). Univeral soil loss equation and revised universal soil loss equation. Handbook of Erosion Modelling, edited by: Morgan R PC and Nearing, MA, 137-67. https://doi.org/10.1002/9781444328455.ch8 De Roo, A. P. J. (1996). The LISEM project: an introduction. Hydrological Processes, 10(8), 1021-1025. https://doi.org/10.1002/(SICI)1099-1085(199608)10:83.0.CO;2-I Rymszewicz, A., Mockler, E., O'Sullivan, J., Bruen, M., Turner, J., Conroy, E., et al. (2015). Assessing the applicability of the Revised Universal Soil Loss Equation (RUSLE) to Irish Catchments. Proceedings of the International Association of Hydrological Sciences, 367, 99-105. https://doi.org/10.5194/piahs-367-99-2015 Piao, S., Zhang, X., Chen, A., Liu, Q., Lian, X., Wang, X., et al. (2019). The impacts of climate extremes on the terrestrial carbon cycle: A review. Science China Earth Sciences, 62, 1551-1563. doi: 10.1007/s11430-018-9363-5 Ma, S., Wang, L. J., Wang, H. Y., Zhao, Y. G., & Jiang, J. (2023). Impacts of land use/land cover and soil property changes on soil erosion in the black soil region, China. Journal of Environmental Management, 328, 117024. https://doi.org/10.1016/j.jenvman.2022.117024 Song, S., & Yan, X. (2022). Evaluation of events of extreme temperature change between neighboring days in CMIP6 models over China. Theoretical and Applied Climatology, 150(1-2), 53-72. doi: 10.21203/rs.3.rs-601223/v1 Su, Y., Wu, X., Guo, C., Chen, X., Liu, L., & Peng, S. (2022). The spatiotemporal pattern and trend of annual precipitation over China from 2001 to 2100. Journal of Lanzhou University (Natural Sciences), 58(05), 641-649. doi: 10.13885/j.issn.0455-2059.2022.05.009 Wang, Y., Li, H., Wang, H., Sun, B., & Chen, H. (2021). Evaluation of CMIP6 model simulations of extreme precipitation in China and comparison with CMIP5. Acta Meteorologica Sinica, 79(03), 369-386. doi: 10.11676/qxxb2021.031 Wang, Z., & Su, Y. (2020). Assessment of soil erosion in the Qinba mountains of the Southern Shaanxi Province in China using the RUSLE model. Sustainability, 12(5), 1733. https://doi.org/10.3390/su12051733 Williams, J. R. (1990). The erosion-productivity impact calculator (EPIC) model: a case history. Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences, 329(1255), 421-428. https://doi.org/10.1098/rstb.1990.0184 Wischmeier, W. H., & Smith, D. D. (1958). Rainfall energy and its relationship to soil loss. Eos, Transactions American Geophysical Union, 39(2), 285-291. doi: 10.1029/TR039I002P00285 Wischmeier, W. H., & Smith, D. D. (1978). Predicting rainfall erosion losses: a guide to conservation planning (No. 537). Department of Agriculture, Science and Education Administration. http://topsoil.nserl.purdue.edu/usle/AH_537.pdf Wu, J., Xia, J., Zeng, S., Liu, X., & Fan, D. (2023). Evaluation of the Performance of CMIP6 Models andFuture Changes Over the Yangtze River Basin. Resources and Environment in the Yangtze Basin , 32(01), 137-150, doi:10.11870/cjlyzyyhj202301013 Yang, J., & Huang, X. (2021). The 30 m annual land cover dataset and its dynamics in China from 1990 to 2019. Earth System Science Data , 13(8), 3907-3925, https://doi.org/10.5194/essd-13-3907-2021 Yang J., Dong J., Xiao X., Dai J. Wu C., Xia J., et al. (2019). Divergent shifts in peak photosynthesis timing of temperate and alpine grasslands in China. Remote Sensing of Environment , 233, 111395, https://doi.org/10.1016/j.rse.2019.111395 Wei, Y., Sun, S., Liang, D., & Jia, Z. (2022). Spatial–temporal variations of NDVI and its response to climate in China from 2001 to 2020. International Journal of Digital Earth , 15(1), 1463-1484, https://doi.org/10.1080/17538947.2022.2116118 Yao, Y., Li, J., Jiang, Y., & Huang, G. (2023). Evaluating the response and adaptation of urban stormwater systems to changed rainfall with the CMIP6 projections. Journal of Environmental Management , 347, 119135, https://doi.org/10.1016/j.jenvman.2023.119135 Zhong, X., Zhang, S., Wu, R., Jing, Y., Men, L., & Zhou, T. (2022). Analysis of Dynamic Changes and Driving Forces of Soil Erosion in Tuojiang River Basin. Research of Soil and Water Conservation , 29(02), 43-49+56, doi:10.13869/j.cnki.rswc.2022.02.003 Supplementary Files Attachment.docx Highlights.docx Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4952297","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":347813406,"identity":"f952da02-e864-4c0c-a6e8-6a6b1739adf1","order_by":0,"name":"Nan Jiang","email":"","orcid":"","institution":"Sichuan Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Nan","middleName":"","lastName":"Jiang","suffix":""},{"id":347813407,"identity":"a6ce21f1-9dbf-46f6-a9c0-34101896195d","order_by":1,"name":"Fuquan Ni","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAvklEQVRIiWNgGAWjYDAC5gMMBz5U/JNjY28/QKQWtgTGgzPOHDDm4zmTQLQW5sO8LQcS50k4GBCnw5yNecNh3oY76W0SDAkMPyq2EdZi2cZWcHDujme5bdKNBxh7ztwmrMXgfo/BgbdnmHPbZA4kMDO2EaPlGI/BAd425nQ2iQQD4rUc5G07nECKFqBfZpxJM2wDBvJB4vxyjHnzhw8VNvLy7e0HH/yoIEILSBecdYAo9ShaRsEoGAWjYBRgBQCsXkOLMbZ4mwAAAABJRU5ErkJggg==","orcid":"https://orcid.org/0000-0002-6362-3008","institution":"Sichuan Agricultural University","correspondingAuthor":true,"prefix":"","firstName":"Fuquan","middleName":"","lastName":"Ni","suffix":""},{"id":347813408,"identity":"87c85029-e761-4663-97e5-c58145a06d05","order_by":2,"name":"Yu Deng","email":"","orcid":"","institution":"Sichuan Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Yu","middleName":"","lastName":"Deng","suffix":""},{"id":347813409,"identity":"3db99f4f-e998-4a94-82e2-d1aa2f8c6805","order_by":3,"name":"Mingyan Wu","email":"","orcid":"","institution":"Sichuan Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Mingyan","middleName":"","lastName":"Wu","suffix":""},{"id":347813410,"identity":"750c6446-2d15-4ed6-bf93-e3ac412fd778","order_by":4,"name":"Mengyu Zhu","email":"","orcid":"","institution":"Sichuan Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Mengyu","middleName":"","lastName":"Zhu","suffix":""},{"id":347813411,"identity":"afe53c10-5581-42ef-ac5f-ccda23db63a6","order_by":5,"name":"Yuxuan Wang","email":"","orcid":"","institution":"Sichuan Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Yuxuan","middleName":"","lastName":"Wang","suffix":""},{"id":347813412,"identity":"d9f76246-77d7-417c-87d1-3214678e7c50","order_by":6,"name":"Huazhun Ren","email":"","orcid":"","institution":"Changjiang Water Resources Commission","correspondingAuthor":false,"prefix":"","firstName":"Huazhun","middleName":"","lastName":"Ren","suffix":""},{"id":347813413,"identity":"8a384ffa-548d-475c-bdbf-e35b6ac05d1a","order_by":7,"name":"Ziying Yue","email":"","orcid":"","institution":"Sichuan Agricultural University","correspondingAuthor":false,"prefix":"","firstName":"Ziying","middleName":"","lastName":"Yue","suffix":""}],"badges":[],"createdAt":"2024-08-21 14:25:09","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4952297/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4952297/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":66940404,"identity":"7fdd87f9-17e4-45f3-bdf1-a746e55dab36","added_by":"auto","created_at":"2024-10-18 08:48:47","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":74084,"visible":true,"origin":"","legend":"\u003cp\u003eMin-Tuo River Basin is a primary tributary of the upper Yangtze River, located between 99° to 106°E and 28° to 34°N. The basin area is approximately 16.30×104km², with elevations ranging from 200 to 6511 meters.\u003c/p\u003e","description":"","filename":"1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/847ba77f4e445a4dd110747c.jpg"},{"id":66940411,"identity":"207f3561-d76f-4366-b8bf-9ee9b0d5380b","added_by":"auto","created_at":"2024-10-18 08:48:48","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":45965,"visible":true,"origin":"","legend":"\u003cp\u003eThe Taylor diagram, the closer the model prediction is to the observation point (Obs), the closer the model prediction is to the actual value. It can be seen that the simulation results of each model differ significantly, but the correlation coefficient of the MME result is 0.99, which is higher than the other six models\u003c/p\u003e","description":"","filename":"2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/54ab631c9696fc36142f11a1.jpg"},{"id":66940409,"identity":"094334ac-a050-4591-9a27-6b9d8b2d25b5","added_by":"auto","created_at":"2024-10-18 08:48:48","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":114332,"visible":true,"origin":"","legend":"\u003cp\u003eBy reclassifying the soil erosion modulus, the spatial distribution characteristics of soil erosion intensity in the Minjiang River Basin can be obtained.\u003c/p\u003e","description":"","filename":"3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/679af3aea1a05487ca296294.jpg"},{"id":66940407,"identity":"12c7ce33-9a4c-4d44-af59-65e0b959257d","added_by":"auto","created_at":"2024-10-18 08:48:47","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":44588,"visible":true,"origin":"","legend":"\u003cp\u003eFrom 2000 to 2020, except for a slight decrease in the slight erosion grade and an unclear change in the extremely severe erosion grade, all other erosion grades increased. The decrease in the slight erosion grade and the increase in the moderate erosion grade were particularly significant from 2010 to 2020. Erosion showed an intensifying trend from 2000 to 2010, and this trend continued to intensify from 2010 to 2020. In comparison to the period from 2000 to 2010, the intensification of erosion was significantly greater from 2010 to 2020, while the alleviation of erosion decreased. Overall, from 2000 to 2020, erosion exhibited an intensifying trend.\u003c/p\u003e","description":"","filename":"4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/050945a4c1b29b27a105dee7.jpg"},{"id":66941612,"identity":"6589f643-9d7a-4a64-9355-10f564436cad","added_by":"auto","created_at":"2024-10-18 08:56:47","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":34210,"visible":true,"origin":"","legend":"\u003cp\u003eFrom a temporal perspective, under the three scenarios, the future precipitation for the Min-Tuo River shows a fluctuating upward trend relative to the historical period (1981-2014). Among them, under the SSP5-8.5 scenario, the future precipitation for the Min-Tuo River exhibits the largest fluctuation amplitude and the greatest increase, followed by the SSP1-2.6 scenario, with the smallest precipitation variation occurring under the SSP2-4.5 scenario. For the Min-Tuo River, under the SSP5-8.5 scenario, future precipitation is expected to increase significantly, and the instability of precipitation factors will also strengthen.\u003c/p\u003e","description":"","filename":"5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/ff4bc921396712176a83ef15.jpg"},{"id":66940412,"identity":"d9dfa549-0b60-4378-ac37-7943d01b7bf0","added_by":"auto","created_at":"2024-10-18 08:48:48","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":78525,"visible":true,"origin":"","legend":"\u003cp\u003eThe various land use types in the Min-Tuo River basin undergo dynamic changes in the future period. In terms of area ranking, the top three land use conversions in the future are: conversion from cropland to forestland, forestland to cropland, and grassland to cropland, with conversion areas of 5339 km², 40874 km², and 27393 km², respectively. The conversion areas for other land use types are relatively small.\u003c/p\u003e","description":"","filename":"6.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/f2d8be7f28f2c09f38b9e13d.jpg"},{"id":66940406,"identity":"88cfc27a-92ae-4529-a703-c34e4265a1c8","added_by":"auto","created_at":"2024-10-18 08:48:47","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":55186,"visible":true,"origin":"","legend":"\u003cp\u003eIn all three scenarios, the ranking of the proportion of each soil erosion level in the Min-TuoRiver basin is as follows: slight \u0026gt; light \u0026gt; moderate \u0026gt; intense \u0026gt; very intense \u0026gt; severe. Looking at the forecast years from 2030 to 2050, there is a slight reduction in slight and light erosion, a slight increase in moderate and intense erosion, and a generally stable situation for very intense and severe erosion. Compared to the current stage, the area affected by soil erosion in the Min-Tuo River basin is increasing.\u003c/p\u003e","description":"","filename":"7.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/f07477baaeba5f8d801af3f2.jpg"},{"id":66941617,"identity":"530ed5f1-e46f-459d-a7a5-1834e9b83c42","added_by":"auto","created_at":"2024-10-18 08:56:48","extension":"jpg","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":83066,"visible":true,"origin":"","legend":"\u003cp\u003eIn terms of the area proportion of soil erosion grade remaining unchanged in the future under the three different scenarios in the future, the stability degree of soil erosion grade under SSP2-4.5 scenario is higher than that under the other two scenarios, and the stability order of soil erosion grade in the future is: SSP2-4.5\u0026gt; SSP1-2.6\u0026gt; SSP5-8.5; The proportion of soil erosion increased in the future showed the order of SSP2-4.5\u0026gt; SSP1-2.6\u0026gt; SSP5-8.5, indicating that the risk of soil erosion increased in SSP5-8.5 scenario was the greatest with the increase of CO2 emission level in the future.\u003c/p\u003e","description":"","filename":"8.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/bad03c8fb40fd344bcaf8821.jpg"},{"id":72205148,"identity":"6d32a2be-b494-4470-a443-887e508aea43","added_by":"auto","created_at":"2024-12-23 16:31:36","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1433289,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/046c34c9-baf2-42ef-94a4-2ad68b2c43c1.pdf"},{"id":66940405,"identity":"4e60ff15-3876-49bf-b717-f68693d5d5b4","added_by":"auto","created_at":"2024-10-18 08:48:47","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":280956,"visible":true,"origin":"","legend":"","description":"","filename":"Attachment.docx","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/7339b74523a4cdd5d0153a67.docx"},{"id":66941616,"identity":"fc43cf95-d19f-4b32-b148-0c667faf00c8","added_by":"auto","created_at":"2024-10-18 08:56:48","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":11780,"visible":true,"origin":"","legend":"","description":"","filename":"Highlights.docx","url":"https://assets-eu.researchsquare.com/files/rs-4952297/v1/f91ee67f55e7b951d7bd42f8.docx"}],"financialInterests":"","formattedTitle":"High greenhouse gas emission scenarios increase soil erosion risk: A case study of the Min-Tuo River Basin, China","fulltext":[{"header":"1. Introduction","content":"\u003cp\u003eSoil erosion is a global environmental threat (Keno \u0026amp; Suryabhagavan, 2015). It is the main cause of soil degradation and is related to soil compaction, soil structure and organic matter loss, and soil acidity (Hessel et al., \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Almost no form of land on our planet is free from soil erosion (Bekele, \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Approximately 30% of the global land area and 80% of ag-ricultural land experience soil erosion, with around 3\u0026nbsp;billion people living on degraded land (Mirzabaev et al., \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Soil erosion is particularly severe in developing regions with high population density and ecological vulnerability (Renard et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). The characteristics of the Min-Tuo River Basin align with these conditions. Additionally, severe forest resource de-struction, the 2008 Wenchuan earthquake (Guo et al., \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2012\u003c/span\u003e), and extreme rainfall events (Piao et al., \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2019\u003c/span\u003e) have further worsened the region's ability to retain soil and resulted in local soil desertification and lithification, leading to serious soil and water loss in the area. Therefore, quantifying soil loss and predicting future soil erosion is of great significance as it can provide scientific basis for watershed soil and water loss management and macro-scale decision-making for soil and water conservation, thus promoting ecological environment protection and restoration in the watershed.\u003c/p\u003e \u003cp\u003eThe issue of soil erosion has become a subject of study for numerous scholars. Among these, the use of soil erosion models has proven to be the most widely effective method for researching soil erosion. Several scholars have successively developed models such as the USLE (Wischmeier \u0026amp; Smith, \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e1958\u003c/span\u003e), WEPP (Woodward, \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e1999\u003c/span\u003e), and LISEM (De Roo, \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e1996\u003c/span\u003e), with the USLE model being the most widely utilized. Subsequently, scholars have made supplementary revisions to the USLE model. In 1997, the United States Department of Agriculture proposed the RUSLE model (Renard \u0026amp; Ferreira, \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e1993\u003c/span\u003e), which comprehensively reflects the influencing factors of soil erosion and has been widely used (Rymszewicz et al., \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Ma et al., \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Wang \u0026amp; Su, \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eWith the development of GIS and RS technologies, their combined use allows for the comprehensive exploration of spatial information and patterns related to soil erosion, thereby providing a technological possibility for studying the spatiotemporal variations in soil erosion. Over the past century, global climate has been trending significantly warmer (Suresh, \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Hairsine \u0026amp; Rose, \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e1992\u003c/span\u003e), leading to increased regional differences in precipitation. Rainfall is one of the important factors influencing soil erosion, and global climate change can cause changes in regional precipitation, thereby affecting soil erosion conditions.Soil erosion is a crucial aspect of sustainable development. Analyzing soil erosion factors and processes using soil erosion models can provide a theoretical basis for formulating regional soil and water conservation measures (Nearing et al., \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e1999\u003c/span\u003e). Therefore, this study is based on the RUSLE model to investigate the current status of soil erosion in the Min-Tuo River Basin and predict future soil erosion conditions based on the CMIP6 model.\u003c/p\u003e \u003cp\u003eIn the past, many related studies have studied the current situation of soil erosion in the Min-Tuo River Basin. For instance, Zhong et al. (\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) used the RUSLE model to assess soil erosion in the Tuojiang River from 2000 to 2018 and explored the spatiotemporal dynamic evolution patterns. Deng \u0026amp; Shi (\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) quantitatively evaluated soil erosion and the spatial distribution of soil and water conservation ecosystem services in the Min-Tuo River Basin using the USLE model. Zhang et al. (\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) generated a soil erosion modulus map for the Min-Tuo River Basin using RUSLE based on supervised classification of the current land use map, indicating that the basin is mainly characterized by mild and moderate erosion areas. Liu et al. (\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2019\u003c/span\u003e) investigated the erosion status and characteristics of the upper reaches of the Min-Tuo River Basin and explored erosion risk in the region using a GIS weighted overlay analysis method. However, there is currently very little research on the overall soil and water loss situation in the entire Min-Tuo River Basin after 2010.\u003c/p\u003e \u003cp\u003eGiven the background mentioned above, this study focuses on the Min-Tuo River Basin as the main area. Through the use of remote sensing images, precipitation records, DEM, soil types, vegetation types, land use, and other data, the study employed the RUSLE model to analyze soil erosion in the Min-Tuo River Basin from 2000 to 2020. Furthermore, based on the CMIP6 model, six models and three different emission scenarios were selected. After bias correction and downscaling using the multi-model ensemble (MME) approach, the study predicted the precipitation in the basin from 2021 to 2050. Additionally, with the assistance of the CA-Markov model, the study forecasted land use types for the years 2030, 2040, and 2050, and subsequently predicted the soil erosion situation in the Min-Tuo River Basin from 2030 to 2050. This research can provide theoretical references for soil erosion control and soil conservation measures in the study area under future climate change scenarios. It can also scientifically predict the soil erosion in large river basins in China at a spatial scale, aiming to provide a scientific and rational research basis for the ecological environment and comprehensive management projects in the Min-Tuo River Basin.\u003c/p\u003e"},{"header":"2. Materials and Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1 Study area\u003c/h2\u003e \u003cp\u003eThe The Min-Tuo River Basin (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e) is mainly composed of four parts: the Minjiang, Dadu River, Qingyijiang, and Tuojiang (Hu et al., \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Particularly, the upper reaches of the Minjiang and the arid valleys of the Dadu River have extremely fragile ecosystems. Under the dual pressure of population and economic development, the soil erosion in the Min-Tuo River Basin is becoming more and more serious. With the acceleration of urbanization, industrialization and agricultural modernization, the expansion of economic scale in the future will have an important impact on soil erosion in the Min-Tuo River Basin. At the same time, with the implementation of national strategies such as the Yangtze River Economic Belt and the Chengdu-Chongqing Economic Circle, higher requirements have been put forward for soil and water conservation in river basins.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.2 RULSE model\u003c/h2\u003e \u003cp\u003eThe In recent years, with the maturity of \"3S\" technology, the difficulty of obtaining geographical information by traditional research methods has been overcome, the efficiency of soil erosion assessment has been improved, and the RUSLE model has been better applied. The RUSLE model expression is as follows:\u003c/p\u003e \u003cp\u003e \u003cspan class=\"InlineEquation\"\u003e \u003cspan class=\"mathinline\"\u003e\\(\\:A=R\\times\\:K\\times\\:L\\times\\:S\\times\\:C\\times\\:\\)\u003c/span\u003e \u003c/span\u003e \u003cem\u003eP\u003c/em\u003e (1)\u003c/p\u003e \u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{A}\\)\u003c/span\u003e\u003c/span\u003e is soil and water conservation quantity (t/hm2\u0026middot;a);\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\:\\text{R}\\)\u003c/span\u003e\u003c/span\u003e is rainfall erosion factor (MJ\u0026middot;mm/hm2\u0026middot;h\u0026middot;a); \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\)\u003c/span\u003e\u003c/span\u003eis soil erodibility factor [t\u0026middot;h/(MJ\u0026middot;mm)];\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\:\\text{L}\\text{S}\\)\u003c/span\u003e\u003c/span\u003e is topographic factor, dimensionless;\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\:\\text{C}\\)\u003c/span\u003e\u003c/span\u003e is vegetation cover factor, dimensionless; P is soil and water conservation measure factor, dimensionless.\u003c/p\u003e \u003cp\u003eThe rainfall data partly comes from the National Science and Technology Infrastructure of China, using the 1km resolution annual precipitation data from 2001 to 2020 in China. This data is derived by synthesizing the yearly cumulative sum based on the monthly precipitation data set at a 1km resolution by Professor Peng Shouzhang from 1901 to 2020, as referenced in Peng et al. (\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2017\u003c/span\u003e) and Su et al. (\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Another part is sourced from the Resource and Environment Data Center of the Chinese Academy of Sciences, which is the spatially interpolated data set of annual mean temperature and annual precipitation since 1980 in China. This data set is based on daily observational data from over 2400 meteorological stations in China, processed through collection, statistics, and spatial interpolation calculations.\u003c/p\u003e \u003cp\u003eThe NDVI data is sourced from the National Ecological Science Data Center of China, as referenced in Yang et al. (\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). The land use/cover data is sourced from the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences. The soil data is sourced from the World Soil Database, where the HWSD (Harmonized World Soil Database) can be queried for soil texture data, including sand, silt, clay, and organic carbon content. The terrain data uses a 30m resolution Digital Elevation Model (DEM) data sourced from the CGIAR-CSI SRTM Elevation Database.\u003c/p\u003e \u003cp\u003eFor the calculation of soil and water conservation, it is necessary to ensure consistency in the resolution and projection coordinate system of each factor when using them collectively. Please refer to Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e for the sources and usages of the relevant data.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eData sources.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"4\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eData classification\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eProcessing\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eResolution\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eData sources\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRainfall data\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEstimated R-value\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1km\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eData Center for Resources and Environment, CAS/National Earth System Science Data Center\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSoil data\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEstimated K value\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1km\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eWorld Soil Database\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTerrain data\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEstimated LS value\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e30m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eCGIAR-CSI SRTM Elevation Database\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNDVI data\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEstimated C value\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e30m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eNational Ecological Science Data Center\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLand use/cover data\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEstimated P value\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e30m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eInstitute of Geographical Sciences and Natural Resources, CAS\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cdiv id=\"Sec5\" class=\"Section3\"\u003e \u003ch2\u003e2.2.1. Rainfall erosivity factor(\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:R\\)\u003c/span\u003e\u003c/span\u003e)\u003c/h2\u003e \u003cp\u003eRainfall erosivity factor is one of the potential forces causing soil erosion, which can be judged by annual rainfall data.Guidelines for Calculation of Soil Loss Amount of Production and Construction Projects (SL773-2018).The formula for calculating \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{R}\\)\u003c/span\u003e\u003c/span\u003e is as follows:\u003cdiv id=\"Equ1\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ1\" name=\"EquationSource\"\u003e\n$$\\:R=0.067{P}_{d}^{1.627}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e2\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:R\\)\u003c/span\u003e\u003c/span\u003e is the erosive factor of annual average rainfall, and the unit is MJ\u0026middot;mm/(hm\u0026sup2;\u0026middot;h); \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{P}_{d}\\)\u003c/span\u003e\u003c/span\u003e is annual average rainfall, unit is mm.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section3\"\u003e \u003ch2\u003e2.2.2. Soil erodibility factor(\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:K\\)\u003c/span\u003e\u003c/span\u003e)\u003c/h2\u003e \u003cp\u003eSoil erodibility factor is one of the essential parameters of RUSLE model, which is used to reflect soil sensitivity to erosion. \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\)\u003c/span\u003e\u003c/span\u003e can be determined by soil texture data by Rao et al.(2020)、Williams (\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e1990\u003c/span\u003e) and Fan et al.(2012). An algorithm for estimating soil erodibility factor in EPIC model is proposed to calculate \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\)\u003c/span\u003e\u003c/span\u003e. The specific formula is as follows:\n\u003c/p\u003e\u003cp\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAjMAAACDCAYAAAB4B5qeAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAACgbSURBVHhe7Z39qxbV+v/H7+9qaj9FSKgHjqgEpRmlgoEPWT8Up4MVEYKSRzsEHUs9PnCIjtn2dBQOWhYFIZVpHpLIPmag4FbBp9hSYmA7RMSftqWeP2B/79dyrru1lzNzz9zPs32/YO2515qZ9XCtp2uutWb2iMEKkRBCCCFESfl/8VEIIYQQopRImRFCCCFEqZEyI4QQQohSI2VGCCGEEKVGyowQQgghSo2UGSGEEEKUGikzQgghhCg1UmaEEEIIUWqkzAghhBCi1EiZEUIIIUSpkTIjhBBCiFIjZUYIIYQQpUbKjGiIY8eORc8++2x07ty5OESI4cGWLVuilStXRr/88kscIoToVqTMiLr5/PPPo1dffTXasWNHdP/998ehopuhzlBAG4VJvhk0K55WsHr16uj555+P5s+fL2VdiJxMmjSp+gCwcOHCaMSIEc75fZ1xqOl9f1CUjt7e3kGqDrdgwYI4dLAahuOaVkL8Y8aMGRwYGIhD6qOvr8+VgTwT37p16+IzyZDeihUrquVcvHhxah56enpS49u9e/fg9OnTU+OxdMiT5Q1/eA1pTJw4sXpNrfz7IEO//vJAepYn8o/80vDbie+S2gbxEF+SLEkzKR4jLR0ccjasnn0Xlt+Pq1NtO4T0qOMk2TRCkbo0suqpv79/SJslz778fQ4cOJDY9vLWpRAhtD/aifVP2qL9tnN+G7Kxs1lImSkpTJoMWv6gxoTMQNeOwZ6BcOfOnbGvPsg7ZaDRgw2kNPI0SNcGczoInSEclG3AT4uLMM5zP1i6lg8gTpx/DWmF1+C3OiDesMNmQZxh3rMgfuRF+UiT+g7bgI+VqxaW77RrOZ+VT9Lx5QLcEw5UxEF4LTrdtpMg7Tx5zwtxFalL4J6seiLc+qTFSZi1YfDDk+o0b10KEUIb8fun3+4gqf/jD9tbvUiZKSk0DLMCMEDRIMKG0ipopOEgWQ9M+sTjD+AMtGkDp6XLU6XB4O3nhSNxEGfeyROIw1fOwnSAuJImAMOUhzxP2MD1WfGFIBff8kM6pJemPFl+siA+ypl1ba1yJ8HEHCq7eeuD6zrVttOgLHkmdPKL4lOLonWZp55CkBn14EPfIK0idZpUl0L45GlPaf2HcNp1o2jPTEn59ttvo0cffTT65ptvoo0bN0avvfaaW+OvBZsabQ0zy2Vx6tQpd5wwYYI71svRo0ejysAfjRs3Lg6JogceeCCqKCSxbyiW7syZM90RpkyZ4o5Xr151R/LEeqwfZxbsH2HttjJZRMuXL49Do6jSMaN//OMf1bVfjh988EH0yiuvOL/PtWvXXD1s3749qkxGLdk/RPrIhTo3LJ3Lly+7Yz1s2rQpWrRoUexrDsgfnnnmGXcsSr1tu5XQzpA/dZ0Ee2psDwBtDz/7BZKopy6L1BPxUweHDx+Ojhw5EofegvAi7bPRuhR3BoyNS5YsiX23Q184dOhQ7BsKe9I+++yz2Fc/UmZKCIM8PPHEE87dd999uRULJoWKElvTZcGAy2TfKAy6aUpH0iZVG+iT7jlx4kT8Kz8obbNnz3adbNq0aXHoLQ4ePBhVnhic4zqO//znP2+bUMjn3Xff7erh+vXr0b333hufaS6mrI0aNcodfZi0khg5cmRUeaqOZsyY4cowduxYp8wWZfz48U6RZGMf8XC0NpjEhg0borVr195WT7RRBj3yQTwMcOHG2kbadju4cOFC/GsoKAhsFt66dauT1Zo1a6LK02p8dij11GVeqF/a6nPPPefkRhtohLS6FMIw5Txp7OMcfZ0HFNqlKcc+c+bMSVV0iiBlpoQcP37cWTRQOlAqmCC6ASZ2Gm4tl+dtmqKD8OjRo+Nf+UF+AwMD0euvv+4Gf3+if//9910HxXFdb2+ve/IOlYFZs2a581zH5IFylFa+0CrGtXRyPwxXlLvuuiv+NRQm2F9//TU6ffq0y+PmzZvdJJs0oGTBq/fE8/PPPzt58SSFosFAFULcyGLp0qVxyO+8++67Lg7i4hqsHH/605/is7fo1radh5s3b0Znz551ytqDDz5YlyKRVpd5sYeVvr4+V+/IMs2aVIusuhTCMOX8nnvucUcfxkQbPwHlOAnaWaNImSkhaLFMKEDjoCFkPSn7hBNqmqtF0gBpE3stx3XAwJ020CaZwk1hSbrHlpuKwhMny0sM+v5T8V/+8pdo27ZtVasAeV6xYkX09ttvO38I17355pvud5qVKLSK0cGZsP0wXBI2MTJhhrA0lwcr5/79++OQ4iAvK+d3333njj60R+RU60keeb3xxhuu7frWmTxtm/rParNmicKhWPjK2/r16124vwSEUobCZvfU88ooZfj666+jTz/91FmuUADS4mlGXdaC/oN8f/vtt1RrUi3y1qUQtWD8ZAm+lUiZKRkMvGfOnKmut9NImBB37drl/LVoxjITyw7koVHmzp3r4vGVk0uXLrkJN4mHH37YHU+ePOmOcOXKFXecPHmyO9ZL0oCdNNkwOaTRykGfyYklIywXhllGpk6d6o55yMp/Xm7cuBH/GgoWKZQP9rjkIZRv3raNnJlkd+7cWW2vXGesWrWqqiR+8sknzupmmJUBpQNoe7Q3lAiuP3DggFvaSiOtnVE/WJ6AOFHWWKpMoll1WYukZay8FK1LcediFhmz0KTBMhSKfhIsQTVMpQOLEsGucarN3/3NWw6EVQafOKS1kE4z0uN+3pTgbRWgTPj9t4hIpzIxxb5br8jieGuE+/ntvxXiwzmL24f47I0nIF3S8d8k4RrutzJa3vy0uIf6MKxu8sqFOP2y1YK0yQNvo5B3ykYejbBtEDdhVk7eSPHP+2S1IcKtTixd8mHxGibXJEizMmBV3/SyuvPLX6RtExdh9maOD/mzMPLI/QZ17OcRmabl2Yf4SLNZ1KpL5OLn20iTBzLzZUeclCstz6SfVZ6suhQihLaU9iaewTVJYw9tH9coUmZKhA1YOH8gsoHPn3xbDWmGk0g9MDhTFitT+Aoo4f6Ex+Bv5TXlImlStThNLlxnEIdNyGnp2gTjXxOWl7wz4PvpJHXWNLjWL1stLE+WHmn79U3+CLc84CdPfv5MKTEYgEyeOMqL35QOsEnRrsHvnwf8nEsrP/kkXpMnR8pidVekbXO0a3Fhmbif+Exevozx+3XNteH9SZB+M9q7UasurdxGnnqiXH59h3EC/cC/hvL78oFadSlECH0jbEfWxsyFbdGw/tooUmZEXdD4GFBtMhLFQYbhACBqw6TNRA0Mov5AaBOxOa7z2yht1h9UuaZWGzaFW21diHTqUUpMSW8GUmZE3dAQecrTIC/aCYOfWVdoe/5gaNaoJBhoUV58GIDNPI7SgrXEB+WIa3wLiBAiGfqX9ada0FfpW81CyoxoCCaIpGUHIVoBCodZXcyZlQbFxsKSnhBtecVfZmLgZUAlHKuNPxBzHcpNmnlcCHE7PFzU6jP0s/DBoVFG8KfSkYUQQgghSolezRZCCCFEqZEyI4QQQohSI2VGCCGEEKVGyowQQgghSo2UGSGEEEKUGikzQgghhCg1UmaEEEIIUWqkzAghhBCi1EiZEUIIIUSpkTIjhBBCiFIjZUYIIYQQpUbKjBBCCCFKjZQZIYQQQpQaKTNCCCGEKDVSZoQQQghRaqTMCCGEEKLUSJkRQgghRKmRMiOEEEKIUiNlRgghhBClRsqMEEIIIUqNlBkhhBBClBopM0IIIYQoNVJmhBBiGHHu3Llo0qRJ0YgRI6L169fHobf45ZdfomeffdadGzt2bPT555/HZ4QoN1JmhBBiGLFs2bLo448/jgYGBqL33nvPKTfG1atX3XFwcDDavHmzu06I4YCUGSFKDBMVT9tl4dixY12b37LIspYMz5w5E82aNSsaN25c9NBDD0X/+9//4jORC8cac+3atejy5cvRK6+8Ep9pDWVrn6IY3dSfpczUyfvvv+9MuVu2bIlDRD0gP0zeuIULF8ahIg9MSjxdjx49Og7pfiZPnhytXLlyiLUAqHtrB36f4veMGTOib775Jg5pDWWSZZoM88IEdPfdd0eHDh2KQ1pDGdunKEajbbGpDIrC7Ny5c3DMmDGDvb29cYiol56ensEFCxbEPpGXvr6+wenTpw8ODAzEIe2B9Gj/tH2OOKAOGU6sT4R+H+KYOHHiYH9/fxzyO9xHm/A5cOCA628cW0EzZJkmF8qCHKxMWXIpQpYMCSd+/xrSxxG+ePFid93u3bvd+VbQqfbZDaxbt87VMbJFDj7IY8WKFbedp14IM+e3D9o919o9SXXeSbLaYjuRMlMHVBwDgWgcBlgpM8Vhomh0QqwHJmkG6xAGXJskIfSH0H+S6j1JmQHCKHMraIYs0+RiA71N6rXkUoQ0GZIGyh+Tn69U4cLJtFXjWKfaZ6dB9lYnyDusa8KsnXAOPyCrJHlZXdq5VvWBRklri+1EykxB0D4ZCLpNOy4rdOZOd4Ky0cmBg3bP4GqDsMEE6U+MoT8JJtPwyZVyhXEDgzn9rtk0S5ZpciF+ZGHkkUsRkmTYaZol0zJC/VsboM1myQGlxuYRrjUFlHtM+UV5QaEpA51ui9ozUxB7G2DChAnuKIYPbIosw74d3kBZsmRJ7GsvtPuLFy9GX3zxxZDXftkfMW/evNh3uz+JypOp23uWBzauQrM3GzZLlmly2b9/f/Tkk0/GvnxyKUIRGbaLTrbPssC+paVLl1bnEdo37acyJ0d33XVXtG/fPhfOZu5Fixa5391Op9ti25QZOrht8LMByfxspC07NhFaebI2LDKg2XcgOLZi85TJm29JZG1S5pyfl1ZvtMxL3vzXknuWrBlQCDfHpsgHH3wwPvs7XGff5ggnUz9+HNclYWmF2L3hplfCkhQryvvtt9+6t1TS4P4smTUKb8ls27YtOnv2rPMjc/JDuPnnz5/v/MgrTSaPP/64k1+nyCPLIoRyIf49e/ZUJ6O8cilCp2UYUkumnKfc9IfhyPjx46PDhw+73+fPn68qK/RlKzMTPhtn2Rhtii/nr1y54n6H0FaQW9dstE2h423xloGm9WA2IzlbwwXWDHFlW7KhHKE5DfMxZQHKiMkwqVyYIDEdIg8c92Cey0Pe6sLMS5zETz65L8lUSV7MpGl5Id/txPLgkzf/kCX3WrLOs3RheUlaGjDTMEfzW15CyGdSu7E19nAtPOla4Pq09kIZiYt7zdTdTJCByYF8WP5ID2dyML+5rLxw3i8n+U+6Pk9dFSVLlkVIk0tYFyYPc82qI+JKaiudIEum5JFz5NfaynDE6p0+bWMRYZSZdmL1j7M2YJuGcYwh9GWwMY3wcIzoRshnp9pi25QZKpWCWiVRib5iUyb8CQwoU1iJXJOnfHR+7s1D3uto9H7adI60SdaHOmnG4F4E0qSj++TNf1G5h7KuNUHSZokvSSmFpLwnYZMdaYWDOOGWDz+dtIErK00UJuLgPNe1AuInv9aPG4X2ZooApOWdNLPaJu0jaRBFJqHMjbz1l4dmy6UIoQybAeVANmF50sKNNJkiHyZs7kvqB2J4kNUW6+mjRWjbMtOpU6eiSiOPbty44cxl/F6+fHl8tlyEJtQLFy644/333++OwDWUtRY3b96MKg0g9jUH1lmnTJkS+6LogQceiK5fvx77bgcTJuZBzKP//e9/49BbmDmcZQ+WfMxUSjhLJHwDBDB/cp547DcmVK6zZaC8S1h5819U7kVl/c4770Rr165N3R+FSZl2Xctkzh6CtCWFo0ePRvfcc49bb/7uu+9cGPGZXEO+//771Py8++67qeeaBfHbB9maAW2Ij7fV4sSJE+7aNP7+979Hc+fOHWKGZ5z59ddfq/ttQrJkWZRmy6UIeWVYBMrx9NNPu77L+AAc8ROeVs40mRK2adOmjshHtI+stlhPHy1C25SZH374wa2lMpmQeX8Caje2JyHLZfHzzz9HI0eOjH3p2NppFhs2bIj+85//xL7WgeyTYOJkr8hzzz3nBpx77703PnNr8Jo+fbprgBXF103ITCrABjU+uoXiQRx8HIsNbFyzd+/e6KuvvnLnP/roo+jTTz+NKk9lDZUzLf9JpMk9lDV1OGbMGNcBqXN/4AYUMxQjO8/R74iUtfJUEc2ePTt1jwrr47YZElmGnD592sn9qaeeir788ksXhoznzJnjfoeg1N13332xb3iQR+lHccyC8QRlzgZL5I5sd+zYEV9xO8NJlnlkWBQeNlGqrV+8/PLLzp/1EDoc26coRlpbrKePFqFtygwT2+7du51jU5w/abSb1atXu8k5yyVBnrE28BSdRxljV3oWaKXElbZb3SwavoLl+zmfFybtJNCIKW9/f79rWAxYxsaNG91kzeCFVYU6nDp1qjuHDJmEsbC9+uqrrkHaUxdPYGxm++2336K//e1vLrzRr4Cm5T+JJLknyZo6RLFGOR0YGKgO2Ab5nzZtWvU8ysiaNWvis7egrAcOHHDhoeWJ+I4cOVK1yiAHf5Mf581ShEUJhY0wnm7ZINgqUD79dpTkalmbmoltmM2CN4BqWcGQsw2WH374YXTw4ME7xhKQR4b1gDxRYP7whz84ayv+VhOOe6ErMu6J9pPVFlvZR9uizDBA8wTPgE1hGMDt1bPhDMsjadiubybDNKjoUMHy/ZzPS7g0FoJiwpsYKJoGystbb73lBhAsGrhQ8XrsscfcYBc2SCxxL730UjUcbb2WcpdFrfz7hHLPI2vy+ec//3lI+cFXRP76178mWoiQCZanXbt2xSG3wCpFfDYIc69vgj158mQ1r8gfZYmlJuTeSsulKbBZLjT7WhkacY2A/D/55BNnBaulaDHpNkJS3rvNtZtGZVqEcNwLXdK4lyQjuda7emlFe2qLMsOgzZM1Aza8/vrrTivzCbVxTPoQPkXyhG3YPTa4YbqyMISVtkej3mUmBlQmRCYof7mBPQ/gD7I8RaaZxrnu3//+d/Tmm2/GIc0FZdGWg6DI3gDfAoK1pre3tzqAhIoM+2FY0sGiE8KE/Mgjj8S+yMmMpZQ85M1/HrkXlXVoAaKMPnktRNzHN0ew6NggjBXL5/jx40NkhDLFd0mKKG7twsrQiEsj6XX4EOT5wgsvRH19fanr6yitjA+MN6tWrXJjAQ9SRUnKe7e5kDwyrAfkieyxTpq/G0mSkVzrXRJZbbFZfTSJtigzPLGaOR2eeeYZZ6nxJwr+e2tPT09VSKbMMHDxxMrEysTAJGmTF/ewdMA1xIXpis1p3M/yyKhRo9x1IfUuMxnkzf9PtEy0TFTbt293frME2MexqDAcUHEsyzC5m9WCp39fFo2CReSDDz5waaF0IbPnn3/enTNFzn6TNtfhyL9ZIoCJ25ZFkLk/kHE91ghMhlaX9s0EIIwNt1zHfcTlx51F3vzXkjv3Z8maOrG9LoShfNBuDNoWm4CJB+fLx98/g2zee++96MUXX3R+4L5ly5ZV0wXisPVkfpMvZGRQFsKwdqWBdevSpUuxLxniLhN5liCxWKHkpVmsqAvaGct6yJx6sv0eaeSRZVloxT9ztD0N7HkDHgjw+/08pJZMy9Y2RXHS2mI9fbQQlYm7pfDKFcng7JVLjvgrSkr1dVRez7Lr/Fe7OM91QFwVpaj6WiD32LW89odrB5WJ+bZXyXjlrDIZufyTR/8VNMJxwGvDVk7f5SHvdcjH5ElefXma7AHZ+nLnt8kWuI/7OVeZ2KtlogyEm5+y+X5kQ52ZPDj68fokvcqZN/+QJfdasuY89xBGOmH5+U257T47j7M0ccTh55Hfdr1h+eZa5G7+sOyc5xXyNJLkZdD+rTwWty+PbiSUHXmmjCHUVVq5wW+fPtRB2FeNLFmWiVCGzcDauN8fIC3cyJIp4TaecBwOshdDyWqL9fTRIuSbHdsAjTsJf0IKB3rusU5Fx2iGQPJAXrp9kugkDGh5FcvhMqG0C9p/Wl8pI2Ffoi0kKTP07WaXe7jIspvGo+HWPuuBuqAN24P6cMEewrLaWifbYtveZsrCPvMNmCF9sxN7ayodBKXLrdvavg3uwezsm/FtSQTzKK4V2BJXKzdolh32uLTC7C2iaObMme4tq2YuSxaBfhUuF7K8xx41nC31AcsRLAmGr7Qb9KXKxJerL/EaPeVuJlmyJL+MQ+Tf339nYTYO+DB2YTb3ZWAkmdK5HnkSlz9mhWnwG5dEERm2g063T5aNqS/kFcocWdmycqugTvlMhb3tybKK5SdcnuNaq1sc1xEGVg5cUntKIuse4rV25fdH7vHz4BPKkvJQLspn+fTpdFvsCmXmiSeecHsFEBrfPLG3XhAmey+Svk/CPbwdYh2eb4jgiONf//pX9PDDD7tw0X7Y44JCI5oPyvuCBQvcRud2YwOk/wFDBjUGbL4xxLeF+A1M/ryiyUMIe6AYAEP+7//+7zbFKA0bIK2/N4M0WTIRp+2/8/fphfApA8j7ATve6ORa4uLzB/b9ljAN9qbt3LnT/Q4pIsN20Mn2ibLAvjfmjIGBgTg0G2u/zBtM8tQ9Ez2TOBM45znHXIQfZRV52wQfsnXr1iF75xgL6Rt9fX3uzVAf9s/1xi9Y9Pf3uz6C/GjjvFjBfYTx0c1a1LqHtsnGXNLiBRy/PxJmzsiSJeVLehu5422xUgBRAExtiC1tzVgUQ8tMxWFN2vaRtRuWe/z6Cv02pFCvtlwUXmOwbByapLnO7vPhulYMV0myzNp/x/p+0p4Awsg7yyxJ9yaVn7GEJZmwvGEayCltySJJhp2mE+2T8bjWuEw7DGXNNgb2cgDnrO64lnrA77dj0qAuSCfc9gBJ9Qykk9SuDdJJyjuyDNOoRdI9lJE8gLVVIE+0QcplZa8lS8qfVM5Ot8WusMyUCUxtlUqrfnpeiHZjTz/NtFK0G6w8lQkvt0maf11Bv2s2SbLEosR/AE6CfNvbcoa9pcGbfVhw7KNhvgkfK7L9trQYS3iS5gnYX4Lw0yBuTPdcG1JUhu2iE+2T9lGZYIdsOzBseYVvFPFxS6sHwDpiqwKcs7ozqxh+llZ8qAvSwZLz448/xqHpEDerBWlQxyzLh3knnI9n2taKcEnInC/n8B6DfyWA43qsNrb6sXTpUvfh0Iri4uTA/VmyhG5ti1Jm6gCzOYNXmScTUW7YS8Zr50lr163EXrc3+NYPSwqEMRCawsG3fuzfSpw/f37IAMi1LAmHZvc06Gf0N8zjrSBJlkn772xvH4M8yxFM2tzDK/jsE6HsTJgWj/8JCCYH++0vUREXH6u0STTcC7h3714XP3COCQ2KyrDddKJ9khaOuvGXgOwjfL29vUM+/2GwjIi/v7/fLfEBbc4+B0Hb88th6bCH076IngYKD4oCy00oS+SNOrR6BJZ8UCp8bI8WygHpQ9onRaw9Jd2DHCgLYSgtlJG88wFQYI+TXzbDyhjKEgjz4bquaIsVYYg6wDyHWS3LdChqg/xohrgk06VIx8zFaebgZuPXFQ6zux+OudrCgLwRzpKDLZOQV8KTlpcsXr9P8bsd/cyXJWUgTfLC0fJq+fPzafm2/Jmf+HwI9+G8XcOSQJgGecCZ3w9Pk2G30e72ybKQySlpqQ/Zhe3IZMk9JlOus7Zsb9Nyn9VHuCzjQ5i/xMN9XI+z6wmzfCTlib5C/ObCtpNE2j0cSQOQD33Rzx9+rvfzB1my5H5bsuqmtjiCP5UMCyFKCNYQ3vRJMv12IzwlYs3pxvy2UpY83Yaf4ecJ9+rVq25ZIM2kn0Q3yzCkm9oncuPL4uGyUV64H8tb1pRJnVYmf7fsMlyhLfMhRdpsN7VFKTNCCCFEDZjE2fvU09OTqRChwHEdX7ovg8KZFxQ13mJasGBB1+3TAikzQgghhCg12gAshBBCiFIjZUYIIYQQpUbKjBBCCCFKjZQZIYQQQpQaKTNCCCGEKDVSZoQQQghRaqTMCCGEEKLUSJkRQgghRKmRMiOEEEKIUiNlRgghhBClRsqMEEIIIUqNlBkhhBBClBopM0IIIYQoNVJmhBBCCFFqpMyIruDzzz+Pjh07FvuEuDOg3eOEEI0hZUZ0nG+++Sb64YcfolmzZsUhQtwZPPvss9H+/ftdHxBC1M+IwQrxbyHazrVr16KFCxdGBw8ejMaNGxeHinZjddAtmLWCyX64Qx+YOXNmdPLkSfUBIepEyozoKCtXrozmzJlzR0xa3cqIESOibhwGaBv33XdftHr16jhk+LJlyxZ3vBPKKkQr0DKT6Bi//PKLewKfN29eHCLysn79eqeETJo0KTp37lwcegs7Z84myiRQGJIUmffff7+Qgol1gXvY98QR50Mex44d6+rcsPwb+LEQGe+++2506dKlIfcMV5YuXRq9/fbbTo5CiOJImREdY9++fW7ClGm9GOyvOHv2rFNCXnrppWjz5s3xmVtMmzbNncP19PS4ibIItsRz/fp1d8wDdXn58mW372n58uXOGUzQa9asccrM1atX49AoWrduXTR9+vTqfpHFixc7Bcbn+eefj06dOhX7hi/0gfnz5zs5CiGKI2VGdIwPPvggevLJJ2OfyMuPP/4YPfbYY+73I488cpvSYRYVs2gUVRa5f8qUKbEvH1jX3nvvvUQr0MaNG11+UVZu3rwZh97ixRdfjHbt2hX7omjChAnxr1vcc889Tkm6E3jqqaeiDz/8MPYJIYogZUZ0BCba/v5+t/FRtIZ33nmnplUm7xIOy0f+0pU5U15QQi5evBh98cUXbpnLwMpz6NAhtxdk9OjRThEDlp3Gjx8fLVq0KDpz5oyzzkycONGd86EMKGx3ApMnT3ay0FKTKDMsHSeNK4wX9vkNxgWWuJuJlBnREX766adozJgxmVYDJsqs/R53KigBhw8fdr/Pnz9ftWYwWBgMGmyerWWVYVnHv8/AguJPqiwf2dKV7/wNq6S1bds2twQGKCwMWCitpMFS040bN9y5vXv3Vq0/LJW98sor0aOPPur8hg12d8or+/fff787XrhwwR2FKBMoMPTzjz/++DYLa6i4YP1lfPL3zDVMZUASou309PQMLliwIPYNZWBgwJ2jeXKduB2Tz/Tp0wcryoIL87vz4sWLnRzz4tcFMicuc729vfGZZHbv3u0cHDhwYLCvr8/9Du83/549e9zR0iSf+H2I706se2SiNi/KyMSJExPHCsJsTAnPE75ixYrY1xh6NVt0BDT1X3/9tbrZ1Idzr732mjuy10Kvq3Y/PJWxuZelkqJ7dMTv8DaX2rwoG1jQsRYnfavKxnOWkSvKzG2WVqwzWHMatcBqmUl0BCa/Bx54IPYNhaWP0Ewpuhvqi8FIikxj3HXXXdUlRCHKAi9zLFmyJPb9jikyWfAW32effRb76kfKjBBCdAlpCr4Q3QoPpuyLu/fee+OQW7Bvjw+i1now5RpeEmgUKTOitGDaZMNZlhNCCNE67NtRfEbBB2tL3g9vogw1ipQZUVrYV8CWrywXkqTwyMl1wgkxXMEqw3enrK3bZxdmz5495CvfzUTKjLijSFJ45OQ64YQYDphFxv+6d/gpB7O8sAE4aZNw0jemiiJlRnQt+niYuNP4/vvv419ClAP2xKCMXLlyJQ4pxtGjR90m4EaRMiM6Ah0gbeDmC7K8rsfXUPnQGmbJ8J8pCjEc4V9T2L+qEKIs8OFLXq+uBzb/8j/YGkXKjOgIfP0x7R8Zbtq0Kfr555+rJkrMkvZ11GaBcsQG4qTPbrcb1pfJSyctUa2UB+UibvuUeachL/bPLdsJsiVtKeZiuMH+RcbstD7OwytjefgtGb4zxoNro9+YcVQSEKLt8KXYMWPGxL72whdn+UKuz86dO28Ly4I4uIcvWnLE+fAVXMpnX+cFuhtfyTTw+1/BTUrfvp7ZSsK0qRvySf7WrVsXh6aHG4RxjnLbF4Eh/BoxMqHc4ddAsyA+rsf5XysG4uZLyH6aXENeTP7mN5Bp3vTtXlyYNqSdTytnKA8f4igiFyG6Cdqv3w+zoA/642GjyDIjOsIf//jH6LfffuuINWLr1q3uvzUb9hXiNEtREvv27XP/zZkniuXLlztnUCaWx8aOHTtkU1xlso8qk27VKlCZ1NwHAoGPzbHunOepnfj5GBVvCdg/deM+0mNJjvOcwwqAf8aMGe4VScKS3iQI5fHCCy84029FIYveeuutODQ9HJAh/5OpMklHR44ccfkD+weS9jE98nrq1KnC9U58fF0YeYcbCPmv3OD/d23+1xPy/eijj5z/6aefjiqDrPsN/APO7du3x75sWPapjJVVK2FI0vmsciJr2k+I1T3lFKKM0AdYbqpl4WW8uHTpkrPmNAspM6Ij2KaxkydPxiHtg0mX/9ZsMNHbPz3My7x589yrhygMIUyuTHBMpvzDRh8msl27dsW+W3Iwpk2bFn377bexLx0mQv4VBAMHa9VM2CzDffXVVy4+lJPKU48z/W7YsMHtPUL5QNFg8AiXWEJ5EDfxjBw50ilgRlo48A8k+XotSgvnUVSB/5L9+OOPu9/A/ci76JeCV6xYkbh3ikGRgfONN96o/hNLOH78eLRjxw5XR8B91JlB+sglD+ztQhFEWTTF1yfpfFY5UeS//PLL2Pc7/INJlN2ishGim0Chp/1nQd+wB7lmIWVGdAwm4q+//jr2dRes/TJBhc6UFzrrxYsXoy+++MJtWDaYzLBeoEiMHj3aTebAZMp/u0ZpYBI1i4WP/wVNJm7S47sMWHksfWDS3rNnj/Nzzv5Lta074w//tw/5ZZLEkmN5ygIFAasS+5d80sKfeeaZ6n/HxorCf0QvgpU3dAbprVq1Kpo7d271qQ+ZkhaD4qhRo6pyAJQpysugaQpPqCQgC6BOk9K29X+UIhRH3+LkU+t8SNpAv3///mjZsmWxTwhRBCkzomMwATLRFF1yaAVYUPx8oBgwQYXOVxKYHLdt21adRG1ytUkdRcOsBXv37q1af1DiWAZ59NFHnd/wX23k6Yb0ent7nZXF0jewVOAnLeICJl97K4B8+OXhNw7LzNSpU+PQZCjHTz/95GTARO8rD2E4RxyyOH36tLP+MFmvXbvW3ZMGefGtVlbe0PmgmDz00ENOmeJ+Jn6UFpRClD4rL4oiG8wBeWCdevDBB53fx0zc1GlS2pQTmb788svuOh+UL86lnTfCcoLJ04frUILpE0KIOqh0WiE6Bv/+Pe+GsWbBRlU2sxpsRKMrmKu1AZP8Wp6Jh82+EN5v/j179rijbQ5l8yf+EPJlcRnERf58uN82nXLkHq5j4y1HNiNzjvvw85tzHEkjJJRHRTlw15ozksJJw/KHTLjGzy/xhmn6cfgbZtPw5WYbra38lpb5yYPFbbA52OrEIC424tYC2Vq5OVq9kx5xpp0HywfOLycyCTeM+3IUQhRHyozoKEwqTDYc2wVvmeSZyLIgDiazZuU7LU+k0cgkZ8pMFs2QRxZM5I3KyRSHZoGCRZydIJQHv1GE2tkHhBhujOBPZbATomOwLMCGzXAfRithyYTNtpj1a21WazUsVZw4ccK9YRPu62gUlkMoZ0Uhum0fjU8r5cESCpuUH3nkkeq+nk7CshhLbf6m53bA8hKbtyvKzJDvJrF8xsbwdudHiOGElBnRFbB3hg2w3TDZCdEu7O0nFBohRP1ImRFCCCFEqdHbTEIIIYQoNVJmhBBCCFFqpMwIIYQQotRImRFCCCFEqZEyI4QQQohSI2VGCCGEEKVGyowQQgghSo2UGSGEEEKUGikzQgghhCgxUfT/Ab21Nm/wKdovAAAAAElFTkSuQmCC\" width=\"563\" height=\"131\"\u003e\u003c/p\u003e \u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\)\u003c/span\u003e\u003c/span\u003e is soil erodibility factor, the unit is t\u0026middot;h/(MJ\u0026middot;mm);\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\:\\text{S}\\text{A}\\text{N}\\:\\)\u003c/span\u003e\u003c/span\u003eis sand content (%); \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{S}\\text{I}\\text{L}\\)\u003c/span\u003e\u003c/span\u003e is silt content (%); \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{C}\\text{L}\\text{A}\\)\u003c/span\u003e\u003c/span\u003e is clay content (%); \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{C}\\)\u003c/span\u003e\u003c/span\u003e is organic carbon content (%).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section3\"\u003e \u003ch2\u003e2.2.3. Topographic factors (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:LS\\)\u003c/span\u003e\u003c/span\u003e)\u003c/h2\u003e \u003cp\u003eTopographic factor is also one of the main causes of soil erosion. Slope length factor (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{L}\\)\u003c/span\u003e\u003c/span\u003e) and slope gradient factor (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{S}\\)\u003c/span\u003e\u003c/span\u003e) are selected in RUSLE model. \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{L}\\text{S}\\)\u003c/span\u003e\u003c/span\u003e reflects the influence of topography fluctuation on soil erosion. The greater the \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{L}\\text{S}\\)\u003c/span\u003e\u003c/span\u003e, the faster the soil erosion rate Although it is difficult to measure \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{L}\\text{S}\\)\u003c/span\u003e\u003c/span\u003e in large-scale research, the influence of topographic factors on soil erosion can still be expressed by topographic fluctuation rate, that is, the maximum height change in a certain height area of the ground (Wischmeier \u0026amp; Smith, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e1978\u003c/span\u003e). The specific formula is as follows:\u003cdiv id=\"Equ4\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ4\" name=\"EquationSource\"\u003e\n$$\\:H={H}_{max}-{H}_{min}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e5\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:H\\)\u003c/span\u003e\u003c/span\u003e topographic relief, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{H}_{max}\\)\u003c/span\u003e\u003c/span\u003e is the highest value of elevation, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{H}_{min}\\)\u003c/span\u003e\u003c/span\u003e is the lowest value of elevation.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section3\"\u003e \u003ch2\u003e2.2.4. Vegetation cover factor(\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:C\\)\u003c/span\u003e\u003c/span\u003e)\u003c/h2\u003e \u003cp\u003eVegetation cover factor plays an important role in soil and water conservation, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{C}\\:\\)\u003c/span\u003e\u003c/span\u003eranges from Between 0 and 1(Yang \u0026amp; Huang, \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). When \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{C}\\)\u003c/span\u003e\u003c/span\u003e=1, the ground is completely bare without vegetation; when \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{C}\\)\u003c/span\u003e\u003c/span\u003e=0 the ground vegetation coverage is good. It is a common way to obtain vegetation coverage from remote sensing images and calculate vegetation coverage factors according to vegetation coverage. Cai et al.(2000)The proposed algorithm is used to estimate the vegetation cover factor through artificial and natural rainfall experiments. The specific formula is as follows:\u003c/p\u003e \u003cp\u003e \u003cem\u003efvc\u003c/em\u003e \u003cspan class=\"InlineEquation\"\u003e \u003cspan class=\"mathinline\"\u003e\\(\\:=\\frac{NDVI-{NDVI}_{min}}{{NDVI}_{max}-{NDVI}_{min}}\\)\u003c/span\u003e \u003c/span\u003e (6)\u003cdiv id=\"Equ5\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ5\" name=\"EquationSource\"\u003e\n$$\\:C=\\left\\{\\begin{array}{c}1,\\:fvc\\le\\:0.095\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\\\\\:0.6508-0.3436Inc,0.095\\:\\:\u0026lt;\\:fvc\u0026lt;0.783\\:\\:\\\\\\:0,\\:fvc\\ge\\:0.783\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\:\\end{array}\\right.\\:\\:$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e7\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eWhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{N}\\text{D}\\text{V}\\text{I}\\)\u003c/span\u003e\u003c/span\u003e is normalized vegetation index(Wei et al., \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2022\u003c/span\u003e); \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{N}\\text{D}\\text{V}\\text{I}}_{\\text{m}\\text{i}\\text{n}}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{N}\\text{D}\\text{V}\\text{I}}_{\\text{m}\\text{a}\\text{x}}\\)\u003c/span\u003e\u003c/span\u003e are the minimum and maximum values of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{N}\\text{D}\\text{V}\\text{I}\\)\u003c/span\u003e\u003c/span\u003e respectively; \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{f}\\text{v}\\text{c}\\)\u003c/span\u003e\u003c/span\u003e is vegetation cover index (%);\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\:\\text{C}\\)\u003c/span\u003e\u003c/span\u003e is vegetation cover factor.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section3\"\u003e \u003ch2\u003e2.2.5. Soil and water conservation measures(\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:P\\)\u003c/span\u003e\u003c/span\u003e)\u003c/h2\u003e \u003cp\u003eSoil and water conservation measures can regulate the amount of soil erosion in a region. The factor of soil and water conservation measures is difficult to calculate directly. Usually, the value of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{P}\\)\u003c/span\u003e\u003c/span\u003e factor of a region is determined according to different land use types. The range of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{P}\\)\u003c/span\u003e\u003c/span\u003e is between 0 and 1. The larger \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{P}\\)\u003c/span\u003e\u003c/span\u003e is, the smaller the amount of soil erosion is. When \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{P}\\)\u003c/span\u003e\u003c/span\u003e=0, it means that no erosion occurs. When \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{P}\\)\u003c/span\u003e\u003c/span\u003e=1, it means that no soil and water conservation measures are taken (Humphrey et al., 2021). According to the research conclusions of relevant scholars (Peng et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). The following values were assigned to the study areas with different land use properties (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Using ArcGIS software to reclassify the assigned data, you can obtain \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{P}\\)\u003c/span\u003e\u003c/span\u003e values and their distribution maps for different land use types within the basin.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eFactor values of soil and water conservation measures in the study area.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"7\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eland use types\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003ecultivated land\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eforest land\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003egrass\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003ewaters\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003econstruction land\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eunused land\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003esoil and water conservation factor\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e0.25\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003e2.3. CMIP6 climate model data\u003c/h2\u003e \u003cp\u003eAccording to the applicability of CMIP6 in the study area (Lu et al., \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Wu et al., \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), six atmospheric circulation models with better performance in CMIP6 (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e) were selected. Due to the differences among the selected six climate models, the bilinear interpolation method was used to unify the spatial resolution of all models. After the interpolation, Delta downscaling (Chen et al., \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2011\u003c/span\u003e) was performed on the interpolated model data at each grid. Finally, the daily meteorological data for the Min-Tuo River Basin under three shared socioeconomic pathways (Table\u0026nbsp;\u003cspan refid=\"Tab4\" class=\"InternalRef\"\u003e4\u003c/span\u003e) were obtained for the future period.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eInformation on the six CMIP6 patterns used in this study.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"4\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\"\u0026times;\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNumber\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSchema name\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eCountries\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eSpatial resolution (\u0026deg;)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eACCESS-ESM1-5(ACC)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eAustralia\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\"\u0026times;\" colname=\"c4\"\u003e \u003cp\u003e1.2\u0026deg;\u0026times; 1.8\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eBCC-CSM2-MR(BCC)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eChina\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\"\u0026times;\" colname=\"c4\"\u003e \u003cp\u003e2.8\u0026deg;\u0026times; 2.8\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCanESM5(CAN)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eCanadian\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\"\u0026times;\" colname=\"c4\"\u003e \u003cp\u003e2.8\u0026deg;\u0026times; 2.8\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eIPSL-CM6A-LR(IPSL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eFrance\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\"\u0026times;\" colname=\"c4\"\u003e \u003cp\u003e1.3\u0026deg;\u0026times; 2.5\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMRI-ESM2-0(MRI)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eJapan\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\"\u0026times;\" colname=\"c4\"\u003e \u003cp\u003e1.1\u0026deg;\u0026times; 1.1\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMIROC6(MIR)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eJapan\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\"\u0026times;\" colname=\"c4\"\u003e \u003cp\u003e1.4\u0026deg;\u0026times; 1.4\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab4\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 4\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eIntroduction to SSPs\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNumber\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePath name\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eForced category\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eShared socio-economic pathways\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003e2100 year radiative forcing/(W/m\u0026sup2;)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSSP 1-2.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003elow forcing\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003esustainable development path\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e2.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSSP 2-4.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003elow forcing\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eintermediate path\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e4.5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSSP 5-8.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003ehigh forcing\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eTraditional fossil fuel-based pathways\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e8.5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cdiv id=\"Sec11\" class=\"Section3\"\u003e \u003ch2\u003e2.2.4. CMIP6 data processing\u003c/h2\u003e \u003cp\u003eThe Delta method (Chen et al, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2011\u003c/span\u003e) is used to downscale and bias correct the CMIP6 data of the selected 6 models. This method compares the differences between GCMs output and observational data, and downscales the climate projection from large scale to smaller spatial scales. It has the advantages of strong interpretability, high accuracy, and high computational efficiency. When describing the changes in climate elements, precipitation is represented using the rate of change, while temperature is represented using the amount of change. The calculation formulas for future precipitation and temperature scenarios at meteorological stations are as follows:\u003cdiv id=\"Equ6\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ6\" name=\"EquationSource\"\u003e\n$$\\:{P}_{f}={P}_{o}\\bullet\\:\\frac{{P}_{Gf}}{{P}_{Go}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e8\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Equ7\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ7\" name=\"EquationSource\"\u003e\n$$\\:{T}_{f}={T}_{o}+\\left({T}_{Gf}-{T}_{Go}\\right)$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e9\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eIn the equations, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{P}}_{\\text{f}}\\)\u003c/span\u003e\u003c/span\u003e、\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{T}}_{\\text{f}}\\)\u003c/span\u003e\u003c/span\u003erepresent the reconstructed future precipitation and temperature sequences; \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{P}}_{\\text{G}\\text{f}}\\)\u003c/span\u003e\u003c/span\u003e、\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{T}}_{\\text{G}\\text{f}}\\)\u003c/span\u003e\u003c/span\u003e represent the future precipitation and temperature sequences estimated by the climate model; \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{P}}_{\\text{G}\\text{o}}\\)\u003c/span\u003e\u003c/span\u003e、\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{T}}_{\\text{G}\\text{o}}\\)\u003c/span\u003e\u003c/span\u003e represent the multi-year average precipitation and temperature simulated by the climate model during the reference period; \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{P}}_{\\text{o}}\\)\u003c/span\u003e\u003c/span\u003e、\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{T}}_{\\text{o}}\\)\u003c/span\u003e\u003c/span\u003e represent the multi-year average precipitation and temperature of the observed field during the reference period.\u003c/p\u003e \u003cp\u003eThe standardized Taylor diagram (Taylor, 2001) is used to assess the model performance of different climate models. It effectively and intuitively illustrates the differences in performance among multiple models and the magnitude of the errors between simulated and actual values, and has been widely applied in climate model evaluation research. Single-model to multi-model ensemble prediction is one of the effective ways to improve model accuracy, and the calculation formula for using MME is as follows:\u003cdiv id=\"Equ8\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ8\" name=\"EquationSource\"\u003e\n$$\\:{F}_{MME}=\\frac{1}{n}\\sum\\:_{i=1}^{n}{F}_{i}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e10\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eIn the equation, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{F}}_{\\text{M}\\text{M}\\text{E}\\:}\\)\u003c/span\u003e\u003c/span\u003erepresents the average result of the multi-model ensemble, and\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\:{\\text{F}}_{\\text{i}\\:}\\:\\)\u003c/span\u003e\u003c/span\u003erepresents the simulation results of each individual model.\u003c/p\u003e \u003cp\u003eQuantitative evaluation was conducted on the future precipitation results of the Min-Tuo River Basin estimated by the six climate models and the multi-model ensemble average, which intuitively shows the relationship between the predicted precipitation and observed precipitation under the six models and the ensemble average(Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e) Therefore, overall, the future precipitation simulation values after the multi-model ensemble average are closer to the observed precipitation values than those of a single model, and the reliability of the simulation results is higher. Therefore, in this study, the multi-model ensemble average was used as the future precipitation data, and the future soil erosion was predicted based on this.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section3\"\u003e \u003ch2\u003e2.2.5. CA-Markov Model\u003c/h2\u003e \u003cp\u003eIn order to estimate the soil erosion modulus of the Min-Tuo River Basin under future climate conditions, it is necessary to predict the land use types in order to obtain the P factor values. Because the Cellular Automata (CA) model and Markov model complement each other well, the CA-Markov model is widely used for future land prediction (Abijith \u0026amp; Saravanan, \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Aqil et al., \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Wang \u0026amp; Wang, \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Therefore, the CA-Markov model is used to predict future land use in the study area.\u003c/p\u003e \u003cp\u003eThe \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\text{a}\\text{p}\\text{p}\\text{a}\\)\u003c/span\u003e\u003c/span\u003e coefficient can be used to test the consistency between the simulation results and the actual data, and is commonly used to assess the accuracy of predicting land use changes. The formula for calculating the \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\text{a}\\text{p}\\text{p}\\text{a}\\)\u003c/span\u003e\u003c/span\u003e coefficient is as follows:\u003cdiv id=\"Equ9\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ9\" name=\"EquationSource\"\u003e\n$$\\:Kappa=\\frac{{p}_{a}-{p}_{c}}{{p}_{p}-{p}_{c}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e11\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eThe formula where \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{p}}_{\\text{a}}\\)\u003c/span\u003e\u003c/span\u003erepresents the probability of correct simulation, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{p}}_{\\text{c}}\\)\u003c/span\u003e\u003c/span\u003e represents the probability of simulated prediction, and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\text{p}}_{\\text{p}}\\)\u003c/span\u003e\u003c/span\u003e represents the probability of ideal simulation. The \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\text{a}\\text{p}\\text{p}\\text{a}\\)\u003c/span\u003e\u003c/span\u003e coefficient ranges between [-1,1], and when 0.6\u0026le;\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\text{a}\\text{p}\\text{p}\\text{a}\\)\u003c/span\u003e\u003c/span\u003e\u0026le;0.8, it indicates significant consistency; when 0.8\u0026lt;\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\text{K}\\text{a}\\text{p}\\text{p}\\text{a}\\)\u003c/span\u003e\u003c/span\u003e\u0026le;1, it indicates very good consistency.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e"},{"header":"3. Results and analysis","content":"\u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Soil Erosion Historical Analysis\u003c/h2\u003e \u003cp\u003eBased on the calculated soil erosion impact factors, the soil erosion modulus for the years 2000, 2010, and 2020 was obtained, and the intensity of soil erosion was analyzed. According to the \"Classification and Grading Standards for Soil Erosion in China\" (SL190-2007) (2008), the soil erosion intensity in the study area was divided into six erosion levels: slight erosion, mild erosion, moderate erosion, severe erosion, very severe erosion, and extreme erosion (Table\u0026nbsp;\u003cspan refid=\"Tab5\" class=\"InternalRef\"\u003e5\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\u003cdiv class=\"gridtable\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e\u003ctable float=\"Yes\" id=\"Tab5\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 5\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eClassification of soil erosion modulus.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e\u003ccolgroup cols=\"7\"\u003e\u003c/colgroup\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003elevel\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003emicrometric\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003emild\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003emoderate\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003estrong\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eextremely strong\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eviolently\u003c/p\u003e \u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eErosion modulus [t/(km²·a)]\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u0026lt;500\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e500 ~ 2500\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e2500 ~ 5000\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e5000 ~ 8000\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e8000 ~ 15000\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u0026gt;15000\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/table\u003e\u003c/div\u003e \u003cp\u003e\u003c/p\u003e \u003cp\u003eAs seen from Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, the erosion in Ya'an, Liangshan Prefecture, Ganzi Prefecture, and the southern part of Aba Prefecture is showing an increasing trend.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eAccording to Table\u0026nbsp;\u003cspan refid=\"Tab6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, it is known that the soil erosion modulus reached its highest value in 2020, at 34939.8t/(km²·a). The overall trend over the 20 years from 2000 to 2020 first decreased and then increased. The increase from 2010 to 2020 was greater than the decrease from 2000 to 2010, at 7.1% and 2.6% respectively.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\u003cdiv class=\"gridtable\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003ctable float=\"Yes\" id=\"Tab6\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 6\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eStatistical table of soil erosion modulus in the study area from 2000 to 2020.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e\u003ccolgroup cols=\"4\"\u003e\u003c/colgroup\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003etime\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMinimum t/(km²·a)\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eMaximum value [t/(km²·a)]\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eRange of change (%)\u003c/p\u003e \u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e2000\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e33493.6\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e2010\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e32621.4\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e-2.6\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e2020\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e34939.8\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.1\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/table\u003e\u003c/div\u003e \u003cp\u003e\u003c/p\u003e \u003cp\u003eFrom Table\u0026nbsp;\u003cspan refid=\"Tab7\" class=\"InternalRef\"\u003e7\u003c/span\u003e, it can be observed that: (1) In 2000, slight erosion and moderate erosion were the main types, covering areas of 112688.41km² and 20007.15km², accounting for 69.72% and 12.38% of the total erosion area, respectively; (2) In 2010, slight erosion and moderate erosion were predominant, covering areas of 97840.53km² and 25076.86km², accounting for 60.54% and 15.52% of the total erosion area, respectively; (3) In 2020, slight erosion and moderate erosion were predominant, covering areas of 54167.59km² and 48477.77km², accounting for 33.52% and 29.99% of the total erosion area, respectively; (4) From 2000 to 2020, the area of slight erosion showed a decreasing trend, reaching its minimum in 2020 at 54167.59km², accounting for 33.52% of the total erosion area. The areas of mild, moderate, and severe erosion showed an increasing trend over the 20 years, reaching their maximum values in 2020 at 28669.73km², 48477.77km², and 24358.86km², accounting for 17.74%, 29.99%, and 15.07% of the total erosion area, respectively. The areas of very severe and extremely severe erosion showed an increasing trend over the 20 years, but with a smaller magnitude, reaching their maximum values in 2020 at 5862.79km² and 84.56km², accounting for 3.63% and 0.05% of the total erosion area, respectively.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\u003cdiv class=\"gridtable\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e\u003ctable float=\"Yes\" id=\"Tab7\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 7\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eStatistical Table of Soil Erosion Area in Study Area from 2000 to 2020.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e\u003ccolgroup cols=\"8\"\u003e\u003c/colgroup\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colspan=\"2\" nameend=\"c2\" namest=\"c1\"\u003e \u003cp\u003eYear\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colspan=\"2\" nameend=\"c4\" namest=\"c3\"\u003e \u003cp\u003e2000\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colspan=\"2\" nameend=\"c6\" namest=\"c5\"\u003e \u003cp\u003e2010\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colspan=\"2\" nameend=\"c8\" namest=\"c7\"\u003e \u003cp\u003e2020\u003c/p\u003e \u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colspan=\"2\" nameend=\"c2\" namest=\"c1\"\u003e \u003cp\u003eArea (proportion)\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eArea/km²\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eRatio/%\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eArea/km²\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003eRatio/%\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003eArea/km²\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e \u003cp\u003eRatio/%\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"5\" rowspan=\"6\"\u003e \u003cp\u003etype\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMicro erosion\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e112688.41\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e69.72\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e97840.53\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e60.54\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e54167.59\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e33.52\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMild erosion\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e17808.95\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e11.02\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e17830.72\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e11.03\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e28669.73\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e17.74\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eModerate erosion\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e20007.15\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e12.38\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e25076.86\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e15.52\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e48477.77\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e29.99\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eIntense erosion\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e8677.3\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e5.37\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e16063.14\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e9.94\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e24358.86\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e15.07\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eExtreme erosion\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e2403.79\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1.49\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e4743.96\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e2.94\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e5862.79\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e3.63\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSevere erosion\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e35.7\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.02\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e66.1\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.04\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e84.56\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e \u003cp\u003e0.05\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/table\u003e\u003c/div\u003e \u003cp\u003e\u003c/p\u003e \u003cp\u003eThe Sankey diagram (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e) provides a more intuitive view of the erosion grade transitions.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003e3.2. Soil Erosion Future Forecast and Analysis\u003c/h2\u003e \u003cdiv id=\"Sec16\" class=\"Section3\"\u003e \u003ch2\u003e3.2.1. Precipitation Forecast\u003c/h2\u003e \u003cp\u003eThe precipitation forecast based on the multi-model ensemble average is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn spatial terms, Wang et al. (\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) used China as the study area, and the results indicated that the spatial estimates of precipitation and extreme precipitation indices from CMIP6 were higher than the actual values. In this study, based on the MME simulations, the future high precipitation area in the Min-Tuo River basin is located in the southwest of the study area. This may be due to the enhanced disturbance of climate in the high-altitude areas of the Min-Tuo River basin under the selected future climate change scenarios, resulting in abnormally high precipitation in these high-altitude regions. The low precipitation area in the future is located in the northern part of the study area, which aligns with the actual spatial distribution of annual precipitation in the study area at present.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section3\"\u003e \u003ch2\u003e3.2.2. Land Use Forecast\u003c/h2\u003e \u003cp\u003eWith the help of IDRISI software, the Kappa index for simulated land use and actual land use in the study area in 2010 and 2020 was calculated to be 0.87 and 0.86, respectively, indicating good consistency. Considering that the forecast period of this study is from 2030 to 2050, we selected representative years 2030, 2040, and 2050 with a 10 year interval to visually display the evolution of future land use types in the Min-Tuo River basin (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e). Under future climate change scenarios, the area of cropland conversion increases. However, if agricultural activities are not managed properly, this could potentially increase the risk of soil erosion in the Min-Tuo River basin.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section3\"\u003e \u003ch2\u003e3.2.3. Soil Erosion Time Distribution Simulation\u003c/h2\u003e \u003cp\u003eUsing the MME method, the future soil erosion values for the Min-Tuo River basin were calculated for three scenarios. Considering the length of the time series, representative years were selected at 10-year intervals to calculate the average soil erosion modulus and total erosion for the Min-Tuo River basin (Table\u0026nbsp;\u003cspan refid=\"Tab8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). Over time, there are differences in soil erosion in the Min-Tuo River basin under different scenarios. Comparatively, in the SSP5-8.5 scenario, the multi-year average soil erosion modulus and total amount in the Minjiang River basin are ranked as follows: SSP5-8.5 \u0026gt; SSP1-2.6 \u0026gt; SSP2-4.5.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\u003cdiv class=\"gridtable\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003ctable float=\"Yes\" id=\"Tab8\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 8\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eSoil erosion statistics of Min-Tuo River in the future.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e\u003ccolgroup cols=\"4\"\u003e\u003c/colgroup\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSSPs\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eYear\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eAverage soil erosion modulus [t/(km²·a)]\u003c/p\u003e \u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eTotal erosion(10⁸t)\u003c/p\u003e \u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"2\" rowspan=\"3\"\u003e \u003cp\u003eSSP1-2.6\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2030\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e442.88\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.22\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2040\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e442.86\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.22\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2050\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e445.46\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.26\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"2\" rowspan=\"3\"\u003e \u003cp\u003eSSP2-4.5\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2030\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e437.39\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.13\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2040\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e437.99\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.14\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2050\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e440.64\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.18\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"2\" rowspan=\"3\"\u003e \u003cp\u003eSSP5-8.5\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2030\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e446.58\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.28\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2040\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e443.36\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.23\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2050\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e447.74\u003c/p\u003e \u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e7.3\u003c/p\u003e \u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/table\u003e\u003c/div\u003e \u003cp\u003e\u003c/p\u003e \u003cp\u003eIn the SSP5-8.5 scenario, the future changes in soil erosion in the Min-Tuo River basin exhibit a \"V\" shape, meaning that the total soil erosion decreases before 2030, but gradually increases after 2030. In terms of the year when the maximum erosion occurs, in all three scenarios, the maximum soil erosion occurs in 2050. This may be due to the fact that the maximum rainfall also occurs in the 2040s in all three scenarios. Over the future time scale, the increase in heavy rainfall or precipitation variability leads to soil erosion and substantial soil loss.\u003c/p\u003e \u003cp\u003eFigure 7 presents the average soil erosion intensity level in the Min-Tuo River basin during the forecast period under the three different scenarios. Therefore, it is necessary to establish an early warning system for soil erosion forecasting to address the future changes in soil erosion.\u003c/p\u003e\u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section3\"\u003e \u003ch2\u003e3.2.4. Prediction of spatial evolution of soil erosion\u003c/h2\u003e \u003cp\u003eFurthermore, the spatial distribution of soil erosion in Min-Tuo River in the future is discussed in three scenarios on the spatial scale (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). Under SSP1-2.6 scenario, the annual average soil erosion intensity grade and intensity change of Min-Tuo River in the future are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e(a). Under SSP1-2.6 scenario, the area of soil erosion in Min-Tuo River reaches 5.35×10⁴km² (above mild erosion), accounting for 32.81% of the total area of the basin. The areas with increased soil erosion intensity are scattered in Ya'an city, Liangsha, Ganzi and the south of Aba in the southwest of the basin. Therefore, under the low emission forcing scenario of SSP1-2.6 in the future, it is necessary to focus on preventing further increase of soil erosion intensity in the corresponding areas mentioned above. As far as the soil erosion grade change under SSP1-2.6 scenario is concerned, the soil erosion grade of 91.76% area of Min-Tuo River will remain stable in the future, 7.13% area will be aggravated, and 1.11% area will be alleviated.\u003c/p\u003e \u003cp\u003eUnder SSP2-4.5 scenario, the spatial characteristics of soil erosion intensity in Min-Tuo River are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e(b). SSP2-4.5 medium emission scenario, the soil erosion grade distribution tends to be consistent with SSP1-2.6 scenario, but in the future, the soil erosion intensity grade is different from SSP1-2.6 scenario, specifically reflected in: 91.89% of the Min-Tuo River soil erosion grade will be stable, 7.00% of the area soil erosion degree will be aggravated, erosion degree will be slightly reduced.\u003c/p\u003e \u003cp\u003eSimilarly, under SSP5-8.5 scenario, the spatial characteristics of soil erosion intensity of Min-Tuo River are similar to those of SSP1-2.6 and SSP2-4.5 scenarios, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e(c), except that the change of soil erosion grade is different from the other two models. The differences are reflected in: In future SSP5-8.5 high emission scenario, the area of soil erosion grade increase in Min-Tuo River will further increase on the basis of SSP2-4.5 scenario, reaching 12.50%, in which the area of soil erosion intensity increase sharply is the same as that in SSP2-4.5 scenario, and the soil erosion grade in 91.39% area of Min-Tuo River will remain stable in the future, and the soil erosion degree will increase by 0.52% compared with SSP2-4.5 scenario. Therefore, it is necessary to identify potential areas with abrupt increase in soil erosion intensity rate under future high emission scenarios, and to prevent slight soil erosion areas from transforming into severe soil erosion areas.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e "},{"header":"4. Discussion","content":"\u003cp\u003eThis study demonstrates the dynamic changes in soil erosion patterns in the Min-Tuo River basin over a 10-year period and forecasts future soil erosion using the recently developed CMIP6 model for approximately 30 years. The selected CMIP6 model is part of the World Climate Research Programme, and the approved climate change simulation and estimation data ensure the scientific credibility of the predictions.\u003c/p\u003e\u003ch2\u003e4.1. Soil erosion aggravation in Min-Tuo River in the future\u003c/h2\u003e\u003cp\u003eThe severity of soil erosion in the Min-Tuo River basin is expected to intensify in the future, leading to increasingly serious issues of soil and water loss. The relationship between future climate and soil erosion is closely intertwined, with climate factors directly or indirectly influencing soil erosion. Among these factors, precipitation has the most significant impact on soil erosion and forms its foundation. This study considers future precipitation conditions, which are projected to increase under three scenarios compared to historical periods, especially with intensified and prolonged precipitation events, including heavy rainfall, which will exacerbate soil and water loss (Yao et al., \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Additionally, as future temperatures are expected to rise continuously (Song \u0026amp; Yan, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), leading to high-temperature droughts and extreme heatwaves, soil cracking will occur, resulting in serious consequences such as debris flows when rainfall infiltrates. Addressing the increasingly severe soil erosion caused by future high-temperature heavy rainfall is a matter of significant concern and requires attention and resolution.\u003c/p\u003e\u003ch2\u003e4.2. Strengthening soil and water conservation measures in farmland\u003c/h2\u003e\u003cp\u003eUnder future climate scenarios, the increased magnitude of precipitation in the basin will exacerbate soil erosion, and changes in land use may further deteriorate soil erosion in certain areas. Considering the compounded effects of land use change, a significant conversion of grassland to cropland may potentially worsen soil erosion. Therefore, in the process of soil and water erosion control in the future Min-Tuo River basin, attention should be paid to maintaining forest and grassland coverage in the watershed and ensuring stable growth of vegetation cover to mitigate or alleviate the current soil and water loss situation. In the context of increased rainfall in the future, it is advisable to increase the proportion of grassland appropriately, leveraging the role of maintaining surface cover to mitigate soil and water loss while reducing ecological water consumption in the region to safeguard water resources.\u003c/p\u003e\u003cp\u003eFurthermore, to prevent the worsening of soil and water loss, the conversion of land from forests, grasslands, and wetlands to cropland should be avoided (Fang, \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). It is recommended to consolidate the achievements of returning farmland to forest and grassland projects, continue to maintain stable growth of vegetation cover to mitigate future soil erosion caused by surface cover factors, and strengthen comprehensive management of small watersheds with slope farmland treatment and accompanying measures for slope surface water systems to enhance watershed water conservation capacity in the upstream areas of the basin.\u003c/p\u003e\u003ch2\u003e4.3. Uncertainty of simulation results\u003c/h2\u003e\u003cp\u003eThe simulation results exhibit a certain degree of uncertainty, stemming from both the uncertainties in climate simulation and soil erosion simulation. Firstly, the uncertainty in climate simulation arises from the inability to access high-resolution climate simulation data, leading to reliance on downscaled climate simulation data for the study area. Enhancing the resolution of climate simulation in the future could improve the accuracy of simulations.\u003c/p\u003e\u003cp\u003eSecondly, the uncertainty in soil erosion simulation is intrinsic to the soil erosion models themselves, affecting the precision of future soil erosion simulations. In this study, the classic RUSLE model was employed to calculate the soil erosion modulus in the Min-Tuo River basin. Compared to other methods for assessing soil loss, the RUSLE model tends to overestimate erosion due to its omission of the sedimentation process in soil erosion calculations, focusing solely on soil erosion quantity. Therefore, it is necessary to consider how to adjust the model to better suit the soil erosion conditions in the Min-Tuo River basin.\u003c/p\u003e\u003cp\u003eIn estimating the soil erosion modulus, it is essential to select more reasonable and applicable methods for calculating the factors of the RUSLE model based on actual conditions in the study area. Specifically, the factors are determined based on land use type, but the allocation values are somewhat subjective due to the lack of consideration for different soil and water conservation measures for different land use categories. Implementing appropriate land management interventions to reverse the trends in land use/land cover change and soil erosion in the study watershed is crucial. Additionally, precision in the calibration and calculation of soil erosion factors needs further improvement. For instance, in calculating the rainfall erosivity factor, using daily or even sub-daily rainfall data when available can enhance accuracy. Moreover, for the calculation of soil erodibility values, field sampling to determine basic soil properties and subsequent verification can enhance precision. Furthermore, while this study only considered land use type as a human factor, factors such as population density and GDP could also be included to further investigate the impact of socio-economic development on soil erosion in the Min-Tuo River basin, enriching the research content.\u003c/p\u003e\u003cp\u003eRegarding the prediction of future soil erosion in the Min-Tuo River basin, this study only calculated the soil erosion modulus and severity at the annual temporal scale, without elucidating the seasonal or monthly variations in future soil erosion conditions. Thus, the seasonal and monthly variations in soil erosion severity under future scenarios remain unclear. Future research could explore the spatiotemporal distribution characteristics of soil erosion at seasonal and monthly scales.\u003c/p\u003e"},{"header":"5. Conclusions","content":"\u003cp\u003eBased on the prediction of future precipitation and land use types of Min-Tuo River, this study predicts the soil erosion modulus in the future period (2021\u0026thinsp;~\u0026thinsp;2050) by using the CMIP6 and scenario data published by the World Climate Research Program, and expounds its temporal and spatial evolution characteristics as follows:\u003c/p\u003e \u003cp\u003eThe analysis of soil erosion in the study area shows that the soil erosion types from 2000 to 2020 are mainly slight erosion, light erosion and moderate erosion, and the area of light erosion decreases especially from 69.72% in 2000 to 33.52% in 2020. The area of soil erosion reduction is less than that of soil erosion intensification in the study area from 2000 to 2010 and 2010 to 2020. It indicates that soil erosion will be aggravated gradually from 2000 to 2020;\u003c/p\u003e \u003cp\u003eSix CMIP6 models and three different emission scenarios were selected to predict the future precipitation of Min-Tuo River. The prediction results showed that the future precipitation of Min-Tuo River fluctuated and increased under three scenarios, and the fluctuation amplitude and increment of future precipitation of Min-Tuo River were the largest under SSP5-8.5 scenario;\u003c/p\u003e \u003cp\u003eAs far as the future land use types are concerned, the top three categories are: cultivated land to forest land, forest land to cultivated land and grassland to cultivated land, with the conversion areas of 5339km\u0026sup2;, 40874km\u0026sup2; and 27393km\u0026sup2; respectively. Under the future climate change scenario, the conversion area of cultivated land will increase, but if the agricultural management activities are improper, it is likely to increase the risk of soil erosion in Min-Tuo River;\u003c/p\u003e \u003cp\u003eCompared with the present stage, the area of soil erosion in Min-Tuo River will increase under the future three climate scenarios. Comparatively speaking, the annual average soil erosion modulus and total amount of Min-Tuo River will be the largest under SSP5-8.5 low emission scenario, with the variation between 446.58\u0026thinsp;~\u0026thinsp;447.74 t/(km\u0026sup2;\u0026middot;a) and 7.28\u0026thinsp;~\u0026thinsp;2.30\u0026times;10⁸t respectively, followed by SSP1-2.6 scenario, and the soil erosion degree of Min-Tuo River will be the lightest under SSP2-4.5 scenari.\u003c/p\u003e \u003cp\u003eSpatially, high-intensity soil erosion will occur in Ya'an, Liangshan, Ganzi and the south of Aba in the southwest of the basin under three different scenarios in the future The stability of soil erosion grade in SSP2-4.5 scenario is higher than that in SSP1-2.6 scenario and SSP5-8.5 scenario; Comparing the changes of soil erosion grade under three scenarios, we can see that the risk of soil erosion degree will become larger and larger with the emission grade of greenhouse gases in the future.\u003c/p\u003e "},{"header":"Declarations","content":"\u003cp\u003eCredit authorship contribution statement\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eNan Jiang:\u0026nbsp;\u003c/strong\u003eData curation, Formal analysis, Investigation, Methodology, Resources,\u0026nbsp;Software, Supervision, Validation, Visualization, Writing\u0026ndash;original draft. \u003cstrong\u003eFuquan Ni:\u0026nbsp;\u003c/strong\u003e Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing\u0026ndash;review \u0026amp; editing. \u003cstrong\u003eYu Deng:\u0026nbsp;\u003c/strong\u003eData curation, Formal analysis, Investigation, Methodology, Resources. \u003cstrong\u003eMingyan Wu:\u003c/strong\u003e Formal analysis, Investigation, Methodology, Resources, Software. \u003cstrong\u003eMengyu Zhu:\u0026nbsp;\u003c/strong\u003eFormal analysis, Investigation, Methodology, Resources, Software. \u003cstrong\u003eYuxuan Wang:\u0026nbsp;\u003c/strong\u003eData curation, Formal analysis, Investigation, Methodology, Software. \u003cstrong\u003eHuazhun Ren:\u0026nbsp;\u003c/strong\u003eFormal analysis, Investigation, Methodology, Resources, Software. \u003cstrong\u003eZiying Yue:\u0026nbsp;\u003c/strong\u003eData curation, Formal analysis, Investigation, Methodology, Software.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eDeclaration of competing interest\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.\u0026nbsp;\u003c/p\u003e\u003cp\u003eThe authors declare the following financial interests/personal relationships which may be considered as potential competing interests:\u003c/p\u003e\n\u003cp\u003eData availability\u003c/p\u003e\n\u003cp\u003eThe data used in this study is as follows: the parameters for the RULSE equation are sourced from Table 1, where the rainfall data is obtained from National Geographic Resource Science SubCenter, National Earth System Science Data Center, National Science \u0026amp; Technology Infrastructure of China, it can be found by the DOI:10.12041/geodata.269728669112188.ver1.db of the data used; Soil data is sourced from the World Soil Database, with the website address http://www.fao.org/soils-portal/so, after entering the URL, select Soil Maps and Databases in the Date Hub to download the relevant data; Terrain data is derived from the global 30 m resolution DEM data set, which is derived from the CGIAR-CSI SRTM Elevation Database(http://srtm.csi.cgiar.org/download); The NDVI utilizes NASA\u0026apos;s MODIS satellite data, visit NASA\u0026apos;s MODIS data download website (https://modis.gsfc.nasa.gov/data/dataprod/mod13.php) and select the NDVI data required; Land use / cover data can be accessed through http://www.globeland30.org, and data download requires registration and application submission, this data is prepared for estimating the P factor, and the value of the P factor is ultimately determined by empirical methods based on relevant literature of the study area(Peng et al. 2017). The data for CMIP6 is obtained from https://esgf-node.llnl.gov/search/cmip6/, the models and scenarios used are listed in Table 2 and Table 3, respectively.\u003c/p\u003e\n\u003cp\u003eFunding\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe authors declare that no funds, grants, or other support were received during the preparation of this manuscript.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eAbijith, D., \u0026amp; Saravanan, S. (2021). Assessment of land use and land cover change detection and prediction using remote sensing and CA Markov in the northern coastal districts of Tamil Nadu, India. \u003cem\u003eEnvironmental Science and Pollution Research, 29(57)\u003c/em\u003e, 86055-86067. https://doi.org/10.1007/s11356-021-15782-6\u003c/li\u003e\n\u003cli\u003eAqil, T., Jianguo, Y., \u0026amp; Faisal, M. (2022). Land change modeler and CA-Markov chain analysis for land use land cover change using satellite data of Peshawar, Pakistan. \u003cem\u003ePhysics and Chemistry of the Earth 128, \u003c/em\u003e103286. https://doi.org/10.1016/j.pce.2022.103286\u003c/li\u003e\n\u003cli\u003eBekele, T. (2019). Effect of Land Use and Land Cover Changes on Soil Erosion in Ethiopia. \u003cem\u003ePeertechz Publications Private Limited, 20(1),\u003c/em\u003e 79-91. https://doi.org/10.1080/15715124.2020.1767636\u003c/li\u003e\n\u003cli\u003eCai, C. F., Ding, S. W., Shi, Z. H., Huang, L., \u0026amp; Zhang, G. Y. (2000). Study of applying USLE and geographical information system IDRISI to predict soil erosion in small watershed. \u003cem\u003eJournal of Soil and Water Conservation, 14,\u003c/em\u003e 19\u0026ndash;24. doi: 10.13870/j.cnki.stbcxb.2000.02.005\u003c/li\u003e\n\u003cli\u003eChen, J., Brissette, F. O. P., \u0026amp; Leconte, R. (2011). Uncertainty of downscaling method in quantifying the impact of climate change on hydrology. Journal of Hydrology 401(3-4), 190-202, https://doi.org/10.1016/j.jhydrol.2011.02.020\u003c/li\u003e\n\u003cli\u003eDeng, B.,\u0026amp; Shi, L. (2022). Remote Sensing Quantitative Research on Soil Erosion in the Upper Reaches of the Minjiang River. \u003cem\u003eFrontiers in Earth Science 10\u003c/em\u003e\u003cem\u003e,\u003c/em\u003e 930535, https://doi.org/10.3389/feart.2022.930535\u003c/li\u003e\n\u003cli\u003eWoodward, D. E. (1999). Method to predict cropland ephemeral gully erosion - ScienceDirect. \u003cem\u003eCatena 37(3-4),\u003c/em\u003e 393-399. https://doi.org/10.1016/S0341-8162(99)00028-4 \u003c/li\u003e\n\u003cli\u003eFan, L., Lu, C., \u0026amp; Chen Z. (2012). A Review of EPIC Model and Its Applications. \u003cem\u003eProgress in Geography, 31(05), \u003c/em\u003e584-592. doi: 10.11820/dlkxjz.2012.05.006\u003c/li\u003e\n\u003cli\u003eSuresh, B. S. (2023). Adaptation in the Scientific Method: An Outline for Mutual Learning and Knowledge Co-production in Climate Science. Science, \u003cem\u003eTechnology and Society 28(4),\u003c/em\u003e 621-638. https://doi.org/10.1177/09717218231197448\u003c/li\u003e\n\u003cli\u003eGuo, B., Tao, H., Liu, B., \u0026amp; Jiang, L. (2012). Characteristics and analysis of soil erosion in Li country after Wenchuan earthquake based on GIS and USLE. \u003cem\u003eTransactions of the Chinese Society of Agricultural Engineering\u003c/em\u003e\u003cem\u003e, 28(14),\u003c/em\u003e 118-126. doi: 10.3969/j.issn.1002-6819.2012.14.01\u003c/li\u003e\n\u003cli\u003eHairsine, P. B.,\u0026amp; Rose, C. W. (1992). Modeling water erosion due to overland flow using physical principles: 2. Rill flow. \u003cem\u003eWater Resources Research 28(1),\u003c/em\u003e 245\u0026ndash;250. https://doi.org/10.1029/91WR02381\u003c/li\u003e\n\u003cli\u003eFang, H. (2021). Changes in Cultivated Land Area and Associated Soil and SOC Losses in Northeastern China: The Role of Land Use Policies. \u003cem\u003eInternational Journal of Environmental Research and Public Health 18(21)\u003c/em\u003e\u003cem\u003e,\u003c/em\u003e 11314. https://doi.org/10.3390/ijerph182111314\u003c/li\u003e\n\u003cli\u003eHessel, R., Jetten, V. G., Liu, B., \u0026amp; Qiu, Y. (2011). Evaluating effects of soil and water management and land use change on the Loess Plateau of China using LISEM. \u003cem\u003eHandbook of Erosion Modelling. Chichester, West Sussex: Wiley-Blackwell,\u003c/em\u003e 223-248. https://doi.org/10.1002/9781444328455.ch12\u003c/li\u003e\n\u003cli\u003eHumphrey, O. S., Osano, O., Aura, C. M., Marriott, A. L., Dowell, S. M., Blake, W. H., et al. (2022). Evaluating spatio-temporal soil erosion dynamics in the Winam Gulf catchment, Kenya for enhanced decision making in the land-lake interface. \u003cem\u003eScience of the Total Environment, 815, \u003c/em\u003e151975. https://doi.org/10.1016/j.scitotenv.2021.151975\u003c/li\u003e\n\u003cli\u003eHu, J., Ma, J., Nie, C., Xue, L., Zhang, Y., Ni, F., et al. (2020). Attribution Analysis of Runoff change in Min-tuo River Basin based on SWAT model simulations, china. \u003cem\u003eScientific reports, 10(1),\u003c/em\u003e 2900. doi: 10.1038/s41598-020-59659-z\u003c/li\u003e\n\u003cli\u003eZhang, J., Li, H., Zhang, X., Yu, H., Liang, F., Qin, Z., et al. (2020). Sensitivity evaluation of soil erosion based on land use types: A case study of Minjiang River Basin. \u003cem\u003eJournal of Intelligent \u0026amp; Fuzzy Systems, 38(5), \u003c/em\u003e5697-5705. https://doi.org/10.3233/JIFS-179658\u003c/li\u003e\n\u003cli\u003eKeno, B., \u0026amp; Suryabhagavan, K. V. (2014). Multitemporal remote sensing of landscape dynamics and pattern change in Dire district, southern Ethiopia. \u003cem\u003eJournal of Geomatics, 8(2),\u003c/em\u003e 189-194. doi: 10.4172/2157-7617.1000226\u003c/li\u003e\n\u003cli\u003eLiu, J., Ni, F., Deng, Y., Li, L., Wang, Z., Zhou, D., et al. (2019). Risk assessment of soil erosion in the upper reaches of Minjiang Rive. \u003cem\u003eSouth-to-North Water Transfers and Water Science \u0026amp; Technology, 17(01),\u003c/em\u003e 105-112. doi: 10.13476/j.cnki.nsbdqk.2019.0015.\u003c/li\u003e\n\u003cli\u003eLu, K., Arshad, M., Ma, X., Ullah, I., Wang, J., \u0026amp; Shao, W. (2022). Evaluating observed and future spatiotemporal changes in precipitation and temperature across China based on CMIP6‐GCMs. \u003cem\u003eInternational Journal of Climatology, 42(15),\u003c/em\u003e 7703-7729. https://doi.org/10.1002/joc.7673\u003c/li\u003e\n\u003cli\u003eMinistry of Water Resources of the People\u0026apos;s Republic of China. (2008). \u003cem\u003eSL190-2007\u003c/em\u003e, classification and grading standards for soil erosion. Beijing: China Water \u0026amp; Power Press. \u003c/li\u003e\n\u003cli\u003eMinistry of Water Resources of the People\u0026apos;s Republic of China. (2018). \u003cem\u003eSL773-2018\u003c/em\u003e, guidance on estimation of soil loss in production and construction projects. Beijing: China Water \u0026amp; Power Press. \u003c/li\u003e\n\u003cli\u003eMirzabaev, A., Goedecke, J., Dubovyk, O., Djanibekov, U., Le, Q. B., \u0026amp; Aw-Hassan, A. (2016). Economics of land degradation in Central Asia. \u003cem\u003eEconomics of land degradation and improvement\u0026ndash;A global assessment for sustainable development,\u003c/em\u003e 261-290. doi: 10.1007/978-3-319-19168-3_10\u003c/li\u003e\n\u003cli\u003eNearing, M. A., Simanton, J. R., Norton, L. D., Bulygin, S. J., \u0026amp; Stone, J. (1999). Soil erosion by surface water flow on a stony, semiarid hillslope. Earth Surface Processes and Landforms: \u003cem\u003eThe Journal of the British Geomorphological Research Group, 24(8),\u003c/em\u003e 677-686. https://doi.org/10.1002/(SICI)1096-9837(199908)24:8\u0026lt;677::AID-ESP981\u0026gt;3.0.CO;2-1\u003c/li\u003e\n\u003cli\u003ePeng, J., Li, D. D., \u0026amp; Zhang, Y. Q. (2007). Analysis of spatial characteristics of soil erosion in mountain areas of northwestern Yunnan based on GIS and RUSLE. \u003cem\u003eJournal of Mountain Science, 25(5),\u003c/em\u003e 548-556. ISSN:1008-2186/CN:51-1516\u003c/li\u003e\n\u003cli\u003ePeng, S., Ding, Y., Wen, Z., Chen, Y., Cao, Y., \u0026amp; Ren, J. (2017). Spatiotemporal change and trend analysis of potential evapotranspiration over the Loess Plateau of China during 2011\u0026ndash;2100. \u003cem\u003eAgricultural and forest meteorology, 233,\u003c/em\u003e 183-194. https://doi.org/10.1016/j.agrformet.2016.11.129\u003c/li\u003e\n\u003cli\u003eWang, Q., \u0026amp; Wang, H. (2022). An integrated approach of logistic-MCE-CA-Markov to predict the land use structure and their micro-spatial characteristics analysis in Wuhan metropolitan area, Central China. \u003cem\u003eEnvironmental Science and Pollution Research,\u003c/em\u003e 1-24. doi: 10.1007/s11356-021-17750-6\u003c/li\u003e\n\u003cli\u003eRao, L., Xu, Y., Hu, J., Wang, Z., Yao, W., \u0026amp; Shen, Z. (2020). Study on Soil Erodibility Factoe K on Soil Cover Area of Piaha Sandstone Region. Journal of Basic Science and Engineering, 28(04), 763-773. doi: 10.16058/j.issn.1005-0930.2020.04.002\u003c/li\u003e\n\u003cli\u003eRenard, K. G., \u0026amp; Ferreira, V. A. (1993). RUSLE model description and database sensitivity. \u003cem\u003eJournal of environmental quality, 22(3),\u003c/em\u003e 458-466. https://doi.org/10.2134/jeq1993.00472425002200030009x\u003c/li\u003e\n\u003cli\u003eRenard, K. G., Yoder, D. C., Lightle, D. T., \u0026amp; Dabney, S. M. (2011). Univeral soil loss equation and revised universal soil loss equation. \u003cem\u003eHandbook of Erosion Modelling, edited by: Morgan R PC and Nearing, MA,\u003c/em\u003e 137-67. https://doi.org/10.1002/9781444328455.ch8\u003c/li\u003e\n\u003cli\u003eDe Roo, A. P. J. (1996). The LISEM project: an introduction. \u003cem\u003eHydrological Processes, 10(8),\u003c/em\u003e 1021-1025. https://doi.org/10.1002/(SICI)1099-1085(199608)10:8\u0026lt;1021::AID-HYP407\u0026gt;3.0.CO;2-I\u003c/li\u003e\n\u003cli\u003eRymszewicz, A., Mockler, E., O\u0026apos;Sullivan, J., Bruen, M., Turner, J., Conroy, E., et al. (2015). Assessing the applicability of the Revised Universal Soil Loss Equation (RUSLE) to Irish Catchments.\u003cem\u003e Proceedings of the International Association of Hydrological Sciences, 367, \u003c/em\u003e99-105. https://doi.org/10.5194/piahs-367-99-2015\u003c/li\u003e\n\u003cli\u003ePiao, S., Zhang, X., Chen, A., Liu, Q., Lian, X., Wang, X., et al. (2019). The impacts of climate extremes on the terrestrial carbon cycle: A review. \u003cem\u003eScience China Earth Sciences, 62,\u003c/em\u003e 1551-1563. doi: 10.1007/s11430-018-9363-5\u003c/li\u003e\n\u003cli\u003eMa, S., Wang, L. J., Wang, H. Y., Zhao, Y. G., \u0026amp; Jiang, J. (2023). Impacts of land use/land cover and soil property changes on soil erosion in the black soil region, China. \u003cem\u003eJournal of Environmental Management, 328,\u003c/em\u003e 117024. https://doi.org/10.1016/j.jenvman.2022.117024\u003c/li\u003e\n\u003cli\u003eSong, S., \u0026amp; Yan, X. (2022). Evaluation of events of extreme temperature change between neighboring days in CMIP6 models over China. \u003cem\u003eTheoretical and Applied Climatology, 150(1-2),\u003c/em\u003e 53-72. doi: 10.21203/rs.3.rs-601223/v1\u003c/li\u003e\n\u003cli\u003eSu, Y., Wu, X., Guo, C., Chen, X., Liu, L., \u0026amp; Peng, S. (2022). The spatiotemporal pattern and trend of annual precipitation over China from 2001 to 2100. \u003cem\u003eJournal of Lanzhou University (Natural Sciences), 58(05),\u003c/em\u003e 641-649. doi: 10.13885/j.issn.0455-2059.2022.05.009\u003c/li\u003e\n\u003cli\u003eWang, Y., Li, H., Wang, H., Sun, B., \u0026amp; Chen, H. (2021). Evaluation of CMIP6 model simulations of extreme precipitation in China and comparison with CMIP5. \u003cem\u003eActa Meteorologica Sinica, 79(03),\u003c/em\u003e 369-386. doi: 10.11676/qxxb2021.031\u003c/li\u003e\n\u003cli\u003eWang, Z., \u0026amp; Su, Y. (2020). Assessment of soil erosion in the Qinba mountains of the Southern Shaanxi Province in China using the RUSLE model. \u003cem\u003eSustainability, 12(5),\u003c/em\u003e 1733. https://doi.org/10.3390/su12051733\u003c/li\u003e\n\u003cli\u003eWilliams, J. R. (1990). The erosion-productivity impact calculator (EPIC) model: a case history. \u003cem\u003ePhilosophical Transactions of the Royal Society of London. Series B: Biological Sciences, 329(1255),\u003c/em\u003e 421-428. https://doi.org/10.1098/rstb.1990.0184\u003c/li\u003e\n\u003cli\u003eWischmeier, W. H., \u0026amp; Smith, D. D. (1958). Rainfall energy and its relationship to soil loss. Eos, \u003cem\u003eTransactions American Geophysical Union, 39(2),\u003c/em\u003e 285-291. doi: 10.1029/TR039I002P00285\u003c/li\u003e\n\u003cli\u003eWischmeier, W. H., \u0026amp; Smith, D. D. (1978). Predicting rainfall erosion losses: a guide to conservation planning (No. 537). \u003cem\u003eDepartment of Agriculture, Science and Education Administration.\u003c/em\u003e http://topsoil.nserl.purdue.edu/usle/AH_537.pdf\u003c/li\u003e\n\u003cli\u003eWu, J., Xia, J., Zeng, S., Liu, X., \u0026amp; Fan, D. (2023). Evaluation of the Performance of CMIP6 Models andFuture Changes Over the Yangtze River Basin. \u003cem\u003eResources and Environment in the Yangtze Basin\u003c/em\u003e, 32(01), 137-150, doi:10.11870/cjlyzyyhj202301013\u003c/li\u003e\n\u003cli\u003eYang, J., \u0026amp; Huang, X. (2021). The 30 m annual land cover dataset and its dynamics in China from 1990 to 2019. \u003cem\u003eEarth System Science Data\u003c/em\u003e, 13(8), 3907-3925, https://doi.org/10.5194/essd-13-3907-2021\u003c/li\u003e\n\u003cli\u003eYang J., Dong J., Xiao X., Dai J. Wu C., Xia J., et al. (2019). Divergent shifts in peak photosynthesis timing of temperate and alpine grasslands in China. \u003cem\u003eRemote Sensing of Environment\u003c/em\u003e, 233, 111395, https://doi.org/10.1016/j.rse.2019.111395\u003c/li\u003e\n\u003cli\u003eWei, Y., Sun, S., Liang, D., \u0026amp; Jia, Z. (2022). Spatial\u0026ndash;temporal variations of NDVI and its response to climate in China from 2001 to 2020. \u003cem\u003eInternational Journal of Digital Earth\u003c/em\u003e, 15(1), 1463-1484, https://doi.org/10.1080/17538947.2022.2116118\u003c/li\u003e\n\u003cli\u003eYao, Y., Li, J., Jiang, Y., \u0026amp; Huang, G. (2023). Evaluating the response and adaptation of urban stormwater systems to changed rainfall with the CMIP6 projections. \u003cem\u003eJournal of Environmental Management\u003c/em\u003e, 347, 119135, https://doi.org/10.1016/j.jenvman.2023.119135\u003c/li\u003e\n\u003cli\u003eZhong, X., Zhang, S., Wu, R., Jing, Y., Men, L., \u0026amp; Zhou, T. (2022). Analysis of Dynamic Changes and Driving Forces of Soil Erosion in Tuojiang River Basin. \u003cem\u003eResearch of Soil and Water Conservation\u003c/em\u003e, 29(02), 43-49+56, doi:10.13869/j.cnki.rswc.2022.02.003\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Soil erosion, CMIP6, RUSLE model, CA-Markov model, Min-Tuo River Basin","lastPublishedDoi":"10.21203/rs.3.rs-4952297/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4952297/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eTo quantify the current status of soil erosion in the Min-Tuo River Basin and predict future soil erosion conditions, the RUSLE model and GIS technology were employed to assess soil erosion in the Min-Tuo River Basin from 2000 to 2020. Based on the CMIP6 model, precipitation in the Min-Tuo River from 2021 to 2050 was predicted, and the CA-Markov model was used to forecast land use types from 2030 to 2050, thereby predicting the spatial distribution trends of soil erosion in 2030, 2040, and 2050. The results indicate that: (1) The overall soil erosion intensity in the study area increased from 2000 to 2020; (2) Future precipitation in the basin, predicted using the CMIP6 model, showed a fluctuating upward trend compared to historical levels. The SSP5-8.5 scenario had the largest fluctuation; (3) In terms of future land use types, an increase in the area of converted cropland may increase the risk of soil erosion; (4) The area affected by soil erosion in the future will increase, with the highest annual average soil erosion modulus and total amount occurring under the SSP5-8.5 high emission scenario; (5) The dynamic evolution of future soil erosion levels indicates that the stability of soil erosion levels under the SSP2-4.5 scenario is higher than under the other scenarios. In the future high CO2 emission scenario, the risk of soil erosion in the Min-Tuo River Basin will increase, necessitating the active implementation of various prevention and control measures to prevent further exacerbation.\u003c/p\u003e","manuscriptTitle":"High greenhouse gas emission scenarios increase soil erosion risk: A case study of the Min-Tuo River Basin, China","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-10-18 08:48:42","doi":"10.21203/rs.3.rs-4952297/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"d63d7664-62eb-4222-a18c-e61a372e0edc","owner":[],"postedDate":"October 18th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2024-12-23T16:23:28+00:00","versionOfRecord":[],"versionCreatedAt":"2024-10-18 08:48:42","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4952297","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4952297","identity":"rs-4952297","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}