A Comparative Study for Provision of Environmental Flows in the Tapi River

: Environmental ﬂows are deﬁned as the ﬂow required into a stream to maintain the river’s ecosystem. The notion of Environmental Flow Allocation (EFA) ensures that a sufﬁcient amount of water is delivered to the stream to maintain ecological integrity. The objective of this study is to examine environmental ﬂows and determine the best acceptable strategy for providing ﬂows into the river in the Lower Tapi Basin. To achieve this objective, daily discharge data from three sites, Ukai (period 1975–2020), Motinaroli (period 1990–2021), and Ghala (period 1995–2005) were collected and analyzed using the Tennant, Tessman, variable monthly ﬂow (VMF), and Smakhtin methodologies. A comparative analysis was carried out on all three sites using the four methodologies. The Tessman and VMF approaches have a strong connection with the computed environmental ﬂow requirements (EFR), according to the results. The calculated EFR was found to be in the range of 30–35% of mean annual ﬂows (MAF). The maximum EFR found at station Ghala is about 54.5% of MAF according to the Tessman method. Such research will help to prevent future degradation of the river by supplying ﬂow in accordance with the EFR, and it will also be used by stakeholders and policymakers to allocate water to preserve the ecosystem.


Introduction
Rivers constitute one of the major parts of the water cycle on the earth and carry basic natural requirements of all the natural life events and ecosystem functions. Since the beginning of time, people have been fascinated by rivers because they rely on them for their water supply, and, as a consequence, the earliest human civilizations developed along the rivers and gradually started to control them, draining and diverting more water as population and water demand increased [1]. In many situations, this led to some rivers having almost no amount of flow in the dry season, causing various forms of damage to the numerous ecosystem services that river water provides. The loss of feeding and breeding habitats in the floodplain water bodies due to the construction of embankments and increased silt load and macrophyte growth are major causes for declining fish resources. In fact, wetlands, rivers, and lakes support a very large share of biodiversity on Earth, even though freshwater ecosystems account for a negligible (0.01%) amount of the planet's water and a small portion of its surface [2].
In this context, the concept of environmental flows emerges as discussions increasingly expanded to cover a wide range of subjects, such as the geomorphic, sediment movement, Earth 2023, 4 572 physically possible [15]. Even the current National Water Policy places "ecology" fourth on the list of water-allocation objectives. In the Indian context, primarily hydrologic methods are used to define environmental flows in rivers. As the progressive destruction of the water environment became apparent, environmental concerns began to gain traction. The issues of EF assessment and management are high on the world agenda at present time. At this time, it remains a new research field. In many countries, including India, there has not even been a nationwide assessment of the water requirements of rivers and their associated aquatic ecosystems. Countries like India, which are also dealing with high population growth nowadays, also have to deal with the allocation of water for irrigation and drinking water purposes making it difficult to allot water to maintain the environmental flow as it is actually needed.
In this context, this study aims to find the environmental flow requirements (EFRs) for the Tapi River (India) by using four different methods, trying to find the best-fitted method among all. The result of the study helped the Ukai Irrigation Circle in releasing the flow from the Ukai dam site. There are significant physical changes in the river that have been observed, which is why the study area was selected. The study is helpful for the Surat Municipal Corporation for the maintenance of the river channel. The flow between the weir and mouth of the river should not have a sufficient amount of flow in the river channel. Flow in this part of the river varies to a great extent between monsoon season and non-monsoon season. This part is also affected by the river meander and sedimentation. After seeing all these problems, an assessment of the EF was required to be carried out at the site.

Methods Used for Assessment of Environmental Flows
Before the 1990s, water management was limited to minimum flow requirements and water quality standards. The last two decades have seen a shift in river management that is more attentive to the seasons and natural circulation regimes of interannual flow variability, as well as the magnitude, timing, and frequency of various flow conditions. According to a thorough investigation of environmental flow techniques conducted by Tharme [16], there are more than 207 significant environmental flow methodologies (EFMs) in use across 44 countries and six regions of the world. These EFMs can be broadly divided into four categories, including [17]: hydrological, hydraulic, habitat-simulation, and holistic (Table 1). Figure 1 shows the graphical representation of the timeseries of the discharge data.
In particular, hydrological methods look at historical runoff series to determine natural flow conditions and advised flow rates. These methods are easy to perform and inexpensive, and they can be used for the preliminary stage of the dam construction. The method is based on the premise that there is a connection between river flow and riverbed ecosystems. The method utilizes time series data of 20-50 years (long term) for the flow of the river, with the ecological flow represented by a proportion of mean yearly flow values. hydrologic index methods include the Tessman method, Tennant method, RVA approach (range of variability), DRM (Desktops Reserve Model), FDC (flow duration curve) analysis method, and the HAT (Hydrologic Assessment Tool). In particular, hydrological methods look at historical runoff series to determine natural flow conditions and advised flow rates. These methods are easy to perform and inexpensive, and they can be used for the preliminary stage of the dam construction. The method is based on the premise that there is a connection between river flow and riverbed ecosystems. The method utilizes time series data of 20-50 years (long term) for the flow of the river, with the ecological flow represented by a proportion of mean yearly flow values. hydrologic index methods include the Tessman method, Tennant method, RVA approach (range of variability), DRM (Desktops Reserve Model), FDC (flow duration curve) analysis method, and the HAT (Hydrologic Assessment Tool).

