Exploring the Role of Reservoir Storage in Enhancing Resilience to Climate Change in Southern Europe

: Recent trends suggest that streamﬂow discharge is diminishing in many rivers of Southern Europe and that interannual variability is increasing. This threatens to aggravate water scarcity problems that periodically arise in this region, because both effects will deteriorate the performance of reservoirs, decreasing their reliable yield. Reservoir storage is the key infrastructure to overcome variability and to enhance water availability in semiarid climates. This paper presents an analysis of the role of reservoir storage in preserving water availability under climate change scenarios. The study is focused on 16 major Southern European basins. Potential water availability was calculated in these basins under current condition and for 35 different climatic projections for the period 2070–2100. The results show that the expected reduction of water availability is comparable to the decrease of the mean annual ﬂow in basins with large storage capacity. For basins with small storage, the expected reduction of water availability is larger than the reduction of mean annual ﬂow. Additionally, a sensitivity analysis was carried out by replicating the analysis assuming variable reservoir volumes from 25% to 175% of current storage. The results show that increasing storage capacity attenuates the reduction of water availability and reduces its uncertainty under climate change projections. This feature would allow water managers to develop suitable policies to mitigate the impacts of climate change, thus enhancing the resilience of the system.


Introduction
Climate change, associated with the recorded rise of average temperatures, which are expected to continue increasing to a greater or lesser extent, may also influence other climatic variables such as precipitation, frost, or evapotranspiration [1]. All these changes may affect, in turn, the hydrological processes and consequently net water resources. This threatens the performance of water resource systems and their capability to supply demand and ecological need as presently planned. Therefore, it is necessary to assess both the impact on water resources and the behavior of water systems under such a scenario [2].
Many authors have devoted significant efforts to evaluate net water resources in climate change projections, on all scales from global to basin [3][4][5][6][7][8][9]. Their results show that climate change will affect, in varying ways and to different extents, each region of the planet. As a global result, it could be synthesized that there will be a reduction of water resources of between 10% and 30% [1]. This is an indicative value, useful for developing macro-policies and for raising awareness in the population.
With regard to Southern Europe, despite the dispersion of the various models, the general trend indicates that net resources will decrease, and that the variability of their distribution will increase, as shown in the results of the Prediction of regional scenarios and uncertainties for defining European climate change risks and effects (PRUDENCE) [10] and Climate change and its impacts at seasonal, decadal and centennial timescales (EN-SEMBLES) [11] projects. In Southern Europe, the prognosis is that traditional water scarcity problems will be aggravated. In many basins in this region, the available resources can hardly meet the existing water demands [12]. These are areas with a benign climate, which favors the implementation of agriculture and the development of tourism and services. All these activities require substantial amounts of water. Although these regions have scarce water resources, they are also resilient as they have long experience in dealing with water scarcity and are well adapted to its management [13].
The present study focuses on understanding the effect of reservoir storage capacity on water availability in Southern European basins under climate change. As above mentioned, these basins are typically characterized by scarce and highly variable water resources. The adopted strategy for water resources development in the last century relied on reservoir storage, as it is necessary to store water during the wet periods for its use in the dry ones. In Spain, for example, the existing 1350 large dams helped to increase water availability from 10% to between 40% and 50% of mean natural flow during the last century [14]. As storage capacity grew in parallel with water use, this water availability is used strictly enough to serve current water needs. Alternative adaptation and mitigation measures are being developed in the current century: controlling irrigation water rights, increasing water use efficiency through localized and drip irrigation, developing non-conventional resources, such as water reutilization and desalination, among others [15,16]. Despite these efforts, projections of climate change suggest less water resources with higher variability, which will negatively affect system performance, so water availability is expected to reverse its growing trend [17].
The analysis of reservoir storage capacity and its relationship with safe yield has been a topic of study since the beginning of the development of large hydraulic systems [18]. Initially, graphical methods were developed to determine the reservoir capacity needed to satisfy a given demand with required reliability, and their use was restricted to single reservoir models. Later methods introduced uncertainty of future inflows and attempted to estimate required reservoir size through statistical analysis of inflows, leading to the concepts of risk of failure and reliability. The development of computing allowed the stochastic generation of synthetic series and the disaggregated analysis of multiple reservoirs in a system [19]. Löf and Hardison [20] provided storage-reliability-yield (SRY) relations for assessing the required storage capacity in the USA. The study was later revisited by Vogel et al. [21], who concluded that areas with lower variability tend to be equipped with within-year storage systems while those with large variability required larger over-year storage facilities. Further developments introduced the concepts of resilience and robustness [22] to complete the reliability-yield analysis [23]. An alternative approach is the simulation of the water resources system behavior, which in conjunction with the power of computers allows the development of complex models that reproduce a simulated operation of the system [24][25][26]. These models are useful as decision-support tools for allocating water among users and assessing the effectiveness of structural and managerial actions [27,28] and their capabilities are even extended to groundwater resources and social and economic considerations [29].
Focus is slowly being placed on the impact of climate change on water availability and the role of reservoir storage to increase resilience. Wurbs et al. [30] highlighted the need to introduce climate change in the analysis of water availability and proposed a methodology to couple climatic and system behavior models. Garrote et al. [31] developed a simulation model specifically suited to account for the role of reservoirs in providing water availability in the context of climate change. Several authors [32][33][34] have argued in favor of adaptive reservoir management as an effective mitigation measure during climate change. Adaptive management requires a good knowledge of the interplay between reservoir storage and the reliability, resilience, and vulnerability of a water supply system subject to uncertain input [35]. Water availability deriving from reservoir systems may become increasingly unstable under climate change [36] and knowledge on how regulated water supply systems react to flow alterations is essential for system managers to design climate adaptation policies. This paper looks beyond the impact of climate change on water availability to provide insight into the performance of reservoir storage systems and their effectiveness of adaptation and mitigation measures. With this purpose we include a regional analysis of the performance of reservoir storage and a sensitivity analysis of the reservoir-yield relations under less abundant resources and larger variability conditions in 16 representative European basins. The objective of the research is to check if reservoir storage enhances resilience to climate change. Given the uncertainty of climate projections, the adopted approach is to evaluate basin response under a large ensemble of plausible future scenarios and to evaluate if reservoir storage plays a role in determining the response to changes in hydrologic forcing. System response is quantified in terms of the elasticity of water availability to climate change, comparing changes in potential water availability with changes in mean annual flow. Elasticity is evaluated with the help of two new indices proposed in this work, which characterize the attenuation of changes and the reduction of uncertainty provided by reservoir storage.

