Water-Energy Nexus for an Italian Storage Hydropower Plant under Multiple Drivers

: Climate change has repercussions on the management of water resources. Particularly, changes in precipitation and temperature impact hydropower generation and revenue by a ﬀ ecting seasonal electricity prices and streamﬂow. This issue exempliﬁes the impact of climate change on the water-energy-nexus, which has raised serious concern. This paper investigates the impact of climate change on hydropower with a multidisciplinary approach. A holistic perspective should be favored as the issue is complex, consequently, we chose to investigate a speciﬁc case study in Italy. It allows grasping the details, which matters in mountainous area. We integrated a hydrological model, hydropower management model, nine climate scenarios, and ﬁve electricity scenarios for a speciﬁc storage hydropower plant. Independently from the scenarios, the results show a glacier volume shrinkage upward of 40% by 2031 and minimum of 50% by 2046. The reservoir mitigates losses of revenue that reach 8% in the worst case, however, are lower compared with run-of-the-river conﬁguration. Changes in price seasonality amplitude also determine modiﬁcations in revenues, while temporal shifts appear to be ine ﬀ ective. For run-of-the-river, any variation in hydrological cycle immediately translates into revenue. Comparing the results of all future scenarios with the base scenario, it can be concluded that an increase in temperature will slightly improve the performances of hydropower.


Introduction
Hydropower has a great importance in global electricity generation, representing the 15% of the worldwide amount [1], and accounting for the 85% of the total production [2] in the world's renewable energy framework. This technology has the great advantage of ensuring operational flexibility and enhancing the possibility of storing energy and managing the production. It increases the reliability of the electric system, acting as an ancillary service. This is a fundamental characteristic especially with the increasing development of intermittent sources like wind and photovoltaic [3][4][5][6].
For these reasons, increasing attention has been lately posed on the quantification of the effects of climate changes on hydropower [7], considering increase in temperature or change in temporal and spatial patterns of liquid and solid precipitation [8,9]. Due to its strong dependency on hydrologic regime, hydropower is among the most vulnerable technologies to climate changes [10]. In addition, as the production of this power source relies not only on inflows but also on electricity prices [11], it becomes important to understand how climate changes will affect these two drivers.
Precipitation and glacier volume will change. In mountainous areas, one can expect a progressive decrease in peak snow accumulation [12][13][14]. Laternser and Schneebeli [12] have shown a decrement

