Next Article in Journal
Temporal Variability of Suspended-Solids Concentration in the Estuarine Channel of Patos Lagoon, Southern Brazil
Previous Article in Journal
Application of LSTM Networks for Water Demand Prediction in Optimal Pump Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sustainable Surface Water Storage Development Pathways and Acceptable Limits for River Basins

1
International Water Management Institute (IWMI), P.O. Box 2075, Colombo 10120, Sri Lanka
2
UNU Institute for Water Environment and Health, 204–175, Longwood Road South, Hamilton, ON L8P 0A1, Canada
3
Faculty of Engineering Technology, Open University of Sri Lanka, P. O. Box 21, Nawala, Nugegoda 10250, Sri Lanka
*
Author to whom correspondence should be addressed.
Water 2021, 13(5), 645; https://doi.org/10.3390/w13050645
Submission received: 22 January 2021 / Revised: 17 February 2021 / Accepted: 18 February 2021 / Published: 28 February 2021
(This article belongs to the Section Water Resources Management, Policy and Governance)

Abstract

:
This paper addresses the questions of acceptable upper limits for storage development and how best to deploy storage capacity in the long-term planning of built surface water storage in river basins. Storage-yield curves are used to establish sustainable storage development pathways and limits for a basin under a range of environmental flow release scenarios. Optimal storage distribution at a sub-basin level, which complies with an identified storage development pathway, can also be estimated. Two new indices are introduced—Water Supply Sustainability and Environmental Flow Sustainability—to help decide which pathways and management strategies are the most appropriate for a basin. Average pathways and conservative and maximum storage limits are illustrated for two example basins. Conservative and maximum withdrawal limits from storage are in the range of 45–50% and 60–65% of the mean annual runoff. The approach can compare the current level of basin storage with an identified pathway and indicate which parts of a basin are over- or under-exploited. A global storage–yield–reliability relationship may also be developed using statistics of annual basin precipitation to facilitate water resource planning in ungauged basins.

1. Introduction

Storing water helps in stabilizing variable water supplies, producing hydropower, and mitigating floods. Water storage options include small and large reservoirs, ponds and tanks, aquifers, soil moisture, and wetlands [1]. Reservoirs, ponds, tanks, and wetlands are generally regarded as surface storage, and the others as subsurface storage. Surface storage has been beneficially exploited for thousands of years. Examples span from dams constructed by Sumerians in the Euphrates and Tigris Rivers in 6500 BC [2] and the ancient tank irrigation systems that flourished in Sri Lanka from around 300 BC [3,4] to the mega dams of the twentieth century in different parts of the world [5,6,7]. However, surface storage is known to create adverse ecological and social consequences, too.
At least 43% of all global water abstractions for irrigation are supplied by surface storage (reservoirs with a storage capacity larger than 0.1 km3). The remainder is split between groundwater supply (22%), aqueduct transfers (8%), and other unspecified sources (27%) [8]. Hydropower accounts for close to 17% of global electricity generation, mostly produced by large dams [9,10]. Reservoirs support a number of other non-consumptive activities such as flood risk reduction, recreation, fishing, and small-scale livelihoods. Apart from the tangible benefits, in certain countries, large-scale water storage infrastructure serves as a strategic tool for promoting domestic peace and stability [11,12] and a symbol of national identity and patriotism [13]. Large water storage infrastructure combined with transboundary water transfers also function as a foreign policy tool to enhance influence at the international level [11,14,15]).
On the other hand, only 23% of rivers longer than 1000 km remain unfragmented by dams [16]. Dams regulate river flow [17,18] and alter downstream flow regimes, which negatively impacts aquatic ecosystems [16,19,20,21] and in turn fisheries, livelihoods, and cultural and recreational activities [22,23,24]. Large dams often displace a significant number of people and affect those living downstream, causing losses to property, livelihoods, and incomes [25,26,27]. The number of people displaced by dams in the last century is at least 80 million [28]. Furthermore, construction of large dams frequently creates conflicts between competing expectations of different groups of stakeholders [15,29].
The current global tally of large dams (15 m high or greater, and between 5 m and 15 m impounding more than 3 million m3) is close to 60,000 and growing [30]. An indicative figure for the number of smaller dams below the International Commission on Large Dams (ICOLD) benchmark is 800,000 and also growing [31]. While surface water storage capacity continues to grow globally, its development remains largely ad-hoc, i.e., driven by the needs of the immediate future, the availability of suitable sites, and short-term political interests. The question of what is an acceptable upper limit for surface storage development in a river basin remains unanswered. Concepts of closed basins [32], sustainability boundaries [33], and water footprints [34] prescribe limits to withdrawals from available river flow but do not suggest an acceptable or optimal storage limit for a particular river basin. Related questions such as how best to deploy the maximum acceptable storage capacity in a river basin in terms of the size, number, and location of reservoirs that maximizes benefits and minimizes environmental and social costs have not been well studied. A few studies, such as [35] and [36], have dealt with the latter question, but not in the context of an upper limit to surface storage development.
There is a significant body of literature on the optimal operation of particular reservoirs or reservoir systems to meet human and ecological water requirements rather than considering them in a basin planning context. Examples include studies [37,38,39,40,41], which seek to optimize operations of existing reservoirs to meet multiple objectives ranging from water supply, flood control, energy generation, and environmental flow. Although there are numerous other examples, a lengthy discussion is beyond the scope of this paper.
This paper moves beyond traditional optimization of existing reservoir operations towards a long-term, multi-decadal vision for surface storage development in a river basin. It recognizes that an acceptable upper limit to storage development is a compromise between storage benefits and ecosystem and social costs [42] and that tools to support the decision-making process for structured storage development are needed. The paper examines how best to determine this compromise between storage benefits and ecosystem costs when developing storage for supplying water for various human needs. It investigates the behavior of a set of indicators measuring water supply benefits and environmental flow releases at different storage capacities, numbers, modes of siting, and spatial distribution of reservoirs (different scenarios of reservoir system configurations).
The paper illustrates how an overarching sustainable storage development pathway and acceptable limits to surface storage development may be identified for a river basin to limit withdrawals from reservoir storage to a sustainable level. It also illustrates the optimal manner in which storage capacity and sustainable withdrawal levels may be distributed among sub-basins so that a river basin as a whole complies with the overarching sustainable storage development pathway, even when the reservoir configuration changes. Finally, it discusses how the results may be used in practice to identify drainage regions having higher surface storage potential, to plan for a more beneficial distribution of storage capacity, and to allocate environmental flows (EF).

2. Materials and Methods

2.1. Example Case Study Basins

Approaches presented here can be applied to any river basin. For illustrative purposes, this paper uses two case study basins in Sri Lanka: the Malwatu Oya in the Dry Zone and the Kalu Ganga in the Wet Zone (Figure 1). The Dry Zone (stretching to the north and the east of the island) has a mean annual precipitation (MAP) of less than 1750 mm and a pronounced dry season with high evaporation rates from May to September. The Wet Zone lies in the southwestern part of the country and receives a MAP of over 2500 mm, which is normally spread throughout the year. This zoning is used in Sri Lanka, which is largely humid overall, hence, the characterizations “wet” and “dry” are relative to the national context. Apart from different climatic zones, the two basins contrast in terms of terrain, existing storage infrastructure, basin area and shape, type of land cover and land use, and dominant water issues.
The Malwatu Oya basin has a drainage area of 3338 km2, a MAP of 1300 mm, and a natural mean annual runoff (MAR) of 0.79 km3 (coefficient of variation 0.47). Its elevation varies from 0 to 740 m amsl. Its gently undulating terrain is dotted with small tank cascade systems that originated in the Third Century BC (Figure 1). This system consists of a series of interconnected small reservoirs that collectively store, regulate, and release incoming flow to paddy cultivations [3]. Apart from 1853 small tanks, this basin has several larger irrigation reservoirs, also built in ancient times but subsequently modified to suit the present-day demands of paddy cultivation. The combination of small and large reservoirs (which store about 43% of the MAR) strives to satisfy the irrigation and domestic needs of the basin during the long dry period from May to September. However, intense water demand and frequent droughts have led to water scarcity in the basin [43].
The Kalu Ganga basin, covering an area of 2296 km2, is about 30% smaller than Malwatu Oya, with a MAP of about 4000 mm. While the natural MAR of Kalu Ganga (8.45 km3) is more than 10 times that of Malwatu Oya, its coefficient of variation of annual flow (0.14) is much less. The basin’s elevation varies from 0 to 2250 m amsl and has a steep gradient in the upper catchment area, but milder low-lying slopes further downstream. Basin agriculture consists mainly of perennial rain-fed crops (as opposed to irrigated paddy in Malwatu Oya) such as tea and rubber. Flooding occurs in the low-lying areas of the basin almost every year [44]. Currently, there is only one human-made reservoir within the basin (Kukuleganga, a run-of-the-river hydropower reservoir storing only 0.02% of the MAR), but there are several proposals for more reservoirs to harness the abundant water resources here for inter-basin transfer and for developing hydropower [45].