Area under Analysis
We present results for 16 major river basins in Europe, which are shown in Figure 1. Basin selection was based on a regional focus on Southern Europe, but including different climates, hydrologic regimes, and storage capacities to allow for a more effective comparison. The selected basins cover a large fraction of the Atlantic and Mediterranean divides of Southern Europe and are representative of the variety of conditions that can be found across the region. The main characteristics of the basins considered in this study are shown in Table 1. Basin areas range from 17,550 km 2 (Segura) to 115,910 km 2 (Loire). Reservoir Storage Volume V includes all reservoirs, except those managed exclusively for hydropower. The basin with largest storage volume is Guadiana, which includes two of the largest reservoirs in Southern Europe: Alqueva (4.15 km 3 ) and La Serena (3.21 km 3 ). Basin hydrology is very variable, with Specific Runoff ranging from 11 mm/year in the Segura basin to 563 mm/year in the Po basin. The most relevant characteristic for this study is specific storage, defined as the ratio of Storage Volume V in km 3 divided by Mean Annual Flow F in km 3 /year for the period 1960-1999. This ratio is usually called Residence Time (in years) and represents the regulation capacity of reservoirs in the basin. In the basins under study, it ranges across three orders of magnitude, from 0.01 years (Arno) to nearly 6 years (Segura).
The spatial support for the analysis is taken from the "Hydro1k" data set [37], derived from the Global 30 Arc-Second Elevation (GTOPO) 30 arc-second digital elevation model of the world. The dataset provides a digital elevation map and a set of topographically derived rasters at 1 km resolution, including streams and drainage basins divided into catchments. The original drainage basins in "Hydro1k" were processed to eliminate catchments which were too small (less than 1000 km 2 ), which were merged to neighboring catchments. The merging was always done with downstream areas and avoiding catchments including reservoirs. The reservoir storage volume in every catchment was obtained from the International Commission on Large Dams (ICOLD) World Register of Dams [38]. We selected dams in the register with more than 0.005 km 3 of storage capacity, excluding dams managed only for hydropower. The reservoirs were georeferenced and linked to the corresponding Hydro1k streams. All dams located in the same Hydro1k subbasin were grouped in an equivalent reservoir adding the storage volume and flooded area (to account for reservoir evaporation losses).  The spatial support for the analysis is taken from the "Hydro1k" data set [37], derived from the Global 30 Arc-Second Elevation (GTOPO) 30 arc-second digital elevation model of the world. The dataset provides a digital elevation map and a set of topographically derived rasters at 1 km resolution, including streams and drainage basins divided into catchments. The original drainage basins in "Hydro1k" were processed to eliminate catchments which were too small (less than 1000 km 2 ), which were merged to neighboring

Methodological Overview
The methodological approach is presented in Figure 2. The analysis is structured in three steps: analysis of the forcing scenarios for water resources systems, analysis of system response in terms of potential water availability, and analysis of the sensitivity of system response to reservoir storage. The analysis of system forcing consists of the compilation of a large name of model runs producing monthly streamflow series in the basins under analysis, both for a historic control period and for a projected future period. Streamflow series for the control period were corrected for bias. Streamflow series for the future period were obtained under different climate scenarios. The scenarios were characterized in terms of the expected changes of mean, standard deviation, and coefficient of variation of annual streamflow. The analysis of system response is focused on the estimation of the potential water availability allowed by current reservoir storage in the basins under analysis. Uncertainty of potential water availability is first characterized for the control and the future periods. Then, the elasticity of water availability to climate changes is explored by comparing changes in potential water availability to changes in mean annual flow, both for individual projections and for the distribution of all projections in each basin. The focus of the analysis is to explore how this elasticity is affected by reservoir storage and streamflow variability. The third step is focused on exploring the sensitivity to reservoir storage. The analyses of the previous step are repeated considering variable storage in each basin. The performance of the system is characterized by two new indices proposed in this study: the attenuation index and the uncertainty index. These indices describe how the performance of the water supply system is affected by changes in streamflow. The main conclusions of the study are obtained by comparing how these indices change as a function of reservoir storage for all basins.
three steps: analysis of the forcing scenarios for water resources systems, analysis of system response in terms of potential water availability, and analysis of the sensitivity of system response to reservoir storage. The analysis of system forcing consists of the compilation of a large name of model runs producing monthly streamflow series in the basins under analysis, both for a historic control period and for a projected future period. Streamflow series for the control period were corrected for bias. Streamflow series for the future period were obtained under different climate scenarios. The scenarios were characterized in terms of the expected changes of mean, standard deviation, and coefficient of variation of annual streamflow. The analysis of system response is focused on the estimation of the potential water availability allowed by current reservoir storage in the basins under analysis. Uncertainty of potential water availability is first characterized for the control and the future periods. Then, the elasticity of water availability to climate changes is explored by comparing changes in potential water availability to changes in mean annual flow, both for individual projections and for the distribution of all projections in each basin. The focus of the analysis is to explore how this elasticity is affected by reservoir storage and streamflow variability. The third step is focused on exploring the sensitivity to reservoir storage. The analyses of the previous step are repeated considering variable storage in each basin. The performance of the system is characterized by two new indices proposed in this study: the attenuation index and the uncertainty index. These indices describe how the performance of the water supply system is affected by changes in streamflow. The main conclusions of the study are obtained by comparing how these indices change as a function of reservoir storage for all basins.

