Using SWAT and Fuzzy TOPSIS to Assess the Impact of Climate Change in the Headwaters of the Segura River Basin ( SE Spain )

The Segura River Basin is one of the most water-stressed basins in Mediterranean Europe. If we add to the actual situation that most climate change projections forecast important decreases in water resource availability in the Mediterranean region, the situation will become totally unsustainable. This study assessed the impact of climate change in the headwaters of the Segura River Basin using the Soil and Water Assessment Tool (SWAT) with bias-corrected precipitation and temperature data from two Regional Climate Models (RCMs) for the medium term (2041–2070) and the long term (2071–2100) under two emission scenarios (RCP4.5 and RCP8.5). Bias correction was performed using the distribution mapping approach. The fuzzy TOPSIS technique was applied to rank a set of nine GCM-RCM combinations, choosing the climate models with a higher relative closeness. The study results show that the SWAT performed satisfactorily for both calibration (NSE = 0.80) and validation (NSE = 0.77) periods. Comparing the long-term and baseline (1971–2000) periods, precipitation showed a negative trend between 6% and 32%, whereas projected annual mean temperatures demonstrated an estimated increase of 1.5–3.3 ◦C. Water resources were estimated to experience a decrease of 2%–54%. These findings provide local water management authorities with very useful information in the face of climate change.


Introduction
Climate change as a result of increased greenhouse gas emissions leads to changes in hydrologic conditions and results in various impacts on the availability of global water resources [1].Spain is one of the countries most vulnerable to the impacts of climate change in Europe due to the high spatial and temporal irregularity of water resources and socio-economic characteristics [2].In addition, future climate tendencies show an increment in the temperature and a significant reduction in total annual rainfall [3].When these impacts occur in regions that already present low water resource availability and frequent droughts, these impacts can be exacerbated.The Segura River Basin (SRB) is situated in SE Spain and is one of the most water-stressed basins in Mediterranean Europe [4].
In global terms, the total water demand for consumption is 1800 hm 3 /year, where 86% corresponds to agricultural use and 10% to urban uses.Against this demand data, annual natural water availability is, on average, around 800 hm 3 /year [5].This deficit is partly covered by water from the Tagus-Segura water transfer and the use of unconventional water resources like treated wastewater and desalinated water, but these resources are not enough and the SRB is still suffering aquifer overexploitation [6].The headwaters of the Segura River Basin (HWSRB) need to be studied thoroughly, as they are the most important sites for water resource generation in the basin [7].The HWSRB have an important relevance in SRB water resources, since they comprise 9% of the water resource contribution, in spite of the fact that they covers only 1.2% of the area of the total watershed.
The SWAT has been successfully and widely used all over the world for different purposes, including the evaluation of climate change impacts on water resources [8].However, SWAT applications assessing the water resources of Spain under changing climate conditions are scarce in scientific literature.Such studies have mostly been conducted in the north of the country, where there is an absence of water scarcity [9][10][11][12][13][14].In the case of Spanish Mediterranean catchments, this model has rarely been used [15].
The Technique for Order Preference by Similarity to an Ideal Solution (TOPSIS) is one of the most used techniques for solving Multi-Criteria Decision Analysis (MCDA) problems and was first developed by Hwang and Yoon [16].In the classic formulation of the TOPSIS method, personal judgements are represented with crisp values.However, crisp data are inadequate to model real-life decision problems under many conditions.That is why the fuzzy TOPSIS method was proposed, whereby the weights of criteria and ratings of alternatives are evaluated by linguistic variables represented by fuzzy numbers to deal with the deficiency in the classic TOPSIS.The fuzzy TOPSIS method has been widely applied in many fields; for example, energy [17], environment [18], industrial processes [19], and climate change [20].However, to the best of our knowledge, the only precedent in the combined use of the SWAT model and the fuzzy TOPSIS method is found in Won et al. [21], wherein the authors assessed the water use vulnerability in 12 basins of South Korea, using SWAT to simulate hydrological components and fuzzy TOPSIS to rank the water use vulnerability in those basins.
The aim of the present study was to evaluate the climate change effect on the water resources of the HWSRB.The specific objectives of this study included: (1) to reduce the uncertainty in climate change projections by applying the fuzzy TOPSIS technique to rank climate models and (2) to explore the water resource response to future climate projections for the HWSRB.To achieve these objectives, we set up a hydrological model, using SWAT for the HWSRB.After calibration and validation of the model, nine different climate models were downloaded from the EURO-CORDEX initiative [22] and the fuzzy TOPSIS technique was applied to select which historical runs had the best fit with the observed climate data during a baseline period .Once those climate models were ranked, some of them were used to evaluate climate change in the study area for the medium term (2041-2070) and long term (2071-2100) using two different representative concentration pathways (RCP4.5 and RCP8.5).The results obtained in this study provide local water management authorities with very useful information for the proper utilisation and management of water resources under climate change conditions in this vulnerable region.

