1. Introduction
The ecological flow regime can provide enough instream flow to protect river health and provide sustainable ecological status in river ecosystems [
1,
2]. According to the literature, four main methods have been developed for ecological flow assessment including hydrological and ecological methods [
3]. However, other combined methods might be proposed based on the complexities in the hydrological and ecological processes in rivers. Hydrologic desktop methods are inexpensive and straightforward for assessing ecological flow. For example, Tennant and flow duration curve analysis methods are simple and well-known methods in this regard [
4]. The most well-known and popular hydraulic rating method is the wetted perimeter method that has been developed based on the relationship between the wetted perimeter and river flow [
5]. Due to complexities in the ecological assessment of river ecosystems, it was essential to improve the assessment methods in this regard for greater focus on the habitats’ requirements [
6]. Instream flow incremental methodology (IFIM) developed a habitat simulation method to assess ecological flow in which hydraulic habitat simulation has been introduced as the core of the method [
7,
8]. Protecting ecological flow is necessary to obtain the long-term purposes of water security [
9].
Hydraulic habitat simulation was originally developed in the structure of the PHABSIM software (V4.1) by applying a univariate habitat model to simulate the suitability of hydraulic habitat parameters including depth, velocity and substrate. In fact, this method simulates the suitability of physical factors separately. Then, a mathematical model such as the multiply method is applied to simulate the combined suitability [
10]. Given the inability of the univariate habitat method to simulate interactions between physical factors, other methods have been introduced in this regard that have been reviewed in the literature [
11,
12,
13]. The fuzzy multivariate method or fuzzy logic approach is a robust method in which interactions between physical parameters could be simulated in hydraulic habitat selection by developing verbal fuzzy rules. Considering expert opinions in the development of physical suitability criteria is one of the requirements in the hydraulic habitat simulation due to complexities of the ecological assessment in the ecosystems. Fuzzy hydraulic habitat simulation has been used as the standalone method to assess and optimize ecological flow in recent studies [
13].
Not only is IFIM not a model, but it also is not a rigid or fixed process to assess ecological flow. In other words, IFIM provides general guidelines to manage river ecosystems and the supply of its requirements in which a combination of simulations and judgements might be applied. Developers encouraged users to have innovation and creativity when generating ecological flow regimes in the rivers [
14]. A significant drawback in all the developed ecological flow methods is the lack of an optimization model in the structure of the assessment methods. In fact, water resource engineers face many complexities in the management of water resources in river basins. In other words, experts need an optimization framework to balance the needs of the environment and water demands. Using optimization frameworks in water resource management has been highlighted in the literature [
15]. Linear programming is the simplest optimization method that has been used in water resource management [
16]. However, it is not robust due to the non-linear nature of water resource management, which needs efficient methods [
17]. Hence, non-linear approaches have been recommended as well [
18,
19]. Two aspects should be considered for selecting the optimization method including robustness and efficiency. Hence, evolutionary algorithms have been proposed for the optimization of water resources as efficient and robust methods [
20]. The classic, as well as the new, generation algorithms have been recommended for finding the optimal solutions [
21]. Moreover, they might be classified based on the utilized methods for finding the best solutions. Some algorithms have been inspired by the social behaviour of animals [
22]. In contrast, other algorithms follow natural laws in the optimization process such as the gravity search algorithm that applies the principals of the gravity law for finding the best solutions [
23]. Many animal- and non-animal-inspired algorithms have been used in the water resource management problem for improving planning [
24,
25,
26,
27].
Some previous studies improved hydrologic methods by adding ecological indicators to the model to assess ecological flow [
28]. However, these methods should be improved due to their lack of focus on the simulation of target species in the study area. Moreover, adding an optimization framework to the ecological flow assessment method is required to balance needs in the river basins, especially in the drought condition. Based on these requirements, this study proposes a novel ecohydraulic method in which field studies are linked to the fuzzy hydraulic habitat simulation and hydrological analysis to assess ecological flow in the structure of the optimization framework at the river basin scale. In fact, the optimization method considers two purposes including minimizing ecological impacts and water supply loss in the river basin. This novel method is developed at the river basin scale; this is a strength point. In fact, previous methods are only assessment tools that are mostly useable at the river reach scale. However, the proposed method is an optimization method that can mitigate the loss of water supply and ecological impacts at the river basin scale simultaneously. This novel method can demonstrate how converting assessment methods to simulation-optimization methods at the basin scale is necessary. In fact, water resource engineers face many ecological complexities in the management of river basins that might need integrated frameworks to balance needs fairly.
3. Results and Discussion
The results of implementing the proposed method in the case study are presented step by step in this section. Then, a full discussion of the results will be presented.
Figure 8 displays a schematic view of the selected river reaches to identify the sensitive ecological reaches in which hydraulic habitat suitability should not be reduced considerably compared with the natural flow. As can be observed, five river reaches are selected for the hydraulic habitat simulation in the proposed method. Four tributaries are the main inflow to the river, which are reproduction habitats for the native fish species. Moreover, the selected reach downstream (simulated reach 5 or SR5) is the reproduction habitat for some species that need to migrate to the river for reproduction. In fact, if hydraulic habitat loss is minimized in the selected reaches, the river habitats can be protected. The simulation of hydraulic habitats in all the streams of the river network is ideal. However, it is not possible practically due to the required time and cost for simulating habitats. Hence, highly suitable areas in the selected river reaches could protect river habitats acceptably. This work carried out field observations in the simulated reach based on the proposed method in the previous section.
Table 1 displays results of the field observations. Six main fish species have been recognized in the river basins, as displayed in
Table 1. The normalized fish population of each species has been computed. The highest normalized fish population was considered as the target species for the hydraulic habitat simulation. It seems logical that the upstream tributaries are the habitats for cold-water fishes. For example, the brown trout (
Salmo trutta) is the main species in the upstream tributaries and is one of the most well-known cold-water fishes. Based on the combination of microhabitat observations and expert opinions, verbal fuzzy rules were developed for the selected target species. For example,
Table 2 shows part of the developed fuzzy rules for the
Copoeta capoeta as the target species for SR4. According to the results of the verbal fuzzy rules, velocity might have a considerable effect on the hydraulic habitat suitability. However, the effect of depth and substrate on the hydraulic habitat suitability is undeniable. The results of this study corroborate the previous studies regarding the fuzzy hydraulic habitat simulation that highlighted the remarkable effect of velocity on the hydraulic habitat suitability.
Figure 9 displays the normalized weighted useable area (NWUA) functions as the main output of the hydraulic habitat simulation in the selected river reaches. The NWUA indicates the ideal ecological flow regime in the river. The regression model to develop final functions was used in the optimization system. The L3 regression model was the best correlated model for all the simulated reaches. The biological response of river habitats to the change in discharge in the rivers is not linear and demonstrates the necessity of using ecological-based methods to assess ecological flow. In other words, some methods such as hydrologic desktop methods consider a linear relationship between an increasing rate of flow and biological response in the habitats. In fact, these methods recommend that increasing ecological flow is appropriate in all hydrologic conditions. However, the developed NWUA functions in this study corroborate that hydrologic desktop methods are not reliable for assessing ecological flow.
Another required input for the optimization model is the monthly flow time series of different hydrological conditions. Based on the stream drought index,
Figure 10 displays a computed index in the long-term period that was used to calculate the average monthly time series in different hydrological condition including drought as well as normal. This study applied a long-term period in the past years in which instream flow was approximately the same as the natural flow. Thus, the outputs of the hydrological analysis might be reliable and robust for further applications. As presented, the mass balance method was utilized to generate flow time series in the selected reaches where hydrometric stations were not available.
Figure 11 shows the flow time series at SR5 as a sample of the results of the hydrological analysis.
Outputs of the optimization model are the main results of the developed model. As presented, this study utilized different evolutionary algorithms to optimize ecological flow in this study.
Figure 12,
Figure 13,
Figure 14,
Figure 15 and
Figure 16 display the proposed ecological flow using different algorithms. The performance of the optimization algorithms might be very different in terms of ecological flow optimization in the developed method. However, the outputs of some algorithms are approximately similar. It demonstrates that utilizing a robust decision-making system is an important requirement for finalizing ecological flow regime in the developed method. As presented in the previous section, two indices were applied in the decision-making system including the RMSE for hydraulic habitat loss and reliability index for water supply in the river basin.
Table 3 displays the RMSE for hydraulic habitat loss in the representative reaches in the case study. The performance of different algorithms is not similar in terms of hydraulic habitat loss. The highest RMSE is more than 0.5; this indicates that reducing the habitat suitability might be a serious concern in some simulated conditions. However, the optimization model can reduce the difference in habitat loss in the natural as well as real condition. In fact, a RMSE less than 0.2 might be considered as a perfect performance for the optimization model. Due to the significant difference between the performances of the optimization algorithms in terms of hydraulic habitat loss, the decision for finalizing the ecological flow regime should be based on the output of the FTOPSIS method. Furthermore,
Table 4 displays the reliability index for water supply as another criterion in the decision-making system. It appears that the performance of the selected algorithms is approximately similar in terms of water supply. In fact, the performance of the algorithms is the same in terms of proposing total offstream flow. A lower reliability index in the dry years could be observed for all the algorithms, as seems logical. Changing the hydrological condition has a significant impact on the water supply. For example, the reliability index in the dry years is less than 40%, which is not favourable for the regional water authority. Conversely, the reliability index in the wet years is more than 75%, so is able to support most of the defined water demand in the case study. The weight of importance for the two criteria including the RMSE and reliability index were considered as very high and high, respectively. In fact, the priority of the optimization system is to protect river habitats. Hence, very high importance for the RMSE is logical and acceptable. However, water supply might be important as well. Thus, considering the reliability index to be of high importance might provide a fair balance between ecological requirements and water supply in the case study.
Table 5 displays a rating of the alternatives in the case study as well. The solution for prioritizing the purposes in the optimization model is to determine CI1 and CI2 based on the weight of importance of each term in the objective functions. However, as highlighted, CI1 and CI2 were defined as one in this study, and the weight of importance was defined in the decision-making system to find the privileged optimization algorithm.
Figure 17 displays the final ranking of the alternatives in the case study based on the closeness coefficient (CC) as the final output of the FTOPSIS method. Either PSO or DE is the privileged algorithm were used to assess the ecological flow regime in the case study. The final proposed ecological flow regime is shown in
Figure 18.
Due to using a complex method, a full discussion of the different aspects of the developed method is required. In fact, each system might have some advantages, drawbacks and limitations that should be noticed for use in practical projects. The fuzzy hydraulic habitat simulation used in the proposed method might be an appropriate option in many cases. However, other novel methods have been proposed in the literature that might be useable. For example, adaptive neuro-fuzzy inference systems (ANFIS)-based hydraulic habitat simulation is highlighted in the previous research works due to its robust performance in some case studies. One of the recommendations in the future studies for investigating technical aspects of the proposed method is to use ANFIS-based models instead of fuzzy hydraulic habitat simulation to compare the performance of the models for assessing ecological flow regimes. Using hydrological desktop methods should be stopped in future studies of ecological flow. In fact, these methods might generate incorrect optimal ecological flow due to their lack of ability to simulate the actual response of the river ecosystem to the changing rate of flow in the river.
The core of the proposed method is an optimization model that should be discussed in terms of technical considerations and computational aspects. This question as to why engineers should utilize decision-making systems in the optimization of ecological flow should be highlighted. The main shortcoming of evolutionary algorithms is their lack of ability to guarantee the global optimization, which means that these methods are not able to provide the optimal solution in the domain globally. This issue is highlighted especially in complex objective functions such as the developed function in this study. Thus, utilizing a robust decision-making system is inevitable in the optimization process. In fact, it is assumed that the optimization model is able to find the global optimization in the optimization process. However, it is not possible in practice. Hence, the decision-making system would come to be used as the post-processing step in the assessment of the ecological flow regime. The final output of the FTOPSIS method is only useable in the case study. In fact, defining the criteria and weight of importance is considerably effective on the final ranking of the methods. This study considered two criteria in this study based on the requirements of the case study. However, developing other criteria might be obligatory in other cases. Furthermore, the weight of importance should be considered based on the priorities in the case study. For instance, water supply might be more important than ecological flow in some cases. Hence, recommended optimization algorithms as the privileged algorithm for assessing ecological flow are only applicable for the case study.
Computational complexities are a critical aspect for utilizing the optimization system in water resources management. This term is defined as the required time and memory for finding the optimal solution in the domain. High computational complexities reduce the applicability and efficiency of the optimization model. In fact, a complex optimization system might need considerable running time to find the best solution. This might be a serious challenge for applying the optimization system in practice because it is needed to carry out numerous simulations or to cover a long-term period in practical projects. Hence, the required time and memory might increase drastically. One of the remarkable advantages of the proposed method is its low computational complexities. In fact, this research work used the final output of the fuzzy hydraulic habitat simulation in the optimization system to reduce the computational complexities. As discussed, using ANFIS-based models might be another option regarding hydraulic habitat simulation. One of the solutions is to apply an ANFIS-based model in the optimization directly. It might be applicable and robust in some cases. However, the computational complexities dramatically rise due to opening data-driven models in each time step. Hence, utilizing the proposed fuzzy hydraulic habitat method is advantageous in this regard.
As presented, the defined objective function in this study contains two main terms including ecological flow and water supply. The objective function tries to minimize water supply loss and ecological impacts simultaneously. Another option for developing the optimization model is to use a multi-objective model in which water supply and ecological flow terms could be defined as two objective functions. Then, an optimal solution could be obtained using multi-objective evolutionary algorithms such as multi-objective particle swarm optimization. In fact, aggregating two terms in a single-objective function might have some significant advantages. First, single-objective optimization algorithms might have lower computational complexities compared with the multi-objective algorithms. This study developed a flexible and upgradable system, which means that other types of data-driven models could be added to the system in future studies. Thus, it is important that complexities are minimized in practical projects. The proposed aggregated form is beneficial in this regard. Moreover, many classic and new generation algorithms have been developed in the form of the single objective. However, the limited number of algorithms has been addressed in the form of multi-objective algorithms. As discussed, one of the challenges in the application of evolutionary algorithms is that there is no guarantee for global optimization that indicates utilizing different algorithms is essential. The aggregated function is advantageous in this regard as well. In fact, applying different evolutionary algorithms is possible using the proposed objective function in this study.
The developed optimization model is applicable for river basins in which significant storage structures such as dams are not available. In the case study, no dam was available in the basin. However, many dams have been constructed on rivers in different countries. Thus, when using the proposed method for these basins, there is a need to revise the optimization model. In fact, the reservoir operation models should be combined with the developed method to optimize ecological flow in the basin scale. Reservoirs might be helpful for regulating ecological flow for time steps in which sufficient water is not available in the river. However, the construction of dams might have other significant impacts for the river ecosystem that should be addressed in the assessment of ecological flow.
Hydraulic habitat suitability should be considered in the assessment of ecological flow. Some previous studies demonstrate the considerable effect of hydraulic habitat suitability in the assessment of the ecological impacts in river basins. However, other factors might be effective as well. Other factors were not remarkably effective in our case study, though they should still be addressed in other case studies. For example, water quality might be a challenge for the management of many river basins which is effective in assessing ecological flow. Thus, water quality simulation should be added to the developed optimization system in these cases. Furthermore, sediment transport might change the geometry of the habitats in some cases. It is recommendable to add sediment transport modelling to the developed model in future studies. This work applied an improved IFIM model to assess ecological flow in which the optimization system and hydrologic analysis were added to the ecological flow model. However, other ecological-based methods have been highlighted in the previous research works as well. For example, building block methodology (BBM) is a known holistic method in which ecology, hydrology, groundwater issues and hydraulics could be considered in the assessment of ecological flow. This method originally could not be linked with optimization models. Utilizing this method might be advantageous in some cases. It is recommendable to focus on linking the BBM method with the developed optimization model in future studies.
The readers should consider that the daily flows are more important for aquatic life because the monthly flows can mask short periods of total lack of water in the riverbed. However, a daily approach is difficult in terms of computational effort as well as practical implementation. Hence, this research work applied the monthly approach. The environmental managers can design daily environmental flow regimes considering the recommended monthly flow regimes through the hydrological trends in each case study.
Climate change could be added to the optimization system as well. In fact, previous studies corroborate that the impact of climate change might alter stream flow significantly. Thus, the mean monthly flow in different hydrological conditions might be changed in future periods. The proposed method is advantageous in this regard. In fact, modified monthly flow due to the impact of climate change could be considered in the optimization model and might be helpful for long-term management of river basins. Moreover, connected river–wetland ecosystems might need other requirements in the assessment of ecological flow that should be added to the proposed system.
Hydrological desktop methods such as Tennant recommend 10% of mean annual flow as the minimum ecological flow in the rivers. However, this study demonstrates that 50% of mean monthly flow is the minimum assessed ecological flow. In contrast, some methods such as hydraulic rating methods might determine ecological flow without considering hydrological requirements in river basins. Thus, the proposed method is robust in terms of ecological and hydrological considerations. In fact, using integrated ecohydraulic methods is a reliable method to assess the ecological flow regime that should be utilized in the optimization model. In other words, assessment of ecological flow should be converted to the optimization of ecological flow. This study proposed a novel concept for assessing ecological flow by linking hydrology analysis, ecological modelling and optimization of water resources in river basins. This method is able to make the decision process easier in the basins. In fact, the lack of an optimization system might escalate conflicts of interest. In the proposed objective function, defining the weight of importance of each term (CI1 and CI2) based on the negotiations among stakeholders can be helpful to minimize conflicts. Finally, it should be discussed that most of the previous studies developed ecological flow models at the river reach scale, which might not be appropriate for use at the river basin scale because it is necessary to use acceptable mean habitat suitability at the watershed scale, which might not be assessable using the previous methods.