Journal topic
Hydrol. Earth Syst. Sci., 22, 6473–6491, 2018
https://doi.org/10.5194/hess-22-6473-2018
Hydrol. Earth Syst. Sci., 22, 6473–6491, 2018
https://doi.org/10.5194/hess-22-6473-2018

Research article 13 Dec 2018

Research article | 13 Dec 2018

# Climate change effects on the hydrology of the headwaters of the Tagus River: implications for the management of the Tagus–Segura transfer

Climate change effects on the hydrology of the headwaters of the Tagus River: implications for the management of the Tagus–Segura transfer
Francisco Pellicer-Martínez1 and José Miguel Martínez-Paz2 Francisco Pellicer-Martínez and José Miguel Martínez-Paz
• 1Department of Civil Engineering, Catholic University of Murcia, Spain
• 2Department of Applied Economics, University of Murcia, Spain

Correspondence: Francisco Pellicer-Martínez (fpellicer@ucam.edu)

Abstract

Currently, climate change is a major concern around the world, especially because of the uncertainty associated with its possible consequences for society. Among them, fluvial alterations can be highlighted in basins whose flows depend on groundwater discharges and snowmelt. This is the case of the headwaters of the Tagus River basin, whose water resources, besides being essential for water uses within this basin, are susceptible to being transferred to the Segura River basin (both basins are in the Iberian Peninsula). This work studies the possible effects that the latest climate change scenarios may have on this transfer, one of the most important ones in southern Europe. In the first place, the possible alterations of the water cycle of the donor basin were estimated. To do this, a hydrological model was calibrated. Then, with this model, three climatic scenarios were simulated, one without climate change and two projections under climate change (Representative Concentration Pathways 4.5 (RCP 4.5) and 8.5 (RCP 8.5)). The results of these three hydrological modelling scenarios were used to determine the possible flows that could be transferred from the Tagus River basin to the Segura River basin, by simulating the water resource exploitation system of the Tagus headwaters. The calibrated hydrological model predicts, for the simulated climate change scenarios, important reductions in the snowfalls and snow covers, the recharge of aquifers, and the available water resources. So, the headwaters of the Tagus River basin would lose part of its natural capacity for regulation. These changes in the water cycle for the climate change scenarios used would imply a reduction of around 70 %–79 % in the possible flows that could be transferred to the Segura basin, with respect to a scenario without climate change. The loss of water resources for the Segura River basin would mean, if no alternative measures were taken, an economic loss of EUR 380–425 million per year, due principally to decreased agricultural production.

1 Introduction

Currently, there are practically no doubts in the scientific community that the Earth is suffering climate change (CC) and that this is due to the anthropic action of greenhouse gas emissions (IPCC2014). At the global level, general circulation models predict a warming of the planet of about 2 C for the year 2050, which will cause a reduction in accumulated ice masses and a rise in the sea level (IPCC2014). These changes in the natural environment, which are already causing alterations in the available resources, have a clear socioeconomic repercussion: decreases in fish stocks, increases in energy consumption, changes in the availability of water resources, and land degradation due to erosion. In areas with greater risk and with low capacity for adaptation, the consequences of CC could become critical, the emigration of their population being the only viable solution .

In the context of CC, water as a resource plays a fundamental role, since – as well as being a basic environmental asset – it is key to human survival and well-being. In general terms, an increase in rainfall in humid areas and a decrease in arid and semi-arid ones are predicted. This situation would be accompanied by an increase in the frequency and intensity of extreme events (droughts and floods), so in areas where already there are shortages of water resources the current situation would be exacerbated (IPCC2014). For areas close to polar regions or mountainous areas, the increase in temperature would reduce the precipitation that falls as snow , as well as the volume of ice in the glaciers and the snow cover on the summits . As a consequence, the fluvial regime for this type of river basin would be modified , with an increase in the risk of floods and the loss of the part of their natural regulatory capacity provided by the ice and snow covers . Thus, in areas where the water reserves derived from snow covers are used during the summer, new reservoirs would have to be built in order to replace the loss of the natural regulation capacity (Özdoǧan2011). Similarly, in areas with important aquifers, part of the natural regulatory capacity could also be lost, since changes in precipitation patterns would affect recharge rates (Smerdon2017). Indeed, a higher intensity of rainfall would favour surface runoff to the detriment of infiltration . In short, the CC predicted for many areas of the planet will suppose an increase in temperature together with a greater availability of surface water – which would increase the water evapotranspiration, accelerating the water cycle and reducing the available water resources that could be used .

The environmental and social repercussions of these physical effects are being studied from multiple perspectives . For example, there is work related to water quality , the effects on ecosystems and the services they provide , and the impact on the food security and water security of the population, although the majority of the studies deal with the impacts on the economic activities which are more sensitive to the availability of water resources, such as agriculture , urban supply , or the hydroelectric sector .

Water resource transfers between basins (inter-basin water transfer: IBWT) are instruments of water resource allocation that, despite the controversy they sometimes provoke, can play an important role in mitigating the effects of CC in many areas of the world . The IBWTs can be an alternative source of supply to basins affected by a decrease in their available resources and/or an increase in their water use demands . In turn, the quantity of available water in the donor basins can change significantly, which makes CC a determining factor that must be analysed when assessing the potential or vulnerability of the IBWT . Although there are some works that specifically studied the effects of CC on IBWTs, basically the focus has been on changes in the fluvial regimes in donor basins . Moreover, there are scarce examples in the specialized literature of the analysis of the effects of CC within a framework of integrated water resource management (Onagi2016). In addition, for an adequate comprehensive study of the effects of CC in an IBWT, and to produce operational indicators for water management plans , including the management of water trade , it is also necessary to consider the effects in the receiving basin.

The Tagus–Segura Aqueduct (TSA) is one of the most important IBWT projects in southern Europe. This hydraulic infrastructure, in operation since 1979, transfers flows from the Tagus Headwaters River Basin (THRB) to the Segura River basin (SRB). The destination of the volumes transferred, which are variable depending on the available water resources in the donor basin, is basically irrigation, but also urban and tourism uses . The regional models of CC forecast, for both basins, an increase in temperature together with a significant diminution in rainfall, so that a decrease in the available water resources in both is foreseen (CEDEX2011a). In addition, the THRB is located in a high, mountainous area where snowfalls are frequent and which extends over important karst aquifers. Then, it is expected that both the precipitation that falls as snow and the aquifer recharge would be reduced. Although there is work in which this problem is explicitly described, with proposals to mitigate the decrease in TSA flows due to CC , no specific modelling of climate scenarios has been made, neither from a hydrological perspective, to determine the water balance, nor by simulation of the water resource exploitation system.

The overall objective of this work was to determine the hypothetical effects that CC may have on the operation of the TSA. So, the transferable flows were estimated considering explicitly the operating rule of this IBWT – which, basically, is based on the available water storage in the main two reservoirs of the donor basin concerned. For this, first of all, the effects on the water cycle of the donor basin were evaluated by means of hydrological modelling in which the precipitation that falls as snow was included. Then, the historical climate data together with two CC scenarios of the fifth assessment report (AR5) of the Intergovernmental Panel on Climate Change (IPCC2013) were recreated in order to analyse the possible alterations in the fluvial regime of the THRB. The results of these three hydrological models were the inputs of the subsequent three simulations of the Tagus Headwaters Water Resources Exploitation System (THWRES), which provided a prediction of the flows that could be transferred to the SRB within a framework of integrated water resource management. Additionally, as another novel contribution of this work, the socioeconomic impacts produced by the climate change on these transferred flows were assessed. Finally, note that the complete methodology was developed by open-source tools and by free software for the scientific community, which facilitates the reproducibility of the work.

Figure 1The Tagus Headwaters River Basin (THRB) and Tagus Headwaters Water Resources Exploitation System (THWRES): the main physical characteristics and water uses (including the TSA). Source data: CHT (2018).

2 The water resource exploitation system of the Tagus Headwaters River Basin

The THRB covers an area of 7000 km2 and is located in the middle of the Iberian Peninsula (Fig. 1). It extends over a high, mountainous area with a continental–Mediterranean (CHJ2016) climate with a marked seasonality between the summer (June–September) and winter (December–March) months . The average annual precipitation is 620 mm and the minimum values occur in summer (June–August). While the average annual temperature is 11 C, in the coldest months there are values less than zero (November–April), so snowfalls are frequent at the higher altitudes . Moreover, much of the THRB extends over karst aquifers, meaning that groundwater exerts an important influence on the surface flows that circulate along the main river streams .

Regarding the water uses within the THRB (urban, industrial, and irrigation), their water necessities represent a very low percentage of the available water resources (around 1000×106 m3 yr−1, on average, in the last 70 years): this area has a low population density, is not conducive to agriculture, and its industrial facilities (hydroelectric power stations and an important thermonuclear power station) do not consume much water.