Description of the Study Area
The HWSRB is located in the southeastern region of Spain and covers an area of about 235 km 2 .The basin is characterised by steep terrain, and the elevation ranges between 898 and 1912 m, as is shown in Figure 1.The drainage network is formed by two main rivers, the Segura River and its tributary, the Madera River.The whole basin drains into the Anchuricas Reservoir.The mean discharge at the reservoir is 1.6 m 3 /s.This reservoir was constructed in 1957; it has a capacity of 6 hm 3 , and the key purpose of the reservoir is to generate electricity.The study area is located mostly on permeable outcrops, limestone, and dolomite of the Upper Cretaceous.Land use in the watershed is highly forested; 61% of the surface is occupied by forests and 19% by Mediterranean shrubland vegetation.The rest of the land use is mainly for range purposes.The predominant soil type is rendzic leptosol, with variable depth always less than 50 cm, good drainage, and abundant stoniness [23].
The climate is typical Mediterranean with clear seasonality, rainy springs and autumns, and dry summers.According to data from 1971 to 2000, the mean annual precipitation ranged from 511 to 1300 mm, with an average value of 878 mm for the HWSRB, and the average annual temperature was 12.4 • C.

Methodology
In order to evaluate the future impacts of climate change on water resources, the hydrological cycle was simulated in the headwaters of the SRB under different climate change scenarios.This process included three main steps: (1)

SWAT Model
The SWAT [24] is a hydrological watershed model to evaluate the land practice water, sediment transport, and agricultural chemical yields in complex watersheds where soils, land use, or management can widely change.
The SWAT is a semi-distributed and physically based model.The balance equation used is where SW is the final water soil content, SW is the initial water soil content, R is the precipitation, Q is the surface runoff, E is the evapotranspiration, W is the percolation, and Q is the amount of baseflow (all in mm).
The basin is divided into sub-basins and those, in turn, are divided into hydrologic response units (HRUs).HRUs are defined by homogeneous regions with the same slope, soil, and land use.Each HRU generates an amount of runoff that is routed to calculate the total runoff.In order to calculate HRUs, the slope was divided into three classes (0%-8%, 8%-30% and >30%) and a threshold level of 10% was established to facilitate model processing and eliminate minor soils, slopes, and land uses for each subbasin.

Input Data for Hydrological Modelling
As for the data used to carry out the hydrological modelling, catchments were defined based on the digital elevation model (DEM), available on the website of the National Center for Geographic Information [25], with an accuracy of 25 m × 25 m.Meteorological data were obtained from the

SWAT Model
The SWAT [24] is a hydrological watershed model to evaluate the land practice water, sediment transport, and agricultural chemical yields in complex watersheds where soils, land use, or management can widely change.
The SWAT is a semi-distributed and physically based model.The balance equation used is where SW t is the final water soil content, SW O is the initial water soil content, R day is the precipitation, Q surf is the surface runoff, E a is the evapotranspiration, W seep is the percolation, and Q gw is the amount of baseflow (all in mm).
The basin is divided into sub-basins and those, in turn, are divided into hydrologic response units (HRUs).HRUs are defined by homogeneous regions with the same slope, soil, and land use.Each HRU generates an amount of runoff that is routed to calculate the total runoff.In order to calculate HRUs, the slope was divided into three classes (0%-8%, 8%-30% and >30%) and a threshold level of 10% was established to facilitate model processing and eliminate minor soils, slopes, and land uses for each subbasin.

