Projection of Climate Change and Consumptive Demands Projections Impacts on Hydropower Generation in the São Francisco River Basin, Brazil

Climate change impacts may influence hydropower generation, especially with the intensification of extreme events and growing demand. In this study, we analyzed future hydroelectric generation using a set of scenarios considering both climate change and consumptive demands in the São Francisco River Basin. This project will increase consumptive demands for the coming decades. Five models from the recently released Coupled Model Intercomparison Project Phase 6 and two scenarios, SSP2-4.5 and SSP5-8.5, were considered to estimate climate change projections. The affluent natural flows, regulated flows, and the hydroelectric energy generated were estimated for four multi-purpose reservoirs considering all existing and new demands. The conjunction of scenarios indicated a possible significant reduction in water availability, increased consumptive demands, especially for irrigation, and reduced power generation. Only at the Sobradinho hydroelectric plant, the decrease ranged from −30% to −50% for the period 2021 to 2050 compared to the historical period (1901 to 2000). The results can provide insights into future energy generation and water resources management in the basin.


Introduction
Extreme events such as severe droughts have increased globally in recent years [1]. Such events significantly influence several sectors such as energy production, agriculture, the economy, and society. Hydropower generation is the most widely used renewable energy resource in the world [2]. Its production is highly connected with the consumptive demands for food production and municipal water supply. With growing demands and modification of the current climate, increasing pressure may arise in hydropower generation and cause important changes in water management, particularly during the occurrence of prolonged droughts.
The Northeast Brazil Region (NBR) suffers from recurrent prolonged droughts, directly affecting this region's water, food, and energy security. Since the 18th century, eight drought events persisted for more than three years [3]. However, the recent multiannual drought of 2012-2018 (seven years) was the most prolonged and severe one [3,4]. Due to this historical coexistence with droughts, both National and State Governments have taken measures to mitigate the impact of these events with the construction of multiyear reservoirs, and lately, cisterns, canals, and social programs to assist the most vulnerable populations [5].
The São Francisco River hydroelectric system produces most of the electric power consumed by the NBR. The basin can generate 16,551 GWh per year, 12% of the nation's affluent natural flows ranging from 60% to 90% for the 21st century compared to the 20th century for various river basins in the NBR [18,19].
The sixth phase of the CMIP project (CMIP6) has recently started to release the first results of the models. It brings advances in CGMs modeling and future scenario considerations [20]. The combination of socio-economic and technological development with future radiative forcing scenarios based on updated data on emissions trends was the innovation of CMIP6 climate projections compared to CMIP5 [20,21].
Although several climate change impacts on the São Francisco River Basin have been reported, the impacts on hydroelectricity generation due to future climate associated with growing consumptive uses are rarely addressed, especially using the most recent CMIP6 models. The increase in water demand for the coming decades, caused by population and wealth growth, enhanced agricultural demand due to higher temperatures and the addition of new demands by a water transfer project may exert significant pressure on the hydroelectric power generation of the NBR [16,18].
To confirm this assumption, we proposed a set of scenarios to consider both climate changes and consumptive demand growth. We evaluated how each demand can increase, and its impact on hydroelectricity production using future affluent natural flows.
Thus, this paper aimed at clarifying the role of climate change and growing demands on hydroelectricity production by estimating the impacts that the future climate projected by the recent CMIP6 models and combined with the growth of projected consumptive demands may cause on the affluent natural flows and hydroelectric power generation on the São Francisco River basin. Results should help decision-makers to better prepare for future scenarios and to reduce climate change impacts on the water-energy-food nexus.

Materials and Methods
The responses of the affluent natural flows and hydroelectric generation were analyzed based on the future climate scenarios SSP2-4.5 and SSP5-8.5 of the CMIP6 models and the possible future scenarios of consumptive demands. The affluent natural flows and future demands for each basin of the reservoir, followed by hydroelectric production for these scenarios were estimated, a process divided into six stages, according to Figure 1. Besides the possible decrease in the average annual precipitation in the NBR, these works also showed an increase in temperature of approximately 4° to 5 °C [14][15][16][17]. Higher temperatures and significantly reduced rainfall may result in reductions in the rates of affluent natural flows ranging from 60% to 90% for the 21st century compared to the 20th century for various river basins in the NBR [18,19].
The sixth phase of the CMIP project (CMIP6) has recently started to release the first results of the models. It brings advances in CGMs modeling and future scenario considerations [20]. The combination of socio-economic and technological development with future radiative forcing scenarios based on updated data on emissions trends was the innovation of CMIP6 climate projections compared to CMIP5 [20,21].
Although several climate change impacts on the São Francisco River Basin have been reported, the impacts on hydroelectricity generation due to future climate associated with growing consumptive uses are rarely addressed, especially using the most recent CMIP6 models. The increase in water demand for the coming decades, caused by population and wealth growth, enhanced agricultural demand due to higher temperatures and the addition of new demands by a water transfer project may exert significant pressure on the hydroelectric power generation of the NBR [16,18].
To confirm this assumption, we proposed a set of scenarios to consider both climate changes and consumptive demand growth. We evaluated how each demand can increase, and its impact on hydroelectricity production using future affluent natural flows.
Thus, this paper aimed at clarifying the role of climate change and growing demands on hydroelectricity production by estimating the impacts that the future climate projected by the recent CMIP6 models and combined with the growth of projected consumptive demands may cause on the affluent natural flows and hydroelectric power generation on the São Francisco River basin. Results should help decision-makers to better prepare for future scenarios and to reduce climate change impacts on the water-energy-food nexus.