2.2. Data and Reservoir System Configuration Scenarios

The main data used in the analyses are simulated time series of natural monthly flows for the period 1961–2013 at the outlets of 148 and 114 minor sub-basins of the Malwatu Oya and Kalu Ganga basins, respectively (Figure 1). Natural or unregulated flows—flows that would occur at a particular location of a river in the absence of any upstream water resources development—were simulated by using the Soil and Water Assessment Tool (SWAT) model. Figure 1 shows the locations of climate stations (6 for Malwatu Oya and 15 for Kalu Ganga, which recorded data for periods ranging from 1948 to 2013) that were used as inputs to the two models. The SWAT model for Malwatu Oya was calibrated and validated for the periods 1971–1975 and 1976–1979, respectively, against observed flows at a heavily regulated location (the only record available; Figure 1). The calibration efficiency was 0.89 for R2 and 0.86 for the Nash–Sutcliffe Coefficient (NSC), while the validation efficiency was 0.82 for R2 and 0.90 for the NSC. The Kalu Ganga model was calibrated and validated for the same periods against three gauging stations (2 regulated and 1 unregulated; Figure 1), obtaining calibration efficiencies of 0.90 for R2 and 0.87 for the NSC, and validation efficiencies of 0.83 for R2 and 0.88 for the NSC.
Subsequently, all existing reservoirs and water management practices were removed from the model setup and the models were executed under current land cover and land use conditions. The output flow time series from the modified models were assumed to represent pseudo-natural or at least unregulated flow. It is acknowledged that the validity of these simulated flows depends on having a high number of sites with observed flow data (both regulated and unregulated) against which to compare model results. However, since these pseudo-natural flows are used to illustrate a generic approach for long-term water storage planning rather than obtain results specific to the study basins, the limitations in the simulated flow time series are of little significance. The SWAT models only play a support function for developing generic methodologies (Section 2.3, Section 2.4, Section 2.5 and Section 2.6), which may be implemented in any other basin using relevant input data.
To develop generic methodologies, this research analyzed three hypothetical reservoir system configuration scenarios (Figure 2 and Figure 3). Configuration 1 consists of a single reservoir placed at the downstream end of the basin through which all runoff is routed. In Configuration 2, hypothetical reservoirs are placed near the outlets of 15 and 17 large sub-basins (those recognized by the Department of Agrarian Development) of the Malwatu Oya and Kalu Ganga basins, respectively. In Configuration 3, the sub-basins are further disaggregated into 148 and 114 minor sub-basins of the Malwatu Oya and Kalu Ganga basins, and hypothetical reservoirs are placed near their outlets. The main features of each configuration are presented in Table 1. The sizes of individual reservoirs decrease gradually from Configuration 1 to Configuration 3.

2.3. Overarching Surface Water Storage Development Pathways

In this paper, the term “water supply” means the supply of water to all consumptive uses including irrigation, industrial, domestic, and other purposes, which together make up a major portion of the water use portfolio of a river basin. We limit our focus to the water supply at this stage to simplify the analysis. The approach presented in this paper may be extended or adjusted to cover other non-consumptive water resource uses.
The economic benefit from a water supply reservoir may be measured in terms of the annual reservoir yield. This is an estimate of how much water can be supplied annually for all those purposes from a reservoir of a given capacity at a given reliability of supply [46] and is usually determined by analyzing the time series of inflows at the location of the reservoir. The annual yield that can be supplied at 100% reliability, for a given flow record by a reservoir of given capacity, which starts full and refills at least once after the worst drought on record, is termed the Safe Yield [47]. This paper refers to Safe Yield as the Water Supply (WS) Yield.
Storage-WS Yield curves (graph of WS Yield with increasing storage capacity) are usually used to design single reservoirs. However, in this research, Storage-WS Yield curves were used as a tool to illustrate the design of overarching surface water storage development pathways for an entire river basin under seven environmental flow (EF) release scenarios, by assuming that a single large reservoir is placed at the outlet of the entire river (reservoir system Configuration 1; Figure 2 and Figure 3).

2.3.1. Environmental Flow Release Scenarios

The seven EF release scenarios consisted of (a) a zero-release scenario where no specific EF releases were allowed except spills from the reservoir, and (b) six EF release scenarios generated by using the method of Smakhtin and Anputhas [48]. This method uses natural or unregulated monthly flow time series at a particular location as input to generate six other monthly flow time series which differ in magnitude but are similar in pattern to the natural flows. The method first estimates a flow duration curve of the natural flows, which is then shifted to the left along the horizontal axis to estimate six environmental flow duration curves that are subsequently converted to time series of environmental flows. These flow time series correspond to six Environmental Management Classes (EMCs), named F through A. The annual and monthly quantities of flow releases increase in magnitude when moving from Class F to Class A. Class F represents the ecological status of a highly degraded river, whereas Class A represents that of a protected healthy river (see [48,49] for further details). When reservoir system Configuration 1 is assumed, these time series of flows represent monthly EF requirements at the most downstream end of the river to maintain it in each of the EMCs. Figure 4 shows extracts of the monthly time series of EF release requirements (corresponding to the six EMCs) for the Malwatu Oya and Kalu Ganga basins that were generated in this manner starting with the natural flows.
In a similar context to WS Yield, the annual quantity of EF released downstream from a given reservoir is referred to as the EF Yield. Table 2 summarizes the minimum EF Yield (as a percentage of the MAR) expected to be released downstream under each of the seven EF release scenarios considered in this exercise.

2.3.2. Reservoir Simulation

The WS Yield corresponding to gradually increasing storage capacities under reservoir system Configuration 1 was estimated with the Sequent Peak Algorithm [50,51] as implemented in the Water Evaluation and Planning (WEAP) model. A varying monthly water demand, closely aligned with the actual temporal distribution of the annual water demand of the Malwatu Oya basin (nearly 65% of which is between October to February, corresponding to the main paddy cultivation season), was imposed on both basins to ensure that the results for the two basins were comparable. For each reservoir capacity, EF releases were varied from zero (hereafter referred to as Class No EF) to Class A so that seven corresponding WS Yield values were obtained. The largest WS Yield was obtained for Class No EF, whereas the smallest WS Yield was obtained for Class A. A range of seven overarching Storage-WS Yield curves were estimated for each basin in this manner (Figure 5). The term “overarching” is used to imply that these Storage-WS Yield curves are for the entire basin.
The basic equations of the Sequent Peak Algorithm are:
For i= 1 to N and k=1 to 12,
S i = S i 1 β k μ E F i 1 + Q i 1   If   S i > 0
Else,
S i = 0
S 0 = S c
where S i = reservoir storage at the beginning of the month i; Q i 1 = reservoir inflow during month i−1; S i 1 = reservoir storage at the beginning of the month i−1; μ = mean annual river flow; β k = fraction of annual demand expected in month k (corresponding to month i−1) out of 12 months of the year;   = total annual demand as a fraction of the mean annual river flow; E F i 1 = environmental flow release during month i−1; N = total number of months in the simulation; S 0 = reservoir storage at the beginning of the simulation, and S c = storage capacity of the reservoir.
The goal of the simulation is to find the maximum for each storage capacity S c that satisfies the condition S 0 = S c and S i = S c for at least one S i   after the worst drought on record.
Gross WS Yields were estimated assuming zero evaporation or seepage losses from reservoirs. Furthermore, reservoir operation did not consider any conservation, buffer, or inactive storage zones. These assumptions were for simplicity, although the omitted parameters may be included if deemed critical. Reservoir storage was allowed to become zero if it was necessary to satisfy the demands but was required to fill up at least once during the total simulation period (53 years) to ensure continuity of storage. While releasing EF, spills from the reservoir were allowed whenever the reservoir reached its full capacity.
Assuming that all streamflow except the WS Yield flows downstream from the reservoir, a set of seven storage-EF Yield curves were generated for each river (Figure 5).

