Dynamic Rule Curves and Streamﬂow under Climate Change for Multipurpose Reservoir Operation Using Honey-Bee Mating Optimization

: Climate change in the watershed above the reservoir has a direct impact on the quantity of streamﬂow that enters the reservoir and the management of water resources. Developing effective reservoir rule curves helps reduce the risk of future failures of water resource management. The purpose of this study was to analyze the inﬂuence of climate change on the volume of streamﬂow entering the Ubolratana Reservoir, Thailand during the years 2020–2049 with climate simulations from the CIMP5 model under RCP4.5 and RCP8.5 scenarios. SWAT models were used to forecast future reservoir streamﬂow quantities. Moreover, suitable reservoir rule curves using the Honey-Bee Mating Optimization (HBMO) were developed and the effectiveness of the new rule curves was assessed. According to the research ﬁndings, the average yearly streamﬂow in the future apparently grew from 32% in the base years (2011–2019) and 65% under the RCP4.5 and RCP8.5 scenarios, respectively. It was discovered that the average monthly streamﬂow was higher in the rainy season than in the dry season. Both of the projected situations have a form compatible with the present rule curves in the section of the new reservoir rule curves generated with the HBMO. Furthermore, the newly constructed rule curves may allow the reservoir to keep more water during the rainy season, thereby assuring that there will be adequate water during the following dry season. Additionally, during the dry season, the reservoir was able to release more water that would be able to reduce the water shortage, indicating that it was able to effectively reduce the amount of water shortage and average overﬂow under RCP4.5 and RCP8.5 situations.


Introduction
Uncertainty has a direct influence on the understanding of hydrology and water resource cycles caused by global climate change, as well as the growing frequency and intensity of droughts and floods throughout the world; these events are jeopardizing the management and development of water resources to meet global demands in all industries, making management more complex and difficult. For the past two decades, climate change has had a global impact on water resource management. Several study groups have sought to create ways for controlling water at its sources in order to deal with the fluctuation of supply sides and demand sides. The majority of such studies have evaluated the consequences of future climate change based on prediction findings from climate models combined with hydrological models to analyze impacts on water allocation efficiency for consumption [1], irrigation [2,3] hydroelectric power generation [4], and procurement of new reservoirs in the future [5].
In Thailand after the Great Flood of 2011, numerous watershed areas experienced drought between 2012 and 2019. The primary reason for this is that rainfall was below normal [6]. Many rivers' average discharge was lower than usual [7]. Government agencies must implement campaign initiatives to encourage consumers and farmers to consume water insufficiency problem. Water intake to reservoirs has been lower than the average amount. In contracts, in certain years, the volume of water flowing into the reservoir surpasses the storage capacity during the rainy season. The water must be drained onto the downstream side, causing floods in residential and agricultural regions. As a result, when Ubolratana Reservoir has to develop suitable and efficient water distribution criteria, taking into account the diversity of hydrological circumstances in the past, present, and future together with the application of evolutionary optimization techniques to create more efficient rule curves. This would be expected to be of great benefit for water resource management.
In the past, the consideration of improving the reservoir rule curves of Ubolratana Reservoir, and the other reservoirs in Northeastern Thailand was a case study based on climate change forecasting from the AR4 models [18,30]. This research draws on climate forecasting data from the CMIP5 model based on the RCP4.5 and RCP8.5 scenarios that use bias correction to be more accurate, including there are various types and different model resolutions. The integration of SWAT hydrological models into the analysis of streamflow conditions has not been previously studied, and the same applies to experiments that link these models to the development of the optimal reservoir rule curves with the HBMO technique. Consequently, the expected outcome of the study is the optimal rule curves, appropriate outcome for the climate change situation and the variation on streamflow in many cases.
The purpose of this research was to use the CMIP5 and SWAT models to examine how global climate change affects the quantity of streamflow input into the Ubolratana Reservoir, as well as to improve the reservoir rule curves by employing the approach of the HBMO and considering the objective function, which is to minimize the quantity of water that is scarce and the amount of water that overflows the reservoir, respectively. The results of this study were predicted to be useful in predicting water scarcity and extreme water circumstances for flexible water management, provided as decision-support information for stakeholders to use as information for climate change policy planning and evaluation of water allocation guidelines to assist future activities.