Current and Future Runoff Scenarios
The focus of the present study is the analysis of the role of reservoir storage to determine how water resources systems react to changes in hydrologic forcing. An effort was made to obtain a wide ensemble of scenarios that would represent the uncertainty linked to climate projections. Therefore, we chose to combine model results obtained under two sets of emission scenarios, the Special Report on Emission Scenarios (SRES) and Representative Concentration Pathways (RCP), in order to increase the size of the ensemble. Current and future runoff scenarios were compiled from three previous studies that include Southern Europe [39][40][41]. These studies were based on results from different climate models developed over the last 15 years under two sets of emission scenarios: SRES and RCP. They jointly describe the uncertainty that is currently challenging water managers.
The first set of scenarios was taken from the output of regional climate models from the PRUDENCE project [10]. The study by González-Zeas et al. [39] was based on the projections of surface runoff made by eight RCMs at 50 km resolution nested in a single global model, referred to as HadAM3H, in emission scenarios A2 and B2. They analyzed current ) and future (2070-2100) time slices. The second set of scenarios was based on the results of the Regional Climate Models (RCMs) of the ENSEMBLES project [11]. The project produced many transient model runs for the time period from 1960 to 2100 using RCMs to characterize model uncertainty. The study by Garrote et al. [40] selected runoff output from four ENSEMBLES models at 25 km resolution under emission scenario A1B to study the major Mediterranean river basins of Europe. They worked with windows of analysis on the transient model runs in three time slices: historical (1960-1990), short term (2020-2050) and long term (2070-2100). In the first two sets of scenarios, monthly runoff time series were directly obtained from the "Total runoff" variable (mrro) produced by RCMs. The values of surface runoff flux available at the nodes of the native grid of the RCMs (50 km resolution in PRUDENCE and 25 km resolution in ENSEMBLES) were used to produce monthly runoff maps by interpolation at the finer grid provided by the Hydro1k dataset (1 km). The center of the RCM grid was taken as a point equal to the average for that cell. Interpolation was based on a weighted mean using the inverse of the distance squared as weight. These runoff maps were combined with the subbasin definitions of Hydro1k to obtain monthly streamflow values for each subbasin. The third set of scenarios was based on the results of the global hydrological model PCRaster GLOBal Water Balance (PCRGLOBWB) model [42] in the Inter-Sectorial Impact Model Intercomparison Project (ISIMIP) [43]. In ISIMIP, the PCRGLOBWB model was forced with five global climate models under historical conditions and climate change projections corresponding to four Representative Concentration Pathways scenarios: RCP-2., RCP-4., RCP-6. and RCP-8., corresponding to radiative forcing in the year 2100 of 2.6, 4.5, 6.0, and 8.5 W/m 2 , respectively. The study by Sordo-Ward et al. [41] used naturalized streamflow from PCRGLOBWB at 50 km resolution to analyze 1261 subbasins covering the entire territory of Western Europe. They considered two time slices in their analysis: historical (1960-1999) and long-term projection (2060-2099). The monthly streamflow time series in the subbasins were also obtained from monthly runoff maps derived from the runoff produced from the PCRGLOBWB model through interpolation at the Hydro1k 1 km grid.
A total of 16 model runs were compiled for the historical period (eight model runs from the PRUDENCE project, three model runs from the ENSEMBLES projects and five model runs from the PCRGLOBWB model). The windows of analysis in this period overlap for years . All these model runs produced different results in the basins under analysis. To assess the quality of these hydrological projections, the results obtained at the working scale of each model run were compared to a reference estimate of mean annual runoff under current conditions. The selected reference was the annual surface runoff layer (Global Composite Runoff Fields) of the University of New Hampshire Global Runoff Data Centre (GRDC) [44]. This data layer was produced by combining a database of observed river discharge information in more than 9900 gauging stations with a climate-driven water balance model to develop consistent runoff fields. The combination of direct readings from gauging stations with the water balance model preserves the spatial distribution of runoff generation and provides the best estimate of observed runoff over large domains. The mean values of the time series compiled for the historical period were compared with mean annual runoff produced by GRDC. The results are presented in Figure 3, which shows the scatterplot resulting from comparing catchment mean annual runoff produced from GRDC with that produced by model runs for the historical period. Model runs corresponding to the Special Report on Emissions Scenarios (SRES) (PRUDENCE and ENSEMBLES projects) show poor agreement. The models that performed best were Universidad de Castilla La Mancha (UCM) and Eidgenössische Technische Hochschule Zürich (ETHZ2), with coefficients of determination slightly lower than 0.6. This poor performance can be explained because runoff was obtained directly from RCM output. Model runs for the RCP scenarios (ISIMIP project) were produced by the hydrological model PCRGLOBWB. They show better performance, with coefficients of determination close to 0.7, but they reveal significant bias for low runoff. The discrepancies obtained in the comparison suggest that bias correction is necessary to overcome this very large model uncertainty. Using the monthly series of individual models without bias correction would imply significant distortion in the regulation provided by reservoirs in each basin. The ratio between reservoir storage capacity and mean annual flow would change for each model run, affecting the evaluation of the regulation capacity provided by the reservoirs. For this reason, runoff derived from RCM results and from PCRGLOBWB was corrected for bias in each location. The chosen method for bias correction was linear scaling [45]. This method is justified by data availability, because GRDC only provides monthly long-term means of runoff. Therefore, all model projections for the historical period have the same mean.
hydrological model PCRGLOBWB. They show better performance, with coefficients of determination close to 0.7, but they reveal significant bias for low runoff. The discrepancies obtained in the comparison suggest that bias correction is necessary to overcome this very large model uncertainty. Using the monthly series of individual models without bias correction would imply significant distortion in the regulation provided by reservoirs in each basin. The ratio between reservoir storage capacity and mean annual flow would change for each model run, affecting the evaluation of the regulation capacity provided by the reservoirs. For this reason, runoff derived from RCM results and from PCR-GLOBWB was corrected for bias in each location. The chosen method for bias correction was linear scaling [45]. This method is justified by data availability, because GRDC only provides monthly long-term means of runoff. Therefore, all model projections for the historical period have the same mean.  The number of model runs compiled for the long-term climate change projection was 35: eight model runs corresponding to the A2 scenario, four model runs corresponding to the B2 scenario, three model runs for the A1B scenario, five model runs for RCP-2 scenario, five model runs for RCP-4 scenario, five model runs for RCP-6 scenario and five model runs for RCP-8 scenario. The windows of analysis in the long-term projection overlap for the years 2070-2099. These projections were corrected for model bias by applying the same correction as in the corresponding model in the historical period. This ensemble of climate projections was put together from different projects developed over a 15-year period, running a range of global climate models under two sets of emission scenarios, and applying different methodologies. It can thus be considered a representative description of the range of scenarios that climate change science is projecting for the region. However, it should be noted that runoff projections derived from climate models are uncertain. Climate models provide a good overall representation of climate, but their performance degrades at the scale of individual grid boxes, indicating that they are not skillful at their smallest scale. The performance of RCMs generally improves after suitably removing bias. However, model errors still remain large, particularly for climatic variables relevant for hydrology, like precipitation or runoff [46]. Given this inherent uncertainty, a basic hypothesis of this work is that water management decisions based on the global analysis of a wide range of projections produces better results than decisions based on a very detailed analysis of a reduced number of projections.
The average annual runoff obtained from GRDC in the period 1960-2000 was also used to characterize the basins under analysis. The relationship between Specific Runoff and reservoir Residence Time is plotted in Figure 4 for all Hydro1k basins in Southern Europe, highlighting the 16 basins under analysis. As can be seen in Figure 4, there is a clear relation between both variables, with larger values of storage corresponding to basins with lower values of specific runoff. The selected basins produce a good coverage of the possible range of behaviors found in the region, from basins with large specific water resources and low storage capacity like Arno, Po or Loire, to others in the opposite situation with very low water resources and large storage volumes as Júcar, Guadiana or Segura.
it should be noted that runoff projections derived from climate models are uncertain. Climate models provide a good overall representation of climate, but their performance degrades at the scale of individual grid boxes, indicating that they are not skillful at their smallest scale. The performance of RCMs generally improves after suitably removing bias. However, model errors still remain large, particularly for climatic variables relevant for hydrology, like precipitation or runoff [46]. Given this inherent uncertainty, a basic hypothesis of this work is that water management decisions based on the global analysis of a wide range of projections produces better results than decisions based on a very detailed analysis of a reduced number of projections.
The average annual runoff obtained from GRDC in the period 1960-2000 was also used to characterize the basins under analysis. The relationship between Specific Runoff and reservoir Residence Time is plotted in Figure 4 for all Hydro1k basins in Southern Europe, highlighting the 16 basins under analysis. As can be seen in Figure 4, there is a clear relation between both variables, with larger values of storage corresponding to basins with lower values of specific runoff. The selected basins produce a good coverage of the possible range of behaviors found in the region, from basins with large specific water resources and low storage capacity like Arno, Po or Loire, to others in the opposite situation with very low water resources and large storage volumes as Júcar, Guadiana or Segura.