2.4. Average Overarching Surface Water Storage Development Pathways

Average overarching Storage-WS Yield curves of the range of seven curves developed earlier for each basin were estimated by considering approaches (a) and (b) below:
(a)
Considering annual WS Yield
Under this approach, the mean of the two extreme WS Yield values (corresponding to classes No EF and A) for each storage capacity of the reservoir system Configuration 1 was estimated. This resulted in an average overarching Storage-WS Yield curve (hereafter referred to as the Mid curve) for each of the two basins under an average EF Yield scenario.
(b)
Considering monthly variation of WS Yield
This approach investigated ways in which a compromise may be reached between the two extreme Storage-WS Yield curves considering the actual monthly variation of the WS and EF Yields. A WS Sustainability Index and an EF Sustainability Index were formulated following [52] as described below.
  • WS Sustainability and EF Sustainability
R e l i a b i l i t y = N s N
R e s i l i e n c e = N u s N N s
V u l n e r a b i l i t y = i = 1 N m i N N s
R e l a t i v e   V u l n e r a b i l i t y = V u l n e r a b i l i t y T a r g e t   D e m a n d
  ( W S / E F )   S u s t a i n a b i l i t y = Reliability + Resilience + ( 1 Relative   Vulnerability ) 3
where Ns = number of months in which the target demand is met; N = total number of months in the period of analysis; Nus = number of months in which the supply has returned to the target demand after failure; and mi = the magnitude of failure in the ith failure month. The value of both WS Sustainability and EF Sustainability indices varies from 0 to 1.
2.
Compromising WS Sustainability and EF Sustainability
As estimated in Section 2.3, for each size of reservoir there is a range of WS Yields and EF Yields corresponding to different levels of ecological protection for the river (Figure 5). The behavior of WS Sustainability against EF Sustainability for each size of reservoir was investigated assuming: (a) water supply health is at its best (i.e., WS Sustainability = 1) when the WS Yield is at the higher end of the spectrum (EF Yield at Class No EF), and deteriorates with gradually increasing release levels of EF and, (b) the health of the river is at its best (i.e., EF Sustainability = 1) if the WS Yield is at the lower end of the spectrum (EF Yield at Class A), and deteriorates when gradually increasing levels of yields are drawn out. In practice, it is impossible to maintain both WS Sustainability and EF Sustainability at 1 and some compromise must be reached. This compromise may be achieved in several ways, for example, (1) maintaining EF Sustainability at 1 and allowing WS Sustainability to decline, (2) maintaining WS Sustainability at 1 and allowing EF Sustainability to decline, or (iii) allow both EF Sustainability and WS Sustainability to decline so that both requirements are partially met. In the latter instance, if both water supply and environmental flow demands are satisfied in equal proportions, EF Sustainability and WS Sustainability indices will become equal and minimum. The WS Yield at which the two indices become equal and minimum was estimated for gradually increasing storage capacities under reservoir system Configuration 1 in both basins. A second average overarching Storage-WS Yield curve (hereafter the “compromise curve”) was constructed for both basins.

2.5. Optimum Distribution of Storage Capacity and Withdrawal Thresholds in Sub-Basins

The overarching Storage-WS Yield curves developed in Section 2.3 and Section 2.4 provide an estimate of expected cumulative yields from a basin at different cumulative storage capacities. To establish the optimum distribution of this cumulative storage capacity and corresponding optimum yields within individual sub-basins, the impact of distributed reservoir system configurations 2 and 3 (Figure 2 and Figure 3) on the indicators WS Yield and EF Yield were investigated. Under Configurations 2 and 3, Storage-WS Yield curves for each sub-basin (large sub-basins under Configuration 2 and minor sub-basins under Configuration (3)) were generated using simulated natural runoff, without making explicit environmental flow releases other than reservoir spills. The individual curves were aggregated under two sub-configurations to generate a cumulative Storage-WS Yield curve for the entire basin. The two sub-configurations are: (a) each cumulative storage capacity is distributed equally among sub-basins, and (b) each cumulative storage capacity is distributed among sub-basins according to their MAR ratio. The intention was to identify the sub-configuration that produced the maximum aggregated WS Yield for the entire basin.

2.6. Regionalized Storage-Yield Curves

All Storage-WS Yield curves were generated by estimating yields at 11 selected storage capacities (0.0, 0.05, 0.1, 0.4, 0.8, 1, 1.4, 2, 3, 4, and 5 MAR units). These were selected to illustrate the shape of the entire curve, including the initial segment and where it levels off. A factor analysis by principal components of all Storage-WS Yield curves of the 148 and 114 minor sub-basins of the Malwatu Oya and Kalu Ganga basins was undertaken to regionalize the curves and identify the factors underlying their shapes.

3. Results and Discussion

3.1. Overarching Surface Water Storage Development Pathways

The Storage-WS Yield curves for both basins under the seven EF release scenarios are shown in Figure 5 (top) and the Storage-EF Yield curves under each EMC in Figure 5 (bottom). The curves have been normalized by dividing both axes by the MAR of each basin. These curves represent different overarching surface storage development pathways that may be adopted in the case study basins. The Storage-WS Yield curves illustrate the theoretical upper thresholds for cumulative annual surface water withdrawals in MAR units corresponding to a given cumulative storage capacity in the two basins under a range of ecological protection categories, from No EF releases to a stringent Environmental Management Class A. The corresponding Storage-EF Yield curves (Figure 5 bottom) illustrate the plausible annual EF releases under each surface storage development pathway.
Similar curves may be developed for other river basins. Any of the curves may be adopted as a potential storage development pathway for an entire river basin. However, the selected pathway will become a sustainable surface water storage development pathway only if it corresponds to the required level of ecosystem protection for a given river. For rivers for which such specific levels of ecosystem protection needs have not been identified, an average storage development pathway may be adopted as explained in Section 3.2.

3.2. Average Overarching Surface Water Storage Development Pathways and Acceptable Limits

