Projection of Future Climate Change and Its Inﬂuence on Surface Runoff of the Upper Yangtze River Basin, China

: Global climate change will modify precipitation and temperatures’ temporal and spatial distribution, trigger more extreme weather events, and impact hydrological processes. The Yangtze River basin is one of the world’s largest basins, and understanding future climate changes is vital for water resource management and supply. Research on predicting future climate change in the upper Yangtze River basin (UYRB) and introducing machine learning algorithms to analyze the impact of climate factors, including extreme weather indicators, on surface runoff is urgently needed. In this study, a statistical downscaling model (SDSM) was used to forecast the future climate in the UYRB, and the Mann–Kendall (MK) or modiﬁed Mann–Kendall (MMK) trend test at a 5% level of signiﬁcance was applied to analyze temporal trends. The Spearman rank correlation (SRC) test at a 5% level of signiﬁcance and random forest regression (RFR) model were employed to identify the key climatic factors affecting surface runoff from annual precipitation, annual temperature, maximum 5-day precipitation (R × 5Day), number of tropical nights (TR), and consecutive dry days (CDD), and the RFR model was also used to predict future runoff. Based on the results, we found that, compared to the selected historical period (1985–2014), the mean annual precipitation (temperature) during the mid-term (2036–2065) increased by 18.93% (12.77%), 17.78% (14.68%), 20.03% (17.03%), and 19.67% (19.29%) under SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5, respectively, and during the long term (2071–2100), increased by 19.44% (12.95%), 22.01% (21.37%), 30.31% (30.32%), and 34.48% (37.97%), respectively. The warming and humidiﬁcation characteristics of the northwestern UYRB were more pronounced. The key climatic factors inﬂuencing surface runoff were annual precipitation, maximum 5-day precipitation (R × 5day), and annual temperature. Because of warming and humidiﬁcation, surface runoff in the UYRB is expected to increase relative to the historical period. The surface runoff during the mid-term (long term) increased by 12.09% (12.58%), 8.15% (6.84%), 8.86% (8.87%), and 5.77% (6.21%) under SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5, respectively. The implementation of sustainable development pathways under the low radiative forcing scenario can be effective in mitigating climate change, but at the same time, it may increase the risk of ﬂoods in the UYRB.


