Estimation of the Climate Change Impact on the Hydrological Balance in Basins of South-Central Chile

: In this study, the SWAT (Soil Water Assessment Tool) hydrological model is implemented to determine the effect of climate change on various hydrological components in two basins located in the foothills of the Andes: the Quino and Muco river basins. The water cycle is analyzed by comparing the model results to climatic data observed in the past (1982–2016) to understand its trend behaviors. Then, the variations and geographical distribution of the components of the hydrological cycle were analyzed using the Representative Concentration Pathway (RCP)8.5 climate scenario to model two periods considering the immediate future (2020–2049) and intermediate future (2050–2079). In this way, in the study area, it is predicted that yearly average temperatures will increase up to 1.7 ◦ C and that annual average precipitation will decrease up to 210 mm for the intermediate future. Obtained results show that the analyzed parameters presented the same trend behavior for both periods of time; however, a greater impact can be expected in the intermediate future. According to the spatial distribution, the impact worsens for all the parameters as the elevation increases in both basins. The model depicted that yearly average evapotranspiration would increase around 5.26% and 5.81% for Quino and Muco basins, respectively, due to the large increase in temperature. This may cause, when combined with the precipitation lessening, a decrease around 9.52% and 9.73% of percolation, 2.38% and 1.76% of surface ﬂow, and 7.44% and 8.14% of groundwater for Quino and Muco basins, respectively, with a consequent decrease of the water yield in 5.25% and 4.98% in the aforementioned watersheds, respectively.


Introduction
Climate change and the effect of anthropic alterations on the hydrological cycle have achieved considerable interest in the scientific community. Several impacts of climate change can be predicted through the variation in the hydrological cycle [1,2]. Strong evidence showing that climate change is altering hydrological cycles at the regional and global levels has been exposed [3,4], indicating consistent predicted effects: changes in precipitation patterns, increased intensity of extreme weather events, glacier retreat and its consequent alterations in the river discharge regimes [4].
The study objectives are first to analyze the annual variation of the hydrometeorological components over the last 35 years. Furthermore, the prediction of the response of the hydrological cycle according to the immediate future (2020-2049) and intermediate future  climate scenarios is assessed.

Study Areas
Quino and Muco basins are located in the Araucanía Region, in the south-central zone of Chile (38 • 10 00 S and 38 • 40 00 S) (Figure 1). This zone supports an important industrial development strongly related to forestry and agricultural production, also linked to hydroelectric generation. It is located on the border of two significant eco-regions for world conservation, the Mediterranean Scrub and the temperate rainforest, within the biodiversity hotspot called Valdivian Rainforest [25].
climatic variations over the last decades. 101 The study objectives are first to analyze the annual variation of the hy-102 dro-meteorological components over the last 35 years. Furthermore, the prediction of the 103 response of the hydrological cycle according to the immediate future (2020-2049) and 104 intermediate future  climate scenarios is assessed. 105