Input Data for Hydrological Modelling
As for the data used to carry out the hydrological modelling, catchments were defined based on the digital elevation model (DEM), available on the website of the National Center for Geographic Information [25], with an accuracy of 25 m × 25 m.Meteorological data were obtained from the high-resolution (approximately 12 km × 12 km) gridded data set called SPAIN02.Detailed documentation of the development and analysis of the SPAIN02 data set can be found in Herrera et al. [26].In this study, potential evapotranspiration was simulated using the Hargreaves method [27] due to the fact that it only requires minimum and maximum temperatures.Oudin et al. [28] checked that water balance models using parsimonious temperature-based methods perform similarly well compared to more data-demanding methods.The discharge data at the catchment outlet were available on the Hydrographical Study Centre website [29].
In addition to DEM, the Geographic Information Systems (GIS) input data required to build the SWAT model setup included a land cover map and a soil map.Land cover data were derived from reclassified Corine Land Cover 2006 [30] (1:50,000).The soil data for the HWSRB were obtained from the Harmonized World Soil Database (HWSD), assembled by the Food and Agriculture Organization of the United Nations (FAO).This database provides data for 16,000 map units containing two different soil layers (0-30 and 30-100 cm deep) [31].

Calibration and Validation of the SWAT Model
Sensitivity analysis was conducted to identify the most influential parameters for streamflow simulation, which were adjusted during calibration.Automatic calibration with the Sequential Uncertainty Fitting programme algorithm (SUFI-2) [32] was run with the sensible parameters and with other relevants in baseflow, groundwater, and runoff to improve the fit.SUFI-2 is a stochastic calibration that provides some relation between calibration and the uncertainty associated with ignorance about natural systems and all other sources, such as driving variables, conceptual model, parameters, and measured data.Detailed documentation of the SUFI-2 algorithm can be found in Abbaspour et al. [33].
Table 1.Evaluation criteria for model performance [34].Moriasi et al. [35] recommended the Nash-Sutcliffe efficiency (NSE), root mean square error to the standard deviation ratio (RSR), and percent bias (PBIAS) as evaluation criteria for model performance.

Climate Scenarios and Statistical Bias Correction Method
The climate simulations used in this study consisted of 9 combinations of General Circulation Models and Regional climate models (GCM-RCM) from the EURO-CORDEX initiative [22] with a grid spacing of about 12.5 km (0.11 • on a rotated grid).The EURO-CORDEX is an international climate downscaling initiative that aims to provide high-resolution climate scenarios for Europe [36].The Water 2017, 9, 149 5 of 14 simulations have been produced assuming concentration pathways RCP4.5 and RCP8.5, described in van Vuuren et al. [37], and are listed in Table 2.For this study, 30 years of data from historical simulation runs  were used as the baseline period.The future climate is represented with two 30-year periods from the scenario simulation runs; medium term (2041-2070) and long term (2071-2100).It is well known that climate model output data contain systematic errors and cannot be used directly in hydrological simulations [38].That is why a bias correction technique was also applied to the downscaled data to increase the accuracy of the results.In this study, the bias correction technique based on distribution mapping of precipitation and temperature was applied.The idea of distribution mapping is to correct the distribution function of climate model values to agree with the observed distribution function.In 2012, Teutschbein and Seibert [39] evaluated different methods for bias correction of regional climate model simulations for hydrological climate change impact studies, and they obtained very good results applying this technique.In order to extract and bias correct data obtained from the climate models, the CMhyd tool was used [40].