Research Area
The research site was Ubolratana Reservoir in Ubolratana District, Khon Kaen Province. The study focused on five watershed areas; Lam Pha Niang, Lam Nam Phue, Upper Lam Nam Phong, Lam Nam Choen, and Lam Nam Phrom, all of which are tributaries of the Chi River Basin in northeast Thailand (Figure 1), with a total water intake area of around 12,000 square kilometers. The reservoir is a rock-fill dam with a clay core with a height of 2 m. The dam crest is 185.00 m above sea level. The basin receives an average of 2470 MCM of water each year. The normal water storage capacity is 2431.3 MCM, with a reservoir area of 370 square kilometers. The main functions of the reservoir are for generating electricity with an annual power generation capacity of approximately 56.1 million kilowatt-hours, irrigation covering an area of approximately 480 square kilometers, flood relief, fisheries, and intercity transportation travel.  [31]. For this study, 10 CMIP5 models were selected by the investigators: MIROC_ESM, BNU, CanESM, MI-ROC5, FGOALS_g2, CESM1_CAM5, GFDL, EC_EARTH, CCSM4, and FGOALS_s2 [32,33]. The data used in the global climate change analysis were supported by the Hydro-Informatics Institute (HII) (Public Organization), which revealed that there are a wide variety of models that can be applied (more than 15 models). However, when comparing the model's data with the measurement stations in the study area, (especially rainfall data) and ranked based on the lowest tolerance. It was found that the models used in this study were among the 10 models with the lowest inaccuracies and were used in this study. Then, in the streamflow analysis, only the climate data from five of the best models were selected. For ease of use, the HII, which has downscaled the data model to a 5 × 5 square kilometer grid. Base year climate data in the study areas used the data for 9 years between 2011-2019, and climate forecasting data from 30-year models between 2020-2049.

Data Bias Correction
The Gamma-Gamma transformation approach was used in this study to correct for rainfall inaccuracy from the GCM. For this study, climate data, particularly precipitation data, courtesy of the Hydro-Informatics Institute (HII), is the agency that produces and distributes data for use in climate change studies in Thailand. This agency has identified the Gamma-Gamma transformation method to mitigate discrepancies in rainfall data. In addition, HII has published a study that applied this method to study the impact of climate change in Thailand on agricultural water demand [34]. In addition, Sharma (2015) has also chosen this method to study rainfall in western Thailand, which found that the   [31]. For this study, 10 CMIP5 models were selected by the investigators: MIROC_ESM, BNU, CanESM, MIROC5, FGOALS_g2, CESM1_CAM5, GFDL, EC_EARTH, CCSM4, and FGOALS_s2 [32,33]. The data used in the global climate change analysis were supported by the Hydro-Informatics Institute (HII) (Public Organization), which revealed that there are a wide variety of models that can be applied (more than 15 models). However, when comparing the model's data with the measurement stations in the study area, (especially rainfall data) and ranked based on the lowest tolerance. It was found that the models used in this study were among the 10 models with the lowest inaccuracies and were used in this study. Then, in the streamflow analysis, only the climate data from five of the best models were selected. For ease of use, the HII, which has downscaled the data model to a 5 × 5 square kilometer grid. Base year climate data in the study areas used the data for 9 years between 2011-2019, and climate forecasting data from 30-year models between 2020-2049.

Data Bias Correction
The Gamma-Gamma transformation approach was used in this study to correct for rainfall inaccuracy from the GCM. For this study, climate data, particularly precipitation data, courtesy of the Hydro-Informatics Institute (HII), is the agency that produces and distributes data for use in climate change studies in Thailand. This agency has identified the Gamma-Gamma transformation method to mitigate discrepancies in rainfall data. In addition, HII has published a study that applied this method to study the impact of climate change in Thailand on agricultural water demand [34]. In addition, Sharma (2015) has also chosen this method to study rainfall in western Thailand, which found that the Gamma-Gamma transformation was more effective in improving rainfall frequency and intensity compared to other methods [35]. The concept of this method is to correct for discrepancies caused by frequency and rainfall between GCM and measurement stations in the base year by creating a cumulative distribution function (CDF). This leads to the creation of appropriate Gamma parameters, with the functionalities and key parameters as shown in Equations (1)-(4).
where α is the shape and β is the size of the data from the GCM and base year monitoring stations at the selected locations to be gamma distribution. xTrunc is the amount of rainfall from CDF treated with the Gamma parameters, which are developed in Equation (2) for Equation (3). The α and β values were calculated by applying the maximum likelihood estimation method to calculate the daily precipitation from the inverse-adjusted GCM as shown in Equation (4).