The water resources generated within the THRB are fundamental for the water uses located downstream of the Entrepeñas and Buendía reservoirs (EBR): irrigated agriculture, urban supply (including the city of Madrid), generation of electric power, and maintenance of environmental flows until the city of Aranjuez. Moreover, a large part of these water resources (up to 650×106 m3 yr−1) is susceptible to being transferred to the neighbouring Guadiana River basin and, further away, to the SRB. The former can receive up to 50×106 m3 yr−1 (BOE2014, 2015), of which 20×106 m3 are for the maintenance of the wetland of Tablas de Daimiel and 30×106 m3 are for urban supply to the populations located in the upper Guadiana River basin (CHG2016). The SRB can receive up to 600×106 m3 yr−1 (gross volume), the maximum monthly flow being 68×106 m3. This IBWT is managed by a complex operating rule that gives priority to the water uses in the Tagus River basin and basically depends on the volume stored in the EBR, which have a total storage capacity of 2494×106 m3. The operating rule (Fig. 2) consists of two conditioning factors (BOE2014, 2015). The first restricts the maximum volume that can be transferred in each hydrological year (October–September) to 650×106 m3. The second establishes the transferable flows for each month according to four levels created from two variables: Vacu, the accumulated volume stored in the EBR at the beginning of the month, and Aacu, the accumulated volume of the flows that entered into the EBR in the previous 12 months. The four levels are the following.

• Level 4. When Vacu is lower than 400×106 m3. Transfers are not allowed (Qtrans=0 m3 month−1).

• Level 3. When Vacu is between 400×106 m3 and the values indicated in Fig. 2, which vary between 586×106 m3 and 688×106 m3, depending on the month. Transfers (Qtrans) of 20×106 m3 month−1 are allowed.

• Level 2. When Vacu is between the volumes established in Level 3 and 1500×106 m3, and in addition Aacu is lower than 1000×106 m3. Transfers (Qtrans) of 38×106 m3 month−1 are allowed.

• Level 1. When Vacu is equal to or greater than 1500×106 m3, or Aacu is equal to or greater than 1000×106 m3. Transfers (Qtrans) of 68×106 m3 month−1 are allowed.

Figure 2Operating rule of the TSA. Based on BOE (2014, 2015).

3 Methodology

The methodology applied to determine the maximum monthly volumes that can be transferred from the THRB to the SRB was structured in two stages (Fig. 3).

Figure 3Methodological framework.

The first stage consisted of modelling the hydrology of the THRB until the EBR. For that, a hydrological model was calibrated for the most recent observed flows. Then, with the calibrated model, three scenarios were recreated. In the first one the historical climate series were used (No CC) without climatic correction coefficients. In the others, the data from two characteristic CC scenarios of AR5 (Representative Concentration Pathways 4.5 (RCP 4.5) and 8.5 (RCP 8.5)) were used (IPCC2013). In the first (RCP 4.5), CO2 emissions increase in the future, until in 2050 they stabilize (stabilization scenario), while in the second (RCP 8.5), a continuous and greater increase in emissions of CO2 is assumed (scenario of very high emissions).

The second stage consisted of simulating the THWRES by means of a decision support system (DSS). This simulation incorporated the future water uses contemplated by the water management board. As water uses downstream of the EBR have priority over possible transferable flows (CHT2015), they were included in this simulation.

The methodological framework stages were developed with open-source tools. QGIS was used in the data processing of spatial information, R was employed for data analysis and hydrological modelling , and the DSS SIMGES was used for the simulation of the water resource exploitation system .

Finally, once the series of transferrable flows were calculated, as a complementary goal of this work, an assessment of the socioeconomic consequences that climate change effects have on the main destiny of these flows, the SRB, was made.

## 3.1 Hydrological modelling

### 3.1.1 abcd water balance model with snowmelt module

The hydrological modelling was carried out using the abcd water balance model (Thomas1981). It was applied in a semi-distributed manner , allowing the use of all the gauging stations in the calibration (and validation) in order to maintain the spatial heterogeneity that defines the parameters and variables. This conceptual model was improved by taking into account the hydrological processes of snow and melting. This water balance model and this structure were selected in order to facilitate the understanding of the developed process, allowing the potential reproducibility of the work.

The hydrological modelling, whose scheme is shown in Fig. 4, began with the snowmelt module proposed by . This module recreates the hydrological processes of precipitation as snow (Sn), snow accumulation on the summits (Snp), and snowmelt (Sm). One equation, with a parameter that depends on the temperature, establishes which part of the precipitation (P) occurs as rain (Rf) and which part occurs as snow (Sn). Snow is accumulated in a storage called snowpack (Snp). Then, another equation controlled by one parameter, which also depends on the temperature, establishes the snowmelt (Sm) of the snowpack when the temperature increases. In the modelling, the melting snow ends up forming part of the storage that represents the soil moisture (S). This snowmelt module has been used in previous work, for example, by in basins of China, and by for the Segura headwaters river basin.

The modelling continued with the incorporation of rainfall (Rf) and melting snow (Sm) into the abcd model. This water balance model simplifies the hydrological cycle in two storages, one that simulates the soil moisture balance (S) and another that represents the groundwater storage (G). The model has four parameters – a, b, c, and d (Thomas1981) – that give the model its name. The parameters “a” and “b” manage the soil moisture balance (S), establishing evapotranspiration (ET) and the water susceptible to being lost by surface run-off and/or percolation (QsG). The third parameter “c” identifies percolation (ΔG) towards aquifers (G) and surface run-off (Qs). The fourth parameter “d” determines the discharge from the aquifer (Qg). The output variables of the model are evapotranspiration (ET) and surface run-off (Q).

Figure 4Conceptual scheme of the abcd model.

### 3.1.2 Calibration–validation process

The parameters were calculated by a cascading calibration process , which consists of determining the parameter values from upstream to downstream. In other words, once the model's parameters are established for upstream catchments, their values become input data in the calibration process of the downstream catchments. The split-sample test, proposed by , is carried out by splitting the data series of observed flows into two periods: calibration and validation.

The objective function used in the calibration is the Nash–Sutcliffe efficiency criterion (ENS). This function quantifies the goodness of fit of the model in order to evaluate the performance of the variables selected for study , which are the flows that go out from each catchment (Q). The ENS varies between $\left]-\mathrm{\infty },\mathrm{1}\right]$ and the closer its value is to 1, the better the performance of the model is. The calibration is developed automatically with the Shuffled Complex Evolution Method (SCE-UA) algorithm . Once the semi-distributed model has been calibrated, another three metrics that evaluate the performance of the model, comparing observed with simulated flows, are calculated: the determination coefficient (R2), the percentage of the bias (PBIAS), and the root mean square error (ERMS) .

## 3.2 Simulation of the water resource exploitation system

The operation of the water exploitation systems for a specific scenario is usually evaluated by a DSS which simulates the part of the water cycle that is anthropically modified. These systems represent the main hydraulic network of an area, generally a river basin, together with its natural and artificial storages (rivers, reservoirs, aquifers, canals, among others), establishing the water uses in them . In this hydraulic network, series of flows are introduced, generally in their natural regime. Then, the DSS simulates this kind of system in an integrated manner, fulfilling the priority criteria among the different water uses and the pre-established operating rules . The DSS provides, as results, flow and volume series related to the main fluxes of the hydraulic network (supplies to uses, water consumptions, returns, evaporation in the reservoirs, etc.). Therefore, it is able to estimate the water uses that are not completely met, which are generally those with lower priority and/or those whose spatial location in the hydraulic network does not make it possible to always guarantee their supply . For example, water uses located upstream of a reservoir have less guarantee than those located downstream.

There are different DSSs for water exploitation systems . But, as was advanced in the Methodology section, in this work the model applied is the simulation module of AQUATOOL called SIMGES (), which is one of the models applied most in Spanish river basins, as well as in other countries (Chile, Italy, Morocco, etc.). SIMGES simulates the water exploitation system on a monthly basis in a conservative flow network, seeking a compatible solution that accomplishes the defined constraints . This DSS was selected since it is able to reproduce complex operating rules such as the TSA operating rule.

## 3.3 Socioeconomic impacts at the Segura River basin (SRB)

The precise quantification of the socioeconomic impact of reductions in the volume of water transferred via the TSA would require new integral simulations of the exploitation system of the receiving basin (SRB), which exceed the scope of this work. However, an initial quantification of this impact has been made, based on the work of , where supply failures were assigned an economic value in irrigation in the SRB, and of , who estimated the economic value of the risk associated with droughts in the Region of Murcia. But, given the order of priority of allocation among the water uses, irrigation would suffer the full brunt of any supply deficit. The almost 270 000 ha of irrigated land in the SRB has a net demand of 1363×106 m3 yr−1 (CHS2015), a large part of which is supplied by the TSA. The whole irrigated area in the SRB is divided into seven irrigation zones (IZs); for each of them the water demand curve is estimated from a linear programming model that optimizes the gross value added (GVA) by the optimal cultivation plans according to the water supply. This crop programming includes the irrigation situations of woody crop maintenance, the change from irrigated to rainfed crops and the abandonment of irrigation plots, as well as the impact on employment.