Fuzzy TOPSIS
Fuzzy TOPSIS is based on the distance of each indicator for each regional climate model from the ideal solution.The steps followed in the application of this technique first include the determination of the fuzzy decision matrix, taking into account the number of climate models used and the indicators evaluated.After that, in order to homogenise the evaluation supplied for all the criteria, their values must be linearly normalised.Finally, the proximity coefficients (D + i ,D − i ) for each alternative are calculated in accordance with ideal and anti-ideal values selected for each indicator.This technique is designed to minimise the distance of a data object from the positive ideal solution (D + i ) and maximise the distance from the negative ideal solution (D − i ) [14].The closeness coefficient (C i ) of each alternative is calculated as: To establish the ranking of climate models, it is sufficient to sort them according to the decreasing values of their closeness coefficient.A clear example of a fuzzy approach to ranking climate models can be found in Raju and Kumar (2015) [20].The climatic variables used were precipitation, minimum temperature, and maximum temperature.The correlation coefficient (CC), normalised root mean square deviation (NRMSD), and skill score (SS) [41] were used as performance indicators.Equal weights were considered for each criterion, and ideal and anti-ideal values for all the indicators were chosen as (1, 1, 1) and (0, 0, 0).

Calibration and Validation
A global sensitivity analysis found the following sensible parameters (Table 3): SOL_AWC, LAT_TTIME, SOL_BD, GW_REVAP, ALPHA_BF, RCHRG_DP, SOL_K, FFCB, OV_N, and GWQMN.The presence of several of these parameters showed the great importance of groundwater (ALPHA_BF, RCHRG_DP, GWQMN, and GW_REVAP) and lateral flow (LAT_TTIME) in this area, as is described in Conan et al. [42] and Galván et al. [43], where shallow aquifers have a relevant role.Some soil properties influence the opposition to the groundwater, and this justifies its presence as a sensible parameter, like SOL_AWC, SOL_K, and SOL_BD, usually listed in other studies [44].In addition to these parameters, CN2 was considered due to its correlation with runoff production, affecting baseflow as well.In addition, in a Mediterranean area where evapotranspiration has high relevance, ESCO was added because of its function in driving the extraction of the evaporative demand from lower soil layers [45,46].
As shown in Table 3, calibration with the SUFI-2 algorithm provided the best fitted values for the parameters.The best iteration was a very good performance based on performance criteria [35], with 0.80 NSE, 1.22% PBIAS, and 0.45 RSR.
The calibrated parameters were similar to previous references in areas with similar Mediterranean, warm, or semi-arid climatic or vegetation characteristics.The GW_REVAP value is close to the upper limit of the range due to the presence of forests in the area, as occurs with OV_N [47]; this allows the transfer of water to the root zone and increases evapotranspiration [46], which affects the baseflow and points out the importance of evapotranspiration in the Mediterranean balance.ALPHA_BF has a value that set aquifers as a medium-high velocity response to recharge [46].The FFCB value is similar to others demonstrated in other Mediterranean watersheds [48].
Validation was required after calibration.The model was run for the 1971-1987 period, including five years for a warm-up period and, as in calibration, comparing monthly streamflow observed values with the simulated values.The model performance in this step was still accurate, with statistics described as good or very good, as shown in Table 4.Only PBIAS was worse than in the calibration period.
Figure 2 shows an accurate global model performance comparing the simulated and observed values; although the calibration period was drier than the validation and the streamflow had maximum peaks in the validation period that the SWAT overestimated, this is not an unusual issue when a model is implemented [7,43].

Selection of RCMs Using TOPSIS
As shown in Table 5, data set grids relating to the period 1971-2000 obtained from SPAIN02 [26] were compared with historical runs obtained from each of the 9 RCMs to assess the CC, NRMSD, and SS under a fuzzy approach.Table 6 presents , , , and a ranking pattern for every regional climate model used.The top three positions are occupied by RCA4_CNRM-CM5, RCA4_EC-EARTH, and HIRHAM5_EC-EARTH, with the relative closeness of 0.6216, 0.6213, and 0.4928.On the contrary, the seventh, eighth, and ninth positions are occupied by CCLM4-8-17_MPI-ESM-LR, WRF331F_IPSL-CM5A-MR, and REMO2009_MPI-ESM-LR, with the relative closeness of 0.1953, 0.1952, and 0.0951.These results suggest that RCA4_CNRM-CM5 and RCA4_EC-EARTH are suitable as input data for the SWAT modelling application in the case study.