Water Availability Analysis
The study is based on the analysis of how climate change affects water availability in the different basins, and how this affect is modified by available reservoir storage. Potential Water Availability (PWA) is defined as the annual water demand that can be satisfied in a point of the drainage network with a given reliability. PWA depends on the mean and variability of the streamflow series, the storage available for flow regulation, the monthly distribution of the demand and the reliability indicator adopted in the analysis. In this study, PWA was estimated with the Water Availability and Adaptation Policy Analysis (WAAPA) model [31,47]. WAAPA simulates the operation of a complex water resources system with many reservoirs. The basic topological unit of WAAPA is the river network. The main components are inflows, reservoirs and demands, all linked to nodes in the network. WAAPA computes the amount of water supplied to demands from a system of reservoirs accounting for ecological flows and evaporation losses. Input data for WAAPA are monthly inflows in relevant points of the river network, monthly demand values, and reservoir data. Reservoirs are described by monthly maximum and minimum capacity, storage-area relationship, monthly rates of evaporation, and monthly required environmental flow. WAAPA applies an algorithm with simple operating rules, where all reservoirs in the basin are jointly managed to satisfy the set of demands, drawing water preferably from reservoirs located upstream. This algorithm is applied to potential demands located in every node in the river network, and therefore water availability is obtained for the entire river network. The main results of WAAPA are time series of monthly volumes supplied to each demand, monthly storage values and monthly values of spills, environmental flows, and evaporation losses in every reservoir. From this output, demand reliability can be computed for the criterion of choice (volume reliability, time reliability at the monthly or annual scale, or more complex criteria).
WAAPA can obtain PWA for a given demand reliability criterion through an iterative scheme that changes local demand values until the reliability criterion is met with a given precision. In this study, PWA is estimated by considering only one type of demand in the system, with constant monthly distribution. This choice was made because the true monthly distribution of demands in each model node is unknown. Results therefore should be considered only approximate and could be fine-tuned if the ratio between urban and irrigation demand was known in every model node. Ecological flows were specified as the 10% percentile of the monthly marginal distribution of natural flows. System performance is evaluated as gross volume reliability. PWA is obtained for 92% volume reliability. This reliability level was chosen as an intermediate value between reliabilities required from urban demands (usually close to 100%) and those required from irrigation demand (usually close to 90%), assuming an approximate distribution of 20% urban demand and 80% irrigation demand, which is typical of Portugal, Spain, and Greece [48].

Results and Discussion
The WAAPA model was run for the European Mediterranean region for the 16 hydrologic scenarios corresponding to the historical period  and for the 35 hydrologic scenarios corresponding to climate projections for the long-term time horizon 2070-2100. The long-term time horizon was chosen for two reasons. Firstly, the results from PRU-DENCE project were only available for this time horizon. Secondly, the changes in the long-term time horizon are usually more accentuated than in the mid-term time horizon and the effects are more apparent. Results were obtained for all catchments in the Hydro1k dataset, but, for the sake of simplicity, we only present global results for the 16 basins under analysis. We first analyze the climate projections, then we present the results obtained for PWA in the basins. Average values of these results are summarized in Table 2 and presented and discussed in detail in the following section. Finally, the role of storage is studied through a sensitivity analysis.