The modelling of the optimal crop plan for each IZ is determined by the following objective function (1) that maximizes the GVA:

$\begin{array}{}\text{(1)}& max\sum _{i}\left({Y}_{i}\cdot {P}_{i}-{C}_{i}\right)\cdot {L}_{i},\end{array}$

where i denotes crop activities under different management options, Yi is the yield of each crop i, Pi the price received by the farmer, Ci the direct costs of production per unit area, and Li the area dedicated to each activity.

The objective function is subject to the following constraints (Eqs. 2–7):

$\begin{array}{}\text{(2)}& & \sum _{i}{L}_{i}\le {L}_{T},\text{(3)}& & \sum _{i}{q}_{i}\cdot {L}_{i}\le {Q}_{T},\text{(4)}& & {L}_{i}^{R}+{L}_{i}^{S}={L}_{i}^{E},\text{(5)}& & {L}_{i}^{R}+{L}_{i}^{M}+{L}_{i}^{P}={L}_{i}^{E},\text{(6)}& & {L}_{i}^{G}\le {L}_{i}^{E},\text{(7)}& & \sum _{i}{\mathit{lf}}_{i}\cdot {L}_{i}\le L{F}_{T}.\end{array}$

LT is the total available surface irrigable in the IZ; qi is the water requirement of each crop per unit area and QT is the availability of water for the entire campaign in the IZ; ${L}_{i}^{R}$ is the surface of irrigated woody crops; ${L}_{i}^{S}$ is the irrigable surface that goes to rainfed; ${L}_{i}^{M}$ is the irrigable surface of woody crops under maintenance irrigation; ${L}_{i}^{P}$ is the surface of the irrigation plots abandoned; ${L}_{i}^{G}$ is the surface of existing irrigated greenhouses; ${L}_{i}^{E}$ is the existing area of each activity in the basin in the reference year; lfi are the labour requirements for each crop, and LFT is the availability of agricultural labour in the IZ.

The first constraint (Eq. 2) prevents each unit of demand (IZ) from cultivating more area than the available net irrigable area. The following constraint (Eq. 3) represents the limitation of water availability for each IZ. The set of constraints (Eq. 4), (Eq. 5), and (Eq. 6) allow simulation of specific management options to certain crop groups. Constraint (Eq. 4) fixes the total area of woody crops such as almond, olive, and wine, distributed between irrigated and rainfed depending on the availability of water. Constraint (Eq. 5) represents citrus and fruit trees, whose total area is equal to the area actually irrigated plus, in situations of scarcity of resources, the surface under maintenance irrigation and/or loss of trees because of not being able to perform the minimum maintenance irrigation. Constraint (Eq. 6) sets the maximum available area for greenhouse crops in the reference year. Finally, constraint (Eq. 7) represents the limitation of the available labour for each IZ. The programme used allows estimation of the gross margin generated under different water availability assumptions, as well as derivation of water demand curves and the marginal value of this resource (Griffin2006).

The necessary data to characterize the technical coefficients of each IZ have been obtained from the sources indicated in and in , updating all the economic figures to EUR of 2017. The programme was solved for each IZ and for the three climatic scenarios (No CC, RCP 4.5, and RCP 8.5), so the availability of water in each of them changes (QT). The differences in the average volume transferred to the SRB in each CC scenario (RCP 4.5 and RCP 8.5) calculated with SIMGES were distributed proportionally for each IZ, taking as a reference the volume transferred in the scenario without CC. Thereby, the socioeconomic impact due to changes in the availability of water was obtained, ceteris paribus the rest of the parameters of the model. Finally, the comparison among the results obtained from GVA and employment for each scenario are presented in the Discussion section.

4 Data source

## 4.1 Hydrological modelling of the Tagus Headwaters River Basin