Selection of RCMs Using TOPSIS
As shown in Table 5, data set grids relating to the period 1971-2000 obtained from SPAIN02 [26] were compared with historical runs obtained from each of the 9 RCMs to assess the CC, NRMSD, and SS under a fuzzy approach.Table 6 presents D + i , D − i , C i , and a ranking pattern for every regional climate model used.The top three positions are occupied by RCA4_CNRM-CM5, RCA4_EC-EARTH, and HIRHAM5_EC-EARTH, with the relative closeness of 0.6216, 0.6213, and 0.4928.On the contrary, the seventh, eighth, and ninth positions are occupied by CCLM4-8-17_MPI-ESM-LR, WRF331F_IPSL-CM5A-MR, and REMO2009_MPI-ESM-LR, with the relative closeness of 0.1953, 0.1952, and 0.0951.These results suggest that RCA4_CNRM-CM5 and RCA4_EC-EARTH are suitable as input data for the SWAT modelling application in the case study.
CCLM4-8-17_CNRM-CM5 0.5268 0.0000 0.0000 0.0000 0.9721 0.7978 0.0000 0.9630 0.7954 CCLM4-8-17_MPI-ESM-LR 0.0713 0.0588 0.9288 0.0766 0.0000 0.0000 0.0764 0.0000 0.0000 HIRHAM5_EC-EARTH 0.9960 0.0000 0.0000 0.9401 0.7408 0.0000 0.9389 0.7338 0.0000 REMO2009_MPI-ESM-LR 0.0000 0.3726 0.0000 0.0000 0.0000 0.1681 0.0000 0.0000 0.1692 RACMO22E_EC-EARTH 0.0000 0.8852 0.3098 0.0700 0.9579 0.0000 0.0699 0.9489 0.0000 RCA4_CNRM-CM5 0.0000 0.9725 0.0000 0.8774 0.8199 0.8003 0.8785 0.8483 0.8059 RCA4_EC-EARTH 0.8780 0.0386 0.7276 0.4320 0.7606 0.9630 0.4314 0.7892 0.9698 RCA4_MPI_ESM_LR 0.5958 0.0344 0.6874 0.0000 0.2168 0.0000 0.0000 0.2850 0.0000 WRF331F_IPSL-CM5A-MR 0.0713 0.0588 0.9288 0.0766 0.0000 0.0000 0.0764 0.0000 0.0000  The mean annual projected precipitation and temperature are displayed in Table 7. Similar patterns are observed, with a general reduction in precipitation and a general increase in temperature.Expected changes in precipitation under RCP4.5 are totally different depending on the model analysed.While RCA4_EC-EARTH projects a slight increase in precipitation, RCA4_CNRM-CM5 projects a reduction that ranges between 13% and 19%.This significant variability was also found in the climate projections published by the Spanish agency of meteorology [49].Under the RCP8.5 scenario, both models show a negative trend in precipitation that ranges between 6% and 17% in the medium term and 32% in the long term.With regards to projected temperature compared to the baseline period, the mean annual temperature suggests a significant and steady increase across the HWSRB in both scenarios.The temperature increase across the HWSRB will range between 0.9 and 1.3 • C in the medium term and 1.3 and 1.8 • C in the long term in the RCP4.5 scenario and between 1.5 and 2.1 • C in the medium term and 2.7 and 3.3 • C in the long term in the RCP8.5 scenario.As for temperature, there were no exceptions; all models showed a higher increase in temperature in the long term compared to the medium term.These results are consistent with other studies in Spanish Mediterranean areas [15,50].
As shown in Figure 3, the temperature increase in winter and autumn is lower than the increase in warmer months, increasing the intra-annual difference of the temperature.The lowest increase in temperature occurs in winter in the medium term in RCP4.5 (0.5 to 1.1 • C), while the largest increases occur during summer in the long term in RCP8.5 (2.9 to 3.8 • C).These results are consistent with other studies in small Mediterranean basins [15,50].With respect to seasonal precipitation change (Figure 4), both models agreed on projecting a decrease in the precipitation for winter, spring, and autumn, while the precipitation will not suffer significant variations in summer.As shown in Figure 3, the temperature increase in winter and autumn is lower than the increase in warmer months, increasing the intra-annual difference of the temperature.The lowest increase in temperature occurs in winter in the medium term in RCP 4.5 (0.5 to 1.1 °C), while the largest increases occur during summer in the long term in RCP 8.5 (2.9 to 3.8 °C).These results are consistent with other studies in small Mediterranean basins [15,50].With respect to seasonal precipitation change (Figure 4), both models agreed on projecting a decrease in the precipitation for winter, spring, and autumn, while the precipitation will not suffer significant variations in summer.

