The Runoff in the Upper Taohe River Basin and Its Responses to Climate Change

Climate change has a significant impact on water resources. Forecasts and simulations of climate runoff processes are essential for assessing the impact of global climate change on runoff variations. This study focuses on the upper Taohe River Basin, which is an important watershed in the semi-arid regions of northwest China. To assess the runoff in the upper Taohe River Basin and the responses to climate change, the SWAT hydrological model was used to analyze future climate change scenarios and their effects on water resources. The results indicate that the minimum temperature would increase gradually in the 21st century and that the minimum temperature change would be more significant than the maximum temperature change, which indicates that minimum temperature changes would make an obvious contribution to future regional warming. Under RCP2.6, the average precipitation would decrease; at the same time, under RCP4.5 and RCP8.5, the average precipitation would increase. In the future, under different climate scenarios, the runoff will exhibit droughts and flood disasters. These research results provide scientific support for water resource utilization and management in the Taohe River Basin.


Introduction
Water resources have become one of the most important global resources and environmental problems in the 21st century and is also an essential aspect in international Earth science development [1][2][3][4]. Globally, rising temperatures due to continuous greenhouse gas emissions have exerted an irreversible impact on climate change, water resources, agriculture, health, energy, and natural systems. In the meantime, it has also added increased risk to the ecological environment as well as economic and social development in the future [5][6][7][8][9][10][11]. How to scientifically and quantitatively assess the impact of future climate change on regional hydrological processes has become a hot topic for meteorological and hydrological researchers [12,13]. Generally, the most scientific approach to predict climate change is based on physical assumptions that describe future global or regional climate change or climate change scenarios [14,15].
The Coupled Model Inter-comparison Project 5 (CMIP5) from the 5th Assessment Report (AR5) is the most successful of the General Circulation Models (GCMs)and includes four future scenarios: Representative Concentration Pathways (RCPs) 2.6, 4.5, 6.0, and 8.5. These scenarios are based on the expected differences in radiation, which have great impact on the future climate. In RCP2.6, scenarios can be found with radiative forcing as low as 3 W/m 2 in the year 2100. RCP4.5 and RCP6.0 assume the stabilization of radiative forcing. However, RCP8.5 assumes increased radiative forcing to 8.5 W·m −2 after 2100 [16]. The upper Taohe River basin is an important water concentration area in northwest China. It is susceptible to the influence of the southwest monsoon originating from the Indian Ocean, and receives abundant precipitation of 400~600 mm per year. Therefore, the natural runoff gradually occupies more than 68.8% of the whole basin, which is the main runoff generation area and an important water supply area of the Taohe River Basin [31][32][33][34]. The water systems in the study area are well developed and have many symmetrical tributaries ( Figure  1c). Generally, the middle and upper reaches of the Taohe River are important ecological barriers in the upper reaches of the Yellow River and play an important role in maintaining the water resources and ecological security in the Yellow River Basin.