Tennant Method
This method, also known as the Montana method, is the most widely used hydrological technique in the world [19]. It was created in 1976 and has been used in a minimum of 25 countries in either its original form or by applying modification to the original form [16]. The method was specially developed to determine the needs of fish in streams of Montana. A different range of flow is given for each different stage of river condition. The flow condition is defined using a proportion of mean annual flows (MAF) for two separate six-month periods. The selection of the mean annual runoff, which only depicts run-off on an annual basis and ignores variations in flow throughout the year, is its main drawback. Many researchers have concluded that, in the monsoon period, 40% of MAF is needed and in the non-monsoon period 20% of MAF is needed to maintain good habitat conditions [16].

Tessman Method
Tennant method gives the EFR for the two different seasons. To overcome this disadvantage of the Tennant method, the Tessman method was developed [20]. The flow of a river varies greatly throughout the seasons (monsoon and non-monsoon). Therefore, the calculation of EFR by considering only MAF is not adequate and the Tessman method that

Tennant Method
This method, also known as the Montana method, is the most widely used hydrological technique in the world [19]. It was created in 1976 and has been used in a minimum of 25 countries in either its original form or by applying modification to the original form [16]. The method was specially developed to determine the needs of fish in streams of Montana. A different range of flow is given for each different stage of river condition. The flow condition is defined using a proportion of mean annual flows (MAF) for two separate six-month periods. The selection of the mean annual runoff, which only depicts run-off on an annual basis and ignores variations in flow throughout the year, is its main drawback. Many researchers have concluded that, in the monsoon period, 40% of MAF is needed and in the non-monsoon period 20% of MAF is needed to maintain good habitat conditions [16].

Tessman Method
Tennant method gives the EFR for the two different seasons. To overcome this disadvantage of the Tennant method, the Tessman method was developed [20]. The flow of a river varies greatly throughout the seasons (monsoon and non-monsoon). Therefore, the calculation of EFR by considering only MAF is not adequate and the Tessman method that considers the mean monthly flow (MMF) to determine the EFR in terms of minimum flow is more appropriate. In particular, the minimum flow is evaluated as:

Variable Monthly Flow (VMF) Method
The variable monthly flow (VMF) method was proposed by Pastor [18] to explain the environmental flow requirements in global water assessments. EFRs are calculated using monthly flow, which follows the totally natural discharge variability. this method. The method was used to increase the protection of the freshwater ecosystem during the low-flow season with a reserve of 60% of the MMF and a minimum considered flow of 30% of the MMF during the high-flow period. This method allows stakeholders to use water during the low-flow season by losing up to 40% of the MMF. The minimum flow is evaluated as:

Smakhtin Method
Smakhtin considered the indices from the flow duration curve for the computation of EFR for the allocation of water into the stream. Smakhtin et al. [21] first applied this method to 128 catchment areas in different parts of the world to assess the status of exploitation of the world's rivers with the consideration of environmental water requirements. According to this method, EFR is roughly represented by a discharge which, on average each year, is exceeded 90% of the time in order to maintain "fair" ecological conditions, and Q90 is the flow that is considered for the allocation of EFR in the stream for 90% of the year. In particular, through reviewing numerous studies on low-flow hydrology, Smakhtin et al. [21] proposed that Q90 varies between 0 and 50% of MAR.

Correlation Analysis
EFR was calculated by taking an average from all four methods and denoted as the calculated EFR. Then, the relationship between the EFR for all the methods and the calculated EFR on a monthly basis was analyzed by means of a simple linear correlation.

Summary of the Paper
In the following schematic working flow diagram, a summary of the paper with the different steps performed to detect the EFR is shown ( Figure 2).