Climate Projections
We first present the characterization of climate projections for the basins under study. Climate projections were taken from the runoff variable of RCM models in the PRUDENCE and ENSEMBLES projects (under SRES emission scenarios) and of the PCRGBLOBWB hydrologic model (under RCP emission scenarios). Mean and coefficient of variation of annual flows were computed for each basin during the historical period and during the long-term projection. Changes in the long-term projection were estimated taking the control period as a reference, applying the following expressions: where F is Mean Annual Flow, SD is the Standard Deviation of the annual time series of streamflow, and CV is the Coefficient of Variation of the annual time series of streamflow (standard deviation of the annual time series divided by mean annual flow). The subindices HIST and PROJ refer to the historical period and to the long-term projection. The results are depicted in Figure 5,  Figure 5 shows negative correlation between changes in Mean Annual Flow ∆F and Coefficient of Variation ∆CV: reduction of F and increase of CV. The general shape of the scatter plot is similar in all basins in Southern Europe. This has clear implications for water management since both factors will negatively impact water availability. This tendency is stronger for basins with larger residence times that, as seen in Figure 4, are located in water scarce regions, already facing strong hydrologic irregularities. The dispersion of results is stronger for basins with larger residence times, presenting an additional challenge for water management. The ensemble of projections, jointly considered, suggests that water managers should be ready to cope with less abundant and more variable water resources in the future. Given the large dispersion of results, water managers should also be ready to deal with greater year-on-year variability or extreme events than in the past. Figure 5 also shows that expected changes in CV are much larger than changes in F, with many basins reaching extreme values close to 2 (a 100% increase). The basins showing more extreme projections are Guadalquivir, Júcar and Guadiana.

Water Availability
The WAAPA model was used to compute Potential Water Availability (PWA) for the historical period and for the long-term projection in the 16 basins under analysis. The re sults are shown in Figure 6, which presents the value of PWA obtained in each basin as a function of the relative rank of the corresponding projection. All 35 projections were used to prepare this figure, thus mixing projections under SRES and RCP emission scenarios An individual plot of each basin is included in the supplementary material, where the joint distribution is compared to the distributions of both sets of emission scenarios. The corresponding emission scenario is identified for each model run available in the long term projection. These plots show that there is no clear correlation between the emission scenario and the projected PWA. Values corresponding to different emission scenarios are mixed and the most extreme scenarios (A2 and RCP-8) do not always produce the mini mum values for PWA.
PWA is expressed as a fraction of Mean Annual Flow (F) in the historical period. Re sults for the historical period are shown in the upper chart (a) and results for the long term projection are shown in the lower chart (b). If all model runs were assumed equi probable, this plot would correspond to the empirical estimation of the probability distri bution function of PWA expected in each basin. The results show that the relative value of PWA to F tends to be larger for basins with larger storage capacity, both in the histor ical and in the projection periods. This fact clearly illustrates the effectiveness of reservoi storage to increase water availability. The plots also show large uncertainty in the estima tion of PWA. For the historical period, this result is remarkable because historical time series were corrected for bias with respect to the GRDC estimation of F and therefore al had the same Mean Annual Flow. The uncertainty in PWA reflects model uncertainty be cause the differences in PWA can only be attributed to the differences in the seasonal and interannual variability of the time series produced by each model run. Unfortunately, the skill of the models to reproduce current hydrological irregularity cannot be evaluated be cause there are no available regional data sets for Southern Europe on interannual natu