Italian Hydropower: Current Situation and Climate Change Impact Projections
With 18.8 GW installed hydropower capacity in 2017, Italy is ranked 11th at the global level and 4th at the European scale [48]. A total of 53.3 GW of electricity is from renewables accounting to 35% of Italian energy mix, where hydropower accounts for 35% of the renewable energy produced [49,50]. Hydropower is strategic since Italy conceals a minimum amount of fossil feedstock and depends on resources from abroad [51]. As of 2016, the import of fossil fuels accounted for 92.6% [52]. In the Italian framework, the National Electricity Strategy (Strategia Elettrica Nazionale, in Italian) [53] and the recent Integrated National Plan for Energy and Climate (Piano Nazionale Integrato per l'Energia e il Clima, in Italian) [54] fix future national energy policies, depicting future pathways in renewable sources sector. They target a progressive phase out of generation from coal by 2025, and 55.4% of national electricity consumption from renewable energy sources (RES) by 2030 [55]. To achieve this goal, new incentivizing policies must double or even triple wind and solar energy at 2030, throughout new constructions as well as revamping and repowering existing installations. It will ask for increasing storage capacity to manage intermittent supply. This will empower the sharing of existing and the building of new pump-and-storage power plants. It nowadays represents the most mature storage technologies, despite new incomers are becoming competitive [56]. In terms of energy generation, the exploitable hydropower potential is becoming scare. It, nevertheless, represents a strategic asset in baseload and peak-load generation. Therefore, the Italian Ministry of Economic Development (Ministero dello Sviluppo Economico) aims at incentivizing the revamping of existing "great hydropower" rather than the construction of new hydropower, with the exception of only "small hydropower plants" (i.e., SHP) [53,54].
The 4268 hydropower plants are heterogeneously spread throughout the country [49]. Small installations are numerous: 3074 (≤1 MW), 886 (1-10 MW). However, the remaining 308 (>10 MW) plants contribute to 75% of the total hydroelectric power generated [49]. Therefore, Northern part of Italy regroups 81.8% of the installations and 59.5% of installed capacity (see Table 1). Following the latest report of the Italian Regulatory Authority [57] the share of hydropower managed by major companies is as such, Enel (37.7%), A2A (9.9%), Erg (3.2%), Edison (5.3%), CVA (6.3%), Hydro Dolomiti Energia (6%), Alperia (4.6%), SEL (4.2%), Iren (2.9%), and others (20.8%). Whereas in central and southern Italy, Abruzzo and Calabria have the maximum hydropower concentration in terms of energy generation. The streamflow regime of Italian hydropower is typically of two types: (a) snow and/or glacier dominated (low flow in winter while high flow in spring and summer); (b) rainfall dominated (depends on the seasonality of rainfall, typically low flow during summer) [18,26]. With the anticipated changes in climate due to global warming, substantial changes are expected for streamflow of snow and/or glacier dominated hydropower catchment [7,58]. With the retreat of glaciers a temporary increase in streamflow is expected for the near future and is completely dependent on the characteristic of the catchment and rate of climate change; with a possibility of change in the regime of hydropower catchment from glacier to snow or snow to rainfall based [18,59]. For mountain catchments, particularly in the Alpine region these shifts in hydrological regime is due to the early onset of melting season and shrinkage of glaciers, resulting in higher streamflow in early summer and spring [18,20,[60][61][62].
In the recent work of Smiraglia et al. [63], the area changes in Italian glaciers in the New Italian Glacier Inventory [64] are compared with both the CGI-CNR Inventory [65] and WGI (World Glacier Inventory) dataset. They observed a change of −30% and −39% respectively in the last five decades. It is therefore necessary to understand the influence of climate change in the Italian hydropower sector. Table 2 gives some insights on the studies conducted in this decade for possible consequences of climate change on Italian hydropower.
Assessment and quantification of climate change impacts on hydropower are quite complex in nature [7,66]. From the existing literature reported in Table 2, it is clear that the implications are site specific depending on the catchment and use of hydropower apart from power generation (e.g., drinking, recreational, irrigation). The evaluation of hydropower at both a regional and national level is quite tiresome due to the need for detailed information such as meteorological variables, streamflow, design specification, inter-basin transfers, water usage among others. Therefore, most of the studies reported in Table 2 evaluate a single hydropower site. From these studies, it was observed for Italian hydropower a reduction in late summer flow and anticipation of peak streamflow for snow dominated catchments. For the comprehensive overview of projected changes refer to Table 2.

Methodology and Case Study
This section presents the models and dataset utilized to build the analysis on a case study, along with the workflow adopted. The hydrological and hydropower management model are described along with the future scenarios of climate and electricity price used for simulations. Finally, the case study will be presented, describing in detail the main features of the basin considered and the data available.

Study Area and Data
This paper focuses on the Sabbione catchment, located in Val d'Ossola, in the northern part of Piedmont region of Italy (see Figure 1a). The main feature of this catchment along with the hydropower plant considered as the closing section is presented in Table 3. The delimitation of the catchment was done considering the Sabbione dam (Morasco hydropower plant) as closure. This catchment has altitude range between 2460-3358 m above sea level (m.a.s.l), comprising four glaciers that cover approximately the 30% of the whole area of the catchment [64]. Recent work by Giaccone et al. [84] observes a glacier retreat in this catchment of 2 km in the time interval 1885-2011.

Methodology and Case Study
This section presents the models and dataset utilized to build the analysis on a case study, along with the workflow adopted. The hydrological and hydropower management model are described along with the future scenarios of climate and electricity price used for simulations. Finally, the case study will be presented, describing in detail the main features of the basin considered and the data available.

Study Area and Data
This paper focuses on the Sabbione catchment, located in Val d'Ossola, in the northern part of Piedmont region of Italy (see Figure 1a). The main feature of this catchment along with the hydropower plant considered as the closing section is presented in Table 3. The delimitation of the catchment was done considering the Sabbione dam (Morasco hydropower plant) as closure. This catchment has altitude range between 2460-3358 m above sea level (m.a.s.l), comprising four glaciers that cover approximately the 30% of the whole area of the catchment [64]. Recent work by Giaccone et al. [84] observes a glacier retreat in this catchment of 2 km in the time interval 1885-2011. The hydropower plants originally built in the glacier valley used to encompass the front of the glacier submerged in Morasco hydropower's reservoir before the glacier retreat. The Sabbione reservoir (see Table 3) is the biggest artificial reservoir of the Ossola and Piedmont Region, the second for dimensions in the whole alpine chain behind Place Moulin in Val d'Aosta [18]. Enel manages this reservoir and dam and exploits the water from snow, precipitation, and glacier fed Toce river basin. This catchment has been studied previously in the literature [11,18] predicting the hydropower future of this basin in terms of power production and glacier coverage only. However, the impact of climate change seasonality on the revenue was neglected. The reservoir is used for hydropower production only, since irrigation is not present and domestic water usage is negligible.  Table 3. Main features of the Sabbione catchment and Morasco power plant [18]. The hydropower plants originally built in the glacier valley used to encompass the front of the glacier submerged in Morasco hydropower's reservoir before the glacier retreat. The Sabbione reservoir (see Table 3) is the biggest artificial reservoir of the Ossola and Piedmont Region, the second for dimensions in the whole alpine chain behind Place Moulin in Val d'Aosta [18]. Enel manages this reservoir and dam and exploits the water from snow, precipitation, and glacier fed Toce river basin. This catchment has been studied previously in the literature [11,18] predicting the hydropower future of this basin in terms of power production and glacier coverage only. However, the impact of climate change seasonality on the revenue was neglected. The reservoir is used for hydropower production only, since irrigation is not present and domestic water usage is negligible.

The Hydrological Model
We used a semi-distributed hydrological model following Bongio et al. [25], Patro et al. [18], and Ranzani et al. [59]. The basin was subdivided into independent sub-domain elevation bands with homogeneous geomorphological characteristics. The model simulates the average daily discharge flowing to the closing section, i.e., the Sabbione hydropower plant, from the topographic information, the local precipitation, and temperature data series. The total discharge is calculated as the sum of four distinct contributions: snow melt, glacial melt, rainfall runoff, and base/groundwater flow. The nivometer and the rain gauge approach [85] separate solid from liquid precipitation, while the degree-day model [86,87] treated snow and glacier melt. The mass balance of the glaciers was calculated at yearly scale considering mass-conservation law and neglecting any motion associated. We computed evapotranspiration with a long-term average temperature criterion [88,89], and effective net precipitation and infiltration with the Curve Number method [90]. Finally, a Nash approach computes all the components at the closing section [91]. A schematic representation of the model can be seen in Figure 2. Rosso [92] formula estimated the Nash parameters for the basin response to liquid precipitation. These are expressed as monomial functions of hortonian ratios RB, RL, RA and the time scale factor LΩ/V.
where LΩ (km) indicates the length of the main stream, V (m/s) the average streamflow velocity of a Rosso [92] formula estimated the Nash parameters for the basin response to liquid precipitation. These are expressed as monomial functions of hortonian ratios R B , R L , R A and the time scale factor L Ω /V.
where L Ω (km) indicates the length of the main stream, V (m/s) the average streamflow velocity of a flood wave in the hydrographic network. From the literature, V assumes values between 2 and 3 m/s, in this work it was assumed equal to 2.5 m/s [92]. Hortonian ratios were evaluated from the quantitative geomorphological description of the hydrographic network, using the Horton-Strahler hierarchical model. The estimation of all the parameters was done using GIS software, to obtain the Horton-Strahler classification for the Sabbione catchment. The results found and the parameters obtained are shown in Table 4. In order to evaluate the maximum and minimum degree day parameters, DDF MAX and DDF MIN respectively, were determined using processing-modelling as reported in Avanzi et al. [93]. We calibrated the model with the meteorological dataset provided by the Regional Agency for Environmental Protection (ARPA Piemonte). The available daily meteorological data include precipitation, temperature, and snow depth from 2000 to 2014 for the Formazza meteorological station located at an altitude of 2453 m.a.s.l. The calibration was done minimizing the Root Mean Square Error (RMSE) (mm w. e.), obtained by comparing simulated and real values of snow water equivalent (SWE) for all the series available.
The values obtained from this phase were DDF MAX = 1.77 mm/d • C and DDF MIN = 0.99 mm/d • C.
The values obtained from this phase were DDFMAX = 1.77 mm/d °C and DDFMIN = 0.99 mm/d °C.  Concerning the temperature lapse rate (TLR), it was calibrated comparing the temperatures recorded by some stations at high altitude in Piedmont region with the ones computed with the following formulation: Concerning the temperature lapse rate (TLR), it was calibrated comparing the temperatures recorded by some stations at high altitude in Piedmont region with the ones computed with the following formulation: where T higher elevation is the temperature observed at the meteorological station located at higher altitude while T lower elevation is the temperature observed at the lower altitude meteorological station, and dH is the difference in the altitude of the two stations under consideration. Therefore, starting from temperature series recorded by lower altitude stations, TLR was computed minimizing the error between simulated and observed ones. The stations of each pair belong to the same catchment and are available in ARPA Piemonte website (http://www.arpa.piemonte.it/) is as follow: The calibration obtained a value of TLR equal to 0.53 • C/100 m, consistent with the values found by Rolland [94]. With all this information, the hydrologic budget was computed for the catchment at daily resolution. This hydrological modelling framework has been tested previously for 42 Italian Alpine catchments of which Sabbione was also considered [18]. The model has proven to provide a sufficiently good performance for the scope of this paper with an average RMSE daily and monthly scale as 2.56 and 1.32 m 3 /s, and average Nash Sutcliffe efficiency (NSE) daily and monthly scales at −0.14 and 0.63. Further information including details of the model and protocol of calibration and validation can be found in Patro et al. [18].

Future Climate Scenarios
The future discharge inflow at Sabbione reservoir was computed according to the future scenarios following Patro et al. [18]. Nine scenarios consider direct assumptions on future climatic patterns of two variables: temperature and precipitation. These "local" scenarios impose trends in the historical series of the climatic variables. The adoption of this approach stems from the fact that (a) several uncertainties affect Representative Concentration Pathways for mountainous areas [7]; (b) it allows a direct evaluation of the uncertainty with a Monte Carlo technique compared to IPCC-like outputs [59]. Bongio et al. [25] analyzed these nine scenarios with IPCC-like scenario, proving comparable results. Each scenario simulates 100 time series of discharges to obtain a 95% confidence interval for the time horizon of 1 January 2016-31 December 2046. These nine scenarios were subdivided into four major categories based on variation in precipitation and temperature time series: • Future-like-present scenarios (scenarios 1 and 2): They evaluate the future evolution of the glaciers at present conditions. Temperature and precipitation assume presents patterns, alternatively considering or neglecting the computation of the mass balance of the glaciers in the basin. Liquid-precipitation scenario (scenario 6). This extreme scenario assumes liquid-only precipitation and neglects temperature warming.

•
The mixed scenarios (scenarios 7-9): These complex scenarios built random series combining liquid-only precipitation and the three above-mentioned trends of temperature.

The Hydropower Management Model
We simulated the energy generation schedule with the model proposed in Ranzani et al. [59]. It maximizes the objective function provided in Equation (4). It provides the hydropower plant operator revenue at an annual scale and hourly time granularity. It also adds the value of the residual water (R T ) to avoid emptying the reservoir at the end of each simulation. Figure 2 shows the general structure of the management model.
where ρ is water density, g is the acceleration of gravity, c is the water flow through the turbines, η is the plant efficiency, ∆t is the time multiplier that corresponds to the unity of time considered for the operations (1 h in the present work), T is the time horizon for the optimization (1 year in this case), h t is the hydraulic head, P t is the electricity day-ahead price, x t is the scheduling multiplier, and R T is the value of the residual volume of water in the reservoir. The latter represents the virtual revenue that would have been generated with the residual water.
The process of management must respect physical constraints linked to the capacity of the reservoir and of the turbines. Therefore, the overall optimization problem can be formulated as: The objective function relies on the parameter x t , regulating the amount of discharge released through the penstock. It embodies the scheduling of the turbines as it prescribes the amount of water to be turbined per unit of time. Φ is the function of the reservoir's height-volume (water level-volume) curve. The model outputs energy generation schedule, revenue, and reservoir water volume behavior.
The model considers a local search algorithm, called Threshold Accepting (TA) [96][97][98] to approximate the objective function maximum. This technique explores the search space by affecting the turbine schedule, so that the choice of a particular setting in the model is mainly linked to the fact that it will permit to easily reach the observed income with low computational weight. It is recorded if the value of the objective function is better. The algorithm accepts a degradation in the objective function if lower to a threshold, which reduces progressively in the search. This approach avoids staying stuck in a local optimum, while staying efficient.
In the absence of data of actual power generated, the model was not validated with real data. Therefore, the configuration found in Ranzani et al. [59] was assumed, as the reliability of the model was demonstrated. We only adapted the parameters of the algorithm to fit with the specific constraints of the Sabbione hydropower installation.

Electricity Price Scenarios
The high uncertainty of long-term electricity price projections determined the necessity of introducing a specific approach. Five electricity price scenarios were taken from Gaudard et al. [26]. Figure 4 shows the average monthly electricity spot price series for Italy (1 January 2005-31 December 2016), taken from Energy Markets Manager (Gestore dei Mercati Energetici, in Italian, see http://www.mercatoelettrico.org/). It shows a double peak, one in winter (heating systems) and the other in summer (cooling system).
introducing a specific approach. Five electricity price scenarios were taken from Gaudard et al. [26]. Figure 4 shows the average monthly electricity spot price series for Italy (1 January 2005-31 December 2016), taken from Energy Markets Manager (Gestore dei Mercati Energetici, in Italian, see http://www.mercatoelettrico.org/). It shows a double peak, one in winter (heating systems) and the other in summer (cooling system). We identified the main periodicity in series via a Fourier transform analysis. The first two orders (annual and biannual cycles) put in evidence the historical shape of the periodicity for Italian prices. From this historical periodicity, five scenarios project the evolution of seasonality pattern of electricity spot prices in Italy. Table 5 details the scenarios considered where the annual mean price is constant to focus on the impact of seasonality. The hydropower management model takes into account the streamflow simulated for the future climate and the price trajectories to provide hydropower plant revenue. We identified the main periodicity in series via a Fourier transform analysis. The first two orders (annual and biannual cycles) put in evidence the historical shape of the periodicity for Italian prices. From this historical periodicity, five scenarios project the evolution of seasonality pattern of electricity spot prices in Italy. Table 5 details the scenarios considered where the annual mean price is constant to focus on the impact of seasonality. The hydropower management model takes into account the streamflow simulated for the future climate and the price trajectories to provide hydropower plant revenue. Shifted Seasonality SS Seasonality is shifted onward by three months. This means that the prices will peak in spring and autumn instead of winter and summer.

Climate Change Effects on Glacier and Flow Regimes
We launched the hydrological model from 2016 to 2046. Each climate scenario considered 100 simulations and extracted the 95% confidence interval for each annual discharge. Figure 5 shows the evolution of glacier's volume over the 30 years.
The rate of glacier retreat is not constant. Until 2030, volumes reduce of almost 50% irrespective of the climate scenario considered. This fast decrease is due to the melting of low altitude glaciers (under 3000 m.a.s.l.), which are more vulnerable to temperature increases [18]. Then, the reduction slows down to reach a new pseudo equilibrium condition. Note that the glacier's mass balance model neglects the dynamic mobility i.e., the transfer of the upper glacier's mass to lower one. This could reduce the decrease of volume at lower altitudes and increase the reduction at higher ones [19]. Figure 5 also shows that the glacier's volume reduces significantly in liquid only precipitation scenarios (6)(7)(8)(9) due to the absence of snow that no longer fuels the glacier. The modifications in the hydrological regime play a relevant role in glacier retreat.

Seasonality
IS The seasonality is shifted by six months Shifted Seasonality SS Seasonality is shifted onward by three months. This means that the prices will peak in spring and autumn instead of winter and summer.

Climate Change Effects on Glacier and Flow Regimes
We launched the hydrological model from 2016 to 2046. Each climate scenario considered 100 simulations and extracted the 95% confidence interval for each annual discharge. Figure 5 shows the evolution of glacier's volume over the 30 years. The rate of glacier retreat is not constant. Until 2030, volumes reduce of almost 50% irrespective of the climate scenario considered. This fast decrease is due to the melting of low altitude glaciers (under 3000 m.a.s.l.), which are more vulnerable to temperature increases [18]. Then, the reduction slows down to reach a new pseudo equilibrium condition. Note that the glacier's mass balance model neglects the dynamic mobility i.e., the transfer of the upper glacier's mass to lower one. This could reduce the decrease of volume at lower altitudes and increase the reduction at higher ones [19]. Figure  5 also shows that the glacier's volume reduces significantly in liquid only precipitation scenarios (6)(7)(8)(9) due to the absence of snow that no longer fuels the glacier. The modifications in the hydrological regime play a relevant role in glacier retreat.  Figure 6 represents the average monthly discharge for the nine climate scenarios in 2030. It also indicates the variation of annual runoff volumes. Presence or absence of glaciers significantly influences the average monthly discharge. The climate scenario 1 (base scenario) neglects the mass balance of glacier, thus observes higher discharge in spring and summer, i.e., when temperature increasing causes the ablation period. In autumn and winter, precipitation governs the discharge. As for climate scenario 2 having meteorological condition as base scenario, the summer peak will decrease due to glacier shrinkage, while the periodicity remains unchanged at 2030. Concerning climate scenarios 3-5, warming affects the hydrological regime with an anticipation of summer peaks due to earlier onset of the ablation of snowpack and glacier. Interestingly, the basin becomes more impulsive and almost equally distributed in the liquid only precipitation scenarios (climate scenarios 6-9). In these scenarios, the curves are almost in phase, differing only in the mean value of discharge. Figures 5 and 6 combined show the glaciers depletion predominance towards the discharge. The decreasing rate of discharge is higher for climate scenario 2 and progressively reduces as higher trends in temperature (climate scenarios 3-5) are implemented, something that occurs also for climate scenarios 6-9. In fact, the greater temperatures induce a faster regression of the glacier, thus increasing runoff. This compensation would continue until the volume of the glacier is enough to balance the losses due to higher temperatures. The results also reflected those found by Ravazzani et al. [11] on the Toce river basin, which comprises our study area, and differ slightly in the percentages of variation of monthly mean discharge. Ravazzani et al. [11] observed a seasonal shift in monthly discharge with a significant increase in the winter period, while the monthly discharge decreases in summer for the simulation period 2011-2050. decrease due to glacier shrinkage, while the periodicity remains unchanged at 2030. Concerning climate scenarios 3-5, warming affects the hydrological regime with an anticipation of summer peaks due to earlier onset of the ablation of snowpack and glacier. Interestingly, the basin becomes more impulsive and almost equally distributed in the liquid only precipitation scenarios (climate scenarios [6][7][8][9]. In these scenarios, the curves are almost in phase, differing only in the mean value of discharge. Figure 6. Average monthly discharge obtained with a 95% confidence interval at 2030 for all scenarios. The mean annual volume variation relative to climate scenario 1 is presented in brackets. Figures 5 and 6 combined show the glaciers depletion predominance towards the discharge. The decreasing rate of discharge is higher for climate scenario 2 and progressively reduces as higher trends in temperature (climate scenarios 3-5) are implemented, something that occurs also for climate scenarios 6-9. In fact, the greater temperatures induce a faster regression of the glacier, thus increasing runoff. This compensation would continue until the volume of the glacier is enough to balance the losses due to higher temperatures. The results also reflected those found by Ravazzani et al. [11] on the Toce river basin, which comprises our study area, and differ slightly in the percentages of variation of monthly mean discharge. Ravazzani et al. [11] observed a seasonal shift in monthly discharge with a significant increase in the winter period, while the monthly discharge decreases in summer for the simulation period 2011-2050.

Climate Change Impact on Reservoir Volumes and Revenue
We simulated the hydropower management with the nine climate scenarios discharge series. It provided the reservoir volumes and annual revenue. Figure 7 shows the reservoir's volume obtained with the historical price setting. The brackets give the relative variations of incomes to base climate scenario 1 in 2030. Please note that the management model tends to keep the upper level of reservoir higher, to exploit the hydraulic head effect. As already observed in Figure 5, the glacier depletion results in a reduction of discharge into the reservoir. However, seasonality matters. In the climate scenarios 1-5, the operator fills the reservoir in spring and summer, when water discharge is high.   Figure 6. Average monthly discharge obtained with a 95% confidence interval at 2030 for all scenarios. The mean annual volume variation relative to climate scenario 1 is presented in brackets.

Climate Change Impact on Reservoir Volumes and Revenue
We simulated the hydropower management with the nine climate scenarios discharge series. It provided the reservoir volumes and annual revenue. Figure 7 shows the reservoir's volume obtained with the historical price setting. The brackets give the relative variations of incomes to base climate scenario 1 in 2030. Please note that the management model tends to keep the upper level of reservoir higher, to exploit the hydraulic head effect. As already observed in Figure 5, the glacier depletion results in a reduction of discharge into the reservoir. However, seasonality matters. In the climate scenarios 1-5, the operator fills the reservoir in spring and summer, when water discharge is high. The warmer the climate is, the earlier the reservoir is filled. The energy is generated when the prices are higher from summer to winter. The reservoir volume largely relies on water seasonality. The warmer the climate is, the earlier the reservoir is filled. The energy is generated when the prices are higher from summer to winter. The reservoir volume largely relies on water seasonality. In contrast, electricity prices largely determine the reservoir volume shape in the climate scenarios 6-9, where the water seasonality is stable. The reservoir is mostly filled in spring when prices are especially low. The water discharge still affects the reservoir volumes, reflecting similarities found in hydrological regimes. Note that the reservoir used volume is much lower. A smaller dam can be built with low amplitude of water seasonality. Climate change should be considered when designing new projects.
Smart management can mitigate the impact of climate change. Water volumes in Figure 6 vary more than income in Figure 7. The management lowers the negative impact thanks to head effect. In addition, the operator keeps generating electricity during super-peak hours. The only loss in revenue is during low profitability hours. On the other hand, the positive impacts are also tempered in all but scenario 7. Reservoir hydropower can smooth climate change impacts in terms of revenue. This result completely aligns with the finding of Ravazzani et al. [11], where they also concluded that optimal regulation of the Toce river basin hydropower for simulation period 2011-2050, anticipates the date of achieving maximum reservoir volume with drawdown in reservoir volume in August and September to capture the increasing autumn inflows.

Impact of Electricity Price Scenarios
We now analyze the impact of the five different electricity price scenarios. We consider the business-like-present hydrological scenario (climate scenario 1) to focus on the impact of the seasonality of the electric demand. Figure 8 shows the evolution of reservoir's volumes according to 10   In contrast, electricity prices largely determine the reservoir volume shape in the climate scenarios 6-9, where the water seasonality is stable. The reservoir is mostly filled in spring when prices are especially low. The water discharge still affects the reservoir volumes, reflecting similarities found in hydrological regimes. Note that the reservoir used volume is much lower. A smaller dam can be built with low amplitude of water seasonality. Climate change should be considered when designing new projects.
Smart management can mitigate the impact of climate change. Water volumes in Figure 6 vary more than income in Figure 7. The management lowers the negative impact thanks to head effect. In addition, the operator keeps generating electricity during super-peak hours. The only loss in revenue is during low profitability hours. On the other hand, the positive impacts are also tempered in all but scenario 7. Reservoir hydropower can smooth climate change impacts in terms of revenue. This result completely aligns with the finding of Ravazzani et al. [11], where they also concluded that optimal regulation of the Toce river basin hydropower for simulation period 2011-2050, anticipates the date of achieving maximum reservoir volume with drawdown in reservoir volume in August and September to capture the increasing autumn inflows.

Impact of Electricity Price Scenarios
We now analyze the impact of the five different electricity price scenarios. We consider the business-like-present hydrological scenario (climate scenario 1) to focus on the impact of the seasonality of the electric demand. Figure 8 shows the evolution of reservoir's volumes according to the five price scenarios. Price seasonality affects the reservoir management. Doubling the amplitude of seasonality (DS) means spring prices are becoming less attractive. The level is flat, has almost no water incomes, and no electricity is being generated. In contrast, summer prices are high, thus the operator generates most of the energy at that period. In the case of inverted seasonality (IS), peak of water income in summer is almost only stored as prices are low. The operator mostly starts generating in winter with a shift of seasonality onward by three months. Finally, without seasonality, the reservoir level is flatter. Energy is generated throughout the year. The reservoir is used to capture the peak hours in the 12 months, as the water mostly intakes in summer. The correlation between water and price seasonality determines reservoir level and optimal design.
The reservoir is large enough to manage a shift of seasonality. IS and SS do not observe a significant variation of revenue (in brackets). In contrast, doubling or removing the amplitude of seasonality affects the revenue. DS and WS keep the mean annual price stable but affect the peak Price seasonality affects the reservoir management. Doubling the amplitude of seasonality (DS) means spring prices are becoming less attractive. The level is flat, has almost no water incomes, and no electricity is being generated. In contrast, summer prices are high, thus the operator generates most of the energy at that period. In the case of inverted seasonality (IS), peak of water income in summer is almost only stored as prices are low. The operator mostly starts generating in winter with a shift of seasonality onward by three months. Finally, without seasonality, the reservoir level is flatter. Energy is generated throughout the year. The reservoir is used to capture the peak hours in the 12 months, as the water mostly intakes in summer. The correlation between water and price seasonality determines reservoir level and optimal design.
The reservoir is large enough to manage a shift of seasonality. IS and SS do not observe a significant variation of revenue (in brackets). In contrast, doubling or removing the amplitude of seasonality affects the revenue. DS and WS keep the mean annual price stable but affect the peak prices. When doubling the amplitude of seasonality, peak prices increase while baseload diminishes. As Sabbione installation can generate during peak price and store during baseload, its revenue significantly increases. Without seasonality, prices are stable over the year and peak prices are smoothed, thus the revenue decreases.

Impact of Mixed Scenarios
We now combine two climate scenarios (2 and 9) with the five price scenarios to compute the normalized incomes with respect to 2015. We selected these climate scenarios because of their significant impact on the reservoir volume and income (seen in Section 4.1). They serve as a general representation for all other climate scenarios. Figure 9 represents the reservoir water volumes obtained from mixed scenarios while Figure 10 outlines the normalized income with respect to the BaU price scenario of climate scenario 1. It can be seen in Figures 9a and 10a, the mixed scenario 2-BaU, to compensate the effect of temperature increase the maximum emptying of the reservoir is regulated to take the advantage of head effect, and it is also in phase with the reference business as the present baseline scenario in dashed lines. While, the mixed scenario 2-DS has a trend similar to that observed in Section 4.3, less storage in winter to accommodate for electricity generation when prices are higher. In the mixed scenario 2-WS, the hydropower generation is aligned to the incoming discharge, resulting in emptying of reservoir in winter while filling in the summer to maximize the revenue. While, in Figures 9b and 10b, it is clear that due to strong temperature increment we have maximum incoming discharge during summer due to glacier ablation with anticipation of the summer peaks. Moreover, in mixed scenario 9-BaU, since the seasonality of prices remains constant (double peak in winter and summer) and the seasonality of incoming discharge is negligible, filling happens exactly in the intermediate seasons and the emptying takes place right in summer peak to maximize the revenue. For mixed scenario 9-DS, the management model takes the advantage of peak both in price and discharge, and reservoir volumes are also higher compared to other mixed scenarios. It is fascinating to note that for the mixed scenario 9-WS, the management model operates the hydropower plant as a run-of-the-river configuration, immediately utilizing the incoming discharge for electricity generation, since prices have no seasonality. Figure 10 highlights the fact that change in one scenario (either climate or price) has a limited impact on the management of the reservoir, while combining two scenarios (climate + price), the situation starts to become complex. It therefore highlights that hydropower should be considered with the nexus perspective. It can be seen in Figure 9a and Figure 10a, the mixed scenario 2-BaU, to compensate the effect of temperature increase the maximum emptying of the reservoir is regulated to take the advantage of head effect, and it is also in phase with the reference business as the present baseline scenario in dashed lines. While, the mixed scenario 2-DS has a trend similar to that observed in Section 4.3, less storage in winter to accommodate for electricity generation when prices are higher. In the mixed

What If the Reservoir Configuration Becomes a Run of River Model?
A run-of-river power plant design allows grasping the mitigation benefit of a reservoir. Energy generation only depends on the seasonality of inflows, as the operator cannot delay water usage to peak price periods. We simulated this configuration assuming that hydrological runoff is immediately used or released if higher than the turbine capacity. If Q F represents the incoming discharge from the river basin and Q P the maximum possible discharge through the penstocks, discharge Q t considered for power production is characterized as follows: We assumed the maximum water head with the reservoir, otherwise it would affect the turbine capacity. It can be viewed as a limit condition to which sedimentation has filled the dam reservoir. Climate change and its impact on glaciers retreat may exacerbate the sedimentation phenomenon [99]. For instance, Ceppo Morelli dam in Italy had an reservoir volume of 0.5 Mm 3 , reduced to zero nowadays [25,100]. It runs as a run-of-river hydropower plant. Table 6 shows a decrease in the incomes for a RoR configuration compared to the reservoir one, because of a non-optimal management of the production in all climate scenarios. It represents a difference of around −2% for climate scenarios 2-5 and around −4% in the other climate scenarios. The shape of historical electricity prices (Figure 4) obviously impacts the variation. The reduction of discharge in summer (where the prices are higher) and an anticipation of the ablation season in spring (when the prices are lower), have an impact on revenues. The hydrological regime tends to shift and go out of phase with prices. It can be concluded that a run of river power plant revenue is dependent not only on the amount of water income, but also on the behavior of prices. Indeed, considering a different electricity market, e.g., the Swiss one, the effects can be different [59]. Anyway, it proves the mitigation benefits of the reservoir. Table 6. Comparison of the income variation between run-of-the-river (RoR) and reservoir configuration in the case of climate scenarios at 2030 w.r.t. baseline climate scenario 1 over historical Italian electricity prices.

Reservoir
Run of River +4.4% −0.4% Table 7 compares the income variation in Sabbione and of the Ceppo Morelli modelled in Gaudard et al. [26], both with RoR configuration and the five electricity price scenarios. Our RoR configuration results are aligned with those obtained by [26]. The increase and decrease of income are comparable. In contrast, the configurations of Sabbione outline substantial differences. The results are aligned with DS prices, but not with the three remaining price scenarios. RoR configuration is highly sensitive to any changes in the correlation (or non-correlation) between the availability of resources (discharge) and the demand (electricity prices) due to the impracticality to manage the water [26]. Even the DS price scenario observes a higher income increasing with a reservoir than RoR configuration. Reservoir permits not only to compensate future temporal shifting in the electricity price series, minimizing the losses, but also to increase revenues in case future climate scenarios are led to this effect.

Conclusions
After a review of the existing literature on the future of Italian hydropower, this study assessed a case study in a multi-driver and multi-model framework, both with presence and absence of storage capacity. The hydrological model utilized nine climate scenarios implementing increasing trends in temperature and modification in the pattern of precipitation, to depict future evolutions of glacier and The rate of glacier depletion accelerates up to 2030, indicating the importance of glacier retreat in water-energy nexus framework of hydropower. At 2030, climate changes are expected to change seasonality of inflows, with an anticipation of spring peaks due to earlier beginning of snow and glacier ablation. Annual water volumes reduce despite glacier depletion, with values ranging from −9.2% to −4.5% for climate scenarios implementing temperature trends. However, an adequate management can mitigate these losses by optimizing the hydraulic head and the scheduling of turbines. For these scenarios the reduction of revenues is of −7.6% and −3.4%. If changes in price seasonality will modify the operation of the reservoir hydropower plant, only the amplitudes are expected to affect revenues. A temporal shift is ineffective, as the peak prices remain stable.
Run-of-river power plants are highly sensitive to changes in discharge and prices. The absence of a storage capacity cannot mitigate evolution. The comparison between reservoir and RoR configuration outlined the higher robustness of the former to face changes in discharge, as the annual losses of income are 40% lower. These confirm the importance of a storage capacity in compensating the losses and increasing the revenues in most favorable scenarios.
The results of this study are site specific, but some conclusions are general since the results are sufficiently detailed for decision-making and allow us to formulate recommendations for future hydropower investments, especially by connecting these results with the electricity market. The case study represents a good benchmark for researchers and dam managers, since the climate and price scenarios implemented are simple, and suitable for both short-term and long-term climate projections, especially for glacier and/or snow fed catchments. In particular, the liquid only precipitation scenario is extreme but suitable for far-future analysis.
There is still room for further improvements in the water-energy nexus modelling of hydropower. One should consider glacier evolution according to albedo and solar radiation, possible alternative electricity markets such as balancing and ancillary services, combined analysis of hydropower with future climate change perspective of other renewable sources, among others. Hydroelectricity represents an interesting asset in deploying a sustainable energy mix. However, existing installations should be designed in order to limit social and environmental impacts. Therefore, there is greater need to study with a multidisciplinary approach for recognizing the various factors influencing hydropower in the long run.