Optimisation of Multipurpose Reservoir Operation by Coupling Soil and Water Assessment Tool ( SWAT ) and Genetic Algorithm for Optimal Operating Policy ( Case Study : Ganga River Basin )

Reservoirs are recognized as one of the most efficient infrastructure components in integrated water resources management. At present, with the ongoing advancement of social economy and requirement of water, the water resources shortage problem has worsened, and the operation of reservoirs, in terms of consumption of flood water, has become significantly important. To achieve optimal reservoirs operating policies, a considerable amount of optimization and simulation models have been introduced in the course of recent years. Subsequently, the assessment and estimation that is associated with the operation of reservoir stays conventional. In the present study, the Soil and Water Assessment Tool (SWAT) models and a Genetic Algorithm model has been employed to two reservoirs in Ganga River basin, India in order to obtain the optimal reservoir operational policies. The objective function has been added to reduce the yearly sum of squared deviation from preferred storage capacity and required release for the irrigation purpose. The rule curves that were estimated via random search have been discovered to be consistent with that of demand requests. Thus, in the present case study, on the basis of the generated result, it has been concluded that GA-derived optimal reservoir operation rules are competitive and promising, and can be efficiently used for the derivation of operation of the reservoir.


Introduction
The apprehensions pertaining to the water and their demands in the existing and foreseeable future will continue to be progressively onerous and would turn out to be increasingly intertwined with several evolving segments, viz.farming, industries, and manufacturing [1][2][3].The impacts of alterations in the global climate has intensified furthermore to these challenges; thus, problems pertaining to water shortage demand more attention and additional actions [4][5][6].There is increasing apprehensions regarding the requisite to formulate the balance between water demands and corresponding rivalling sectors, especially in river basins in which the majority of the water is allotted or exhausted.As a consequence of the increase in population and various socio-economic advancement in the recent decades within the Ganga River Basin, the demand as well as the claim for the available water has swollen, and assessing the demand and supply of available water turns out to be essentially urgent [7][8][9].Most tropical and sub-tropical areas, such as India, which is categorized by considerable annual and seasonal variations in precipitation.Such variability steps up the need for the development of infrastructure and the necessity for the management of available water resources in a sustainable and economic manner.The system analysis methodology stands out as the most significant enhancement in the field of management of water resources, and has vast prospective in stipulating the applicable provision for potent and efficient administration in the evolving circumstance [10][11][12].The downstream area of the Ganga River has undergone a tremendous amount of increase in population and intense agricultural and industrial expansion has all induced a considerable increment in the requirement of need of water.The demand for water from the regions outside of this basin is also escalating.Hence, defining techniques that can accomplish a safe, assured, secure, and sustainable supply of water in this region has caught the attention of both society and government [13,14].
In recent decades, considerable study has been conducted exercising optimization of reservoir to ascertain the optimal and efficient operating policies [15][16][17][18], and the stated area has developed into a major topic in examining the sustainability of existing water resources [16,19,20].To curb the problem of inadequate supply of water during lean season, focus has been given on improving management of existing water resources, especially in the optimization of operations and the policies of reservoir [18,21,22].Reservoirs has the ability to boost the dependability of water supply in fostering livelihoods, increasing agricultural yield, and reducing the susceptibility of farmers to droughts.A reservoir could suffice different objectives through water impounding and streamflow regulatory.Regular reservoir operational objectives comprise of the control of flood, generation of hydropower, agricultural, industrial and domestic supplies of water, navigation, recreation, and so on [10,22,23].Moreover, the ecological and environmental issues that can result from the operation of reservoir and diversion of water have also appealed to stakeholders and decision makers [24][25][26].Usually, most of the reservoirs all around the world are accomplished by employing predefined operating policies [17,20,27,28].Operating policies of reservoirs have been aimed for regulating the release of water by keeping in mind the interests of the reservoir stakeholders and decision makers, volume of water impounded, inflows, demands of water, release capacity, and downstream constraints [29][30][31].Loucks and van Beek, (2017) [32] mentioned two types of operating policies in their study; one distinguishes the target water storage for a specific time, which is typically a distribution rule based on single objective storage curve, whereas the second one categorizes the reservoirs storage zones exercising various rule curves that are correlated by specific distribution policies.In addition to an exemplary operating policy should not only enhance the performance of the corresponding reservoir (e.g., reduce deficit in supply of water and boost monetary profits), although should also be able to accomplish the preservation of hoarded water inside the reservoir [17,[33][34][35].However, more often than not there is the circumstance that a predefined policy could not able to accomplish the objectives of a best policy.Thus, the enhancement of the operational policy of the reservoirs is essential [27,36], and the simulation of reservoirs along with the optimization algorithms can be an efficient and effective methodology for fine-tuning the rule curves [16,35].The operation of the reservoir is a large scale multi-objective optimization problem, involving agriculture, hydrology, hydropower, reliability, and environment flows [37][38][39].
The development of operation of reservoir is a complicated problem, which includes multiple decision variables, several conflicting and contradictory objectives together with significant uncertainty and risk [16,28,36].Designing and attaining operating rules for multi-reservoir structures is a challenging job and has been extensively generated during the scientific past of researches that are associated with the management of water resources [30,31,40].Conventionally, operation of reservoir is accomplished based on heuristic measures, comprising storage curves and skewed decisions made by the decision makers to release water from the reservoir conferring upon the existing stage of the reservoir, prevailing water demands, current hydrological conditions, and the period of the year.The expression rule curve has been predominantly utilized to represent the operational policies that outline the best or optimum level of storage and render a tool for release policies i.e., recognised as a function of stored water [41,42].The problem becomes even more complex and intricate when the constraints or the goals that are associated to quality of water are taken into consideration [28,43].In the last decade, an array of approaches has been formulated based not only on conventional stochastic and probabilistic investigations, but also on the data and the estimation of critical hydrologic events, keeping in mind the end goal to build the efficiency and efficacy of the reservoir [44,45].The operation of the reservoir is a non-convex and nonlinear problem, however few of the academicians have employed a linear method to incorporate the methodology of linear programming for recommending their solution [23,46].Reliability is the major concern among the researchers and necessitates that a complete stochastic approach is applied to decipher the problems of reservoir operation, and yet, only some of the researchers have undertaken to addressed this issue [17,47].
A numerous amount of computer-based analytical techniques, from computation to optimization, are propounded for articulating, analysing, investigating, evaluating, and explaining problems pertaining to water resources planning and management [11,13,32].In the initial years, the majority of the optimization methods for the analysis of water resources systems incorporated dynamic programming (DP), linear programming (LP), and several nonlinear programming techniques.Several studies that are associated with water resources have previously applied these techniques complemented by auxiliary methods for a diverse range of problems [32,48,49].Rapid advancements have transpired in these subjects, aided by high-speed computers, which have attributed to its development.During last decades, evolutionary optimization methods, viz.genetic algorithms, have become popular in the application of global optimization in the field of planning and management of reservoirs [16,18,20,22,23,27,35,36,50].While, almost all of the prevalent algorithms are accomplished in undertaking complicated and complex planning and management problems, evolutionary algorithms methods have been able to attract a vast amount of consideration vis-à-vis, their capability in optimizing complex systems [31,34,44].Advancement in high speed computing techniques have led to a vast number of researches, which employed evolutionary algorithms (EAs) to different disciplines of water resources planning, operational and management problems [13,23], involving the calibration of surface and ground-water [51,52], treatment of water [53], and operation policy of reservoirs [16,30,31].Moreover, Abraham and Jain, (2005) [54] indicated that there are numerous advantages in employing EAs since they demand comparatively less understanding regarding the prevailing questions being answered, are little susceptible to the Pareto front's continuity or their shape, are robust and could be easily implemented.The competence of EAs on anonymous objective functions involving several numerical and scientific characteristics, viz.non-differentiable, discontinuous, non-convex, and multimodal [23,31,46], have been accounted in numerous researches.Evolutionary algorithms (EAs) have been comprehensively employed to a number of extents of water resources over the last decade, such as urban drainage and sewage schemes [55], water distribution systems [39,56], hydrologic and fluvial models [1], and water supply and sewage treatment systems [13,44,55].
GAs, which were initially ascribed by Holland, (1975) [57] and further developed by Goldberg, (1989) [58], have been one of those developments that has offered a potent tool for optimization.GAs, analogous to natural selection based on the renowned Darwinian theory, take into account the native genetic operators and the continued existence of the fittest in an effort to obtain the optimal solution in a considered set of solution [10,17,32,56].They have become a handful tool for answering complicated problems of optimization in several different disciplines.Following the introduction of GA and its comprehensive utilisation in the different disciplines of water resources planning and management problems.Over the years, GAs have gained attention in the field of global optimization application to address storage reservoir optimization operation, planning and management [18,30,31,41,43,47,59].Chang et al. (2005) [18] employed a GA-based methodology to explore the optimal operational policy for the releases of water from reservoir, and then utilised these results in a developed ANFIS model.Oliveira and Loucks, (1997) [60], who applied GA to assess operating policies for multi-reservoir systems, establish that it could be successfully utilised in determining the efficient and the successful operational rules.Wardlaw and Sharif, (1999) [37] also applied GAs to answer the reservoir operation problems optimally.Yuan et al., (2008) [61] applied an advanced differential EA for optimal scheduling of hydro generation on a daily basis.Numerous researches have used EAs to manage the reservoir management problems realistically, particularly the scheduling of release from the reservoir and the planning for hydropower generation [10,16,20,23,31,33,[62][63][64].Today, especially for reservoirs operation problems, EAs are acknowledged as suitable tools to support decision-making due to their multiple advantages.
The reality of an ever increasing urban population coupled with socio-economic advancement during the last few years encompassing the Ganga River Basin (GRB), has substantially increased the requirement of water, besides, the boundless need of the water management practices necessitates the study evaluating the water supply security [65,66].Moreover, anthropogenic activities, for example, river operations and extensive regulations have tremendously significantly altered the water budget of the basin.In addition, the ecological and environmental problems that are associated with the operation of reservoir and the diversion of water have likewise appealed to the concern of the society and government [6,67].Moreover, areas in the downstream part of the Ganga River have witnessed ever increasing growth in the population and intense agricultural and industrial advancement, leading to a considerable surge in the demand of water.Therefore, finding modes of water supply to attain a safe, economical, sustainable, and secure supply in this area, have enticed government, stakeholders, policy makers, decision makers, and society.
As explained before, the water demand of this basin is also rising.Consequently, discovering techniques that can accomplish a safe, assured, secure, and sustainable supply of water in this region has catch the attention of both society and government.Nevertheless, the simulation of credible flows for GRB is the primary step before heading towards the optimizing of reservoir operation problem, as it is necessary not only to manage the available water resources of the basin, but also to have a enhanced knowledge of the influence of possible future alterations on the availability of water [9].Hydrological models have been necessary for examining hydrologic processes and their reactions to both endogenous and human influences, although because of the complexity involving the exemplification of complicated environmental conditions and processes, the hydrologic models require to be calibrated prior to its application in order to attain a realistic match [68][69][70].Especially, water yields, inflows, irrigation water demand, and overall water balance of a watershed can be simulated through the hydrologic simulations by the Soil and Water Assessment tool (SWAT) [71,72].
The objectives of this research have been to (1) compute a range of demands of water (e.g., irrigation requirement) in the GRB; and, afterwards, (2) figure out different optimal approaches for two distinctive reservoirs for different priorities to warrant the river's environmental health, while minimalizing the deficit in supply of water and amplifying the generation of hydropower as well.

Study Area
A multi-objective optimization model has been introduced to assist the operation of the reservoir in this study, for which two reservoirs, namely Rihand dam and Tehri Dam, in the Ganga River Basin, India has been chosen to examine the new model.

The Ganga River Basin
GRB encompasses an area of about 1.08 Mkm 2 and it canvases a distance of 1200 km, before confluencing into the Bay of Bengal in the east.The Ganga River basin is located in the northern side of the Indian sub-continent, which is encompassed within the latitudes of 22 • 30 and 31 • 30 North and the longitude of 73 • 30 and 89 • 00 East (Figure 1).The basin area of 0.86 Mkm 2 in India, which is almost equal to the 26 percent of overall geographical area of the India, is distributed between eleven states.The average annual rainfall in the GRB varies from 350 mm at the western end to 2000 mm near the delta at the eastern part.The amount of rainfall received by the basin not only changes throughout the region, but it is also limited to only a few months of the year mainly during the months of monsoon extending from June to October, thus instigating low flow situations in the Ganga as well as its tributaries throughout the dry season of November through May.Therefore, because of both intra and inter annual unevenness in precipitation and consequently in surface runoff, the realistic distribution and allotment of water from the river through reservoirs is necessary for preserving, sustaining, exploiting, and managing water resources in this area.

The Reservoirs
Rihand dam is one of the Asia's largest man-made reservoirs developed.Reservoir's total area is roughly 457 km 2 .Total area of submergence is approximately 46,600 ha (Table 1).The reservoir is a multipurpose and multiyear structure, which commenced its operation in 1962.Thus, the reservoir is of utmost importance to the people of the area.The released water courses towards north across the industrial regions and meets the Son River near Chopan.The water of reservoir is primarily utilised for: The Tehri dam is situated near Tehri town in the region of Kumaon-Garhwal Himalaya in Uttarakhand, and is the highest constructed reservoir in India having a height of 260.5 m (Table 1).The dam creates a reservoir of 4.0 cubic kilometres with a surface area of 52 km 2 .The water of reservoir is primarily utilised for drinking, irrigation, generation of hydropower, and industrial purposes.
The reservoirs operation policy involved some random operational factors to fulfil various demands from different reservoir stakeholders viz., meeting irrigation demand, meeting emergent demand for power, and so on.The water that is released from the Rihand reservoir is recommended on the basis of water availability in the reservoir at the end of monsoon season, irrigation water demand of state of Bihar, and the generation of hydropower and several releases proposed by Uttar Pradesh Power Corporation Limited (UPPCL) (Lucknow, UP, India)/Uttar Pradesh Jal Vidyut Nigam Limited (UPJVNL) (Sonebhadra, India).

The Reservoirs
Rihand dam is one of the Asia's largest man-made reservoirs developed.Reservoir's total area is roughly 457 km 2 .Total area of submergence is approximately 46,600 ha (Table 1).The reservoir is a multipurpose and multiyear structure, which commenced its operation in 1962.Thus, the reservoir is of utmost importance to the people of the area.The released water courses towards north across the industrial regions and meets the Son River near Chopan.The water of reservoir is primarily utilised for: The Tehri dam is situated near Tehri town in the region of Kumaon-Garhwal Himalaya in Uttarakhand, and is the highest constructed reservoir in India having a height of 260.5 m (Table 1).The dam creates a reservoir of 4.0 cubic kilometres with a surface area of 52 km 2 .The water of reservoir is primarily utilised for drinking, irrigation, generation of hydropower, and industrial purposes.
The reservoirs operation policy involved some random operational factors to fulfil various demands from different reservoir stakeholders viz., meeting irrigation demand, meeting emergent demand for power, and so on.The water that is released from the Rihand reservoir is recommended on the basis of water availability in the reservoir at the end of monsoon season, irrigation water demand of state of Bihar, and the generation of hydropower and several releases proposed by Uttar Pradesh Power Corporation Limited (UPPCL) (Lucknow, UP, India)/Uttar Pradesh Jal Vidyut Nigam Limited (UPJVNL) (Sonebhadra, India).

Data and Methodology
The water quantity the optimal operation model for reservoirs based on SWAT includes two segments.The first segment is the coupling of the distributed hydrologic model (SWAT) with the function of floodgates and reservoirs.This segment is the important part and is effective in estimating hydrologic deviation in the available water resources within the basin scale and specifies the essential set of conditions for the optimal reservoir function.The duration of study ranges from January 1960 to December 2014.The data employed involve the observed hydro-meteorology data for the study region, the reservoirs water levels, quantity of the generated hydropower, and surface runoff at different stations.The second segment is the optimal operation model for reservoirs, which would be employed to achieve two purposes, viz.maximizing both the hydropower potential and water resources availability.The framework in Figure 2 demonstrates the schematic figure of the approach that was adopted in this research in order to attain the objectives.

Data and Methodology
The water quantity the optimal operation model for reservoirs based on SWAT includes two segments.The first segment is the coupling of the distributed hydrologic model (SWAT) with the function of floodgates and reservoirs.This segment is the important part and is effective in estimating hydrologic deviation in the available water resources within the basin scale and specifies the essential set of conditions for the optimal reservoir function.The duration of study ranges from January 1960 to December 2014.The data employed involve the observed hydro-meteorology data for the study region, the reservoirs water levels, quantity of the generated hydropower, and surface runoff at different stations.The second segment is the optimal operation model for reservoirs, which would be employed to achieve two purposes, viz.maximizing both the hydropower potential and water resources availability.The framework in Figure 2 demonstrates the schematic figure of the approach that was adopted in this research in order to attain the objectives.For enhancing the reservoirs function, the prerequisite knowledge of the total demand of water is requisite.The collective demand of water comprises of the fundamental demands (e.g., for satisfying ecological and environmental needs, sustaining navigation, etc.), industrial and domestic water needs, and irrigation water demands.For the computation of water that is needed to fulfil irrigation demand, a traditional technique for assessing water demand for crop is employed.For enhancing the reservoirs function, the prerequisite knowledge of the total demand of water is requisite.The collective demand of water comprises of the fundamental demands (e.g., for satisfying ecological and environmental needs, sustaining navigation, etc.), industrial and domestic water needs, and irrigation water demands.For the computation of water that is needed to fulfil irrigation demand, a traditional technique for assessing water demand for crop is employed.

SWAT Model
The SWAT model is a physically based semi-distributed continuous time model that is developed by the United States Department of Agriculture (USDA) to explore the effect of climate, land use changes on water, and sediment in large un-gauged basins.SWAT can operate on a large basin and is able to replicate numerous processes viz.surface runoff in rivers, sediment and so on, on a daily/sub-daily time step [72,73].SWAT simulates several hydrologic procedures, comprising the generation of streamflow, exercising either the Green and Ampt infiltration equation or Soil Conservation Service (SCS) curve number method.For the quantification of evapotranspiration, the SWAT model offers several approaches, namely, Hargreaves, Priestley-Taylor, or Penman-Monteith methods.Groundwater flow, percolation, and lateral flow are evaluated using the mass balance of the underlying system.
SWAT model has been structured for the entire Ganga river basin [9].The basin has been apportioned into 1045 sub-basins and 20,628 Hydrological Response Units (HRUs) employing Digital Elevation Model (DEM) data from Shuttle Radar Topographic Mission (SRTM), land cover data from National Remote Sensing Centre (NRSC), and soil attributes dataset from the Food and Agriculture Organization (FAO) [9].The daily reanalysis and re-gridded weather data from IMD (rainfall, temperature) has been used.Daily rainfall data are at a resolution of 0.25 • × 0.25 • latitude by longitude grid points, whereas temperature data are at 0.5 • × 0.5 • resolution [74][75][76].
Various SWAT parameters that has been employed for calibration purpose, namely, Snow parameters (SMFMX, SMTMP, SMFMN, SFTMP), Elevation band parameters (PLAPS, TLAPS), and Hydrological parameters (CN2, SOL_AWC, ESCO, GW_DELAY, and GW_REVAP) [9,72,73].The CN2 is the curve number parameter, which controls the fraction of water that infiltrates into the soil or comes as surface runoff through overland flow [72,73].SOL_AWC is the available water capacity of the soil layer and is referred to as the plant available water [72,73].ESCO is the soil evaporation compensation factor and it controls the soil evaporative demand that is to be met [72,73].When the value of ESCO gets closer to 0, the model will receive more water from the lower soil level to fulfil the evaporative demand.GW_REVAP is the ground water revap coefficient, which is a SWAT term to describe the movement of water into overlying unsaturated layers as a function of water demand for ET.The GW_REVAP parameter controls the amount of water in the capillary fringe that separates the unsaturated zone and the saturated zone to move upward (to fulfill evaporative demand) [72,73].A larger GW_REVAP value results in a larger transfer rate from the shallow aquifer to the unsaturated zone.GW_DELAY is the groundwater delay time and increase in the value of GW_DELAY, increases the time between water exits the soil profile and enters the shallow aquifer [72,73].SMFMX and SMFMN are the melt factor for snow on 21 June and 21 December, respectively, which account for the influence of density of snow pack on snow melt.SFTMP and SMTMP are the snowfall temperature and the snow melt base temperature respectively.PLAPS and TLAPS are the precipitation lapse rate and temperature lapse rate respectively, which accounts for the change in the precipitation and temperature on account of elevation [72,73].SWAT model simulation has been compared with the observed data at Uttarkashi and Rihand gauge stations provided by Central Water Commission (on a daily basis), employing commonly used statistical parameters (Figure 3).

Reservoir Inflow
The estimation of stream flow coming to the reservoir system is the important factor; hence, the estimation of credible flows for the basin is essential.The inflow coming to reservoir is employed to compute the net irrigation demand and constructing mass balance constraint for the reservoir.Hence, the computation of credible flows for the basin by employing the SWAT model has been necessary.Reservoir inflow data for each month t, It, has been estimated from the calibrated and validated SWAT model [9].The surface runoff simulated for the two reservoirs sub-basin for the 54-year period (1960-2013) has been derived from the SWAT model [9].It should be highlighted that the simulation of SWAT model has been done at a daily time step, which is then aggregated to attain data at a monthly time step to employ in this research (Table 2).

Irrigation Water Demand
The water collected in the Rihand dam irrigates about 2500 km 2 of the agricultural land.The power house has the installed capacity of 300 MW.
The water that is needed to meet the irrigation demand is generally computed by utilising the methodology that was proposed by Brouwer and Heibloem, (1986) [77].The water necessary to meet the irrigation demand leans upon the water demand of the crops, which has been outlined as "the amount of water required to satisfy the loss of water due to evapotranspiration (ETcrop) of a crop that

Reservoir Inflow
The estimation of stream flow coming to the reservoir system is the important factor; hence, the estimation of credible flows for the basin is essential.The inflow coming to reservoir is employed to compute the net irrigation demand and constructing mass balance constraint for the reservoir.Hence, the computation of credible flows for the basin by employing the SWAT model has been necessary.Reservoir inflow data for each month t, I t , has been estimated from the calibrated and validated SWAT model [9].The surface runoff simulated for the two reservoirs sub-basin for the 54-year period (1960-2013) has been derived from the SWAT model [9].It should be highlighted that the simulation of SWAT model has been done at a daily time step, which is then aggregated to attain data at a monthly time step to employ in this research (Table 2).

Irrigation Water Demand
The water collected in the Rihand dam irrigates about 2500 km 2 of the agricultural land.The power house has the installed capacity of 300 MW.
The water that is needed to meet the irrigation demand is generally computed by utilising the methodology that was proposed by Brouwer and Heibloem, (1986) [77].The water necessary to meet the irrigation demand leans upon the water demand of the crops, which has been outlined as "the amount of water required to satisfy the loss of water due to evapotranspiration (ET crop ) of a crop that is free from diseases, raising in fields under non-restricting soil conditions embracing fertility of the soil and availability of water in the soil, thus accomplishing high level of production in the specified cultivating atmosphere" [78].Therefore, in order to compute the water that is needed to meet the irrigation demand aimed at the intended farming land, the subsequent parameters has been computed: crop evapotranspiration, reference crop evapotranspiration, net irrigation, and actual demand of water required for irrigation with attention being given to both the manner of irrigation and the efficiency of application [79].

Reference Crop Evapotranspiration
Brouwer and Heibloem, (1986) [77] proposed a methodology to compute the reference crop evapotranspiration, ET 0 (mm/mon), by employing the available meteorology data of temperature (Equation ( 1)), and the sunshine duration or radiation (Equation ( 1)) from Brouwer and Heibloem, (1986) [77]: where p = Average Daily Percentage of Yearly Daytime, T mean = (T max + T min )/2.T mean = the mean daily temperature, T max = (Sum of all T max values during the month/number of days of the month), and T min = (Sum of all T min = values during the month/number of days of the month).

Crop Evapotranspiration
After the estimation of reference crop evapotranspiration (ET 0 ), the equation for computing crop evapotranspiration (ET crop ) is given below: where K c is the crop coefficient that is associated with the attributes of the crop and the stages of its growth.

Net Irrigation Water Demand
As the difference between effective rainfall and the evapotranspiration of crop has been adopted to compute the required amount of water that is needed for irrigation, it is imperative to assess the effective rainfall first.Furthermore, the effective rainfall indicates the quantity of the rainfall that gets infiltrated inside the soil and subsequently gets evapotranspirated by the crops into the air [78].Moreover, not the entire rainfall that precipitates has been efficient, as some portion of the rainfall might be lost because of surface runoff, canopy intervention, or evaporation [78].However, these hydrological processes could be computed by utilising physically based hydrological models (e.g., SWAT) [9].In the present research, stream flow, lateral flow, and seepage simulated by employing the SWAT model has been adopted, and the effective rainfall could be estimated from the equation below: where R eff is the effective precipitation, and PREC is rainfall reaching the surface of the soil, Q is stream flow, LATQ is subsurface lateral flow (coming from the root zone and meeting the main channel), and SEEP is seepage that is coming out of the bottom of the soil profile.
The net water needed to meet the irrigation demand, Ir net , is computed by using the following equation: Ir net = ET crop − R eff (4)

Irrigation Water Demand
The water that is needed to meet the irrigation demand has been computed based on the requirement of crop for water, effective rainfall, and the efficiency of the irrigation.To compute the irrigation water demand, an efficiency factor has been employed in order to account for the loss in the water that is incurred during the conveyance and the application on the field [78].The efficiency generally comprises of three sections, that symbolise the loss in the water through conveyance, field canal, and field application [78].In this study, for conveyance and canal efficiencies, the value of 0.8 has been used.While, for field application efficiency, the value of 0.6 (corresponding to furrow irrigation) has been adopted.Therefore, the final project efficiency for our study has been taken as 0.4 (0.8 × 0.8 × 0.6).

Estimation of Other Water Demands
Apart from the water that is needed to meet irrigation demands, the water of the reservoirs is used to meet domestic viz.drinking and industrial demands, navigation, and the maintenance of minimum water in the downstream of the reservoir to sustain ecology, environmental, and aquatic habitat.For the present study, other demands, such as domestic and industrial demands, have been contemplated as the constant requirements.The target release for the reservoirs is deliberated to be equal to the irrigation demand, while the storage is considered to be equal to the summation of the dead storage and the other essential demands for the next three months.That is to say, that this would ensure that the three months of other necessary demands is achieved by the reservoir.The ecological streamflow for sustaining aquatic habitat is to provide the minimum streamflow regimen for aquatic.For the present study, the 10% of the annual average flow for the river that was suggested by Tenant (1976) [80] has been taken as the minimum instantaneous flow that is required to sustain short-term survival for aquatic life.

Multi-Objective Optimization Model
Multi-objective optimization problems for the reservoirs involve the simultaneous optimization of numerous non-commensurable and often contradictory objectives.Often, the reservoirs have two primary uses, generation of hydropower and supply sufficient quantity of water to meet irrigation demand.As can be seen that these objectives are generally contradictory and conflicting in nature with each other, an increase in the level of generation of hydroelectric power demands a higher level of water, which means a lesser amount of water releasing to meet irrigation demand.Likewise, the policy makers and decision makers should consider the probable trade-off between maximizing generation of hydroelectric power and gratifying irrigation and navigation requirements before reaching out for the most optimal policy.

Genetic Algorithm Model Formulation
Genetic Algorithms (GAs), as proposed by John Holland [57,58], are a section of evolutionary algorithms (EAs), which mimic courses that are observed in natural evolution to optimize a certain objective function [43,81,82].GAs, analogous to Darwinian natural selection, is a search technique centred around the methodology of natural selection, evolution, and natural genetics, through generations in an endeavour to obtain an optimal solution.
GA lets a gamut comprised of different entities (solution points) and steered random searches to evolve under stated selection policies to the point, which maximizes/minimizes the considered objective problem [57,58].Similar to the other EAs, no derivatives are required by GAs to assess the best solutions.GAs could also be employed to accomplish both local and global investigation within a big search space [43,83].The fundamentals that are necessary for GAs to commence are the representation of parameter that is needed to be maximised, the genetic operator, and the objective function.The decision variables that are necessary for the problem of optimization has been programmed to a single artificial gene within the whole GA set, which can be outlined as one of the many likely answers for the given optimisation problem.The aptness of a chromosome to be considered as a prime solution to a certain optimisation problem is relied upon the degree of the objective function that is represented by it, in order to assess the likelihood of existence.The investigation of GA begins with an arbitrarily generated population of chromosomes initially, and is then combined through genetic operators, namely selection, crossover, and mutation, to refine the assessment of solutions by generating fitter chromosomes through iterations.The selection process is essentially to select chromosomes in the gene pool with a superior fitness values for survival or generating offspring within the population of the subsequent production.Following the evaluation of fitness function, selection has been accomplished via roulette wheel selection [58], which is analogous to the essential tool of the "Survival of the fittest", to form the next generation, with every portion of a roulette wheel being fragmented to fitness values.Therefore, the strings having higher fitness value would possess a greater possibility of being transferred to the next generation.The information of the two likely solutions of the given optimisation problem is reciprocated through crossover between two random positions.Single-point crossover or multi-point crossover is elected randomly.Mutation is a key method, which allows for the introduction of new genetic chromosome to a population.Mutation allows for the modification of a string to a stipulated extent, which could be both negative or positive.Following the mutation, the function's fitness is again computed to obtain the prime solution for the particular set of decision variables.

Genetic Algorithm Model Development
In the present study, the primary aims of the optimization of reservoir model are to minimalize the deficit of water supply, and to maximize the generation of hydroelectric power.The aforementioned objectives of the optimisation problem are mutually conflicting/competitive, as hydropower generation demands higher head in the reservoir to maximize power production and subsequently generate more energy.However, supply of water to meet irrigation and other demands necessitates release of water to meet these demands.
For the present study, three scenarios coupled with three conditions of priority has been devised to simulate the model optimization.For the first scenario, the objective function is to minimize the deficit of water supply only, designated as the supply priority scenario.Whereas for the second scenario, designated as the power priority the objective function is maximization of generation of hydropower only.While the third scenario intents to maximize the advantages of both supply priority scenario and the power priority scenario simultaneously.Therefore, while the first two scenarios have been intended to regulate the reservoir with a specific and solitary objective, the third scenario involves two conflicting/competitive goals having equal priorities.However, basic water demand must be met to warrant the river's environmental well-being, irrespective of the three scenarios.The genetic algorithm toolbox in MATLAB has been used to design the optimal operating policies for managing water reservoir systems.MATLAB has been utilised due to its ability to provide numerous built in auxiliary functions is helpful in function optimization, is totally portable, and it is methodical in numerical computations.

Objective Function for Rihand Dam and Tehri Dam
In this study, the objective function for the proposed GA model (the fitness function) for the Rihand Dam and Tehri Dam is minimalizing the squared deviation of monthly demand for irrigation deficit and the aberration in the objective storage and maximization of hydropower generation.The following equation represents the mathematical equation of the optimization model: To minimalize the deficit in supply of water; where, D t = water demand for irrigation during the month "t', R t = release of water to irrigation demand during the month "t", S T = target storage in Mm 3 , and S t = initial storage during the month "t".
To maximize the energy production, the following equation represents the mathematical equation of the optimization model: where E is the overall energy generated in M kWh (million kWh); p is the coefficient of power generation; R t is the total discharge into the riverbed turbine during the period t in Mm 3 ; and, H t are the net head available to turbine in meters in period t.
The optimization of the reservoir optimization problem has been subjected to the subsequent constraints:

Constraints
The solution of the multiobjective optimisation problem is constrained by several restrictions that are imposed on the decision variable.Constraints may define the gamut of the defined objective function, or on the other hand, can add restrictions on the solution for a given problem.

Mass Balance Constraint
The correlation between the month to month storage is driven through the equation of continuity.The equation of continuity has been specified as, where, Q t = inflow in the reservoir for the month "t", S t+1 =final storage of the reservoir for the month "t", E t = Loss due to evaporation from the reservoir for the time period "t", and O t = Overflow for the reservoir during month "t".

Release Constraint
The water that is released for the irrigation purpose during the contemplated month ought to be equal or less than the net demand of irrigation during the stated month and this constraint is represented by the following equation, Storage Limits The storage of the reservoir during any month must never be larger than the underlined reservoir capacity, and it must never be lower than the defined reservoir dead storage.Arithmetically, the constraint has been represented by the following equation: Over Flow Constraint Over flow constraint manages the condition, once the final storage surpasses the maximum defined limit of the reservoir.Mathematically, the constraint is represented by the following equation: where, E is the maximum quantity of hydropower in M kWh (million kWh) that could be generated; p is coefficient of power production; R t is the discharge into the riverbed turbine for the period t in Mm 3 ; and, H t are the net head that are available to turbine in meters during period t.

Results and Discussion
The reservoir optimization/simulation model that is characterised by Equations ( 5)-( 12) is adopted in the Rihand Dam, and Tehri Dam to obtain the operative policies for simultaneously optimizing several objectives of the generation of hydropower and the supply of water.In the present study of the GA model, stream flow coming into the reservoir system on a monthly basis, and, the monthly demands for irrigation purposes are the two critical input variables.The primary goal has been to estimate the amount of water that is needed to be released to satisfy the overall demand.That is to say that the objective function for the optimisation of reservoir depends upon the releases from the reservoir at every time step, therefore, water discharges must be the decision variable whereupon the GA revolves around.Additionally, the single variable that is unknown in the optimisation of reservoir operation model is the target storage (S t ), and it is computed by applying the GA for every month to minimalize the deficit of supply of water and the maximization of hydropower generation.
For the present study, three scenarios coupled with three conditions of priority (namely, supply priority, power priority, and equal priority) have been devised to simulate the model optimization.Employing the aforementioned three optimisation scenarios and the GA optimization technique, the stipulated reservoir storages are accomplished for the different operating policies.
Figures 4 and 5 pictures the associated target reservoir storage.From Figures 4 and 5, it has been noted that, to meet the demand in the first scenario (supply priority), it requires larger water release from the Rihand and Tehri dam, and consequently lower water level, resulting in the smallest attained final storages.On the other hand, the second scenario demands the impounding of water inside the reservoir at a higher stage to achieve added head and thus generating more energy.That is to say, the reservoir managers and associated authority are inclined to hold the water in the reservoir, subsequently ensuring the highest attained target storages.It can be noted that releases for the purpose of irrigation has reduced considerably in this scenario.Nevertheless, for the scenario of equal priority, the target storages that were achieved after optimized are in between the two scenarios.Although the derived operating rule curves for the reservoirs for three different scenarios serve distinctive objectives, in spite of that, the constraint ensuring that the attainment of basic surface runoff demands enumerated in the reservoir optimization model, can ensure the minimum streamflow in the river, and thus, should be able to warrant the river's environmental health.The other possible and optional operational, storage, and release policies would be able to support the reservoir authority and policy makers in formulating an efficient and sustainable guideline for several conflicting/competing priorities rendered and for diverse inflow scenarios.The approach of multi-objective GA has been thus very much beneficial, in fabricating a set of well-defined solution for the objectives, which are contradictory/conflicting and would ultimately serve for better operation involving a small time of computation.
The average monthly results of reservoir storage that were obtained by the genetic algorithm and historical operations are shown in Figures 4 and 5.It demonstrates that (1) the maximum storage volumes for both the reservoirs for all the three priority scenarios are much smaller than those of historical operations, thus providing enough water in the river to warrant river health; (2) most of the storage volumes obtained from the GA for the supply priority are close to the average value for the historical operations and much greater than the lowest values, thus providing more water to meet irrigation demand; (3) for the power priority scenario, the storage volumes are close to the maximum value for the historical operations.In summary, the results based on optimal solutions that are computed by GA and historical operations implies that by means of the GA optimisation approach, both the water demand and the generation of hydropower could be optimised and should not cause any significant shortage of water in most of the years (Figures 4 and 5).Based on the hydrologic modelling using SWAT and Genetic Algorithm, this study provides reasonable operating policies for the two important reservoirs for the basin.Moreover, the optimized results that were obtained should be competent in delivering better operation choices in response to wet and dry circumstances in water supply.Whenever a drought condition has been anticipated in the coming future, the storage-governed priority measure should be able to strengthen the potential of the reservoir to take care of the emergent requirements of water or prospective water lend.The reservoirs concerned might be holding a greater volume of water by cutting out the release quantity   Based on the hydrologic modelling using SWAT and Genetic Algorithm, this study provides reasonable operating policies for the two important reservoirs for the basin.Moreover, the optimized results that were obtained should be competent in delivering better operation choices in response to wet and dry circumstances in water supply.Whenever a drought condition has been anticipated in the coming future, the storage-governed priority measure should be able to strengthen the potential of the reservoir to take care of the emergent requirements of water or prospective water lend.The reservoirs concerned might be holding a greater volume of water by cutting out the release quantity  Based on the hydrologic modelling using SWAT and Genetic Algorithm, this study provides reasonable operating policies for the two important reservoirs for the basin.Moreover, the optimized results that were obtained should be competent in delivering better operation choices in response to wet and dry circumstances in water supply.Whenever a drought condition has been anticipated in the coming future, the storage-governed priority measure should be able to strengthen the potential of the reservoir to take care of the emergent requirements of water or prospective water lend.The reservoirs concerned might be holding a greater volume of water by cutting out the release quantity for certain dates, however complementing on other days.The higher storage volume should be able to achieve developing drought-response jobs.However, if the available water is above average and no scarcities have been anticipated, hydropower generation-governed priority results in the added generation of hydropower [23,28].The increase in the production of environmental friendly energy from hydropower sectors can also help in mitigating the greenhouse gas emitted from the power supply from thermal power plants and different sources of energy can be substituted by the additional hydro-power generation.Nevertheless, the third scenario of equal priority can be adopted, as both of the approaches for the mitigation of drought and satisfying irrigation demand and increased hydropower generation become important.Several studies have investigated the optimal reservoir operating rules for the reservoirs using genetic algorithm [16,30,31] and have advocated that the optimal reservoir policies that were obtained by employing the genetic algorithm should be better than the currently employed random operational policy for both meeting total water demand for the irrigation purposes and maximizing the annual production of hydropower.
Nowadays, with the rapid growth in the population, and consequently, the increase in the demands, the ecological health of the river and the shortage of water supply have progressively developed into a prominent problem and thus, have been able to attract the concerns of both government and the public.Predominantly, the revamp of the operating policies of the reservoirs could influence the way in which the reservoirs are being operated to meet the demands.However, added unforeseen and serious complications, such as mitigating flood, might stop the reservoir operators from attaining the desired target or final storages.Yet, the endeavour to accomplish the optimal operating policies should help in warranting the sustainability of river ecosystem and diminish water supply deficit, while undertaking hydropower generation and serving added demands.Hence, the optimal final storage that was obtained in this research may provide decision makers and reservoir operators with informative and valuable information that is needed for efficient operation of the Rihand and Tehri reservoir.

Summary and Conclusions
The conventional optimization model offers merely single optimal solution in a particular simulation, and applies a preference vector that necessitates extremely careful alterations to attain the optimal solutions.Nevertheless, the Multi-objective Genetic Algorithm could attain non-dominated solutions to multi-objective problems in a single simulation employing the population-by-population methodology.Genetic algorithms propose a dominant optimization methodology and have popularly established itself in reservoir operation planning and management; however, rising numbers of variables demands large amount of computation, and ends in a very time-consuming evolution amid successive generations, therefore abruptly reducing the likelihood of converging into an optimum solution.In addition, the efficacy and effectiveness of the solution does not depend on a function that is outlined by the operator and decision makers, subsequently these objective functions could be employed precisely without demanding any transformations.
Evolutionary computation methods exhibit huge prospective in optimizing complicated problems.In this research, a Multi-Objective GA approach has been successfully employed for the problem of the optimization of the operation of the multi-objective reservoir.The foremost benefit of the Multi-Objective Genetic Algorithm methodology is its ability in able to discover several pareto optimal solutions in a particular simulation, which has been both efficient and effective at the same time and it supports the operators and decision makers to make fitting assessments and decisions in diverse situations.This research suggests a constrained genetic algorithm for the planning and management of multi-purpose reservoir, which includes the meeting of human requirement and the generation of hydropower requirements.A multi-objective model has been articulated by incorporating hydropower and irrigation as the two conflicting/competing objectives and the multi-objective reservoir optimization model, exercising a combination of the SWAT simulation model and Genetic Optimisation Algorithm, has been applied to develop efficient operation policies for the Rihand reservoir system and the Tehri Reservoir system, in India, in order to minimalize the annual sum of squared deviation of release to meet irrigation demand and target storage, and to maximize the reservoir capacity to generate hydropower.The operating policies of reservoirs (i.e., final storage) conforming to the three priority situations (supply of water only, hydropower generation only, and equal priority) have been developed.Under the new operation policy, optimised by GA, the trade-off concerning meeting the irrigation and other demands and generation of hydropower has been resolved.An optimal solution has been determined, which is proficient enough in sufficing the numerous demands and generation of hydropower, as compared to the contemporary procedures.The balanced optimum solution could uphold a higher level of water in the reservoir, and at the same time, meet the several downstream demands.Additionally, the optimised policy may amplify the generation of hydropower during the flood season, as well as increase the likelihood of accessibility of water for the subsequent dry season.In this study, optimal solutions for the regulations reservoir have been accomplished.The assessment of the results implies that the operation of the reservoir for hydropower generation and water supply can be enhanced considerably.The optimization model that was proposed for the two reservoirs system has been robust and adaptable in serving multi-objectives under different priorities.Thus, this research has been successfully able to exhibit the efficiency and the effectiveness of Multi-Objective Genetic Algorithm for varying multi-objective reservoir operation strategies.It has been exhibited that the technique has been efficient in optimising the rule curves managing the operation of the reservoirs in a multi-objective framework.

Sustainability 2018 , 20 Figure 1 .
Figure 1.Location of the Ganga River Basin in India.

Figure 1 .
Figure 1.Location of the Ganga River Basin in India.

Figure 2 .
Figure 2. Schematic diagram of the methodology adopted in the study.
Irrigation water need (Brouwer and Heibloem, 1986) Determine Objective Function Define Constraints Apply Genetic Algorithm to derive Optimal Operation policy Determine monthly inflow from SWAT model

Figure 2 .
Figure 2. Schematic diagram of the methodology adopted in the study.

Figure 5 .
Figure 5. Target Storage for Tehri Dam for three scenarios.

Figure 4 .
Figure 4. Target Storage for Rihand Dam for three scenarios.

Figure 5 .
Figure 5. Target Storage for Tehri Dam for three scenarios.

Figure 5 .
Figure 5. Target Storage for Tehri Dam for three scenarios.

Table 1 .
Characteristic of the Rihand and Tehri Dam.

Table 1 .
Characteristic of the Rihand and Tehri Dam.

Table 2 .
Analysis of Inflow for the Reservoirs simulated from the Soil and Water Assessment Tool (SWAT) model.

Table 2 .
Analysis of Inflow for the Reservoirs simulated from the Soil and Water Assessment Tool (SWAT) model.