Assessment of Rainfall And Temperature Trends in The Yellow River Basin, China from 2023 to 2100

: China's Yellow River Basin (YRB) is sensitive to climate change due to 12 its delicate ecosystem and complex geography. Water scarcity, soil erosion, and 13 desertification are major challenges. To mitigate the YRB's ecological difficulties, 14 climate change must be predicted. Based on the analysis of the evolution features of 15 hydro-meteorological elements, the CMIP6 climate model dataset with Delta 16 downscaling and Empirical Orthogonal Function (EOF) is utilized to quantitatively 17 explore the future variations of precipitation and temperature in the YRB. The 18 following results are drawn: The spatial resolution of the CMIP6 climate model is less 19 than 0.5°×0.5° (i.e., about 55 km×55 km), which is improved to 1 km×1 km by the 20 downscaling of Delta, and has outstanding applicability to precipitation and 21 temperature in the YRB. The most accurate models for monthly mean temperature are 22 CESM2-WACCM, NorESM2-LM, and ACCESS-CM2, and for precipitation are 23


Introduction
The world nowadays is extremely concerned about global warming.The Intergovernmental Panel on Climate Change (IPCC) noted in its Sixth Assessment Report that global warming has been significant since the industrial revolution (Masson-Delmotte et al., 2018).Climate warming has become a widely acknowledged fact over the past 70 years as evidenced by the rise in the average global temperature, the melting of snow glaciers, and the rise in sea level.Nearly every area of human life is impacted by climate change, which is directly tied to human society.As a result, the future trend of climate change is a crucial area for research for all of humanity.One of the most severely impacted regions in China is the YRB, which is one of the sensitive areas and has a warming rate that is higher than the global average for the same time period (Omer et al., 2020).Droughts, floods, and heat waves have all significantly risen in the YRB as just a result of global warming (Li et al., 2022;Zhou et al., 2019), making this region more severely affected by climate change than most other regions of the globe.Consequently, it is critical that the YRB undertake future climate change research.
Presently, future climate change is predicted using global climate models (GCMs).
Before making predictions, it is vital to evaluate the models' simulation capabilities due to their complexity.Based on CMIP5 (Coupled Model Intercomparison Project Phase 5), the capacity of climate models to simulate elements of the global or regional climate has been assessed.The findings demonstrate that practically all models predict an increase in future temperatures in various places (Zhou et al., 2019;Zhang et al., 2019).

In comparison to CMIP3 and CMIP5, the International Coupled Model Comparison
Program, which is already in its sixth phase (CMIP6), compared to CMIP3 and CMIP5, has the most models participating, the best-designed scientific experiments, the most comprehensive simulation data provided, and more advanced advancements in quantity, resolution, and accuracy.The Shared Socioeconomic Pathway and Representative Concentration Pathway (SSP-RCP) combined scenarios, a new generation of combined emission scenarios introduced by the CMIP6 model, will serve as the foundation for future study on the global climate during the next 5-10 years.It has been demonstrated that using the CMIP6 model dataset will result in simulation results that are more in line with observations than CMIP5 (Chen et al., 2021;Fu et al., 2021).The current research on the application of CMIP6 in the YRB is relatively insufficient, so the study of future climate change in this basin can develop better related water conservation strategies for the basin.
Meteorological data can be used to study climate trends, but gridded climate datasets are currently the most popular and practical option.These datasets can better simulate large-scale climate characteristics, but their highest spatial resolution is only 0.5° (roughly 55 km), making it impossible to depict the impact of topography and landscape on regional climate.Improving the spatial resolution of the gridded climate dataset by spatial downscaling is the current remedy for the aforementioned issue (Peng et al., 2018).Statistical downscaling and dynamical downscaling are two common spatial downscaling methods for climate data.Most research use statistical downscaling, which doesn't require complicated physical formulae (Wang et al., 2022).Dynamic downscaling, in contrast, necessitates a large number of parameters, is driven by higher computer power, and occasionally fails to faithfully represent the characteristics of climate change at tiny geographic scales (Gebrechorkos et al., 2019).
Utilizing temperature and precipitation as logical indicators can help scientists assess the dangers of climate change and develop climate change responses that are successful (Mohammad et al., 2018).On the basis of the CMIP5 model, studies have been done to simulate the temperature and precipitation in the Chinese region (Zhou et al., 2019;Zhang et al., 2019).The outcomes of CMIP6's climate prediction differ from those of CMIP5 (Eyring et al., 2019;Wang et al., 2022).Particularly, due to the physical upgrade of the climate model, the warming in the SSPs scenario in CMIP6 is often greater than that in the RCPs scenario in CMIP5, and the expected increases in regional total and extreme precipitation in China tend to be higher than those in CMIP5 (Zhu et al., 2021).Additionally, it has been demonstrated that CMIP6 models generally perform better when reproducing mean temperature and precipitation, as evidenced by lower mode differences and less cold bias (Jiang et al., 2020;Yu et al., 2022).The multi-model ensemble (MME), which is based on the synthesis of various CMIP6 models, has been used in numerous studies for climate models since it can successfully minimize the uncertainty of individual models and has superior simulation capabilities than most models (Xu et al., 2021;Yan et al., 2021).
In this study, we optimize the CMIP6 climate model data using historical data, and we compare various spatial interpolation methods in downscaling analysis based on different evaluation indicators, in order to select several climate models with the best performance results for arithmetic averaging, and ultimately obtain a comprehensive MME.The EOF is used to explore the temporal distribution of the annual temperature and based on the ideal MME in order to predict and analyze the trends of climate change in the YRB over the coming years (2023-2100).

Study area
The YRB is situated in northern China (32°-42°N, 95°-120°E), and is approximately 1900 km long from east to west and 1100 km wide from north to south, with a basin area of 7.95×10 5 km² (Guo et al., 2022) (Fig. 1).The area has a continental climate, vast territory and numerous mountain ranges.With a semi-humid climate in the south, a semi-arid climate in the middle, and an arid climate in the northwest, the variations in topography and terrain cause an unequal distribution of precipitation in the basin.The annual precipitation in the basin ranges from 200 to 650 mm and decreases from southeast to northwest.The East Asian monsoon influences regular floods in the region, and the annual temperature ranges from 4 to 12 °C, with a decreasing trend from east to west.The YRB's land use types are primarily grassland (49.6%) and farmland (27.2%), with forest land (9.3%) and unused land following (13.9%).Vegetation types range from east to west and include the deciduous broadleaved forest belt, grassland belt, desert belt, and Tibetan Plateau vegetation belt.(3) The CMIP6 climate model data is obtained from the Earth System Grid Federation (ESGF) (https://esgf-node.llnl.gov/projects/cmip6/)and include four future periods (2021-2040, 2041-2060, 2061-2080, 2081-2100) for four Shared Socioeconomic Pathways (ssps), which are different emission scenarios of ssp1-2.6 (ssp126), ssp2-4.5 (ssp245), ssp3-7.0(ssp370) and ssp5-8.5 (ssp585).The simulated datasets of monthly mean temperature (tas) and monthly rainfall (pr) for 39 climate models are selected (Table 1).

Model performance evaluation indicators
As part of the assessment of model simulation capabilities, we compare the spatial distribution of elements including temperature and precipitation over climate averages, and analyze the annual average variability of each model meteorological element over the selected time period.
(1) Mean Absolute Error (MAE) MAE is the average for the sum of absolute errors between all simulations and corresponding observations, which is better than the average error and effectively prevents errors to offset each other, thus correctly indicating the extent of the overall deviation of the prediction (Wang et al., 2014).It can be mathematically written as： where: Pi is the downscaled or raw value, Oi is the observed value, and n is the number of records of all validated stations.
(2) Time Skills Score (TS) SS can be used to quantitatively evaluate the effect of individual climate models on the simulation of temporal variability, which measures the difference between the simulated and observed interannual variability of the time series at each image element in the study area as a ratio of standard deviations (Zhu et al., 2021).It can be mathematically written as： where: STD m 、 STD o are the standard deviation values of the time series at each image point in the simulated and observed fields, respectively; according to the evaluation criteria, the closer the value of TS is to 0, the smaller the difference between the standard deviation of the observed and simulated fields is, indicating that the simulation ability of the climate model is better, and vice versa, the worse the simulation performance of the climate model is.
(3) Taylor diagram-based evaluation method (S) The three indicators, correlation coefficient, standard deviation and central root mean square error of the combined model and measured series, are plotted in polar coordinates as a scatter plot to visually judge the simulation effectiveness of the climate model (Zhu et al., 2021).In order to quantify the evaluation results of Taylor diagram, the standard deviation and correlation coefficient of the simulated and observed fields are introduced as indicator S to further quantitatively evaluate the simulation capability of the climate model.It can be mathematically written as： where, R is the correlation coefficient between each climate model and the observed data; Ro is the maximum of the correlation coefficients of all the evaluated objects; σ f = STD m /STD  ，STD m and STD o are the standard deviations of the simulated and observed fields, respectively.It can be seen that the more consistent the simulated results are with the observed data, the more S tends to be close to 1, and the better the simulation effect of the climate model.
(4) Spatial Skills Score (SS) The spatial skill score is used to measure the ability of each climate model to simulate the spatial distribution characteristics of the evaluation variables in the same period, taking into account the spatial mean deviation and correlation coefficient between the simulated and observed fields.It can be mathematically written as： ( ) where, SS is the dimensionless index; m and o are the data of simulated and observed fields, respectively; MSE(m,o) and MSE(o ‾,o) denotes the mean square error of simulated and observed fields, respectively; r m,o is the correlation coefficient of simulated and observed fields; s m and s o are the mean squared errors of the simulated and observed fields, respectively.It can be seen that the closer the SS is to 1, the better the simulation effect of the climate model.

Statistical downscaling
In this study, observations from 93 ground-based meteorological stations in the YRB area are used as independent variables, climate model simulations of monthly precipitation and temperature at the corresponding stations are used as dependent variables by the application of Delta downscaling using a spatial raster of multi-year averages of January-December precipitation and temperature with a resolution of 1 km for the base period of 1971-2000.Delta downscaling is a method to obtain future climate scenarios based on the series of measured climate elements in the base period and the characteristic values of variables (e.g., absolute increase in temperature, relative rate of change in precipitation, etc.) of regional future climate models, which are superimposed on each other.This method has the advantage of being relatively simple and less computationally intensive, and is capable of downscaling GCM simulations to specific stations.Different downscaling methods are used for different variables.It can be mathematically written as： where, Tf is the temperature downscaling data for the prediction period, T Gf is the simulated temperature data for each grid point in the study area during the prediction period, TGo is the simulated average temperature data for the grid points in the study area during the base period, and To is the measured average temperature data for the corresponding grid points during the base period, Pf is the precipitation data after the Delta downscaling, PGf is the simulated precipitation data at each grid point in the study area during the predict period, PGo is the simulated precipitation data at the grid point in the study area during the base period, and Po is the measured precipitation data at the corresponding grid point during the base period.The specific downscaling process is shown in the figure (Fig. 2).

Space field analysis
The empirical orthogonal function (EOF) is a method for analyzing the structural features in matrix data and extracting the main data feature quantities (Sang et al,. 2021).
When the PCA method was introduced into the analysis of meteorological problems, the meteorological field referred to the method as EOF decomposition.It can be mathematically written as： where, V is the spatial eigenvector corresponding to the EOF input variable X, Λ is the diagonal matrix of eigenvalues, and the contribution rate of the corresponding V is calculated based on Λ. T is the time coefficient matrix, which is obtained by multiplying the spatial eigenvector V with the input variable X, and finally decomposing X into the product of the spatial eigenvector and the time coefficient matrix.It can be mathematically written as：

Downscaling results evaluation and multimodal ensemble
For temperature, the BI shows minimum error among the four different evaluation indicators, so this interpolation method is used to further evaluate the effect of   there are some points with poor fitting.For temperature, the R 2 for the CESM2-WACCM, NorESM2-LM, and ACCESS-CM2 observed and simulated temperatures are 0.934, 0.933, and 0.932, respectively, and the regression coefficients are 0.967, 0.971, and 0.931, with most points near the 1:1 line.The fit results of the three climate models are similar, and the R 2 and regression coefficients are greater than 0.9, all passing the 99% confidence level test.Arithmetic mean ensemble of the preferred temperature and precipitation climate models is used to obtain the MME, which performs better than the other single models in the evaluation test, with its MAE dropping to 21.874 mm and both SS and S improving significantly, indicating that MME provides a better fit.The simulated and observed values of the MME dataset are more concentrated than the three independent climate models, which can reduce the errors caused by outliers to a certain extent and improve the fit of the simulated data (Fig. 5 and 6).

Trends in precipitation
The multi-year mean precipitation from 1901 to 2014 is 465.78 mm, with a standard deviation of 31.24mm.The annual precipitation fluctuates widely up and down, with an average variation rate of -1.45 mm/10a.The annual precipitation in different scenarios under the future scenario has similar interannual variation trends, which mainly shows that the fluctuations of precipitation phases in different scenarios are generally consistent, and there is a lagging or early relationship in the appearance of precipitation peaks in some years.In the four ssp126, ssp245, ssp370, and ssp585 scenarios, the increase rates of annual precipitation are 6.89 mm/10a, 5.31 mm/10a, 7.02 mm/10a, and 10.18 mm/10a, respectively.All four scenarios exhibit a strong growing tendency.Overall, there are no discernible changes between the annual precipitation fluctuation trends in the past and the future (Fig. 7).