Average overarching surface water storage development pathways were developed following the two approaches in Section 2.4. While a detailed discussion on the results of approach (a) is not required, the results of approach (b) are discussed in detail in this section.
Figure 6 illustrates the WS Sustainability versus EF Sustainability space for a reservoir of storage capacity with 1.0 MAR units under Configuration 1. The graphs at the top show the behavior of the WS Sustainability and EF Sustainability indices when an attempt is made to release gradually increasing EF Yields starting from its lowest level at Class No EF (where WS Sustainability=1). The graphs at the bottom show the behavior of the two indices when an attempt is made to draw out gradually increasing WS Yields starting from its lowest level at Class A (where EF Sustainability = 1). The arrows illustrate different pathways through which a compromise between the two sustainability indices may be made. Each path corresponds to a unique reservoir management strategy (in this case an entire river basin) and defines the priority assigned to meeting WS demands, meeting EF demands, and storing water (Table 3). On the path Priority: EF, the value of the EF Sustainability index is maintained at 1 while that of the WS Sustainability index is allowed to decline. On the path Priority: WS and Water Storage, the value of the WS Sustainability Index is maintained at 1 while that of the EF Sustainability index is allowed to decline. The path Priority: Equal EF and WS maintains two equal indices (but less than 1) by sharing water deficits in equal proportions between WS and EF requirements. The Compromise pathway was constructed by estimating the WS Yield at which the two indices reach a minimum on the path Priority: Equal EF and WS for gradually increasing storage capacities (Figure 7).
Along the fourth pathway, Priority: WS, values of both the EF Sustainability and WS Sustainability indices are less than 1, but the value of the WS Sustainability index is always higher than the EF Sustainability index since water supply demands are prioritized over EF demands and filling the reservoir. Although not exactly an average pathway, it is possible to construct another Compromise pathway (Compromise WS) by estimating the WS Yield at which the two indices reach a minimum on the path Priority: WS. In this case, meeting WS demands are prioritized over meeting EF demands. The exact coordinates of where a river basin is located in Figure 6 depend on the water supply demands, the level of ecological protection required by the river, and the river basin management strategy adopted in terms of priorities assigned to WS demands, EF demands, and reservoir filling. A Compromise pathway appropriate for the river in question can be formulated based on these coordinates.
The alignment of the Mid, Compromise, and Compromise WS pathways for the two basins compared to overarching pathways described in Section 3.1 is shown in Figure 7.
In rivers where no specific ecological protection requirement has been identified, an average pathway such as the Mid or Compromise curve (Figure 7) may be adopted. For example, for the Malwatu Oya basin, the Mid and Compromise curves envelop the Environmental Management Class C curve up to a storage capacity of 0.7 MAR units and this envelope moves to the Class B curve beyond a storage capacity of 1.3 MAR units (Figure 7 top). It can be inferred that up to a cumulative storage capacity of approximately 1.0 MAR units, the Class C curve represents a reasonable middle-ground storage development pathway unless there is a specific need for the river to be in a more protected state. Figure 7 (top) illustrates that marginal gains in WS Yield from storage capacities above 1.0 MAR units along Mid and Compromise pathways are minimal and that beyond this point accomplishing a satisfactory trade-off between water supply and environmental needs will be challenging. The mean of Mid and Compromise yields at this storage capacity is 45% of the MAR. Both curves level off at an average yield of 65%, corresponding to a storage capacity of 5.0 MAR units. Accordingly, it is possible to define a conservative limit to storage capacity expansion in the Malwatu Oya basin at 1.0 MAR units and 45% of MAR water withdrawals. A maximum withdrawal limit from surface storage may be set at 65% of the MAR, assuming that more stringent levels of ecosystem protection are not necessary.
Similarly, in the case of Kalu Ganga, considering the alignment of Mid and Compromise curves, either the Class C or the Class D curve may be adopted depending on ecosystem protection requirements. A conservative limit to storage capacity expansion in Kalu Ganga may be set at 0.5 MAR units and 50% of MAR water withdrawals and a maximum limit may be set at 60% of MAR, assuming that more stringent levels of ecosystem protection are not necessary.

3.3. Optimal Distribution of Storage Capacity and Withdrawal Thresholds in Sub-Basins

Cumulative WS Yield estimates for the entire basin under reservoir system Configurations 2 and 3 (Section 2.5) revealed that the cumulative yield from the entire basin is maximized when the cumulative storage capacity of the basin is distributed in sub-basins according to their MAR ratios (also observed by Pitman in 1995 [53]). However, the maximum WS Yield achieved from Configuration 1 was higher than the maximum aggregated WS Yield obtained from Configurations 2 and 3 for the same cumulative storage capacity (Figure 7). This illustrates that the maximum achievable WS Yield from a lumped reservoir of a certain capacity is higher than the cumulative yield obtained from a network of distributed reservoirs of the same aggregated capacity. The cumulative yield for the entire basin declines with subsequent levels of reservoir distribution but may stabilize at high levels of distribution since homogeneous regions are likely to have similar (normalized) storage versus WS Yield curves. The results imply that more releases are available in downstream parts of the basin as EF when small reservoirs are distributed in the river network, even if no explicit EF releases are made. Distributing storage capacity within sub-basins according to their MAR ratios and limiting abstractions to an appropriate fraction of the local runoff generated within each sub-basin is one strategy to ensure that cumulative water supply yields are maximized, the river basin as a whole complies with a pre-selected storage development pathway, and corresponding EF releases are maintained. Figure 8 and Figure 9 illustrate how such a strategy may be practically formulated in river basin planning.
Figure 8 shows the feasible storage capacity (S) and withdrawals (Y) in each large sub-basin, for the entire Malwatu Oya basin to stay in the Mid scenario if the cumulative storage capacity in the basin were to be increased to 1 MAR unit from its existing capacity of 0.43 MAR units. Figure 9 shows the same distribution for Kalu Ganga if the storage capacity in the basin were increased to 0.5 MAR units. The figures have been produced by distributing the cumulative storage capacity (1.0 MAR units in Malwatu Oya and 0.5 MAR units in Kalu Ganga) in large sub-basins according to their MAR ratios and estimating individual WS Yields corresponding to the storage capacity assigned to each large sub-basin using Storage-WS Yield curves developed under Configuration 2. To estimate the WS-Yield for each large sub-basin under the Mid scenario, individual WS Yields obtained under Configuration 2 were multiplied by the ratio between the ordinates of the Configuration 2 curve and the Mid curve for the entire basin (Figure 7) at 1.0 MAR units in Malwatu Oya and 0.5 MAR units in Kalu Ganga. Similar spatially disaggregated estimates may be made for other cumulative storage capacities and overarching pathways. It is also possible to adopt variable pathways and corresponding storage capacities and withdrawal thresholds for different sub-basins depending on ecological protection requirements.

3.4. Comparing Results with the Current Situation in the Case Study Basins and Practical Implications for Planning

The availability of more in-stream flow releases when reservoirs are distributed over the river network than when they are lumped together (Section 3.3) provides a scientific basis for the high level of reservoir distribution seen in the ancient tank irrigation system in the Malwatu Oya basin, although EF is probably not the only reason for this distribution. One aspect not covered in this research is the impact of evaporation on actual WS Yields when storage capacity is distributed over a large number of small reservoirs. Beyond a certain level of disaggregation of reservoir capacity over the river network, losses from evaporation may outweigh the gains in EF. Future research could investigate the relative impact of evaporation versus disaggregation of storage capacity on WS and EF Yields.
If the existing storage capacities in each minor sub-basin of the Malwatu Oya were lumped together (as in Configuration 3) and were drawing their maximum yields irrespective of EF releases, the cumulative annual 100% temporally reliable WS Yield would be 0.31 of the MAR of Malwatu Oya. Based on a comparison of current basin outflow with estimated natural outflow and actual demands, gross annual withdrawals in the basin were found to be as high as 0.55 MAR units with a time reliability of only 75–80%. The current annual in-stream flow release of 0.45 MAR units implies that the river is in Environmental Management Class B (Table 2). However, a close examination of the actual monthly time series of in-stream flow revealed that although the in-stream flow release is above that for Class B during the wet season, it drops to below that for Class F during the dry season. Hence, the river is in Class F under current practices.
Figure 10 classifies large sub-basins of Malwatu Oya by their current state of exploitation based on a comparison of their actual storage capacity and level of withdrawals with where they should be to place the basin as a whole in the Mid scenario (assuming the Mid scenario is adopted). It illustrates that four of the 15 major sub-basins are over-exploited in terms of storage, while five others are under-exploited. Table 4 compares the current spatial distribution of storage capacity among large sub-basins with the expected distribution to maximize yield gains, i.e., a distribution which follows the ratio among the MAR of large sub-basins. Sub-basins over-/under-exploited by more than 20% are highlighted in Table 4. Positive percentages in the last column imply over-exploitation, while negative percentages imply under-exploitation.
Any new surface water storage development in the Malwatu Oya basin needs to focus on under-exploited areas to maximize yield gains and ensure environmental protection. Since the basin is already saturated with surface storage, the development of sustainable groundwater may be the next best option to release the stress on surface water. In this context, regularization of informal groundwater development currently underway in the basin is a high priority. Any transfers from other basins need to ensure they do not unduly compromise the ecological integrity of the source basins. Since Kalu Ganga is relatively unregulated at present, these results may be useful for development planning.

3.5. Use of Regionalized Storage-Yield Curves as a Measure of the Surface Storage Potential of River Basins and as a Basin Characteristic