SWAT Hydrological Model
The SWAT (Soil and Water Assessment Tool) model was created to aid in the management of water resources, and it was utilized in the evaluations for estimating the impact of water resource management and water pollution in watersheds and large basins [36], the quantity of streamflow that has changed, the amount of sediment and water quality in streams affected by changes in land use and climate in both past, present and future projections [37], which could be divided into distinct stages of watershed processing. For example, in the main watershed, sub-watershed zones are being created. Calculations that demonstrate outcomes daily and at extended intervals are also included. This considers variables from hydrological processes with the water balance equation as in Equation (5).
where SW t is the final soil water content; SW 0 is the initial soil water content, t is the time (days), R day is the precipitation (mm) on the day i, Q surf is the surface water content on the day I, E a is the evaporative transpiration amount on the day I, W seep is the amount of water seeping into the basement on the day i, and Q gw is the amount of groundwater returning to the stream on the day i.

Data Input
In the implementation process, the SWAT method requires the import of basic physical data, including a digital elevation model (DEM) with elevation values between 90 to 1596 m (MSL). The watershed area has a slope from the west (mainly mountains and upstream forests) to the eastern lowland area where the Ubolratana Reservoir is located (see Figure 1). As for the soil type map (Figure 2a), it indicates that more than 50% of the soil is clay, which is in the eastern lowland, followed by clay loam soil, which is mainly in the eastern lowland of the study area. The types of land use in the study area were mostly agricultural areas. It was found that the use of land for rice farming which is most distributed in the eastern lowland area, combined with sugarcane and cassava plantation in the central area of the basin. In the west, most areas are watershed forests. The land use spatial distribution map is illustrated in Figure 2b.
Station E68A (Lam Pha Niang Basin), E29 Station (Upper Phong Basin), and E85 Station (Lam Chuan River). Basin). These data are from 2011-2019 supported by the Royal Irrigation Department of Thailand. The data used for evaluating the effectiveness of the SWATcomputed results for the different types, intervals, scales, and data sources used in this study are summarized and shown in Table 1.  SWAT-CUP (SWAT Calibration and Uncertainty Procedure) is a SWAT-compatible model. When compared to the old approach of manual correction by trial and error, the SWAT model's sensitive variable analysis, calibration, and validation procedures have more flexibility and take less time. The outcome of altering the sensitivity variable will Daily climate data includes rainfall, temperature, humidity, wind speed, and solar intensity. Daily rainfall data were collected from 9 rain gauge stations distributed in the study area and 1 climate station (Khon Kaen station) located in the southeastern part of the watershed, as shown in Figure 1. There are 4 stations of streamflow and sediment data, of which 3 stations are located in the watershed areas above the Ubolratana Reservoir, are Station E68A (Lam Pha Niang Basin), E29 Station (Upper Phong Basin), and E85 Station (Lam Chuan River). Basin). These data are from 2011-2019 supported by the Royal Irrigation Department of Thailand. The data used for evaluating the effectiveness of the SWAT-computed results for the different types, intervals, scales, and data sources used in this study are summarized and shown in Table 1. SWAT-CUP (SWAT Calibration and Uncertainty Procedure) is a SWAT-compatible model. When compared to the old approach of manual correction by trial and error, the SWAT model's sensitive variable analysis, calibration, and validation procedures have more flexibility and take less time. The outcome of altering the sensitivity variable will serve as a guide for the best calibration and adjustment of the solution(s) between the SWAT generated results and the station data.  [38]. For this study, the use of the SUFI-2 technique was selected to apply in the operation. The SUFI-2 technique is uncertainty analysis consisting of predictive P-factors representing the actual measured values that appear in the simulation results for 95% of the uncertainty of the simulation. The prediction (95% prediction uncertainty; 95PPU) and R-factor are calculated as the ratio of the mean amplitude range of the 95PPU to the standard variance of the actual data. The calculated 95PPU values were positioned at 2.5% and 97.5% of the cumulative probability distribution of the variables considered. Using Latin hypercube sampling [38] as this technique requires the least number of sensitivity variables but can produce the best results compared to other methods [39]. Eight parameters from the most vulnerable model types were chosen for examination in this study. Eight parameters from the most vulnerable model types were chosen for examination in this study. The results of the modification of the parameters that calculated streamflow from the model closest to the data from the measurement station are shown in Table 2. Then, the results were compared with the data from the measurement station, and the efficiency was assessed using two statistical indices to check the accuracy of the results [40], which showed the level of accuracy of the monthly streamflow comparison results. It is divided into four levels as shown in Table 3 [41].  (6), is between 0-1, with values greater than 0.6 indicating that the two data are correlated at a level of reliability.