Trends in temperature
The multi-year mean temperature from 1901 to 2014 is 6.59 °C, with a standard deviation of 0.59 °C and an average change rate at 0.0038 °C/10a.Similar to precipitation, the annual temperature in the YRB under the four future scenarios shows a more obvious increasing trend compared with the history, and a clear linear tendency.
The temperature increase rates under the four scenarios of ssp126, ssp245, ssp370, and ssp585 are 0.1 °C/10a, 0.3 °C/10a, 0.5 °C/10a, and 0.7 °C/10a, respectively, and the temperature increases with the emissions concentration in different scenarios.Under the ssp126 scenario, the average change rate is 0.1 ℃/10a, which will exceed 8 ℃ in 2100.Under the ssp585 scenario, the interannual average change rate of annual temperature is 0.7 °C/10a (P < 0.01), and the R 2 exceeds 0.98, indicating a significant increase in annual temperature with a mainly linear tendency, reaching a maximum of 12.87 °C within 100 years in 2100.The ssp245 and ssp370 have the closest temperature trends to the 1980-2014 period (Fig. 10).
The EOF1 under all four climate scenarios shows positive region-wide eigenvalues.
In the ssp126 scenario, the variance contribution of the EOF1 in the YRB temperature is 82.43%, which plays a major role.The variance contribution of EOF2 is 9.49%, while the contribution of the EOF1 in the other three scenarios is above 95%.The EOF1 eigenvalues of different scenarios are positive, indicating that the basin-wide temperature shows an increasing trend.The EOF1 in the two low emission scenarios shows a relatively obvious east-west distribution, and the positive center of the ssp245 is obviously south compared with the ssp126, and the eigenvalues are smaller at the lower reaches of the Yellow River, while the EOF1 in the two high emission scenarios shows a more similar north-south distribution.The spatial trend of the EOF2 of the ssp126 shows a more obvious east-west difference, reflecting a secondary feature of opposite temperature trends in the basin above the XiaHeYan and in the middle and lower reaches (Fig. 11).
The PC1 of different scenarios shows a linear increasing trend, and the temperature increases faster with increasing carbon emission concentration from the low emission scenario ssp126 and ssp245 to the high emission scenarios ssp370 and ssp585.Each EOF1 show negative values until the 2050s and 2060s, indicating that the temperature of the YRB is relatively low to the multi-year mean temperature, while it is relatively high after the 2060s, with distinctly cooler and warmer periods.In contrast to PC1, PC2 of ssp126 has no obvious trend and alternates frequently between warm and cold periods (Fig. 12).

Improvements in climate model preference
The regional terrain, subsoil, and atmospheric conditions-all of which have clear regional characteristics-have a substantial impact on temperature and precipitation.In In addition to meeting the needs of this study, the Delta method can be used to further downscale the climate data to obtain higher resolution (30 m×30 m) for the study of climate trends in the YRB, and these high-resolution climate data can be used as driving data for other fields like hydrological simulation, vegetation succession simulation, crop growth simulation, etc.