Water Availability
The WAAPA model was used to compute Potential Water Availability (PWA) for the historical period and for the long-term projection in the 16 basins under analysis. The results are shown in Figure 6, which presents the value of PWA obtained in each basin as a function of the relative rank of the corresponding projection. All 35 projections were used to prepare this figure, thus mixing projections under SRES and RCP emission scenarios. An individual plot of each basin is included in the Supplementary Materials, where the joint distribution is compared to the distributions of both sets of emission scenarios. The corresponding emission scenario is identified for each model run available in the long-term projection. These plots show that there is no clear correlation between the emission scenario and the projected PWA. Values corresponding to different emission scenarios are mixed and the most extreme scenarios (A2 and RCP-8) do not always produce the minimum values for PWA.
PWA is expressed as a fraction of Mean Annual Flow (F) in the historical period. Results for the historical period are shown in the upper chart (a) and results for the long-term projection are shown in the lower chart (b). If all model runs were assumed equiprobable, this plot would correspond to the empirical estimation of the probability distribution function of PWA expected in each basin. The results show that the relative value of PWA to F tends to be larger for basins with larger storage capacity, both in the historical and in the projection periods. This fact clearly illustrates the effectiveness of reservoir storage to increase water availability. The plots also show large uncertainty in the estimation of PWA. For the historical period, this result is remarkable because historical time series were corrected for bias with respect to the GRDC estimation of F and therefore all had the same Mean Annual Flow. The uncertainty in PWA reflects model uncertainty because the differences in PWA can only be attributed to the differences in the seasonal and interannual variability of the time series produced by each model run. Unfortunately, the skill of the models to reproduce current hydrological irregularity cannot be evaluated because there are no available regional data sets for Southern Europe on interannual naturalized streamflow variability. Except in the Arno basin, PWA is expected to decrease significantly in the long-term projection with respect to the historical period, with average reductions between 15% and 35%. These reductions are the consequence of reduced F and increased CV. The most significant reductions are projected for the basins of South Western Europe: Guadiana and Guadalquivir (35% on average) and Duero (28% on average). The uncertainty of PWA in the long-term projection is larger than that in the historical period due to the additional variability introduced by emission scenarios. However, the large model uncertainty hinders the interpretation of results obtained for different emission scenarios.
The estimated changes in PWA are compared to estimated changes in F in Figures  7 and 8. Figure 7    Except in the Arno basin, PWA is expected to decrease significantly in the long-term projection with respect to the historical period, with average reductions between 15% and 35%. These reductions are the consequence of reduced F and increased CV. The most significant reductions are projected for the basins of South Western Europe: Guadiana and Guadalquivir (35% on average) and Duero (28% on average). The uncertainty of PWA in the long-term projection is larger than that in the historical period due to the additional variability introduced by emission scenarios. However, the large model uncertainty hinders the interpretation of results obtained for different emission scenarios.
The estimated changes in PWA are compared to estimated changes in F in  In Figures 7 and 8, the changes of F are estimated from the first expression shown in Equation (1). The changes of PWA are similarly estimated from the comparison of values obtained in the long-term projection and the control period for the same model, applying the following expression: where PWA is Potential Water Availability and the sub-indices HIST and PROJ refer to the historical period and to the long-term projection. The mean values plotted in chart (b) of Figure 7 reveal that basins with small storage capacity (Po, Loire, Tiber, Garonne and Rhône) show a larger reduction of PWA than the reduction of F. The Arno basin is the exception, with no reduction of PWA despite a small average reduction of F. The basins with larger storage capacity, in general, show a smaller reduction of PWA than the reduction of F. Struma-Strymon, Vardar-Axios, Ebro, Maritsa-Evros, Guadalquivir, Tajo-Tejo, Júcar, Guadiana and Segura belong to this group. Duero-Douro is an exception, with larger reduction of PWA than of F. This may be explained because Duero-Douro shows the largest difference between change in F and change in SD. The wide scatter of changes in F and PWA in chart (a) of Figure 7 shows that there is no exact relation between changes in Mean Annual Flow (∆F) and changes in Potential Water Availability (∆PWA). For individual projections, changes in PWA may be larger, equal or smaller than changes in F. This is, in part, a consequence of changes in hydrologic variability, which may explain why negative changes in F produce positive changes in PWA and vice versa. However, the comparison of Figure 5 and Figure 7 shows that changes in hydrologic variability alone cannot explain the diversity of behaviors seen in Figure 7. Hydrologic variability is measured in terms of Coefficient of Variation of annual flows and is therefore referred to interannual variability. Basins with small storage capacity show a behavior more exposed to changes in CV because, for them, water availability is almost directly determined by short-duration dry periods of the streamflow series. These dry periods show a large variability among model runs, which explains the variability observed in values of PWA. As basin storage grows larger, the reservoirs attenuate the effect of short-duration dry periods and the interannual variability becomes A remarkable effect shown in Figure 8 is that the uncertainty regarding changes in F (gray line) seems to grow as specific storage grows. The larger spreads of the estimated probability distributions appear in basins with larger specific storage, in the bottom row. Basins with comparatively smaller storage capacity, in the first row, show much less uncertainty on changes in F. This shows that reservoir storage was developed where it was required: in basins with large hydrologic variability. Furthermore, the difference in uncertainty between changes in F and changes in PWA, which is large in basins with small storage, is progressively reduced as specific storage increases. obtained in the long-term projection and the control period for the same model, applying the following expression: where PWA is Potential Water Availability and the sub-indices HIST and PROJ refer to the historical period and to the long-term projection.
The mean values plotted in chart (b) of Figure 7 reveal that basins with small storage capacity (Po, Loire, Tiber, Garonne and Rhône) show a larger reduction of PWA than the reduction of F. The Arno basin is the exception, with no reduction of PWA despite a small average reduction of F. The basins with larger storage capacity, in general, show a smaller reduction of PWA than the reduction of F. Struma-Strymon, Vardar-Axios, Ebro, Maritsa-Evros, Guadalquivir, Tajo-Tejo, Júcar, Guadiana and Segura belong to this group. Duero-Douro is an exception, with larger reduction of PWA than of F. This may be explained because Duero-Douro shows the largest difference between change in F and change in SD. The wide scatter of changes in F and PWA in chart (a) of Figure 7 shows that there is no exact relation between changes in Mean Annual Flow (∆F) and changes in Potential Water Availability (∆PWA). For individual projections, changes in PWA may be larger, equal or smaller than changes in F. This is, in part, a consequence of changes in hydrologic variability, which may explain why negative changes in F produce positive changes in PWA and vice versa. However, the comparison of Figures 5 and 7 shows that changes in hydrologic variability alone cannot explain the diversity of behaviors seen in Figure 7. Hydrologic variability is measured in terms of Coefficient of Variation of annual flows and is therefore referred to interannual variability. Basins with small storage capacity show a behavior more exposed to changes in CV because, for them, water availability is almost directly determined by short-duration dry periods of the streamflow series. These dry periods show a large variability among model runs, which explains the variability observed in values of PWA. As basin storage grows larger, the reservoirs attenuate the effect of short-duration dry periods and the interannual variability becomes less important. Basins with storage capacity larger than mean annual flow show a much less sensitivity to changes in the coefficient of variation of mean annual flows. Figure 8 is useful to analyze the effect of the uncertainty on emission scenarios. The estimated probability distributions shown in Figure 8 reveal a wide range of behaviors. The basins were classified in five groups (A1, A2, A3, B1 and B2), according to the relative value of the distributions of changes in F and PWA. Group A1 is integrated by basins where the distribution of expected reductions in PWA is to the left of the distribution of expected reductions in F, suggesting that the availability of reservoir storage tends to dampen the effect of climate change. Struma-Strymon, Vardar-Axios, Guadalquivir and Tajo belong to this group. In the second group, A2, the distributions of expected changes in F and PWA are very similar. This group is formed by Ebro, Júcar and Segura. The only basin in Group A3 (Douro-Duero) presents larger expected reductions in PWA than in F. In group B, the probability distributions of F and PWA cross each other. In group B1, the distribution of changes in PWA is to the left of the distribution of changes in F for low probability values. For high probability values, the distribution of changes in PWA is to the right of the distribution of changes in F. This results in larger uncertainty for changes in PWA than in F. This effect may be due to increased exposure to changes in variability due to lack of regulation storage. Arno, Po, Loire, Tiber, Garonne, Rhône and Maritsa-Evros belong to group B1. The only basin in group B2 is Guadiana, where the distribution of changes in PWA is to the right of the distribution of changes in F for low probability values and to the left for high probability values. Guadiana shows less uncertainty in changes of PWA than in changes in F, due to its large reservoir storage.
A remarkable effect shown in Figure 8 is that the uncertainty regarding changes in F (gray line) seems to grow as specific storage grows. The larger spreads of the estimated probability distributions appear in basins with larger specific storage, in the bottom row.
Basins with comparatively smaller storage capacity, in the first row, show much less uncertainty on changes in F. This shows that reservoir storage was developed where it was required: in basins with large hydrologic variability. Furthermore, the difference in uncertainty between changes in F and changes in PWA, which is large in basins with small storage, is progressively reduced as specific storage increases.