2.
The Nash Sutcliffe efficiency (NSE) coefficient, as shown in Equation (7), is between −∞ and 1, with values greater than 0.5 indicating that the two data are correlated at a level of reliability.
where n is the total number of data. Q oi is the i-order value, Q oa is the mean from all measurements, Q si is the i-order model, Q sa is the i-order value from all models, Q s is the calculated value from the model, and Q o is the measurement value. The HBMO Algorithm is a hybrid search algorithm based on bee mating behavior. The biological bee breeding process is transformed into a mathematical modeling program. As a result, the phases in the adjustment process were properly outlined. Mating is the first step in algorithm development. Every queen bee makes a flight based on her power and speed throughout each mating flight. Equation (8) determines the likelihood of mating between individual male bees and queen bees. The likelihood of mating is high during the start of the mating flight when the queen bee's velocity is high, or when a male bee is sufficiently numerous to mate, the probability of mating is high.
After the movement of the queen bees or after mating, energy, and speed decrease according to Equations (9) and (10). When all queen bees have completed a pairing flight, they begin to breed to achieve the required number of embryos. The queen bees are selected in proportion to the queen bee's fitness and are artificially inseminated with sperm randomly selected from the queen bee's sperm sac. The worker bees would be selected in proportion to their fitness to be used to improve larval outcomes. After the embryos were born, they would be sorted according to their fitness. The best larvae replace the worst queen bees until there are no better embryos than any queen bees. The remaining larvae are then killed and new matings begin until there is a perfect mating. All predetermined will be completed or meet converging criteria [42].
where Prob (Q, D) is the probability of mating between the male bee D and the queen bee Q or the probability of successful mating; ∆( f ) is the difference between the male bee's fitness (f (D)) and the fitness of the queen bee (f (Q)); S(t) is the speed of the queen bee at the time.
where E(t) is the queen's energy; S(t) is the queen's speed; α is a factor ∈ [0, 1] and γ is the amount of energy reduction after each transition.

Water Equilibrium Simulation Model
The models HEC-3, HEC-5, and HEC-RAS were used in a simulation study of the reservoir system in each watershed [43]. Water balance principles were used. In this study, a simulation model of the reservoir system was created by using the same principles as in the above model, to facilitate connection with the Honey Bee Mating Optimization and begin calculating the water balance of each reservoir. To begin calculating the water balance of each reservoir from the rule curves, the initial storage volume of the reservoir was set at full capacity or the maximum storage level; the discharge volume could be calculated following the Standard Operating Rule as shown in Figure 3 and Equation (11). Then, the available water cost of the reservoir could be calculated for the next month with the principles of the water balance equation as shown in Equation (12).
where R υ,τ is the amount of water discharged from the reservoir during the year υ in the month τ (τ is 1 to 12 referring to January to December); D τ is the demand for water at the bottom of the basin during month τ; x τ is the lower boundary of the rule curves of the month τ; y τ the upper boundary of the rule curves of the month τ; and W υ,τ is the amount of original water level available in the basin of the month τ.
where S υ,τ is the amount of water stored in the reservoir at the end of the month τ; Q υ,τ is the average streamflow in the month τ; Eτ is the evaporation loss in the month τ; and DS (dead storage) is unused storage volume.
where Rυ,τ is the amount of water discharged from the reservoir during the year υ in the month τ (τ is 1 to 12 referring to January to December); Dτ is the demand for water at the bottom of the basin during month τ; xτ is the lower boundary of the rule curves of the month τ; yτ the upper boundary of the rule curves of the month τ; and Wυ,τ is the amount of original water level available in the basin of the month τ.
where Sυ,τ is the amount of water stored in the reservoir at the end of the month τ; Qυ,τ is the average streamflow in the month τ; Eτ is the evaporation loss in the month τ; and DS (dead storage) is unused storage volume. The reservoir rule curves were generated using the HBMO Algorithm Optimal Solution in this study. In the instance of shortage frequency, the target function for determining the solution was the least average shortage, as illustrated in Equation (13).
where is the length of the original water quantity data set; ℎ is the amount of water shortage in the year (The amount of water released is less than the water demand target) The reservoir rule curves were generated using the HBMO Algorithm Optimal Solution in this study. In the instance of shortage frequency, the target function for determining the solution was the least average shortage, as illustrated in Equation (13).