Study Areas 107
Quino and Muco basins are located in the Araucanía Region, in the south-central 108 zone of Chile (38°10'00" S and 38°40'00" S) ( Figure 1). This zone supports an important 109 industrial development strongly related to forestry and agricultural production, also 110 linked to hydroelectric generation. It is located on the border of two significant 111 eco-regions for world conservation, the Mediterranean Scrub and the temperate rainfor-112 est, within the biodiversity hotspot called Valdivian Rainforest [25].
113 Selected pre-cordilleran hydrographic systems, Quino and Muco basins, have 299 114 and 651 km 2 of area and elevations ranging from 305 to 1724 mamsl and from 189 to 1469 115 mamsl, respectively. The basins are located in a transition area between Mediterranean 116 dominance and temperate humid conditions [26]. Annual average precipitations in the 117 basins vary between 1253 mm and 2693 mm and annual average temperatures between 118 10 °C and 12 °C. Selected pre-cordilleran hydrographic systems, Quino and Muco basins, have 299 and 651 km 2 of area and elevations ranging from 305 to 1724 mamsl and from 189 to 1469 mamsl, respectively. The basins are located in a transition area between Mediterranean dominance and temperate humid conditions [26]. Annual average precipitations in the basins vary between 1253 mm and 2693 mm and annual average temperatures between 10 • C and 12 • C.

Hydrological Modelling 2.2.1. Model Description
The SWAT model (Texas A&M University, Texas, USA, Version 2012. 10.4.21) allows simulating several physical processes concerning the hydrological cycle at different time scales. The model presents three different ways to discretize the basin: grids, slopes and sub-basins schemes [27]. In this work, the monthly time scale and the sub-basin scheme derived from the digital elevation model (DEM) were approached. Thereby, sub-basin configuration considers the preservation of the channels and the natural flow path within the model. Firstly, the waterways were defined, opting to create them from the minimum area threshold automatically established by SWAT, achieving an accurate representation of the channels.
Next, a hydrologic response units (HRUs) discretization method was also applied in order to improve the simulation accuracy. HRUs were obtained by integrating the soil type, land use and slope. This method allows reflecting more accurately the spatial variability of the process through which rain transforms into runoff, as well as the routing of the water in each HRU [28].
Additionally, daily climatic information on precipitation and maximum and minimum temperatures was incorporated to calculate the main processes of the hydrological cycle, using the Hargreaves method [29] to determine the potential evapotranspiration.

Databases
The input data for SWAT hydrological model were obtained using the Geographic Information System ArcGIS, ArcMap 10.4. (ESRI Geographic Information System Company, California, USA) The DEM was developed using Alos-1 Palsar sensor images with a spatial resolution of 12.5 m (https://vertex.daac.asf.alaska.edu/, accessed on 12 November 2020). Soil type information was obtained from the agrological studies of the IX Region performed by the Natural Resources Information Center (CIREN) in 2002 for the Araucanía region [30]. The aforementioned study provides information about the soil types present on the study area and their associated physical-chemical properties. This information was included in the database used in the SWAT model for the hydrological simulations [28]. Land use and coverage data for Muco and Quino basins for 2011 were acquired from satellite image processing carried out by Heilmayr et al. [31].
On the one hand, meteorological data used to analyze the climatic variability in the past were extracted from the gridded product CR2MET developed by the Center for Climate and Resilience Research (CR) 2 (http://www.cr2.cl/, accessed on 12 November 2020) [32]. This product is based on actual observations and atmospheric reanalysis on a 0.05 • latitude-longitude grid (~5 km spatial resolution) for the continental Chilean territory. In this study, daily meteorological information (precipitation and extreme temperatures) was used for a 35-year time period .
On the other hand, future climate information was obtained from the (CR) 2 climate data generator [33] at a spatial resolution of 10 km (http://www.cr2.cl/, accessed on 12 November 2020).The GCMs from the Fifth Phase of the Coupled Model Intercomparison Project (CMIP5) report of the Intergovernmental Panel on Climate Change (IPCC) were evaluated and categorized according to several performance indices to find the global models correctly representing the climate in the southeastern Pacific region [33]. Based on this analysis, and considering the data availability, the global model MPI-ESM-MR (Max Planck Institute for Meteorology, Hamburg, Germany) was used together with the Regional Climate Model System version 4 (RegCM4, National Center for Atmospheric Research, Colorado, USA.). In this way, historical data from 1976-2005 and a period of climate simulations from 2007-2100 were used. It is important to highlight that the database used in this work as supplied by (CR) 2 is the only available climatic product for the study area and that there are two Representative Concentration Pathways (RCP) simulated scenarios: RCP2.6 and RCP8.5. In the case of the RCP2.6 scenario, it could only be possible if greenhouse gas emissions become extremely controlled [4]. Even while the climate system is in a constant change state, achieving RCP2.6 would certainly present the smallest impacts on natural ecosystems and human activities [11]. To reach such radiation levels, greenhouse gas emissions should be significantly reduced [34]. Moreover, the other scenarios consider the increase of greenhouse gas emissions and radiative forcing in different tendencies, with RCP8.5 the most critical scenario [11]. According to the aforementioned reasons, the RCP8.5 climate scenario from the local product RegCM4-MPI-ESM-MR ((CR) 2 , Santiago, Chile) was used in this work [33].
It is worth noting that the available data from CR2MET (CR) 2 , Santiago, Chile) database for past climate and data from RegCM4-MPI-ESM-MR used for future climate were extracted at the same geographical points (latitude-longitude) (Figure 2d). Point extraction was carried out to avoid possible biases induced by changes in the spatial positioning of the meteorological data provided to the SWAT model between calibration and projection. The observational data from CR2MET and the historical data of RegCM4-MPI-ESM-MR climate product have been validated by (CR) 2 in its report "Regional Climate Simulations" [33]. A validation procedure was performed to evaluate the capacity of the RegCM4 to simulate the current climate. According to their report [33], excellent results were obtained in the model behavior for simulating precipitation and extreme temperatures in the south-central zone of Chile.  Calibration, sensitivity and uncertainty analysis were performed by using the SU-210 FI-2 algorithm (Sequential Uncertainty Fitting Version 2) [35,36] included in SWAT-CUP 211 (SWAT Calibration, Uncertainty Procedures) [37]. This is a semi-automated reverse 212 modeling procedure that combines and optimizes the modeling parameters for an accu-213 rate result [35]. 214 Abbaspour et al. [35] have suggested that the sensitivity analysis can be determined 215 through the Latin Hypercube SUFI-2 sampling, performing "n" parameter combinations. 216 In the present investigation, the number of combinations was set at 500 simulations. To 217 Additionally, observed flow data were obtained from (CR) 2 at the closure point of both basins, corresponding to the fluviometric stations Rio Quino en Longitudinal and Rio Muco en Puente Muco (Figure 2d). In order to select the gauging points, a consistency analysis was performed. For this purpose, the quality, coherence and the time extension of the databases were evaluated using the fluviometric data from (CR) 2 [33]. Both hydrometric stations have 95% of the information related to daily average flow. For the Quino basin, we obtained information from 1980-2012 and, in the case of Muco basin, from 1980-2016.
Information of the cover/use, soil types, slope, points of climatic data and flows for Quino and Muco basin are represented in Figure 2.

Calibration, Sensitivity and Uncertainty Analysis
Calibration, sensitivity and uncertainty analysis were performed by using the SUFI-2 algorithm (Sequential Uncertainty Fitting Version 2) [35,36] included in SWAT-CUP (SWAT Calibration, Uncertainty Procedures) [37]. This is a semi-automated reverse modeling procedure that combines and optimizes the modeling parameters for an accurate result [35].
Abbaspour et al. [35] have suggested that the sensitivity analysis can be determined through the Latin Hypercube SUFI-2 sampling, performing "n" parameter combinations. In the present investigation, the number of combinations was set at 500 simulations. To start evaluating the simulations, SUFI-2 determines the relative significance of each parameter (t test), as a measure of their sensitivity [35]. In this study, the global sensitivity method of the SUFI-2 algorithm was implemented.
Calibration and sensitivity analysis were initiated using 24 parameters where 9 of them were selected according to the highest sensitivity to perform the calibration ( Table 1). The definition of the parameters to be calibrated and their ranges were performed considering studies from several reports [28,37]. In the first instance, the calibration was performed on a monthly scale considering the measured flows at the closing point of Quino and Muco river basin corresponding to the fluviometric stations "Rio Quino en Longitudinal" and "Rio Muco in Puente Muco", respectively, using the land cover/use corresponding to 1986.The model was run considering three warm-up years (1979)(1980)(1981) for the calibration period from 1982-1992. Validation for both basins was determined applying the same adjusted values used in the calibration but considering a different set of observed data corresponding to 2011, in order to ensure that the model has a satisfactory precision range within its applicability domain.
Then, the validation procedure for Quino basin was conducted using a fluviometric data set from 2000-2013; furthermore, a data set from 2006-2016 was used for the Muco basin model validation. Modeled flow representation after calibration and validation procedures were evaluated using statistical tests such as the efficiency Nash-Sutcliffe (NSE) index, percentage bias (PBIAS) and determination coefficient (R 2 ) criteria, according to Moriasi et al. [38].
Finally, after being calibrated and validated with CR2MET data, the model was executed again, with information from the climatic model RegCM4-MPI-ESM-MR in the

Parameterization
The parameterization was carried out to apply the adjustments suggested in the calibration procedure for the study basins. The global modification terms relative (r) and replace (v), were considered based on soil type, use and cover, location, slope, subbasin number, soil layer, combinations of related parameters and others as suggested elsewhere [21,39]. Table 1 shows the global modification terms applied to the most sensitive parameters for Quino and Muco basins parameterization.

Climatic Variability Analysis
Climatic variability analysis was conducted once the hydrological model was calibrated and validated for both basins also with the CR2MET product. In this way, meteorological elements and hydrological components such as precipitation, temperature, ET (evapotranspiration), PERC (percolation), SURQ (surface runoff), GW-Q (underground flow) and WYLD (discharge) were analyzed. A Mann-Kendall (MK) test [40,41] was applied in order to determine the parameter trends in the last 35 years for the Quino and Muco basins.
The MK test is a non-parametric method [42] widely used in several hydroclimatic studies [43][44][45][46]. The advantage of the MK test lies in the fact that trend analysis is almost not influenced by a small number of outliers for sequence analysis. In this study, the significance level of the test was set as 0.05.
The magnitude of predicted changes in the trend of the meteorological variables was estimated using the Sen slope method [47,48]. In this technique, the gradient was calculated as a change in the measurements, correlated to time change units. The Theil-Sen estimator provides additional trend information on the variation intensity [47,49].

Climatic Change Analysis
In order to perform climatic studies, the SWAT model allows incorporating weather anomalies to simulate climate change as well as the prediction of future weather patterns [27,50], and in this way the modifications of the water balance components can be calculated.
The basins' hydrological responses were simulated using the local climate model RegCM4-MPI-ESM-MR [33]. The future climate projection RCP8.5 was selected to consider the most critical scenario in CMIP5.
Yearly averages of the water cycle parameters (ET, PERC, SURQ, GW_Q and WYLD) were obtained considering periods of thirty years. In such analyses, the "historical period" (1976-2005) from the climate model was set as the baseline while the period from 2020-2049 was considered "immediate future" and 2050-2079 as "intermediate future".
Both future scenarios were compared to baseline using Equations (1) and (2), respectively, at a sub-basin scale. In this way, future changes in the absolute values of the hydrological parameters predicted could be spatially analyzed, improving the comprehension of the water cycle within the basin for the next decades under extreme climate scenarios. Future relative changes at a basin scale were also reported in this study. For such analyses, projected scenarios were compared to historical periods according to Equations (3) and (4) The analysis was carried out on an annual scale comparing projected scenarios to historical periods in order to quantify the expected impact in the next decades in Quino and Muco basins. Additionally, a Student's t test for related samples was performed to determine the existence of significant differences for the projected climatic data and the projected hydrological cycle components.
In order to determine the existence of significant differences for every component of the water cycle among simulation results for the historical, immediate and intermediate periods, a Student's t test statistic test for paired samples was applied. The test was performed once it was checked that every sample could be properly fitted to a normal distribution via the Kolmogorov-Smirnov test. Additionally, in order to comply with the test requirements, all the simulations were performed under similar conditions but changing the climatic scenarios. In this way, the impact of climatic change over the water cycle components distribution was unveiled. This analysis was conducted on an annual scale, at the sub-basin level, for a 30-year period with a 95% significance level.

Model Calibration and Validation
According to the hydrological network, Quino and Muco basins were divided into 99 and 91 sub-basins with average surfaces of 3.02 km 2 and 7.16 km 2 , respectively. In this way, the surface range of the sub-basins oscillates between 0.006 km 2 and 14.33 km 2 in Quino and between 0.01 km 2 and 45.93 km 2 in the Muco basin. The SWAT model was run to obtain the water discharge in all the sub-basins for both watersheds.
The model performance was calculated in the lowest basin points corresponding to the fluviometric stations Rio Quino en Longitudinal and Rio Muco en Puente Muco ( Figure 2). In the Quino basin, the determination coefficient and the Nash-Sutcliffe efficiency parameter presented performance classification of "very good", showing R 2 values of 0.88 and 0.89 and NSE 0.84 and 0.79 in the calibration and validation, respectively, on a monthly scale. The percentage of bias in the calibration maintained a "very good" index (−3.11%), and the statistic obtained in the validation (−20.68) was classified as satisfactory by Moriasi et al. [38] (Figure 3a). On the other hand, in the Muco basin, statistics obtained showed "very good" classification for R 2 (0.88 and 0.92) and NSE (0.88 and 0.89) in the calibration and validation process, respectively. The PBIAS showed a "very good" performance in the calibration phase (5.92%) and "good" in the validation phase (−11.85%) (Figure 3b). Figure 3 depicts the performance of the discharge simulations during the calibration and validation periods for both basins. Performance values recorded "very good" to "satisfactory" grades, according to the classification proposed by Moriasi et al. [38], considering monthly scale verification. ues of 0.88 and 0.89 and NSE 0.84 and 0.79 in the calibration and validation, respectively, 314 on a monthly scale. The percentage of bias in the calibration maintained a "very good" 315 index (−3.11%), and the statistic obtained in the validation (−20.68) was classified as sat-316 isfactory by Moriasi et al. [38] (Figure 3a). On the other hand, in the Muco basin, statistics 317 obtained showed "very good" classification for R 2 (0.88 and 0.92) and NSE (0.88 and 0.89) 318 in the calibration and validation process, respectively. The PBIAS showed a "very good" 319 performance in the calibration phase (5.92%) and "good" in the validation phase 320 (−11.85%) (Figure 3b

Annual Variation of the Hydrometeorological Components in the Last 35 Years
According to the results reported in Figure 4, the yearly variation trend of the hydrological components is directly related to the meteorological parameters ( Figure 4). In order to evaluate such changes in the climate and hydrology of the study basins, 35  The precipitation in the study basins presented a significant decreasing trend during the period, with p-value of 0.002 (Quino) and 0.005 (Muco) (Figure 4b,i). The mean annual precipitation value of the period for the Quino basin was 1811 mm, while for the Muco basin it was 1781 mm. Such a decreasing trend in precipitation, in addition to the temperature increase during these years, could cause a decrease in hydrological components such as ET, PERC, SURQ, GW_Q and WYLD, as shown in the model results.
Additionally, during the 35-year period, the evapotranspiration presented an annual average of 643 mm for Quino and values around 646 mm for Muco basin (Figure 4c,j). This parameter does not show a significant trend in its behavior; however, the remaining hydrological components presented significant trends according to the Mann-Kendall test.
Yearly mean values obtained for PERC, SURQ, GW_Q and WYLD were 644 mm, 377 mm, 599 mm and 1111 mm, respectively, for Quino basin (Figure 4d,e,f,g). Meanwhile, for Muco basin the annual averages for the period were 626 mm (PERC), 378 mm (SURQ), 590 mm (GW_Q) and 1081 mm (WYLD) (Figure 4k,l,m,n), showing similar behavior in both basins.  (Figure 5a). A significant and increased temperature rise can be observed as spatially represented from downstream to upper stream both in the Quino and Muco basins. Furthermore, predicted precipitation changes (Figure 5b) indicate a significant decrease for both study areas with a p-value <0.0001 for immediate and intermediate future climatic periods. As shown, the higher changes in precipitation for both basins are expected to occur upstream. Raw p-values of the tests for input data and water cycle parameters modeled are reported as information in Appendix A, Table A1. Additionally, during the 35-year period, the evapotranspiration presented an annual 349 average of 643 mm for Quino and values around 646 mm for Muco basin (Figure 4c,j). 350 This parameter does not show a significant trend in its behavior; however, the remaining 351 hydrological components presented significant trends according to the Mann-Kendall 352 test. Yearly mean values obtained for PERC, SURQ, GW_Q and WYLD were 644mm, 377 353 mm, 599 mm and 1111 mm, respectively, for Quino basin (Figure 4d,e,f,g). Meanwhile, 354 for Muco basin the annual averages for the period were 626 mm (PERC), 378 mm 355 (SURQ), 590 mm (GW_Q) and 1081 mm (WYLD) (Figure 4k,l,m,n), showing similar be-356 havior in both basins. 357 The aforementioned increase in temperature and a decrease in precipitation could be even expected for both climatic periods. The greatest changes are predicted to occur in the intermediate future, with an increase up to 1.7 • C in the mean yearly temperature and an average yearly precipitation decrease up to 200 mm (Figure 5a,b). Such trends agree with the behavior observed during the historical period.

Evapotranspiration (ET)
The geographic distribution of the projected changes in evapotranspiration is represented in Figure 6. According to the Student's t test, for both periods significant differences exist with p-values <0.0001 for both basins. The higher increase in evapotranspiration is expected to occur in the mountain area of both basins.

Evapotranspiration (ET) 380
The geographic distribution of the projected changes in evapotranspiration is rep-381 resented in Figure 6.According to the Student's t test, for both periods significant differ-382 ences exist with p-values <0.0001 for both basins. The higher increase in evapotranspira-383 tion is expected to occur in the mountain area of both basins.

Percolation (PERC)
Similarly, percolation projections in the study basins depict significant decreasing trends during the climatic periods analyzed. Results represented in Figure 7 show the geographical spread of the decreasing trend in this water cycle parameter for both basins during the next decades, with the highest changes in the mountain area.  The decreasing trend is consistent in both river basins, indicating an intensification of the trend observed since the historical period analysis (Figure 7).

Surface Runoff (SURQ)
Following the performed analysis, the surface runoff shows a decreasing trend in both basins during the projected climatic periods, with a greater decrease in the intermediate future ( Figure 8). According to Student's t test, significant differences could be predicted with p-values lower than 0.0001. The variation of this water cycle parameter during the immediate future seems to be limited downstream, changing less than 10 mm. However it would be geographically balanced in the middle part of the basins where the largest number of sub-basins with projected annual variation reaching 59 mm can be found. However, in the intermediate future, the Quino basin does not present a clearly homogeneous correlation to height. Meanwhile, in the Muco basin, a marked decreasing variation is observed as the height increases in almost all its sub-basins.

Surface Runoff (SURQ) 407
Following the performed analysis, the surface runoff shows a decreasing trend in 408 both basins during the projected climatic periods, with a greater decrease in the inter-409 mediate future (Figure 8). According to Student's t test, significant differences could be 410 predicted with p-values lower than 0.0001. The variation of this water cycle parameter 411 during the immediate future seems to be limited downstream, changing less than 10 mm.  Decreasing trends could also be predicted during both periods studied, with signif-425 icant values for Quino and Muco basins for the groundwater flow. If analyzed as spa-426 tially distributed, a generalized change strongly influenced by the sub-basins altitude can 427 be observed in Figure 9.

431
The groundwater flow projection worsens in the intermediate future with maximum 432 variation values obtained up to 204mm. This downward trend seems to be consistent 433 with the behavior that occurred during the historical period. 434 For both climatic periods, the statistical tests highlight a generalized and significant decrease in the surface runoff; such a trend also agrees with the changes observed during the historical period.

Groundwater Flow (GW_Q)
Decreasing trends could also be predicted during both periods studied, with significant values for Quino and Muco basins for the groundwater flow. If analyzed as spatially distributed, a generalized change strongly influenced by the sub-basins altitude can be observed in Figure 9. Following the performed analysis, the surface runoff shows a decreasing trend in 408 both basins during the projected climatic periods, with a greater decrease in the inter-409 mediate future (Figure 8). According to Student's t test, significant differences could be 410 predicted with p-values lower than 0.0001. The variation of this water cycle parameter 411 during the immediate future seems to be limited downstream, changing less than 10 mm. Decreasing trends could also be predicted during both periods studied, with signif-425 icant values for Quino and Muco basins for the groundwater flow. If analyzed as spa-426 tially distributed, a generalized change strongly influenced by the sub-basins altitude can 427 be observed in Figure 9.

431
The groundwater flow projection worsens in the intermediate future with maximum 432 variation values obtained up to 204mm. This downward trend seems to be consistent 433 with the behavior that occurred during the historical period. 434  Finally, according to the Student's t test, the projected water discharge suggests a significant decrease in the coming decades for both study basins ( Figure 10) in all the sub-basins. It can be spatially distinguished that, as in the previous parameters, there is a marked relationship between the predicted changes and the sub-basin's altitude, with a decrease in projected discharge greater as the elevation increases for both hydrographic basins.
Finally, according to the Student's t test, the projected water discharge suggests a 436 significant decrease in the coming decades for both study basins ( Figure 10) in all the 437 sub-basins. It can be spatially distinguished that, as in the previous parameters, there is a 438 marked relationship between the predicted changes and the sub-basin's altitude, with a 439 decrease in projected discharge greater as the elevation increases for both hydrographic 440 basins. 441 It is worth noting that the trends obtained for both climatic periods seem homoge-442 neous, showing a discharge decrease. However, this tendency seems to intensify during 443 the intermediate future if compared to the immediate future ( Figure 10   It is worth noting that the trends obtained for both climatic periods seem homogeneous, showing a discharge decrease. However, this tendency seems to intensify during the intermediate future if compared to the immediate future ( Figure 10).
The yearly average decrease in discharge for the intermediate future in the lower and middle areas of the basins could reach 80 mm to 200 mm, while in the mountainous area a drop of up to 304 mm could be expected. This significant reduction projected in the Quino and Muco basins for each of their sub-basins is also consistent with the trend obtained in the historical period ( Figure 10). Table 2 shows the yearly relative and absolute changes in Quino and Muco basins for both periods modeled at a basin scale. As it was expected from previous sections, higher average change percentage could be expected in almost every parameter for the intermediate future if compared to the immediate future.

Relative and Absolute Changes of the Hydrological Balance Components at a Basin Scale
In both climatic periods, temperature is the parameter with the highest relative change both for Quino and for Muco basin, projecting a maximum increase of 14.78% and 14.16%, respectively. Ordered from highest to lowest change are the PERC, GW_Q, ET and WYLD parameters, while the SURQ presented the lowest relative change with a decrease of 0.72% and 0.58% (immediate future and intermediate future) for Quino basin and 1.31% and 0.84% (immediate future and intermediate future) for Muco basin.

Discussion
It could be expected that, by incorporating the precipitation and temperature anomaly from a future climate scenario (RCP 8.5), the remaining water cycle components become modified, since precipitation is the main component of change in the hydrological cycle and basically confers it mobility [51][52][53]. On the other hand, temperature is one of the components governing evapotranspiration [54], and the transpiration fraction is the process causing greatest water loss in the basin system [55]. This is why, when these drivers become altered, the components of the water balance are also modified.
These temperature and rainfall changes both in the Quino and Muco basins seems to be induced from moderate to severe changes in the regime of the hydrological cycle components evaluated. Obtained results imply that the average annual temperature increase up to 1.5 • C, representing a percentage change of 14.78% and 14.16% for Quino and Muco, respectively, leads to a decrease in PERC (−9.52 mm and −9.73 mm), SURQ (−2.38 mm and −1.76 mm), GW_Q (−7.44 mm and -8.14 mm) and WYLD (−5.25 mm and −4.98 mm) and an increase in ET (+5.26 mm and +5.81 mm) for Quino and Muco basins, respectively, in the intermediate future scenario. The potential loss of water resources through the deficit of rainfall and greater evaporation has been quantified in different studies, reaffirming that this geographical area of Chile is particularly vulnerable to climate change [16,[56][57][58][59].
On the other hand, in both study basins, on a spatial scale, it can be observed that the annual average changes for both study periods (immediate future and intermediate future) present the same behavioral trend following a staggered pattern. It is relevant that the greatest changes in ET, PERC, SURQ, GW_Q and WYLD take place in the mountainous zone (pre-mountain range) and minor changes in the alluvial plain (central depression). The behavior of these parameters is directly related to the temperature increase behavior between 0.5 • C and 1.7 • C from the lower part to the upper part of the basin, respectively, and a decrease in precipitation of 10 mm to 200 mm between the two time periods for these geographical areas.
A similar behavior, different in geographical area, was also reported by Falvey and Garreaud [60]. They observed a cooling of coastal temperatures and a gradual warming over the Andes mountain range for the central and northern zone of Chile. They reported a decrease in temperatures (−0.15 • /decade) along the coast, slight increases both in the central valley and the mountain range stations, with a significant increase of almost +0.25 • C per decade. The orographic effect induced by the Andes Mountains on precipitation defines a remarkable character of the hydrological regime in the central-southern region of the country, where precipitation on the slopes of the mountains can be from two to three times greater than in the lower elevations [59,61]. This effect, in addition to the climatic projections, could cause a greater change for the windward Andes slope.
Several studies suggest variable hydrological responses of the water cycle parameter as the effect of future climate change scenarios around the world.  2009) indicated that increasing 2.8 • C in temperature causes a change in real evapotranspiration and water production of +6.2% and −13.9%, respectively. Andreini et al. [62] concluded that the large changes that occur in surface runoff correspond to a lesser change in precipitation. Additionally, in India, Pandey et al. [63] quantified the climate change impact (2071-2100) in the Armur basin using SWAT, suggesting that evapotranspiration and water production will increase 28% and 49%, respectively.
Grusson et al. [19] also compared different climatic ensembles in the south zone of France during 2010-2050. In their study, a strong impact on green water was shown, specifically a reduction in soil water content (SWC) and a substantial increase in evapotranspiration (ET) in winter. In summer, however, the ET faced lower flows due to the lack of SWC, which highlights future deficits of green water reserves. Besides, Yu et al. [23] showed an increase in annual runoff depth at the northeast zone of China by 18.1% (RCP 2.6), 11.8% (RCP 4.5), 23.6% (RCP 6.0) and 11.5% (RCP 8.5), compared to the base years. Therefore, worldwide, differentiated behavior can be observed depending on the physical-geographical conditions of each study area.
Specifically in Chile, Stehr et al. [24] modeled the hydrological response to climate change in two sub-basins of the Biobío River, depicting that a 30% decrease in precipitation could conduct a 45% and 32% reduction in the water flow of the study sub-basins. In addition, Barrientos et al. [64] recently reported an study also conducted in experimental sub-basins of the BioBio river obtaining similar results in the trends. Even when these forecasts suggest more dramatic reductions in precipitation and runoff than the present study, the observed trends become similar. Generally, as it has been reported by the (CR) 2 group studies [33] and as the aforementioned studies suggest, the warming effects and consequences seem to be greater for northern than for southern zones in Chile [33].
Finally, the present study showed that future climate projections maintain the same behavioral trend that was evidenced in the simulated past. This similarity shows the quality of the data used for the study and the correct performance of the calibration and validation for both study basins (Quino and Muco). In this way, it can be stated that the performance metrics calculated after the configuration of the SWAT model by the calibration and validation procedures are functional for long-term hydro-climatic evaluation in the south-central zone of Chile.
However, as temperatures rise and precipitation decreases, the land cover/use is unlikely to remain unchanged until 2079 even without human intervention. It is necessary to develop studies addressing the individual and combined effects of the land use/cover change and future climate change in basins of the south-central zone of Chile. The importance of treating the data in this way is to determine which of the alterations can affect the hydrological components to a greater and lesser extent in order to determine their individual and combined effect. This could help on the one hand understanding the effect of these two stressors on the components of the hydrological cycle and, on the other hand, to address most probable future scenarios. The present manuscript provides the first results to determine the importance of the individual effect of climate change in the south-central zone of Chile. Future work by this group will intend to study the problems described above.

Conclusions
This study determined the possible effects of future climate on the water balance of the Quino and Muco basins for land use/cover by using the SWAT hydrological model. The hydrological response of the basin was simulated using the local climate model RegCM4-MPI-ESM-MR, for the future climate projection RCP 8.5, as it is the most critical scenario in CMIP5. The future projection was divided into two time periods (immediate future and intermediate future). Calibration and sensitivity analyses were carried out to ensure that the hydrological model output was reliable. The results show for both analysis periods an increase in temperature and a decrease in precipitation. The analyzed parameters show the same trend behavior for both time periods; however, a greater impact in the intermediate future could be expected.
For the intermediate future, on an annual average, temperatures are expected to increase up to 1.7 • C and rainfall will decrease by around 210 mm. As a consequence, as temperature rises, ET would increase causing, combined with the precipitation effect, a decrease in PERC, SURQ, GW_Q and WYLD. Additionally, the study results suggest that the expected changes for both periods of time would present the same behavior that has occurred since the 1980s, with a greater effect on the intermediate future.
The present investigation results could be useful for future studies in areas such as native forest cover, forest plantation and agriculture. Furthermore, information is provided to investigate how future climate changes could affect hydrological processes in watersheds located in the central south zone of Chile. The data and results shown here could help in the integrated and sustainable basin management within the region.