The Influence of Storage
The results obtained in Section 2.2 suggest that reservoir storage plays a relevant role in controlling how projected changes in Mean Annual Flow may be translated into changes in Potential Water Availability. However, the large variability of local conditions in the studied basins introduces uncertainties in the analysis. In this section we further explore the influence of reservoir storage on changes in water availability through a sensitivity analysis that discounts for local conditions. We repeated the water availability analysis but considering different storage volumes in each basin. Potential Water Availability was computed in current and future scenarios in the 16 basins, assuming changing reservoir volumes of 25%, 50%, 75%, 100%, 125%, 150% and 175% of current storage. Storage was proportionally reduced or increased in the same location of existing reservoirs. This choice was made for convenience, without any implications for projected future evolution of storage in the region. In fact, the most likely scenario in the future for European basins is a progressive reduction of available storage due to reservoir sedimentation, with very little additional storage being built. Figure 9 shows the scatter plots of changes in  Figure 9 reveal that the dispersion of the scatter plot gets reduced as the storage capacity is increased. This effect is more marked for the Guadiana basin, which has the largest reservoir storage.
In order to assess the global behavior, the values obtained for changes in F and PWA were classified according to the relative rank of the corresponding projection, obtaining an empirical estimate of their probability distributions, under the assumption that all scenarios analyzed are equally likely. The results are shown in Figure 10, which presents the estimate of the probability distribution of changes in Mean Annual Flow (∆F, blue line) and changes in Potential Water Availability (∆PWA) for different storage values in colored lines from brown (25% of current storage volume) to green (175% of current reservoir storage). The basins analyzed showed variable sensitivity to storage. Some basins, like Arno, Ebro, Maritsa-Evros or Guadiana, showed very little sensitivity to storage capacity because the distributions of expected changes of PWA are very similar. Other basins, like Po, Tiber, Rhône or Struma-Strymon, presented significant differences in behavior depending on the storage volume assumed. In some of these basins, the estimated probability distributions of changes in PWA for high storage values (green color) were located to the right, indicating less reductions of PWA. Po, Loire, Rhône, Duero-Douro, Ebro, Júcar, Guadiana and Segura show this behavior. For other basins, however, the probability distributions for high storage values were located to the left. Tiber, Struma-Strymon and Guadalquivir belong to this group.
This range of behaviors illustrates the complex relations between hydrologic variability and reservoir storage in determining water availability in climate change scenarios, suggesting that a specific analysis for local conditions is required to translate projections of changes in mean annual flow into projections of changes in water availability in basins with significant storage capacity. RCP scenarios are represented as circles. Results shown in Figure 9 reveal that the dispersion of the scatter plot gets reduced as the storage capacity is increased. This effect is more marked for the Guadiana basin, which has the largest reservoir storage. In order to assess the global behavior, the values obtained for changes in F and PWA were classified according to the relative rank of the corresponding projection, obtaining an empirical estimate of their probability distributions, under the assumption that all scenarios analyzed are equally likely. The results are shown in Figure 10, which presents the estimate of the probability distribution of changes in Mean Annual Flow (∆F, blue line) and changes in Potential Water Availability (∆PWA) for different storage values in colored lines from brown (25% of current storage volume) to green (175% of current reservoir storage). The basins analyzed showed variable sensitivity to storage. Some basins, like This range of behaviors illustrates the complex relations between hydrologic var bility and reservoir storage in determining water availability in climate change scenari suggesting that a specific analysis for local conditions is required to translate projectio of changes in mean annual flow into projections of changes in water availability in bas with significant storage capacity. The influence of reservoir storage on the elasticity of water availability was analyz by computing an attenuation index I A , defined in the following expression: where ∆PWA is the change in Potential Water Availability and ∆F is the change in Me Annual Flow. As seen in the previous sections, most changes in PWA and F are redu tions and therefore ∆PWA and ∆F are negative. A positive value of this index indica an attenuation of the effect of climate change: the absolute value of ∆PWA is smaller th the absolute value ∆F.
We explore how the attenuation index I A changes with reservoir storage. The resu are shown in Figure 11, which presents the value of the attenuation index as a function reservoir storage in each basin for all available projections (thin grey lines), the avera values (solid lines in the color code corresponding to the basin) and average values pl and minus one standard deviation (dotted lines in the color code corresponding to t The influence of reservoir storage on the elasticity of water availability was analyzed by computing an attenuation index I A , defined in the following expression: where ∆PWA is the change in Potential Water Availability and ∆F is the change in Mean Annual Flow. As seen in the previous sections, most changes in PWA and F are reductions and therefore ∆PWA and ∆F are negative. A positive value of this index indicates an attenuation of the effect of climate change: the absolute value of ∆PWA is smaller than the absolute value ∆F. We explore how the attenuation index I A changes with reservoir storage. The results are shown in Figure 11, which presents the value of the attenuation index as a function of reservoir storage in each basin for all available projections (thin grey lines), the average values (solid lines in the color code corresponding to the basin) and average values plus and minus one standard deviation (dotted lines in the color code corresponding to the basin). The results show a large variability for individual projections, which translates into large uncertainty for water managers. The variability of the attenuation index appears to be progressively reduced as specific storage grows across basin locations (from top row to bottom row). This suggests that reservoir storage plays a relevant role in reducing uncertainty on the effects of climate change projections on water availability. basin). The results show a large variability for individual projections, which translates into large uncertainty for water managers. The variability of the attenuation index appears to be progressively reduced as specific storage grows across basin locations (from top row to bottom row). This suggests that reservoir storage plays a relevant role in reducing uncertainty on the effects of climate change projections on water availability. The effect of reservoir storage on the variability of the attenuation index I A is further explored by analyzing the uncertainty index I U , defined as: where σ(∆PWA) is the standard deviation of the changes in Potential Water Availability for all projections and σ(∆F) is the standard deviation of the changes in Mean Annual Flow for all projections. I U index compares the variability of the projections of changes in water availability to that of the projections of mean flow. A negative value of this index indicates a reduction of the uncertainty of climate change projections: the variability of ∆PWA is smaller than the variability ∆F.
The summary of results found in the analysis of the I A and I U indices is shown in Figure 12. Chart (a) of Figure 12 compares the average of the values of the I A index obtained in the sensitivity analyses of reservoir storage for all basins. Chart (b) of Figure 12 represents the corresponding values of the I U index. The effect of reservoir storage on the variability of the attenuation index I A is further explored by analyzing the uncertainty index I U , defined as: where σ(∆PWA) is the standard deviation of the changes in Potential Water Availability for all projections and σ(∆F) is the standard deviation of the changes in Mean Annual Flow for all projections. I U index compares the variability of the projections of changes in water availability to that of the projections of mean flow. A negative value of this index indicates a reduction of the uncertainty of climate change projections: the variability of ∆PWA is smaller than the variability ∆F.
The summary of results found in the analysis of the I A and I U indices is shown in Figure 12. Chart (a) of Figure 12 compares the average of the values of the I A index obtained in the sensitivity analyses of reservoir storage for all basins. Chart (b) of Figure 12 represents the corresponding values of the I U index.
The plots shown in Figure 12 indicate that increased reservoir storage results in larger values of the attenuation index I A in most basins and smaller values of the uncertainty index I U in all basins. The results for the attenuation index are less conclusive than those for the uncertainty index. Out of the 16 basins analyzed, I A is observed to decrease with increasing reservoir storage in four basins: Tiber, Garonne, Struma-Strymon and Guadalquivir. In the case of Garonne and Struma-Strymon, the decrease only covers the range from 25% to 100% of current reservoir storage. For storage volumes between 100% and 175% of current reservoir storage I A is increasing with increasing reservoir storage. The case of Tiber basin may be explained because reservoirs only cover 12% of the contributing area and therefore 88% of the flow is unregulated. Guadalquivir basin is exposed to the most extreme reduction of Mean Annual Flow (43% on average) and this may have an influence on the observed behavior. In the case of the uncertainty index, the reduction with increasing storage is observed for all basins. These results are valid for the range of storage volumes explored in the sensitivity analysis in each basin individually and for all basins as a whole, regardless of basin size and location in Southern Europe, and therefore show a clear picture of the role played by reservoir storage in attenuating the impact of reduced streamflow on water availability and on reducing the uncertainty of climate change projections. It is unlikely that reservoir storage will be further increased in Southern Europe. Most basins already have an adequate amount of storage and additional storage capacity would not increase water availability in a scenario of decreasing resources. However, water managers should be aware that proper management of currently available storage will be helpful to address the challenges posed by climate change. The plots shown in Figure 12 indicate that increased reservoir storage results in larger values of the attenuation index I A in most basins and smaller values of the uncertainty index I U in all basins. The results for the attenuation index are less conclusive than those for the uncertainty index. Out of the 16 basins analyzed, I A is observed to decrease with increasing reservoir storage in four basins: Tiber, Garonne, Struma-Strymon and Guadalquivir. In the case of Garonne and Struma-Strymon, the decrease only covers the range from 25% to 100% of current reservoir storage. For storage volumes between 100% and 175% of current reservoir storage I A is increasing with increasing reservoir storage. The case of Tiber basin may be explained because reservoirs only cover 12% of the contributing area and therefore 88% of the flow is unregulated. Guadalquivir basin is exposed to the most extreme reduction of Mean Annual Flow (43% on average) and this may have an influence on the observed behavior. In the case of the uncertainty index, the reduction with increasing storage is observed for all basins. These results are valid for the range of storage volumes explored in the sensitivity analysis in each basin individually and for all basins as a whole, regardless of basin size and location in Southern Europe, and therefore show a clear picture of the role played by reservoir storage in attenuating the impact of reduced streamflow on water availability and on reducing the uncertainty of climate change projections. It is unlikely that reservoir storage will be further increased in Southern Europe. Most basins already have an adequate amount of storage and additional storage capacity would not increase water availability in a scenario of decreasing resources. However, water managers should be aware that proper management of currently available storage will be helpful to address the challenges posed by climate change.