Future trends of precipitation and temperature in the YRB
Numerous studies have indicated that future average precipitation and temperature in China will increase (Tian et al., 2021;Zhu et al., 2021) and temperature relative to 1900-2014, with annual precipitation growth rates of 0.689 mm/a and 1.018 mm/a and temperature growth rates of 0.01 °C/a and 0.07 °C/a for the ssp126 and ssp585 scenarios, respectively.According to these findings, the entire YRB will grow more humid and warmer in the future.
The YRB region of China has a continental monsoon climate with contemporaneous rain and heat, and yearly precipitation is governed mostly by summer precipitation, which is influenced by both the East Asian monsoon and the Indian monsoon (Zhu et al., 2016).Therefore, the unique properties of each mode can be understood in terms of precipitation formation.On the basis of the EOF1 of precipitation, it is evident that all scenarios exhibit a trend of increasing precipitation across the region, as indicated by consistently higher or lower precipitation across the region, with larger amplitudes in the southeast and smaller amplitudes in the northwest and northeast.The EOF1 of several scenarios demonstrates a constant increase in temperature across the YRB, and all northern regions have higher typical values, which is consistent with the global warming trend process.For the low emission scenarios ssp126 and ssp245, the positive center of the EOF1 is located in the northern portion of the basin; for the high emission scenarios ssp370 and ssp585, it is located in the northern portion of the basin and the Yellow River source region.The second mode EOF2 of ssp126 reflects the secondary feature of the opposite temperature trend in the basin above the XiaHeYan and the middle and lower reaches of the basin, which may be related to the fact that the basin above the XiaHeYan is distributed with mountains and plateaus, having a higher elevation and being more sensitive to temperature fluctuations.
Comparing the predicted and observed data of the CMIP6 downscaled climate model reveals that the simulated values are generally smaller than the observed values, and the simulation performance of different climate models for different climate variables in the same region varies greatly, with a good fit for temperature but a poor fit for precipitation, necessitating separate verification of the regional applicability under each variable.In addition, it has excellent simulation capabilities for modeling interannual variability, although there is still a discrepancy between its simulated and actual values due to the climate model's atmospheric physical framework and climate variability.
In general, the applicability of the CMIP6 downscaled climate model over the YRB is higher than that of the CMIP6 unscaled and CMIP5 downscaled climate models, and its uncertainty still needs to be further examined.According to the CMIP6 downscaled climate model, the temperature of the YRB shows a significant increase in the future period, while precipitation mainly shows a significant increase, but with large interannual variability, indicating that the YRB is affected by global warming and tends to be warm and dry, which may exacerbate the water scarcity issue in the YRB.The In general, the MME is superior than the single climate simulation, and the CMIP6 climate model is relatively better for the monthly temperature simulation of the YRB with MME, while the precipitation fitting is poor.
Precipitation and temperature in the YRB during the 21st century are divided into two periods according to the trends, namely, the period of low rainfall and colder temperatures before 2060 and the period of high rainfall and warmer temperatures after 2060.The EOF1 for all scenarios of precipitation indicates a trend of increasing precipitation over the region, which is consistent with either high or low rainfall across the region, with larger amplitudes in the southeast and smaller amplitudes in the northwest and northeast.The EOF1 for each scenario of temperature shows a significant rising trend, and the warming amplitude grow as the emission concentration rises.The ssp126 and ssp245 regional warming centers are concentrated in the northern portion of the basin, whereas the ssp370 and ssp585 regional warming centers are concentrated in the northern portion of the basin and the Yellow River source region.
The climate model of CMIP6 after Delta downscaling has uncertainties in the simulation of precipitation and temperature in the YRB and performs well in interannual trends, but there are still discrepancies between specific annual numerical simulations and reality.Therefore, in future research, additional downscaling or bias correction methods and alternative high-resolution data sources can be utilized to further minimize the climate model's uncertainty.

Fig. 1
Fig.1 Location and DEM of the YRB and 93 hydrological stations distributed across the rectangular

Fig. 2
Fig.2 Delta downscaling process of climate models in the YRB temperature simulation.The ACCESS-CM2, ACCESS-ESM1-5, CASMS-CSM1-0, CESM2, CESM2-FV, E3SM-1-0, INM-CM-5-0, CanESM5, CESM2-WACCM, FIO-ESM-2-0, INM-CM5-0, NorESM2-LM, etc., rank higher among all climate models for each evaluation indicators, with simulated and observed MAE in the range of 23-26 mm, SS in the range of 0.16-0.39,S in the range of 0.82-1, and TS in the range of 0-0.005.It can be seen from the evaluation criteria of the four evaluation indicators that each higher-ranked climate model has a better performance in the temperature simulation performance level of the YRB, indicating that these models can better reflect the temperature characteristics of this region.Considering the data integrity of the climate models and the performance of the actual fits, three climate models, CESM2-WACCM, NorESM2-LM, and ACCESS-CM2, are used to study the future temperature trends in the YRB (Fig.3).

Fig. 3
Fig.3The fitting performance ranking of the monthly precipitation series in the YRB from 1995 to 2014

Fig. 4
Fig.4The fitting performance ranking of the monthly temperature series in the YRB from 1995 to 2014

Fig. 5 Fig. 6
Fig.5 Scatter plots of monthly precipitation (pr) sequences in the YRB from 1995 to 2014 based on

Fig. 7 Fig. 8
Fig.7 Changes in annual precipitation in the YRB from the historical period (1901-2014) to the future

Fig. 9
Fig.9The main modes of EOF of rainfall spatial field in the YRB for 2022-2100 under scenarios