Annual Streamflow Change
Under the RCP 4.5 emissions scenario, the variability in the projected streamflow is very high, ranging from +4% to −35% in the medium term and projecting in the long term an important decrease that is expected to range between 2% and 23% (Figure 5).The possibility of a slight increment in the streamflow, as can be seen for the RCA4_EC-EARTH model, agrees with estimated projections by the Spanish government for the SRB [51], as does the prediction of a higher streamflow reduction in the medium term compared with the long term.These results can also be compared to those obtained by Estrela et al. [3], who also projected a reduction in mean annual runoff for the SRB between 21% and 33%.Overall, the increase in temperature and the projected decrease in precipitation will result in increased evapotranspiration, which will interact to reduce streamflow significantly [52].Comparing precipitation and streamflow results, it can be seen that, due to higher actual evapotranspiration, the decreases obtained in streamflow exceed those in precipitation by 20%.These results are consistent with other studies in Mediterranean climates [53], in which the streamflow is very sensitive to a decrease in precipitation in basins with high evapotranspiration rates.

Annual Streamflow Change
Under the RCP4.5 emissions scenario, the variability in the projected streamflow is very high, ranging from +4% to −35% in the medium term and projecting in the long term an important decrease that is expected to range between 2% and 23% (Figure 5).The possibility of a slight increment in the streamflow, as can be seen for the RCA4_EC-EARTH model, agrees with estimated projections by the Spanish government for the SRB [51], as does the prediction of a higher streamflow reduction in the medium term compared with the long term.These results can also be compared to those obtained by Estrela et al. [3], who also projected a reduction in mean annual runoff for the SRB between 21% and 33%.Overall, the increase in temperature and the projected decrease in precipitation will result in increased evapotranspiration, which will interact to reduce streamflow significantly [52].Comparing precipitation and streamflow results, it can be seen that, due to higher actual evapotranspiration, the decreases obtained in streamflow exceed those in precipitation by 20%.These results are consistent with other studies in Mediterranean climates [53], in which the streamflow is very sensitive to a decrease in precipitation in basins with high evapotranspiration rates.
streamflow reduction in the medium term compared with the long term.These results can also be compared to those obtained by Estrela et al. [3], who also projected a reduction in mean annual runoff for the SRB between 21% and 33%.Overall, the increase in temperature and the projected decrease in precipitation will result in increased evapotranspiration, which will interact to reduce streamflow significantly [52].Comparing precipitation and streamflow results, it can be seen that, due to higher actual evapotranspiration, the decreases obtained in streamflow exceed those in precipitation by 20%.These results are consistent with other studies in Mediterranean climates [53], in which the streamflow is very sensitive to a decrease in precipitation in basins with high evapotranspiration rates.

Seasonal Streamflow Change
Figure 6 shows the seasonal streamflow changes resulting from the estimated scenarios.Overall, a general seasonal streamflow decrease is expected for both scenarios and models.Only in winter and spring does the RCA4_EC-EARTH model estimate an increase of streamflow, which is consistent with the increment of precipitation estimated by this model.In summer, despite a lack of a clear decrease in the precipitation, an important decrease in streamflow is estimated due to the increase in the temperature, which would cause an increase in the evapotranspiration.Figure 6 shows the seasonal streamflow changes resulting from the estimated scenarios.Overall, a general seasonal streamflow decrease is expected for both scenarios and models.Only in winter and spring does the RCA4_EC-EARTH model estimate an increase of streamflow, which is consistent with the increment of precipitation estimated by this model.In summer, despite a lack of a clear decrease in the precipitation, an important decrease in streamflow is estimated due to the increase in the temperature, which would cause an increase in the evapotranspiration.

