Development of an Evaluation System for Sustaining Reservoir Functions — A Case Study of Shiwen Reservoir in Taiwan

Reservoirs serve as important facilities, providing a stable source of public water in Taiwan. As construction of new reservoirs becomes more difficult, it is essential to ensure sustainable functionality of reservoirs in the future. To this end, this study proposes a system for reservoir sustainability evaluation. The evaluation system consists of social justice, environmental protection, and economic development containing 12 indicators which are grouped into six categories: flood control, sediment management, water resources allocation, river ecology, water quality, and benefit and fairness. Moreover, evaluation system operational procedures to supplement planning and decision-making processes are proposed, and applied in a case study of the Shiwen reservoir planning in Taiwan. The planned reservoir in this case study is rated as “Good”, nearly “Excellent”, in sustainability as evaluated with the Sustainability Confidence Index (SCI). Additionally, Analytic Network Process (ANP) results indicate that the flood control capacity and sediment management are the first and second most important indicators for the reservoir. If desilting operations had been conducted, the SCI values would have increased from 3.3 to 3.7, warranting an “Excellent” rating for the reservoir. The case study demonstrates that decision-makers can apply the proposed system when managing reservoir evaluations.


Introduction
In general, human activities subsequently lead to many environmental impacts on biodiversity, ecosystem functioning, and water cycle conservation efforts [1], yet reservoirs comprise an effective water infrastructure for integrated water resources development and management [2][3][4].Changes in the natural environment are thus an inevitable consequence of reservoir construction and operation.Sustainability is one popular concept delineated in Deciding the Future of Dams and Research [5].For a sustainable approach to reservoirs, a sediment management plan is necessary [6].Recently, numerous studies have indicated that Integrated Water Resources Management (IWRM) should consider and seek a precautionary approach to manage aquatic habitats [1,[7][8][9][10] as well as sedimentation and flood controls.Reservoir sedimentation reduces the demand of a reservoir and can be mitigated by actions such as sediment bypassing, sediment sluicing, or sediment flushing to avoid sediment accumulation [11].Flood potential is a significant parameter when considering the tradeoff between flood control and conservation for sustaining reservoir functions in integrated water resources development [3].Concerning water resources management, as sediment accumulation and potential flood problems not only decrease a reservoir's demand but also influence riverine ecosystems [1,[7][8][9][10]12] (e.g., fish habitats and downstream water quality), an effective evaluation system for sustaining multi-objective reservoir functions is essential to integrate water resources development and management.Accordingly, there are many factors which influence the functions of a reservoir, such as flood, sedimentation, water quality, and surrounding habitats [4,11,13,14].Those factors with additional indices should be considered as multiple indicators for sustainable reservoir management.

Flood Control Indicator
Flood control systems include the reservoir itself and its downstream river channel, the primary purpose of which is to reduce flood flow due to considerations for the safety of both the reservoir and river channel [15].Geng et al. [16] used flood control capability indicators to perform assessments that included a river level regulating index and flood control capability rate of change indices.Hsieh [17] integrated current flood status data with a wetness index to create a flood early warning index featuring signal lights, called the flood alert indicator (FAI).Recently, Chen [18] included minimum and maximum reservoir water levels of flood control for multiple objective operations.Some studies on flood control operations take into account inflow, non-damaging discharge, and flood peak [19,20], while other studies on flood control capability assessment mostly focus on reservoir safety and river flood flow standards.Scarrott et al. [21], Montaldo et al. [22], and Gioia [23] used reservoir routing approach in flood control.

Sediment Management Indicator
In sediment management, however, the service life of a dammed reservoir and downstream fish habitat quality is highly influenced by sedimentation [24].Numerous indicators have been used for basin-wide sediment management-a management concept that involves reducing sediment yield in the upper catchment areas, enhancing in-stream sediment transport capacity, and desilting reservoir sediment.A drawdown ratio defined by Atkinson [25] has been used as an indicator to evaluate the elevation of river outlets or desilting outlets (including flushing, sluicing, and bypassing) in Japan.Graf et al. [26] employed a mean annual sedimentation rate as an indicator to investigate and analyze life spans of reservoirs across the United States, and then graded reservoir sedimentation states according to a five-level rating scheme based upon annual loss of reservoir capacity.

Water Resources Allocation Indicator
Water resource management concentrates on social and economic benefits and equality as well as on maintaining the ecosystem sustainability [27,28].Water resource allocation is a key function of reservoirs.Many indicators for stress, allocation, or water resource shortage are currently used to evaluate the possibility of water shortage [2,29,30].The shortage index (SI) is one of the established water shortage indices [2,29,30] that reflects the state of water resource allocation.Water shortage level results can then be used to revise water resource allocation strategies.In 2004, the Taiwan Water Resources Agency [31] established a three-level water supply stability indicator to evaluate water resource allocation.Huang and Chou [32] proposed regulations for adjusting water resource allocation based on Shihmen Reservoir's operational rule curves, and further established a drought index for reservoir storage in Taiwan.Zeng et al. [2] applied SI to verify the effectiveness of the proposed water transfer triggering mechanism for multi-reservoir operation performances.Tseng et al. [29] used SI to evaluate drought events and their possible impacts to the water resources system of the Tsengwen Reservoir of Taiwan.Lian et al. [30] applied SI for optimal Danjiangkou Reservoir operations in China.For domestic water use purposes, reservoir water quality is the most important part of sustaining water supply functionality, and is the topic of many studies producing various related indicators such as the water quality index (WQI) for reservoir operations and evaluations [12,33,34].Moreover, water quality is one popular ecological indicator in sustainable development [35].

River Ecological Indicator
For river ecology conservation, adequate reservoir flow release is necessary to mitigate adverse effects on aquatic organisms' natural habitats caused by changes in downstream flow.Riverine ecosystem protection needs to maintain natural flow and sediment regimes downstream from dams [36].This is needed to integrate environmental flows into multi-objective reservoir management [9].Only recently have river ecosystems, including fish habitats and environmental flows, been considered in reservoir operations [18], though the valuable insights of stream habitat ecological modeling and assessment have been widely identified in previous studies [37][38][39].Recently, numerous studies applying the IWRM concept have incorporated precautionary measures for aquatic habitats [1,[7][8][9][10].For such aquatic habitat studies, ecological flows can be determined using statistical methods, such as the Indicators of Hydrologic Alteration (IHA) method [9].Xu et al. [15] used an eco-friendly reservoir operating rule for restoring a fish migration passage into their reservoir operating rules.To manage reservoir-related fish habitat issues, river habitat models can be employed to simulate changes in river flow and habitat availability under various flow scenarios.As such, weighted usable area (WUA) indices have been widely applied in habitat assessments including reservoir operations [33,[38][39][40][41][42][43].Nikghalb et al. [44] determined ecological flow using WUA, and compared these results with Q 95 results.Aguilar and Polo [45] tried using various approaches to determine minimum environmental flow values including using WUA with various flows.Similar to the studies of Li et al. [40], Pasha et al. [41], Yao et al. [42], and Xu et al. [15], this study used WUA as an indicator with Q 95 to represent the ecological condition of the study reservoir.However, in Taiwan, the concept of a minimum discharge required by fish has been reformulated into a relation between volume and habitat [46].A WUA compatible with the living requirements of fish within the research area is then calculated.

Gini Coefficient
Equity, stability, and economic efficiency are key components in decision-and policymaking for sustainable development [28].The Gini coefficient is a ratio of difference of ideal equity and real uses [47].The coefficient can be one of the effective indices employed to explore equity in water allocation strategies [28].The Gini coefficient here was calculated based on research on the fairness of water allocation in South Africa [48], future changes on water discharge in the coming 30 years (e.g., Masaki et al. [49]), and the inequality in water supply and demand from 1999 to 2006 in the Yellow River basin (e.g., Wang et al. [50]).Wang et al. [50] used the Gini coefficient in assessing equity in domestic water supply.Masaki et al. [49] used the Gini and Lorenz-asymmetry coefficients to evaluate the characteristics and applicability in river flow regimes under future climate change.Rajah et al. [51] applied a Gini index to evaluate the precipitation distributions for wet-day frequency and decreasing wet-day frequency.Peipoch et al. [52] applied the Gini coefficient to evaluate the effects of ecological simplification on flow variability.Hu et al. [28] used the Gini coefficient to find equity in the water allocation strategies in a river basin.Golusin and Ivanović [35] used the Gini coefficient, water quality, and water supply as indicators for Southeastern Europe countries.

Sustainable Indicators
Numerous researchers in various fields of expertise have established several indicators to determine the sustainability of a planned multipurpose reservoir.To analyze water supply, Xu et al. [53] defined a system dynamic model using a set of sustainability confidence indices (SCI) to evaluate the level of sustainability according to water supply and consumption of various water resource systems along the Yellow River in China.Huang [54] studied the application of SCI when conducting sustainability evaluations of important reservoir catchment areas in Taiwan.Chen et al. [55] also established a set of combined indicators for dam safety, reservoir management, water ecological system stability, reservoir benefits, and reservoir sustainability to assess the level of sustainable reservoir development.Petersson [56] investigated reservoir and dam sustainability issues and used changes in flow to rate reservoir and dam sustainability levels.Kumambala et al. [57] proposed the use of a hydrologic indicator, a human health indicator, and an environment indicator for water resource utilization in Malawi.Chu et al. [58] studied the impact on overall reservoir performance according to changes in the characteristics of local ecology, water quality, hydrology, and habitat environment.Recent attempts to solve the abovementioned environmental problems [12] and maintain reservoir functionality (i.e., power generation, flood control, habitat protection, water allocation, water quality, and sediment control) include using multi-objective reservoir operation modeling where water utilization and flood control targets are constrained by ecological flow or water quality.
Proposed SCI indicators in the reviewed literature have only been used for single purpose analyses such as water quality, water quantity, and flooding or sediment control.A multi-objective reservoir evaluation, however, could be used to solve environmental and ecological problems while ensuring the sustainable functions of reservoirs [11,12,18,41].The main focus of this study is to incorporate sustainability concepts when making planning considerations and during operations of multi-purpose reservoirs with the aim of improving reservoir sustainability for reservoirs in Taiwan.Since sustainability concepts apply to the environment, economy, and society, it is therefore necessary to establish a multi-factor indicator-based evaluation system to determine reservoir sustainability levels.For Taiwan reservoir management, this study proposes an indicator system that includes calculations and standardizations employed for each indicator, weighting factors for each indicator, integrated indicator analysis, and a reservoir sustainability assessment.The proposed system was applied in a case study to demonstrate its functional application for sediment management in the reservoir planning process.

Materials and Methods
In order to identify indicators, 19 experts and scholars within hydraulic engineering, hydrology, and ecohydrology fields discussed 12 indicators from three different perspectives, including that of academia, non-governmental organizations, and government authorities in Taiwan (i.e., Water Resources Planning Institute, Water Resources Agency, Ministry of Economic Affairs [59]) (Table S1 in Supplementary Material), by way of research projects (i.e., Water Resources Planning Institute, Water Resources Agency, Ministry of Economic Affairs [60,61]) and a case study using a national reservoir calculus (i.e., Water Resources Planning Institute, Water Resources Agency, Ministry of Economic Affairs [62]).With the indicators and their restrictions defined, index calculations, rating schemes, and other operational procedures were integrated to formulate a reservoir SCI system.The evaluation system framework (Figure 1) incorporates "sustainable management" concepts within the planning processes so that newly constructed reservoirs can meet sustainability requirements.To evaluate the sustainability of existing reservoirs and, if necessary, to facilitate improvements so as to enhance their sustainability, the system consists of 12 indicators in six categories related to social justice, environmental protection, and economic development.The code I 11 could be considered as intergenerational equity.This study considered the code by benefit sharing.

Indicators for Reservoir Function Evaluation
Indicators for assessing reservoir function were grouped into six categories.Index assessments and standards must also be considered in addition to establishing analytical tools, since each index requires extensive computational and analytical work and too many indicators increase the system complexity and user accessibility.Therefore, this study selected a total of 12 indicators from published literature, taking into consideration Taiwan's unique hydrological and geological conditions as well as the specifications of indicator applicability.The following sections further describe the indicators listed in Table 1.
1  The flood control system includes the reservoir itself and downstream river channels, the primary purpose of which is to reduce flood flow for reservoir and river channel safety.The indicators selected for this category are a flood reduction ratio (code: I 1 ) and a safe flood flow ratio (I 2 ).The former is defined as the ratio of peak inflow reduction (Q pi − Q po ) to peak inflow (Q pi ), while the latter is defined as the ratio of peak outflow from the reservoir (Q po ) to the designed maximum flow of the downstream river channel (Q r ).The indicator I 1 represents the reservoir's flood reduction capability.The higher the I 1 value, the better the flood reduction capability of the reservoir.However, since reservoir functions include water supply, I 1 is subject to a number of restrictions as well.Q r indicates whether downstream river channels have the capacity needed to accommodate peak outflow from the reservoir.Hence, a smaller I 2 suggests a decreased potential of disastrous downstream flooding.

Sediment Management (C 2 )
Two indicators selected for this category are the effective reservoir capacity ratio (I 3 ) and the adequacy of sediment supply to downstream and coastal areas (I 4 ).The indicator I 3 is defined as the ratio of effective reservoir capacity at a specific point in time to the capacity at the time of its completion (or renovation).Effective reservoir capacity at a given time is equal to the initial effective reservoir capacity (V 0 ) minus the volume of sediment (V s ) accumulated from years before.The indicator I 4 is defined as the ratio of annual sediment supply to downstream areas after completion of the reservoir (Q d ) to the annual sediment supply prior to completion (Q s ).Q s can be acquired from historical data, while Q d can be estimated by simulation or determined as the reservoir becomes operational.
Changes in I 3 reflect the effects of desilting and sedimentation, and therefore can be used as a basis to decide whether to take desilting or soil conservation measures in upper catchment areas.The indicator I 4 reflects the downstream river channel and coastal area, as well as river ecology, and protection needs, and therefore can be used as a basis for adjusting reservoir desilting efforts.

Water Resource Allocation (C 3 )
The water supply ratio and the water allocation ratio which are included in water resource allocation (C 3 ) were simulated and analyzed based on annual records of daily water supply and demand, and have been set by the reliability standard frequently adopted in Taiwan (SI = 1.0 was adopted in the cases in this research).Therefore, no further descriptions were considered or calculations were performed regarding this index.The indicators selected for this category are water supply stability (I 5 ) and adequacy of water resource allocation (I 6 ).The indicator I 5 is defined as the ratio of the amount of water the reservoir actually supplied (SA) to the planned water supply (S), and is used to determine whether the reservoir is capable of providing the planned water supply.The indicator I 6 is defined as the ratio of actual water allocated to various targets (SA) or planned water supply (S) to the corresponding demand of the targets (D).It is used to evaluate the adequacy of water allocation, and to determine if the reservoir must adjust its water supply to various targets.For reservoirs currently in planning phases, the SA value is obtained by simulating reservoir operations to obtain a planned water supply value versus the actual supply.In terms of stability or adequacy, any water resource supplied that exceeds demand is considered poor water resource allocation.

River Ecology (C 4 )
This study selected an ecological base flow ratio (I 7 ) and a measurement of change in river habitat area (I 8 ) as the two indicators for the conservation of river ecology.The indicator I 7 is defined as the ratio of the ecological base flow released from the reservoir (Q sb ) to the dry season flow with a 95% probability of exceedance (Q 95 ).It is used to evaluate the effects of the base flow on downstream ecology and to act as a reference to determine if discharge adjustments are needed during the dry season.The indicator I 8 is defined as the percentage of change in WUA before and after reservoir construction.It is used to assess whether the reservoir discharge meets ecological requirements, and as a reference when reviewing river ecology conservation plans.
The WUA represents the equivalent area suitable for the survival of indicator organisms, and can be calculated using the following equation: in which, F[ . . .] is the combined suitability index of habitat mesh i and is used as the weighting factor; f (V i ), f (D i ), and f (d i ) are, respectively, the suitability indices of flow velocity V i , water depth D i , and bed sediment particle size d i of mesh i; A i is the habitat area of mesh i; and n represents the total number of meshes in the habitat.The suitability index can be obtained from the habitat suitability curve [63].

Water Quality (C 5 )
Good water quality is a key factor for ensuring reservoir sustainability.Reservoir water quality directly affects water purification and treatment costs, drinking water quality, and public hygiene.Water quality maintenance indicators can be divided into two types.The first type, used to assess the degree of water body pollution, includes a river pollution index (RPI) and a water quality index (WQI).The second type, used to assess the degree of eutrophication in a reservoir water body, is the total phosphorus (TP) under normal conditions.
RPI is composed of the four parameters Dissolved Oxygen (DO), Biochemical Oxygen Demand (BOD), Suspended Solids (SS), and Ammoniacal Nitrogen (NH 3 -N), and graded according to their concentrations.WQI includes the seven parameters DO, BOD, pH, NH 3 -N, E. coli colonies, SS, and TP, which are collectively referred to as WQI 7 .WQI 7 is calculated using the following equation [64]: in which, w i is the weighting factor of each water quality parameter, based on the outcomes of expert surveys [64]; q i are the grades of each water quality parameter, and can be obtained through model simulations or on-site surveys.
In this study, WQI 7 and TP are selected for a water body pollution indicator (I 9 ) and for a eutrophication indicator (I 10 ), respectively.I 9 can be calculated upon obtaining the seven water quality parameters included in WQI 7 .Results are then used to categorize the water quality and as a reference for assessing water body pollution status.For the indicator I 10 , the TP concentration can be obtained from field survey or from a Vollenweider model [65] simulation.Parameters required of this model include the total volume of water in reservoir inflow of TP, settling rate of TP, area of inundation, and outflow velocity.

Benefit and Fairness (C 6 )
As mentioned earlier, the sustainable development concept requires satisfying the needs of the present generation without threatening the needs of future generations [66].Considering that the main objective of this research is to understand and analyze the future economic benefit of water and the reduction of actual water supply, "benefit" and "fairness" were taken as indices to demonstrate how the indices change over time from the present to future, and to further estimate and detect variances between time periods that will affect future generations [50,66].Hence, reservoir planning or renovation must consider the issue of water supply benefits and fairness.Indicators for benefits and fairness include cost-benefit ratios, economic benefits, benefit indices, Gini coefficients, a fairness index, and public acceptance.The Gini coefficient is the ratio of area A to the sum of area B and area A [28], and can be used for water allocation equity [28].In this study, industrial revenue derived from water will change over time, and shortages of water supply will worsen.Since the two evaluation results (e.g., benefit and fairness) are both categorized as changes affecting generations, both are considered to be in the same category [50,66].
After considering data accessibility and indicator representativeness, this study selected the mean annual production value per unit water use (GWR i ) as a benefit indicator (I 11 ) and the Gini coefficient (G ini ) as a fairness indicator (I 12 ).In other words, I 11 = ∑ n i GWR i /n while I 12 = G ini .GWR i is calculated by using the following equation: In the equation above, GW ti is the production value per unit water use (NT$/m 3 ) in the year t i ; T is the time interval (five years); i is the index of year census taken; GWR i is the mean annual change of GW ti within the interval of T years.Hence, ΣGWR i /n is the mean over the total number of n intervals.
The Gini coefficient is defined by the following equation [48]: in which, X i is the dimensionless cumulated time, Y i is the dimensionless cumulated volume of water supplied, i = 1, 2, . . ., n, with n as the total number of time segments.As shown in Figure 2, (X i − X i−1 )(Y i + Y i−1 ) = 2∆B i is equal to the increase in the dimensionless water supply volume from the time X i−1 to X i .Hence, the second part on the right side of Equation (4) becomes 2 Under such conditions, the Y i = f (X i ) curve shown in Figure 2 coincides with the 45 • -line, and all the increments ∆Y i are the same for each equal segment of ∆X i .This is why the 45 • -line is referred to as the line of equality or complete fairness.The larger the deviation of the Y i = f (X i ) curve from the 45 • -line, the less fairness.Hence, the area A in Figure 2 between the curve and the 45 • -line represents the degree of inequality.The curve represents real water allocation [52].Therefore, the area A (dark color) represent differences between the perfect water allocation and real water allocation during various cumulated time periods.In the worst case of extreme inequality, the Y i = f (X i ) curve goes along the horizontal axis and suddenly rises to 1 at X i = 1.
The equity of water allocation can be calculated by the equitable sharing of the used water quantity for economic benefits [28].Evaluation outcomes of the benefit and fairness indicators could be used to assess current and future water supply situations.When poor benefit or fairness is observed, managerial actions to improve the situation can be taken.
In the equation above, GWti is the production value per unit water use (NT$/m 3 ) in the year ti; T is the time interval (five years); i is the index of year census taken; GWRi is the mean annual change of GWti within the interval of T years.Hence, ΣGWRi/n is the mean over the total number of n intervals.
The Gini coefficient is defined by the following equation [48]: in which, Xi is the dimensionless cumulated time, Yi is the dimensionless cumulated volume of water supplied, i = 1, 2, …, n, with n as the total number of time segments.As shown in Figure 2, (  −  − )(  +  − ) = ∆  is equal to the increase in the dimensionless water supply volume from the time Xi−1 to Xi.Hence, the second part on the right side of Equation ( 4) becomes  ∑ ∆ . If  ∑ ∆    = , then Gini = 0.Under such conditions, the Yi = f(Xi) curve shown in Figure 2 coincides with the 45°-line, and all the increments ΔYi are the same for each equal segment of ΔXi.This is why the 45°-line is referred to as the line of equality or complete fairness.The larger the deviation of the Yi = f(Xi) curve from the 45°-line, the less fairness.Hence, the area A in Figure 2 between the curve and the 45°-line represents the degree of inequality.The curve represents real water allocation [52].Therefore, the area A (dark color) represent differences between the perfect water allocation and real water allocation during various cumulated time periods.In the worst case of extreme inequality, the Yi = f(Xi) curve goes along the horizontal axis and suddenly rises to 1 at Xi = 1.
The equity of water allocation can be calculated by the equitable sharing of the used water quantity for economic benefits [28].Evaluation outcomes of the benefit and fairness indicators could be used to assess current and future water supply situations.When poor benefit or fairness is observed, managerial actions to improve the situation can be taken.

Reservoir Sustainability Rating
Literature points to isometric rating schemes, which are used to identify the relative levels of different indicators as the most used evaluation indicators.In this study, each indicator is rated by using a four-level scheme of Excellent, Good, Fair, and Poor with index scores (Si) of 4, 3, 2, and 1,

Reservoir Sustainability Rating
Literature points to isometric rating schemes, which are used to identify the relative levels of different indicators as the most used evaluation indicators.In this study, each indicator is rated by using a four-level scheme of Excellent, Good, Fair, and Poor with index scores (S i ) of 4, 3, 2, and 1, respectively.This scheme helps resolve difficulties posed by differences in index units or values when calculating overall assessment ratings.For flood control capability, I 1 represents the reservoir's flood reduction capacity.Hence, I 1 ≥ 30% is regarded as Excellent, indicating that the reservoir can perform robust flood reduction functions; I 1 = 21~30% is rated as Good, meaning that the reservoir is still capable of providing good flood reduction; I 1 = 11~20% is rated as Fair, meaning that the reservoir may be losing its capability to provide adequate flood reduction; I 1 ≤ 10% is rated as Poor, indicating that the reservoir is losing its flood reduction capability.Rating schemes for the other indicators follow the same principle.
To assess the overall reservoir sustainability, each indicator weighting factor was analyzed with the analytic network process (ANP) (i.e., Water Resources Planning Institute, Water Resources Agency, Ministry of Economic Affairs [68]).ANP survey questionnaires were distributed to three groups of experts in academic institutions, non-governmental organizations (NGOs), and public agencies with hydrology, hydraulic engineering, and ecohydrology backgrounds.The weighting factor for each category of indicators was then determined by averaging the returned questionnaire responses (Section S3 in Supplementary Material).To simplify the ANP surveys, the two indicators in each category were given the same weighting.Table 3 lists the weighting factors of the 12 indicators (δ i , i = 1, 2, . . ., 12).The sum of the 12 weighting factors Σδ i is equal to 1.For the sustainability evaluation, index scores of each indicator are multiplied by their corresponding weighting factor and then summed.Since the sums range from 1 to 4, the sustainability rating was divided into four levels with index scores as shown in Table 4. Reservoirs rated as "Excellent" are in very good condition with high sustainability, and should be monitored for any change in their sustainability; reservoirs rated as "Good" are considered sustainable with some indices degrading, and measures should be taken to upgrade them; reservoirs rated as "Average" are those with deteriorating sustainability, and in-depth investigations should be carried out to formulate improvement strategies; those rated as "Poor" should undergo drastic improvements to restore their sustainability.Operational procedures for applying evaluation indicators in reservoir planning is shown in Figure 3. Required data are first collected while compatible models are established for each of the indicators.ANP is adopted as the tool for decision analysis in Figure 3.The collected data and models were then used to compute indices while index scores (S i ) and weighting factors were obtained from Tables 2 and 3 (δ i ), respectively.Finally, the total score of reservoir sustainability was calculated as Σ(δ i S i ) which is the sustainability confidence index (SCI = Σ(S i ) × (δ i )).To evaluate changes in reservoir sustainability due to environmental, economic, or social factors, analyses are carried out to identify key factors affecting the reservoir's sustainability (such as flood control, desilting, or water quality).Following this, relevant indices and weighting factors are revised, and sustainability is re-evaluated accordingly.Operational procedures for applying evaluation indicators in reservoir planning is shown in Figure 3. Required data are first collected while compatible models are established for each of the indicators.ANP is adopted as the tool for decision analysis in Figure 3.The collected data and models were then used to compute indices while index scores (Si) and weighting factors were obtained from Tables 2 and 3 (δi), respectively.Finally, the total score of reservoir sustainability was calculated as Σ(δiSi) which is the sustainability confidence index (SCI = Σ(Si) × (δi)).To evaluate changes in reservoir sustainability due to environmental, economic, or social factors, analyses are carried out to identify key factors affecting the reservoir's sustainability (such as flood control, desilting, or water quality).Following this, relevant indices and weighting factors are revised, and sustainability is reevaluated accordingly.

Study Area
This study selected the Shiwen Reservoir, currently in its planning phase, to demonstrate an application of the proposed evaluation system.The following section will briefly describe the

Study Area
This study selected the Shiwen Reservoir, currently in its planning phase, to demonstrate an application of the proposed evaluation system.The following section will briefly describe the background of the reservoir, calculation and analysis procedures of each index (S i & δ i ), and evaluation of the sustainability confidence index (SCI) where SCI = Σ(S i ) × (δ i ).Details are available in relevant reports [67] (Tables S4-S6 in Supplementary Material).Shiwen Creek, in Taiwan's Pingtung County, is the site of the Shiwen Reservoir, as shown in Figure 4.The reservoir has a catchment area of 65.53 km 2 and an effective capacity of 65 × 10 6 m 3 and a planned daily water supply of 2 × 10 5 m 3 , which satisfies Pingtung area water requirements until the target year (2031).The Shiwen dam projected height is 103 m, with a normal full pool level at elevation (EL).170 m, a dam crest at EL. 178 m, and a total capacity of 70.13 × 10 6 m 3 .Major flood discharge facilities include a spillway and a sluiceway with designed discharges of 1439 m 3 /s and 947 m 3 /s, respectively.The spillway has two bays, 21 m wide each, and its crest is at EL. 170 m.The sluiceway has three bays, 4 m in diameter each, and its outlet is located at the bottom of the dam.

Sustainability Analysis
Table 5 lists the calculated results of indices and SCIs.Though marginally rated as "Excellent", the sum of sustainability confidence indices (ΣSCI) for the case excluding desilting facilities is equal to 3.3, while indices I3 and I4 are rated as "Poor" and "Fair", respectively.The inclusion of desilting facilities raises indices I3 and I4 to the levels of "Good" and "Excellent", respectively, and consequently, the ΣSCI increases from 3.3 points to 3.7 points, making Shiwen Reservoir even more sustainable.The index scores and SI rating ranges are converted into a radar chart, shown in Figure 5, to facilitate the planning and decision-making processes.The following describes the results of the analysis for each indicator.

Sustainability Analysis
Table 5 lists the calculated results of indices and SCIs.Though marginally rated as "Excellent", the sum of sustainability confidence indices (ΣSCI) for the case excluding desilting facilities is equal to 3.3, while indices I 3 and I 4 are rated as "Poor" and "Fair", respectively.The inclusion of desilting facilities raises indices I 3 and I 4 to the levels of "Good" and "Excellent", respectively, and consequently, the ΣSCI increases from 3.3 points to 3.7 points, making Shiwen Reservoir even more sustainable.The index scores and SI rating ranges are converted into a radar chart, shown in Figure 5, to facilitate the planning and decision-making processes.The following describes the results of the analysis for each indicator.

Flood Control Capability
The simulated peak inflow and outflow for a 100-year return for the Shiwen Reservoir are 1136 m 3 /s and 769 m 3 /s, respectively.Hence, the flood reduction ratio I 1 = 32.3% and is rated as "Excellent".The designed flood flows for a 100-year return period for river channels downstream are 1180 m 3 /s for Shiwen Creek and 1460 m 3 /s for Shuaimang River.Therefore, the safe flood flow ratios are 65.2% for Shiwen Creek and 52.7% (=I 2 ) for Shuaimang River.Hence, the Shiwen Reservoir is rated as "Excellent" for its flood control functionality.

Sediment Management
The effective capacity ratio I 3 for the Shiwen Reservoir excluding desilting operations is 49.61% and is rated as "Poor".When including desilting operations, I 3 becomes 80.11% and is rated as "Good".The index I 4 for adequacy of sediment supply to downstream channels and coastal areas is 13.5%, and rated as "Fair" when excluding desilting operations.With desilting operations, I 4 becomes 65.25% and is rated as "Excellent".

Water Resource Allocation
The simulated mean annual water supply value (70.2 × 10 6 m 3 ) results in a water supply stability index (I 5 ) of 96.16%, which differs from the "Excellent" rating of the planned estimate value (73.0 × 10 6 m 3 ).Increased public water consumption in the Pingtung Area by the target year (2031) is projected to exceed the current supply by 87,000 m 3 /day.Since the Shiwen Reservoir is capable of supplying an additional 120,000 m 3 /day, the adequacy of the water allocation index is 12/8.7 × 100% = 137.93%.Water resources for agricultural purposes are the priority, and its allocation index is assumed to be 100%.Therefore, the simple mean of the two, (100 + 137.93)/2 = 118.97%, is taken as the overall allocation index (I 6 ) and is rated as "Good".

River Ecology
The planned ecology baseflow for Shiwen Reservoir is 0.086 m 3 /s, while the dry season flow with a 95% probability of exceedance is 0.0336 m 3 /s.Therefore, the ecological baseflow ratio (I 7 ) is 2.56 and is rated as "Good".The WUA of downstream river habitats before and after completion of the reservoir are calculated by using a 2D river simulation model and are reduced from 613 m 2 during the wet season and 321 m 2 during the dry season to 245 m 2 .The percentage of reduction (I 8 ) ranges from 60% to 23%, with a simple mean of 41.5%, and is rated as "Fair".

Water Quality
Sampling stations in the Shiwen Reservoir catchment area collected water quality data.Data analysis shows that index I 9 (=WQI 7 ) ranges from 77 to 92, with a simple mean of 87, and is rated as "Excellent".Vollenweider Model simulation results yield a total phosphorus content (TP) of 4.35 to 15.47 g/L.Average TP of the reservoir water body is 11 g/L.Therefore, the eutrophication index I 10 is rated as "Excellent".

Benefits and Fairness
The Shiwen Reservoir is expected to be a major source of water for the Pingtung Area.The area's total industrial and commercial production for 2006 and 2011, as well as the public water supply, provided the basis for calculating unit water use and its effect on production value for these two years.The results are 825 and 1301 $TWD/m 3 for unit water use and production effect, respectively.The Directorate-General of Budget, Accounting, and Statistics of the Taiwan Executive Yuan conducts industrial, commercial, and service censuses once every five years, while the Water Resources Agency has compiled annual water use data for each administrative district since 2001.A time step that correlates with the five-year census frequency is used in the first dataset (n = 1) from 2006 with T = 5.The mean annual change in the effect of unit water use on production value is 2 × (1301 − 825)/((1301 + 825) × 5) × 100% = 9%.The benefits index I 11 is thus 9.0% and is rated as "Excellent".
Daily water supply volume estimates of reservoir operations are simulated using 46 years of historical hydrological data from 1962 to 2007 (Table S7 in Supplementary).The simulated results are then used to calculate the cumulated daily water supply ratios (simulated water supply/planned water supply).The yearly minimum water supply ratios are then ranked in ascending order.A normalized cumulative curve is then plotted and compared with the 45-degree line equality curve.The water supply inequality level, which is represented by the area between the two lines, is then calculated.The result obtained is the Gini coefficient, I 12 = 0.2, and is rated as "Excellent".

Discussion
In reservoir operations, there is a tendency to increase the maximum water level, and therefore capacity, to meet downstream water demands.Yet, the ability to meet downstream demands is a function of a range of factors including infrastructure, operations, and water supply volume/timing.Increasing maximum water levels does not necessarily increase capacity to meet downstream demands since, in doing so, flood storage capacity is reduced and reservoir flood control functionality is possibly affected, impacting the overall reservoir operations.Moreover, reservoir desilting can lead to excessive turbidity, and result in water supply issues and downstream ecological impacts.Accordingly, integrated reservoir management aims to maximize and maintain reservoir functions such as power generation, flood control, habitat protection, and water allocation, as well as water quality and sediment controls [11,18].During flooding periods, however, inundation mitigation of downstream areas might counter the natural flow of supply for ecological needs, which is deprioritized to increase upstream flood control [9].Thus, this study applied an evaluation system for sustaining reservoir functions to a case study in Taiwan by taking into account reservoir desilting and catchment area sediment yield reduction, as well as considering water quantity and quality, habitat protection, and operational benefit and fairness.
Numerous recommendations for using statistical approaches in environmental flow studies have been reported.For example, Morrison et al. [69] used IHA to assess the impacts of environmental flow on reservoir and recreational operations.Nikghalb et al. [44] as well as Alexandre et al. [70] used Q 95 to evaluate impacts of flows on movement patterns and habitat.They also determined ecological flow using WUA, and compared these results with Q 95 results.Aguilar and Polo [45], Li et al. [40], Pasha et al. [41], Yao et al. [42], and Xu et al. [15] used WUA as an indicator with Q 95 to represent the ecological condition of their study reservoir.However, Nikghalb et al. [44] indicated that different criteria are needed to determine the environmental flow requirement from the WUA−discharge relationship.This study used WUA change ratio values to avoid the above issue, but this value should be thoroughly evaluated when used in other applications.For further tests, the IHA change ratio should be included when evaluating the results of WUA and Q 95 .
As recommended by recent studies in relevant fields [2,3,11,12,18,29,30,38,39,[41][42][43], the proposed indicators were established to evaluate environmental, water resource-related, and ecological sustainability in Taiwan, with the exception of fairness which is represented by the benefit index and Gini Coefficient.Therefore, in contrast with recent studies that have evaluated multi-purpose reservoirs, the proposed evaluation system takes the benefit index and Gini coefficient into account to ensure future water supply and reservoir sustainability.Although the Gini coefficient has been applied in the evaluation of the characteristics and applicability in river flow regimes under future climate change, it can also be used to evaluate the distribution of precipitation of wet-day frequency and decreasing wet-day, to assess effects of ecological simplification on flow variability, as well as to find equity in the water allocation strategies in a river basin.Moreover, ANP was used to obtain weighted factors for each proposed indicator to deal with trade-offs within the six categories, namely, flood control, sediment management, water resource allocation, river ecology, water quality, and benefit and fairness.The proposed approach can be universally applied to Taiwan reservoirs if more case studies are conducted using this study's methodology.
This study applied the proposed evaluation system to a reservoir in its planning phase to achieve the reservoir's flood control capacity and sediment management aims.The evaluation SCI value for the reservoir in this case study is rated as "Good" and "Excellent" in sustainability.Additionally, the ANP results indicated that the flood control capacity and sediment management are the first and second most important indicators for the reservoir.If desilting operations had been conducted, the SCI values would have increased from 3.3 to 3.7, and a SCI value rated as "Excellent" achieved.However, reservoirs in Taiwan must be equipped with adequate desilting facilities due to the unique hydrological and geological conditions characterized by the rapid rise and recession of flood flow, and the potential for extremely severe landslides in the catchment area.Reservoir sustainability evaluations must consider the fact that environmental, economic, and social factors play an important role in the overall planning of sustainable water resource systems.

Conclusions
This study proposed a reservoir evaluation system that includes 12 indicators grouped into six categories of flood control, sediment management, water resource allocation, river ecology, water quality, as well as benefit and fairness.The system was successfully applied in a case study.In this study, a framework, concepts, processes, index calculations, and rating schemes were presented along with a proposed reservoir sustainability evaluation procedure.This procedure could supplement current reservoir planning and decision-making.A four-level scheme of Excellent, Good, Fair, and Poor with index scores of 4, 3, 2, and 1 points, respectively, was established to rate the index obtained for each indicator.This rating scheme can be used to assess reservoir sustainability and for planning revisions if needed.A weighting factor for each evaluation indicator was determined using the ANP method.Index scores of each indicator are multiplied by its weighting factor and then summed to yield the reservoir sustainability confidence index.
Since reservoir sustainability may involve historical and environmental factors, the suitability of its evaluation indicators requires thorough examination.Lack of relevant data, which is often the case, may require additional investigation and analyses.Indicator weighting factors, used for the planning of sustainable reservoirs, may be adjusted, if necessary, given that different reservoirs face varied issues.This study recommends further studies on this issue to be conducted, since the specific indicators for this study may not be generally applicable to all reservoirs in Taiwan.Indicators need to be more soundly justified and tested prior to further applications.Accordingly, feedback on the applications of the proposed system will improve its models.Future work may include an evaluation of the sustainability of all of Taiwan's reservoirs using the proposed system, so as to a form a decision-making support platform.The proposed procedure may be integrated with a Decision Support System for sustainable reservoir management in further study.It is also suggested that the impact of climate change on reservoir sustainability is investigated to identify any required adjustment plans.

Figure 2 .
Figure 2. Definition of the Gini coefficient.

Figure 2 .
Figure 2. Definition of the Gini coefficient.

Figure 3 .
Figure 3. Reservoir sustainability evaluation system flowchart for planning and decision-making support.

Figure 3 .
Figure 3. Reservoir sustainability evaluation system flowchart for planning and decision-making support.
reservoir, calculation and analysis procedures of each index (Si & δi), and evaluation of the sustainability confidence index (SCI) where SCI = Σ(Si) × (δi).Details are available in relevant reports[67] (Tables S4-S6 in Supplementary Material).Shiwen Creek, in Taiwan's Pingtung County, is the site of the Shiwen Reservoir, as shown in Figure4.The reservoir has a catchment area of 65.53 km 2 and an effective capacity of 65 × 10 6 m 3 and a planned daily water supply of 2 × 10 5 m 3 , which satisfies Pingtung area water requirements until the target year (2031).The Shiwen dam projected height is 103 m, with a normal full pool level at elevation (EL).170 m, a dam crest at EL. 178 m, and a total capacity of 70.13 × 10 6 m 3 .Major flood discharge facilities include a spillway and a sluiceway with designed discharges of 1439 m 3 /s and 947 m 3 /s, respectively.The spillway has two bays, 21 m wide each, and its crest is at EL. 170 m.The sluiceway has three bays, 4 m in diameter each, and its outlet is located at the bottom of the dam.

Figure 5 .
Figure 5. Radar chart of indices (a) and bar chart of SCI (b) for the Shiwen Reservoir.I 1 : Flood reduction ratio; I 2 : Safe flooding flow ratio; I 3 : Effective reservoir capacity ratio; I 4 : Sediment supply ratio; I 5 : Water supply ratio; I 6 : Water allocation ratio; I 7 : Ecological baseflow ratio; I 8 : Change in WUA; I 9 : Water pollution index; I 10 : Eutrophication index; I 11 : Benefits index; I 12 : Fairness index.

Table 1 .
A summary of various evaluation indicators for reservoir evaluation.
lists the ranges of index values, corresponding levels, and scores of the indicators selected in this study.The indicator ranges in Table 2 (upper and lower limit) were taken from relevant published literature, and further revised in this study for the national reservoir

Table 2 .
Rating scheme for reservoir sustainability indicators.

Table 3 .
Weighting factors of indicators.

Table 4 .
Reservoir sustainability rating and grade.

Table 4 .
Reservoir sustainability rating and grade.

Table 5 .
Indicator weighting factors and Shiwen Reservoir sustainability rating during desilting operations.

Table 5 .
Indicator weighting factors and Shiwen Reservoir sustainability rating during desilting operations.
* Results for the absence of desilting operations are provided in parentheses.ΣSCI = 3.7 (3.3).