Study Area
The Tapi River is the tenth largest on India's peninsula. After the Narmada, it is the second largest westward river. Rivers drain from such parts of Madhya Pradesh, as well as parts of Maharashtra and Gujarat. The Tapi River Basin (TRB) is located between longitudes 72°33′ to 78°17′ and latitudes 20°9′ to 21°50′. It begins at an altitude of 752 m near the Multai reserve forest in Betul, Madhya Pradesh state. The river is 724 km long from its source to its mouth of the Arabian Sea. The left-bank tributaries are the Vaghur, the Am-

Study Area
The Tapi River is the tenth largest on India's peninsula. After the Narmada, it is the second largest westward river. Rivers drain from such parts of Madhya Pradesh, as well as parts of Maharashtra and Gujarat. The Tapi River Basin (TRB) is located between longitudes 72 • 33 to 78 • 17 and latitudes 20 • 9 to 21 • 50 . It begins at an altitude of 752 m near the Multai reserve forest in Betul, Madhya Pradesh state. The river is 724 km long from its source to its mouth of the Arabian Sea. The left-bank tributaries are the Vaghur, the Amravati, the Bori, the Girna, the Purna, and the Sipna, and the right-bank tributaries are the Suki, the Arunavati, the Gomai, and the Aner. The Tapi River Basin covers 65,145 square kilometers, which is equivalent to 2% of India's total land area. The Tapi River basin receives a maximum annual rainfall of up to 2550 mm, with an average of about 835 mm. Madhya Pradesh has a catchment area of 9804 square kilometers, Maharashtra has a catchment area of 51,504 square kilometers, and Gujarat has a catchment area of 3837 square kilometers, these being the three states that make up the catchment area of the river. A catchment is a basin-shaped area of land, bounded by natural features such as hills or mountains, from which surface and sub-surface water flow into streams, rivers, and wetlands. Water flows into, and collects in, the lowest areas in the landscape. The main Tapi basin is separated into three sub-basins: the Lower Tapi Basin (LTB) from Sarangkheda to the Sea, the Middle Tapi Basin from Hathnur to the influence of Purna with the main Tapi, and the Upper Tapi Basin up to Hathnur (a combined area of 29,430 square km) (6745 square km). The entire basin has a significant amount of of flow obstruction structures like dams, weirs, barrages, etc., on the stream. About 90 percent of total rainfall is received during the monsoon months, of which 50% is received during July and August. The Tapi River basin shows different climatic characteristics due to the variation in topography from the upper to lower part of the basin. The average rainfall in the Tapi basin is 888.0 mm. The hydrological response of a river basin can be interrelated with the physiographic characteristics of the drainage basin, such as size, shape, slope, drainage density and size, length of the streams, etc. Specifically, 358 small and large dams, mostly earthen type, can be identified in the basin. Ukai dam is the highest dam in the basin (81 m), and construction was completed in the year 1972. The study area was selected in the LTB, between Ukai dam to the mouth of the Arabian Sea, primarily with Surat and Hazira twin cities, as well as tens of small towns and villages along the river's course. The downstream reach of the Ukai dam is affected by the degradation of the river and alteration of flow. The flow of water into a stream fluctuates significantly due to the storage of water and diversion of water through dams for different water supply needs; lower basins mostly covered the plain of Surat city. The twin cities of Surat and Hazira are nearly 106 km downstream of the Ukai Dam and are affected by recurring floods at regular intervals. Figure 3 shows the location of the Tapi River Basin in India and estimated EF sites in the lower Tapi Basin.
The entire basin is divided into three sub-basins from which the Lower Tapi basin was selected for the study. A stream map of the Lower Tapi basin is also presented. The red-color line indicates the location of the Ukai Dam. The TRB has a total of 13 CWC gauging stations across the basin and 4 of these are located in Gujarat while the others are in Maharashtra and Madhya Pradesh. After obtaining flow data, three stations were selected for the analysis: Ukai (EF1), Motinaroli (EF2), and Ghala (EF3). The following table lists the features of the sites, as well as their names, locations, and characteristics. The section of the river from Ukai Dam to the mouth of the river at the Arabian Sea was considered in LTB, primarily occupying Surat along with Hazira twin cities, as well as tens of small towns and villages along the river's course. The feature of the site and period of the hydrologic datasets is represented in Table 2. Ukai dam is affected by the degradation of the river and alteration of flow. The flow of water into a stream fluctuates significantly due to the storage of water and diversion of water through dams for different water supply needs; lower basins mostly covered the plain of Surat city. The twin cities of Surat and Hazira are nearly 106 km downstream of the Ukai Dam and are affected by recurring floods at regular intervals. Figure 3 shows the location of the Tapi River Basin in India and estimated EF sites in the lower Tapi Basin. The entire basin is divided into three sub-basins from which the Lower Tapi basin was selected for the study. A stream map of the Lower Tapi basin is also presented. The