Materials and Methods
The responses of the affluent natural flows and hydroelectric generation were analyzed based on the future climate scenarios SSP2-4.5 and SSP5-8.5 of the CMIP6 models and the possible future scenarios of consumptive demands. The affluent natural flows and future demands for each basin of the reservoir, followed by hydroelectric production for these scenarios were estimated, a process divided into six stages, according to Figure 1.  In Step 1, the annual time series of precipitation variables and maximum, minimum and mean upper air temperatures (UAT) observed, as well as those projected by CMIP6 were extracted with the classic Thiessen method [22] for the basin of the reservoirs of São Francisco River Basin (SFRB). Still in Step 1, potential evapotranspiration was estimated with the maximum, minimum, and mean UAT using the Hargreaves-Sammani method.
In Step 2, the statistical corrections to remove bias from the precipitation and Potential Evapotranspiration (PET) time series of the CMIP6 models were performed with the observed data. For this purpose, the gamma cumulative distribution function (CDF) was used.
In Step 3, the precipitation and PET data with bias removed from the CMIP6 were used as input by the hydrological rainfall-runoff model soil moisture accounting procedure (SMAP) to generate the natural flow data from the basins of Retiro Baixo, Sobradinho, Itaparica, and Três Marias. Still in Step 3, the affluent natural flows were estimated for the remaining reservoirs (Xingó, Moxotó, Complexo Paulo Afonso and Queimado). For this purpose, we used the natural flow data generated with the SMAP model and the monthly series of naturalized flows made available by the Brazilian National Electrical System Operator (ONS in Portuguese), dividing them in two groups with four reservoirs each. The four reservoirs with the calibrated SMAP model were used to predict the affluent natural flows for the other four remaining reservoirs.
In Step 4, the municipal demands were quantified with geoprocessing techniques and aggregated for each main reservoir. Then, their projections for the period 2017 to 2050 were estimated using the exponential smoothing model. In Step 5, the operation of the reservoir was simulated, focusing on meeting water demands through the Information System for Water Allocation Management (SIGA in Portuguese) software [23]. For this purpose, the natural flow data obtained with the SMAP model, data from the CMIP6 models in the SSP2-4.5 and SSP5-8.5 scenarios, and the consumptive demands scenarios were used. Still in Step 5, with the SIGA simulation results, it was possible to estimate the hydroelectric energy generation.
In Step 6, the regularized flows with 90% duration (Q90) and hydroelectric power analysis were performed through the Mann-Kendall Sen test and percentage anomaly. The next topics discuss the study region, database, and methodology applied.

Study Area
The study region of this work was the SFRB, as shown in Figure 2. It drains an area of 639,219 km 2 , which is equivalent to about 8% of the Brazilian territory. Starting in Serra da Canastra in Minas Gerais, the São Francisco River is 2700 km long. It flows through the states of Bahia and Pernambuco in a south-north direction, reaching the Atlantic Ocean through the border of Alagoas and Sergipe. SFRB traverses seven states: Bahia, Minas Gerais, Pernambuco, Alagoas, Sergipe, Goiás, and Distrito Federal. Besides, it covers 507 municipalities (almost 9% of all municipalities in Brazil), with a resident population that reaches 20 million people [24].
The main SFRB reservoirs are Três Marias, Sobradinho, and Itaparica (also known as Luíz Gonzaga), which have useful volumes of 15,278 hm 3 , 28,669 hm 3 and 3549 hm 3 , and surface areas of 1040 km 2 , 4214 km 2 and 828 km 2 , respectively [25]. The water transfer project will divert water in the São Francisco River between Sobradinho and Itaparica dams, in the State of Pernambuco (see Figure 2). The project comprises two independent canals, water pumping stations, small reservoirs, and a hydroelectric power plant (HPP). It will meet the water supply needs for multiple uses in the states of Pernambuco, Paraíba, Rio Grande do Norte, and Ceará [24].
The project region is known for long periods of drought. The area that suffers the most from the effects of these prolonged droughts is located in the semi-arid part north of the São Francisco River. This area partially covers the states of Pernambuco, Ceará, Paraíba, and Rio Grande do Norte (regions that will be comprised with the transposition). The enterprise will enable water supply for various purposes (human, irrigation, and animal supply, fish, and shrimp farming) in an area with about 12 million inhabitants [24]. Before completing the transposition axes, the current total demand in the basin accounts for 77% of irrigation, 11% of urban demand, and 7% of industrial demand [24]. The main SFRB reservoirs are Três Marias, Sobradinho, and Itaparica (also known as Luíz Gonzaga), which have useful volumes of 15,278 hm 3 , 28,669 hm 3 and 3549 hm 3 , and surface areas of 1040 km 2 , 4214 km 2 and 828 km 2 , respectively [25]. The water transfer project will divert water in the São Francisco River between Sobradinho and Itaparica dams, in the State of Pernambuco (see Figure 2). The project comprises two independent canals, water pumping stations, small reservoirs, and a hydroelectric power plant (HPP). It will meet the water supply needs for multiple uses in the states of Pernambuco, Paraíba, Rio Grande do Norte, and Ceará [24].
The project region is known for long periods of drought. The area that suffers the most from the effects of these prolonged droughts is located in the semi-arid part north of the São Francisco River. This area partially covers the states of Pernambuco, Ceará, Paraíba, and Rio Grande do Norte (regions that will be comprised with the transposition). The enterprise will enable water supply for various purposes (human, irrigation, and animal supply, fish, and shrimp farming) in an area with about 12 million inhabitants [24]. Before completing the transposition axes, the current total demand in the basin accounts for 77% of irrigation, 11% of urban demand, and 7% of industrial demand [24].
Besides meeting the consumptive demands, SFRB plays an essential role in electricity generation, with installed potential corresponding to 12% of the nation's total [25]. However, due to the drought in the region that began in 2012, hydroelectric generation supplied only 56%, 42%, 34%, 30%, 21%, and 21% of the total electrical energy demand by the NBR from 2013 to 2018. This production was far below previous years. The decline in hydroelectricity production is illustrated in Figure 3. This shortfall was mainly replaced by fossil fuel energy (thermal) with 39%, 47%, 42%, 33% and 25%, and the growth of wind energy with an increase from 5% to 52% during this period. Besides meeting the consumptive demands, SFRB plays an essential role in electricity generation, with installed potential corresponding to 12% of the nation's total [25]. However, due to the drought in the region that began in 2012, hydroelectric generation supplied only 56%, 42%, 34%, 30%, 21%, and 21% of the total electrical energy demand by the NBR from 2013 to 2018. This production was far below previous years. The decline in hydroelectricity production is illustrated in Figure 3. This shortfall was mainly replaced by fossil fuel energy (thermal) with 39%, 47%, 42%, 33% and 25%, and the growth of wind energy with an increase from 5% to 52% during this period.

Observed Data
For the calibration of the parameters of the SMAP rainfall-runoff model, monthly historical precipitation data obtained from the Brazilian National Water Agency (in Portuguese, ANA) rainfall stations from 1961 to 2016 [26] and maximum, minimum, and mean temperatures from the Climate Research Unit (CRU) with spatial resolution of 0.5°