Conclusions
Projected changes in hydrologic regime and water availability were analyzed in 16

Conclusions
Projected changes in hydrologic regime and water availability were analyzed in 16 basins in Southern Europe applying the WAAPA model to streamflow time series obtained from 35 climate projections under 7 emission scenarios. The analysis of climate projections concluded that a significant reduction of mean annual flow can be expected in most basins. The reduction in the mean is supplemented by a strong increase in the coefficient of variation, due to an increase of the variability of the projected series. This analysis is uncertain due to the very large variability introduced by the different models and emission scenarios examined. The overall result implies a corresponding reduction in potential water availability, with variable results across basins depending on hydrologic regime and reservoir storage. Basins with large storage values showed reductions of water availability comparable to the reductions of mean annual flow. Basins with small storage capacity showed a larger reduction of water availability than the reduction of mean annual flow. Although model and emission scenario uncertainties are larger than the expected reduction of water availability, a consistent picture emerges from the joint analysis of all projections, requiring significant adaptation measures to compensate for the projected reduction of water availability.
The influence of reservoir storage on basin response to climate change was studied through a sensitivity analysis where current reservoir storage was modified to examine its effects on water availability with values ranging from 25% to 175% of current storage values. The results showed very large variability, which illustrates the complex interplay between hydrologic regime, reservoir storage and water availability. Two indices were introduced to clarify the overall behavior: the attenuation index and the uncertainty index. The attenuation index compares the changes in water availability to the changes in mean annual flow. Positive values of this index indicate an attenuation of the impact of climate change projection on water availability. The uncertainty index compares the variability of changes in water availability and in mean annual flow. Positive values of this index indicate a reduction of the uncertainty of climate change projections on water availability. The results of the sensitivity analysis showed that increasing reservoir storage attenuates the reduction of water availability and reduces the uncertainty of climate projection. The results are valid for each individual basin within the range of storage volumes examined and for the set of 16 Southern European basins analyzed in this work. The effect observed for reservoir storage is a positive factor for system managers since decisions become harder as uncertainty grows. This feature would allow water managers to develop suitable policies to mitigate the impacts of climate change, thus enhancing the resilience of the system.