Spatial Assessment of Expected Changes
As shown in Table 8, a significant agreement for the whole basin area was found due to the reduced extension and the homogeneous hydrological characteristics of the HWSRB.However, sub-basin 1 will suffer slightly higher reductions in its water resources.

Spatial Assessment of Expected Changes
As shown in Table 8, a significant agreement for the whole basin area was found due to the reduced extension and the homogeneous hydrological characteristics of the HWSRB.However, sub-basin 1 will suffer slightly higher reductions in its water resources.

Conclusions
In this study the impacts of projected climate change on water resources in the HWSRB were assessed.The Soil and Water Assessment Tool was used to simulate watershed hydrological processes, and the fuzzy TOPSIS technique was applied in order to select suitable RCM-GCM combinations and reduce uncertainties associated with climate modelling.Simulations with the calibrated model were then conducted for the medium term (2041-2070) and the long term (2071-2100) under two different representative concentration pathways, RCP4.5 and RCP8.5, based on CMIP5.The main findings can be summarised as follow:

•
The SWAT model was able to reproduce the current hydrological conditions of the basin.The statistical results of calibration were NSE = 0.80, RSR = 0.45, and PBIAS = 1.22.The validation results were NSE = 0.77, RSR = 0.48, and PBIAS = −12.68.These results are indicative of the SWAT model's good performance.

•
NRMS, CC, and SS were used to rank nine coupled runs of the GCM-RCM model, applying the fuzzy TOPSIS technique.Higher relative closeness was obtained by RCA4_CNRM-CM5 and RCA4_EC-EARTH.That is why these combinations were suggested for the assessment of the impact of climate change in the HWSRB.

•
Based on the future projections, average annual temperature will increase about 3 • C and precipitation will decrease by 32% by the end of the century.

•
Compared with the baseline period , future water resources in the HWSRB will experience a considerable change as the result of the changing temperature and precipitation.
In the medium term (2041-2070), streamflow presents a high variability under the RCP4.5 and RCP8.5 scenarios, respectively.The largest alterations to streamflow are projected under RCP8.5 for 2071-2100, when they will decline between 33% and 54%.

•
The results obtained from this modelling study may have strong implications in a basin that is already suffering from high water stress.The simulated impacts of climate change should be incorporated into water resource management plans to develop sustainable strategies.Future strategies should be focussed on decreasing demands and increasing the amount of unconventional water resources, but if the magnitude of the climate change renders these strategies insufficient, the need for new water transfer from another basin could arise.

Figure 1 .
Figure 1.Location of the headwaters of the Segura River Basin (HWSRB).
setting up the hydrological model with observed stream flow and climate data; (2) selecting climate projections based on the fuzzy TOPSIS technique; and (3) incorporating climate scenarios into the hydrological model to evaluate the impact of the climate change in the headwaters of the SRB in the medium term (2041-2070) and long term (2071-2100).

Figure 1 .
Figure 1.Location of the headwaters of the Segura River Basin (HWSRB).
In order to evaluate the future impacts of climate change on water resources, the hydrological cycle was simulated in the headwaters of the SRB under different climate change scenarios.This process included three main steps: (1) setting up the hydrological model with observed stream flow and climate data; (2) selecting climate projections based on the fuzzy TOPSIS technique; and (3) incorporating climate scenarios into the hydrological model to evaluate the impact of the climate change in the headwaters of the SRB in the medium term (2041-2070) and long term (2071-2100).

Table 2 .
Overview of the Regional Climate Models (RCMs) considered.

Table 3 .
Range and final parameter values after calibration.

Table 4 .
Calibration and validation period performance.

Table 6 .
Ranking pattern of global climate models.

Table 6 .
Ranking pattern of global climate models.

Table 7 .
Precipitation (mm) and temperature ( • C) means with their variations.

Table 8 .
Spatial streamflow means distribution (m 3 /s) and relative changes of runoff.

Table 8 .
Spatial streamflow means distribution (m 3 /s) and relative changes of runoff.