Observed Data
For the calibration of the parameters of the SMAP rainfall-runoff model, monthly historical precipitation data obtained from the Brazilian National Water Agency (in Portuguese, ANA) rainfall stations from 1961 to 2016 [26] and maximum, minimum, and mean temperatures from the Climate Research Unit (CRU) with spatial resolution of 0.5 • × 0.5 • were used [27]. General information on the Level-Area-Volume curve of the reservoirs for the SIGA program usage was obtained from the operating reservoir companies, the Minas Gerais Energy Company (CEMIG) regarding the Três Marias reservoir, and the São Francisco Hydroelectric Company (CHESF) for the other reservoirs [25,28]. Finally, the monthly hydroelectric energy observed from 1999 to 2016 from the HPP of the São Francisco River was obtained from the ONS and used to validate the hydroelectric energy simulated by the SIGA software [9].   [33] In this study, the precipitation variables, monthly mean, maximum and minimum UAT for the historical scenario (1901 to 2000, 20th century) and SSP2-4.5 and SSP5-8.5 projection scenarios for the 21st century were used (2016-2100).
There are five different shared socioeconomic pathways (SSPs). The future sustainable potential and fossil fuel-driven growth are the most extreme scenarios, SSP1 and SSP5, respectively [34]; high inequality between or within countries, are given by SSP3 and SSP4, respectively [35] and the SSP2 scenario is the "middle term" of the aforementioned scenarios [36]. For each SSP, different Representative Concentration Pathways (RCPs) represent different radiative forcing projections leading to an extensive range of global warming levels throughout the century [37]. This new structure allows the standardization of all socio-economic assumptions (e.g., population, gross domestic product and poverty, among others) in each scenario's modeled representations. Plus, it allows a more subtle investigation of the variety of paths that climate outcomes can follow.

PET Estimation through Hargreaves-Samani Method
For PET estimation, the mean, maximum and minimum UAT data (in Celsius degrees) of the CMIP6 models were used in the SSP2-4.5 and SSP5-8.5 projection scenarios using the Hargreaves-Samani method [38,39], according to Equation (1): where PET is given in mm/month; and the Average External Radiation (Ra) was estimated from the latitude and month of the year, according to the work of [38].

Bias Correction
The results of climate models, such as the CMIP6 models, show systematic errors, especially related to bias [40]. Thus, the statistical correction via CDF with the mapping probability was performed in a monthly time series of precipitation and PET to reduce such errors. CDF is widely used in several studies that use climate models, for bias correction of precipitation and PET variables [17][18][19]41].
The methodology is based on the study by [42] in which the accumulated probability distribution curves of the monthly modeled and observed series are generated, obtaining 12 gamma adjustments, one for each month of the year. After that, it was possible to check the probabilistic behavior of the modeled data in relation to those observed. Thus, the corrected precipitation/PET value of the CMIP6 models was obtained by taking the precipitation/PET value corresponding to the same probability of occurrence in the observed values curve.

The Soil Moisture Accounting Procedure (SMAP) Rainfall-Runoff Hydrological Model
In this study, the monthly SMAP model was used. Developed by [43], it uses little observational data and low computational demands, making it very popular. SMAP covers the process of hydrological balance in two hypothetical reservoirs, one of soil (R soil ) and one underground (R sub ).
The water stored in each reservoir is updated at each time interval due to the incorporation of a new average rainfall data through the following transfer Equations (2)-(5): where E s is the surface flow (mm), P is the precipitation (mm), T u is the initial moisture content (adimensional), E r is the actual evapotranspiration, E p is the potential evapotranspiration (mm), P es is the surface flow parameter (adimensional), C rec is the recharge coefficient (adimensional), E b is the base flow (mm), K the recession constant (month −1 ), R soil is the soil reservoir, R sub is the underground reservoir and R ec is the underground recharge. Figure 4 shows how the monthly SMAP model works.  T uin is the initial soil moisture rate that determines the initial level of the second reservoir R soil and the initial base flow (E bin ) that defines the initial value of R sub . In this work, the calibration took place in a single step, in which the parameters Sat, P es and C rec and K of level R sub that generates the E b underwent an optimization through the solver tool present in the Microsoft Excel spreadsheet, followed by a manual refinement by trial T uin is the initial soil moisture rate that determines the initial level of the second reservoir R soil and the initial base flow (E bin ) that defines the initial value of R sub . In this work, the calibration took place in a single step, in which the parameters Sat, P es and C rec and K of level R sub that generates the E b underwent an optimization through the solver tool present in the Microsoft Excel spreadsheet, followed by a manual refinement by trial and error. Additionally, T uin and E bin were manually adjusted with the Nash-Sutcliffe efficiency coefficient (NASH) for each simulation.
According to [44], NASH is one of the most used criteria in Hydrology and is an objective function that measures the adjustment efficiency between a model and observed data. According to the formulation presented by [45], NASH is calculated through Equation (6): where N the samples number, o i the observed data, m i the data being modeled and o i the average of the observed data.
The reference values for NASH coefficient, suggested by [44,46], classify as acceptable the values from 0.36 to 0.75 and excellent above 0.8. Its maximum value corresponding to 1, which indicates a perfect adjustment between the modeled and observed data. Table  2 shows the parameters of the SMAP model calibrated for the SFRB reservoirs, using the estimated natural flow data from the ONS. The different periods in the calibration of the reservoirs took place according to the availability of data from ANA stations. From those, only the Itaparica reservoir was not classified as excellent calibration, presenting an acceptable NASH coefficient.

Linear Regression Model
The linear regression model was used in this study to generate the flow rates of the HPP that do not have the calibrated SMAP model (Queimado, Xingó, Moxotó, Complexo Paulo Afonso 1, 2, 3 and 4). This model was used with good performance at representing the flows in several reservoirs [17,19]. The first methodological stage occurred with the standardization of the monthly series of naturalized flows using Equation (7): where Z is the normalized flow, i the number of months (ranging from 1 to 12), k the number of years (ranging from 1931 to 2005), j the number of reservoirs (a total of 8), q i,j,k naturalized flow from station j in month i in year k, q i,j is the matrix that represents the average of all the months and ranks, and σ i,j is the matrix that represents the standard deviation of the monthly series of all ranks. The second step was to perform the linear regression of ONS naturalized flows to obtain the parameters of each reservoir that does not have the SMAP calibrated, considering all other reservoirs as explanatory variables (Itaparica, Sobradinho, Três Marias and Retiro Baixo). The linear regression was given by Equation (8): where pk are the flows of the four basins obtained through the SMAP model (varying from 1 to 4), pj are the reservoirs that do not have the SMAP calibrated (varying from 1 to 4) and β i,pk the regression coefficients. Figure 5 shows the linear regression model performance in estimating the flows of the four reservoirs that do not have the SMAP model calibrated. It is possible to observe that the Queimado and Xingó reservoirs presented acceptable results. In contrast, in the Moxotó and Complexo Paulo Afonso reservoirs 1, 2, 3 and 4, results were a little below the acceptable value of 0.36, showing 0.34 and 0.35, respectively. Thus, with the coefficients obtained from the ONS naturalized flow series and the flows obtained through the SMAP model with the data from the CMIP6 models, the affluent natural flows for the climate change scenarios for the other reservoirs were obtained using the linear regression parameters.
number of years (ranging from 1931 to 2005), j the number of reservoirs (a total of 8), q i,j,k naturalized flow from station j in month i in year k, q i,j is the matrix that represents the average of all the months and ranks, and σ i,j is the matrix that represents the standard deviation of the monthly series of all ranks.
The second step was to perform the linear regression of ONS naturalized flows to obtain the parameters of each reservoir that does not have the SMAP calibrated, considering all other reservoirs as explanatory variables (Itaparica, Sobradinho, Três Marias and Retiro Baixo). The linear regression was given by Equation (8): where pk are the flows of the four basins obtained through the SMAP model (varying from 1 to 4), pj are the reservoirs that do not have the SMAP calibrated (varying from 1 to 4) and β i,pk the regression coefficients. Figure 5 shows the linear regression model performance in estimating the flows of the four reservoirs that do not have the SMAP model calibrated. It is possible to observe that the Queimado and Xingó reservoirs presented acceptable results. In contrast, in the Moxotó and Complexo Paulo Afonso reservoirs 1, 2, 3 and 4, results were a little below the acceptable value of 0.36, showing 0.34 and 0.35, respectively. Thus, with the coefficients obtained from the ONS naturalized flow series and the flows obtained through the SMAP model with the data from the CMIP6 models, the affluent natural flows for the climate change scenarios for the other reservoirs were obtained using the linear regression parameters.

