Combined Effects of Land Use/Cover Change and Climate Change on Runoff in the Jinghe River Basin, China

: In the context of global warming and intensiﬁed human activities, the quantitative assessment of the combined effects of land use/cover change (LUCC) and climate change on the hydrological cycle is crucial. This study was based on the simulation results of future climate and LUCC in the Jinghe River Basin (JRB) using the GFDL–ESM2M and CA–Markov combined with the SWAT models to simulate the runoff changes under different scenarios. The results revealed that the future annual precipitation and average temperature in the JRB are on the increase, and the future LUCC changes are mainly reﬂected in the increase in forest and urban lands and decrease in farmlands. Changes in runoff in the JRB are dominated by precipitation, and the frequency of extreme events increases with the increase in the concentration of CO 2 emissions. Under four climate scenarios, the contribution of future climate change to runoff changes in the JRB is − 8.06%, − 27.30%, − 8.12%, and +1.10%, respectively, whereas the inﬂuence of future LUCC changes is smaller, ranging from 1.14–1.64%. In response to the future risk of increasing water-resources stress in the JRB, the results of this study can provide a scientiﬁc basis for ecological protection and water-resources management and development.


Introduction
The mechanisms of runoff generation and confluence in river basins have changed, thereby leading to changes in river-runoff dynamics [1,2].The fifth assessment report of the Intergovernmental Panel on Climate Change indicated that the global mean temperature will increase by about 1.1 • C and may reach or exceed 1.5 • C in the next 20 years [3].This temperature increase is expected to enhance the water-vapor-holding capacity of the atmospheric boundary layer, resulting in increased intensity and frequency of extreme precipitation events, especially in the middle-low latitudes [4,5].Simultaneously, changes in meteorological factors, such as temperature and precipitation, will also alter the hydrological cycle in river basins globally [6].In recent years, urbanization, land desertification, and deforestation have intensified under the influence of human activities, resulting in global land use/cover change (LUCC), which in turn alters the underlying surface conditions [7][8][9].Changes in the underlying surface conditions can affect the hydrological cycle in river basins by altering soil evaporation, canopy interception, runoff frequency and intensity, and eventually, river dynamics [10].Future climate and LUCC-especially in the next few decades-will have a substantial impact on hydrological changes in most river basins globally, with the combined effects of the two on runoff being notably different in distinct river basins [11,12].Therefore, simulating and predicting the hydrological processes in river basins under different future climate and LUCC scenarios have great theoretical and practical importance for elucidating the drivers of future runoff changes as well as realizing scientific management and rational planning of water resources.
Currently, the projections for future LUCC typically have generalized hypotheses about future conversions and are modeled based on socioeconomic and geographical driving factors [13,14].Commonly used models to simulate LUCC include CA-Markov, CLUE-S, Markov-FLUS and other coupled models [15,16].Among them, the CA-Markov model is widely used in regional LUCC simulation and prediction [17].Chen et al. [18] predicted the future LUCC in the Luanhe River Basin based on the CA-Markov model, and used the SWAT model to assess the combined impacts of future climate and LUCC on hydrological drought in conjunction with the results of RCP2.6, RCP4.5, and RCP8.5 future climate scenarios.Daba and You [19] projected the 2038 LUCC along the upper stream of the Awash River based on the CA-Markov model and made recommendations for sustainable water-resource planning and management in the basin.Li et al. [20] quantitatively predicted the LUCC at different periods using the CA-Markov model and investigated the impact of LUCC on runoff and evapotranspiration in the upper reach of the Han River in China.Overall, LUCC primarily affects hydrology and water resources by influencing the hydrological cycle process in river basins, and the variability between different LUCC scenarios will affect the prediction of the hydrological response [21,22].Therefore, considering the future, LUCC is essential to improve the simulation accuracy of future hydrological processes.
Hydrological processes typically involve complex nonlinear, abrupt, and random processes [23,24].The quantitative identification of the runoff response to future climate and LUCC has become a pressing issue in relation to watershed hydrological cycles [25].Increasing evidence has proposed a watershed test, statistical analysis, and model simulation methods [26][27][28].Among these, hydrological model simulation has been widely used because it comprehensively considers the spatial heterogeneity of river basins and physical processes [11,29].The calibrated hydrological model can have different climate and LUCC scenarios input to assess the hydrological response patterns in a changing environment [30,31].However, previous studies have focused more on the impacts of future climate change on hydrological processes and less on the impacts of future LUCC, which will increase the uncertainty of runoff simulation results [32,33].
In this study, a semi-distributed model was used to simulate the hydrological processes in a river basin in China by inputting meteorological and LUCC data describing different future scenarios.Future runoff changes were predicted, and driving factors were quantitatively identified.The Jinghe River Basin (JRB) was selected specifically because of the large number of ecological restoration projects, such as the Grain for Green Project and Protection of Forest System in the Middle Yellow River [34,35], and the pronounced impact of the LUCC on the hydrological cycle.We analyzed the characteristics of future climate and LUCC under different future climate scenarios and periods in the study area.Based on this analysis, we predicted the runoff change under the combined action of the two abovementioned factors and quantitatively identified their contributions to the runoff change.We make recommendations that can serve as a scientific basis for watershed management in the context of climate change and LUCC.

Study Area
The Jing River emanates from the foothills of Liupan Mountain and flows through Shaanxi, Gansu, and Ningxia Provinces into the Wei River in Chenjiatan, Xi'an City, Shaanxi Province.With a total length of 455.1 km, the Jing River is the largest tributary of the Wei River and a secondary tributary of the Yellow River (Figure 1a,b).The JRB is located in a high-intensity soil and water erosion area of the Loess Plateau (34   1c).In the JRB, areas with slopes of less than 30 • account for 91.13% of the total drainage area, whereas areas with slopes between 15 • and 30 • account for 37.35% of the total drainage area (Figure 1d).The JRB is characterized by different area ratios of slope-aspect grades (flat ground (2.54%) < ubac (10.24%) < adret (11.09%) < semi-adret (37.46%) < semi-ubac (38.68%); Figure 1e) and dominated by primitive soil, accounting for approximately 79.95% of the total drainage area (Figure 1f).The JRB has a temperate continental climate, with a rainy and hot period as well as distinct seasons; the annual mean precipitation and temperature are 517.9mm and 8 • C, respectively [36].The JRB is dominated by farm-and grasslands, with the sum of the two areas accounting for more than 85% of the total area of the basin [37].There are no large reservoirs in the JRB; therefore, the influence of reservoirs on natural runoff is limited.
Atmosphere 2023, 14, 1237 3 of 20 located in a high-intensity soil and water erosion area of the Loess Plateau (34°46'-37°19' N, 106°14'-108°42' E; 356-2919 m elevation; 45,421 km 2 drainage area; Figure 1c).In the JRB, areas with slopes of less than 30° account for 91.13% of the total drainage area, whereas areas with slopes between 15° and 30° account for 37.35% of the total drainage area (Figure 1d).The JRB is characterized by different area ratios of slope-aspect grades (flat ground (2.54%) < ubac (10.24%) < adret (11.09%) < semi-adret (37.46%) < semi-ubac (38.68%); Figure 1e) and dominated by primitive soil, accounting for approximately 79.95% of the total drainage area (Figure 1f).The JRB has a temperate continental climate, with a rainy and hot period as well as distinct seasons; the annual mean precipitation and temperature are 517.9mm and 8 °C, respectively [36].The JRB is dominated by farm-and grasslands, with the sum of the two areas accounting for more than 85% of the total area of the basin [37].There are no large reservoirs in the JRB; therefore, the influence of reservoirs on natural runoff is limited.

Data Sources
Daily measurements from 14 national meteorological stations (China Meteorological Data Service Centre, http://data.cma.cn,accessed on 15 September 2022) spanning from 1951 to 2019 in and around the JRB (Figure 1c) were used.This study was based on the standard-sequence method proposed by DeGaetano et al. [25], with the sites with the best correlations being used to interpolate missing data.On this basis, the inverse-distanceweighting (IDW) method was used to obtain the spatial change characteristics of each meteorological indicator in the JRB.In addition, the Penman-Monteith formula, which is recommended by the World Food and Agriculture Organization, was used for calculating the potential evapotranspiration at each site [38].Monthly runoff data, also spanning from 1951 to 2019, came from the Zhangjiashan Hydrological Station (that is, the control station of the JRB) and was primarily sourced from the Yellow River Conservancy Commission.
This study used the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios of the geophysical fluid dynamics laboratory Earth-system model with modular ocean model

Data Sources
Daily measurements from 14 national meteorological stations (China Meteorological Data Service Centre, http://data.cma.cn,accessed on 15 September 2022) spanning from 1951 to 2019 in and around the JRB (Figure 1c) were used.This study was based on the standard-sequence method proposed by DeGaetano et al. [25], with the sites with the best correlations being used to interpolate missing data.On this basis, the inverse-distanceweighting (IDW) method was used to obtain the spatial change characteristics of each meteorological indicator in the JRB.In addition, the Penman-Monteith formula, which is recommended by the World Food and Agriculture Organization, was used for calculating the potential evapotranspiration at each site [38].Monthly runoff data, also spanning from 1951 to 2019, came from the Zhangjiashan Hydrological Station (that is, the control station of the JRB) and was primarily sourced from the Yellow River Conservancy Commission.
This study used the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios of the geophysical fluid dynamics laboratory Earth-system model with modular ocean model (GFDL-ESM2M) version 4, proposed by the National Oceanic and Atmospheric Administration.Datasets, obtained from the National Tibetan Plateau Data Center (http: //data.tpdc.ac.cn, accessed on 3 October 2022).Previous studies have revealed that the widely used GFDL-ESM2M climate model is better suited to a variety of climatic conditions in different regions.In addition, the adaptability of the precipitation and temperature datasets contained in the GFDL-ESM2M climate model in China was tested and verified to better reflect the future climate-change characteristics of China [39].
The 30 × 30 m resolution digital elevation model (DEM) of the JRB was derived from the Geospatial Data Cloud (http://www.gscloud.cn,accessed on 5 January 2023).The 1 × 1 km resolution LUCC and soil datasets for each period were obtained from the Resources and Environment Science and Data Center of the Chinese Academy of Sciences (http: //www.resdc.cn,accessed on 20 November 2022).Here, we used Krasovsky_1940_Albersprojected coordinates.

Cellular Automata-Markov Model
The Cellular Automata-Markov (CA-Markov) model can use local simple transformations to simulate and predict the spatial structure of complex systems and can compensate for the shortcomings of Markov in spatial prediction [40].The CA-Markov model, through the combination of the two techniques, can simulate the spatiotemporal variation characteristics of complex geographic systems [19].In this study, the CA-Markov model was used to predict future LUCC.Initially, a suitability atlas of LUCC transfer under the influence of various factors was established and was subsequently combined with the area-transfer matrix output using the Markov model, with future changes predicted based on the LUCC data of the baseline year 2018 [41].The CA model is formulated as follows: where S i is the cellular state at the i-th moment, N is the cellular neighborhood, and f is the cellular transition rule.In this study, the process of LUCC type is regarded as a Markov process, and the calculation formula of transition probability matrix is as follows [42]: where P i,j is the probability that state i transfers to state j, and satisfies that P i,j > 0 and ∑ n i=1 P i,j = 1.
The LUCC status at each moment can be calculated using the following Equation (3): The Kappa coefficient (K a ) examines the consistency between the simulated image results and measured image data in their entirety and is widely used in accuracy evaluations of LUCC simulations and remote sensing image interpretations [43]: where P 0 is the comparison probability between the predicted and measured images, and P c is the expected value of the correct ratio of the predicted image.According to the results of K a , the prediction accuracy can be divided into the following five levels: almost inconsistent [0.0, 0.2), low consistency [0.2, 0.4), moderately consistent [0.4,0.6), highly consistent [0.6, 0.8), and almost completely consistent [0.8, 1.0].

SWAT Model
In this study, the SWAT model was used in simulating the runoff change process in the JRB from 1951 to 1996, with the years 1951-1960, 1961-1980, and 1981-1996 being used as the warm-up, calibration, and verification periods, respectively, and the SUFI-2 algorithm of SWAT-CUP was used for parameter optimization [44].In addition, the correlation coefficient (R 2 ), Nash-Sutcliffe efficiency coefficient (NS), and relative error (RE) were used to verify the accuracy of the simulation results and evaluate the applicability of the SWAT model in the JRB [45,46]: where Q and S are the observation and simulated values of annual-scale runoff, respectively, in m 3 /s; Q and S are the means of the observation and simulated values, respectively, in m 3 /s; m is the number of simulated fields.

Scenario Settings
To analyze the law governing the responses of future runoff changes to climate and LUCC in the JRB, hydrometeorological data from 1951 to 2019 and LUCC in 2018 were used as the baseline scenario (A0).The following two set-ups for future climate and LUCC combinations were established: A1: Set-up considering future climate change alone, that is, the combination of meteorological data from 2020 to 2050 under four future scenarios and LUCC in 2018 (historical).
A2: Set-up considering the combined effect of future climate and LUCC, that is, future runoff changes under the combinations of 2020-2029 meteorological data and LUCC in 2025, 2030-2039 meteorological data and LUCC in 2035, and 2040-2050 meteorological data and LUCC in 2045.
Through a comparative analysis of the runoff simulation results under different scenarios, the influences of climate and LUCC on runoff change can be quantitatively analyzed (Figure 2).

Historical Climate Change Characteristics
Based on meteorological data from 1951 to 2019 from the JRB, this study used the IDW interpolation method to obtain the spatial-variation characteristics of each meteorological element in the river basin (Figure 3).The multi-annual mean precipitation in the JRB lies between 339.6 and 590.7 mm, with its spatial variation being similar to that of the multi-annual mean relative humidity (52.7-69.6%),thereby gradually decreasing from south to north.The multi-annual mean potential evapotranspiration in the JRB lies between 869.9 and 1,060.6 mm, with its spatial variation being similar to the multi-annual mean sunshine duration (5.4-7.7 h), thereby gradually increasing from southwest to northeast.The multi-annual mean maximum temperature (13.2-17.5 °C), minimum temperature (0.6-8.1 °C), and mean temperature (7.0-11.8°C) in the JRB exhibit similar spatial variation patterns, thereby gradually decreasing from southeast to northwest.The multiannual mean wind speed in the JRB lies between 1.6 and 5.9 m/s, with its spatial variation generally characterized as low in the east and high in the west.

Historical Climate Change Characteristics
Based on meteorological data from 1951 to 2019 from the JRB, this study used the IDW interpolation method to obtain the spatial-variation characteristics of each meteorological element in the river basin (Figure 3).The multi-annual mean precipitation in the JRB lies between 339.6 and 590.7 mm, with its spatial variation being similar to that of the multi-annual mean relative humidity (52.7-69.6%),thereby gradually decreasing from south to north.The multi-annual mean potential evapotranspiration in the JRB lies between 869.9 and 1,060.6 mm, with its spatial variation being similar to the multi-annual mean sunshine duration (5.4-7.7 h), thereby gradually increasing from southwest to northeast.The multi-annual mean maximum temperature (13.2-17.5 • C), minimum temperature (0.6-8.1 • C), and mean temperature (7.0-11.8• C) in the JRB exhibit similar spatial variation patterns, thereby gradually decreasing from southeast to northwest.The multi-annual mean wind speed in the JRB lies between 1.6 and 5.9 m/s, with its spatial variation generally characterized as low in the east and high in the west.

Future Climate Change Characteristics
Based on 23 grid points with spatial resolutions of 0.5 × 0.5 • in and around the JRB, the meteorological elements in the river basin from 2020 to 2050 under four climate scenarios were calculated, and the variation characteristics of annual precipitation and mean values of temperature under different climate scenarios were analyzed (Figure 4).Under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the multi-annual mean precipitations in the JRB are 505.2,474.4,515.2, and 521.9 mm, respectively (Figure 4a-d).Among them, the multi-annual mean precipitations under the RCP6.0 and RCP8.5 climate scenarios are relatively close to the multi-annual mean precipitation in the 1951-2019 historical period (517.9 mm), showing differences of −2.7 and 4.0 mm, respectively.Under the RCP4.5 climate scenario, the multi-annual mean precipitation exhibits the greatest difference (−43.5 mm) from that of the historical period.Under the four climate scenarios, the annual precipitations in the JRB exhibit increasing trends with rates of 1.360, 3.360, 4.420, and 0.114 mm/a, respectively.Under the RCP6.0 climate scenario alone, the annual precipitation shows a significant (p < 0.05) increasing trend.Under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the multi-annual mean temperatures in the JRB from 2020 to 2050 are 7.2, 7.3, 7.2, and 7.6 • C, respectively.Except for the RCP2.6 climate scenario, the annual mean temperature shows a significant (p < 0.05) increasing trend (Figure 4e-h), with the rate of increase in the RCP8.5 scenario being the highest (0.054 • C/a).Under the RCP2.6 scenario, the increase rate is the lowest (0.008 • C/a).Under the RCP4.5 and RCP6.0 climate scenarios, the increase rates are 0.029 and 0.021 • C/a, respectively.

Future Climate Change Characteristics
Based on 23 grid points with spatial resolutions of 0.5 × 0.5° in and around the JRB, the meteorological elements in the river basin from 2020 to 2050 under four climate scenarios were calculated, and the variation characteristics of annual precipitation and mean values of temperature under different climate scenarios were analyzed (Figure 4).Under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the multi-annual mean precipitations in the JRB are 505.2,474.4,515.2, and 521.9 mm, respectively (Figure 4a-d).Among them, the multi-annual mean precipitations under the RCP6.0 and RCP8.5 climate scenarios are relatively close to the multi-annual mean precipitation in the 1951-2019 historical period (517.9 mm), showing differences of -2.7 and 4.0 mm, respectively.Under the RCP4.5 climate scenario, the multi-annual mean precipitation exhibits the greatest difference (-43.5 mm) from that of the historical period.Under the four climate scenarios, the annual precipitations in the JRB exhibit increasing trends with rates of 1.360, 3.360, 4.420, and 0.114 mm/a, respectively.Under the RCP6.0 climate scenario alone, the annual precipitation shows a significant (p < 0.05) increasing trend.Under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the multi-annual mean temperatures in the JRB from 2020 to

Historical LUCC Characteristics
The spatial changes and area proportions of various types of LUCC in the JRB from 1980 to 2018 are shown in Figure 5.In the past 40 years, the farmland area in the JRB first increased and then decreased, reaching the maximum in 1990 and becoming smaller than the grassland area after 2010.The forest land area showed an oscillating pattern, reaching the minimum value in 2000, with a maximum increase of 10.59% from 2000 to 2010.The grassland area did not change before 2010, with a maximum increase of 0.62%, but increased rapidly after 2010.Water bodies have exhibited a decreasing trend in the past 40 years, especially in the period from 2010 to 2018, with a maximum decrease of 17.67%.In the past 40 years, the urban land area has shown an increasing trend, with a change of greater than 10% after 1990; of which the increase in 2018 was 34.91% compared to that in 2010.Overall, LUCC in the JRB over the past 40 years was characterized by a decrease in farmland and an increase in urban and forest land.

CA-Markov Model Applicability Analysis
To verify the applicability of the CA-Markov model in predicting LUCC in the JRB, this study constructed a suitability atlas of different LUCC types based on the DEM, slope, aspect, soil (Figure 1c-f), and historical climate (Figure 3) data.At the same time, based on the LUCC transition matrix in 1980-1990 and 1990-2000, combined with the measured LUCC in 1990 and 2000, the LUCC in 2000 and 2010 was predicted.Finally, the 2000 and 2010 forecasts were compared with the measured LUCC using the Crosstab module of the IDRISI 17 software, with the calculated Kappa coefficients being 0.948 and 0.939, respectively (Figure 6).The results showed that the CA-Markov model had high prediction accuracy for the LUCC in the JRB in 2000 and 2010, with predictions being consistent with observations, thereby meeting the accuracy requirements for LUCC predictions in the JRB.

Future LUCC Scenarios
The results of the forecasted LUCC in the JRB for 2025, 2035, and 2045 are shown in Figure 7 and Table 1.The LUCC of each type in each future period was indistinguishable from 2018, and the area-size order was as follows: grassland > farmland > forest land > urban land > water.Farmland in the JRB showed a decreasing trend.Compared with 2018, farmland area in 2025, 2035, and 2045 decreased by 56, 170, and 162 km 2 , respectively.Forest area increased by 3.05%, 6.30%, and 6.47% in 2025, 2035, and 2045, respectively.Grassland decreased by ~1%.Compared to 2018, urban land increased significantly by 11.40%.Except for 2045, water bodies remained unchanged.In addition, the characteristics of future LUCC in the river basin are as follows: the forest area in the eastern mountainous region increased, the farmland area in the middle and lower reaches decreased, and urban land expanded outward based on the original extent.In several LUCC types with a greater degree of change, an increase in the area of urban land favors an increase in the rate of runoff production; a decrease in the area of farmland favors a decrease in agricultural water use and an increase in the area of forest land, which has the effect of storing water and intercepting rainfall to effectively conserve water, resulting in a decrease in the amount of runoff.However, the magnitude of future changes in urban land is greater compared to forest land, and there is a greater likelihood that future LUCC will favor an increase in runoff.
from 2018, and the area-size order was as follows: grassland > farmland > forest land > urban land > water.Farmland in the JRB showed a decreasing trend.Compared with 2018, farmland area in 2025, 2035, and 2045 decreased by 56, 170, and 162 km 2 , respectively.Forest area increased by 3.05%, 6.30%, and 6.47% in 2025, 2035, and 2045, respectively.Grassland decreased by ~1%.Compared to 2018, urban land increased significantly by 11.40%.Except for 2045, water bodies remained unchanged.In addition, the characteristics of future LUCC in the river basin are as follows: the forest area in the eastern mountainous region increased, the farmland area in the middle and lower reaches decreased, and urban land expanded outward based on the original extent.In several LUCC types with a greater degree of change, an increase in the area of urban land favors an increase in the rate of runoff production; a decrease in the area of farmland favors a decrease in agricultural water use and an increase in the area of forest land, which has the effect of storing water and intercepting rainfall to effectively conserve water, resulting in a decrease in the amount of runoff.However, the magnitude of future changes in urban land is greater compared to forest land, and there is a greater likelihood that future LUCC will favor an increase in runoff.

Applicability Analysis of the SWAT Model
The simulation results of the calibration and verification periods in relation to the calibration and verification of the applicability of the SWAT model to the JRB are shown in Figure 8.The R 2 , NS, and RE of the observed and simulated annual runoff were 0.90, 0.89, and 3.33% in the calibration period and 0.73, 0.56, and -7.68% in the verification period, respectively.According to the results of each evaluation index, the SWAT model was more effective in simulating the annual runoff in the JRB.The simulation results of the calibration and verification periods in relation to the calibration and verification of the applicability of the SWAT model to the JRB are shown in Figure 8.The R 2 , NS, and RE of the observed and simulated annual runoff were 0.90, 0.89, and 3.33% in the calibration period and 0.73, 0.56, and −7.68% in the verification period, respectively.According to the results of each evaluation index, the SWAT model was more effective in simulating the annual runoff in the JRB.

Runoff Response to Future Climate Change
Based on the SWAT model, annual runoff changes in the JRB under the four climate scenarios were simulated and predicted (Figure 9).Under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the multi-annual mean runoff-flow rate from 2020 to 2050 was 45.80, 36.22,45.77, and 50.36 m³/s, respectively.Under the RCP2.6,RCP4.5, and RCP6.0 climate scenarios, all annual runoff-flow rates exhibited increasing trends, with increasing rates of 0.006, 0.672, and 1.135 m³/(s•a), respectively.Under the RCP4.5 and RCP6.0 climate scenarios, the annual runoff-flow rate showed significant (p < 0.05) increasing trends.Under the RCP8.5 climate scenario, the annual runoff-flow rate in the JRB exhibited a nonsignificant decreasing trend at a rate of −0.103 m³/(s•a).To further analyze the relationship between precipitation and runoff change in the JRB from 2020 to 2050, the future change in precipitation and runoff-flow rate under different scenarios were studied.The RCP4.5 climate scenario was selected since its precipitation was the lowest.The baseline analysis involved three other scenarios (Table 2).The results showed that the runoff-flow rate is consistent with the change in precipitation,

Runoff Response to Future Climate Change
Based on the SWAT model, annual runoff changes in the JRB under the four climate scenarios were simulated and predicted (Figure 9).Under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the multi-annual mean runoff-flow rate from 2020 to 2050 was 45.80, 36.22,45.77, and 50.36 m³/s, respectively.Under the RCP2.6,RCP4.5, and RCP6.0 climate scenarios, all annual runoff-flow rates exhibited increasing trends, with increasing rates of 0.006, 0.672, and 1.135 m³/(s•a), respectively.Under the RCP4.5 and RCP6.0 climate scenarios, the annual runoff-flow rate showed significant (p < 0.05) increasing trends.Under the RCP8.5 climate scenario, the annual runoff-flow rate in the JRB exhibited a nonsignificant decreasing trend at a rate of −0.103 m³/(s•a).To further analyze the relationship between precipitation and runoff change in the JRB from 2020 to 2050, the future change in precipitation and runoff-flow rate under different scenarios were studied.The RCP4.5 climate scenario was selected since its precipitation was the lowest.The baseline analysis involved three other scenarios (Table 2).The results showed that the runoff-flow rate is consistent with the change in precipitation, with the ratio of ΔR/ΔP being positive.In addition, under the RCP2.6,RCP6.0, and RCP8.5 climate scenarios, ΔR/ΔP was 0.31, 0.23, and 0.30 m³/(mm•s), respectively, indicating that the sensitivity of runoff-flow rate changes to precipitation changes was different under To further analyze the relationship between precipitation and runoff change in the JRB from 2020 to 2050, the future change in precipitation and runoff-flow rate under different scenarios were studied.The RCP4.5 climate scenario was selected since its precipitation was the lowest.The baseline analysis involved three other scenarios (Table 2).The results showed that the runoff-flow rate is consistent with the change in precipitation, with the ratio of ∆R/∆P being positive.In addition, under the RCP2.6,RCP6.0, and RCP8.5 climate scenarios, ∆R/∆P was 0.31, 0.23, and 0.30 m 3 /(mm•s), respectively, indicating that the sensitivity of runoff-flow rate changes to precipitation changes was different under different scenarios; that is, it was also affected by changes in meteorological factors, such as temperature and evaporation.The characteristics of hydrometeorological change in the JRB in different future periods under the four climate scenarios are shown in Figure 10.Under the combined influences of climate and LUCC in the future, the multi-annual mean runoff-flow rate from 2020 to 2050 under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios was 46.38, 36.81, 46.29, and 51.04 m 3 /s, respectively; that is, with the increase of carbon dioxide emission concentration, the runoff-flow rate first decreased and then increased, which is entirely consistent with the change in precipitation.In addition, with the increase in carbon dioxide emission concentration, the variation trends of runoff-flow rate and precipitation are completely consistent in the early (2020-2029) and late (2040-2050) periods and are consistent in the middle (2030-2039) period.Except for the RCP8.5 climate scenario, the maximum values of the multi-annual mean runoff-flow rate and precipitation of the other three climate scenarios all appeared in the later period; moreover, the maximum values of annual mean temperature of the four climate scenarios all appeared in the later period.This shows that the runoff change in the JRB is dominated by precipitation, and with the increase in carbon dioxide emission concentration, precipitation and runoff would increase in the future.
To further elucidate the characteristics of hydrometeorological changes in different periods, the discreteness of the hydrometeorological elements in the JRB was calculated (Table 3).Under the RCP6.0 climate scenario, the dispersion degree of the runoff-flow rate in the JRB was the highest, with the extreme-value ratio and coefficient of variation being 13.41 and 0.63, respectively.Under the RCP2.6 climate scenario, the dispersion degree of runoff-flow rate was the lowest, with the extreme-value ratio and coefficient of variation being 4.92 and 0.37, respectively.In addition, under the RCP2.6 climate scenario, the dispersion degree of the runoff-flow rate in the JRB was the highest in the middle period, with the coefficient of variation being 0.46.Under the RCP4.5 climate scenario, the dispersion degree of the runoff-flow rate was the highest in the earlier period, with the coefficient of variation being 0.43.Under the RCP6.0 and RCP8.5 climate scenarios, the dispersion degree of the runoff-flow rate was the highest in the later period, with the coefficients of variation being 0.61 and 0.73, respectively.that the change in runoff-flow rate in the JRB was not only affected by meteorological factors, such as precipitation and air temperature, but also by changes in LUCC and other factors.Further, with the increase in carbon dioxide emission concentration, the possibility of extreme precipitation and flood disasters in the future increases.Under the four climate scenarios, precipitation and mean temperature in the JRB in different periods were more concentrated than the runoff-flow rate.In particular, mean temperature had the lowest degree of dispersion.Under the RCP6.0 climate scenario, the maximum extreme ratios of precipitation and mean temperature from 2020 to 2050 were 2.44 and 1.34, respectively, whereas the coefficients of variation were less than those of the RCP8.5 climate scenario; that is, 0.20 and 0.07, respectively.Under the RCP2.6 climate scenario, the minimum coefficient of variation of precipitation and mean temperature were 0.16 and 0.06, respectively.The extreme-value ratio of precipitation was slightly higher than that of the RCP4.5 climate scenario; that is, 2.04, with the minimum extreme value ratio of mean temperature being 1.26.In addition, under the four climate scenarios, precipitation had the largest dispersion during the later period.The dispersion degree of mean temperature was the highest in the later period under the RCP2.6 and RCP4.5 climate scenarios, but the lowest under the RCP6.0 and RCP8.5 climate scenarios.This shows that the change in runoff-flow rate in the JRB was not only affected by meteorological factors, such as precipitation and air temperature, but also by changes in LUCC and other factors.Further, with the increase in carbon dioxide emission concentration, the possibility of extreme precipitation and flood disasters in the future increases.

Combined Effects of Climate Change and LUCC on Runoff Change
Under the influence of future climate alone, except for the RCP8.5 climate scenario, the multi-annual mean runoff-flow rate in the JRB decreased to a certain extent compared with the historical period (1961-2019).In addition, under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the contribution rate of climate change to the change in runoffflow rate in the JRB was −8.06%, −27.30%, −8.12%, and +1.10%, respectively (Table 4).Under the combined influence of future climate change and LUCC, the multi-annual mean runoff-flow rate in the JRB decreased by 6.79%, 25.66%, and 6.98% under the RCP2.6,RCP4.5, and RCP6.0 climate scenarios, respectively.Under the RCP8.5 climate scenario, the multi-annual mean runoff-flow rate in the JRB increased by 2.45%.Considering the runoffresponse processes in the A1 and A2 set-ups, the results showed that the contribution rate of the LUCC to the change in runoff-flow rate in the JRB was +1.27%, +1.64%, +1.14%, and +1.35% under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, respectively.A1 is the rate of change compared with the baseline period, which reflects the impact of climate change on runoff change; A2 is the rate of change compared with A1, which reflects the impact of LUCC on runoff change.
In general, influenced by the increasing trend of future precipitation, future runoff in the JRB showed increasing trends under the RCP2.6,RCP4.5, and RCP6.0 climate scenarios.However, under the RCP8.5 climate scenario, future precipitation in the JRB showed a nonsignificant increasing trend, with the temperature increasing significantly, thereby resulting in a decreasing trend in runoff.In addition, except for the RCP8.5 climate scenario, the multi-annual mean precipitation in the JRB was smaller than that in the historical period, with future climate change likely leading to a decrease in the multi-annual mean runoff in the river basin.Future LUCC will be conducive to the increase in the multi-annual mean runoff in the JRB, albeit with a relatively small impact; therefore, under the combined influence of future climate change and LUCC, it is more likely that the multi-annual mean runoff in the JRB will decrease in the future.

Uncertainty Analysis
The existing studies do not consider LUCC and vegetation species changes that influence the hydrological cycle processes in watersheds.It is also difficult to accurately strip out the effects of LUCC and vegetation-species changes on hydrological cycle processes, especially the effects of LUCC, in practical applications [47].In previous studies, most of the drivers of hydrological change were divided into two major factors: climate change and human activities, and first-order approximations were used, such as, LUCC was approximated as a result of human activities [48].However, changes in LUCC and vegetation species alter hydrological processes such as filling, infiltration, and soil water-storage capacity during the runoff-generation process in the watershed and surface roughness, surface water storage, and river catchment path during confluence process, which lead to increased canopy interception and increased conversion of precipitation to evapotranspiration [49].Therefore, human activities play an important role in altering the hydrological cycle of a watershed [50].However, there are relatively few studies on the mechanism of change in hydrological processes in watersheds driven by LUCC, which need to be further strengthened at a later stage.
The SWAT and CA-Markov model used in this study all have certain internal systematic errors, and although they act as reference values for future studies of runoff change and its driving factors, their accuracy still needs to be further improved [51,52].Therefore, the next step needs to comprehensively and accurately identify the depressions characteristics of DEM, develop and improve the underlying code of the SWAT model, and combine with more powerful computational capabilities and high-precision LUCC data, so as to obtain the runoff-change characteristics and the reliable future-prediction results that better reflect the actual watershed [53,54].
In addition, research on climate change and its impacts belongs to the 'if-then-what' type.Future emission scenarios and hydrological model selection and future LUCC will increase the uncertainty of climate-change impact on runoff [55,56].Despite the uncertainties associated with modeling predictions of future climate impacts on runoff changes based on climate-model outputs, specific scenarios remain fundamental to helping identify possible future climate and hydrologic changes [57].Previous studies have shown that the choice of GCM relative to emission scenarios and downscaling methods is the largest contributor to uncertainty in climate-change impact evaluation [58,59].Therefore, when evaluating the impact of climate change in a specific region on future hydrological and water resources, it is crucial to select a suitable future climate model and to simulate future hydrological processes based on multiple climate scenarios to make the resulting conclusions more reliable.

Strategies for Rational Utilization of Water Resources in the Future
Under the four future scenarios, the future runoff alone in the JRB under the RCP8.5 scenario is slightly higher than that of the historical period, whereas the remaining scenarios are smaller than those of the historical period, especially the multi-year-average runoff in the basin under the RCP4.5 scenario decreasing by 27.30%, which indicates that climate change in the future is likely to exacerbate water scarcity in the basin.However, this study shows that the future LUCC will be beneficial to alleviate the pressure on water resources.If policy favors solving the contradiction between the supply and demand of water resources in the watershed, it should continue with the current policy; that is, it should increase the area of future urban land in accordance with the development trend of the present day, and the land retired from farmland should be changed into grassland more often.If policy favors alleviating the pressure on flood control in the watershed, it should control the continuous increase in the area of future urban land, and the retired farmland should be changed into forest land more often [22].In addition, water resources in the JRB are concentrated in the flood season, accounting for more than 60% of the annual runoff.In the future, under the influence of changing climatic conditions, the distribution of water resources within the year will be increasingly uneven, and it is recommended to build large reservoirs and groundwater reservoirs in and around the basin so as to realize the optimal allocation of water resources through the joint storage of the two, to feed the dryness with the abundance and to improve the utilization rate of water resources [60].

Limitations of the Study
This study predicts the first level-LUCC types in the JRB alone, and its spatial resolution is 1 × 1 km, which needs to be improved in the detailed portrayal of LUCC.However, the actual individual LUCC is characterized by scattered and small patches of change, which may have implications for predicting future LUCC in the watershed [61,62].In addition, this study quantitatively distinguishes the contribution of LUCC and climate change based on the results of runoff simulation under different scenarios but does not further investigate the mechanism of their influence; nor does it quantify the sensitivity of runoff to each level of LUCC type, which still requires further research.In addition, the influence of human activities on runoff change in the JRB has gradually increased in recent years.According to the statistics, the surface water supply in the JRB is 270 million m 3 , accounting for 56.72% of the total water consumption, and 65.83% of the surface water supply is supplied in the form of water diversion, and more than 80% of the diverted water is used for agricultural irrigation (Figure 11).However, the impact on future water use owing to the reduction in farmland area was not considered in this study, which may have an impact on future runoff changes.
of the present day, and the land retired from farmland should be changed into grassland more often.If policy favors alleviating the pressure on flood control in the watershed, it should control the continuous increase in the area of future urban land, and the retired farmland should be changed into forest land more often [22].In addition, water resources in the JRB are concentrated in the flood season, accounting for more than 60% of the annual runoff.In the future, under the influence of changing climatic conditions, the distribution of water resources within the year will be increasingly uneven, and it is recommended to build large reservoirs and groundwater reservoirs in and around the basin so as to realize the optimal allocation of water resources through the joint storage of the two, to feed the dryness with the abundance and to improve the utilization rate of water resources [60].

Limitations of the Study
This study predicts the first level-LUCC types in the JRB alone, and its spatial resolution is 1 × 1 km, which needs to be improved in the detailed portrayal of LUCC.However, the actual individual LUCC is characterized by scattered and small patches of change, which may have implications for predicting future LUCC in the watershed [61,62].In addition, this study quantitatively distinguishes the contribution of LUCC and climate change based on the results of runoff simulation under different scenarios but does not further investigate the mechanism of their influence; nor does it quantify the sensitivity of runoff to each level of LUCC type, which still requires further research.In addition, the influence of human activities on runoff change in the JRB has gradually increased in recent years.According to the statistics, the surface water supply in the JRB is 270 million m 3 , accounting for 56.72% of the total water consumption, and 65.83% of the surface water supply is supplied in the form of water diversion, and more than 80% of the diverted water is used for agricultural irrigation (Figure 11).However, the impact on future water use owing to the reduction in farmland area was not considered in this study, which may have an impact on future runoff changes.

Conclusions
In this study, the JRB was selected as the study area, the historical (1951-2019) and future (2020-2050) climate change characteristics were analyzed, and the LUCC in the JRB in 2025, 2035, and 2045 was predicted.The law and attribution of runoff change under different combination scenarios were studied, and the main conclusions are as follows: (1) Under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the annual precipitation and mean temperature in the JRB from 2020 to 2050 increases at a rate of 1.36 and 0.008, 3.36 and 0.029, 4.42 and 0.021, and 0.114 mm/a and 0.054 • C/a, respectively.Under the four scenarios, the change of precipitation was not significant, whereas the temperature change was significant (p < 0.05), except under the RCP2.6 scenario.Compared with 2018, forest and urban land areas in the JRB in 2025, 2035, and 2045 increase by 3.05% and 9.46%, 6.30% and 10.05%, and 6.47% and 11.40%, respectively, whereas grassland, farmland, and water bodies slightly decrease.
(2) Based on the results of each evaluation index, the SWAT model was more efficient in simulating the annual runoff change in the JRB.Under the RCP2.6,RCP4.5, RCP6.0, and RCP8.5 climate scenarios, the contribution rates of future climate change to runoff change in the JRB were −8.06%, −27.30%, −8.12%, and +1.10%, whereas the contribution rates of the future LUCC were +1.27%, +1.64%, +1.14%, and +1.35%, respectively.In addition, with the increase in carbon dioxide emission concentration, future precipitation and runoff also increased, thereby increasing the possibility of extreme precipitation and flood disaster events.
(3) There are many uncertainties in the impacts of LUCC and climate change on runoff changes in the JRB, which need to be explored by continued in-depth studies in terms of quantification of uncertainties and different LUCC types of sensitivities.In addition, the future impacts of LUCC and climate change on runoff vary with different climate characteristics and watersheds.Therefore, there is a need to apply the research framework used in this study to other watersheds, thus contributing to future waterresource management in a wider range of watersheds.

Data Availability Statement:
The digital elevation model, soil raster maps, land use, and future climate data used in this study are freely available and can be accessed from the website provided in the data section of the manuscript.Historical climate and runoff data were obtained from the China Meteorological Data Service Centre and Yellow River Conservancy Commission and can be obtained from these departments through official channels.

Figure 2 .
Figure 2. Schematic of the framework used to determine the impacts of climate change and LUCC on runoff.

Figure 2 .
Figure 2. Schematic of the framework used to determine the impacts of climate change and LUCC on runoff.

3. 2 .
Characterization and Prediction of LUCC 3.2.1.Historical LUCC Characteristics The spatial changes and area proportions of various types of LUCC in the JRB from 1980 to 2018 are shown in Figure 5.In the past 40 years, the farmland area in the JRB first increased and then decreased, reaching the maximum in 1990 and becoming smaller than the grassland area after 2010.The forest land area showed an oscillating pattern, reaching the minimum value in 2000, with a maximum increase of 10.59% from 2000 to 2010.The grassland area did not change before 2010, with a maximum increase of 0.62%, but increased rapidly after 2010.Water bodies have exhibited a decreasing trend in the past 40 years, especially in the period from 2010 to 2018, with a maximum decrease of 17.67%.In the past 40 years, the urban land area has shown an increasing trend, with a change of greater than 10% after 1990; of which the increase in 2018 was 34.91% compared to that in 2010.Overall, LUCC in the JRB over the past 40 years was characterized by a decrease in farmland and an increase in urban and forest land.Atmosphere 2023, 14, 1237 8 of 20

Figure 6 .
Figure 6.Results of LUCC projections and their differences in the JRB: (a) and (b) in 2000; (c) and

Figure 5 .
Figure 5. Historical LUCC characteristics and statistics in the JRB: (a) 1980, (b) 1990, (c) 2000, (d) 2010, (e) 2018, and (f) characteristics of change statistics.3.2.2.CA-Markov Model Applicability AnalysisTo verify the applicability of the CA-Markov model in predicting LUCC in the JRB, this study constructed a suitability atlas of different LUCC types based on the DEM, slope, aspect, soil (Figure1c-f), and historical climate (Figure3) data.At the same time, based on the LUCC transition matrix in 1980-1990 and 1990-2000, combined with the measured LUCC in 1990 and 2000, the LUCC in 2000 and 2010 was predicted.Finally, the 2000 and 2010 forecasts were compared with the measured LUCC using the Crosstab module of the IDRISI 17 software, with the calculated Kappa coefficients being 0.948 and 0.939, respectively (Figure6).The results showed that the CA-Markov model had high prediction accuracy for the LUCC in the JRB in 2000 and 2010, with predictions being consistent with observations, thereby meeting the accuracy requirements for LUCC predictions in the JRB.

Figure 6 .
Figure 6.Results of LUCC projections and their differences in the JRB: (a) and (b) in 2000; (c) and (d) in 2010.Figure 6. Results of LUCC projections and their differences in the JRB: (a,b) in 2000; (c,d) in 2010.

Figure 6 .
Figure 6.Results of LUCC projections and their differences in the JRB: (a) and (b) in 2000; (c) and (d) in 2010.Figure 6. Results of LUCC projections and their differences in the JRB: (a,b) in 2000; (c,d) in 2010.

20 Figure 8 .
Figure 8.Comparison of the observed and simulated runoff in the calibration and verification periods in the JRB: (a) simulation effect of runoff flow rate and (b) correlation between observed and simulated runoff flow rates.

Figure 8 .
Figure 8.Comparison of the observed and simulated runoff in the calibration and verification periods in the JRB: (a) simulation effect of runoff flow rate and (b) correlation between observed and simulated runoff flow rates.

Figure 8 .
Figure 8.Comparison of the observed and simulated runoff in the calibration and verification periods in the JRB: (a) simulation effect of runoff flow rate and (b) correlation between observed and simulated runoff flow rates.

Figure 11 .
Figure 11.Statistics of water use in the JRB: (a) surface water supply and (b) water consumption.Figure 11.Statistics of water use in the JRB: (a) surface water supply and (b) water consumption.

Figure 11 .
Figure 11.Statistics of water use in the JRB: (a) surface water supply and (b) water consumption.Figure 11.Statistics of water use in the JRB: (a) surface water supply and (b) water consumption.

Table 1 .
LUCC area statistics in the JRB during different periods.

Table 1 .
LUCC area statistics in the JRB during different periods.
1-5 represent farmland, forest land, grassland, water, and urban land, respectively.3.3.Quantitative Evaluation of the Impact of Climate Change and LUCC on Runoff Change 3.3.1.Applicability Analysis of the SWAT Model

Table 2 .
Analysis of the future relationship between precipitation and runoff-flow rate changes.

Table 3 .
Discrete characteristics of hydrometeorological elements in different periods under different climate scenarios.

Table 3 .
Cont. is the extreme value ratio and C v is the coefficient of variation. K

Table 4 .
Quantitative identification of the contribution rates of future climate change and LUCC to runoff.