Introduction
Global climate change has become a significant challenge facing humanity and has received extensive attention from society and academia [1,2].The sixth assessment report of the Intergovernmental Panel on Climate Change (IPCC) pointed out that climate change has caused tremendous damage and increasingly irreversible losses to terrestrial, freshwater, coastal, and marine ecosystems [3][4][5][6][7][8][9][10]. As climatic conditions are closely related to the hydrological cycle, it is important to study climate change for regional and global water resource management and utilization [11][12][13][14].
Temperature and precipitation are the basic meteorological factors [15].Climate change has changed the spatial and temporal distribution patterns of global temperature and precipitation and increased the probability of extreme weather events in more than 80% of the world [16].Climate change exhibits global and regional characteristics influenced by various factors.Different regions have distinct climate change features.In areas with complex topographies and volatile local climates, the regional natural conditions are more sensitive to climate change, making the resulting impacts worthy of further research.
Climate change affects the hydrological cycle [17,18].Global warming is accelerating the hydrological cycle and changing the intensity and frequency of precipitation.Moreover, precipitation is directly involved in the hydrological cycle, changing the runoff of the basin [16,19,20].The increase in temperature has changed the water resources' temporal and spatial distribution and increased the probability of extreme climate events [21,22].
Studies on rainfall-runoff modeling show that precipitation is an important recharge source for surface runoff, and regional surface runoff is highly sensitive to changes in precipitation [23][24][25].The increase in temperature leads to an increase in the evaporation of soil moisture and surface runoff and a decrease in streamflow [26][27][28].Climate change has enhanced the temporal and spatial variability of precipitation and temperature, leading to the frequent occurrence of extreme weather events.Extreme weather events such as hot extreme rainstorms are more likely to cause drought and flood disasters [29], which greatly reduce the utilization capacity of regional water resources and cause significant damage to the ecological environment and human social economy [11,16,30].
Extreme climatic factors need more attention in studies on surface runoff and climate.Existing studies are usually conducted at the administrative scale, and few studies have been conducted at the large watershed scale where climatic conditions are complex [31,32].In addition, most studies have utilized statistical methods, such as linear regression, multiple linear regression, or nonparametric tests, to identify correlations between climatic factors and surface runoff and to identify factors that influence surface runoff [27,[33][34][35].However, these methods make it difficult to accurately capture the relationship between climatic factors and surface runoff.This results in low confidence in the identification of key factors affecting surface runoff.
Machine learning is a technique that utilizes algorithms to learn from data and automatically improve the performance of a system.In machine learning, algorithms learn laws and patterns from input data and make predictions or decisions about new data based on these laws and patterns [36].Machine learning techniques can be categorized into three types: supervised learning, unsupervised learning, and reinforcement learning.With the enhancement of the computing power of devices and the increase in data volume, machine-learning algorithms have been widely applied in the field of hydrometeorology, including the establishment of prediction models, decision-making classifications, and data mining [37][38][39][40][41][42].In the study of capturing the correlation between surface runoff and climatic factors and deciding the key factors affecting surface runoff, the application of machine-learning algorithms is expected to enhance the credibility of the conclusions.
The Yangtze River, the longest river in Asia and the third longest in the world, spans 6300 km [8,43].The upper Yangtze River basin (UYRB) is characterized by complex terrain and is significantly influenced by the Qinghai-Tibet Plateau, southwest monsoon, and southeast monsoon.The region is sensitive and vulnerable to climate change [44].Against the backdrop of global warming, the frequency and intensity of droughts and floods in the upstream Yangtze River region have increased.Aquatic and terrestrial ecosystems are deteriorating, leading to frequent natural disasters.These factors severely impede the sustainable socioeconomic development of the Yangtze River basin [45].In the future, climate change will likely alter the hydrological cycle and runoff processes in the UYRB.This has a more pronounced impact on the distribution of water resources in the UYRB.It also adversely affects the water resources for hydropower generation and agricultural irrigation.
At present, research on predicting future climate change in the UYRB and introducing machine learning algorithms to analyze the impact of climate factors, including extreme weather indicators, on surface runoff is urgently needed.To achieve these goals, this study utilized meteorological observation data from 1985 to 2014 in the UYRB and general circulation model (GCM) data under shared socioeconomic pathways (SSP) scenarios.A statistical downscaling model, a random forest regression model, and a multiple regression model were employed to predict future precipitation, temperature, and runoff in the UYRB.The specific research objectives were as follows: (1) to analyze the spatiotemporal precipitation and temperature patterns in the UYRB; (2) to forecast spatiotemporal changes in precipitation and temperature in the UYRB under different SSP scenarios; (3) to identify the key climatic factors influencing runoff in the UYRB and predict trends of surface runoff changes.

Research Area
The Yangtze River originates in the Tanggula Mountains in the Qinghai-Tibet Plateau, with a total length of 6300 km.Yichang, Hubei Province, serves as the dividing point between the upper and lower reaches of the Yangtze River.The upper Yangtze River basin (UYRB) is located between 90 • 32 ~111 • 33 E and 24 • 42 ~35 • 55 N (Figure 1).It spans approximately 4500 km in length and covers an area of 100 km 2 , which accounts for over half of the total Yangtze River basin area.The UYRB encompasses nine provinces, municipalities, and autonomous regions, including Yunnan, Sichuan, and Guizhou [11,30].The topography of the UYRB is characterized by significant variations, with the land elevation gradually decreasing from west to east.The regional elevation difference can reach up to 7000 m.The river flows through the Qinghai-Tibet Plateau, Yungui Plateau, eastern extension of the Qinling Mountains, and Sichuan basin.The terrain and landforms within the UYRB are diverse and complex, with the largest areas occupied by mountains and plateaus, accounting for 50% and 30% of the total area of the UYRB, respectively.Hills account for 18% of the area, while plains make up only 2% [46].The climate in the UYRB exhibits distinctive patterns on both temporal and spatial scales due to the influences of the East Asian monsoon, the South Asian monsoon, and the topography of the Qinghai-Tibet Plateau.At the temporal scale, the climate in the UYRB exhibits significant interannual variations, uneven seasonal distribution, and notable rainy and hot periods.The annual precipitation over multiple years in the UYRB ranges from 723 to 1134 mm, with an annual temperature of 8.6 to 16.8 • C. Precipitation in the basin is concentrated from April to October (Figure 2a), accounting for over 90% of the total annual precipitation.Moreover, the monthly average temperatures during this period remain above 10 • C. Throughout the year, the precipitation (temperature) in the UYRB gradually increases from January, reaches its peak in July at 162.81 mm (20.31 • C), and then gradually decreases.
At the spatial scale, the climate conditions in the UYRB also exhibit considerable variability.Specifically, the western high-altitude areas have a cold, temperate climate characterized by long winters.The majority of the remaining areas, with elevations below 3000 m, are influenced by the East Asian monsoon and the West Asian monsoon, resulting in abundant rainfall and considerable heat.Annual precipitation can reflect the moisture level of a region.Typically, areas with annual precipitation above 800 mm are classified as "humid regions", while areas with precipitation between 400 and 800 mm fall into the category of "semi-humid regions."Areas with precipitation ranging from 200 to 400 mm are classified as "semi-arid regions", and areas with precipitation below 200 mm are categorized as "arid regions."As shown in Figure 2b, the mean annual precipitation in the UYRB increases from northwest to southeast and is sequentially distributed in semi-arid, semi-humid, and humid regions.The area of humid regions accounts for about 50% of the total area of the UYRB, mainly distributed in the Wu River basin, Yibin-Yichang River basin, and downstream of the other three sub-basins.Figure 2c

Data
The observed precipitation, temperature data, and NCEP-DOE Reanalysis 2 (R-2) dataset were employed in this study to establish a statistical downscaling model for the UYRB.The CMIP6 climate data was utilized for predicting future climate.A comprehensive description of the data is provided below: The daily precipitation and temperature data were sourced from the National Meteorological Science Data Center (http://data.cma.cn/,accessed on 15 June 2023).A total of 58 meteorological stations within and around the UYRB were included, as depicted in Figure 1.The information for all stations can be found in Table A1.The observed data covers the period from 1985 to 2014 (1999-2014 for Dongxing station).Missing values were represented by −999, enabling proper identification of data sequences when applying the statistical downscaling model (SDSM).
The annual surface runoff records of Yichang hydrologic station, the outlet of the UYRB, were collected from the literature [8].Surface runoff data covers the period of 1985-2014.The surface runoff data were converted into surface runoff depth (mm) for subsequent analysis.
NCEP-DOE Reanalysis 2 (R-2) was utilized for constructing the SDSM as the input dataset [47].It should be noted that the dataset used is an improved version of its predecessor, NCEP/NCAR Reanalysis 1, as it includes updated parameterizations of physical processes and error fixes [47].NCEP-DOE Reanalysis 2 (R-2) dataset covers the same time period as the observed data, which is from 1985 to 2014.
The CMIP6 daily predictor variables from the Canadian Earth System Model version 5 (CanESM5), the Max Planck Institute for Meteorology Earth System Model version 1.2 (MPI-ESM1.2),and the Norwegian Earth System Model version 2 (NorESM2) experiments were used to project climate change in the UYRB over a medium future period (2036-2065) and long-term future period (2071-2100).Four SSP datasets were used, including SSP1-2.6,SSP2-4.5, SSP3-7.0, and SSP5-8.5.The R-2 and CMIP6 daily predictors are all grid datasets.The grid has a uniform longitudinal/latitudinal resolution of 2.8125 • [47].
All predictor variables were standardized.The 1985-2014 date range was selected as the reference period for standardization of the predictor variables [48,49].Standardized values (n) are produced from predictor values (x) utilizing the mean (µ) and standard deviation (σ) over the 1985-2014 reference period for each data source and according to individual grid boxes using the following expression:

Methods
To carry out the purpose of this paper, the following methods were used: (1) Statistical downscaling model (SDSM) was employed to project future precipitation and temperature in UYRB; (2) Mann-Kendall (MK)/modified Mann-Kendall (MMK) trend test and Theil-Sen's slope were utilized to assess the trends and the trend magnitudes; (3) extreme weather indicators recommended by ECTTDI, such as number of tropical nights (TR), maximum consecutive 5-day precipitation (R×5day), and consecutive dry days (CDD), were used to measure extreme climate events.The definitions of extreme climate indices can be seen in Table 1 [43,44]; (4) the Spearman rank correlation test (SRC) and random forest regression (RFR) model were employed to identify the key climatic factors that influence runoff; (5) multiple regression model was applied to predicting runoff; (6) some other common statistical analysis methods were used in this research.The general circulation model (GCM) is currently the most feasible method for largescale climate simulation and climate change prediction research.However, the characteristics of GCM data are large-scale and low-resolution, which makes GCM data unable to meet the needs of regional studies for high-resolution meteorological data [50].The downscaling methods of large-scale climate data usually include statistical downscaling methods (SDSM), dynamic downscaling methods, and downscaling methods combining statistics and dynamics [29,51,52].Among the three types of methods, SDSM is easy to use and understand and has a clear physical meaning.Compared with the other two, SDSM has a simple calculation process and a small calculation load.It is a common method for downscaling GCM data [53,54].In this paper, SDSM is used for downscaling and forecasting climate change.
The basic process of SDSM is shown in Figure 3. First, the statistical relationship is determined between observed data and NCEP reanalysis data, f (x).Then a statistical model is established, and the accuracy of the model is tested.Finally, the statistical model is used to downscale the GCM data to obtain future meteorological data in the UYRB.NCEP reanalysis data and GCM data used by SDSM are both large-scale predictor variables.There are 23 predictor variables in total, and these are listed in Table 2.The Mann-Kendall trend test is applied at 5% significance level to detect the significant trends in hydrometeorological time series in this study.Mann-Kendall (MK) trend test is a widely used non-parametric statistical test method.It is particularly suitable for analyzing meteorological data sequences that vary over time [56][57][58].The MK trend test is usually employed to assess the presence of a monotonic trend in a time series.A monotonically increasing/decreasing trend implies that the variable consistently increases/decreases over time, but the trend may be linear or non-linear.The advantages of MK trend test over parametric tests are as follows: (1) it is applicable to various distributions and does not require the tested variable to follow a normal distribution assumption; (2) it requires a small sample size, with a minimum data sample size of 8 to 10.
Null hypothesis, H 0 , of MK test assumes that there is no trend in the data (x 1 , x 2 , . .., x n ).Alternative hypothesis H a is that the data have a monotonic trend [59].
Test statistic S: where x j and x i are the annual data values in j-th year and i-th year, respectively, n is the length of the data.The expression for the sgn() function is as follows: For large n, the statistic S is approximately normally distributed with mean E(S) = 0 and variance Var(S) as follows: First, the statistical indicator, Z MK , is calculated: If |Z MK | > Z 1−α/2 , the null hypothesis H 0 is not acceptable, i.e., at the confidence level α, there is a significant monotonic trend in the time series data.A positive value of Z Mk indicates an increasing trend, while a negative value indicates a decreasing trend.In this study, the significance level for the MK trend test is 5%, Z 1−α/2 = 1.96.
Then, Theil-Sen's slope, β, is calculated: where Median() is the median function.Theil-Sen's slope represents the trend magnitude of a variable per unit time [26].

Modified Mann-Kendall Trend Test
Hydrometeorological time series normally display statistically significant autocorrelation, which affects the power of MK trend test to detect trends correctly in time series data.Prewhitening is a procedure for removing autocorrelation within a given time series by adding white noise series to the original series [60].But it has been reported that removal of positive/negative autocorrelation by prewhitening removes/adds a portion of trend, then reduces/increases the detection rate of significant trend in MK test [61,62].
Hamed and Rao proposed the modified Mann-Kendall (MMK) trend test to address autocorrelation issues with a variance correction approach [8,56,63].The MMK trend test is applied at 5% significance level in this study.
The variance Var*(S) in the case of autocorrelated series is as follows: where n n * S represents a correction due to the autocorrelation in the data.Var(S) is the same as the MK test.The expression of n n * S is as follows: where ρ S (i) is the autocorrelation function of the ranks of the observations.The statistical indicator, Z MK is calculated: If |Z MK | > Z 1−α/2 , the null hypothesis H 0 is not acceptable, i.e., at the confidence level α, there is a significant monotonic trend in the time series data.A positive value of Z Mk indicates an increasing trend, while a negative value indicates a decreasing trend.In this study, the significance level for the MMK trend test is 5%, Z 1−α/2 = 1.96.

Random Forest Regression Algorithm
Random forest regression model is used to identify the key climatic factors affecting surface runoff and project future surface runoff in the UYRB.Random forest is an ensemble machine learning algorithm.The random forest algorithm combines multiple weak classifiers and finally obtains the results by voting or averaging, which makes the classification or regression results of the overall model have high accuracy and generalizing performance [64].
The weak classifier used in random forest is CART decision tree, which is also called classification regression tree.When the dependent variable of the dataset is a continuous variable, the tree algorithm is a regression tree with the mean of the leaf nodes as the predicted value [65].The random forest algorithm composed of regression trees is the random forest regression (RFR) algorithm.The algorithmic process is as follows: (1) Randomly sample the training dataset with dropout and construct a sub-training dataset with the same capacity as the training dataset; (2) Use the sub-training sample set to train a CART regression tree model.In the training process, it is necessary to randomly select some features from all feature sets and then select the optimal features according to the minimum mean square error principle; (3) Repeat steps ( 1) and (2) to generate multiple CART regression trees to form a forest; (4) The mean value of the prediction results of all CART regression trees is taken as the final prediction result of random forest.
It can be seen that in the process of establishing the random forest regression model, the algorithm itself has completed the screening of important features.Therefore, in this paper, the annual average temperature, annual precipitation, CDD, TR, and R×5day in the UYRB from 1985 to 2014 are taken as independent variables, and the annual surface runoff depth of the basin is taken as dependent variable.Using Python's sklearn tool library, a random forest regression model was established to calculate the importance of the independent variables to the dependent variable and identify the key climatic factors that affect the runoff of the basin.
Using 80% of the meteorological and runoff data from 1985 to 2014 in the UYRB as the training dataset, and the remaining 20% as the testing dataset.The number of weak classifiers, n_estimators, is an important parameter of an RFR model, which has a significant impact on the fitting accuracy of the model.The model was optimized using cyclic iterative method to determine the value of n_estimators.In this study, n_estimators = 100, as it was found through pre-experiments that R 2 and Pearson's correlation coefficient were both greater than 0.90 with n_estimators = 100.All other parameters of the RFR model are default values.The established RFR model was applied, and the importance score of variables was calculated with function feature_importances_.
Before using the RFR model, we conducted a multicollinearity test on the dataset.If the variance inflation factor (VIF) value of the variable was greater than 10, it was considered that there was collinearity between the variables [66,67].In this case, data preprocessing was necessary to remove the influence of collinearity on the identification of important variables in the RFR model.
After obtaining the key climatic factors affecting surface runoff, a new RFR model was established using the key climatic factors and surface runoff in the UYRB from 1985 to 2014.The new RFR model and future climate data were applied to predict surface runoff series in the UYRB.

Multiple Regression Model
Multiple regression (MR) model typically refers to a linear regression model that involves two or more independent variables, aiming to explain the linear relationship between the dependent variable and multiple independent variables.Its mathematical model is expressed as follows: In the equation, y represents the dependent variable, E(y) denotes the estimated value of dependent variable, β 0 is a constant term, n represents the number of independent variables, x i represents the i-th independent variable, and β i represents the coefficient of the i-th independent variable.
The response of surface runoff to climate change is nonlinear, so the accuracy of applying multiple linear regression model to predict runoff is usually not as good as that of other methods, such as neural networks, random forests, support vector machines, etc.However, the prediction results of multiple linear regression model are still able to reflect the overall trend of runoff series [68,69].Meanwhile, because the multiple linear regression model is characterized by rapid calculation, simple operation, and intuitive results, multiple linear regression is still applied in the research of runoff prediction [70][71][72].
In this study, two sets of key climatic factors affecting surface runoff were identified using the Spearman rank correlation (SRC) test and random forest regression (RFR), so we needed to build a brand-new prediction model to fairly evaluate the credibility of the two sets of key climatic factors.The performance of this brand-new prediction model does not need to be particularly good; what is important is that it can reflect the trend of surface runoff.To achieve this, a multiple regression model was constructed with climatic factors in UYRB during the historical period (1985-2014) as independent variables and surface runoff depth as the dependent variable.Two multiple regression models with different input datasets were established.The multiple regression model with higher R 2 , higher adjusted R 2 , and smaller RMSE is considered to have used more reasonable input variables.These variables are the key climatic factors that influence surface runoff.

Temporal Variation Patterns of Precipitation and Temperature
Table A2 provides the Z MK values and trend magnitudes obtained from the MK/MMK trend test for annual precipitation and annual average temperature from 1985 to 2014 at the 58 meteorological stations in the UYRB.Out of the 58 meteorological stations, the annual precipitation series of 26 stations exhibited decreasing trends, with significant (at a 5% level) decreasing trends observed at 2 stations, namely Zhanyi and Huili.Additionally, the annual precipitation of 32 stations showed increasing trends, with significant increasing trends observed at the Tuotuohe, Wudaoliang, Dari, and Suoxian stations.The trend slopes of the annual precipitation series for all stations are represented using Theil-Sen's slope estimator.From 1985 to 2014, the trend slopes of annual precipitation in the UYRB ranged from −101.83 mm/decade (Huili) to 87.94 mm/decade (Dongxing).The annual precipitation of the whole UYRB was calculated using the Thiessen polygon method.Figure 4 illustrates the annual precipitation in the historical period (1985-2014) of the UYRB.The multi-year mean annual precipitation is 807.03 mm, with the maximum precipitation recorded at 923.26 mm in 1998 and the minimum at 707.62 mm in 2006.The annual precipitation in the historical period of UYRB exhibits a slightly increasing trend, with a trend slope of 10.44 mm/decade.
Regarding the annual average temperature, out of the 58 meteorological stations, except for Guiyang station, where the temperature shows a significant (at a 5% level) decrease, the temperature at the remaining 57 stations exhibited increasing trends.Among them, 53 stations showed significant (at a 5% level) warming trends.The trend slopes of the annual average temperature ranged from −0.42 • C/decade (Guiyang) to 0.99 • C/decade (Leibo).Figure 5 illustrates the annual average temperature in the historical period  in the UYRB, which showed an increasing trend.The maximum temperature recorded was 12.45 • C in 2006, which coincided with the year of the minimum precipitation.The minimum temperature occurred in 1992, with a value of 10.64 • C. The Z MK value for the annual average temperature in the UYRB is 4.14, indicating a significant increasing trend at a 5% level.The warming trend slope is 0.42 • C/decade.

Spatial Variation Patterns of Precipitation and Temperature
The spatial patterns of the trend slopes in annual precipitation  in the UYRB are shown in Figure 6.In the western high-altitude regions, there was a general increasing trend in annual precipitation, with an increase ranging from 0 to 60 mm/decade.Moreover, the closer to the source of the Jinsha River basin, the greater the magnitude of the increase in annual precipitation.Precipitation in the Jialing River basin and Mintuo River basin also increased.Annual precipitation decreased in the southwestern part of the UYRB.In particular, the lower reaches of the Jinsha River basin showed the most pronounced decrease in precipitation.The annual precipitation series of Zhanyi and Huili decreased significantly (at a 5% level), with slopes of −90.75 mm/decade and −101.83mm/decade, respectively.
The spatial patterns of the trend slopes in annual average temperature from 1985 to 2014 in the UYRB are depicted in Figure 7.It can be observed that almost the entire UYRB exhibits significant warming trends, with localized areas in the southeast showing decreasing temperature trends.Across the whole UYRB, the warming trend slope decreases from northwest to southeast, and the significance of the trend gradually diminishes.In the high-altitude regions of the northwest, the warming magnitude ranges from 0.4 to 0.8 • C/decade, which is higher than that in the eastern and southeastern parts of the UYRB.Near Guiyang, the trend of the annual average temperature shifts from an increasing trend to a decreasing trend, with a slope of −0.42 • C/decade for Guiyang.Note: "Sig."means significant at 5% level.Upward solid triangles represent significant increasing trends.Upward hollow triangles represent an increasing trend that is not significant.A downward solid triangle represents a significant downward trend.Downward hollow triangles represent a non-significant downward trend.

Identification and Verification of Statistical Downscaling Model
Figure 8 presents the simulation results of the statistical downscaling model (SDSM) during the identification period (1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004) and verification period (2005-2014).It can be observed that SDSM performs better at simulating the temperature compared to the precipitation.In each month, the simulated temperature values by SDSM closely match the observed values, capturing the dynamic characteristics of the temperature in the UYRB.However, the simulated rainfall tends to be overestimated, especially during the rainy season (from May to October).In terms of dynamic characteristics, the simulated results can also capture the variability of precipitation in the UYRB well.The simulation results of the annual precipitation and temperature by SDSM during the identification and verification periods are shown in Figure 9.At the annual time scale, the SDSM demonstrates better accuracy in simulating the temperature than the precipitation, and the simulated precipitation values exceed the observed values, which is consistent with the monthly scale simulation results (Figure 8). Figure 9 illustrates the linear trends of precipitation and temperature.It can be observed that, both in the identification and verification periods, the trends of the simulated values by the SDSM align with the observed values.In the identification period, both precipitation and temperature show an increasing trend.In the verification period, the precipitation exhibits an increasing trend, while the temperature shows a decreasing trend.This discrepancy with the MK/MMK trend test of the annual average temperature in the UYRB from 1985 to 2014 may be attributed to the fact that the verification period is too short, at only 10 years (2005-2014).Hydro-meteorological series are usually highly random over a short period, which greatly affects the accuracy of trend analysis.The trend of temperature we obtained is a shortperiod trend.From Figure 10a, it can be observed that the mean annual precipitation (MAP) for all four future SSP scenarios is higher than that in the historical period.In the mid-term (2036-2065), the MAP values of the four scenarios are close, as is shown in Figure 10a, with the four curves close together and overlapping with each other.The MAP of SSP1-2.6 in the mid-term (2036-2065) is 959.84 mm, which is an increase of 18.93% compared to the historical period .The MAP of SSP2-4.5 in the mid-term is 950.54 mm (an increase of 17.78%), SSP3-7.0 is 968.68 mm (an increase of 20.03%), and SSP5-8.5 is 965.79 mm (an increase of 19.67%) (Table 3).In the mid-term (2036-2065), SSP1-2.6 shows a decreasing trend in precipitation with a slope of −5.30 mm/decade.The other three scenarios show an increasing trend in precipitation.Under SSP3-7.0,precipitation increases significantly (at a 5% level) with a slope of 17.06 mm/decade.
In the long term (2071-2100), the MAP values of the four scenarios appear significantly different.In Figure 10a, the four curves in the long term are spread out, and the precipitation curve is more upward for a high radiative forcing value.The maximum MAP is 1085.30mm under SSP5-8.5, and the minimum is 963.92 mm under SSP1-2.6.In the long term (2071-2100), the precipitation curves of SSP2-4.5, SSP3-7.0, and SSP5-8.5 have a clear increasing trend, except for SSP1-2.6.The slope increases as the radiative forcing increases.3 show that the precipitation increases significantly (at a 5% level) under all four future scenarios.
From Figure 10b, it can be seen that under all four future scenarios, the mean annual average temperature (MAT) values in the UYRB are higher than those in the historical period.The MAT values increase with increasing radiative forcing.The temperature curves in Figure 10b can be divided into three categories.
The first category: Temperature curves continue to increase in both the mid-term (2036-2065) and the long term (2071-2100), and the slopes are relatively stable.The temperature curves of SSP3-7.0 and SSP5-8.5 belong to the first category.As can be seen in Table 4, the temperature of SSP3-7.0 increases significantly (at a 5% level) in both the mid-term and the long term.The MAT of the mid-term is 13.47 • C, which is an increase of 17.03% compared with the historical period.The long term MAT is 15.00 • C, which is an increase of 30.32% compared with the historical period.The trend slopes for the mid-term and long term under SSP3-7.0 are close to each other, with values of 0.47 and 0.45 • C/decade, respectively.The temperature of SSP5-8.5 shows a significant increasing trend (at a 5% level) both in the mid-term and the long term.The MAT values of the mid-term and long term are 13.73 • C (an increase of 19.29%) and 15.88 • C (an increase of 37.97%), respectively.The temperature slopes for the mid-term and long term under SSP5-8.5 are 0.54 • C/decade and 0.67 • C/decade.The second category: Temperature curves increase in both the mid-term (2036-2065) and long term (2071-2100), but the slope in the long term is obviously smaller than that in the mid-term.The temperature curve of SSP2-4.5 belongs to the second category.Table 4 shows that the temperature under SSP2-4.5 increases significantly (at a 5% level) in both the mid-term and the long term.The MAT of the mid-term is 13.20 • C, which is an increase of 14.68% compared with the historical period.The long term MAT is 13.97 • C, which is an increase of 21.37% compared with the historical period.The slope is 0.29 • C/decade in the mid-term and 0.10 • C/decade in the long term, which is a decrease of 65.52% compared to the mid-term.
The third category: Temperature shows an increasing trend in the mid-term (2036-2065) and a decreasing trend in the long term (2071-2100).The temperature curve of SSP1-2.6 belongs to the third category.Table 4 shows that the temperature under SSP1-2.6 increases significantly (at a 5% level) in the mid-term with a MAT value of 12.98 • C, which is an increase of 12.77% compared to the historical period.The slope in the mid-term is 0.16 • C/decade.In the long term, the temperature shows a decreasing trend.The MAT in the long term is 13.00 • C, which is an increase of 12.95% compared to the historical period.The temperature slope in the long term is −0.05 • C/decade.

Spatial Variation Patterns of Future Precipitation and Temperature
The spatial distributions of the mean annual precipitation in the future periods over the UYRB are shown in Figure 11.From northwest to southeast, the annual precipitation increases under all four scenarios.With increasing radiative forcing, the area of the humid region in UYRB expands, and the semi-humid and semi-arid regions decrease.
In the mid-term (2036-2065), the spatial patterns of precipitation under the four scenarios are the same: see the left column of Figure 11 (Figure 11a-d).Here, SSP1-2.6 is used as an example for the analysis (Figure 11a), and the analysis of the other three scenarios is omitted.Compared with the historical period (1985-2014), the upstream region of the Mintuo River basin changes from a semi-humid region to a humid region, which increases the total area of the humid region in the UYRB.Within the humid region, the area of the heavy precipitation region with precipitation greater than 1200 mm increases.Almost all of the Yibin-Yichang River basin becomes a heavy precipitation region, as do the lower Jialing River basin, the lower Wu River basin, and the lower Jinsha River basin.Semi-arid regions are distributed in the upstream region of the Jinsha River basin, and the area of semi-arid regions in the mid-term is smaller than that in the historical period.
In the long term (2071-2100), the total area of the humid regions in the UYRB is greater than that in the mid-term (2036-2065), as shown in the right column of Figure 11 (Figure 11e-h).The spatial distribution patterns of the mean annual precipitation are similar under the SSP3-7.0 and SSP5-8.5 scenarios.The area of humid regions under the high radiative forcing scenario is larger than that under the low radiative forcing scenario.The middle reaches of the Jinsha River basin change from a semi-humid region to a humid region, and the area of the heavy precipitation region increases with the increasing radiative forcing.In contrast to SSP1-2.6 and SSP2-4.5 scenarios (Figure 11e,f), under SSP3-7.0 and SSP5-8.5 (Figure 11g,h), almost the whole Wu River basin turns into a heavy precipitation region.
Figure 12 shows the spatial distribution of the trend slopes for the future annual precipitation of the UYRB.The results of the MK/MMK trend test and trend slopes can be found in Table A3.In the mid-term (2036-2065), the spatial pattern of precipitation slopes under SSP1-2.6 is different from the other three scenarios, as shown in the left column of Figure 12 (Figure 12a-d).Under SSP1-2.6,precipitation decreases in the lower Jinsha River basin, the lower Mintuo River basin, the Jialing River basin, and the Yibin-Yichang River basin.Precipitation increases in the other regions of the UYRB.The decreasing regions are indicated in yellow and orange, and the increasing regions are blue in Figure 12a.Under SSP2-4.5, SSP3-7.0, and SSP5-8.5 (Figure 12b-d), precipitation increases significantly (at 5e% level) in most parts of the UYRB, including the upper and lower parts of the Jinsha River basin, the Yibin-Yichang River basin, and the Wu River basin.In the Jialing River basin and the upper Mintuo River basin, precipitation increases under SSP3-7.0 and SSP5-8.5 (Figure 12c,d) and decreases under SSP2-4.5 (Figure 12b).
In the long term (2071-2100), precipitation increases in most parts of the UYRB under all four scenarios, as shown in the right column of Figure 12 (Figure 12e-h).Under SSP1-2.6,precipitation decreases only in the upper Jinsha River basin, the Wu River basin, and the Yibin-Yichang River basin and increases in the remaining regions.The spatial pattern of the precipitation slope under SSP1-2.6 in the long term is the opposite of that in the mid-term.Under SSP2-4.5, precipitation decreases in the middle reaches of the Jinsha River basin and the middle reaches of the Wu River basin and increases in the remainder of the UYRB.Under SSP3-7.0 and SSP5-8.5, precipitation increases significantly (at a 5% level) over the vast majority of the UYRB.The regions with decreasing precipitation are small in size, and they are only distributed within the lower Jinsha River basin, the lower Mintuo River basin, and the Wu River basin.
The spatial distribution of the mean annual average temperature in the future over the UYRB is depicted in Figure 13.Under the four scenarios, the mean annual average temperature in the UYRB increases along the northwest to southeast direction, and this spatial pattern is shown in Figure 13 as a change in color from blue to red.In the mid-term (2036-2065), the spatial patterns of temperature under the four scenarios are the same, as seen in the left column of Figure 13 (Figure 13a-d).The area of the cold regions with temperatures less than 0 • C (indicated in blue color) in the northwestern part of the UYRB shrinks slightly as the radiative forcing increases.The area of the hightemperature regions with temperatures greater than 20 • C increases slightly but remains mainly within the Yibin-Yichang River basin.
In the long term (2071-2100), the area of the cold regions decreases and the area of the high-temperature regions increases dramatically as the radiative forcing increases.The area of high-temperature regions under SSP5-8.5 is the largest and is distributed in the Mintuo River basin, Jialing River basin, Yibin-Yichang River basin, and Wu River basin.Under SSP1-2.6, the spatial pattern of temperature in the long term is almost the same as in the mid-term (Figure 13a,e).Under SSP2-4.5, SSP3-7.0, and SSP5-8.5 (Figure 13f-h), the area of the cold regions in the long term is smaller than that in the mid-term, and the area of the high-temperature regions is larger.
The trend slopes of annual average temperature in the future for the four scenarios in the UYRB can be seen in Figure 14 and Table A4.In different future periods, temperature over the entire UYRB under all four scenarios increases significantly (at a 5% level), except for SSP1-2.6 in the long term (Figure 14e).Under SSP2-4.5, the slope in the long term is smaller than that in the mid-term.In the mid-term, the temperature increases significantly (at a 5% level) in the whole UYRB under SSP2-4.5 (Figure 14b), while in the long term (Figure 14f), the increasing temperature trends are not significant in the Mintuo River basin, Jialing River basin, or Yibin-Yichang River basin.Under SSP3-7.0, the slopes in the Jialing River basin, Mintuo River basin, and Yibin-Yichang River basin in the long term (Figure 14g) are smaller than those in the mid-term (Figure 14c).Under SSP5-8.5, the slope in the long term (Figure 14h) is greater than that in the mid-term (Figure 14d).

Key Climatic Factors Affecting Surface Runoff
In the historical period (1985-2014), the annual surface runoff depth in UYRB showed a decreasing trend (Figure 15).The trend test yielded a result of Z MK = −1.41,indicating a non-significant decreasing trend with a slope of −14.43 mm/decade.In the historical period, the mean annual surface runoff (MAR) depth in the UYRB was 418.62 mm, with a maximum value of 517.63 mm observed in 1998, which coincided with a year of heavy rainfall and resulted in a severe flood disaster in the Yangtze River basin.The minimum value of 284.60 mm was recorded in 2006, which was a year with the lowest precipitation and highest temperature.This indicates that extreme weather events may have an impact on the surface runoff in the UYRB.
We considered annual precipitation, annual average temperature, number of tropical nights (TR), maximum consecutive 5-day precipitation (R×5day), and consecutive dry days (CDD) as the climatic factors to be identified, using the Spearman rank correlation (SRC) test and the random forest regression (RFR) model to identify the key climatic factors influencing surface runoff depth.The Spearman rank correlation test is applied at a 5% level of significance, and it assesses the strength of the relationship between different factors and surface runoff by examining the magnitude of the correlation coefficients.The random forest regression model calculates the importance score of each climatic factor in relation to surface runoff.The sum of the importance scores of all factors is 1, with higher importance score values indicating greater significance of the factor in influencing surface runoff.Before the climatic factor datasets were input into the RFR model, the datasets were analyzed for multicollinearity.If the variance inflation factor (VIF) value of the variable was greater than 10, it was considered that there was collinearity between the variables.The results of the correlation test and importance analysis can be found in Table 5.The results of the SRC test indicate a strong positive correlation between precipitation and surface runoff, with the correlation coefficient being the highest at 0.75.The R×5day shows a moderate positive correlation with surface runoff (correlation coefficient of 0.46), while the CDD exhibits a very weak positive correlation with surface runoff (correlation coefficient of 0.15).The mean annual average temperature and TR are negatively correlated with surface runoff, with correlation coefficients of −0.31 (weak correlation) and −0.44 (moderate correlation), respectively.At a significance level of 5%, the relationships between annual precipitation, R×5day, TR, and surface runoff are significant and are the key climatic factors influencing surface runoff.
The VIF values for all variables are less than 10, indicating that there is no multicollinearity affecting the importance analysis results of the RFR model.The importance analysis results reveal that precipitation has the greatest importance in explaining the variation in surface runoff, with an importance score of 0.57.In descending order of importance, the factors are ranked as follows: annual precipitation > R×5day > annual average temperature > CDD > TR.The combined importance of precipitation, R×5day, and temperature reaches 0.85, indicating a high explanatory power for surface runoff.On the other hand, CDD and TR have a minor impact on surface runoff.
The two methods, SRC and RFR, identify slightly different key climatic factors that influence surface runoff.Both the SRC and RFR results confirm that precipitation and R×5day are crucial factors affecting surface runoff.The SRC identifies the third key factor as TR, whereas the RFR identifies annual temperature as the third key factor.These variations in results may arise from differences in the underlying algorithms and statistical approaches used by the two methods.Nonetheless, the consensus on the importance of precipitation and R×5day highlights their significant role in determining surface runoff patterns.

Changes in Future Surface Runoff
Since there are differences in the conclusions drawn from the Spearman rank correlation (SRC) test and the random forest regression (RFR) model, two MR models, MR SRC and MR RFR, were constructed, and the model with the higher simulation accuracy was considered to have used more reasonable input variables.The two MR models' parameters and indicators can be found in Table 6, and the simulation results are illustrated in Figure 16.Table 6 shows that the R 2 and adjusted R 2 values of MR RFR are higher than those of MR SRC , and the RMSE is smaller.Figure 16  The predictions of surface runoff in the UYRB under the four future scenarios are shown in Figure 17.Table 7 provides the statistical characteristics and trends of surface runoff in the future.It can be seen that the mean annual surface runoff (MAR) values of the UYRB in the mid-term and long term are greater than those in the historical period, which indicates an increase in future surface runoff in the UYRB.
In the mid-term (2036-2065), SSP1-2.6 has the largest MAR of 467.22 mm, which is an increase of 12.09% compared with the historical period, followed by SSP3-7.0 at 453.75 mm (an increase of 8.86%), SSP2-4.5 at 450.79 mm (an increase of 8.15%), and SSP5-8.5 at 440.87 mm (an increase of 5.77%).Surface runoff tends to decrease under all four future scenarios.Surface runoff decreases significantly (at a 5% level) under SSP1-2.6,with a slope of -9.46 mm/decade.SSP5-8.5 has the largest slope of −11.06 mm/decade.
In the long term (2071-2100), the MAR values of the four future scenarios do not differ significantly from those in the mid-term (2036-2065).The MAR is the largest under SSP1-2.6 with 469.27 mm.Surface runoff under all four scenarios shows an increasing trend, which is the opposite of the mid-term.Surface runoff increased significantly (at a 5% level) under SSP1-2.6,with a trend slope of 7.99 mm/decade.SSP2-4.5 has the largest trend slope of 9.15 mm/decade.

Discussion
During the historical period (1985-2014), the climate in the UYRB exhibited a trend towards warming and humidification.Both annual precipitation and annual average temperature increased over time, with the warming and humidification trend being more pronounced in the northwest and central regions compared to the other areas.The increase in precipitation across the entire basin was not significant.Spatially, the areas showing an increasing trend in precipitation were mainly located in the northwest and in the Jialing River basin.The closer to the source of the Yangtze River, the greater the slope increase in annual precipitation.The temperature series in the whole UYRB showed a significant increasing trend.Spatially, the slope of temperature decreased from northwest to southeast, with the source area and central region of the UYRB experiencing a larger increase in temperature.
In the future period, the MAP and MAT values of the UYRB were shown to be greater than those in the historical period.The climate of the UYRB will experience warming and humidification in the future.The trend slopes of precipitation and temperature increase with increasing radiative forcing.These findings are consistent with the conclusions of the study conducted by Li  Our study reveals that the temporal and spatial variation patterns of precipitation and temperature under SSP1-2.6 in the UYRB are different from the other three scenarios.The precipitation series under SSP1-2.6 first decreases and then increases from the mid-term (2036-2065) to the long term (2071-2100), while the temperature under SSP1-2.6 first increases and then decreases.Precipitation and temperature under the other three scenarios show increasing trends in both the mid-term (2036-2065) and long term (2071-2100).The reason for this is that SSP1-2.6 represents the sustainable development pathway with a low radiative forcing level.The SSP1-2.6 scenario involves sustainable development practices and proactive measures to mitigate climate change, resulting in lower fossil fuel consumption and greenhouse gas emissions compared to the other scenarios [74,75].Another low radiative forcing scenario, SSP2-4.5, shows an overall increasing trend in precipitation and temperature.Comparatively, the high radiative forcing scenarios, SSP3-7.0 and SSP5-8.5, characterized by intensive social and economic development patterns with substantial fossil fuel consumption and greenhouse gas emissions, lead to significant increases in precipitation and temperature across all future periods.At the spatial scale, the mean annual precipitation and mean annual temperature in the UYRB under the four future scenarios exhibit similar spatial distribution patterns to the historical period.The upstream region generally experiences lower precipitation and temperatures compared to the downstream region, primarily due to variations in elevation and topography.At the same time, the warming and humidification of the climate are more pronounced in the high-altitude regions.
In the process of identifying key climatic factors affecting surface runoff, the results of the random forest regression model are more reasonable than those of the Spearman rank correlation test.This may be because the relationship between climatic factors and surface runoff is very complex, and traditional statistical analysis methods have weaker recognition and processing capabilities than machine learning algorithms for this complex relationship.According to the results of the random forest regression model, the surface runoff in the UYRB is mainly affected by annual precipitation, R×5day, and annual average temperature.The sum of the three factors' interpretation degree of surface runoff reaches 85%.Annual precipitation is the most important influencing factor of surface runoff, followed by R×5day and temperature.From the parameters of the multiple linear regression model, it can be seen that precipitation and R×5day are positively correlated with surface runoff, while temperature is negatively correlated with surface runoff.This is because precipitation is an important source of surface runoff, and changes in temperature affect evaporative conditions in the region.An increase in air temperature leads to stronger evaporation, which in turn increases the amount of surface runoff loss [76].
During the historical period (1985-2014), the surface runoff in the UYRB exhibited a decreasing trend [8,30].The MAR values of the UYRB were higher than the historical period under all four scenarios, and the streamflow in the UYRB increased [11,46].Compared with the historical period, the SSP1-2.6 scenario showed the highest increase in surface runoff, with an increase of 12.09% in the mid-term (2036-2065) and 12.58% in the long term (2071-2100).Surface runoff under all four scenarios showed a decreasing trend in the mid-term (2036-2065) and an increasing trend in the long term (2071-2100).Among the four scenarios, only the surface runoff trends under SSP1-2.6 were significant (at a 5% level) in both the mid-term and long term.This is due to the unique temporal and spatial variation patterns of climate under SSP1-2.6.In the mid-term (2036-2065), decreasing precipitation and increasing temperatures under SSP1-2.6 resulted in a significant (at a 5% level) decrease in surface runoff due to decreased recharge and increased loss of surface runoff.Similarly, when precipitation increased and temperature decreased in the long term (2071-2100), surface runoff increased significantly (at a 5% level).
The trends of surface runoff are not significant under SSP2-4.5,SSP3-7.0, or SSP5-8.5.This is due to the fact that precipitation and temperature act in opposite directions on surface runoff under these three scenarios.Under SSP2-4.5, SSP3-7.0, and SSP5-8.5, both precipitation and temperature increase in all future periods.However, the trends of surface runoff in two future periods are different.This indicates that the intensity of the influence of precipitation and temperature on surface runoff will change in different future periods.In the mid-term (2036-2065), the influence of temperature on surface runoff dominates, causing surface runoff to decrease.In the long term (2071-2100), precipitation has a stronger influence on surface runoff than temperature.

Limitations of This Study
It should be acknowledged that this study contains some limitations.First, in this study, the mean of a model ensemble is used as the future climate projection in the UYRB [77].However, no further analysis was conducted on the uncertainty of the predicted results.This may lead to an incomplete analysis when discussing the impact of climate change on surface runoff [1].
Second, no physically based hydrologic models (e.g., SWAT) were used to predict future surface runoff in this study, but rather a data-driven climate-runoff prediction model was utilized.The reason for this is that building a reliable hydrologic model usually requires a large amount of observational data, including daily precipitation, daily maximum temperature, daily minimum temperature, runoff records, and some other spatial data [78].In this study, the data we collected were daily precipitation, daily mean temperature, and annual runoff at the watershed outlet, and it is difficult to build a reliable hydrological model using these data, so a data-driven climate-runoff prediction model is a better choice.It must also be acknowledged that the climate-runoff prediction model developed in this paper has limited accuracy and can only be used to analyze surface runoff trends.This is because the model only considers climate change, and other factors that affect surface runoff, such as land use and human activities, are not considered.

Conclusions
This study analyzed the temporal and spatial patterns of precipitation and temperature in the upper Yangtze River basin (UYRB) during the historical period (1985-2014) and two future periods: mid-term (2036-2065) and long term (2071-2100).The key climatic factors influencing surface runoff were identified, and the surface runoff sequence under climate change conditions was predicted.The main conclusions are as follows: (1) During the historical period (1985-2014), the UYRB experienced warming and humidification trends, with both annual precipitation and temperature increasing over time.
The warming and humidification trends were more pronounced in the northwestern and central regions.
(2) Under SSP1-2.6,SSP2-4.5, SSP3-7.0, and SSP5-8.5, relative to those in the historical period , the mean annual precipitation (MAP) and mean annual average temperature (MAT) over the UYRB increase in the future.The MAP and MAT values increase more in the long term (2071-2100) than in the mid-term (2036-2065).
(3) Under SSP2-4.5, SSP3-7.0, and SSP5-8.5, precipitation and temperature show increasing trends in both the mid-term (2036-2065) and long term (2071-2100).Under SSP1-2.6, the precipitation series first decreases and then increases from the mid-term (2036-2065) to the long term (2071-2100), while the temperature first increases and then decreases.The trend slopes of precipitation and temperature are minimized under SSP1-2.6.This indicates that the sustainable development pathway is conducive to reducing the impact of climate change in the UYRB.(4) At the spatial scale, the humid regions (annual precipitation > 800 mm) and hightemperature regions (annual temperature > 20
illustrates that the mean annual average temperature also increases from northwest to southeast.Cold regions with mean annual temperatures below 0 • C are concentrated in the high-altitude areas in the northwest, with elevations exceeding 4000 m.The mean annual temperature in the central and eastern parts of the UYRB ranges from 14 to 20 • C.

Figure 2 .
Figure 2. The temporal and spatial patterns of climate in the UYRB (1985-2014), (a) monthly precipitation and temperature, (b) mean annual precipitation, (c) mean annual average temperature.

Figure 3 .
Figure 3. Flow chart of basic ideas of statistical downscaling method (SDSM).
SDSM 4.2 software is used in this paper [55].In the process of establishing the statistical model, 1985-2004 is the identification period of the model, and 2005-2014 is the verification period.The established SDSM model is applied to downscale the large-scale data of CanESM5, MPI-ESM1.2,and NorESM2.The mean value of the model ensemble is used as the future climate projection in the UYRB.

Figure 6 .
Figure 6.Spatial distribution of annual precipitation trend slope in the UYRB (1985-2014).Note: "Sig."means significant at 5% level.Upward solid triangles represent significant increasing trends.Upward hollow triangles represent an increasing trend that is not significant.A downward solid triangle represents a significant downward trend.Downward hollow triangles represent a non-significant downward trend.

Figure 7 .
Figure 7. Spatial distribution of annual average temperature trend slope in the UYRB (1985-2014).Note: "Sig."means significant at 5% level.Upward solid triangles represent significant increasing trends.Upward hollow triangles represent an increasing trend that is not significant.A downward solid triangle represents a significant downward trend.Downward hollow triangles represent a non-significant downward trend.

Figure 8 .
Figure 8. Simulation results of SDSM during identification (a) and verification (b) periods.Note: Obs_Pre is the observed precipitation, Sim_Pre is the simulated precipitation, Obs_Tem is the observed temperature, and Sim_Tem is the simulated temperature.

3. 3 .
Figure 10 represents the annual precipitation and annual average temperature series for the historical (1985-2014) and future (mid-term: 2036-2065, long term: 2071-2100) periods in the UYRB.The results of the MK trend test and trend slopes for future precipitation are provided in Table3.

Figure 9 .
Figure 9. Simulation results and linear trends of annual precipitation and temperature over the UYRB during 1985-2014.Note: Obs_Pre is the observed precipitation, Sim_Pre is the simulated precipitation, Obs_Tem is the observed temperature, and Sim_Tem is the simulated temperature.

Figure 10 .
Figure 10.Prediction results of (a) future annual precipitation and (b) annual average temperature for the whole UYRB.

Figure 16 .
Figure 16.Comparison of prediction results between MR SRC and MR RFR models.
also demonstrates that the predicted values of MR RFR are closer to the 1:1 line.This means that the input variables of MR RFR are more reasonable.The key climatic factors affecting surface runoff are annual precipitation, R×5day, and annual average temperature.The key climatic factors and surface runoff data in the UYRB from 1985 to 2014 were used to build a new RFR model.The new RFR model was utilized to predict future surface runoff by inputting climate data from the four future scenarios.

Table 2 .
List of the 23 predictor variable IDs and corresponding names.

Table 3 .
Projected mean annual precipitation (MAP) change in different future periods in the UYRB.

Table 4 .
Projected mean annual temperature (MAT) change in different future periods in the UYRB.

Table 5 .
Results of correlation test and importance analysis.
Note: The "Bold" represents a significant result at the 5% level.

Table 6 .
Parameters and indicators of two multiple regression models.

Table 7 .
Projected mean annual surface runoff (MAR) change in different future periods in the UYRB.
Note:The "Bold" represents a significant result at the 5% level.The "italics" represent the MMK values of time series with autocorrelation.

Table A1 .
(5)) in the future are projected to increase.Conversely, the semi-arid regions (annual precipitation between 200 mm and 400 mm) and cold regions (annual temperature below 0 • C) are expected to decrease.These changes are enhanced with increasing values of radiative forcing.The slopes of precipitation and temperature are greater in high-altitude regions than those in low-altitude regions.(5)Theintroduction of the random forest regression algorithm can effectively improve the credibility of the identification of key climatic factors influencing surface runoff.Apart from annual precipitation and annual average temperature, the key climatic factors influencing surface runoff also include the extreme climate indicator, R×5day.Precip-itation and R×5day are positively correlated with surface runoff, and temperature is negatively correlated with surface runoff.(6) Under SSP1-2.6,SSP2-4.5, SSP3-7.0, and SSP5-8.5, the mean annual surface runoff (MAR) values in the UYRB are greater than those in the historical period in both the mid-term (2036-2065) and the long term (2071-2100).Surface runoff tends to decrease in the mid-term (2036-2065) due to temperature increases and to increase in the long term (2071-2100) due to precipitation increases.Under SSP1-2.6, the surface runoff increases significantly (at a 5% level), which means the risk of flood disaster in the future in the UYRB may increase under the scenario of lower radiative forcing.Cont.

Table A2 .
The MK/MMK test (at 5% level) results and trend slope of annual mean precipitation and temperature over the UYRB during 1985-2014.The "*" indicates passing the 5% significance level.The italics represent the MMK values of time series with autocorrelation.Positive/negative Z MK values indicate increasing/decreasing trend.

Table A3 .
The MK/MMK test (at 5% level) results and trend slope (mm/decade) of annual mean precipitation over the UYRB during 2036-2065.The bold indicates passing the 5% significance level.The italics represent the MMK values of time series with autocorrelation.Positive/negative Z MK values indicate increasing/decreasing trend.