Climate Change Scenario Data
The future climate change scenario data used in this paper were derived from the RegCM4.6, which is based on the future climate prediction dataset of the National Tibetan plateau Science Data Center (http://data.tpdc.ac.cn/zh-hans/data/ accessed on 18 March 2022) for northwest China. This dataset is based on four different greenhouse gas emission concentrations (RCP2.6, RCP4.5, RCP6.0, and RCP8.5) determined by the HadGEM2-ES of the regional climate model RegCM4.6, which can simulate and predict the average temperature and precipitation in northwest China from 2007 to 2099 [35]. The future climate change output data set has a spatial resolution of 0.25° × 0.25° and 3 h, daily, and yearly temporal intervals that cover northwest China (Gansu, Ningxia, Qinghai, Xinjiang, and Shaanxi Provinces). This study selected climate variables during the period from 2007 to 2099 under three different emission scenarios (RCP2.6, RCP4.5 and RCP8.5) in the middle and upper reaches of the Taohe River basin/UTB as the forcing data for the SWAT model to project future runoff changes. The upper Taohe River basin is an important water concentration area in northwest China. It is susceptible to the influence of the southwest monsoon originating from the Indian Ocean, and receives abundant precipitation of 400~600 mm per year. Therefore, the natural runoff gradually occupies more than 68.8% of the whole basin, which is the main runoff generation area and an important water supply area of the Taohe River Basin [31][32][33][34]. The water systems in the study area are well developed and have many symmetrical tributaries ( Figure 1c). Generally, the middle and upper reaches of the Taohe River are important ecological barriers in the upper reaches of the Yellow River and play an important role in maintaining the water resources and ecological security in the Yellow River Basin.

Climate Change Scenario Data
The future climate change scenario data used in this paper were derived from the RegCM4.6, which is based on the future climate prediction dataset of the National Tibetan plateau Science Data Center (http://data.tpdc.ac.cn/zh-hans/data/ accessed on 18 March 2022) for northwest China. This dataset is based on four different greenhouse gas emission concentrations (RCP2.6, RCP4.5, RCP6.0, and RCP8.5) determined by the HadGEM2-ES of the regional climate model RegCM4.6, which can simulate and predict the average temperature and precipitation in northwest China from 2007 to 2099 [35]. The future climate change output data set has a spatial resolution of 0.25 • × 0.25 • and 3 h, daily, and yearly temporal intervals that cover northwest China (Gansu, Ningxia, Qinghai, Xinjiang, and Shaanxi Provinces). This study selected climate variables during the period from 2007 to 2099 under three different emission scenarios (RCP2.6, RCP4.5 and RCP8.5) in the middle and upper reaches of the Taohe River basin/UTB as the forcing data for the SWAT model to project future runoff changes.

Meteorological and Hydrological Data
Data from eight major meteorological stations in the middle and upper reaches of the Taohe River and its surrounding areas were selected as the input data for the hydrological model. Information on major weather stations is shown in Table 1. The meteorological station data during 1986-2018 were obtained from the China National Meteorological Data Sharing Website (http://data.cma.cn/ accessed on 14 August 2020). The meteorological elements include wind speed, evaporation, atmospheric pressure, temperature (average temperature, maximum temperature, minimum temperature), average vapor pressure, average relative humidity, precipitation, and sunshine hours. The monthly runoff data of the representative hydrological stations (Luqu, Xiabagou, and Minxian) in the middle and upper reaches of the Tao River from 1986 to 2014 were provided by the Gansu Hydrology and Water Resources Bureau (Table 2).

SWAT Model Construction Data
The soil and water assessment tool (SWAT) is a well-known hydrological modeling tool that has been applied in various hydrologic simulations. It has been applied to analyze soil water conservation structures and their impacts on reducing runoff. The SWAT model has been demonstrated to be a powerful tool when selecting the most technically effective management strategies to reduce soil degradation [36].
The upper Taohe River Basin can be extracted and divided into 28 subbasins and 388 hydrological response units (HRUs) from DEM. Monthly runoff data from three hydrological stations from 1986 to 2010 are regarded as the calibration period, and data from the period from 2011 to 2014 are used for verification. SWAT-CUP is used for parameter sensitivity analysis and uncertainty analysis.
The land use data used in this research are from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn accessed on 10 may 2021). There are some differences between China's land use classification system and the SWAT model, and it is necessary to reclassify regional land use data. After reclassification, there are 6 first-level types and 19 s-level types of land use that can be produced in the upper Taohe River Basin. They are divided into nine new land types (Table 3) to highlight the main land use distribution characteristics of the basin. from the Nanjing Soil Research Institute, Chinese Academy of Sciences, v1.1 with a 1 km resolution was used to provide reliable soil parameters for SWAT. The SWAT model has collected a large number of physical and mathematical equations to describe the hydrological processes and material transport processes, and many parameters have been applied, most of which have clear hydrological physical significance, main parameters are shown in Table 4. Note: Calibration method: r_ indicates that the existing parameters will be multiplied by 1 (given value), a_ indicates that the given value will be added to the existing parameters, v_ indicates that the existing parameters will be replaced by the given value.

Evaluation Indicators
The calibration and verification of the SWAT model was carried out using the Nash efficiency coefficient (NSE), certainty coefficient R 2 , and relative deviation (RE) as evaluation indexes to evaluate the simulation results with the measured values. The NSE coefficient reflect the fitting degree between the observed and the simulated values. The closer the NSE coefficient is to 1, the better the simulation effect is. If the NSE coefficient is greater than 0.5, the simulation of the model is considered successful. R 2 is used to characterize the correlation degree of the variables, which is used to evaluate the consistency of the change trend between the simulated value and the measured value. The calculated value tends to 1, indicating that the simulation effect is better. It is generally considered that R 2 > 0.6 can be used as a criterion for evaluating the correlation between simulated and measured runoff values. Re represents the relative deviation between the simulated value and the measured value. The more the calculated value tends to 0, the better the effect is. It is generally believed that Re < 20%, and the simulation results can be accepted [37]. The formulas for the evaluation indicators are presented as Equations (1)-(3).
n represents the length of the simulation time; Q obs and Q sim represent the average values of the runoff observation and simulation values in the simulation period, respectively; Q obs and Q sim represent the runoff observation and simulation values in the research time, respectively.  [35] found that the temperature bias of HadGEM2-ES is generally within ±2.5 • C in the southeast and in the south during the historical period of 1985-2004. This article has evaluated and revised the temperature output values under the RCP4.5 scenario based on the observed meteorological patterns in the upper Tao River Basin based on the historical period of 2007-2018. The deviation between the simulated value and the measured value is exhibited in Figure 2. The deviation between the simulated value and the measured value is small in summer and large in winter. The maximum temperature deviation can reach 7.5 • C, but the temperature variation is consistent with the observed value ( Figure 2). According to the relationship between the altitude and temperature in each site, the average observed temperature from 2007 to 2018 is about 0.8 • C. In comparison, the simulated annual average temperature and the observed multiyear average are largely consistent under the different climate change scenarios. Corrections through inverse deductions between the linear equation and the vertical decline rate of the temperature, the average temperature, and the maximum/minimum temperature are adjusted separately. Therefore, the corrected data are more in line with the actual situation of the study area and meet the research needs of future climate change scenarios.

Result Analysis
2007 to 2018 is about 0.8 °C. In comparison, the simulated annual average temperature and the observed multiyear average are largely consistent under the different climate change scenarios. Corrections through inverse deductions between the linear equation and the vertical decline rate of the temperature, the average temperature, and the maximum/minimum temperature are adjusted separately. Therefore, the corrected data are more in line with the actual situation of the study area and meet the research needs of future climate change scenarios. Under the three greenhouse gas emission scenarios, the average annual temperature in the study area shows a consistent warming trend in the future period ( Figure 3). The average temperature in the future study area would be about 2.83 °C, 3.32 °C, and 4.24 °C under RCP2.6, RCP4.5, and RCP8.5 during 2007-2100, and the change rates may be 0.10 °C/10a, 0.20 °C/10a, and 0.52 °C/10a. The average temperature in the 2080s would be about 0.54 °C, 1.14 °C, and 3.44 °C higher than the average temperature in 1956-1997 [38], which may be consistent with the global warming trends. Overall, the temperature in the upper Taohe River Basin increases with the increase in emission scenarios. Under the three greenhouse gas emission scenarios, the average annual temperature in the study area shows a consistent warming trend in the future period ( Figure 3). The average temperature in the future study area would be about 2.83 • C, 3.32 • C, and 4.24 • C under RCP2.6, RCP4.5, and RCP8.5 during 2007-2100, and the change rates may be 0.10 • C/10a, 0.20 • C/10a, and 0.52 • C/10a. The average temperature in the 2080s would be about 0.54 • C, 1.14 • C, and 3.44 • C higher than the average temperature in 1956-1997 [38], which may be consistent with the global warming trends. Overall, the temperature in the upper Taohe River Basin increases with the increase in emission scenarios. At the same time, the possible future temperature changes in this basin are also analyzed from the two aspects of maximum and minimum temperatures. The annual average maximum and minimum temperature changes in the upper and middle reaches of the Taohe River under the three greenhouse gas emission scenarios were set for four stages:   The maximum temperature and minimum temperature show increasing trends under the three scenarios. The change in the maximum temperature under the RCP2.6 scenario increases steadily, and the increments in the three stages would be 0.13 • C, 0.44 • C, and 0.44 • C higher than those from 2007 to 2018, which is consistent with the change trends in the average temperature. Under the RCP4.5 scenario, the maximum temperature would increase significantly, and the temperature could be 0.63 • C, 1.02 • C, and 1.42 • C higher than in 2007-2018, under the three scenarios, with a maximum change range of 13-17%. The future maximum temperature shows significant changes under the RCP8.5 scenario, which would increase to 0.49 • C, 1.64 • C, and 3.34 • C, which could be higher than the temperature in 2007-2018, with a maximum range of 33%. The results indicate that the highest temperature in the 21st century would increase gradually.
Compared to the maximum temperature variations, the change range of the minimum temperature in different scenarios is consistent with the maximum temperature in the future. Under the RCP2.6 scenario, the minimum temperature in the three stages would increase to 0.26 • C, 0.61 • C, and 0.60 • C, which would be slightly higher than that in 2007-2018, respectively. The change under the RCP4.5 scenario would be significantly enhanced, and the temperatures of the three stages would increase to 0.71 • C, 1.12 • C, and 1.59 • C, which would be higher than those from 2007 to 2018. Under the RCP8.5 scenario, the minimum temperature changed significantly, the minimum temperature would be 0.64 • C, 1.79 • C, and 3.57 • C higher than those from 2007 to 2018. The results indicate that the minimum temperature will increase gradually in the 21st century and that the minimum temperature change would be more significant than the maximum temperature change, which indicates that the minimum temperature changes make an obvious contribution to future regional warming.

Projection of Future Precipitation Change in Upper Taohe River Basin
According to the precipitation data output by the climate model in the historical period (1985~2015) in the northwest region, the simulation effect in the eastern of Qinghai-Tibet Plateau is poor, which may be due to the influence of monsoon circulation on the Qinghai-Tibet Plateau, resulting in a false high-value precipitation center in the climate model. Compared to the observed precipitation from 2007 to 2018 [35], it was found that the simulated precipitation from the climate model is similarly overestimated in the eastern part of the middle and upper reaches of the Taohe River. Therefore, the seriously overestimated stations were removed and compared with the precipitation changes observed by all of the stations in the region. Figure 5 shows the future multiyear precipitation changes in the middle and upper reaches of the Taohe River from 2007 to 2099 and compares the not-removed and removed gridded precipitation. The moving average is a reasonable and practical way to analyze the meteorological data. The 5-year moving average has been selected as a default estimator for the annual survey, partially because it seems easy to understand and compute. The results show that the 5-year moving average curve could be presented as the precipitation change trend.

Applicability Evaluation of SWAT Model
The results of hydrological process curves during the calibration period and verification period showed that SWAT model can capture the time and flow of the flood peaks in the three hydrological stations well ( Figure 6). The simulated value in the dry season is also consistent with the basic flow of the basin, but the simulation abilities among the three stations are different. The Luqu, Xiabagou, and Minxian hydrological stations are the main outlets of the source area in the Taohe River. Observational data showed that from 1986 to 2014, the runoff from the Luqu, Xiabagou, and Minxian hydrological stations in the Taohe River Basin showed a significant response to precipitation. From 1986 to 2003, the runoff showed a slight decrease along with precipitation. During 2003 to 2011, the runoff decreased as the precipitation decreased; then, from 2011 to 2014, the runoff increased as the precipitation increased ( Figure 6). If we take 1986 to 2010 as the model calibration period and 2011 to 2014 as the model validation period, the SWAT model can capture the variation trends in the three hydrological stations well. The overall runoff simulation values of the determination coefficient R 2 , the Nash efficiency coefficient NSE, and the relative error Re of the monthly runoff simulation at regular rates fall within a small uncertainty interval (Table 5). Under the RCP2.6 scenario, the average precipitation would be about 658 mm in the study area in 2020-2099, demonstrating an insignificant decreasing trend (−3.69 mm/10a), and the average precipitation after excluding abnormal stations would be 620 mm, which is closer to the measured precipitation. The precipitation would fluctuate greatly in the 2030s, while the precipitation would achieve its lowest level in the 2080s. Under the RCP4.5 scenario, the average precipitation in the future would be 677.5 mm, and the average precipitation after excluding abnormal stations would be 638.5 mm, demonstrating an insignificant increasing trend (4.97 mm/10a). The precipitation would be the highest in the 2070s (689 mm) and would reach its lowest value in the 2020s (591.3 mm). The future average precipitation under the RCP8.5 scenario would reach 693.2 mm, and the average precipitation after excluding abnormal stations would be 653.9 mm, demonstrating an increasing trend (12.28 mm/10a). The average precipitation would experience an insignificant change before the 2080s, but significant fluctuations would appear after the 2080s.

Applicability Evaluation of SWAT Model
The results of hydrological process curves during the calibration period and verification period showed that SWAT model can capture the time and flow of the flood peaks in the three hydrological stations well ( Figure 6). The simulated value in the dry season is also consistent with the basic flow of the basin, but the simulation abilities among the three stations are different. The Luqu, Xiabagou, and Minxian hydrological stations are the main outlets of the source area in the Taohe River. Observational data showed that from 1986 to 2014, the runoff from the Luqu, Xiabagou, and Minxian hydrological stations in the Taohe River Basin showed a significant response to precipitation. From 1986 to 2003, the runoff showed a slight decrease along with precipitation. During 2003 to 2011, the runoff decreased as the precipitation decreased; then, from 2011 to 2014, the runoff increased as the precipitation increased ( Figure 6). If we take 1986 to 2010 as the model calibration period and 2011 to 2014 as the model validation period, the SWAT model can capture the variation trends in the three hydrological stations well. The overall runoff simulation values of the determination coefficient R 2 , the Nash efficiency coefficient NSE, and the relative error Re of the monthly runoff simulation at regular rates fall within a small uncertainty interval (Table 5). Luqu station, which is located at the source of the Taohe River Basin, the simula performance of this station can be easily observed, with R 2 , NSE, and Re showing mont runoff simulation rates of 0.79, 0.89, and 2.39%, respectively. In the verification period, NSE and Re also reached 0.89, 0.95, and −8.9%, respectively. The R 2 , NSE, and Re of monthly runoff simulation at the Xiabagou station were 0.77, 0.88, and −3.86%, resp  Luqu station, which is located at the source of the Taohe River Basin, the simulated performance of this station can be easily observed, with R 2 , NSE, and Re showing monthly runoff simulation rates of 0.79, 0.89, and 2.39%, respectively. In the verification period, R 2 , NSE and Re also reached 0.89, 0.95, and −8.9%, respectively. The R 2 , NSE, and Re of the monthly runoff simulation at the Xiabagou station were 0.77, 0.88, and −3.86%, respectively, and 0.89, 0.96, −5.07% during the verification period. The Minxian station serves as the total outlet of the upper Taohe River Basin.
The results showed that monthly runoff simulation rates of R 2 , NSE, and Re were 0.83, 0.91, −14.6%, respectively, and during the verification period, R 2 , NSE, and Re also reached 0.87, 0.94, −8.6%, respectively. The above assessment results indicate that the distributed hydrological model SWAT is feasible to simulate runoff in the middle and upper reaches of Taohe River, which lays a foundation for the subsequent study, which is focus on the response of water resources to climate change in the Taohe River Basin.
The calibration and validation results of three on-site observations showed that the SWAT hydrological model is able to produce an acceptable simulation of runoff at a monthly time step, producing reliable results and meeting the research requirements.

Projection of Future Runoff Change in the Middle and Upper Reaches of the Tao River
Based on the good application of the SWAT model, the annual runoff changes under the three greenhouse gas emission scenarios in the middle and upper reaches of the Taohe River could be predicted from 2020 to 2099 by inputting the corrected RCP temperatures and eliminating abnormal precipitation data from the grids. As shown in Figure 7, to compare the long-term runoff changes, this research takes the average runoff trends from 1956 to 2014 as the historical period, allowing the runoff changes in different future periods to be analyzed intuitively.
Under the RCP2.6 scenario, the annual average runoff would be about 30.9 × 10 8 m 3 , and the overall change trend is similar to the runoff from 2003 to 2014. Over the whole period, the relative maximum and minimum runoff would alternately appear in the 1930s, while the overall minimum runoff would appear in the mid-1980s and would be as low as 17.1 × 10 8 m 3 . Comparably, the overall average runoff would be the highest in the 1950s, with an annual average runoff of 35.2 × 10 8 m 3 , and would be the lowest in the 1980s at 24.6 × 10 8 m 3 . Drought risk might be estimated in the future. Under the RCP4.5 scenario, the annual average runoff would be about 32.5 × 10 8 m 3 , and the overall change trend is 15% lower than during the period of 1956−1985. The highest runoff values might be observed during the 2040s and 2070s. Under the RCP8.5 scenario, the annual average runoff would be about 32.2 × 10 8 m 3 , with no significant increasing trend being observed. The maximum runoff would appear in the mid-2080s, and the minimum would appear in the mid-2060s and would be as low as 14.8 × 10 8 m 3 . 24.6 × 10 8 m 3 . Drought risk might be estimated in the future. Under the RCP4.5 scenario, the annual average runoff would be about 32.5 × 10 8 m 3 , and the overall change trend is 15% lower than during the period of 1956−1985. The highest runoff values might be observed during the 2040s and 2070s. Under the RCP8.5 scenario, the annual average runoff would be about 32.2 × 10 8 m 3 , with no significant increasing trend being observed. The maximum runoff would appear in the mid-2080s, and the minimum would appear in the mid-2060s and would be as low as 14.8 × 10 8 m 3 .  Table 6 shows the projected runoff changes in each season. Under the RCP2.6 scenario, runoff would decrease significantly in summer, and insignificant changes would be  Table 6 shows the projected runoff changes in each season. Under the RCP2.6 scenario, runoff would decrease significantly in summer, and insignificant changes would be observed during the other seasons. Under the RCP4.5 and RCP8.5 scenarios, the runoff in all seasons would show a fluctuating trend. In general, the future runoff in the upper Taohe River Basin would show a decreasing trend in summer and increasing in autumn. In other seasons, significant fluctuations can be observed the future runoff.

Uncertainty Analysis of SWAT Model
The SWAT model has a complex structure and involves many equations and variables [39]. Due to the complexity and randomness of hydrological processes, there are many uncertain factors in the processes of hydrological model simulation that could cause interference to the simulation effect. According to the observed situation in the upper Taohe River Basin and in the relevant literature, 23 parameters in the SWAT model were selected for the preliminary overall sensitivity analysis. The calibration number was set to 500 times, and the iteration number was set to 10 times. The parameters suitable for the hydrological process changes in the study area were selected carefully and thoroughly. The contribution of different parameters to the calibration and validation of the model in different study areas are also different and can have a varying degree of impact on the model results [40,41].

Uncertainty Analysis of Future Climate Change Scenarios
In this study, RegCM4.6 driven by the global climate model HadGEM2-ES was used to simulate future climate change in northwest China. Pan et al. [35] compared the historical period  and found that the output results of the model had relatively large simulation errors in the middle and upper reaches of the Taohe River, regardless of the average temperature or precipitation. Because the Taohe River is located at the junction of the Tibetan Plateau and the Loess Plateau, it is jointly affected by the southwest monsoon and the East Asian monsoon as well as by the climate fluctuations of the Tibetan Plateau, and the atmospheric circulation changes are also complex. Together with the rugged terrain and inhomogeneous underlying surface characteristics, it is difficult to obtain accurate model simulations [42]. Moreover, the meteorological stations in the upper and middle reaches of the Taohe River and its surrounding areas are relatively scarce, and the measured data are limited, greatly affecting the reliability of the simulated values. Under the RCP4.5 scenario, the precipitation in the middle and upper reaches of the Taohe River showed a significant decreasing trend from the late 2070s to the late 2080s, with a change rate of −8.97 mm/10a; under RCP8.5, the precipitation showed a significant upward trend, with a change rate of 6.50 mm/10a during the same period. This uncertainty stems from the differences in the scenarios, which will increase the uncertainty of future precipitation trends.
Future climate prediction is complex and highly uncertain, increasing the uncertainty in the exploration of the impact of climate change on watershed hydrological processes. There are many sources of uncertainty in climate prediction, including different future greenhouse gas emission scenarios, natural variability within the climate system, and climate process description methods [43,44]. Even though scenarios and models cannot completely overcome these uncertainties, to deepen the understanding of various physical processes in the climate system, the output results of climate models will become increasingly accurate.

Research Limitation and Implication
We have deployed instruments in regions (such as the Qilian Mountains) containing the cryosphere, and the time series of observation data are too short to support a complete watershed runoff simulation. It needs to be emphasized that the study area selected in this article is in a semi-arid area, where the distribution of cryosphere elements is negligible and their impact on runoff can be basically ignored, so runoff in this area is mainly affected by climate (temperature and precipitation). In the future, it is necessary to select areas covered by cryospheric elements for further comparative analysis to separate the contributions of cryospheric elements and climate on watershed runoff.
The precipitation results of other models have some common results and some noncommon ones, but the precipitation results of the RCMs are quite heterogeneous. Different simulation ranges, physical parameterization schemes, dynamic frameworks, and atmospheric boundary condition data will have a greater impact on precipitation [19]. Conditional data will have a greater impact on precipitation. This paper is an attempt to use the simulation results of RCM to drive SWAT for runoff simulation, deepening our understanding of climate change impact on runoff changes in the Taohe River basin. In the future, multiple RCMs need to be used for comprehensive evaluation to provide more reasonable runoff simulation results.

Conclusions
In this paper, based on the on-site observations of meteorological and hydrological stations along with the output data from future climate change scenarios and the SWAT hydrological model, future regional climate change and its impact on runoff changes were comprehensively analyzed, and the following conclusions were obtained: Due to the differences in the emission scenarios, the warming amplitude is different. From 2020 to 2099, the average temperature change rates under the RCP2.6, RCP4.5, and RCP8.5 scenarios were 0.10 • C/10a, 0.20 • C/10a, and 0.54 • C/10a, respectively. In the RCP2.6 emission scenario, the average precipitation shows a decreasing trend in the future, with a reduction rate of 3.69 mm/10a. Under the RCP4.5 and RCP8.5 scenarios, the future precipitation shows an increasing trend, with increasing rates of 4.97 mm/10a and 12.28 mm/10a, respectively.
The calibration and validation results of the three on-site observations (Luqu, Xiabagou, and Minxian) in the upper Taohe River Basin showed that the SWAT hydrological model is able to produce an acceptable runoff simulation at a monthly time step. The

Data Availability Statement:
The data in this study is confidential. Please contact the corresponding author if necessary.

Conflicts of Interest:
The authors declare no conflict of interest.