Consumptive Demands
The withdrawal flows of consumptive demands (irrigation, human supply and industry) from 1961 to 2017 were obtained from ANA for the state municipalities in Brazil and were aggregated in the basin of the reservoirs, see Figure 2. For this purpose, they were quantified through geoprocessing techniques, being expressed through Equation (9): where QD b is the average flow of demand located inside the basin b (m 3 /s); QD m is the average flow of demand inside the municipality m (m 3 /s); AB mb is the area of the municipality m inside the basin b (km 2 ); and AT m is the total area of the municipality m (km 2 ). The ANA methodology to estimate the consumptive demands was based on investigations of the technical coefficients available in the literature and their compatibility with the data already available [47]. More details on the estimates of the demands are available in the Manual of Consumptive Uses of Water in Brazil in [47], as well as the metadata used in this study.
Four scenarios were considered to analyze the influence of consumptive demands in the future (2018 to 2050) for the river basins that make up the SFRB: Scenario 1 with static demands, i.e., considering the last observed withdrawal flow; Scenario 2 with consumptive demands projected with the exponential smoothing (ETS) model; Scenario 3 with consumptive demands projected with the ETS model considering values below 95% of the confidence interval; and Scenario 4 with consumptive demands projected with the ETS model considering values above 95% of the confidence interval.
The ETS model was proposed in the late 1950s [48][49][50], and the predictions generated by this model are produced through weighted averages of past observations, with the weights decreasing exponentially as the age of the observations, i.e., the more recent the observation, the greater the associated weight. Thus, being y T an observation in time T of a time series, the forecast in time T + 1 will be given by: where α is the parameter for smoothing the growth rate that decreases with time. However, the process has to start at T = 1 with some adjusted value. This adjusted value is called level and is represented by l 0 in Equation (11): Moreover, the ETS models are defined in terms of components, making it easier to add extensions, such as trend and seasonality [51]. The name ETS is due to its three parameter extensions: the Error, the Trend and the Seasonal components (E, T, S). These components, in turn, can be: in the Error, Additive (A) or Multiplicative (M) component; in the Seasonal, A or M or none (N) components; and in the Trend, A or M or N component or Dampened Additive (Ad) or Dampened Multiplicative (Md) [51].
Since its creation, with the variation of these components, a total of 30 models have been developed [51][52][53][54]. Each model consists of an equation that will describe the data for l 0 , trend (b) and seasonality (s) as a function of time t (l t , b t , s t , s t−1 , . . . , s t−m + 1 ) using likelihood and space methods for the calculations. In general, these equations are expressed as Equations (12) and (13): where y t is the forecast as a function of a state vector x t , which will have a transition that describes the states evolution and ε t is the error associated with white noise. For an ETS model (A, M, M), that is, with additive error, with tendency and multiplicative seasonal component, the expressions for y t , l 0 , b t and s t (Equations (14)-(17)), are: where β and γ are, respectively, the growth parameter of b t and the component parameter of s t . For further details on the ETS models equations, please refer to [53]. For this study, ETS was based on the methods described by [53], available through the "forecast" package in the R environment. The best ETS model was chosen with Akaike's information criterion (AIC) (Equation (18)), Schwarz's Bayesian information criterion (BIC) (Equation (19)) and AIC with removed bias (AICc) (Equation (20)).
where LL is the log likelihood, T p is the total of parameters and T is the number of observations.
where k is the estimation of the parameters of the models obtained through the least squares method.
where the expression is the bias correction.

The Decision Support System (SIGA)
The reservoir operation was done through a flow network decision support system named SIGA [23]. It enables the input of different demands and priorities. Furthermore, it provides the supplied durations inherent to the permanence curve of the monthly affluent flows of the reservoirs. The four multiple-use reservoirs of the São Francisco River basin were modeled as a cascade series, as shown in Figure 6. The new demands of the two channels were considered to be met by the system by unifying its uses: human, industrial and irrigation. These uses were also modeled on the downstream stretch of each basin reservoir, through demands. Besides meeting the demands, the modeled system must also need to meet the maximum and minimum discharge limits in the river sections as regulated by the system operator [6], according to Table 3. The values of the minimum limits were used in this study. Besides meeting the demands, the modeled system must also need to meet the maximum and minimum discharge limits in the river sections as regulated by the system operator [6], according to Table 3. The values of the minimum limits were used in this study. The simulation of the flow network is performed by means of the equation of the water mass balance in reservoirs expressed in Equation (21): where V t is the volume stored at the beginning of time period t (hm 3 ); V t+1 is the volume of the reservoir at the beginning of the next time period t + 1 (hm 3 ); I t is the volume affluent to the reservoir during time period t (natural and transfer) (hm 3 ); E t * is the evaporated slide during the time period t, assumed constant along the period (hm); A t is the water surface area at the beginning of period t, assumed constant for short time intervals and function of V t (hm 2 ); R t are the operational withdrawals (hm 3 ); and S t is the overflow (withdrawal of uncontrolled water) (hm 3 ). The flow network simulations had the operational withdrawals optimized to guarantee the fulfillment of the operational limits, maximize the fulfillment of demands and minimize the evaporated volume, during the chosen simulation periods, through the SIGA software [23]. In order to do so, SIGA uses the multiple objective particle swarm optimization (MOPSO) as an optimization algorithm [23].
During a period of failure to meet the demands, SIGA adopts a system of priorities, from the lowest to the highest, to determine the system nodes that will suffer shortages. This study considered the priority system defined in Table 4, which aims to prioritize human supply over irrigation and industry. The second highest priority was attributed to the new water transfer project, supposedly used to assist other human supply systems of basins.

