Next Article in Journal
What Determines Innovative Performance of International Joint Ventures? Assessing the Effects of Foreign Managerial Control
Next Article in Special Issue
Artificial Neural Network Optimized with a Genetic Algorithm for Seasonal Groundwater Table Depth Prediction in Uttar Pradesh, India
Previous Article in Journal
Urban Freight Last Mile Logistics—Challenges and Opportunities to Improve Sustainability: A Literature Review
Previous Article in Special Issue
New Insights on Flood Mapping Procedure: Two Case Studies in Poland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Environmental Flows Assessment in Nepal: The Case of Kaligandaki River

1
Department of Civil Engineering, Khwopa College of Engineering, Bhaktapur 44800, Nepal
2
CERIS, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais 1, 1049-001 Lisbon, Portugal
3
College of Water Conservancy and Hydropower Engineering, Hohai University, Gulou District, Nanjing 210098, China
4
Department of Sanitary Engineering and Water Management, University of Agriculture in Krakow, St. Mickiewicza 24–28, 30-059 Krakow, Poland
*
Authors to whom correspondence should be addressed.
Sustainability 2020, 12(21), 8766; https://0-doi-org.brum.beds.ac.uk/10.3390/su12218766
Submission received: 10 September 2020 / Revised: 13 October 2020 / Accepted: 19 October 2020 / Published: 22 October 2020

Abstract

:
Environmental flow assessments (e-flows) are relatively new practices, especially in developing countries such as Nepal. This study presents a comprehensive analysis of the influence of hydrologically based e-flow methods in the natural flow regime. The study used different hydrological-based methods, namely, the Global Environmental Flow Calculator, the Tennant method, the flow duration curve method, the dynamic method, the mean annual flow method, and the annual distribution method to allocate e-flows in the Kaligandaki River. The most common practice for setting e-flows consists of allocating a specific percentage of mean annual flow or portion of flow derived from specific percentiles of the flow duration curve. However, e-flow releases should mimic the river’s intra-annual variability to meet the specific ecological function at different river trophic levels and in different periods over a year covering biotas life stages. The suitability of the methods was analyzed using the Indicators of Hydrological Alterations and e-flows components. The annual distribution method and the 30%Q-D (30% of daily discharge) methods showed a low alteration at the five global indexes for each group of Indicators of Hydrological Alterations and e-flows components, which allowed us to conclude that these methods are superior to the other methods. Hence, the study results concluded that 30%Q-D and annual distribution methods are more suitable for the e-flows implementation to meet the riverine ecosystem’s annual dynamic demand to maintain the river’s health. This case study can be used as a guideline to allocate e-flows in the Kaligandaki River, particularly for small hydropower plants.

1. Introduction

1.1. Overview of E-Flow Concept

Climate change, rapid globalization, economic boost, and ever-increasing populations demand better management of water resources. To meet humankind’s demands (i.e., energy, water), many water conservancy and hydropower projects have been built [1]. Several hydropower plants are under construction or planned to be constructed, which will modify the rivers’ natural flow regime either through impoundments such as dams and barrages or through diversion work for agriculture or urban supplies [2]. These projects contribute to tackling the climate change effect (drought, floods) and improve the living standard of citizens (rent standards). However, its adverse effects, such as reducing total flow discharge left to the river, changing the seasonal flow regime [3,4], altering the magnitude and frequency of floods, and modifying the groundwater table are a severe threat to the aquatic ecosystem existence [5,6]. Most rivers with dams and diversion projects are at the ecological tipping point, which means act now or face the projects’ worst effects soon. Hence, along with humankind’s benefits, it is required to protect the rivers, natural lakes, and groundwater of the watershed by mitigating the riverine ecosystem [7].
The escalating hydrological alterations of the rivers flow regime, and its resulting severe impacts on the riverine ecosystem’s health is recognized globally [8,9,10]. With the advent of growing public consciousness in the river health and its hydrological alteration causing adverse impacts, river scientists developed the science of environmental flows (e-flows) assessments, which aid in determining the quality and quantity of water required for the protection of the riverine ecosystem and its inhabitants. The US clean water act 1977, 1992 European Commission (EC), Habitat Directive, and 2000 EC Water Framework Directives are examples that stressed mandatory river restoration to improve rivers’ ecological status. The Water Sustainability Act (2016) mandated that streamflow should not reduce below the environmental flow needed due to groundwater pumping in British Columbia and Canada. Environmental flow has been mandated in British Columbia, Canada, and California, USA [11]. An “Instream Flow Requirement (IFR),” “environmental flows,” “Environmental Flows Requirement (EFR),” or “environmental water demand (EWD)” are used as interchangeable terminology for the e-flows [2]. E-flows represent the flow regime provided within a river downstream to maintain a certain acceptable river health level and are widely used as environmental protection measures in many water conservancies projects [12,13]. E-flows is defined as “the quantity, timing, and quality of freshwater flows and levels necessary to sustain aquatic ecosystems which, in turn, support human cultures, economies, sustainable livelihoods, and well-being,” according to the Brisbane Declaration [14]. The process of defining the e-flows is known as “environmental flow assessment (EFA)” [15]. E-flows represent the water needed for the river ecosystem. However, there are other demands as well, such as residential, agriculture, and farming. Hence, all respective sectors’ demand should be considered to understand the relationship between water availability and water stress on the river ecosystem. For instance, Xu et al. [16] integrated the e-flow requirement in water-stress impact analysis to inform energy system deployment. This study did not consider any water needs except for hydropower production. All of the basin’s existing facilities’ water needs should be integrated to conduct water-stress analysis, giving us an idea about the basin’s energy system deployment. Groundwater plays a vital role in river hydrology and ecology at various scales (spatial and temporal), evaluating the spatial and temporal pumping effects on streamflow depletion while allocating environmental flows [17,18]. Gleeson and Richter [18] suggested that “high levels of ecological protection will be provided if groundwater pumping decreases monthly natural baseflow by less than 10% through time”. However, the discussed e-flows methodologies do not consider the impact of groundwater pumping. The study recommends integrating the impact of groundwater pumping and the impact of water infrastructure on environmental flows.
Environmental flow methodologies stem from a need to conserve mostly rivers and wetlands sustainably with the appropriate ecological balance in the water system close to the natural flow regime. This concern led to the development of different approaches in the environmental flows’ methodologies worldwide. Tharme [8] reviewed more than 200 methodologies spread across different geographical spectra. Though all of these methodologies, in principle, have the same goal to achieve a suitable environmental flow regime of water bodies, they differ in their working principle and assumptions made during derivation. These vast methodologies can be classified into four major categories: (1) hydrological, (2) hydraulic rating, (3) habitat simulation, and (4) holistic methodologies. The first two categories, hydrological and hydraulic rating, are based on the assumption that water bodies’ habitats/ecosystem functions (e.g., rivers) degrade with reducing water availability in the river. Whereas, habitat simulation assumes an “optimal” flow in the river, sustaining the river ecology in balance [9,19]. Holistic methods are more comprehensive methods with environmental flows designed to mimic the natural hydrograph [20]. Table 1 shows a short overview of all types of methods. The different factors, namely, types of the river (e.g., perennial, seasonal, high base flow, flashy); ecological importance of the river; stakeholders involved; cost; and difficulty of obtaining data decide the preference of the EFA methods [21]. For brevity, the concepts and details of these methods are described in the literature [19,20,22,23,24,25].
Hydrological methods use historical flow data records (e.g., daily, monthly, seasonal, yearly flow) to recommend an e-flows setting for maintaining the desired level of river health [24]. Rather than focusing on the optimized environment for single species, this method encompasses the conservation of rivers’ overall ecological integrity. The hydrological methods assume that there exists a relationship between flow parameters and biological attributes [27]. Based on this, different parameters are developed to optimize the flow regime. However, hydrological methods have some limitations. They can be used only in gauged catchments, are sensitive to hydrological data, and assume that all the aquatic organisms need the same quantity of water for survival. Due to its easiness of application, rapid assessment, low cost, and fewer field visits with most of the data readily arrived or can be simulated; hydrological methods are widely used to calculate environmental flow [28]. About 30% of all methods are hydrological-based [29]. Most of the preliminary and planning studies’ hydrological methods are used [30]. Suwal et al. [31] used the desktop hydrological method “Global Environmental Flow Calculator” (GEFC) to calculate the e-flow for different classes. Pastor et al. [32] recommended the suitability of the Variable Monthly Flow (VMF) and Tessman methods, especially for the variable flow regimes river.

1.2. Environmental Flows Practices in Nepal