The digital elevation model (DEM) employed has a 25 m resolution and is available on the website of the National Geographic Institute of Spain (http://www.cnig.es/, last access: 15 December 2017). This DEM was used as an auxiliary variable in the interpolation models of the climatic variables, and to delimit the main streams and catchments using the D8 algorithm . The locations of the 12 gauging stations, which have observed flows in the same period, were used to establish the outlet points of the 12 catchments into which the THRB was divided (Fig. 5). The data series of the observed flows are available in the gauging yearbook of the Official Gauging Station Network of Spain (MITECO2018) and cover the period from September 1985 to December 2009. These observed flows have been previously naturalized to be used in the calibration–validation process (Wurbs2006). For that, the main human alterations located upstream of each gauging station were undone: regulation and evaporation in the reservoirs, as well as the derivations for urban, agricultural, and industrial uses (also, the returns of these uses were considered). As available data for the 12 gauging stations exist for the same period, it is possible to calibrate the model jointly for all the catchments. Since the objective of this work is to obtain a statistically significant calibration and data for its testing, each series of observed flows was divided into two periods (Klemeš1986), the first (from September 1985 to July 1995) being used in the validation and the second (from August 1995 to December 2009) being used in the calibration. Thus, the parameters to be used for the CC projections were determined with the most recent data.

Figure 5Spatial distribution of climate stations with AR5 information and the grid obtained from Spain02v5.

The historical climate series used in the calibration–validation comprise the period from September 1980 to December 2009. Thus, the 5 years before the observed flow series were used to warm up the hydrological model. Two information sources were used to obtain the monthly climatic series. The first was the Spain02v5 dataset , the historical series used for the calibration–validation of the hydrological model and for its simulation in the scenario without climate change (No CC). For the No CC scenario, these series of data were extended from October 1940 to September 2010 (70 consecutive years), and they were used in the simulation without climatic correction coefficients. The advantages of using Spain02v5 are that the daily series of precipitation and temperature are refined (without outliers and/or in homogeneities), and that they include the spatial variability of the climatic variables in a grid with a 12.5 km resolution (Fig. 5). Since these data are daily, they were aggregated on a monthly basis to apply them in the model. The second source of data was the State Meteorological Agency of Spain (AEMET) and was used for the two CC scenarios of AR5 (IPCC2013). This source (AEMET) provides the regionalized projections of 27 models (13 for RCP 4.5 and 14 for RCP 8.5), using the statistical method of analogues . The daily temperature and precipitation series of the 6 thermometric and 48 precipitation stations closest to the THRB were used (Fig. 5). These daily data were also aggregated to monthly series. Next, the reference historical data series of each model were compared with Spain02v5 data, using as a control period the 1971–2005 interval. The 10 models that best fitted for both temperature and precipitation were assembled using the Simple Average Forecast Combination (SA) and Bias-Corrected Eigenvector Forecast Combination (EIG2) , available in the GeomComb R-CRAN package . Then, both ensembles were also compared with Spain02v5 data using the same control period (1971–2005). Finally, the series obtained by EIG2 were used, with a lower prediction error compared to the rest of the series. For example, the ENS values obtained with EIG2 for the temperatures were 0.87, while in the separate models they never exceeded 0.77. For the precipitations, the ENS value was 0.30, while in the models it was always lower than 0. In addition, another advantage of the EIG2 method is that it allows one to correct the bias produced by predictive models. The period used in the simulation of CC scenarios (RCPs 4.5 and 8.5) is from October 2020 to September 2090 (also 70 consecutive years).

Based on the average monthly temperature data of the three climatic scenarios, the potential evapotranspiration series were estimated using the Thornthwaite method . As this method tends to underestimate the potential evapotranspiration, the series generated were corrected from a linear regression between the estimated series and those used by the SIMPA hydrological model (BOE2007). These series used by SIMPA have already regionalized for the Iberian Peninsula based on the Penman–Monteith method , correcting the underestimation of the Thornthwaite method.

Once the monthly series of precipitation, temperature, and potential evapotranspiration had been calculated, they were spatially interpolated on the cells using the thin-plate splines method (Wahba1990), which is based on local interpolation from polynomials. This method turns out to be relatively robust against non-compliance with the statistical assumptions necessary in methods such as kriging, being used with good results for the interpolation of climatic variables such as rainfall and temperatures . Finally, each catchment was assigned the average value of the cells over which it extends.

Figure 6Scheme of the water resource exploitation system of the Tagus Headwaters, as far as Aranjuez.

## 4.2 Water resource exploitation system of the Tagus Headwaters Basin (THWRES)

The THWRES covers the basin upstream of the EBR and the water uses located in the basin downstream of these reservoirs, up to the city of Aranjuez (Figs. 1 and 6). In its design, the possible future water uses for the years 2016, 2021, and 2033 contemplated by the water management board (CHT2015) were analysed. And, as there are no relevant differences between them, the uses predicted for the year 2033 were taken as the reference for the three climatic scenarios (No CC, RCP 4.5, and RCP 8.5) (Table 1). In addition, as the available water resources in the catchments downstream of the EBR are so low with respect to those generated in the EBR drainage basin, they were neglected in the water resource exploitation system simulation.

Table 1Main water uses in the THWRES (unit: 106 m3 yr−1).

Regarding water uses, the majority of the urban and irrigation uses are concentrated downstream of the EBR, requiring more than 85 % of the total volume demanded. For urban use, a water return of 80 % was considered with respect to the volume supplied. In 2033 irrigation will have a low water return into the system (<1 %) due to improvements in irrigation facilities and in the application systems which will be used, so this low return was implemented in the DSS. The main industrial use of water occurs at the Trillo Nuclear Plant. This industrial use requires a constant flow over time that returns 46 % of the water supplied. Finally, an environmental flow of around 11 m3 s−1 must circulate in the Tagus River from the EBR to the city of Aranjuez. The order of priority in the water resource allocation among the different water uses is urban supplies, Trillo Nuclear Plant, and irrigation. Once the water uses in the THWRES have been supplied, the possibility of authorizing transfers is evaluated (BOE2014, 2015). Among the four uses of the transfer there is no stipulated clear criterion in the operating rule that governs it, since sometimes it is done discretionally depending on the needs or level of urgency of the uses. In fact, the transfer to Tablas de Daimiel has occurred only once, to avoid serious damage to this wetland. Therefore, in order to accomplish the general criteria of Spanish legislation, the order of priority followed is urban supply to Segura (TSA), urban supply to the populations located in the upper Guadiana River basin, and irrigation supply to Segura (TSA) and Tablas de Daimiel.

5 Results

## 5.1 Climate change effects on the hydrology of the Tagus Headwaters River Basin

The values of the criterion coefficients calculated in the hydrological modelling show that the model employed reproduced properly the surface flows in the THRB in the calibration period: high values of ENS and R2, together with low relative errors (ERMS) and volume errors (PBIAS). However, in the validation period, there are some low values for the goodness-of-fit coefficients calculated, indicating that the results of these catchments have greater uncertainty. These results can be explained by the fact that the validation period was used just after the warming-up process. Thereby, it could cause, in some catchments, the warming to extend to a part of the validation period. This would entail the calibration process using a part of the validation to adjust the initial parameters, obtaining worse adjustments in the validation of the model. However, it is important to highlight that the parameters used in the simulations are adjusted with the more recent data, providing a good performance of the surface flows in the THRB. In addition, the best performance in calibration corresponded to the outlets of the catchments of Entrepeñas and Buendía, which were the flows used as the input in the subsequent simulation of the water resource exploitation system. In fact, both catchments had NSE values around 0.80 and low PBIAS (Table 2).

Table 2Results of the calibration–validation process for the 12 gauging stations in the hydrological modelling.

The simulation of the historical climate series of 1940–2010 with the calibrated model provided an average annual resource of 954.6×106 m3 yr−1 (Q). This series was temporally moved to the 2020–2090 time period in order to reproduce a future climate scenario without climate change (No CC). The simulations for climate change scenarios RCP 4.5 and RCP 8.5 with the same calibrated model indicated that the THRB could suffer a considerable loss of its natural water resources. The RCP 4.5 scenario forecasted a value of 575.6×106 m3 yr−1, representing a decrease of 39.7 %, while the RCP 8.5 scenario predicted a 46.6 % decline in resources to 508.9×106 m3 yr−1, on average (Fig. 7). This is due to a combination of a reduction in precipitation (15 % and 20 % for each scenario, respectively) and an increase in potential evapotranspiration resulting from an increase in temperature of 2.2 and 3.4 C, respectively, for CC scenarios RCP 4.5 and RCP 8.5.

Figure 7Surface flows that enter into the Entrepeñas and Buendía reservoirs (EBR).

Regarding the snow-melting modelling, the scenario without CC (historical series) provided an average precipitation as snow (Sn) of 185.4×106 m3 month−1 (considering the whole year), reaching a maximum value of 406.9×106 m3 month−1. The hydrological modelling considered the snow for the months between October and May to be relevant, especially that of December, January, and February. The snow accumulated on the summits (Snp) represented an average reserve for each winter of about 112.7×106 m3, with a maximum value of 481.5×106 m3 (Fig. 8a). For the CC scenarios, the snowfall would significantly decrease; in fact, the models did not detect relevant snow-melting processes for the months of October and May. The snowfall would drop by 68 % for RCP 4.5, the average value being around 59.3×106 m3 month−1, while for RCP 8.5 the snowfall would drop by 90 %, giving an average value of about 19.2×106 m3 month−1. These values would decrease the snow covers to a similar extent. In the RCP 4.5 scenario the average snow cover of each winter would be 29.5×106 m3 (Fig. 8b), which represents a reduction of 74 %, whereas the RCP 8.5 snow cover would suffer a reduction of 80 %, reaching an average value of 22.1×106 m3 (Fig. 8c).

Figure 8Monthly distribution of the snow cover accumulated on the summits (Snp) in the Tagus Headwaters River Basin (unit: 106 m3 month−1): (a) No CC. (b) RCP 4.5. (c) RCP 8.5.

The aquifers' recharge estimated by the hydrological modelling (ΔG) had an average value of about 771.0×106 m3 yr−1 for the historical climate series, with high monthly variability. This value indicates that almost 80 % of the surface flows (Q) have previously passed through aquifers, which underlines the relevance of the groundwater in this river basin. The maximum recharge occurs during the winter and spring months (January–May), the minimum values being found in the summer months (Fig. 9a). For the CC scenarios, the recharge is significantly reduced and concentrated in the months of February and April. In the RCP 4.5 scenario the aquifers' recharge is reduced by 54 % (to 357.7×106 m3 yr−1) and represents 62 % of the surface flow (Q) (Fig. 9b). For RCP 8.5 the aquifers' recharge is reduced by 62 % (to 290.0×106 m3 yr−1) and would represent only 57 % of the surface flow (Q) (Fig. 9c). The snowfall reduction is one of the reasons for this decrease in aquifer recharge, since the melting snow in the model becomes soil moisture (S), as usually happens in nature. These changes in the aquifers' recharge alter the pattern of the groundwater discharge (Qg), which is part of the surface flows (Q). In fact, the relative differences in the possible groundwater discharges between the summer and spring months would be greater (Fig. 10a). Overall, the CC scenarios predict a significant reduction in the water resources with respect to the historical climate series, together with a loss of the natural capacity of regulation of the river basin itself, manifested as an increase in the relative values of the maximum flows (Fig. 10b).

Figure 9Monthly distribution of the recharge (ΔG) in the Tagus Headwaters River Basin (unit: 106 m3 month−1): (a) No CC. (b) RCP 4.5. (c) RCP 8.5.

Figure 10Alteration of the intra-annual pattern of the groundwater discharges from aquifers (a) and surface flows (b), in relative terms (%).

## 5.2 Climate change effects on the Tagus Headwaters Water Resources Exploitation System

The simulation of the THWRES was carried out with the flows obtained in the previous hydrological modelling. The results obtained using the historical climate series (No CC) and the stabilization scenario (RCP 4.5) indicate that the supply to the water uses in the Tagus River basin is guaranteed, in both cases, for the whole of the time period simulated. However, for the scenario of very high emissions (RCP 8.5), the simulation indicates that there would be supply deficits in the THWRES in some years towards the end of the time period considered.

Once the water uses in the Tagus River basin have been supplied, the TSA operating rule comes into effect, providing the monthly volumes that would be possible to transfer in each climatic scenario evaluated. The results are presented in Fig. 11, which shows the volume stored in each month in the EBR (Vacu), the transferred volume (Qtrans) that leaves the Tagus River basin, and the flow that would go through the TSA to the SRB.

Figure 11The storage volume in the EBR (Vacu), the total transferred flow (Qtrans), and the flow to the SRB (TSA) (unit: 106 m3 month−1): (a) No CC. (b) RCP 4.5. (c) RCP 8.5.

Table 3 summarizes the main statistics of the annual transfer series, specifying which part goes to the SRB. If the future climate were similar to that of past years (No CC), the average transferable volume would be about 450×106 m3 yr−1 (Fig. 11a). This volume would break down, following the priority criterion established, into about 40×106 m3 yr−1 for the Guadiana River basin (29×106 m3 yr−1 for urban supply and about 10×106 m3 yr−1 for Tablas de Daimiel) and 411×106 m3 yr−1 for the TSA. Thus, given the average losses by infiltration and evaporation from the transfer channel, which are around 10 % (CHS2015), the net volume that would reach the SRB would be 370×106 m3 yr−1, a value close to those of the actual series of transfers that have occurred since it came into operation . So, if there were No CC, despite the fact that the volume reaching the SRB via the TSA would be, on average, much lower than the planned 540×106 m3 yr−1 (600×106 m3 yr−1 minus 10 % losses), it would be possible to transfer an average of 411×106 m3 yr−1 to the SRB.

Table 3Water volumes transferred (Qtrans and TSA) in each climate scenario (unit: 106 m3 yr−1).

In the RCP 4.5 scenario, the transferable volumes drop to 143×106 m3 year−1 on average, 17×106 m3 yr−1 being for the urban supply in the Guadiana River basin, 2.5×106 m3 yr−1 for Tablas de Daimiel, and 123×106 m3 yr−1 for the SRB by means of the TSA infrastructure. So, considering losses of 10 %, this would mean net water resources of 111×106 m3 yr−1 being transferred to the SRB, barely 20 % of the maximum transferable volume. As worrisome as this important decrease in the average value is the existence of consecutive periods of 3 and 4 years in which no transfer would occur (Fig. 11b).

This situation would be aggravated for the climatic scenario RCP 8.5 since the transferable volume would be reduced to about 100×106 m3 yr−1, on average, distributed as 12×106 m3 yr−1 for the urban supply in the Guadiana River basin, 1.5×106 m3 yr−1 for the maintenance of the Tablas de Daimiel wetland, and 86×106 m3 yr−1 for the TSA. Thus, the SRB would receive approximately throughout the year the volume planned for just 1 month in the operating rule. This scenario worsens the duration and frequency of the no-transfer periods, which intensify over time, reaching a situation of total cessation of the TSA from the year 2067 due to a lack of accumulated volumes (Vacu) in the EBR (Fig. 11c).

6 Discussion

The hydrological modelling carried out with the AR5 CC scenarios predicts a sharp decrease in the water resources of the THRB. These results are in line with the previous simulations made on the same scale , as well as on a larger scale . The hydrological modelling indicated that the increase in temperature would generate a decline in snowfall, which would lead to a reduction in the snow cover period of 2 months. The relevant snowfalls would have a delay of 1 month and the snowmelt would start a month earlier. The decreases provided by the simulations are similar to those published by for the same mountainous area, in which a decrease of up to 70 % was predicted for the interval 2070–2100 in the RCP 4.5 scenario, and greater than 90 % for the RCP 8.5 scenario (same interval). In addition, also predicted a decrease in the 2-month snow cover period. Regarding the aquifer recharge, the hydrological modelling indicated a decline greater than 50 %. Although these values are slightly higher than those published by for the whole Tajo River basin, they are in line with the provisions in the THRB, since the majority of the projections used for these two scenarios (RCP 4.5 and RCP 8.5) indicated diminishments of around 50 % in aquifer recharge within this area. These alterations of the hydrological processes will probably generate a change in the fluvial regime. In relative terms, the months of autumn (October, November, December) would suffer the major decrease, this change being consistent with the results of . At the annual level, these alterations would result in greater variability of the surface flows of the basin due to increases in the relative difference between the maximum and minimum flows of spring and summer, respectively. In short, these changes are going to diminish the natural capacity of the THRB to regulate its own water resources. However, the results of the THWRES simulations indicate that the EBR will not reach their maximum volume in any situation, and no new infrastructure will be necessary to overcome this loss of regulation.

The effects of these physical changes on the THWRES will translate into average decreases around 70 %–79 % in the volume expected to be available for the TSA in both scenarios evaluated (RCP 4.5 and RCP 8.5). In addition, there would be long periods without transfers for the RCP 4.5 scenario, and this IBWT could even stop operating for the last third of the simulation period under the RCP 8.5 scenario, when there would be deficits regarding the water uses in the Tagus River basin itself.

At this point it is necessary to consider two reflections. The first is the uncertainty of regionalized projections, whose origin is associated with the uncertainties inherited during the different stages of their generation. To this uncertainty must be added the influence of uncontrolled local factors in the regionalization, which have a great relevance in the case of precipitation. In this sense, in Spain, while the regionalized series for temperature are in line with the historical reference series, those of precipitation have strong uncertainties associated with the models . Even so, the use of these data sources is the best way to understand the repercussions of CC for water resources and hence to adapt strategies to CC. The second is that the simulations carried out for the water resource exploitation system represent scenarios in which the uses of the water and the operating rule do not change with time. The flows that enter in the model are the only input data that vary. So, they are not tools to predict future behaviour, since water needs, infrastructure, and water policy are not static. However, they do serve to evaluate the effects that CC could have on the TSA if measures are not taken regarding the management of the demand for water or the incorporation of alternative sources of water resources into this system.

The socioeconomic impact of the decreases in the TSA supply were calculated by the methodology presented in Sect. 3.3. First, the global decrease in TSA flows to the SRB was distributed among the seven IZs, according to the results obtained by , which were conditioned by the water demand and by the topology of the exploitation system of the SRB. Once these deficits were determined, they were valued using the linear programming presented, calculating GVA and agricultural employment for each IZ according to water allocations available in each scenario. The results are presented in Table 4, which shows the differences in relation to a scenario with no climate change (No CC).

Table 4Differences of the RCP 4.5 and RCP 8.5 scenarios with respect to No CC (these values represent decreases with respect to the No CC scenario).

By calculating the ratio between the GVA and TSA the marginal value of water used in irrigation was obtained: 1.29 EUR m−3 for the RCP 4.5 scenario and 1.31 EUR m−3 for the RCP 8.5 scenario (average values for the whole SRB). The marginal value establishes the upper limit of the marginal productivity of water, and therefore it represents the maximum payment capacity of the sector for water in each scenario. The figures obtained are in line with those presented in other studies for the area .

Thus, the decrease in the TSA supply for the RCP 4.5 scenario would mean a direct loss of 380 million EUR yr−1 (at the prices of 2017), while for RCP 8.5 it would amount to 425 million EUR yr−1. Given that the GVA of irrigation in the SRB is around 1870 million EUR yr−1 , these figures represent, respectively, direct losses of 20 % and 23 % in terms of GVA. In addition, the decrease in the volumes transferred in the CC scenarios means that the irrigated area is occupied by more labour-extensive crops, causing the loss of around 6690 and 7369 direct full-time agrarian jobs per year in the sector (Table 4). To all these direct effects we should add the drag effects of the primary sector on other economic sectors that depend directly on the use of irrigation in the area (food industry, marketing, transport, inputs, etc.) and make up the well-known agro-industrial cluster of the Region of Murcia . In this sense, some authors propose a multiplier of no less than 3 to calculate the total economic impact of agricultural production on other related activities. To this economic impact it would be possible to add the environmental one , as it would be, for example, the reduction of the flows circulating in the Segura River, since it is part of the distribution network of the TSA volumes, which would also have effects on the ecosystems associated with this river .

7 Conclusions

In this work, the possible effects of CC on the Tagus Headwaters Water Resources Exploitation System have been evaluated. In particular, the work has focused on the hydrological changes that would arise in the Tagus Headwaters River Basin (THRB) if the predicted climatic scenarios arose, and how these would affect the uses of the basin itself and the IBWT (the TSA) that transfers flows to the Segura River basin. The main effects on the water cycle would be the reduction of snowfalls and their accumulation on the summits, and a decrease in aquifer recharge. These changes would generate a significant loss of the natural regulation capacity of the THRB itself, which could not be corrected with new infrastructure since there would also be a significant reduction in the water resources.

The hypothetical flows that can be transferred by the TSA would suffer a significant decrease if the simulated CC scenarios came into being. For RCP 4.5 there would be an average reduction of 70 %, while for RCP 8.5 it would be 79 %, both figures being with respect to the series without CC. These reductions would result in direct economic losses in the irrigation sector of around 20 %. Beyond these average figures, the increase in the zero-transfer periods would have a high impact, since it would further increase the uncertainty associated with this source of supply, preventing adequate planning of water uses in the receiving basin.

Finally, in addition to the undoubted interest that is presented by the case study analysed, this work shows a complete sequential methodological framework which should serve as a guide for the comprehensive evaluation of IBWT within a framework of integrated water resource management. Moreover, the entire methodology has been developed with open-source tools, facilitating its reproducibility in other areas.

Data availability
Data availability.

The climatic historical data were taken from the Spain02v5 dataset available at http://www.meteo.unican.es/datasets/spain02 (Spain02v5, 2018).

Author contributions
Author contributions.

Both authors contributed equally to this work.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

This article is part of the special issue “Assessing impacts and adaptation to global change in water resource systems depending on natural storage from groundwater and/or snowpacks”. It is not associated with a conference.

Acknowledgements
Acknowledgements.

This research work has been supported by project 19342/PI/14 funded by “Fundación Séneca-Agencia de Ciencia y Tecnología de la Región de Murcia” in the framework of PCTIRM 2011–2014. Moreover, the authors thank AEMET and UC for the data provided for this work. We appreciate the valuable comments and suggestions provided by the editor and two anonymous referees.

Edited by: David Pulido-Velazquez
Reviewed by: two anonymous referees

References

Albaladejo-García, J., Martínez-Paz, J., and Colino, J.: Financial evaluation of the feasibility of using desalinated water in the greenhouse agriculture of Campo de Níjar (Almería, Spain), ITEA-Inf. Tec. Econ. Ag., 114, 6, https://doi.org/10.12706/itea.2018.024, 2018 (in Spanish). a

Allen, R. G., Pereira, L., Raes, D., and Smith, M.: Crop evapotranspiration: Guidelines for computing crop requirements, FAO Irrigation and drainage paper 56, Rome, 15 pp., available at: https://appgeodb.nancy.inra.fr/biljou/pdf/Allen_FAO1998.pdf (last access: 18 January 2018), 1998. a

Amblar-Francés, P., Casado-Calle, M., Saavedra, A., Ramos, P., and Rodríguez, E.: Guía de Escenarios Regionalizados de Cambio Climático sobre España a Partir de los Resultados del IPCC-AR5, Agencia Estatal de Meteorología, Gobierno de España, available at: http://hdl.handle.net/20.500.11765/7956, 17 December 2017. a, b

Bajracharya, A., Bajracharya, S., Shrestha, A., and Maharjan, S.: Climate change impact assessment on the hydrological regime of the Kaligandaki Basin, Nepal, Sci. Total Environ., 625, 837–848, https://doi.org/10.1016/j.scitotenv.2017.12.332, 2018. a

BOE: Acuerdo para encomienda de gestión por el Ministerio de Agricultura, Alimentación y Medio Ambiente (Dirección General del Agua) al CEDEX, del Ministerio de Fomento, para la realización de asistencia técnica, investigación y desarrollo tecnológico en materias competencia de la Dirección General del Agua, Boletín Oficial del Estado (BOE), 287, 49436–49458, available at: https://www.boe.es/diario_boe/txt.php?id=BOE-A-2007-20623 (last access: 15 August 2018), 2007 (in Spanish). a

BOE: Real Decreto 773/2014, de 12 de septiembre, por el que se aprueban diversas normas reguladoras del trasvase por el acueducto Tajo-Segura, Ministerio de Agricultura, Alimentación y Medio Ambiente, Boletín Oficial del Estado (BOE), 223, 71634–71639, available at: https://www.boe.es/buscar/doc.php?id=BOE-A-2014-9336 (15 August 2018), 2014 (in Spanish). a, b, c, d

BOE: Ley 21/2015, de 20 de julio, por la que se modifica la Ley 43/2003, de 21 de noviembre, de Montes, Jefatura del Estado, Boletín Oficial del Estado (BOE), 173, 60234–600272, available at: https://www.boe.es/diario_boe/txt.php?id=BOE-A-2015-8146 (15 August 2018), 2015 (in Spanish). a, b, c, d

CEDEX: Estudio de los Impactos del Cambio Climático en los Recursos Hídricos y las Masas de Agua, Informe Final, Centro de Estudios y Experimentación de Obras Públicas, Ministerio de Agricultura, Alimentación y Medio Ambiente, available at: http://www.mapama.gob.es/es/agua/temas/planificacion-hidrologica/planificacion-hidrologica/EGest_CC_RH.aspx (15 August 2018), 2011a. a

CEDEX: Estudio de los Impactos del Cambio Climático en los Recursos Hídricos y las Masas de Agua. Ficha 1: Evaluación del Impacto del Cambio Climático en los Recursos Hídricos en Régimen Natural, Centro de Estudios y Experimentación de Obras Públicas, Ministerio de Agricultura, Alimentación y Medio Ambiente, available at: http://www.cedex.es (15 August 2018), 2011b (in Spanish). a

CEDEX: Evaluación del Impacto del Cambio Climático en los Recursos Hídricos y Sequías en España, Informe Final, Centro de Estudios Hidrográficos, Centro de Estudios y Experimentación de Obras Públicas, Ministerio de Agricultura, Alimentación y Medio Ambiente, available at: http://www.cedex.es/CEDEX/LANG_CASTELLANO/ORGANISMO/CENTYLAB/CEH/Documentos_Descargas/EvaluacionimpactoCCsequiasEspana2017.htm (15 August 2018), 2017 (in Spanish). a, b, c

Chavez-Jimenez, A., Granados, A., Garrote, L., and Martín-Carrasco, F.: Adapting water allocation to irrigation demands to constraints in water availability imposed by Climate Change, Water Resour. Manag., 29, 1413–1430, https://doi.org/10.1007/s11269-014-0882-x, 2015. a

CHG: Plan Hidrológico de la Parte Española de la Demarcación Hidrográfica del Guadiana, Memoria (Parte I), Confederación Hidrográfica del Guadiana, Ministerio de Agricultura, Alimentación y Medio Ambiente, available at: http://planhidrologico2015.chguadiana.es/?corp=planhidrologico2015&url=61 (15 August 2018), 2016 (in Spanish). a

CHJ: Plan Hidrológico de la Demarcación Hidrográfica del Júcar. Memoria, Ciclo de Planificación Hidrológica 2015–2021, Confederación Hidrográfica del Júcar, Ministerio de Agricultura, Alimentación y Medio Ambiente, available at: https://www.chj.es/es-es/medioambiente/planificacionhidrologica/Paginas/PHC-2015-2021-Plan-Hidrologico-cuenca.aspx (15 August 2018), 2016. a

CHS: Proyecto de Plan Hidrológico de la Cuenca del Segura, Anejo 02: Recursos Hídricos, Confederación Hidrográfica del Segura, Ministerio de Agricultura, Alimentación y Medio Ambiente, available at: http://www.chsegura.es/chs/planificacionydma/planificacion/, 2015 (in Spanish). a, b

CHT: Plan Hidrológico de la Parte Española de la Demarcación Hidrográfica del Tajo, Memoria, Confederación Hidrográfica del Tajo, Ministerio de Agricultura, Alimentación y Medio Ambiente, available at: http://www.chtajo.es/LaCuenca/Planes/PlanHidrologico/Planif_2015-2021/Paginas/Plan_2015-2021.aspx (last access: 15 August 2018), 2015. a, b

CHT: Centro de descargas de capas en formato shape de la Cuenca Hidrográfica del Tajo, available at: http://www.chtajo.es/LaCuenca/Paginas/DescargaDCapas.aspx (last access: 17 December 2017), 2018 (in Spanish). a

Colino, J., Martínez-Carrasco, F., and Martínez-Paz, J. M.: El impacto de la PAC renovada sobre el sector agrario de la Región de Murcia, de la Región de Murcia, available at: https://www.cesmurcia.es/cesmurcia/DocumentServlet?docid=/publicaciones/ficheros/564.pdf (last access: 15 September 2018), 2014. a

Díaz, P., Morley, K. M., and Yeh, D. H.: Resilient urban water supply: preparing for the slow-moving consequences of climate change, Water Practice and Technology, 12, 123–138, https://doi.org/10.2166/wpt.2017.016, 2017. a

Duan, Q., Sorooshian, S., and Gupta, V. K.: Optimal use of the SCE-UA global optimization method for calibrating watershed models, J. Hydrol., 158, 265–284, https://doi.org/10.1016/0022-1694(94)90057-4, 1994. a

Flörke, M., Schneider, C., and McDonald, R. I.: Water competition between cities and agriculture driven by climate change and urban growth, Nature Sustainability, 1, 51–58, https://doi.org/10.1038/s41893-017-0006-8, 2018. a

Giupponi, C. and Gain, A. K.: Integrated water resources management (IWRM) for climate change adaptation, Reg. Environ. Change, 17, 1865–1867, https://doi.org/10.1007/s10113-017-1173-x, 2017. a

Gomariz-Castillo, F., Alonso-Sarría, F., and Cabezas-Calvo-Rubio, F.: Calibration and spatial modelling of daily ET0 in semiarid areas using Hargreaves equation, Earth Sci. Inform., 11, 325–340, https://doi.org/10.1007/s12145-017-0327-1, 2018. a

Griffin, R. C.: Water Resource Economics: The Analysis of Scarcity, Policies and Projects, The MIT Press Cambridge, available at: https://mitpress.mit.edu/books/water-resource-economics (last access: 15 September 2018), 2006. a

Grindlay, A., Zamorano, M., Rodríguez, M., Molero, E., and Urrea, M.: Implementation of the European Water Framework Directive: Integration of hydrological and regional planning at the Segura River Basin, southeast Spain, Land Use Policy, 28, 242–256, https://doi.org/10.1016/j.landusepol.2010.06.005, 2011. a

Guerreiro, S. B., Birkinshaw, S., Kilsby, C., Fowler, H. J., and Lewis, E.: Dry getting drier – The future of transnational river basins in Iberia, J. Hydrol., 12, 238–252, https://doi.org/10.1016/j.ejrh.2017.05.009, 2017. a

Gupta, H. V. and Kling, H.: On typical range, sensitivity, and normalization of Mean Squared Error and Nash-Sutcliffe Efficiency type metrics, Water Resour. Res., 47, w10601, https://doi.org/10.1029/2011WR010962, 2011. a

Herrera, S., Fernández, J., and Gutiérrez, J. M.: Update of the Spain02 gridded observational dataset for EURO-CORDEX evaluation: assessing the effect of the interpolation methodology, Int. J. Climatol., 36, 900–908, https://doi.org/10.1002/joc.4391, 2016. a

Hsiao, C. and Wan, S. K.: Is there an optimal forecast combination?, J. Econometrics, 178, 294–309, https://doi.org/10.1016/j.jeconom.2013.11.003, 2014. a

Hutchinson, M. F.: Interpolating mean rainfall using thin plate smoothing splines, Int. J. Geogr. Inf. Syst., 9, 385–403, https://doi.org/10.1080/02693799508902045, 1995. a

IPCC: Climate Change 2013 – The Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, UK, https://doi.org/10.1017/CBO9781107415324, 2013. a, b, c

IPCC: Climate Change 2014 – Impacts, Adaptation and Vulnerability: Part B: Regional Aspects: Working Group II Contribution to the IPCC Fifth Assessment Report, vol. 2, Cambridge University Press, Cambridge, UK, https://doi.org/10.1017/CBO9781107415386, 2014. a, b, c

Jha, C. K., Gupta, V., Chattopadhyay, U., and Sreeraman, B. A.: Migration as adaptation strategy to cope with climate change: A study of farmers' migration in rural India, Int. J Clim. Chang. Str., 10, 121–141, https://doi.org/10.1108/IJCCSM-03-2017-0059, 2018. a

Kahil, M. T., Connor, J. D., and Albiac, J.: Efficient water management policies for irrigation adaptation to climate change in Southern Europe, Ecol. Econ., 120, 226–233, https://doi.org/10.1016/j.ecolecon.2015.11.004, 2015. a

Klemeš, V.: Operational testing of hydrological simulation models, Hydrolog. Sci. J., 31, 13–24, https://doi.org/10.1080/02626668609491024, 1986. a, b

Li, L., Zhang, L., Xia, J., Gippel, C. J., Wang, R., and Zeng, S.: Implications of modelled climate and land cover changes on runoff in the middle route of the south to north water transfer project in China, Water Resour. Manag., 29, 2563–2579, https://doi.org/10.1007/s11269-015-0957-3, 2015. a

Li, Z. L., Xu, Z. X., and Li, Z. J.: Performance of WASMOD and SWAT on hydrological simulation in Yingluoxia watershed in northwest of China, Hydrol. Process., 25, 2001–2008, https://doi.org/10.1002/hyp.7944, 2011. a

Li, Z. L., Shao, Q. X., Xu, Z. X., and Xu, C. Y.: Uncertainty issues of a conceptual water balance model for a semi-arid watershed in north-west of China, Hydrol. Process., 27, 304–312, https://doi.org/10.1002/hyp.9258, 2013. a

Lobanova, A., Koch, H., Liersch, S., Hattermann, F. F., and Krysanova, V.: Impacts of changing climate on the hydrology and hydropower production of the Tagus River basin, Hydrol. Process., 30, 5039–5052, https://doi.org/10.1002/hyp.10966, 2016. a

Lobanova, A., Liersch, S., Nunes, J. P., Didovets, I., Stagl, J., Huang, S., Koch, H., del Rocío Rivas López, M., Maule, C. F., Hattermann, F., and Krysanova, V.: Hydrological impacts of moderate and high-end climate change across European river basins, J. Hydrol., 18, 15–30, https://doi.org/10.1016/j.ejrh.2018.05.003, 2018. a, b, c

Lorenzo-Lacruz, J., Vicente-Serrano, S., López-Moreno, J., Beguería, S., García-Ruiz, J., and Cuadrat, J.: The impact of droughts and water management on various hydrological systems in the headwaters of the Tagus River (central Spain), J. Hydrol., 386, 13–26, https://doi.org/10.1016/j.jhydrol.2010.01.001, 2010. a

Martínez-Paz, J. M. and Pellicer-Martínez, F.: Valoración económica del riesgo asociado a las sequías en la Región de Murcia, in: Riesgos Ambientales en la Región de Murcia, edited by: Conesa, C. and Pérez, P., 269–294, EDITUM, Murcia, 2018 (in Spanish). a, b, c

Martínez-Paz, J. M., Perni, A., Ruiz-Campuzano, P., and Pellicer-Martínez, F.: Valoración económica de los fallos de suministro en los regadíos de la cuenca del Segura, Revista Española de Estudios Agrosociales y Pesqueros REEAP, 244, 35–67, 2016 (in Spanish). a, b, c, d

Martínez-Paz, J. M., Banos-González, I., Martínez-Fernández, J., and Esteve-Selma, M. A.: Assessment of management measures for the conservation of traditional irrigated lands: the case of the Huerta of Murcia (Spain), Land Use Policy, 81, 382–391, https://doi.org/10.1016/j.landusepol.2018.10.050, 2018. a

McKenney, D. W., Pedlar, J. H., Papadopol, P., and Hutchinson, M. F.: The development of 1901–2000 historical monthly climate models for Canada and the United States, Agr. Forest Meteorol., 138, 69–81, https://doi.org/10.1016/j.agrformet.2006.03.012, 2006. a

Meza, F. J., Wilks, D. S., Gurovich, L., and Bambach, N.: Impacts of Climate Change on irrigated agriculture in the Maipo basin, Chile: Reliability of water rights and changes in the demand for irrigation, J. Water Res. Plan. Man., 138, 421–430, https://doi.org/10.1061/(ASCE)WR.1943-5452.0000216, 2012. a

Mitchell, T. D. and Hulme, M.: Predicting regional climate change: living with uncertainty, Prog. Phys. Geog., 23, 57–78, https://doi.org/10.1177/030913339902300103, 1999. a

MITECO: Anuario de Aforos de la Red Oficial de Estaciones de Aforos, Ministerio para la Transición Ecológica, available at: https://www.miteco.gob.es/es/agua/temas/evaluacion-de-los-recursos-hidricos/sistema-informacion-anuario-aforos/default.aspx, last access: 15 August 2018. a

Molina-Navarro, E., Trolle, D., Martínez-Perez, S., Sastre-Merlin, A., and Jeppesen, E.: Hydrological and water quality impact assessment of a Mediterranean limno-reservoir under climate change and land use management scenarios, J. Hydrol., 509, 354–366, https://doi.org/10.1016/j.jhydrol.2013.11.053, 2014. a, b

Morán-Tejeda, E., Lorenzo-Lacruz, J., Ignacio López-Moreno, J., Rahman, K., and Beniston, M.: Streamflow timing of mountain rivers in Spain: Recent changes and future projections, J. Hydrol., 517, 1114–1127, https://doi.org/10.1016/j.jhydrol.2014.06.053, 2014. a

Morote, A. F., Olcina, J., and Rico, A. M.: Challenges and proposals for socio-ecological sustainability of the Tagus-Segura Aqueduct (Spain) under Climate Change, Sustainability, 9, 2058, https://doi.org/10.3390/su9112058, 2017. a, b

Nash, J. and Sutcliffe, J.: River flow forecasting through conceptual models part I – A discussion of principles, J. Hydrol., 10, 282–290, https://doi.org/10.1016/0022-1694(70)90255-6, 1970. a

O'Callaghan, J. F. and Mark, D. M.: The extraction of drainage networks from digital elevation data, Comput. Vision Graph., 28, 328–344, 1984. a

Olmstead, S. M.: Climate Change adaptation and water resource management: A review of the literature, Energ. Econ., 46, 500–509, https://doi.org/10.1016/j.eneco.2013.09.005, 2014. a

Onagi, E.: Climate Change and integrated approach to water resource management in the Murray-Darling basin, in: Sustainable Development and Disaster Risk Reduction, edited by: Uitto, J. I. and Shaw, R., 173–187, Springer Japan, Tokyo, https://doi.org/10.1007/978-4-431-55078-5_11, 2016. a

Özdoǧan, M.: Climate change impacts on snow water availability in the Euphrates-Tigris basin, Hydrol. Earth Syst. Sci., 15, 2789–2803, https://doi.org/10.5194/hess-15-2789-2011, 2011 a

Pedro-Monzonís, M., Jiménez-Fernández, P., Solera, A., and Jiménez-Gavilán, P.: The use of AQUATOOL DSS applied to the System of Environmental-Economic Accounting for Water (SEEAW), J. Hydrol., 533, 1–14, https://doi.org/10.1016/j.jhydrol.2015.11.034, 2016a. a, b

Pedro-Monzonís, M., Solera, A., Ferrer, J., Andreu, J., and Estrela, T.: Water accounting for stressed river basins based on water resources management models, Sci. Total Environ., 565, 181–190, https://doi.org/10.1016/j.scitotenv.2016.04.161, 2016b. a, b

Pellicer-Martínez, F. and Martínez-Paz, J. M.: Assessment of interbasin groundwater flows between catchments using a semi-distributed water balance model, J. Hydrol., 519, 1848–1858, https://doi.org/10.1016/j.jhydrol.2014.09.067, 2014. a

Pellicer-Martínez, F. and Martínez-Paz, J. M.: Contrast and transferability of parameters of lumped water balance models in the Segura River Basin (Spain), Water Environ. J., 29, 43–50, https://doi.org/10.1111/wej.12091, 2015. a

Pellicer-Martínez, F., González-Soto, I., and Martínez-Paz, J. M.: Analysis of incorporating groundwater exchanges in hydrological models, Hydrol. Process., 29, 4361–4366, https://doi.org/10.1002/hyp.10586, 2015. a

Perni, A. and Martínez-Paz, J. M.: Measuring conflicts in the management of anthropized ecosystems: Evidence from a choice experiment in a human-created Mediterranean wetland, J. Environ. Manage., 203, 40–50, https://doi.org/10.1016/j.jenvman.2017.07.049, 2017. a

Pulido-Velazquez, D., García-Aróstegui, J. L., Molina, J. L., and Pulido-Velazquez, M.: Assessment of future groundwater recharge in semi-arid regions under climate change scenarios (Serral-Salinas aquifer, SE Spain). Could increased rainfall variability increase the recharge rate?, Hydrol. Process., 29, 828–844, https://doi.org/10.1002/hyp.10191, 2014. a

Pulido-Velazquez, M., Alvarez-Mendiola, E., and Andreu, J.: Design of efficient water pricing policies integrating basinwide resource opportunity costs, J. Water Res. Plan. Man., 139, 583–592, https://doi.org/10.1061/(ASCE)WR.1943-5452.0000262, 2013. a

QGIS Development Team: QGIS Geographic Information System, available at: https://www.osgeo.org/projects/qgis/ (last access: 15 August 2018), 2016. a

R Core Team: R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, available at: https://www.r-project.org/ (last access: 15 August 2018), 2016. a

Shevnina, E., Kourzeneva, E., Kovalenko, V., and Vihma, T.: Assessment of extreme flood events in a changing climate for a long-term planning of socio-economic infrastructure in the Russian Arctic, Hydrol. Earth Syst. Sci., 21, 2559–2578, https://doi.org/10.5194/hess-21-2559-2017, 2017. a

Shrestha, S., Shrestha, M., and Babel, M. S.: Assessment of climate change impact on water diversion strategies of Melamchi Water Supply Project in Nepal, Theor. Appl. Climatol., 128, 311–323, https://doi.org/10.1007/s00704-015-1713-6, 2017. a

Skøien, J., Blöschl, G., Laaha, G., Pebesma, E., Parajka, J., and Viglione, A.: rtop: An R package for interpolation of data with a variable spatial support, with an example from river networks, Comput. Geosci., 67, 180–190, https://doi.org/10.1016/j.cageo.2014.02.009, 2014. a

Smerdon, B. D.: A synopsis of climate change effects on groundwater recharge, J. Hydrol., 555, 125–128, https://doi.org/10.1016/j.jhydrol.2017.09.047, 2017. a

Solaun, K. and Cerdá, E.: The impact of Climate Change on the generation of hydroelectric power-A Case study in Southern Spain, Energies, 10, 1343, https://doi.org/10.3390/en10091343, 2017. a

Spain02v5: Climatic historical data, available at: http://www.meteo.unican.es/datasets/spain02, last access: 18 January 2018.

Szczypta, C., Gascoin, S., Houet, T., Hagolle, O., Dejoux, J.-F., Vigneau, C., and Fanise, P.: Impact of climate and land cover changes on snow cover in a small Pyrenean catchment, J. Hydrol., 521, 84–99, https://doi.org/10.1016/j.jhydrol.2014.11.060, 2015. a

Thomas, H.: Improved Methods for National Water Assessment, Report, contract WR 15249270, U.S. Water Resources Council, available at: https://pubs.er.usgs.gov/publication/70046351 (last access: 15 August 2018), 1981. a, b

Thornthwaite, C. W.: An approach toward a rational classification of climate, Geogr. Rev., 38, 55–94, 1948. a

Tumushabe, J. T.: Climate Change, food security and sustainable development in Africa, in: The Palgrave Handbook of African Politics, Governance and Development, edited by: Oloruntoba, S. O. and Falola, T., 853–868, Palgrave Macmillan US, New York, https://doi.org/10.1057/978-1-349-95232-8_53, 2018. a

Wahba, G.: Spline Models for Observational Data, CBMS-NSF Regional Conference Series in Applied Mathematics, Society for Industrial and Applied Mathematics, Philadelphia, USA, 1990. a

Wang, G. Q., Zhang, J. Y., Xuan, Y. Q., Liu, J. F., Jin, J. L., Bao, Z. X., He, R. M., Liu, C. S., Liu, Y. L., and Yan, X. L.: Simulating the impact of Climate Change on runoff in a typical river catchment of the Loess Plateau, China, J. Hydrometeorol., 14, 1553–1561, https://doi.org/10.1175/JHM-D-12-081.1, 2013. a

Warziniack, T., Lawson, M., and Karen Dante-Wood, S.: Effects of Climate Change on ecosystem services in the Northern Rockies, in: Climate Change and Rocky Mountain Ecosystems, edited by: Halofsky, J. E. and Peterson, D. L., 189–208, Springer International Publishing, Cham, https://doi.org/10.1007/978-3-319-56928-4_10, 2018. a

Weiss, C. E. and Roetzer, G. R.: GeomComb: (Geometric) Forecast Combination Methods, r package version 1.0, available at: https://CRAN.R-project.org/package=GeomComb (last access: 15 August 2018), 2016. a

World Bank Group: High and Dry: Climate Change, Water, and the Economy, World Bank Group, available at: https://openknowledge.worldbank.org/handle/10986/23665 (last access: 15 September 2018), 2016. a

Wurbs, R. A.: Methods for developing naturalized monthly flows at gaged and ungagged sites, J. Hydrol. Eng., 11, 55–64, https://doi.org/10.1061/(ASCE)1084-0699(2006)11:1(55), 2006.  a

Xu, C.-Y. and Singh, V. P.: A Review on monthly water balance models for water resources investigations, Water Resour. Manag., 12, 20–50, https://doi.org/10.1023/A:1007916816469, 1998. a

Xue, X., Zhang, K., Hong, Y., Gourley, J. J., Kellogg, W., McPherson, R. A., Wan, Z., and Austin, B. N.: New multisite cascading calibration approach for hydrological models: Case Study in the Red River Basin using the VIC model, J. Hydrol. Eng., 21, 05015019, https://doi.org/10.1061/(ASCE)HE.1943-5584.0001282, 2016. a

Zare, F., Elsawah, S., Iwanaga, T., Jakeman, A. J., and Pierce, S. A.: Integrated water assessment and modelling: A bibliometric analysis of trends in the water resource sector, J. Hydrol., 552, 765–778, https://doi.org/10.1016/j.jhydrol.2017.07.031, 2017. a

Zhang, E., Yin, X., Xu, Z., and Yang, Z.: Bottom-up quantification of inter-basin water transfer vulnerability to climate change, Ecol. Indic., 92, 195–206, https://doi.org/10.1016/j.ecolind.2017.04.019, 2018. a

Zhang, L., Li, S., Loáiciga, H. A., Zhuang, Y., and Du, Y.: Opportunities and challenges of interbasin water transfers: a literature review with bibliometric analysis, Scientometrics, 105, 279–294, https://doi.org/10.1007/s11192-015-1656-9, 2015. a

Zhang, L. P., Qin, L. L., Yang, Z., Xia, J., and Zeng, S.: Climate change impacts on hydrological processes in the water source area of the Middle Route of the South-to-North Water Diversion Project, Water Int., 37, 564–584, https://doi.org/10.1080/02508060.2012.692108, 2012. a