The shape of the storage-yield curve for a particular river basin depends on the characteristics of the inflow time series, which in turn depend on the climate, geomorphology, and geological features of the basin. For example, the curves for the Malwatu Oya and its sub-basins have a gradual gradient in the initial segment before leveling off, whereas those for Kalu Ganga have a much steeper gradient in the initial segment. Yield gains from storage are reached much quicker in the latter than in the former. The maximum potential yield (No EF) curve for Malwatu Oya levels off at a storage capacity of 3.8 MAR units (3001 million m3), whereas in Kalu Ganga this point is reached much earlier at a capacity of 0.4 MAR units (3381 million m3) (Figure 7). The corresponding yields are 772 and 7312 million m3, respectively. Although the storage capacities at which the maximum yields are reached are of the same order of magnitude, the corresponding yields differ by an order of magnitude 10.
A factor analysis by principal components of the set of curves for the 148 minor sub-basins of the Malwatu Oya revealed that 88% of the variance in the shape of the curves (especially for storage capacities above 0.1 MAR units) may be explained by one common factor, which is highly positively correlated with the mean and standard deviation, and negatively correlated with the skewness of annual precipitation of individual sub-basins. Therefore, 88% of the variance in runoff generated by the minor sub-basins in the Malwatu Oya and their storage-yield curves can be explained by their annual precipitation statistics. Considering the variation of factor values across the basin, Malwatu Oya may be represented by six regional storage-yield curves (Figure 11; Table 5).
In the Kalu Ganga, 92% of the variance in the curves could be explained by two common factors, with the first factor accounting for 77%. The first factor explains the variance in the flatter part of the curves beyond 0.4 MAR units, and the second explains the variance in the initial gradient up to 0.4 MAR units. Variance in the flatter part of the curves beyond 0.4 MAR units (explained by factor 1) is negatively correlated with the standard deviation, and positively correlated with the skewness of annual precipitation. It was also found that variance in the initial gradient of the curves up to 0.4 MAR units (explained by factor 2) is positively correlated with the mean annual precipitation. Therefore, 92% of the variance in runoff generated by individual minor sub-basins in the Kalu Ganga and their storage-yield curves can be explained by their annual precipitation statistics. Considering the variation of the common factors across the basin, Kalu Ganga is represented by nine regional storage-yield curves (Figure 12, Table 5).
The differences in storage-yield curves among basins and sub-basins within each basin can be explained by differences in the statistics of annual precipitation. However, since simulated runoff was used, the characteristics of the storage-yield curves depend on the climate and catchment characteristics input during modeling. Hence, the actual extent to which storage yield curves are influenced by characteristics of annual precipitation versus other catchment characteristics needs to be verified at places where substantial measured flow data are available.
According to the generalized global storage–yield–reliability relationship developed by Kuria and Vogel [54,55], standardized yields (yields as a fraction of the MAR) in the range 0.3–0.8 from a given reservoir capacity can be expressed by a combination of the statistics of the inflow time series (a combination of mean, standard deviation, and skewness coefficient of annual inflows; Equation (9)). The possibility of replacing the statistics of annual inflows with the statistics of annual precipitation needs to be further explored and would be beneficial for storage planning in regions where measured and modeled runoff data are not readily available. Emerging methods for generating synthetic flow records such as continuous simulation modeling (e.g., [56,57]) may also prove useful to establish storage-yield curves in ungauged basins.
Y = 0.651 S 0.203 Z r 0.306 μ 1.135 σ 0.342 γ 0.017
where Y = yield; S = storage capacity in million cubic metres; Z r = standard normal variate with r equal to reliability; μ , σ , γ =   the mean, standard deviation, and skewness coefficient of annual inflows.

4. Conclusions

The results of this research illustrate that Storage-WS yield curves for entire river basins under different EF release scenarios are a useful approach for designing overarching sustainable storage development pathways and acceptable limits to surface storage development in river basins. Storage-WS Yield curves and the indices of WS Sustainability and EF Sustainability may be conjunctively used to identify the most appropriate storage development pathway and acceptable storage limits for a given river basin so that aquatic ecosystem health and related services are not unduly compromised as the basin’s cumulative storage capacity increases. The research demonstrates that the optimal arrangement of storage capacity is to have the cumulative capacity distributed among sub-basins according to the MAR ratio among them, and that withdrawals from each sub-basin may be appropriately limited to maintain the entire basin on the identified pathway.
While this approach may help plan long-term storage development in relatively under-developed river basins, it may also be used to compare the current status of developed basins against what is acceptable, and on that basis to propose remedies to any deviations. Alternative approaches for establishing storage-yield curves including the use of synthetic flow records may prove a fruitful future direction in research to implement this paper’s recommendations in regions with limited river flow data.
The research fills an important knowledge gap on the questions of acceptable storage limits and optimal reservoir arrangements for river basins. This generic approach can be adapted to individual river basins by: (1) using reliabilities of supply less than 100%; (2) using alternative EF estimation methods; (3) formulating an alternative EF Sustainability index, measuring the extent to which instream flow releases differ from natural flows, instead of its current formulation; (4) considering evaporation and percolation losses and inactive storage zones for reservoirs; (5) using daily flow data instead of monthly flow data especially in the case of smaller capacity reservoirs; (6) including indicators on economic sustainability to compare different reservoir system configurations; (7) expanding the approach to cover other water use sectors including non-consumptive ones and, (8) considering future climate and socioeconomic scenarios.

Author Contributions

Conceptualization, N.E., V.S. and L.U.; methodology, N.E.; formal analysis, N.E.; investigation, N.E.; writing—original draft preparation, N.E.; writing—review & editing, V.S. and L.U.; supervision, V.S. and L.U. All authors have approved the content of the submitted manuscript.

Funding

This research was supported by the International Water Management Institute as part of the CGIAR Research Programs on Climate Change Agriculture and Food Security and Water Land and Ecosystems.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Simulated runoff data used in the analysis may be made available on request from the authors.

Acknowledgments