The world is focusing and giving prime importance to the ecological impact assessment of water-related projects like hydropower and developing a robust e-flows setting methodology to ensure rivers’ ecological integrity [27,33,34]. The concept of e-flows is burgeoning on a global scale. However, in Nepal, the e-flows concept is still in its infancy as investors and the government are not serious about the rivers’ e-flows requirements while constructing water conservancy projects on the rivers in Nepal. Many hydropower projects (storage, run of the river, and peaking run of the river) are at the construction phase. Most of the projects are located on most of Nepal’s major rivers’, namely in Kaligandaki, Karnali, and Mahakali. Developed countries have used different advanced and robust methodological approaches for e-flows regulation. However, developing nations like Nepal usually use hydrological methods. Though hydropower development goes back a century, environmental flow concern was not significant until the Water Resource Act, 1992 in Nepal. This act gave the basic idea about the “environmental study” but did not encompass the broader environmental flows aspect. The Environmental Protection Act (EPA), 1997, was the turning point for e-flows in Nepal as it gave basic guidelines and highlighted the need for Initial Environmental Examination (IEE) and Environmental Impact Assessment (EIA) based on the installed capacity of projects. The EPA made the Environmental Protection Rules (EPR, 1997), making it mandatory to conduct an EIA for projects above 50 MW, and an IEE study for projects below 50 MW [35]. The introduction of the Hydropower Development Policy, 2001, became the paradigm shift in the e-flows setting in Nepal as it exclusively defined the minimum flow requirements for hydropower projects built onwards. Ecological, downstream, or environmental flows are the terms used in Nepal for the e-flows. The policy stated that “Provision should be made to release such quantum of water which is higher of either at least 10% of the minimum monthly average discharge of the river/stream or the minimum required quantum as identified in the environmental impact assessment study report” [35]. The working Policy for Construction and Operation of Physical Infrastructure within Protected Area (2009) further defined different provisions of e-flows that is: If the headworks is within the conservation areas, at least 50% of the monthly flow is considered as e-flows and if the headworks are not within conservation areas, but the downstream flow through the conservation areas, at least 10% of the monthly flow is considered as e-flows [35]. With the support of the International Water Management Institute (IWMI) and the World-Wide Fund for Nature (WWF), e-flows assessment in Nepal is undergoing in many rivers [36,37]. About 50% of all large rivers are affected by dam construction. Studies have identified that cold-water fish in Nepal are threatened due to block connectivity by dam blockages of hydropower projects [38]. Snow trout and gold mahseer are critically endangered fish species, and dark mahseer and Gangetic ailia are the vulnerable fish species according to the International Union for Conservation of Nature (IUCN).
Furthermore, the irrigation policy (2014) [39] gave the working direction and decided to use only the residual water for irrigation purposes after maintaining the minimum required water flow in the river and creeks. During the last two decades, many hydropower projects had been constructed. More projects are in the construction phase and the pipeline [40]. The EIA reports show the requirement of e-flows in the projects; however, due to a lack of monitoring resources, none of Nepal’s projects has been following the policy regarding the e-flows implementation, which is a severe threat to the downstream ecosystem of the projects. For instance, the Modi River did not get environmental flows as prescribed in EIA and IEE reports except for the wet season [41,42].
The case study on the environmental flow assessment discussed in this paper is located in Nepal’s Himalayan region. Table 1 shows that hydrological methods are best suited for Nepal on the primary investigation of e-flows. The main objectives of the present study are: (i) study different existing hydrological e-flow methods (EFMs) and allocate e-flows using six different e-flow methods; (ii) comparison of 6 EFMs; (iii) compute flow alteration using the Indicators of Hydrological Alterations (IHA) indicators and e-flows components (EFC) and application of global indexes; and (iv) to suggest better e-flows assessment methods (hydrological). The Kaligandaki River is considered as a case study for this investigation. The study is organized as follows: Section 2 presents details about the study area and applied methodologies. Section 3 shows the main results of the study. The discussion and conclusion are developed in Section 4 and Section 5, respectively.

2. Materials and Method

2.1. Study Area