Hydroelectric Energy Production
Conventionally, the power produced in a HPP is given by: where G is a constant with a value of 9.81.10 −3 (kg/m 2 s 2 ) that represents the product between the density of the water (1000 Kg/m 3 ), the gravity (g) and constant with a value of 10 6 that converts the energy of (J) into (MJ); η t is the hydraulic performance of the turbine; hl is the penstock head loss (m); q is the turbine flow in only one generating unit (m 3 /s).
Usually, the net drop height has been represented using fourth-order polynomials depending on the stored volume (in the upstream level) and the discharged flow (in the downstream level). The following polynomial gives the value of the upstream level: where fcm is the value of the upstream level (m); a 0 , . . . ,a 4 are the coefficients of the polynomial that represents the upstream level for the reservoir; V is the stored volume (hm 3 ).
On the other hand, the downstream level of the plant is the river level after the turbine. The same can be given through the polynomial that relates the discharged flow (turbine flow-Q plus spilling flow-S into the HPP) with the respective coefficients: where fcj is the value of the downstream elevation (m); S is the spilled flow rate from the HPP (m 3 /s); b0 . . . b4 are the coefficients of the polynomial that represents the downstream elevation to the reservoir. Thus, based on Equations (23) and (24), the forebay elevation (m) becomes: However, not all of this forebay elevation is available at the turbine entrance, as part of it is lost when the water passes through the forced conduit and other through penstock head loss. The net fall height of the HPP is defined as the difference between hb and the hydraulic loss, expressed by: where pl is the hydraulic net loss (m). The hydraulic losses, in turn, were modeled as a quadratic function of the turbine flow and a certain constant, as follows: where k p is a constant that depends on the physical characteristics of the forced duct that connects the reservoir to a certain HPP (s 2 /m 5 ).
As for the turbine performance, to build a model closer to the reality of the hydraulic performance, two variables must be taken into consideration: the net drop height it is submitted to and its turbine flow [55]. Thus, this data set is represented through the following concave quadratic function: where r 0 , . . . , r 5 are the coefficients of the polynomial that represent the yield of a given hydroelectric unit. The values of the coefficients of the polynomials that are upstream (variable "a", Equation (23)), downstream (variable "b", Equation (24)) and the efficiency of the turbine (variable "r", Equation (28)) for each reservoir, as well as the constants used to obtain the hydraulic loss (variable "kp", Equation (27)), are in Table 5.
These levels were represented through fourth-degree polynomials, obtained through polynomial regression. The turbine-generator yield (%) of HPP Itaparica, provided by [55], equivalent to η equal to 91%, was used. The kp constant to obtain the hydraulic losses, belonging to Equation (28), considered the kp values of 2.53 × 10 −5 , for HPP Três Marias; 9.84 × 10 −7 , for HPP Sobradinho; and 0 for HPP Itaparica, which made pl equal to zero, that is, without losses.   These levels were represented through fourth-degree polynomials, obtained through polynomial regression. The turbine-generator yield (%) of HPP Itaparica, provided by [55], equivalent to η equal to 91%, was used. The kp constant to obtain the hydraulic losses, belonging to Equation (28), considered the kp values of 2.53 × 10 −5 , for HPP Três Marias; 9.84 × 10 −7 , for HPP Sobradinho; and 0 for HPP Itaparica, which made pl equal to zero, that is, without losses. Figure 7 shows the validation with the NASH coefficient obtained comparing the estimation of the energy generated by Itaparica, Sobradinho and Três Marias HPP obtained through the SIGA software with the observed series obtained from ONS in the period 1999 to 2016. As NASH prioritizes the maximum values of the time series [56], even the hydroelectric power generated with the model presenting discrepancies in Três Marias when As NASH prioritizes the maximum values of the time series [56], even the hydroelectric power generated with the model presenting discrepancies in Três Marias when compared to the observed hydroelectric power, the maximum values were well represented, resulting in NASH of approximately 0.64.

Projection Analysis
The analysis of the projections of the affluent natural flows and hydroelectric energy were carried out for the four reservoirs and the three main HPP of the São Francisco River, through the Mann-Kendall Sen test, Q90 and percentage anomaly.

The Mann-Kendall Sen Test
The non-parametric Mann-Kendall Sen test was used to check trends in the flow series. According to [57], the nonparametric Mann-Kendall-Sen test has been widely used for hydrological studies and is recommended by the WMO.
Therefore, according to [58], the statistics of the Mann-Kendall test for a series (Z1, Z2, . . . , Zn) coming from a sample of n independent random variables and identically distributed, is given by Equation (29): where S are the values of the series in annual time intervals; i and j are the time indexes, and n is the number of series elements [57]. The term sign (Z j − Z i ) is the result of Equation (30): For the null hypothesis test H0 (which considers not having a positive trend), the TAU variables and p-value were used. The TAU variable is related to the correlation classification coefficient and quantifies the monotonic association, being given by Equation (31): where n is the size of the series. H0 is accepted in the Mann-Kendall TAU when the test is less than a critical value named α, which for this study was α = 0.05 (for a statistical significance of 95%), that is, for TAU < α, the series has no positive trend. Otherwise, i.e., for TAU ≥ α, the time series has a positive tendency. Otherwise, the p-value of the statistic S considers the H0 true, for p-value > α and false for p-value ≤ α [57].
The Sen estimator provides the magnitude of the detected trends. According to [57], it is estimated through the statistic Q, given by Equation (32): where X i and X j are related to the values of the variable under study in times i and j [57]. The positive or negative value for Q indicates increasing or decreasing tendency, respectively.

Permanence Curve and Regularized Flow (Q90)
The performance of the system was measured in terms of meeting the consumptive demands during the observed period and in two future scenarios: (i) the possible increase of gas emissions that intensify the study effect (SSP2-4.5 and SSP5-8.5); and (ii) possible future annual consumptive demands (D1, D2, D3 and D4). The indicator used to measure this degree of service was the permanence curve (PC).
The permanence curve or duration of flows is equivalent to a hydrogram of the accumulated frequency of flows in a river [59]. The same is expressed through the relationship between the flow and the frequency (empirical estimate) with which this flow is exceeded or equalized, such as: where m is the order of the ordered value and n is the number of the series values. This study used monthly flows. PC helps in the flow data analysis when the percentage of time the river has flows in a given range is required or what percentage of time a river has sufficient flow to meet a given demand [59]. In this study, PC was used to obtain the flow that is exceeded 90% of the time, which is called regularized flow or Q90. This has been used as a reference for legislation in environment and water resources in most Brazilian states and provides a flow duration that supplies 90% of the time to end-users.

Anomaly Percentage
For the statistical analysis of the annual energy projections, the calculation of anomaly percentage was used. It is expressed in Equation (27) and used the CMIP6 model projections related to scenarios SSP2-4.5 and SSP5-8.5 and the demands D1, D1, D3 and D4 in the period from 2021 to 2050 (21st century). A comparison is made with the representation of the 20th century (1901 to 2000) simulated by CMIP6 models using historical simulated data.
where E a XXI is the annual average energy of the projections for the scenarios of the 21st century and E a XX is the annual average energy of the of the 20th century.

Streamflow
The observed historical time series of affluent natural flows for the period 1931 to 2019 (see Figure 8) show high inter-annual variability, with maximum values close to 5000 m 3 /s in the reservoirs of Itaparica and Sobradinho, above 1500 m 3 /s in Três Marias and close to 350 m 3 /s in Retiro Baixo. The minimum values were below 1000 m 3 /s in the reservoirs of Itaparica and Sobradinho, 500 m 3 /s in Três Marias and 100 m 3 /s in Retiro Baixo, and they were all registered in the last years. This fact reflects the prolonged water scarcity that has reached the SFRB in recent years, especially in those basins that are totally or partially inserted in the semi-arid region, such as the Itaparica and Sobradinho. In those basins, the Mann-Kendall Sen test indicated a negative trend. to [15,61], the last positive phase of AMO occurred from 1995 onwards and still continues. Therefore, the low values of affluent natural flows observed may be partly due to AMO's positive phase and the expected higher Sea Surface Temperatures (SST) in the tropical Atlantic Ocean. However, even with the decadal cycle indicating a phase with lower affluent natural flows since 1995, this positive phase of AMO is unprecedented and may be attributed to the extreme increase in global SST associated with anthropogenic global warming [15,61].  For the set of models that showed a significant trend, the increase in greenhouse gas emissions suggests a more significant impact on flows. In most cases, the slope module is larger for the SSP5-8.5 scenario than for the SSP2-4.5 scenario. This fact may be related to the considerable temperature increase projected by the SSP5-8.5 scenario, resulting in increased PET projected for the 21st century. This result was also found with the CMIP5 data by [17]. The MIROC6 and MRI-ESM2-0 models indicated a positive trend. In contrast, the other models projected a negative trend, with the declivity of the CanESM5 model presenting a value below −30 for the Queimado, Sobradinho, Itaparica, Moxotó, Xingó and Complexo Paulo Afonso reservoirs. In addition, a decadal variability with well-defined periods was observed, with higher values of affluent natural flows between 1930 and 1960 and between 1980 and 1995 and lower flows between 1960 and 1990 in all reservoirs, which coincided with the phases of the Atlantic Multidecadal Oscillation (AMO) and the Pacific Decadal Oscillation (PDO) [15,32,60,61]. According to [60], the period from 1930 to 1960 adjusted to the positive phase of AMO (with higher temperatures in the tropical Atlantic Ocean) with the drought period in the NBR, while the lower temperatures of the tropical Atlantic Ocean, i.e., the negative phase of AMO, coincided with greater rainfall in the region from 1960 to 1990. According to [15,61], the last positive phase of AMO occurred from 1995 onwards and still continues. Therefore, the low values of affluent natural flows observed may be partly due to AMO's positive phase and the expected higher Sea Surface Temperatures (SST) in the tropical Atlantic Ocean.
However, even with the decadal cycle indicating a phase with lower affluent natural flows since 1995, this positive phase of AMO is unprecedented and may be attributed to the extreme increase in global SST associated with anthropogenic global warming [15,61]. Figure 9 shows the trend declivities of the CMIP6 models for the SSP2-4.5 and SSP5-8.5 scenarios (from 2017 to 2100) concerning the history (period from 1901 to 2000) of the affluent natural flows of the studied reservoirs, according to the Mann-Kendall Sen test. For the set of models that showed a significant trend, the increase in greenhouse gas emissions suggests a more significant impact on flows. In most cases, the slope module is larger for the SSP5-8.5 scenario than for the SSP2-4.5 scenario. This fact may be related to the considerable temperature increase projected by the SSP5-8.5 scenario, resulting in increased PET projected for the 21st century. This result was also found with the CMIP5 data by [17]. The MIROC6 and MRI-ESM2-0 models indicated a positive trend. In contrast, the other models projected a negative trend, with the declivity of the CanESM5 model presenting a value below −30 for the Queimado, Sobradinho, Itaparica, Moxotó, Xingó and Complexo Paulo Afonso reservoirs.  Figure 10 presents the box-plots of the ratio between Q90 for the SSP2-4.5 and SSP5-8.5 scenarios varying in 30 years (2021-2050 and 2051-2080), respectively, and Q90 in the historical period (1901-2000) for the SFRB reservoirs, corresponding to the five models analyzed. The models indicated reductions in Q90 for most of the SFRB reservoirs and the two analyzed periods. The exception was the Itaparica reservoir, which for the SSP5-8.5 scenario from 2021 to 2050, the model median indicated an increase in Q90. The increase in Q90 presented by the Itaparica basin may be related to the low calibration of the SMAP  Figure 10 presents the box-plots of the ratio between Q90 for the SSP2-4.5 and SSP5-8.5 scenarios varying in 30 years (2021-2050 and 2051-2080), respectively, and Q90 in the historical period  for the SFRB reservoirs, corresponding to the five models analyzed. The models indicated reductions in Q90 for most of the SFRB reservoirs and the two analyzed periods. The exception was the Itaparica reservoir, which for the SSP5-8.5 scenario from 2021 to 2050, the model median indicated an increase in Q90. The increase in Q90 presented by the Itaparica basin may be related to the low calibration of the SMAP model (with a low NASH index of 0.36 according to Table 2), since this basin is located below Sobradinho and the water transfer project for the semi-arid region of the NRB, for which the expectation would be a decrease in Q90.