Results
The discharge data on a daily basis of the Tapi River were utilized to evaluate the environmental flows for water allocation. The study area was analyzed by considering several hydrological indicators, i.e., mean annual flow (MAF) at a particular site, mean monthly flow (MMF) at a particular site, base flow index, which is computed by BFI = Q90/MAF, and hydrological variability index, which is computed by HVI = Q25 − Q75/Q50, where Q90, Q75, Q50, and Q25 are annual flows that are equal to or exceeded 90%, 75%, 50%, and 25% of the time, respectively ( Table 3).
The EFRs in the Tapi River were calculated by means of the Tennant, Tessman, VMF, and Smakhtin methods for determining a minimum flow for the lower reach of the Tapi River. The results and the comparisons of the four methods are represented in Table 4 and Figures 3 and 4.  (HFRs). The environmental flows requirement (EFR) for all three stations was obtained at a monthly scale and the comparison of the values evaluated with the different methods is shown in Figure 4. As a result, for the EF1 (Ukai) station, the lowest EFR values were obtained by means of the Smakhtin method, while the VMF method shows the highest value of EFR during the monsoon season (Figure 4a).  At this station, August, September, and October were identified as the months with the highest flows. Moreover, the calculated EFR is highest during the high flow months, while the rest of the months show an average EFR value. As regards the station EF2 (Motinaroli), the Smakhtin and the Tessaman methods showed the lowest and highest values of EFR, respectively (Figure 4b) and the highest flows were detected in July, August, and September. The calculated EFR showed the highest values during the monsoon season (close to the Tessman method ones), while it showed moderate value in the rest of the months. Finally, also at station EF3 (Ghala), the lowest and highest values of EFR were detected with the Smakhtin and Tessaman methods, respectively (Figure 4c). Moreover, August, September, and October were identified as the highest flow months. And for October to March, the Smakhtin and the Tennant methods showed almost equal values of EFR.
For the different methods and for the different stations, Figure 5 shows the ratio between the annual EFR and the mean annual flow (MAF). Generally, for all the stations, the Tessman method showed the highest percentage values while the Smakhtin method shows the lowest ones. In particular, considering all the values obtained by all four methods and for all the stations, the lowest values were observed at the Motinaroli station for the Smakhtin method while the highest values were observed at the Ghala station for the Tessman method. In fact, results showed that, for the Ukai station, the EFR values obtained by means of the Tennant, Tessman, VMF, and Smakhtin methods resulted in 33%, 49%, 37%, and 20% of the MAF, respectively. Similar behavior was detected at the Ghala station, in which the highest percentage was evaluated for the Tessman method (about 53%), followed by the VMF method, the Tennant method, and the Smakhtin method (16.2%). Regarding the Motinaroli station, in this case the Tessman and the Smakhtin methods also showed the highest (about 47%) and the lowest (about 13%) percentages but, differently from the previous stations, the Tennant method showed percentage values higher than the VMF one. For the different methods and for the different stations, Figure 5 shows the ratio tween the annual EFR and the mean annual flow (MAF). Generally, for all the statio the Tessman method showed the highest percentage values while the Smakhtin meth shows the lowest ones. In particular, considering all the values obtained by all four me ods and for all the stations, the lowest values were observed at the Motinaroli station the Smakhtin method while the highest values were observed at the Ghala station for Tessman method. In fact, results showed that, for the Ukai station, the EFR values obtain by means of the Tennant, Tessman, VMF, and Smakhtin methods resulted in 33%, 49 37%, and 20% of the MAF, respectively. Similar behavior was detected at the Ghala s tion, in which the highest percentage was evaluated for the Tessman method (about 53 followed by the VMF method, the Tennant method, and the Smakhtin method (16.2 Regarding the Motinaroli station, in this case the Tessman and the Smakhtin methods a showed the highest (about 47%) and the lowest (about 13%) percentages but, differen from the previous stations, the Tennant method showed percentage values higher th the VMF one.  Figure 6 illustrates the correlation between the EFRs determined using the four d ferent methods and the calculated EFR values. Results evidenced a very good correlat between the calculated EFR and the Tessman and the variable monthly flow (VMF) me ods, with correlation coefficient (R 2 ) values higher than 0.9. A good correlation was a  Figure 6 illustrates the correlation between the EFRs determined using the four different methods and the calculated EFR values. Results evidenced a very good correlation between the calculated EFR and the Tessman and the variable monthly flow (VMF) methods, with correlation coefficient (R 2 ) values higher than 0.9. A good correlation was also identified for the Smakhtin method, which presented correlation coefficients of about 0.63. Finally, the lowest values were detected for the Tennant method with a R 2 lower than 0.5.
Earth 2023, 4, FOR PEER REVIEW 12 Figure 6. Linear relationship between the EFR for all the methods and the calculated EFR on a monthly basis.

Conclusions and Discussion
In this study, the EFR was estimated using the Tennant, Tessman, variable monthly flow (VMF), and Smakhtin methods, starting from daily data measured in three-gauge stations (Ukai, Motinaroli, Ghala) in the lower reach of the Tapi River.
Generally, data on living organisms, the geometric data of the river, and real knowledge of the ecological effects of flow alternation are needed to correctly define EFs, but these data on the river ecosystem are not always available, an issue identified in past studies in the majority of the world's river basins [18Error! Reference source not found.]. For this reason, the choice of EF methods for this study was limited to hydrological meth-

Conclusions and Discussion
In this study, the EFR was estimated using the Tennant, Tessman, variable monthly flow (VMF), and Smakhtin methods, starting from daily data measured in three-gauge stations (Ukai, Motinaroli, Ghala) in the lower reach of the Tapi River.
Generally, data on living organisms, the geometric data of the river, and real knowledge of the ecological effects of flow alternation are needed to correctly define EFs, but these data on the river ecosystem are not always available, an issue identified in past studies in the majority of the world's river basins [18]. For this reason, the choice of EF methods for this study was limited to hydrological methods that were largely applied in past studies in which these methods were compared. For example, Hairan et al. [22] carried out a study on the evolution of the environmental flow of the Selangor River (Malaysia) using the hydrological index approach and the software Global Environmental Flow Calculator (GEFC). They analyzed daily flow data collected from Malaysia's Department of Irrigation and Drainage (DID) over a 60-year period  and demonstrated that, from the overall comparison of the pre-and post-impact periods, the minimum flow in the dry seasons of the year increased, while the maximum flow decreased in the monsoon seasons during the post-impact period. Kumar et al. [23] applied three separate hydrological methodologies for EFA to calculate the minimal environmental flow requirements for the Srisailam dam (India) for the period 1968-2015. They identified that the reservoir operation policy was causing a severe hydrological alteration in the high flow season, especially in the month of July. Moreover, they concluded that in the case of the non-availability of environmental information, hydrological indicators can be used to provide the basic assessment of environmental flow requirements. Sahoo et al. [24] carried out an environmental assessment on the Mahanadi River Basin (India) by using the hydrological approach of environmental flow assessment. Using GEFC software, they evaluated the usefulness of numerous desktop hydrology-based environmental flow assessment methods such as Tennant, Tessman, variable monthly flow (VMF), range of variability approach (RVA), and FDC shifting methodology to evaluate environmental flows from daily discharge data.
Following these approaches, in this paper, the EFR values evaluated with the different methods were compared and, for each gauge station, the combination of the four methods, given by the average value of the EFR, was evaluated together with the low flow requirement (LFR) and the high flow requirement (HFR). Moreover, the results obtained with the four different methods were compared with the calculated EFRs in three different stations, also considering monthly variability, which indicates a significant improvement over previous studies based only on an annual scale [21,25]. Finally, the correlation between the EFRs determined using the four different methods and the calculated EFR values was investigated. In this study, the Tessman and the VMF methods gave the highest correlation with the calculated EFR. This result has also been evinced in past studies, such as that of Pastor et al. [18], in which these two methods showed good correlation with locally calculated EFRs in several case studies from a wide range of climates, flow regimes, and freshwater ecosystems. The timing and duration of low flows are critical to the health of the river ecosystem. In this context, the present study ensures the minimum value of flow necessary to maintain the sustainability of the watercourse. In any case, because of the lack of data on ecosystem response, the EF approaches in this study were limited to hydrological methods. In order to better define EFs, biological data, and river morphology data, an understanding of the ecological implications of flow alternation could be included in future studies and, thus, more information on ecological data should be gathered, as well as data from fisheries, as the timing and duration of low flows are essential to aquatic ecosystems' health. Moreover, this study could be further improved by considering the whole Tapi River Basin.