Fig. 10
Fig.10 Time coefficients (PC1-PC3) corresponding to the main modes of EOF of the spatial field of

Fig. 11 Fig. 12
Fig.11 The main modes of EOF of the temperature spatial field in the YRB for 2022-2100 under scenarios ssp126 (a, c), ssp245 (e), ssp370 (b) and ssp585 (d); where the EOF1_tas_ssp126 (contribution %), EOF1_tas_ssp245 (contribution %), EOF1_tas_ssp370 (contribution %), EOF1_tas _ssp585 (contribution %) denote the 1st mode of temperature (tas) in ssp126, ssp245, ssp370, and ssp585 scenarios, respectively, and EOF2_tas_ssp126 (% contribution) denotes the 2nd mode of temperature (tas) in ssp126 scenario, and the contribution rate of each mode is shown in parentheses order to obtain high-resolution grid map data for illustrating the spatio-temporal characteristics of climate change,Shiru et al. (2020)  andTian et al. (2021) applied direct spatial interpolation of meteorological station data; however, the impact of topographic features on climate change is still not taken into consideration.For spatial downscaling of the gridded climate dataset, the Delta method with a high spatial resolution reference climate dataset that incorporates topography effects is employed, such that the resulting high spatial resolution gridded map data has accurate geographical location and area information.Consequently, the spatial downscaling dataset obtained in this work comprises topographic features and gives very precise geographical information, such as geographic location, area ratio, and other statistical indicators.Based on this dataset, the results about an examination of the spatial distribution characteristics of climate change in the YRB are more trustworthy.Compared to other studies using Delta spatial downscaling(Wang et Chen., 2014), this study obtained lower MAE values for precipitation and temperature downscaling in the ranges of 23-26 mm and 2.1-2.4 °C, indicating the feasibility of downscaling GCM meteorological data to a 1 km × 1 km grid in the YRB using the Delta method.In addition, five interpolation methods are utilized in the downscaling process of this study, which decreases the errors created by the Delta method in the interpolation process and provides more fair and accurate downscaling results.The BI and Kriging are the most applicable interpolation methods for the downscaling of precipitation and temperature in the research area, and they offer promising application possibilities in the climate study of the YRB.Wang at al. (2021) synthetically evaluated the applicability of 19 climate models to daily precipitation in the YRB by using the rank score (RS) method with eight indicators like MAE and correlation coefficient, and the results showed that the climate models with better simulation results are in the order of the MRI-ESM2-0, ACCESS-CM2, CNRM-CM6-1, CNRM-ESM2-1, FGOALS-f3-L, and MPI-ESMI-2-HR.Peng et al. (2018) investigated the simulation performance of 28 climate models for climate change in the Loess Plateau region of China and determined that the GFDL-ESM2 and NorESM1-M had the best monthly precipitation and temperature simulation performance.The inconsistency between these conclusions and the paper's findings is likely owing to the variable evaluation indicators, temporal resolution, and downscaling or lack thereof of the datasets, all of which might influence the evaluation outcomes.This research uses 39 CMIP6 climate models, five spatial interpolations, Delta downscaling, and four separate assessment metrics to analyze the application of watershed precipitation; the conclusions are exhaustive and trustworthy.
findings ofChun et al. (2022)  indicate that the YRB will experience more precipitation during the carbon-neutral period than during the carbon-peak period(2015)(2016)(2017)(2018)(2019)(2020)(2021)(2022)(2023)(2024)(2025)(2026)(2027)(2028)(2029)(2030), which may mitigate the drought-related disasters caused by global warming-induced increases in evapotranspiration.In light of this, China's ability to reach carbon neutrality by the end of this century is crucial to resolving the water crisis in the YRB.5.ConclusionsThis study combines different interpolation methods with Delta downscaling to downscale the climate models of monthly temperature and precipitation in the YRB, and verifies their simulation accuracy by four different evaluation indicators, and finally, the spatio-temporal characteristics of monthly temperature and precipitation are analyzed according to the MME with the best simulation effect.The Delta downscaling method not only generates more detailed climate information, but also minimizes the uncertainty of the climate model dataset.According to the four indicators, MAE, TS, S and SS, the best fits of each climate model for the monthly mean temperature of the YRB are CESM2-WACCM, NorESM2-LM, and ACCESS-CM2, and the best fits for monthly precipitation are ACCESS-ESM1-5, CESM2-WACCM, and IPSL-CM6A-LR.

Table 1 .
Summary of the 39 climate models of CMIP6