Projection of Consumptive Demands
In the historical time series, annual consumptive demands withdrawn from 1961 to 2017 increased over time, as shown in Figures 11-14 and Table 6. The demand for irrigation had the highest average annual growth rates: 6.80%, 7.42%, 10.99% and 9.29% in Itaparica, Sobradinho, Três Marias and Retiro Baixo, respectively. It is also the irrigation demand that has the highest values in three of the four analyzed reservoirs. The exception was the Retiro Baixo reservoir, where the human demand exceeded the others.
The ETS model projected growth scenarios for the period 2018 to 2050 for most of the consumptive demands. The exception was the D2 scenario (more optimistic). In this scenario, industrial demands reduced its withdrawals with an average annual rate of −3.73%, −0.70%, −0.08% and −1.53% for the reservoirs of Itaparica, Sobradinho, Três Marias and Retiro Baixo, respectively. Human supply also showed a reduction on the average annual rate of −1.02% for the reservoir of Três Marias in the D2 scenario.

Projection of Consumptive Demands
In the historical time series, annual consumptive demands withdrawn from 1961 to 2017 increased over time, as shown in Figures 11-14 and Table 6. The demand for irrigation had the highest average annual growth rates: 6.80%, 7.42%, 10.99% and 9.29% in Itaparica, Sobradinho, Três Marias and Retiro Baixo, respectively. It is also the irrigation demand that has the highest values in three of the four analyzed reservoirs. The exception was the Retiro Baixo reservoir, where the human demand exceeded the others.

Projection of Consumptive Demands
In the historical time series, annual consumptive demands withdrawn from 1961 to 2017 increased over time, as shown in Figures 11-14 and Table 6. The demand for irrigation had the highest average annual growth rates: 6.80%, 7.42%, 10.99% and 9.29% in Itaparica, Sobradinho, Três Marias and Retiro Baixo, respectively. It is also the irrigation demand that has the highest values in three of the four analyzed reservoirs. The exception was the Retiro Baixo reservoir, where the human demand exceeded the others.
The ETS model projected growth scenarios for the period 2018 to 2050 for most of the consumptive demands. The exception was the D2 scenario (more optimistic). In this scenario, industrial demands reduced its withdrawals with an average annual rate of −3.73%, −0.70%, −0.08% and −1.53% for the reservoirs of Itaparica, Sobradinho, Três Marias and Retiro Baixo, respectively. Human supply also showed a reduction on the average annual rate of −1.02% for the reservoir of Três Marias in the D2 scenario.            The ETS model projected growth scenarios for the period 2018 to 2050 for most of the consumptive demands. The exception was the D2 scenario (more optimistic). In this scenario, industrial demands reduced its withdrawals with an average annual rate of −3.73%, −0.70%, −0.08% and −1.53% for the reservoirs of Itaparica, Sobradinho, Três Marias and Retiro Baixo, respectively. Human supply also showed a reduction on the average annual rate of −1.02% for the reservoir of Três Marias in the D2 scenario.

Hydroelectric Energy Production
Future hydroelectric energy production considered CMIP6 climatic scenarios SSP2-4.5 and SSP5-8.5 (in the period from 2021 to 2050), and projected demands of D1, D2, D3 and D4 scenarios relative to the historical period from 1901 to 2000. Figure 15 indicates the anomalies behavior in the average percent of the hydroelectric energy produced at Itaparica, Sobradinho and Três Marias HPP.

Hydroelectric Energy Production
Future hydroelectric energy production considered CMIP6 climatic scenarios SSP2-4.5 and SSP5-8.5 (in the period from 2021 to 2050), and projected demands of D1, D2, D3 and D4 scenarios relative to the historical period from 1901 to 2000. Figure 15 indicates the anomalies behavior in the average percent of the hydroelectric energy produced at Itaparica, Sobradinho and Três Marias HPP. The SSP5-8.5 scenario models indicated more intense reductions of hydropower than the SSP2-4.5 scenario for two, out of the three HPP analyzed in this study. Itaparica HPP was the only exception, with the models indicating, for the two scenarios of future climate and most scenarios of consumptive demands, an increase of hydropower generation.
For the consumptive demands, Sobradinho's HPP was the one that showed the highest sensitivity to future scenarios of these demands. It presented a more significant decrease in the hydroelectricity generation when the demands became more pessimistic. The D4 scenario presented greater intensity of the decrease in hydroelectric power with a median of less than −50%. For Três Marias and Itaparica HPP, there was little change of the consumptive demands, with the medians of the anomalies close to −10% and 10%, respec- The SSP5-8.5 scenario models indicated more intense reductions of hydropower than the SSP2-4.5 scenario for two, out of the three HPP analyzed in this study. Itaparica HPP was the only exception, with the models indicating, for the two scenarios of future climate and most scenarios of consumptive demands, an increase of hydropower generation.
For the consumptive demands, Sobradinho's HPP was the one that showed the highest sensitivity to future scenarios of these demands. It presented a more significant decrease in the hydroelectricity generation when the demands became more pessimistic. The D4 scenario presented greater intensity of the decrease in hydroelectric power with a median of less than −50%. For Três Marias and Itaparica HPP, there was little change of the consumptive demands, with the medians of the anomalies close to −10% and 10%, respectively, in the two climate scenarios (SSP2-4.5 and SSP5-8.5). Figure 16 shows the trend slopes of the CMIP6 models for the future climate scenarios SSP2-4.5 and SSP5-8.5 and the consumptive demands D1, D2, D3 and D4 using the Mann-Kendall Sen test of the energies of the Itaparica, Sobradinho and Três Marias HPP.  Figure 16 shows the trend slopes of the CMIP6 models for the future climate scenarios SSP2-4.5 and SSP5-8.5 and the consumptive demands D1, D2, D3 and D4 using the Mann-Kendall Sen test of the energies of the Itaparica, Sobradinho and Três Marias HPP. From the set of models with a significant trend, there is clear evidence that at Sobradinho HPP, the increase in greenhouse gas emissions suggests a more significant impact on electricity generation since the slope module is higher for the SSP5-8.5 scenario than for the SSP2-4.5. Sobradinho HPP also had the greater magnitude of the slope modules for the most pessimistic consumptive demands scenarios.
Among the models that indicated reductions in electricity generation, the IPSL-CM6A-LR model stands out, reaching the −0.7 slope module in scenario D4. In the Três Marias UH, for the SSP2-4.5 scenario, only the MRI-ESM2-0 model presented a significant trend with a negative slope in all scenarios of consumptive demands. In contrast, in the SSP5-8.5 scenario, only the IPSL-CM6A-LR model presented a significant trend, which was also negative. In the Itaparica HPP, only in the SSP5-8.5 climate scenario and in the D4 demand scenario, the BCC-CSM2-MR and CanESM5 models presented a significant trend with negative declivities close to −0.05.