The authors are grateful to Matthew McCartney, Research Group Leader, Sustainable Water Infrastructure and Ecosystems (IWMI) for his constructive feedback on the manuscript and would also like to thank Lal Mutuwatte (IWMI) and Madusanka Thilakarathne (Consultant, IWMI) for their help in developing the SWAT models.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McCartney, M.; Smakhtin, V. Water storage in an era of climate change: Addressing the challenge of increasing rainfall variability. In IWMI Blue Paper; International Water Management Institute: Colombo, Sri Lanka, 2010. [Google Scholar]
  2. McCully, P. Silenced Rivers: The Ecology and Politics of Large Dams, 2nd ed.; Zed Books: London, UK; New York, NY, USA, 1996. [Google Scholar] [CrossRef]
  3. Panabokke, C.R.; Sakthivadivel, R.; Weerasinghe, A.D. Small Tanks in Sri Lanka: Evolution, Present Status and Issues; International Water Management Institute (IWMI): Colombo, Sri Lanka, 2002. [Google Scholar]
  4. Kekulandala, B.; Jacobs, B.; Cunningham, R. Management of small irrigation tank cascade systems (STCS) in Sri Lanka: Past, present and future. Clim. Dev. 2020. [Google Scholar] [CrossRef]
  5. Wada, Y.; Reager, J.T.; Chao, B.F.; Wang, J.; Lo, M.H.; Song, C.; Li, Y.; Gardner, A.S. Recent changes in land water storage and its contribution to sea level variations. Surv. Geophys. 2017. [Google Scholar] [CrossRef] [Green Version]
  6. Hogeboom, R.J.; Knook, L.; Hoekstra, A.Y. The blue water footprint of the world’s artificial reservoirs for hydroelectricity, irrigation, residential and industrial water supply, flood protection, fishing and recreation. Adv. Water Resour. 2018, 113, 285–294. [Google Scholar] [CrossRef]
  7. Shi, H.; Chen, J.; Liu, S.; Sivakumar, B. The role of large dams in promoting economic development under the pressure of population growth. Sustainability 2019, 11, 2965. [Google Scholar] [CrossRef] [Green Version]
  8. Hanasaki, N.; Yoshikawa, S.; Pokhrel, Y.; Kanae, S. A global hydrological simulation to specify the sources of water used by humans. Hydrol. Earth Syst. Sci. 2018, 22, 789–817. [Google Scholar] [CrossRef] [Green Version]
  9. International Energy Agency (IEA). Key World Energy Statistics; International Energy Agency: Paris, France, 2018. [Google Scholar]
  10. International Hydropower Association (IHA). Hydropower Status Report; International Hydropower Association: London, UK, 2018. [Google Scholar]
  11. Conker, A.; Hussein, H. Hydraulic mission at home, hydraulic mission abroad? Examining Turkey’s regional “Pax-Aquarum” and its limits. Sustainability 2019, 11, 228. [Google Scholar] [CrossRef] [Green Version]
  12. Hussein, H.; Conker, A.; Grandi, M. Small is beautiful but not trendy: Understanding the allure of big hydraulic works in the Euphrates-Tigris and Nile waterscapes. Mediterr. Polit. 2020. [Google Scholar] [CrossRef]
  13. Menga, F. Building a nation through a dam: The case of Rogun in Tajikistan. Natl. Pap. 2015. [Google Scholar] [CrossRef]
  14. Conker, A.; Hussein, H. Hydropolitics and issue-linkage along the Orontes River Basin: An analysis of the Lebanon–Syria and Syria–Turkey hydropolitical relations. Int. Environ. Agreem. Polit. Law Econ. 2020. [Google Scholar] [CrossRef] [Green Version]
  15. Vij, S.; Warner, J.; Barua, A. Power in water diplomacy. Water Int. 2020. [Google Scholar] [CrossRef]
  16. Grill, G.; Lehner, B.; Thieme, M.; Geenen, B.; Tickner, D.; Antonelli, F.; Babu, S.; Borrelli, P.; Cheng, L.; Crochetiere, H.; et al. Mapping the world’s free-flowing rivers. Nature 2019, 569, 215–221. [Google Scholar] [CrossRef]
  17. Anderson, E.P.; Jenkins, C.N.; Heilpern, S.; Maldonado-Ocampo, J.A.; Carvajal-Vallejos, F.M.; Encalada, A.C.; Rivadeneira, J.F.; Hidalgo, M.; Cañas, C.M.; Ortega, H.; et al. Fragmentation of Andes-to-Amazon connectivity by hydropower dams. Sci. Adv. 2018, 4. [Google Scholar] [CrossRef] [Green Version]
  18. Puijenbroek, P.J.T.M.; Buijse, A.D.; Kraak, M.H.S.; Verdonschot, P.F.M. Species and river specific effects of river fragmentation on European anadromous fish species. River Res. Appl. 2019, 35, 68–77. [Google Scholar] [CrossRef]
  19. Lehner, B.; Reidy Liermann, C.; Revenga, C.; Vörösmarty, C.; Fekete, B.; Crouzet, P.; Döll, P.; Endejan, M.; Frenken, K.; Magome, J.; et al. High resolution mapping of the world’s reservoirs and dams for sustainable river-flow management. Front. Ecol. Environ. 2011, 9, 494–502. [Google Scholar] [CrossRef] [Green Version]
  20. Bejarano, M.D.; Sordo-Ward, Á.; Alonso, C.; Nilsson, C. Characterizing effects of hydropower plants on sub-daily flow regimes. J. Hydrol. 2017, 550, 186–200. [Google Scholar] [CrossRef] [Green Version]
  21. Dong, N.; Yu, Z.; Yang, C.; Yang, M.; Wang, W. Hydrological impact of a reservoir network in the upper Gan River Basin, China. Hydrol. Process. 2019, 33, 1709–1723. [Google Scholar] [CrossRef]
  22. Beck, M.W.; Claassen, A.H.; Hundt, P.J. Environmental and livelihood impacts of dams: Common lessons across development gradients that challenge sustainability. Int. J. River Basin Manag. 2012. [Google Scholar] [CrossRef]
  23. Poff, N.L.R.; Olden, J.D. Can dams be designed for sustainability? Science 2017, 358, 1252–1253. [Google Scholar] [CrossRef] [PubMed]
  24. Soukhaphon, A.; Baird, I.G.; Hogan, Z.S. The Impacts of Hydropower Dams in the Mekong River Basin: A Review. Water 2021, 13, 265. [Google Scholar] [CrossRef]
  25. Van der Zaag, P.; Gupta, J. Scale issues in the governance of water storage projects. Water Resour. Res. 2008, 44. [Google Scholar] [CrossRef] [Green Version]
  26. Richter, B.D.; Postel, S.; Revenga, C.; Scudder, T.; Lehner, B.; Churchill, A.; Chow, M. Lost in development’s shadow: The downstream human consequences of dams. Water Altern. 2010, 3, 14–42. [Google Scholar]
  27. Zarfl, C.L.; Lumsdon, A.E.; Berlekamp, J.; Tydecks, L.; Tockner, K. A global boom in hydropower dam construction. Aquat. Sci. 2014, 77, 161–170. [Google Scholar] [CrossRef]
  28. Scudder, T. Development-induced community resettlement. In New Directions in Social Impact Assessment: Conceptual and Methodological Advances; Vanclay, F., Esteves, A.M., Eds.; Edward Elgar: Cheltanham, UK; Northhampton, MA, USA, 2011. [Google Scholar] [CrossRef]
  29. Shah, E.; Boelens, R.; Bruins, B. Reflections: Contested epistemologies on large dams and mega-hydraulic development. Water 2019, 11, 417. [Google Scholar] [CrossRef] [Green Version]
  30. International Commission on Large Dams (ICOLD). Available online: https://www.icold-cigb.org/GB/world_register/general_synthesis.asp (accessed on 22 February 2021).
  31. United Nations Economic and Social Commission for Asia and the Pacific (UN ESCAP). Enhancing Regional Cooperation in Infrastructure Development Including that Related to Disaster Management; United Nations: New York, NY, USA, 2006. [Google Scholar] [CrossRef]
  32. Seckler, D. The New Era of Water Resources Management: From ‘Dry’ To ‘Wet’ Water Savings; IIMI Research Report 1. International Irrigation Management Institute: Colombo, Sri Lanka, 1996. [Google Scholar]
  33. Richter, B.D. Re-thinking environmental flows: From allocations and reserves to sustainability boundaries. River Res. Appl. 2009, 26, 1052–1063. [Google Scholar] [CrossRef]
  34. Hoekstra, A.Y. Sustainable, efficient, and equitable water use: The three pillars under wise freshwater allocation. Wiley Interdiscip. Rev. 2014, 1, 31–40. [Google Scholar] [CrossRef]
  35. Jager, H.I.; Efroymson, R.A.; Opperman, J.J.; Kelly, M.R. Spatial design principles for sustainable hydropower development in river basins. Renew. Sustain. Energy Rev. 2015. [Google Scholar] [CrossRef] [Green Version]
  36. Roozbahani, R.; Abbasi, B.; Schreider, S.; Hosseinifard, Z. A basin-wide approach for water allocation and dams location-allocation. Ann. Oper. Res. 2019. [Google Scholar] [CrossRef]
  37. González-Zeas, D.; Garrote, L.; Iglesias, A. Assessing maximum potential water withdrawal for food production under climate change: An application in Spain. J. Water Clim. Chang. 2014, 5, 633–651. [Google Scholar] [CrossRef]
  38. Porse, E.C.; Sandoval-Solis, S.; Lane, B.A. Integrating Environmental Flows into Multi-Objective Reservoir Management for a Transboundary, Water-Scarce River Basin: Rio Grande/Bravo. Water Resour. Manag. 2015. [Google Scholar] [CrossRef]
  39. Rheinheimer, D.E.; Liu, P.; Guo, S. Re-operating the Three Gorges reservoir for environmental flows: A preliminary assessment of trade-offs. River Res. Appl. 2016, 32, 257–266. [Google Scholar] [CrossRef]
  40. Chen, M.; Dong, Z.; Jia, W.; Ni, X.; Yao, H. Multi-objective joint optimal operation of reservoir system and analysis of objectives competition mechanism: A case study in the upper reach of the Yangtze River. Water 2019, 11, 2542. [Google Scholar] [CrossRef] [Green Version]
  41. Jiang, Z.; Liu, P.; Ji, C.; Zhang, H.; Chen, Y. Ecological flow considered multi-objective storage energy operation chart optimization of large-scale mixed reservoirs. J. Hydrol. 2019. [Google Scholar] [CrossRef]
  42. Eriyagama, N.; Smakhtin, V.; Udamulla, L. How much artificial surface storage is acceptable in a river basin and where should it be located: A review. Earth-Sci. Rev. 2020, 208, 103294. [Google Scholar] [CrossRef]
  43. Meegastenna, T.J. Assessment of climate change impacts and adaptation measures to Malwatu Oya river basin in North Central province of Sri Lanka. In International Workshop on Innovation of Developing the Strategy for Impact Assessment of and Adaptation to the Climate Change as the “New Normal” (CLIMATE). In Proceedings of the 3rd World Irrigation Forum, Bali, Indonesia, 1–7 September 2019; International Commission on Irrigation and Drainage (ICID). 2019. Available online: http://www.icid.org/wif3_bali_2019/wif3_ws_climate.pdf (accessed on 22 February 2021).
  44. Nandalal, H.K.; Ratnayake, U.R. Setting up of indices to measure vulnerability of structures during a flood. In Proceedings of the International Conference on Sustainable Built Environment (ICSBE-2010), Kandy, Sri Lanka, 13–14 December 2010; Available online: https://core.ac.uk/download/pdf/141733401.pdf (accessed on 22 February 2021).
  45. Dharmasena, G.T. Integrated Water Resources Management and Sound Information System–Sri Lankan Experience. In Sustainable Development of Water Resources, Water Supply and Environmental Sanitation. In Proceedings of the 32nd WEDC International Conference, Colombo, Sri Lanka, 13–17 November 2006; Fisher, J., Ed.; WEDC, Loughborough University: Loughborough, UK, 2006; pp. 239–242. Available online: https://repository.lboro.ac.uk/articles/conference%20contribution/Integrated_water_resources_management_and_sound_information_system_-_Sri_Lankan_experience/9596771/1 (accessed on 22 February 2021).
  46. Vogel, R.M.; Stedinger, J.R. Generalized storage-reliability-yield relationships. J. Hydrol. 1987, 89, 303–327. [Google Scholar] [CrossRef]
  47. Vogel, R.; Seiber, J.; Archfield, S.; Smith, M.; Apse, C.; Huber-Lee, A. Relations among storage, yield, and instream flow. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef]
  48. Smakhtin, V.U.; Anputhas, M.A. Assessment of Environmental Flow Requirements of Indian River Basins; IWMI Research Report 107; International Water Management Institute: Colombo, Sri Lanka, 2006; p. 36. [Google Scholar]
  49. Smakhtin, V.U.; Eriyagama, N. Developing a software package for global desktop assessment of environmental flows. Environ. Model. Softw. 2008, 23, 1396–1406. [Google Scholar] [CrossRef]
  50. Thomas, H.A., Jr.; Burden, R.P. Operations Research in Water Quality Management. Final Report to Bureau of State Services, Division of Water Supply and Pollution Control, Public Health Service, Department of Health, Education and Welfare [under] Contract PH86-62-140; Division of Engineering and Applied Physics, Harvard University: Cambridge, MA, USA, 1963. [Google Scholar]
  51. Rippl, W. The capacity of storage-reservoirs for water supply. In Minutes of the Proceedings of the Institute of Civil Engineers 71, 1883; Paper Number 1864, pp. 270–278. Available online: https://www.icevirtuallibrary.com/doi/abs/10.1680/imotp.1883.21797 (accessed on 22 February 2021).
  52. McCartney, M.; Rebelo, L.; Xenarios, S.; Smakhtin, V. Agricultural water storage in an era of climate change: Assessing need and effectiveness in Africa. In IWMI Research Report; 152. 37p; International Water Management Institute: Colombo, Sri Lanka, 2013. [Google Scholar] [CrossRef] [Green Version]
  53. Pitman, W.V. Towards Improved Utilization of South Africa’s Water Resources. In Proceedings of the Seventh South African National Hydrology Symposium, Grahamstown, South Africa, 4–6 September 1995; Institute for Water Research, Rhodes University: Grahamstown, South Africa, 1995. [Google Scholar]
  54. Kuria, F.; Vogel, R. Global storage-reliability-yield relationships for water supply reservoirs. Water Resour. Manag. 2014, 29, 1591–1605. [Google Scholar] [CrossRef]
  55. Kuria, F.W.; Vogel, R.M. A global water supply reservoir yield model with uncertainty analysis. Environ. Res. Lett. 2014. [Google Scholar] [CrossRef] [Green Version]
  56. Grimaldi, S.; Nardi, F.; Piscopia, R.; Petroselli, A.; Apollonio, C. Continuous hydrologic modelling for design simulation in small and ungauged basins: A step forward and some tests for its practical use. J. Hydrol. 2020. [Google Scholar] [CrossRef]
  57. Winter, B.; Schneeberger, K.; Dung, N.V.; Huttenlau, M.; Achleitner, S.; Stötter, J.; Merz, B.; Vorogushyn, S. A continuous modelling approach for design flood estimation on sub-daily time scale. Hydrol. Sci. J. 2019. [Google Scholar] [CrossRef]
Figure 1. Study basins in climatic zones of Sri Lanka (center), the sub-basin structure, and hydrography of the Malwatu Oya (top), and the Kalu Ganga (bottom) basins.
Figure 1. Study basins in climatic zones of Sri Lanka (center), the sub-basin structure, and hydrography of the Malwatu Oya (top), and the Kalu Ganga (bottom) basins.
Water 13 00645 g001
Figure 2. Three reservoir system configurations in Malwatu Oya ranging from one large hypothetical reservoir at the downstream end to 15 medium-size reservoirs in large sub-basins and 148 smaller reservoirs in minor sub-basins.
Figure 2. Three reservoir system configurations in Malwatu Oya ranging from one large hypothetical reservoir at the downstream end to 15 medium-size reservoirs in large sub-basins and 148 smaller reservoirs in minor sub-basins.
Water 13 00645 g002
Figure 3. Three reservoir system configurations in Kalu Ganga ranging from one large hypothetical reservoir at the downstream end to 17 medium-size reservoirs in large sub-basins and 114 smaller reservoirs in minor sub-basins.
Figure 3. Three reservoir system configurations in Kalu Ganga ranging from one large hypothetical reservoir at the downstream end to 17 medium-size reservoirs in large sub-basins and 114 smaller reservoirs in minor sub-basins.
Water 13 00645 g003
Figure 4. Extracts from Environmental Flow time series required at the outlets of Malwatu Oya (top) and Kalu Ganga (bottom) for the rivers to be maintained in six Environmental Management Classes (EMCs).
Figure 4. Extracts from Environmental Flow time series required at the outlets of Malwatu Oya (top) and Kalu Ganga (bottom) for the rivers to be maintained in six Environmental Management Classes (EMCs).
Water 13 00645 g004
Figure 5. Behavior of Water Supply Yield (top) and Environmental Flow Yield (bottom) under different Environmental Management Classes (EMC) for a hypothetical reservoir through which all river runoff is routed for Malwatu Oya (left) and Kalu Ganga (right).
Figure 5. Behavior of Water Supply Yield (top) and Environmental Flow Yield (bottom) under different Environmental Management Classes (EMC) for a hypothetical reservoir through which all river runoff is routed for Malwatu Oya (left) and Kalu Ganga (right).
Water 13 00645 g005
Figure 6. Variation of WS Sustainability and EF Sustainability for Malwatu Oya (left) and Kalu Ganga (right) basins when gradually increasing levels of EF are released (top) and gradually increasing levels of yields are withdrawn (bottom) from a reservoir with a storage capacity of 1.0 MAR units, with different priorities assigned to water supply, environmental flow provision, and water storage.
Figure 6. Variation of WS Sustainability and EF Sustainability for Malwatu Oya (left) and Kalu Ganga (right) basins when gradually increasing levels of EF are released (top) and gradually increasing levels of yields are withdrawn (bottom) from a reservoir with a storage capacity of 1.0 MAR units, with different priorities assigned to water supply, environmental flow provision, and water storage.
Water 13 00645 g006
Figure 7. Behavior of WS Yield with storage capacity for different Environmental Management Classes, curves of Mid, Compromise, and Compromise WS, and different scenarios of reservoir system configurations for Malwatu Oya (top) and Kalu Ganga (bottom).
Figure 7. Behavior of WS Yield with storage capacity for different Environmental Management Classes, curves of Mid, Compromise, and Compromise WS, and different scenarios of reservoir system configurations for Malwatu Oya (top) and Kalu Ganga (bottom).
Water 13 00645 g007
Figure 8. Feasible storage capacity (S) and withdrawals (Y) in each large sub-basin of Malwatu Oya for the entire basin to stay on the Mid scenario if the cumulative storage capacity were increased to 1.0 MAR unit. S and Y values given in terms of MAR units of the entire basin.
Figure 8. Feasible storage capacity (S) and withdrawals (Y) in each large sub-basin of Malwatu Oya for the entire basin to stay on the Mid scenario if the cumulative storage capacity were increased to 1.0 MAR unit. S and Y values given in terms of MAR units of the entire basin.
Water 13 00645 g008
Figure 9. Feasible storage capacity (S) and withdrawals (Y) in each large sub-basin of Kalu Ganga for the entire basin to stay on the Mid scenario if the cumulative storage capacity were increased to 0.5 MAR units. S and Y values given in terms of MAR units of the entire basin.
Figure 9. Feasible storage capacity (S) and withdrawals (Y) in each large sub-basin of Kalu Ganga for the entire basin to stay on the Mid scenario if the cumulative storage capacity were increased to 0.5 MAR units. S and Y values given in terms of MAR units of the entire basin.
Water 13 00645 g009
Figure 10. Classification of large sub-basins of Malwatu Oya by their current state of exploitation of surface storage. Sub-basin numbers shown in black. Sub-basin details provided in Table 4.
Figure 10. Classification of large sub-basins of Malwatu Oya by their current state of exploitation of surface storage. Sub-basin numbers shown in black. Sub-basin details provided in Table 4.
Water 13 00645 g010
Figure 11. Regional Storage-Yield Curves of Malwatu Oya with minor sub-basins having similar curves indicated by color and number. Graphs show storage capacity on the horizontal axis; WS Yield on the vertical axis in MAR units.
Figure 11. Regional Storage-Yield Curves of Malwatu Oya with minor sub-basins having similar curves indicated by color and number. Graphs show storage capacity on the horizontal axis; WS Yield on the vertical axis in MAR units.
Water 13 00645 g011
Figure 12. Regional Storage-Yield Curves of Kalu Ganga with minor sub-basins having similar curves indicated by color and number. Graphs show storage capacity on the horizontal axis; WS Yield on the vertical axis in MAR units.
Figure 12. Regional Storage-Yield Curves of Kalu Ganga with minor sub-basins having similar curves indicated by color and number. Graphs show storage capacity on the horizontal axis; WS Yield on the vertical axis in MAR units.
Water 13 00645 g012
Table 1. Main features of reservoir system Configurations 1, 2, and 3.
Table 1. Main features of reservoir system Configurations 1, 2, and 3.
BasinMalwatu OyaKalu Ganga
Configuration123123
Number of reservoirs115148117114
Average catchment area per reservoir (km2)3338 22322.6229617225.7
Average catchment area as a fraction of total area 10.070.00710.060.009
MAR of average catchment area (km3)0.790.060.0068.450.510.075
MAR of average catchment area as a fraction of total MAR10.070.00710.060.009
Note: MAR stands for Mean Annual Runoff.
Table 2. Environmental Flow (EF) Yield requirement at the downstream end of Malwatu Oya and Kalu Ganga basins under seven EF release scenarios.
Table 2. Environmental Flow (EF) Yield requirement at the downstream end of Malwatu Oya and Kalu Ganga basins under seven EF release scenarios.
EF Release ScenarioEnvironmental Management Class (EMC)EF Releases as Percentage of MAR
Malwatu OyaKalu Ganga
1No EF00
2F9.725
3E12.633.3
4D16.643.1
5C23.254.6
6B35.567.8
7A58.382.8
Note: The zero-EF release scenario is shown as EMC No EF.
Table 3. Reservoir and river basin management strategies through which water supply (WS) and environmental flow (EF) demands can be fully or partially met (see Figure 6).
Table 3. Reservoir and river basin management strategies through which water supply (WS) and environmental flow (EF) demands can be fully or partially met (see Figure 6).
Path (Reservoir/River Basin Management Strategy)Explanation
Priority: EF EF releases are made first before satisfying WS demands; any excess water is stored in the reservoir/river basin after satisfying both demands.
Priority: Equal EF and WSEqual and higher priority is assigned to meeting both EF and WS demands than storing in the reservoir/river basin. Shortages are shared in equal proportion between the two demands. Water is stored in the reservoir/river basin only after meeting both demands.
Priority: WSWS demands are met first, EF releases next, and excess water is stored in the reservoir/river basin after satisfying both demands.
Priority: WS and Water StorageWS demands are met first, storage is allowed next; EF needs are only satisfied through storage spills.
Table 4. Degree of over-/under-exploitation of large sub-basins of Malwatu Oya. Colors blue, red and green correspond to under, over and correct exploitation levels.
Table 4. Degree of over-/under-exploitation of large sub-basins of Malwatu Oya. Colors blue, red and green correspond to under, over and correct exploitation levels.
No.Sub-basinMAR
(mill. m3)
Existing Storage Capacity (mill. m3) Existing Storage Capacity as % of MARMAR as % of the Total Basin MAR Existing Storage Capacity as % of Total Storage CapacityDegree of over-/ under-Exploitation
1Kadahatu Oya26.3610.7540.793.223.11−3.29
2Narivili Aru27.2010.2737.743.322.97−10.53
3Maha Kanadara Oya40.1931.9279.444.909.23+88.35
4Maminiya Oya35.1116.8548.004.284.87+13.81
5Sangili Kanadara Oya35.6517.9550.344.355.19+19.36
6Upper Kanadara Oya49.0335.7172.835.9810.33+72.69
7Upper Malwatu Oya71.9135.9750.028.7710.40+18.60
8Boo Oya76.3542.8356.099.3112.39+33.00
9Middle Sangili Kanadara Oya14.185.4538.451.731.58−8.82
10Lower Weli Oya28.322.779.773.450.80−76.83
11Kal Aru41.673.287.865.080.95−81.35
12Upper Middle Malwatu Oya161.94114.2470.5419.7533.04+67.26
13Lower Kanadara Oya41.966.6515.845.121.92−62.43
14Lower Middle Malwatu Oya59.377.2612.237.242.10−70.99
15Lower Malwatu Oya110.563.863.4913.491.12−91.72
Total789.70345.7543100100
Table 5. Regional storage-yield curves represented as flow tables for Malwatu Oya and Kalu Ganga basins (Figure 11, Figure 12).
Table 5. Regional storage-yield curves represented as flow tables for Malwatu Oya and Kalu Ganga basins (Figure 11, Figure 12).
Storage Capacity (MAR units) WS Yield (MAR units)
Malwatu OyaKalu Ganga
Region Number
123456123456789
0.000.00.000.010.010.010.020.000.010.000.000.000.010.020.030.03
0.050.070.070.090.110.120.130.060.160.080.160.170.190.220.230.23
0.10.100.110.130.160.170.210.100.230.150.280.290.310.330.370.40
0.40.230.260.300.360.390.460.220.370.420.670.630.660.740.760.84
0.80.300.370.440.510.600.630.300.480.690.850.780.750.810.880.94
1.00.330.400.490.560.660.690.330.540.800.920.830.780.850.900.95
1.40.390.470.550.630.730.780.410.650.840.960.890.810.860.930.97
20.480.580.640.730.820.890.520.820.890.990.970.870.880.960.99
30.620.710.780.850.920.970.670.950.951.021.000.940.911.001.01
40.730.830.870.950.991.000.761.001.001.051.031.000.941.031.04
50.810.910.931.001.011.030.831.041.031.071.051.040.971.061.07
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eriyagama, N.; Smakhtin, V.; Udamulla, L. Sustainable Surface Water Storage Development Pathways and Acceptable Limits for River Basins. Water 2021, 13, 645. https://doi.org/10.3390/w13050645

AMA Style

Eriyagama N, Smakhtin V, Udamulla L. Sustainable Surface Water Storage Development Pathways and Acceptable Limits for River Basins. Water. 2021; 13(5):645. https://doi.org/10.3390/w13050645

Chicago/Turabian Style

Eriyagama, Nishadi, Vladimir Smakhtin, and Lakshika Udamulla. 2021. "Sustainable Surface Water Storage Development Pathways and Acceptable Limits for River Basins" Water 13, no. 5: 645. https://doi.org/10.3390/w13050645

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