Min(Aver Sh
where n is the length of the original water quantity data set; Sh v is the amount of water shortage in the year v (The amount of water released is less than the water demand target).

Reservoir Rule Curves Efficiency Evaluation
By analyzing the frequency of occurrence of an incident, the rule curves assessment was set to evaluate two parts: water scarcity and excess release water with mean and maximum values of Magnitude and Duration through the performances of the test rule curves with future monthly streamflow scenarios from 2020 to 2049. Changes in greenhouse gas emissions are RCP4.5 and RCP8.5, which are two different types of RCP.

Model Performance Assessment
Evaluation of model performance assessed the accuracy between the calculation of streamflow from the SWAT model calculated from the average monthly streamflow volume during 2011-2019 and the streamflow data from 4 measurement stations in the study areas, namely E68A Station (Lam Pha Niang River Basin), E29 Station (Upper Phong River Basin), Ubolratana Dam Station, and E85 Station (Lam Nam Choen River Basin) in the same period. The model's performance was evaluated using an index of R 2 ranging from 0.62-0.88 and NSE between 0.50-0.81, which were both within the acceptable accuracy range as shown in Table 4. Comparative results of streamflow volumes from the SWAT model and streamflow data from Ubolratana Dam Station are shown in Figure 4. The average annual streamflow from the SWAT model is 5147.34 MCM and that of the measurement station is 2385.56 MCM.

Model Performance Assessment
Evaluation of model performance assessed the accuracy between the calcul streamflow from the SWAT model calculated from the average monthly streamfl ume during 2011-2019 and the streamflow data from 4 measurement stations in th areas, namely E68A Station (Lam Pha Niang River Basin), E29 Station (Upper Phon Basin), Ubolratana Dam Station, and E85 Station (Lam Nam Choen River Basin same period. The model's performance was evaluated using an index of R 2 rangin 0.62-0.88 and NSE between 0.50-0.81, which were both within the acceptable a range as shown in Table 4.  Figure 4. The average annual stre from the SWAT model is 5147.34 MCM and that of the measurement station is MCM.

Forecasting of Future Streamflow Volumes
Forecasted future streamflow from 2020 to 2049 were expected to be impa climate change based on the CIMP5 model under the RCP4.5 projection case. In tot was a 32% increase in the average annual streamflow in the future. With the MIRO model, the streamflow volume was likely to increase to a maximum of 4734.9  (Table 5).
Climate change was projected to influence future streamflow levels between 2020 and 2049, according to the CIMP5 model under the RCP8.5 forecast. The results showed that the average annual streamflow across all models tended to increase.  (Figure 7). Looking at the seasonal average monthly streamflow volumes, the trend of change in average monthly water volume was similar under the RCP4.5 projection case but had a greater proportion of increase. Overall, the average monthly streamflow volume increased during the rainy season by 3551.80 MCM (57.32%) and by 401.32 MCM (213.81%) in the dry season. The increase was significant in both the rainy and dry seasons compared to the other models (Table 5), with a significant increase in percentage in October (Figure 8).   (Table  5).       (Figure 7). Looking at the seasonal average monthly streamflow volumes, the trend of change in average monthly water volume was similar under the RCP4.5 projection case but had a greater proportion of increase. Overall, the average monthly streamflow volume increased during the rainy season by 3551.80 MCM (57.32%) and by 401.32 MCM (213.81%) in the dry season. The increase was significant in both the rainy and dry seasons compared to the other models (Table 5), with a significant increase in percentage in October ( Figure 8).

Optimal Reservoir Rule Curves by HBMO Algorithm
The findings of the Ubolratana Reservoir rule curves generated with the HBMO Algorithm approach based on the CIMP5 climate change impacts of 5 models under RCP4.5 and RCP8.5 projection cases were compared to the present Ubolratana Reservoir rule curves. The rule curves in both predicted situations were discovered to be identical to the existing rule curves. However, from July to September, the newly developed upper rule curves were higher than the current rule curves. This effected an increase in the amount of water stored in the reservoir, resulting in a sufficient water supply for the next dry season. In the upper rule curves of the two forecast cases, the shape corresponded to the current rule curves, but the lower rule curves developed lower than the current ones during the dry season from December to April. This means that the reservoir can release more water than with the current rule curves. It can reduce water scarcity, making it possible to respond to water users in irrigated areas (Figures 9 and 10). According to recent study, applying the Harris Hawks Optimization (HHO) technique for searching in the Ubolratana reservoir, Thailand, the optimal rule curves with the HHO technique was similar to the current rule curves. The upper rule curves developed were higher than the current rule curves throughout the rainy season, allowing for additional water storage at the end of the rainy season [44].

Reservoir Rule Curves Efficiency Evaluation
The purpose of evaluating the efficiency of reservoir rule curves is to test the func-

Reservoir Rule Curves Efficiency Evaluation
The purpose of evaluating the efficiency of reservoir rule curves is to test the functions of the rule curves in order to know the results that could support the changing water situations due to various uncertainties, whether in past periods or for scenarios that may occur in the future. The assessment of rule curves had two parts, namely, water shortage and excess release water by assessing the frequency of occurrence of an incident through mean and maximum values of Magnitude and Duration.
We evaluated the efficiency of the current reservoir rule curves and the reservoir rule curves obtained from future streamflow during 2020-2049, which yielded five CIMP5 models of climate change under the RCP4.5 scenario. In all models except the MIROC5 model, the reservoir rule curves were able to lower the mean water deficit and mean overflow when compared to the present rule curves. Under the RCP4.5 scenario, the reservoir rule curves from the MIROC_ESM model were the most efficient ones in reducing mean water deficit and mean overflow when compared to the reservoir rule curves in other models (Table 6). Under the RCP8.5 scenario, the results showed that the reservoir rule curves in all models were able to reduce the average water shortage compared to the current rule curves. Moreover, the reservoir rule curves from the MIROC5 model could also help reduce the over-average water flow. The efficiency evaluation indicated that the reservoir rule curves from the MIROC5 model were able to reduce the average water shortage and average overflow the best when compared to the reservoir rule curves of all models (Table 7).

Conclusions
There were two primary objectives of this research. The first was to investigate how global climate change has affected the quantity of streamflow that flows into the Ubolratana Reservoir in the years 2020-2592. Second, these modifications will be utilized as data for improving the suitable reservoir rule curves using the HBMO algorithm approach, as well as evaluating the effectiveness of the newly designed reservoir rule curves.
The results of this study showed that future streamflow data are based on the SWAT model. The forecast years 2020-2049 were projected to be influenced by climate change from the CIMP5 model, according to the findings of this study. Both RCP4.5 and RCP8.5 were expected to rise under the anticipated conditions. Under RCP4.5 and RCP8.5, the future overall average annual streamflow will rise by 32% and 65%, respectively. The MIROC_ESM model had the highest average annual streamflow compared to other models. However, there is a different study (FGOALS_g2 model, under the RCP4.5 forecast case), which indicates that the future annual mean streamflow tends to decline. When we considered the average monthly streamflow volume in the future according to seasons, it was found that the trend of change in streamflow volume was consistent with both under the forecasting cases. The average monthly streamflow volume was expected to increase markedly during the wet season (August to November) and at the beginning of the dry season (December).
Then, the Ubolratana Reservoir rule curves developed by HBMO Algorithm was created. There were five CIMP5 climate models under the RCP4.5 and 8.5 forecast cases, for which the developed rule curves were shaped in accordance with the current rule curves. Moreover, the developed rule curves could also allow the reservoir to hold more water during the rainy season. This should ensure that there will be enough water in the next dry season. In addition, during the dry season, reservoirs will be able to release more water, thereby reducing water scarcity. Finally, the future rule curves in the reservoir as a result of the climate change examined in this study would be able to answer the objective functions, which is to acquire the least average water scarcity amount. The rule curves will also be rated for their efficiency in reducing water scarcity and overflow compared to the current rule curves.