Conclusions
The impacts of climate change in the coming decades may influence the intensification of extreme events such as floods, droughts and temperature increases. Climate change, added to an increase in consumptive demands caused by population and wealth growth, can put significant pressure on the hydrosystems, especially when amplified by additional new demands from water transfer projects. Such a situation can significantly influence hydroelectric generation. Thus, this study aimed at generating information on the impact of these changes (climate and demand) upon SFRB's electricity production. Decision-makers can use such information to adopt energy policies and measures to mitigate the possible impacts of such scenarios.
In the time series of affluent natural flows observed in the period from 1931 to 2016, the minimum annual flows of the Itaparica, Sobradinho, Três Marias and Retiro Baixo reservoirs were presented during the last six years. Negative trends, according to the Mann-Kendall Sen test, were found for the first two reservoirs. This period of drought over the NBR has already been reported in several studies, in which they relate to the low- From the set of models with a significant trend, there is clear evidence that at Sobradinho HPP, the increase in greenhouse gas emissions suggests a more significant impact on electricity generation since the slope module is higher for the SSP5-8.5 scenario than for the SSP2-4.5. Sobradinho HPP also had the greater magnitude of the slope modules for the most pessimistic consumptive demands scenarios.
Among the models that indicated reductions in electricity generation, the IPSL-CM6A-LR model stands out, reaching the −0.7 slope module in scenario D4. In the Três Marias UH, for the SSP2-4.5 scenario, only the MRI-ESM2-0 model presented a significant trend with a negative slope in all scenarios of consumptive demands. In contrast, in the SSP5-8.5 scenario, only the IPSL-CM6A-LR model presented a significant trend, which was also negative. In the Itaparica HPP, only in the SSP5-8.5 climate scenario and in the D4 demand scenario, the BCC-CSM2-MR and CanESM5 models presented a significant trend with negative declivities close to −0.05.

Conclusions
The impacts of climate change in the coming decades may influence the intensification of extreme events such as floods, droughts and temperature increases. Climate change, added to an increase in consumptive demands caused by population and wealth growth, can put significant pressure on the hydrosystems, especially when amplified by additional new demands from water transfer projects. Such a situation can significantly influence hydroelectric generation. Thus, this study aimed at generating information on the impact of these changes (climate and demand) upon SFRB's electricity production. Decision-makers can use such information to adopt energy policies and measures to mitigate the possible impacts of such scenarios.
In the time series of affluent natural flows observed in the period from 1931 to 2016, the minimum annual flows of the Itaparica, Sobradinho, Três Marias and Retiro Baixo reservoirs were presented during the last six years. Negative trends, according to the Mann-Kendall Sen test, were found for the first two reservoirs. This period of drought over the NBR has already been reported in several studies, in which they relate to the low-frequency natural variability of the positive phase of the PDO and AMO atmospheric phenomena [15,60,61]. However, according to [15] and [61], this positive phase of AMO and PDO is unprecedented and may be attributed to the extreme increase in global SSM due to anthropogenic global warming.
As a consequence of the prolonged drought, in 2016 ANA reduced the flow of the São Francisco River from 800 m 3 /s to 700 m 3 /s [6] and in May 2017 to 600 m 3 /s [25].
This reduction in affluent natural flows continues in both climate scenarios, SSP2-4.5 and SSP5-8.5, from 2017 to 2100 compared to the historical period (from 1901 to 2000) of the CMIP6 models. The set of models indicated reductions in Q90 for most SFRB reservoirs and in the two periods analyzed (2021-2050 and 2051-2100). Except for the Itaparica reservoir, which for the SSP5-8.5 scenario, from 2021 to 2050, the median of models indicated an increase in Q90. For the set of models that presented a significant trend according to the Mann-Kendall Sen test, an increase in greenhouse gas emissions suggests a more significant impact on the affluent natural flows because the slope module is greater for scenario SSP5-8.5 than for scenario SSP2-4.5.
The analyses of the historical time series of flows taken for the consumptive demands from 1961 to 2016 showed an exponential increase in these demands. Irrigation demands presented the highest average annual growth rates, with values of 6.80%, 7.42%, 10.99% and 9.29% for the reservoirs of Itaparica, Sobradinho, Três Marias and Retiro Baixo, respectively. According to [24], about 77% of the current total demand from the São Francisco River is for irrigation, 11% for urban demand and 7% for industrial demand.
Four possible scenarios of these consumptive demands were projected for 2017 to 2051, coupled with the future climate scenario (SSP2-4.5 and SSP5-8.5). For the SSP5-8.5 scenario, the models indicated projections with more intense hydropower reductions than the SSP2-4.5 scenario at the Sobradinho and Três Marias HPP. The exception was the Itaparica HPP, in which the models indicated, for the two future climate scenarios and for most scenarios of consumptive demands, an increase in hydroelectric power. The Sobradinho HPP was the one that showed the highest sensitivity to the future scenarios of consumptive demands, presenting greater intensity in the decrease of hydroelectric energy. In contrast, the demand scenarios were more pessimistic, i.e., scenario D4 presented greater intensity in the decrease of hydroelectric generation with a median of less than -50%, to the detriment of the other scenarios of consumptive demands D1 (no change in demands from 2016) with approximately −30%, D2 with approximately −40% and D3 with approximately −50% in scenario SSP5-8.5.
The possible reduction in water and energy availability and the increase in consumptive demands, especially irrigation demand, reveals that future conflict may intensify among multiple uses. Moreover, possible economic downturn due to water reduction for agriculture and industry, shortage of municipal water supply and crisis in the region's electricity sector are possible future outcomes.
The need for new policies to restrict irrigation growth in SFRB and to prevent and oppose illegal water abstraction is evident. As far as hydroelectric power generation is concerned, with the increase in the need for water to meet consumptive demands, it is necessary to increasingly invest in alternative forms of electricity generation to supply this reduction in hydroelectric generation.
Two possibilities to increase future energy generation can be implemented, resulting in positive and negative climate change feedback. Positive feedback for climate change would be an investment in non-renewable energy, such as generating thermal energy that uses fossil fuels and emits CO 2 into the atmosphere. In 2015, this energy type accounted for 47% of all energy generated in the NEB region. The other possibility facing these scenarios and with negative feedback for climate change would be the massive investment in renewable sources such as wind and solar. In 2018 the wind energy share was already approximately 52% in the region [9]. This type of energy in Brazil, according to [62], is already cheaper than electricity generated from fossil fuel, biomass and nuclear sources.
Overall, this study showed that climate change associated with growing consumptive demands both from existing and new users may significantly affect hydroelectricity production in SFRB. In an uncertain future, proactively planning how the water-energy-food nexus will interact in the next decades is key to promote a sustainable future for our society. Our results support a necessity to invest in renewable energy to account for the hydroelectric production losses and invest in irrigation regulations as the main consumptive user in the region.  Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: The study did not report any data.