Kaligandaki is one of Nepal’s major rivers and the Narayani River (Figure 1). Narayani River joins the Ganges River in India as a left-bank tributary, eventually draining at the Bay of Bengal.
This paper’s study station is Kotagaon Shringe hydrometric station located at 27°45′00″ N latitude and 84°20′50″ E longitude at an elevation of 198 m [6]. The station lies downstream of the powerhouses of the Kaligandaki-A hydropower station. The daily flow values recorded at the Kotagaon station were obtained from the Department of Hydrology and Meteorology (DHM), Nepal (http://dhm.gov.np/). The mean daily flows data from the year 1964–2015 were used in the study. The details of the Kaligandaki River are shown in Table 2. These flow values were processed in MS-Excel.

2.2. Methodology

In this study, six different hydrologic-based EFA methods were used to evaluate e-flows in the Kaligandaki River and later compared its effectiveness and influence on the natural flow regime using IHA indicators. The leading cause of using hydrological methods is the lack of ecological information of the basin, which is a must in other advanced methods, such as holistic habitat simulation methods. The figurative workflow of the study is shown in Figure 2. The key features of the used methods are discussed in the following sections.

2.2.1. Annual Distribution Method

The Annual Distribution Method (ADM) method is based on the mean ratio index, a ratio between annual average runoff and minimum annual average runoff. The ADM assumes that the minimum monthly average runoff of rivers could meet the essential water needs to maintain essential ecological, environmental function, give a favourable environment for the survival of the aquatic organism and riverine ecosystem will not suffer severe irreversible damage; however, it cannot reflect the hydrological characteristics of the river. Nevertheless, for many years, the average monthly discharge process can better reflect the overall historical river discharge process, such as timing, duration, frequency, and flow rate change. The mean ratio index is determined and quantified in the same period. According to the long series data of the natural average monthly flow of hydrological station, the average annual discharge Q ¯ and the minimum annual discharge Q ¯ min are calculated respectively [43]. The calculation formulas are given below:
Q ¯ = 1 12 i = 1 12 q ¯ i
q ¯ i = 1 n j = 1 12 q i j
Q ¯ min = 1 12 i = 1 12 q min ( i )
q min ( i ) = min ( q i j ) , j = 1 , 2 , , n
where q ¯ i is the average monthly discharge of ith month during series of n years (m3/s); q min ( i ) is the minimum monthly discharge of ith month for n years (m3/s); q i j is the average discharge of ith month jth year; n is the number of years data available.
Calculate the mean ratio index ( η ) using an annual average discharge ( Q ¯ ) and minimum annual average discharge ( Q ¯ min ) at the same time.
η = Q ¯ min Q ¯
Q i = q ¯ i × η
Calculate the basic ecological flow for each month ( Q i ) using Equation (6), where i is the month, i.e., i = 1, 2–12. The ADM method calculates monthly environmental flow; however, we needed daily flow for a year to compute IHA flow alteration. The study used Equation (7) to convert monthly environmental flow into daily flow [44].
Q e n v , j i = Q n , j i Q n i × Q e n v i
where Q e n v , j i is the minimum environmental flow on the jth day of the ith month, Q n , j i is the natural flow on the jth day of the ith month, Q n i is the average value of the natural flow of the ith month, and Q e n v i is the minimum environmental flow of the ith month.

2.2.2. Global Environmental Flow Calculator

The study used the “Global Environmental Flow Calculator” (GEFC) software [13] to calculate the environmental flow of different Environmental Management Classes (EMC). The GEFC software implements the “FDC Shifting” method. The details about the software are described in Smakhtin and Eriyagama [13]. The method describes ‘Environmental Management Classes’ (EMC). Table 3 shows the details about EMC and its corresponding ecological description with a management perspective. It classified EMC into six classes giving six similar environmental flow levels to each class. The software gives the percentage of mean annual flow (MAF) in each EMC. The study considers only EMC from B to F.

2.2.3. Flow Duration Curve Analysis

The Flow Duration Curve Analysis (FDCA) method is another popular method used in environmental flow calculations. A typical flow duration curve represents the proportion of flow exceeded for a particular time in the river section. Based on this exceedance curve, the minimum threshold is defined to preserve the ecological integrity of rivers. Generally, indices related to flow duration curves are developed. The maximum abstraction level of water from the river can be subsequently calculated. This method is useful in setting environmental flows downstream of hydropower. The equation used to compute the exceedance probability, which also is referred to as the flow-duration percentile, is given as:
P = 100 × ( m / ( n + 1 ) )
where P is the exceedance probability, m is the ranking, from highest to lowest, of all daily mean flows for the specified period of record, and n is the total number of daily means flows in the recorded period.
In this study, Q80%, Q85%, and Q90% are used as low flow indices.

2.2.4. Tennant Method

The Tennant method [45,46] suggested that specific percentages of the average annual flow (AAF) are necessary to maintain a river ecosystem’s biological integrity. It assumed that the aquatic organisms’ water requirement depends on different life cycles, such as the reproductive stage, global growth stage, etc. These stages have different water requirements. Hence, the method divided the whole year into a spawning period (April–September) and the general growth period (October–March). Table 4 shows the percentage standards of an aquatic organism’s water requirement in different life cycle stages.

2.2.5. Dynamic Methods

30% of Mean Daily Flow (30%Q-D).
This is based on the concept of minimum daily flow. It releases 30% of mean daily flow, considering a long series of interannual mean daily flow data, allowing for dynamic e-flows releases [47,48,49].

2.2.6. Mean Annual Flow

At least 10 or 25% of the Mean Annual Flow (MAF) must be released to the downstream depending upon the degree of environmental protection in the river reach. For satisfactory results, at least five years of continuous daily flow data are needed. For more representative results, a long time series of interannual mean daily flow data is required for this method [47].

2.2.7. Indicators of Hydrological Alteration (IHA) and Global Indexes

Richter et al. [50] used 32 “ecologically relevant” hydrological indexes to develop a set of Indicators of Hydrologic Alteration (IHA). The IHA indexes represent five essential parameters of the natural flow regime: magnitude, frequency, duration, timing, and rate of change [51]. The IHA indexes captured most of the variation and information described by 171 indexes, which makes IHA indexes the best choice to represent alteration on rivers [52]. The study used the Richter et al. [53] approach to assess the flow regime and e-flows component (EFC) alteration using relative mean deviation between the natural flow regime (NFR) and e-flows release from the 6 EFMs.
Kuriqi et al. [47] proposed a global index for a separate group of IHA indicators to simplify the analysis. Initially, the alteration of each indicator was computed using Equation (9), which is the relative mean difference between natural flow regime (NFR) and altered flow regimes (AFR) divided by nfr. Here, afr is the regime that will be obtained after releasing the e-flows instead of natural flows. After that, each group’s new global index was computed using Equation (10), each group’s average. The details about indexes are shown in Table 5.
H I i , j = | H I i , j ( n f r ) H I i , j ( a f r ) H I i , j ( n f r ) |
I j = H I i , j N
where H I i , j is the relative mean difference of i hydrological indicator of j group. I j : global alteration index of j group, H I i , j ( n f r ) : hydrological indicator related to nfr, H I i , j ( a f r ) : hydrological indicator related to the altered flow regime, and N ; the total number of IHA indexes for each group. Here, i = 1,2–12 depend upon the group and j = 1,2,3,4,5. The five indices are classified, as shown in Table 6.

2.2.8. Environmental Flow Components

The five essential components of flow, namely: low flows, extreme low flows, high flow pulses, small floods, and large floods, have been identified ecologically essential and have been incorporated in IHA software as “environmental flow components” [51]. Each of the flow components has a respective role in the life of an organism. For instance, extreme low flows reduce water connectivity, restricting organisms’ movement; high flow pulses help aquatic mobile organisms move upstream and downstream of rivers. Different species, different geographic conditions, and different rivers (perennial, ephemeral) could be life-threatening, causing the death of organisms or may provide favorable conditions for aquatic organisms for their life stages [51]. The study used the e-flows component (EFC) to assess the alteration using each parameter’s relative mean difference.

2.2.9. Limitations of the Methodology

The hydrological methods are the most simple, straightforward, and data-friendly methods that have been used extensively for the preliminary study of the e-flows allocation [48]. So, when there are several available solutions for the same problem, there is no guarantee of the best solution from both fixed and scientific perspectives. Single flow indices and other hydrological methods (such as Tennant, GEFC, FDC, ADM) have been applied globally, having many advantages. These methods simplify river basin planning work, needless time, and money. They require a low-level of knowledge related to the eco-hydrology of the basins [29,57]. The study’s analysis was limited to hydrological-based e-flows methods, which required only time-series data of mean daily flows [47]. Other sophisticated methods are also available such as habitat simulation and holistic approach, which are more directly involved in ecological concepts. They were not used in the analysis due to a lack of data in the present case study. However, the study can be continued further using those methods if the required data are available.

3. Results

3.1. E-Fows Allocation

Different hydrological methods are used to compute the potential e-flows in the river Table 7 presents the estimated e-flows from different methods. Daily flow data collected over 51 years (between 1964 and 2015) were used to calculate e-flows. The GEFC was used to calculate e-flows for different management classes.
This case study only considered Class B classes to Class F, as Class A flows could not be considered as e-flows for any project. The lowest class, ‘F’, shows 15.7% of MAF, that is, 70.44 m3/s, which indicates that below this amount would characterize the river as a dead environment. An average annual e-flows allocation of 147.16 m3/s (32.8% of MAF) is expected to maintain the essential ecosystem functions. Figure 3 and Figure 4 show the interannual discharge and monthly mean flow of all EFMs.
The study considers the Tennant method’s fair and degrading condition criteria, which allocates 10% of the MAF (equivalent to 44.87 m3/s) for the flow from October to March and 30% of the MAF (equivalent to 134.6 m3/s) for the flow from April to September. The method is simple and easy to implement, which lets it become more user friendly. It must be noted that 10% of the 6 months’ flow is the lowest among all the methods studied. Hence, it may create critical conditions for the river ecosystem during a 10% period.
The mean annual flow method that is 10% MAF allocated 44.97 m3/s to maintain the health of river reach, which is the second-lowest allocation showing its vulnerability of implementation. The dynamic methods of 30%Q-D methods follow the natural hydrograph pattern showing its applicability in the EFs allocation.
The FDC curve could be determined using whole multi-year data, or the FDC curve could be determined for each year separately. Młyński et al. [58] showed significant differences between the Qp (‘p’ percentage exceedance discharge) for the multi-year curve and Qp for mean or median annual curves. The study used multi-year curves to determine the e-flows of the Kaligandaki River. FDC or percentile methods such as Q80, Q85, and Q90 was considered for the study. This method allocated e-flows as the mean daily discharge that is equaled or exceeded by 80% (Q80), 85% (Q85), and 90% (Q90). The e-flows suggested by these methods are Q80= 49.04 m3/s, Q85 = 45.65 m3/s, and Q90 = 43.05 m3/s. The EF estimations from the FDCA methods considered here are less than Class F (equivalent to 70.44 m3/s), which critically modified the river ecosystem. Hence, it can be concluded that FDCA methods are not suitable for the EF estimations.
The ADM allocated the lowest e-flows at 45.07 m3/s for March and the highest for August at 785.59 m3/s, which is shown in Table 5. The method considered the intra-annual variation of the flow regimes, which can meet the actual need to sustain the ecological function of the river instead of considering the specific percentages of average annual runoff or specific guarantee rate of frequency curve of natural runoff as e-flows.

3.2. Interannual and Seasonal E-Flows Characterization

Figure 5 and Figure 6 show that the constructed FDC’s shape for different EF methods differs from NFR for mean annual flow and seasonal flow (autumn, spring, summer, and winter). The FDC curve of ADM and 30%Q-D methods show a slope, which means this method tries to maintain the rivers’ variability, but not by other methods. From the figure, we can see all methods on the FDC curve are below the mean NFR.
The FDC curve of ADM and 30%Q-D show a similar pattern as NFR but with lower values. Simultaneously, all methods are low flow fixed values methods, represented in a straight line rather than normal FDC curves. The annual mean FDC curves show that all methods value is less than low flows of the NFR, which gives us an idea about the insufficiency of the flow in the river. The seasonal FDC of all methods shows that autumn and summer FDC is the same as the mean annual FDC; however, FDC of spring and winter shows that Class B and Class C e-flows are higher than mean seasonal flows, which demonstrates suitability for the low flow seasons. In contrast, all other methods are below the mean seasonal flows. The ADM and 30%Q-D methods show a similar pattern as NFR FDC; however, they give very low e-flows during spring and winter seasons. Most of the e-flows methods give a straight line FDC, which means it does not consider the river’s flow variability.
The methods which mimic the shape of NFR will be the better methods to sustain the health of the river. In this way, we can choose the different methods for a different season or even different months to take each method’s strength and to give one robust method for e-flows calculation in different periods.

3.3. Flow Regime Alteration

3.3.1. IHA Alteration

Figure 7 shows the e-flows regime alteration degree using five indices developed by Kuriqi et al. [47] for all e-flow methods used in the study. The detailed information about the e-flows regime alteration regarding each e-flow method in Appendix A (Table A1). Each method has its strengths and weaknesses, allowing a certain percentage of flow into the rivers. Figure 7 shows the flow alteration due to different EFMs, which may have a high impact on the hydro-ecosystem and related hydro-ecological process.
For instance, 10% MAF, Q85%, and Q90% show a high Mean Monthly Flow Alteration Index (IMM). Whereas, 30%Q-D, Q80%, Class E, Class F, and Tennant show moderate IMM. While ADM, Class B, Class C, and Class D show mild IMM. For the second global index, Magnitude and Duration of Extreme Flow Alteration Index (IMDE), only Class B shows great alteration. ADM and Class D show mild alteration while all remaining methods show moderate alteration. Looking at the Timing of Extreme Flow Alteration Index (ITE), 10% MAF, Q80%, Q85%, Q90%, and Tennant show moderate alteration. Only Class B shows mild alteration, but all remaining EFMs methods show low alteration. Frequency and Duration Alteration Index (IFD) shows that environmental management classes Class C, Class D, and Class E show moderate. In contrast, 30%Q-D, ADM, and Class B show low alteration, but other EFMs, for instance, 10%MAF, Q80%, Q85%, Q90%, Class F, and Tennant show high alteration. The Rate and Frequency Alteration Index (IRF), all except ADM and 30%Q-D, show high alteration while ADM and 30% Q-D show low alteration showing these methods effectiveness in e-flows determination.

3.3.2. E-flows Components (EFC)

The study calculated relative changes in the mean of e-flows component (EFC) regarding six EFMs against NFR using the values obtained from IHA software. The results are shown in Table 6. The results show that most of the EFC low flow changes 100% negatively except for 30% Q-D, ADM, and Tennant methods. For 30% Q-D, all EFC Low Flows changes −70%, while in ADM, it varies between −31% during June and 65% during August. The Tennant method varies between −6% during April and −100% for the winter season, extending up to March month. Table 6 shows the relative mean change of the EFC parameters against NFR. Most of the methods, namely 10% MAF, Q80%, Q85%, Q90%, Class B to Class F, show all EFC parameters changes of −100%. While 30%Q-D method shows duration, timing, and frequency of extreme low, high flow, small flood, and enormous flood vary by 0%, and remaining parameters show −70% relative changes against NFR. The ADM method shows a variation of relative mean change within a range between −2% of high flow frequency and +260% of significant flood frequency.
In comparison, the Tennant method shows a variation between −39% extremely low peak and 1393% extremely low duration. The variability of flow is high when the positive and relative change of Cv is high. In contrast, flow variability is low when Cv’s negative and relative change is of low value.

4. Discussion

The results show that MAF, Tennant (Oct–Mar), and FDC methods allocate e-flows less than the GEFC class F (critically modified) method; hence these methods are not recommended for the e-flows assessment for the present case-study. Further, Class B (Slightly modified) and Class C (moderately modified) allocated a considerable amount of e-flows; however, they gave a fixed value which is a certain percentage of MAF. Nevertheless, river flow is dynamic with interannual variability. Hence, ADM and dynamic method (30%Q-D) are recommended for e-flows allocation. These methods consider the changing characteristics of natural runoff. Hence, to ensure the sustainability of the river health, we must choose the methods that maintain the river’s flow dynamics rather than an absolute fixed percentage of MAF of the river.
NFR of many rivers worldwide had been altered by anthropogenic activities such as diversions and impoundment work [5,6,44]. For the rivers with limited data, time, and funding, the hydrological methods are the best method to allocate e-flows. The applied five alteration indexes showed that there is a considerable difference in alteration regarding each EFM. The degree of alteration was varied from EFM to EFM. The EFM, such as ADM and 30%Q-D, appeared to be less altered. This may be due to the consideration of the flow variability of the river [48]. Here, the ITE, IFD, and IRF, which have an essential role in sustaining the different ecological processes [54,59], were preserved near NFR, which means they are less altered only for two EFMs, that is, ADM and 30%Q-D. Overall, the results obtained by taking an average of five global indexes showed that the 10% MAF, Q80%, 85%, 90%, and Tennant methods gave high alteration, class D, E, and F showed moderate alterations, indicating their unsuitability in the e-flows allocation. The remaining methods showed mild alterations, indicating the suitability of e-flows allocation. The ADM and 30%Q-D methods showed the lowest and second-lowest alteration among all considered methods with values of 0.380 and 0.406, respectively. This is because the ADM and dynamic methods were the only methods that considered the concept of a dynamic pattern of the river like the NFR.
All six EFMs showed a high alteration of an e-flows component (EFC). Among those methods, except ADM, 30%Q-D, and the Tennant method, all others showed a −100% alteration, confirming the methods’ applicability. ADM gave a lower EFC alteration; it showed a monthly low flow alteration between −31% to −65%, resulting from less water in the river than the naturally available water in the riverine ecosystem. The low alteration of monthly flows means it maintains the temperature, flow velocity, and connectivity needed for most of the aquatic habitat than other methods [51]. Extreme low flows, high flow pulses, small floods, and large floods low alterations further give strong supports to the suitability of the method, which is shown in Table 8, because these components play a crucial role in maintaining the health of the rivers [51].
The FDC is a plot of the observed historical river flow data collected with different temporal resolutions daily, weekly, monthly, seasonally, and annually at the gauge station and when that flow is equaled or exceeded [60,61,62]. The FDC is an informative method that shows the characteristics of flow. The slope of the flow-duration curve is a quantitative measure of the flow regime variability [62]. The NFR FDC generally showed low flows exceeded most of the time while high flows are exceeded infrequently. However, the paper’s studied e-flows methods showed low flow exceedance 100% of the time, except for the ADM and 30%Q-D. The FDC curve plotted for annual mean and mean seasonal flows showed a drastic reduction in river flows, especially the high flood. However, the ADM and 30%Q-D methods try to mimic the FDC curve of NFR, which might help sustain the riverine ecosystems’ health. Other methods that release fixed minimum flow rather than dynamic release gives a straight line FDC, removing high floods, small floods and it is lower than NFR over the year. The FDC is useful to evaluate the relationship between magnitude and frequency of the river flow; however, it does not maintain temporal sequences of flows and so is unsuccessful in meeting the criteria of the timing or duration of the e-flows [63]. Nevertheless, as recommend by Kuriqi et al. [9,49], to guarantee suitable habitat conditions during low flow periods, AMD and 30%Q-D should be combined with other methods by setting a minimum of e-flows to be released downstream of the water intake. The results of this case study are aligned with the recommendations of the investigations above.

5. Conclusions

In developing countries like Nepal, where ecological information is insufficient and has inadequate baseline data regarding e-flows assessment, the hydrological methods can estimate e-flows requirements for planning and study phases. This study aimed to discuss e-flows calculation methodologies (hydrological) and to discuss the present status of the e-flows in Nepal. We compared six hydrological-based EFMs to allocate e-flows, to evaluate flow alteration, to estimate relative change (%) of the EFC against NFR, and to characterize the interannual and seasonal e-flows of the Kaligandaki River.
The results of the study showed that the global indexes such as Frequency and Duration Alteration Index (IFD) and the Rate and Frequency Alteration Index (IRF) showed a high alteration for all methods, except for the ADM and dynamic method (30%Q-D), which in turn showed a low alteration. The remaining three indexes, namely the Mean Monthly Flow Alteration Index (IMM), the Magnitude and Duration of Extreme Flow Alteration Index (IMDE), and the Timing of Extreme Flow Alteration Index (ITE) showed moderate and mild alteration for all hydrological-based EFMs investigated in this case study. In the overall analysis, it can be seen that the flow alteration of five indexes and the e-flows component (EFC) is lower for ADM and dynamic methods compared to other hydrological methods considered in this study. Furthermore, the FDC of annual mean flows and annual seasonal mean flows showed a dramatic decrease in river flows, especially the high flows in most e-flow methods except for ADM and 30%Q-D methods. This concludes the practicability of the ADM and dynamic methods; it reflects the interannual variability of the river to meet the specific ecological function of the different sections of the river in different periods. Nevertheless, we suggest that the application of those methods should be made under biota requirements at a given river.
The ADM method used in the study was specially designed for large and medium-size perennial rivers, not for temporary or seasonal rivers. The runoff process of these rivers may be intervened more often, which can create a large error in the calculation procedure. Furthermore, many researchers highlighted that the interannual and intra-annual variability of flow in the river must be maintained to sustain the ecological biodiversity. Hence, within hydrological methods, a method that considers the dynamic nature of the flow regime of the river is recommended for the e-flows allocation.

Author Contributions

Conceptualization, methodology, software, N.S. and A.K.; formal analysis, validation, data curation, writing—original draft preparation, N.S. and A.K.; writing—review and editing, visualization, N.S., A.K., X.H., J.D., D.M., and A.W.; funding acquisition and supervision, A.K. and A.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The first author would like to thank Kamal Prasad Pandey for his invaluable information on Nepal’s hydrology and data acquisition. The first author is highly indebted to Khwopa College of Engineering for rendering a warm environment for the research work. Alban Kuriqi was supported by a Ph.D. scholarship granted by Fundação para a Ciência e a Tecnologia, I.P. (FCT), Portugal, under the Ph.D. Program FLUVIO–River Restoration and Management, grant number PD/BD/114558/2016. Thanks to the anonymous reviewers for their comments and suggestions to further improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. IHA indicators are showing flow alteration regarding all applied EFMs against NFR.
Table A1. IHA indicators are showing flow alteration regarding all applied EFMs against NFR.
IHA ParametersMean
30% Q-D10% MAFQ80%Q85%Q90%ADMClass BClass CClass DClass EClass FTennant
Group #1
July0.700.960.960.960.960.460.820.870.910.930.940.89
August0.700.970.970.970.970.460.850.900.920.940.950.91
September0.700.950.950.950.960.540.790.850.890.920.930.87
October0.700.900.890.900.910.500.540.680.770.820.850.79
November0.700.800.780.800.810.500.120.350.530.630.690.79
December0.700.700.680.700.720.500.030.100.300.450.540.70
January0.700.620.580.610.630.500.050.060.140.290.400.62
February0.700.550.510.540.570.520.050.060.080.180.290.55
March0.700.500.450.490.520.490.040.050.060.140.230.50
April0.700.540.500.540.560.460.040.060.100.200.300.35
May0.700.690.660.680.700.340.050.140.290.420.510.19
June0.700.890.880.880.890.430.510.640.730.790.820.67
Group #2
1-day minimum0.701.100.090.070.060.020.690.000.000.020.060.32
3-day minimum0.700.390.340.380.420.421.001.000.000.030.100.39
7-day minimum0.700.410.350.400.430.421.000.000.000.040.110.41
30-day minimum0.700.460.410.450.480.461.000.020.030.090.180.46
90-day minimum0.700.510.470.500.530.471.000.030.050.140.240.51
1-day maximum0.700.990.990.990.990.641.000.960.970.980.980.97
3-day maximum0.700.990.980.980.990.590.930.950.960.970.980.96
7-day maximum0.700.980.980.980.980.530.910.940.950.960.970.94
30-day maximum0.700.970.970.970.970.510.870.910.940.950.960.92
90-day maximum0.700.960.960.960.960.490.830.880.910.930.940.89
Number of zero days0.000.000.000.000.000.000.000.000.000.000.000.00
Base flow index04.5803570.00080.00080100.80711.4556770.2556470.2531940.1848670.0801630.452621
Group #3
Date of minimum0.0000.940.900.940.940.050.000.0000.0000.1480.4751.000
Date of maximum0.000.180.180.180.180.030.180.180.1850.1780.1780.183
Group #4
Low pulse count010.988495110.70969200.0674540.2726590.5806160.7360161
Low pulse duration0.001.000.911.001.000.410.000.090.560.610.751.00
High pulse count0.001.001.001.001.000.291.001.001.001.001.000.44
High pulse duration0.001.001.001.001.000.461.001.001.001.001.004.56
Low Pulse Threshold0.700.570.530.560.590.450.010.050.140.250.341.00
High Pulse Threshold0.700.960.950.960.960.520.790.850.890.920.930.88
Group #5
Rise rate0.701.000.991.001.000.670.910.930.960.970.980.86
Fall rate0.701.000.991.001.000.630.900.920.940.950.970.68
Number of reversals0.001.001.001.001.000.080.440.530.690.860.940.87

References

  1. Couto, T.B.; Olden, J.D. Global Proliferation of Small Hydropower Plants—Science and Policy. Front. Ecol. Environ. 2018, 16, 91–100. [Google Scholar] [CrossRef]
  2. Karimi, S.S.; Yasi, M.; Eslamian, S. Use of Hydrological Methods for Assessment of Environmental Flow in a River Reach. Int. J. Environ. Sci. Technol. 2012, 9, 549–558. [Google Scholar] [CrossRef] [Green Version]
  3. Kuriqi, A.; Ali, R.; Pham, Q.B.; Gambini, J.M.; Gupta, V.; Malik, A.; Linh, N.T.T.; Joshi, Y.; Anh, D.T.; Nam, V.T.; et al. Seasonality Shift and Streamflow Flow Variability Trends in Central India. Acta Geophys. 2020, 68, 1461–1475. [Google Scholar] [CrossRef]
  4. Ali, R.; Kuriqi, A.; Abubaker, S.; Kisi, O. Long-Term Trends and Seasonality Detection of the Observed Flow in Yangtze River Using Mann-Kendall and Sen’s Innovative Trend Method. Water 2019, 11, 1855. [Google Scholar] [CrossRef] [Green Version]
  5. Huang, X.; Suwal, N.; Fan, J.; Pandey, K.P.; Jia, Y. Hydrological Alteration Assessment by Histogram Comparison Approach: A Case Study of Erdu River Basin, China. J. Coast. Res. 2019, 93, 139–145. [Google Scholar] [CrossRef]
  6. Gao, Y.; Pandey, K.P.; Huang, X.; Suwal, N.; Bhattarai, K.P. Estimation of Hydrologic Alteration in Kaligandaki River Using Representative Hydrologic Indices. Water 2019, 11, 688. [Google Scholar]
  7. Suwal, N.; Huang, X.; Pandey, K.P.; Bhattarai, K.P. Assessment of Hydrological Alteration and Selection of Representative Hydrological Indicators in Erdu River. In Proceedings of the ICWRER 2019, Nanjing, China, 14–18 June 2019. [Google Scholar]
  8. Tharme, R.E. A Global Perspective on Environmental Flow Assessment: Emerging Trends in the Development and Application of Environmental Flow Methodologies for Rivers. River Res. Appl. 2003, 19, 397–441. [Google Scholar] [CrossRef]
  9. Kuriqi, A.; Pinheiro, A.N.; Sordo-Ward, A.; Garrote, L. Water-Energy-Ecosystem Nexus: Balancing Competing Interests at a Run-of-River Hydropower Plant Coupling a hydrologic–ecohydraulic Approach. Energy Convers. Manag. 2020, 223, 113267. [Google Scholar] [CrossRef]
  10. Ali, R.; Kuriqi, A.; Abubaker, S.; Kisi, O. Hydrologic Alteration at the Upper and Middle Part of the Yangtze River, China: Towards Sustainable Water Resource Management Under Increasing Water Exploitation. Sustainability 2019, 11, 5176. [Google Scholar] [CrossRef] [Green Version]
  11. Li, Q.; Gleeson, T.; Zipper, S.C.; Kerr, B. Too Many Streams and Not Enough Time or Money? New Analytical Depletion Functions for Rapid and Accurate Streamflow Depletion Estimates. OSF Preprints 2020. Available online: https://osf.io/gfhym (accessed on 21 October 2020).
  12. Dyson, M.; Bergkamp, G.; Scanlon, J. Flow: The Essentials of Environmental Flows; IUCN: Gland, Switzerland; Cambridge, UK, 2003; pp. 20–87. [Google Scholar]
  13. Smakhtin, V.; Eriyagama, N. Developing a Software Package for Global Desktop Assessment of Environmental Flows. Environ. Model. Softw. 2008, 23, 1396–1406. [Google Scholar] [CrossRef]
  14. Arthington, A.H.; Bhaduri, A.; Bunn, S.E.; Jackson, S.E.; Tharme, R.E.; Tickner, D.; Young, B.; Acreman, M.; Baker, N.; Capon, S.; et al. The Brisbane Declaration and Global Action Agenda on Environmental Flows. Front. Environ. Sci. 2018, 6, 6. [Google Scholar] [CrossRef] [Green Version]
  15. Pittock, J.; Lankford, B.A. Environmental Water Requirements: Demand Management in an Era of Water Scarcity. J. Integr. Environ. Sci. 2010, 7, 75–93. [Google Scholar] [CrossRef]
  16. Xu, H.; Lee, U.; Coleman, A.M.; Wigmosta, M.S.; Sun, N.; Hawkins, T.R.; Wang, M.Q. Balancing Water Sustainability and Productivity Objectives in Microalgae Cultivation: Siting Open Ponds by Considering Seasonal Water-Stress Impact Using AWARE-US. Environ. Sci. Technol. 2020, 54, 2091–2102. [Google Scholar] [CrossRef] [PubMed]
  17. De Graaf, I.E.M.; Gleeson, T.; Van Beek, L.P.H.R.; Sutanudjaja, E.H.; Bierkens, M.F.P. Environmental Flow Limits to Global Groundwater Pumping. Nat. Cell Biol. 2019, 574, 90–94. [Google Scholar] [CrossRef]
  18. Gleeson, T.; Richter, B. How Much Groundwater Can We Pump and Protect Environmental Flows through Time? Presumptive Standards for Conjunctive Management of Aquifers and Rivers. River Res. Appl. 2017, 34, 83–92. [Google Scholar] [CrossRef]
  19. Jowett, I.G. Instream Flow Methods: A Comparison of Approaches. Regul. Rivers Res. Manag. 1997, 13, 115–127. [Google Scholar] [CrossRef]
  20. Williams, J.G.; Moyle, P.B.; Webb, J.A.; Kondolf, G.M. Environmental Flow Assessment: Methods and Applications; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  21. Lumbroso, D.M.; Sakamoto, D.; Johnstone, W.M.; Tagg, A.F.; Lence, B.J. Development of a Life Safety Model to Estimate the Risk Posed to People by Dam Failures and Floods. Dams Reserv. 2011, 21, 31–43. [Google Scholar] [CrossRef] [Green Version]
  22. Acreman, M.C.; Dunbar, M.J. Defining Environmental River Flow Requirements—A Review. Hydrol. Earth Syst. Sci. 2004, 8, 861–876. [Google Scholar] [CrossRef]
  23. Shokoohi, A.; Hong, Y. Using Hydrologic and Hydraulically Derived Geometric Parameters of Perennial Rivers to Determine Minimum Water Requirements of Ecological Habitats (case Study: Mazandaran Sea Basin-Iran). Hydrol. Process. 2011, 25, 3490–3498. [Google Scholar] [CrossRef]
  24. Fuladipanah, M.; Jorabloo, M. Hydrological Method to Evaluate Environmental Flow (case Study: Gharasou River, Ardabil). Int. J. Environ. Ecol. Eng. 2015, 9, 62–65. [Google Scholar]
  25. Dubey, A.; Singh, O.; Shekhar, S.; Pohshna, C. Assessment of Environmental Flow Requirement Using Environmental Management Classes-Flow Duration Curve for Narmada River. Int. J. Curr. Microbiol. Appl. Sci. 2019, 8, 891–897. [Google Scholar] [CrossRef]
  26. Pandey, K.P. Study on Hydrologic Alteration and Alteration Parameter Reduction Methods. Master’s Dissertation, Hohai University, Nanjing, China, 2019. [Google Scholar]
  27. Smakhtin, V.U.; Shilpakar, R.L.; Hughes, D.A. Hydrology-Based Assessment of Environmental Flows: An Example from Nepal. Hydrol. Sci. J. 2006, 51, 207–222. [Google Scholar] [CrossRef] [Green Version]
  28. Suwal, N. Research on Optimal Operation of Cascade Hydropower Stations Considering Ecological Flows. Master’s Dissertation, Hohai University, Nanjing, China, 2019. [Google Scholar]
  29. Młyński, D.; Operacz, A.; Walega, A. Sensitivity of Methods for Calculating Environmental Flows Based on Hydrological Characteristics of Watercourses Regarding the Hydropower Potential of Rivers. J. Clean. Prod. 2020, 250, 119527. [Google Scholar] [CrossRef]
  30. Operacz, A.; Wałęga, A.; Cupak, A.; Tomaszewska, B. The Comparison of Environmental Flow Assessment—The Barrier for Investment in Poland or River Protection? J. Clean. Prod. 2018, 193, 575–592. [Google Scholar] [CrossRef]
  31. Suwal, N.; Huang, X.; Kuriqi, A.; Chen, Y.; Pandey, K.P.; Bhattarai, K.P. Optimisation of Cascade Reservoir Operation Considering Environmental Flows for Different Environmental Management Classes. Renew. Energy 2020, 158, 453–464. [Google Scholar] [CrossRef]
  32. Pastor, A.V.; Ludwig, F.; Biemans, H.; Hoff, H.; Kabat, P. Accounting for Environmental Flow Requirements in Global Water Assessments. Hydrol. Earth Syst. Sci. 2014, 18, 5041–5059. [Google Scholar] [CrossRef] [Green Version]
  33. Smakhtin, V.; Anputhas, M. An Assessment of Environmental Flow Requirements of Indian River Basins; IWMI: Colombo, Sri Lanka, 2006; Volume 107. [Google Scholar]
  34. Poff, N.L.; Tharme, R.E.; Arthington, A.H. Evolution of Environmental Flows Assessment Science, Principles, and Methodologies. In Water for the Environment; Academia Press: Cambridge, MA, USA, 2017; pp. 203–236. [Google Scholar]
  35. Gaudel, P. Environmental Assessment of Hydropower Development in Nepal: Current Practices and Emerging Challenges. Vidyut. Febraury 2015. Available online: https://www.researchgate.net/publication/316080737_Environmental_Assessment_of_Hydropower_Development_in_Nepal_Current_Practices_and_Emerging_Challenges (accessed on 21 October 2020).
  36. Doody, T.; Cuddy, S.; Bhatta. Connecting Flow and Ecology in Nepal: Current State of Knowledge for the Koshi Basin; Sustainable Development Investment Portfolio (SDIP) Project; CSIRO: Canberra, Australia, 2016. [Google Scholar]
  37. Oglethorpe, J.; Regmi, S.; Bartlett, R.; Dongol, B.S.; Wikramanayake, E.; Freeman, S.J.O. The Value of a River Basin Approach in Climate Adaptation. In Proceedings of the International Conference on Climate Change Innovation and Resilience for Sustainable Livelihoods, Kathmandu, Nepal, 12–14 January 2015. [Google Scholar]
  38. Gubhaju, S.R. Impact of Damming on the Aquatic Fauna in Nepalese Rivers. In Cold Water Fisheries in Thetrans-Himalayan Countries; Petr, T., Swar, S.B., Eds.; Food and Agriculture Organization of the United Nations: Rome, Italy, 2002; pp. 129–145. [Google Scholar]
  39. Panta, S.K.; Resurrección, B.P. Gender and Caste Relations Amidst a Changing Political Situation in Nepal: Insights from a Farmer-Managed Irrigation System. Gender Technol. Dev. 2014, 18, 219–247. [Google Scholar] [CrossRef]
  40. International Hydropower, A. Hydropower Status Report: Sector Trends and Insights; IHA: London, UK, 2018. [Google Scholar]
  41. Jalsrot Vikas Sanstha; GWP Nepal. Assessment of the Environmental Flow in the Gandaki River Basin: A Case of Modi Khola; Jalsrot Vikas Sanstha: Kathmandu, Nepal; GWP Nepal: Hetauda, Nepal, 2016. [Google Scholar]
  42. Rijal, N.; Shrestha, H.K.; Bruins, B. Environmental Flow Assessment of Hewa Khola A and Lower Hewa Khola Hydropower Projects in Nepal. Hydro Nepal J. Water Energy Environ. 2018, 23, 71–78. [Google Scholar] [CrossRef]
  43. Zha-Rong, P.A.N.; Xiao-Hong, R.U.A.N.; Jing, X.U. A New Calculation Method of Instream Basic Ecological Water Demand. J. Hydraul. Eng. 2013, 44, 119–126. [Google Scholar]
  44. Zhang, H.; Chang, J.; Gao, C.; Wu, H.; Wang, Y.; Lei, K.; Long, R.; Zhang, L. Cascade Hydropower Plants Operation Considering Comprehensive Ecological Water Demands. Energy Convers. Manag. 2019, 180, 119–133. [Google Scholar] [CrossRef]
  45. Tennant, D.L. Instream Flow Regimens for Fish, Wildlife, Recreation and Related Environmental Resources. Fish 1976, 1, 6–10. [Google Scholar] [CrossRef]
  46. Wałęga, A.; Mlynski, D.; Kokoszka, R.; Miernik, W. Possibilities of Applying Hydrological Methods for Determining Environmental Flows in Select Catchments of the Upper Dunajec Basin. Pol. J. Environ. Stud. 2015, 24, 2663–2676. [Google Scholar] [CrossRef]
  47. Kuriqi, A.; Pinheiro, A.N.; Sordo-Ward, A.; Garrote, L. Influence of Hydrologically Based Environmental Flow Methods on Flow Alteration and Energy Production in a Run-of-River Hydropower Plant. J. Clean. Prod. 2019, 232, 1028–1042. [Google Scholar] [CrossRef]
  48. Kuriqi, A.; Pinheiro, A.N.; Sordo-Ward, A.; Garrote, L. Flow Regime Aspects in Determining Environmental Flows and Maximizing Energy Production at Run-of-River Hydropower Plants. Appl. Energy 2019, 256, 113980. [Google Scholar] [CrossRef]
  49. Bejarano, M.; Sordo-Ward, A.; Gabriel-Martin, I.; Garrote, L. Tradeoff Between Economic and Environmental Costs and Benefits of Hydropower Production at Run-of-River-Diversion Schemes under Different Environmental Flows Scenarios. J. Hydrol. 2019, 572, 790–804. [Google Scholar] [CrossRef]
  50. Richter, B.D.; Baumgartner, J.V.; Powell, J.; Braun, D.P. A Method for Assessing Hydrologic Alteration Within Ecosystems. Conserv. Biol. 1996, 10, 1163–1174. [Google Scholar] [CrossRef] [Green Version]
  51. Mathews, R.; Richter, B.D. Application of the Indicators of Hydrologic Alteration Software in Environmental Flow Setting1. JAWRA J. Am. Water Resour. Assoc. 2007, 43, 1400–1413. [Google Scholar] [CrossRef]
  52. Olden, J.D.; Poff, N.L. Redundancy and the Choice of Hydrologic Indices for Characterizing Streamflow Regimes. River Res. Appl. 2003, 19, 101–121. [Google Scholar] [CrossRef]
  53. Fausch, K.D.; Bestgen, K.R. Ecology of Fishes Indigenous to the Central and Southwestern Great Plains. In Ecological Studies; Springer: New York, NY, USA, 1997; pp. 131–166. [Google Scholar]
  54. Rood, S.B.; Mahoney, J.M. River Damming and Riparian Cottonwoods along the Marias River, Montana. Rivers 1995, 5, 195–207. [Google Scholar]
  55. Richter, B.; Baumgartner, J.; Wigington, R.; Braun, D. How Much Water Does a River Need? Freshw. Biol. 1997, 37, 231–249. [Google Scholar] [CrossRef] [Green Version]
  56. Richter, B.D.; Baumgartner, J.V.; Braun, D.P.; Powell, J. A Spatial Assessment of Hydrologic Alteration within a River Network. Regul. Rivers Res. Manag. 1998, 14, 329–340. [Google Scholar] [CrossRef]
  57. Książek, L.; Woś, A.; Florek, J.; Wyrębek, M.; Młyński, D.; Wałęga, A. Combined Use of the Hydraulic and Hydrological Methods to Calculate the Environmental Flow: Wisloka River, Poland: Case Study. Environ. Monit. Assess. 2019, 191, 254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Młyński, D.; Wałęga, A.; Ozga-Zielinski, B.; Ciupak, M.; Petroselli, A. New Approach for Determining the Quantiles of Maximum Annual Flows in Ungauged Catchments Using the EBA4SUB Model. J. Hydrol. 2020, 589, 125198. [Google Scholar] [CrossRef]
  59. Cushman, R.M. Review of Ecological Effects of Rapidly Varying Flows Downstream from Hydroelectric Facilities. North Am. J. Fish. Manag. 1985, 5, 330–339. [Google Scholar] [CrossRef]
  60. Verma, R.K.; Murthy, S.; Verma, S.; Mishra, S.K. Design Flow Duration Curves for Environmental Flows Estimation in Damodar River Basin, India. Appl. Water Sci. 2016, 7, 1283–1293. [Google Scholar] [CrossRef] [Green Version]
  61. Vogel, R.M.; Fennessey, N.M. Flow-Duration Curves. I: New Interpretation and Confidence Intervals. J. Water Resour. Plan. Manag. 1994, 120, 485–504. [Google Scholar] [CrossRef]
  62. Searcy, J.K. Flow-Duration Curves; manual of hydrology. Part 2. US Geol. Survey Water Supply Paper 1542-A., Low flow techniques; United States Government Printing Office: Washington, DC, USA, 1959. [Google Scholar]
  63. Jain, S.K.; Kumar, P. Environmental Flows in India: Towards Sustainable Water Management. Hydrol. Sci. J. 2014, 59, 751–769. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Hydropower and hydrometric station locations within the Kaligandaki River Basin, Nepal. The upper figure shows the mean daily time series of flow discharge from the year 1964–2015.
Figure 1. Hydropower and hydrometric station locations within the Kaligandaki River Basin, Nepal. The upper figure shows the mean daily time series of flow discharge from the year 1964–2015.
Sustainability 12 08766 g001
Figure 2. The flow chart shows the main steps and e-flow methods applied in this study.
Figure 2. The flow chart shows the main steps and e-flow methods applied in this study.
Sustainability 12 08766 g002
Figure 3. Interannual discharge regarding all applied EFMs. The dashed line connects the mean values of e-flows discharge for each EFMs.
Figure 3. Interannual discharge regarding all applied EFMs. The dashed line connects the mean values of e-flows discharge for each EFMs.
Sustainability 12 08766 g003
Figure 4. Mean monthly e-flows regarding all EFMs. The year is divided into two parts; the first part refers to high flow seasons (July–December). The second part refers to the low flow seasons (January–June). The reading of the left-vertical axis corresponds to the first part. The reading of the right-vertical axis corresponds to the second part of the year.
Figure 4. Mean monthly e-flows regarding all EFMs. The year is divided into two parts; the first part refers to high flow seasons (July–December). The second part refers to the low flow seasons (January–June). The reading of the left-vertical axis corresponds to the first part. The reading of the right-vertical axis corresponds to the second part of the year.
Sustainability 12 08766 g004
Figure 5. Flow duration curve of interannual mean daily e-flows from 1964 to 2015, regarding of all applied EFMs.
Figure 5. Flow duration curve of interannual mean daily e-flows from 1964 to 2015, regarding of all applied EFMs.
Sustainability 12 08766 g005
Figure 6. Flow duration curve of seasonal mean daily e-flows regarding of all applied EFMs. Summer (June, July, August), autumn (September, October, November), winter (December, January, February), and spring (March, April, May).
Figure 6. Flow duration curve of seasonal mean daily e-flows regarding of all applied EFMs. Summer (June, July, August), autumn (September, October, November), winter (December, January, February), and spring (March, April, May).
Sustainability 12 08766 g006
Figure 7. Global indices derived from indicators of hydrological alteration regarding all applied EFMs. Namely, Mean Monthly Flow Alteration Index (IMM), Magnitude and Duration of Extreme Flow Alteration Index (IMDE), Timing of Extreme Flow Alteration Index (ITE), Frequency and Duration Alteration Index (IFD), and Rate and Frequency Alteration Index (IRF). The five indexes are classified into four categories, low (0–0.25), mild (0.25–0.5), moderate (0.5–0.75), and high alteration (0.75–1).
Figure 7. Global indices derived from indicators of hydrological alteration regarding all applied EFMs. Namely, Mean Monthly Flow Alteration Index (IMM), Magnitude and Duration of Extreme Flow Alteration Index (IMDE), Timing of Extreme Flow Alteration Index (ITE), Frequency and Duration Alteration Index (IFD), and Rate and Frequency Alteration Index (IRF). The five indexes are classified into four categories, low (0–0.25), mild (0.25–0.5), moderate (0.5–0.75), and high alteration (0.75–1).
Sustainability 12 08766 g007
Table 1. Selection criteria for different environmental flow calculation methodologies [26].
Table 1. Selection criteria for different environmental flow calculation methodologies [26].
Method CategoryResolution LevelEcosystemTimeCost
HydrologicVery Low/LowRiverShortLess
Hydraulic ratingLowRiverShort/LongLess/Medium
Habit simulationMedium/HighRiverMedium/LongMedium/High
HolisticHighWetland, floodplains, LongHigh
Table 2. Flow data and watershed description of Kaligandaki River, Nepal.
Table 2. Flow data and watershed description of Kaligandaki River, Nepal.
NameDetails
Elevation190 m to 8168 m
Total catchment area11.851 km2
Location82°52.8′ E to 84°26.3′ E, 27°43.2′ N to 29°19.8′ N
Mean annual precipitation1396 mm
Flow data Series1 January 1964–31 December 2015
Min flow (m3/s)46
Mean flow (m3/s)449.7
Max flow (m3/s)6840
Min average monthly flow (m3/s)
10% of min average monthly flow (m3/s)
90
9
Table 3. Description of Environmental Management Classes (EMC) used in the Global Environmental Flow Calculator [2,13].
Table 3. Description of Environmental Management Classes (EMC) used in the Global Environmental Flow Calculator [2,13].
EMCMost likely Ecological ConditionManagement Perspective
A
(Natural)
Same as natural rivers with insignificant modification of instream and riparian habitatProtected rivers and basins. Reserves and national parks. No new water projects (dams, diversions, etc.) allowed.
B
(Slightly modified)
Largely intact biodiversity and habitats despite anthropogenic activities (dam, diversion, basin modifications)Water supply schemes or irrigation development present and/or allowed.
C
(Moderately modified)
The biota’s habitats and movement have been impacted, but essential ecosystem functions are still unmodified; some sensitive species are vanished and/or reduced in extent; alien species survived.Multiple disturbances (for instance, dams, diversions, habitat modification, and reduced water quality) related to the need for socio-economic development
D
(Largely modified)
Substantial changes in natural habitat, biota, and essential ecosystem functions have occurred; a lower than expected species richness; the much-lowered presence of intolerant species; alien species prevail.Significant and precise visible disturbances (such as dams, diversions, transfers, habitat modification, and water quality degradation) associated with basin and water resources development
E
(Seriously modified)
Habitat diversity and availability have declined; species richness is strikingly lower than expected; only tolerant species remain; indigenous species can no longer breed; alien species have invaded the ecosystem.High human population density and extensive water resources exploitation. This class is not suitable as a management goal. The management team should move to a higher class to restore the flow pattern of the river.
F
(Critically modified)
Modifications have reached a tipping point; the ecosystem has been completely modified with an almost complete loss of natural habitat and biota; in the worst case, the underlying ecosystem functions have been destroyed, and changes are irreversible.This status is not acceptable from the management perspective. Management interventions are necessary to restore flow patterns and river habitats (if still possible/feasible) to ‘move’ a river to a higher management category.
Table 4. Flow recommendations as per the Tennant method, based on Tennant [45].
Table 4. Flow recommendations as per the Tennant method, based on Tennant [45].
Aquatic-Habitat Condition for Small StreamRecommended Base Flow (% of MAF)
General Period (October–March)Fish Spawning Period (April–September)
Flushing or maximum200% of the average flow
Optimum range60–10060–100
Outstanding4060
Excellent3050
Good2040
Fair or degrading1030
Poor or minimum1010
Severe degradation<10% of average flow to zero flow
Table 5. List of the Indicators of Hydrologic Alteration (IHA) parameters, their ecological significance, and regime characteristics and global indexes for each IHA group.
Table 5. List of the Indicators of Hydrologic Alteration (IHA) parameters, their ecological significance, and regime characteristics and global indexes for each IHA group.
Global Index for Each GroupIHA ParametersRegime Characteristic (Specific Alteration)Ecological Significance
Mean Monthly Flow Alteration Index (Imm)Group 1:
Mean value of each calendar month
Magnitude
(increased variation)
Guaranteed favourable habitat conditions and flow regime (quantity, quality, and temperature) for aquatic and terrestrial organisms. Availability of food and cover for fur-bearing mammals.
Magnitude and Duration of Extreme Flow Alteration Index (IMDE)Group 2:
Annual minima, 1, 3, 7, 30, 90 day means
Annual maxima, 1,3,7,30,90 day means
Number of zero-flow days
Baseflow index: 7 day minimum flow/mean flow for the year
Magnitude and Duration (prolonged low flows; altered inundation duration; prolonged inundation)Structuring of aquatic ecosystems by abiotic and biotic factors. The shaping of river channel morphology and physical habitat conditions.
Timing of Extreme Flow Alteration Index (ITE)Group 3:
Julian date of each annual 1 day maximum
Julian date of each annual 1 day minimum
Timing (oss of seasonal flow peaks)Disrupt cues for fish: (spawning, egg hatching, migration) [54]. Evolution of the life history and behaviour mechanism of the aquatic organisms [48].
Frequency and Duration Alteration Index (IFD)Group 4:
No. of high pulses each year
No. of low pulses each year
Mean duration of high pulses within each year (days)
Mean duration of low pulses within each year (days)
Frequency and Duration
(flow stabilization)
Availability of floodplain habitats for aquatic organisms. Influences bedload transport, channel sediment textures, and duration of substrate disturbance (high pulses). Nutrient and organic matter exchanges between river and floodplain.
Rate and Frequency Alteration Index (IRF)Group 5:
Means of all positive differences between consecutive daily values Means of all negative differences between consecutive daily values
Reversals
Rate of change and Frequency (rapid changes in river stage; accelerated flood recession)Wash out and stranding of aquatic species [55]. Failure of seedling establishment [56].
Note: All the ecological significance of the IHA parameters is not listed in the table. A few are listed for more ecological significance; the authors recommend literature related to IHA parameters [51,53,57,58].
Table 6. Range of alteration for global indexes. Where ‘0’ means no alteration while ‘1’ implies the highest alteration.
Table 6. Range of alteration for global indexes. Where ‘0’ means no alteration while ‘1’ implies the highest alteration.
Range0.00–0.250.25–0.500.50–0.750.75–1.00
AlterationLowMildModerateHigh
Table 7. E-flow allocation of Kaligandaki River regarding all applied e-flows methods (EFMs).
Table 7. E-flow allocation of Kaligandaki River regarding all applied e-flows methods (EFMs).
MethodClasses(%of MAF)E-flows (m3/s)
GEFCClass B47.8214.46
Class C32.8147.16
Class D23.7106.33
Class E18.683.45
Class F15.770.44
TennantOct–Mar1044.87
Apr–Sept30134.6
FDCQ80% FDCA, 49.04
Q85% FDCA, 45.65
Q90% FDCA, 43.05
Mean annual flow10%MAF 44.97
Dynamic methods30%Q-D 30% of daily flow
Annual Distribution Method
MonthJanFebMarchAprilMayJuneJulyAugustSepOctNovDec
E-flow (m3/s)60.5750.0545.0749.7573.54204.58646.33785.59584.67274.83127.6183.71
Table 8. Relative change (%) of the EFC regarding six EFMs against the natural flow regime (NFR), represented by the mean values. The sign ( + ) symbolizes an increase and (−) a decrease.
Table 8. Relative change (%) of the EFC regarding six EFMs against the natural flow regime (NFR), represented by the mean values. The sign ( + ) symbolizes an increase and (−) a decrease.
E-Flows Components (EFC)Dynamic E-FlowsMinimum AnnualFDC CurveADMGlobal Environmental Flow CalculatorTennant
EFC Low Flows30%Q-D10%MAFQ80%Q85%Q90%Class BClass CClass DClass EClass F
July—Low Flow−70−100−100−100−100−46−100−100−100−100−100−70
August—Low Flow−70−100−100−100−100−65−100−100−100−100−100−61
September—Low Flow−70−100−100−100−100−35−100−100−100−100−100−72
October—Low Flow−70−100−100−100−100−40−100−100−100−100−100−64
November—Low Flow−70−100−100−100−100−49−100−100−100−100−100−41
December—Low Flow−70−100−100−100−100−49−100−100−100−100−100−100
January—Low Flow−70−100−100−100−100−48−100−100−100−100−100−100
February—Low Flow−70−100−100−100−100−49−100−100−100−100−100−100
March—Low Flow−70−100−100−100−100−50−100−100−100−100−100−100
April—Low Flow−70−100−100−100−100−48−100−100−100−100−100−6
May—Low Flow−70−100−100−100−100−34−100−100−100−100−100−15
June—Low Flow−70−100−100−100−100−31−100−100−100−100−100−52
EFC Parameters
Extreme low peak−70−100−100−100−100−41−100−100−100−100−100−39
Extreme low duration0−100−100−100−100−55−100−100−100−100−1001393
Extreme low timing0−100−100−100−100−3−100−100−100−100−100234
Extreme low freq.0−100−100−100−10082%−100−100−100−100−100−74
High flow peak−70−100−100−100−100−62−100−100−100−100−100−100
High flow duration0−100−100−100−100−30−100−100−100−100−100−100
High flow timing0−100−100−100−10026−100−100−100−100−100−100
High flow frequency0−100−100−100−100−2−100−100−100−100−100−100
High flow rise rate−70−100−100−100−100−74−100−100−100−100−100−100
High flow fall rate−70−100−100−100−100−81−100−100−100−100−100−100
Small Flood peak−70−100−100−100−100−68−100−100−100−100−100−100
Small Flood duration0−100−100−100−10011−100−100−100−100−100−100
Small Flood timing0−100−100−100−1006−100−100−100−100−100−100
Small Flood frequency0−100−100−100−100−64−100−100−100−100−100−100
Small Flood rise rate−70−100−100−100−100−96−100−100−100−100−100−100
Small Flood fall rate−70−100−100−100−100−85−100−100−100−100−100−100
Large flood peak−70−100−100−100−100−74−100−100−100−100−100−100
Large flood duration0−100−100−100−100−18−100−100−100−100−100−100
Large flood timing0−100−100−100−100−2−100−100−100−100−100−100
Large flood frequency0−100−100−100−100260−100−100−100−100−100−100
Large flood rise rate−70−100−100−100−100−88−100−100−100−100−100−100
Large flood fall rate−70−100−100−100−100−71−100−100−100−100−100−100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Suwal, N.; Kuriqi, A.; Huang, X.; Delgado, J.; Młyński, D.; Walega, A. Environmental Flows Assessment in Nepal: The Case of Kaligandaki River. Sustainability 2020, 12, 8766. https://0-doi-org.brum.beds.ac.uk/10.3390/su12218766

AMA Style

Suwal N, Kuriqi A, Huang X, Delgado J, Młyński D, Walega A. Environmental Flows Assessment in Nepal: The Case of Kaligandaki River. Sustainability. 2020; 12(21):8766. https://0-doi-org.brum.beds.ac.uk/10.3390/su12218766

Chicago/Turabian Style

Suwal, Naresh, Alban Kuriqi, Xianfeng Huang, João Delgado, Dariusz Młyński, and Andrzej Walega. 2020. "Environmental Flows Assessment in Nepal: The Case of Kaligandaki River" Sustainability 12, no. 21: 8766. https://0-doi-org.brum.beds.ac.uk/10.3390/su12218766

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop