Next Article in Journal
Research Advances on Financial Inclusion: A Bibliometric Analysis
Next Article in Special Issue
On Sustainable Consumption: The Implications of Trade in Virtual Water for the EU’s Food Security
Previous Article in Journal
Pollen and Fungal Spores Evaluation in Relation to Occupants and Microclimate in Indoor Workplaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Info-Gap Models for Optimal Multi-Year Management of Regional Water Resources Systems under Uncertainty

Department of Natural Resources and Environmental Management, University of Haifa, Haifa 3498838, Israel
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(6), 3152; https://doi.org/10.3390/su13063152
Submission received: 11 February 2021 / Revised: 2 March 2021 / Accepted: 8 March 2021 / Published: 13 March 2021
(This article belongs to the Special Issue Sustainable Water Resource Management in a Changing Climate)

Abstract

:
The common practices for the planning and management of Water Resources Systems (WSSs) have been challenged in the last few decades by global climate change processes, which are observed around the world in increasing frequencies. Climate change is manifested by climate variability, temperature increase, and extreme events such as droughts and floods, which have a decisive effect on natural resource availability and in turn on water quality. Historical records may not be sufficient to reliably account for uncertain future predictions under climate change conditions. While such highly uncertain situations become the “normal” case worldwide, the traditional practices of probabilistic risk measures cannot be used to appropriately quantify the uncertain phenomena under non-stationarity conditions. To better account for uncertain future conditions, the objective of this study is to develop a water management model based on Info-Gap Decision Theory (IGDT) using optimization under deep uncertainty conditions. The Info-Gap theory is a framework that measures the confidence in the operational decisions by quantifying the robustness to uncertainty without accounting for any probabilistic data. To demonstrate the method as a tool to better guide the long-term sustainable operation of the water supply system under uncertain future conditions, we applied the Info-Gap model to the Sea of Galilee (SoG) regional WSS, which is a significant part of the Israeli National Water System (INWS). For Israel, which is, like other Middle East semi-arid regions, prone to dry conditions and limited water availability, there are well-founded concerns that prolonged periods of drought lie ahead, as a consequence of the global climate change processes. This study contributes a management tool for decision making under deep uncertainty to improve the decision-making process and better adapt to unpredictable uncertain future conditions. We demonstrate how the IGDT could be formulated and used to analyze WSSs under different settings and demonstrate how decisions could be derived from the IGDT formulation. We also show a sensitivity analysis for the obtained solutions.

1. Introduction

The optimal planning and management of Water Supply Systems (WSSs) require the consideration of multiple aspects to assure that the needs of all stakeholders are satisfied [1]. The complexity and variety of system components that compose the WSS, the size of these components and the large sums of money involved in the regular operation, modification and future expansion of WSSs intensify the importance of optimal decision-making [2].
Traditionally, decision-making in WSSs deals with operation and investment decisions for spatial planning and construction efforts to utilize the available water resources in the most effective ways [3]. More recently, sustainability considerations were incorporated to assure that the trade-off between satisfying the increasing demands (i.e., population growth) and the availability of the adequate water resources will be kept in balance in the long term (for example, because of climate change, a decline in water quality in the sources, etc.) [1]. That is, we want to guarantee that the system will continue to function in the future.
Global climate change processes manifested by global warming effects and extreme events such as droughts and floods are observed around the world in increasing frequencies [2]. These factors challenge the accepted practices of the planning and management of WSSs. Climate change has a decisive effect on the availability of natural resources, thus there is a need to address uncertain climate effects in the decisions made today for the benefit of meeting sustainability goals in the future.
There are a variety of approaches to account for uncertainty in the planning and management of WSSs. A common way to deal with uncertainty is by a deterministic model, where the random variables are assumed to be known and predefined. As such, they could be replaced by values that represent the best estimates of the random variables. However, when the resulted problem parameters significantly deviate from the predicted values, replacing the random variables with poorly predicted values may lead to a solution that is not feasible with respect to the realized random variables [4].
To overcome the limited ability of deterministic models to reflect the real world uncertain behavior adequately, a wide range of stochastic modeling frameworks were developed to better account for different realizations of the uncertain factors [5]. In stochastic optimization models, the uncertain variables are replaced by random variables which are quantified using probabilistic or fuzzy logic models [6]. Both approaches use normalized mathematical functions to quantify uncertainty, but each aims at providing a different view of the incomplete information [1].
A fundamental aspect of stochastic programming is to cast the uncertain variables as Probability Density Functions (PDFs) or scenarios (discrete PDFs) [7]. When relevant data are extensive, a wide range of statistically based approaches (i.e., Monte-Carlo simulations) are available to construct the PDFs, or at least they can be reliably estimated. However, in many cases, the available data are not sufficient to reliably characterize a representative distribution [8]. In the water resources management domain, this is emphasized when the past occurrences of extreme events such as droughts or floods are rare, thus the data to construct the tails of the distribution are sparse, which may lead to unpredictable behavior [2]. Such highly uncertain situations are considered as Knightian (“deep” or “severe”) uncertainty, where the traditional practice of probability is confronted since the uncertainty cannot be quantified appropriately using PDFs [9].
In cases where a decision-maker is confident with the underlying probability distributions, a trustable optimal outcome can be achieved [10]. However, optimizing performance under nominal conditions will often be vulnerable when conditions depart from expectations [5]. Hence, relying on optimal predicted outcomes under Knightian uncertainty is unrealistic and may lead to poor decisions [3]. Ben-Haim [7] argues that in such cases when “true” optimality is complicated to achieve, it is preferred to look only for the outcomes that are critical and must be satisfied, despite the fact that these outcomes are not necessarily optimal. For such situations of severe uncertainty, other types of decision-making methods are available. These methods strive for the decision’s robustness to uncertainty instead of the decision’s optimality [11]. Identifying the minimum information required to satisfy a wide range of possible conditions will lead to the selection of a robust course of action [9,12]. More specifically, a robust option is one that performs well even when future conditions deviate from the predicted best estimate [6]. That is, a robust outcome reflects the immunity degree of a decision to uncertainty. Robustness to uncertainty analysis is a valuable tool as it enables the decision-maker to identify management options that perform acceptably well under a wide range of plausible future conditions [5].
In this study, we proposed the Info-Gap Decision Theory (IGDT), which is a quantified theory of robustness [6] that looks for the maximum robustness to failure instead of optimizing a predefined objective function. Unlike stochastic programming, IGDT is a non-probabilistic method which quantifies uncertainty using the concept of uncertainty sets, and attempts to optimize the robustness of decisions to uncertainty instead of optimizing an objective for a given probabilistic quantification of uncertainty [2,12]. The Info-Gap approach enables the decision-maker to investigate the impact of different decisions on the size of the uncertainty set, and therefore, to be able to quantify the robustness to failure of the system [1]. The size of the uncertainty set defines the cloud of possibilities consistent with the limited available information, where robustness is defined as the largest uncertainty set in which no failure will occur [10].
In the WSSs context, the severe lack of knowledge may be due to the unknown impact of global climate change trends, where future predictions are no longer realistic if taken based on past behavior. Generally, there is an agreement on a global warming trend [2], where models such as the IPCC global circulation models (IPCC 2007) predict incidences of extremely wet and extremely dry events. Specifically to the Middle East, there is a high likelihood of a dry scenario before the end of the 21st century, where models predict a 15% decrease in water availability in the region [13].
In Israel, the rainfall regimes in the last few decades reflect a consistent decreasing trend in average precipitation [14,15,16]. As Israel’s climate and hydrology patterns appear to be shifting toward dryer conditions with yet unfamiliar extent, we aim in this study to apply the Info-Gap Decision Theory to the Sea of Galilee (SoG) regional WSS, which is a significant part of the Israeli National Water System (INWS). Two weather-associated factors, rising temperatures and reduced rainfall, are the main drivers for the reduction in the SoG replenishment in recent years (see Section 2.2), where a 3.3% drop in rainfall per decade was observed during the last 40 years [14]. Thus, the SoG regional WSS management decisions will be investigated with IGDT in light of the uncertain behavior of future hydrological conditions by assessing the robustness of the system to various operation policies.
The remainder of this paper is organized as follows: Section 2 sets the background by briefly reviewing the available literature on IGDT applications for WSSs management problems under uncertainty, and gives an overview of the SoG WSS hydrological regime. Section 3 describes the concepts and structure of an Information Gap Decision Theory model. In Section 4, the formulation of the proposed regional WSS model is presented. The case study and the results are presented in Section 5. Finally, Section 6 concludes the study.

2. Background

2.1. Literature Review

In the past two decades, the Info-Gap decision theory gained a significant amount of applications in a variety of fields, such as engineering, economics, ecology, environmental protection, biology, medicine, project management, and national security [6,7]. However, among the multiple applications of IGDT, only a few were related to the domain of water resource systems, while only a handful addressed the planning and management of WSS problems, as we detail below.
Hipel and Ben-Haim [1] introduced IGDT to WSSs management problems by describing Info-Gap modeling to represent different sources of hydrological uncertainty. They illustrated decision making in watershed management problems using a few different Info-Gap models, each relevant to a different possible amount of prior information. They considered a lumped model of a WSS without modeling the spatial and temporal resolution of the problem (e.g., no conveyance network was considered). Woods et al. [11] used IGDT to assess the relative performance of various supply infrastructure portfolios meant to maintain the supply-demand balance in the Thames water resource system, UK. By comparing the robustness curves of the infrastructure options, they were able to prioritize the infrastructure options according to the projections of water availability, demand, and energy costs. A similar attempt to apply IGDT to quantitatively assess the robustness of various supply-side and demand-side management options was made by Korteling et al. [2] using a quasi-realistic water resources system based on a drift reservoir and demand node, located in southwest Cornwall, UK. In conjunction with IGDT, they included a Multi-Criteria Decision Analysis (MCDA) assessment for the management options to weight preferred performance results.
Matrosov et al. [9] demonstrated a joint application of Robust Decision Making (RDM) and IGDT frameworks to identify robust water supply portfolios for the planned capacity expansion of London’s water supply system in the Thames basin, UK. They showed that the methods are complementary since they provide different insights about the performance and vulnerabilities of the analyzed strategies. Matrosov et al. [9] concluded that a decision-maker would gain insights by deploying both frameworks in conjugation. Similar to Matrosov et al. [9], Roach et al. [3] compared the IGDT and Robust Optimization (RO) methods by assessing their applicability for a water resources management problem. They stressed that applying a wider range of decision methodologies will enable more encompassing investigation into potential future uncertainties.
In the broader context of water resources problems, IGDT was mainly applied to flood risk management problems. Hine and Hall [17] introduced IGDT to the field of flood management as an alternative to probability and fuzzy sets methods, to account for the uncertainty more adequately in highly uncertain fluvial rivers. They applied IGDT to a simple flood management problem to cost-effectively reduce the flood risk for surrounding residential properties. In their study, they evaluated the robustness of flood management decisions under uncertain upstream discharges into a floodplain. In a later work performed by Hine and Hall [5] for the river Trent catchment in the UK, an Info-Gap model was constructed around a parametric statistical flood discharges frequency analysis, resulting in a hybrid model of uncertainty. The robustness of a set of options for flood risk management accounted for two major sources of uncertainty: the flood depth model and the flow frequency. The analysis is considered a hybrid approach, in the sense that an Info-Gap model was used to analyze uncertainties in an underlying probabilistic model of the variability in fluvial flows. An additional application of Info-Gap analysis to a flood risk management problem was performed in [6], where few defense rising options have been considered to deal with potential damage caused by a given depth of floodwater. Therefore, the management decisions were related to the timing of interventions that should be taken. IGDT was tested as an alternative to probability and fuzzy sets theories while highlighting their unsuitability in cases of severe uncertainty.
In light of the above, there is a lack of studies that utilize the Info-Gap theory for the management of WSSs. The objective of this research is to develop a water management model based on the Info-Gap theory using optimization under uncertainty conditions to guide the long-term operation of the Sea of Galilee (SoG) WSS.

2.2. Sea of Galilee Hydrology and Climate Change Indications

The SoG watershed is divided between three Mediterranean countries, Israel, Lebanon, and Syria, and spans over 2730 km squared, where 25% of its area resides outside Israel’s borders [18]. It is situated in the northern part of Israel, which is, like other Middle East semi-arid regions, prone to dry conditions and limited water availability. The Sea of Galilee (Lake Kinneret) is the only natural surface freshwater reservoir in Israel and it is among the major water resources of the state [13,15]. The lake traditionally provided about a third of the country’s domestic, agricultural and industrial water requirements. Since the 1990s, extraction from Lake Kinneret exceeded the annual replenishment and its water levels declined to the lowest allowed level, endangering the ecological stability and the water quality of the lake [18]. In recent years, the Israeli government has undertaken extensive investments in water reclamation and desalination infrastructure in the country [14,15], while promoting water conservation. This has relieved the country’s reliance on the lake and significantly reduced the amount of water pumping. As a result, the lake’s ecological and environmental conditions have been improved, despite the persistent severe droughts in recent years [16]. Nowadays, the lake’s main uses are local water supply (mainly for agriculture), fishery, and tourism.
Climate change indicators are manifested in the SoG drainage basin. Using the Standard Precipitation Index (SPI), Givati et al. [19] indicated negative SPI (high level of aridity) of 11 and 17 during the 1927–1970 and 1970–2014 periods, respectively. This result indicates climate change as a higher level of aridity in the watershed. More specifically, the annual available water (inflows minus evaporation) in the SoG has decreased significantly during the last 40 years [14,15,16,19,20]. In Figure 1a, the linear regression line between the availability of SoG waters and Years (1976–2017) indicates a decline from 527 to 202 million cubic meters per year (MCM/Y). This result agrees with the conclusion of Gophen et al. [16] that indicated a decline from 470 to 225 MCM/Y in the period 1985–2016. Generally, there is a consensus among the Israeli government experts and researchers that climate change has an adverse effect on the SoG water level due to a combination of reduced precipitation and increased temperature/evaporation, which reduces the annual inflows. The flow in Israel’s streams and rivers is predominantly supplied by rainfall, where a general flow depletion trend has been measured across the SoG watershed over the past decade [16]. The extremeness of climate change in the watershed is highlighted when examining the Jordan River discharges during the last 40 years. Figure 3 in Gophen et al. [16] shows a decrease of about 53%. The increase in evaporation from the lake is associated with an increase in temperature in the watershed [16]. Figure 1b shows linear regressions between years and annual means of daily air maximum and minimum temperature measured in the KEFAR BLUM metrological station during 1976–2017. The annual maximum and minimum were elevated by about 2 °C during this period. This result agrees with Gophen et al. [16], which indicated increases of 2.7 °C and 1.5 °C in maximum and minimum temperature when analyzing data from DAFNA metrological station in the last 40 years.
These climate change indicators (Figure 1a,b) in the SoG watershed necessitate addressing severe uncertainty conditions in the planning and management of the SoG WSS. This is demonstrated by analyzing the recent years’ data and the SoG water level (Figure 1c), which shows five consecutive years of drought between years 2013–2018 [15]. The lake’s level reached −214.62 m in November 2018, which is defined as a black (catastrophic) line in the Israeli Water Authority (IWA). This drought was then followed by a wet winter (2019–2020) where the level of the lake increased by 5.73 m (a total volume increase of about 950 million cubic meters), leaving the lake almost full at level −208.89 m in May 2020. Although the water level of the lake did rise significantly, which ended the drought, the prospect for the coming future indicates that such extremism and high variability will challenge the water sector in the future [15,19]. Therefore, there are planned interventions in which desalination water will be “imported” from outside the SoG watershed to support the SoG WSS when needed.
The IWA employs a set of models for the planning and management of the regional WSS [13,21]. Although these models have been developed during the years to enable higher resolutions of various aspects of the system, support a multi-year time horizon, and account for both water quantity and quality, the models are essentially deterministic [21]. Attempts to account for uncertainty usually rely on a scenario-based approach to accommodate various possible realizations of the uncertain parameters. However, as explained previously, using deterministic models to make decisions under uncertainty might result in poor performance. Thus, in this study we present a decision support system based on the Info-Gap theory, which enables the decision-maker to come up with informed decisions by quantifying their confidence in the evaluated management policies through the evaluation of robustness to uncertainty [10,11]. By doing so, the Info-gap method allows the decision-maker to derive a trade-off between robustness and performance requirements. This trade-off enables the decision-maker to rank alternative policies based on their subjective attitude for balancing robustness and performance.

3. The Information Gap Decision Theory

A fundamental concept governing the IGDT approach is an expending uncertainty set relative to a nominal value (best estimate) of the uncertain parameter. When the uncertain parameter is equal to the nominal value, then no uncertainty is involved. Thus, as one step toward constructing the Info-Gap model, a nominal value needs to be obtained. A general representation of an optimization problem can be formulated as follows:
f ^ =   min x f x , θ s u b j e c t   t o
h i x , θ 0   i = 1 m
x 0
θ   U α
where x is the vector of decision variables, θ is the vector of uncertain parameters and f x , θ is the objective function. To achieve a feasible solution for the objective function, the constraints h i need to be fully satisfied. Note that we only include less-than-or-equal constraints because equality and larger-than constraints can be manipulated to be written as inequalities in a less-than-or-equal form.
When uncertainty is introduced, the uncertain parameter θ deviates from the nominal values θ ^ by an unknown amount α. The behavior of the uncertain parameter is described by an uncertainty set U α , which can be mathematically defined as follows:
U α = θ : θ θ ^ σ α α 0
In Equation (5), the deviation is set relative to the standard deviation σ of the uncertain parameters. That is, the expansion of the uncertainty set U α farther away from the nominal value θ ^ is relative to the standard deviation.
Using Equations (1)–(4), we can determine the Base Run (BR) value of the objective function by setting the uncertain parameters θ = θ ^ , that is, assuming there is no uncertainty involved:
f ^ B R =   min x f x , θ ^ s u b j e c t   t o
h i x , θ ^ 0   i = 1 m
x 0
As Info-Gap analysis begins with the best estimate of the uncertain parameter, and solving Equations (6)–(8) will obtain the objective function value at that nominal value. This base run value ( f ^ B R ) could be used as a reference to quantify the price-of-robustness.
Given a predefined B u d g e t , we can search for the maximum allowed deviation α that can be tolerated, whilst still leading to performance within requirements. Therefore, the Info-Gap optimization problem is to find a decision variable x which maximizes the allowed deviation α while still being feasible to all realizations in the uncertain set, as given in Equations (9)–(12).
α ^ =   max x α s u b j e c t   t o
max   θ   U α f x ,   θ B u d g e t
h i x ,   θ 0
x 0
where B u d g e t represents the threshold for deviation from the BR value. When no uncertainty is involved, the total annual operating cost is equal to the base run value of the objective function, hence when B u d g e t = f ^ B R we obtain the deterministic solution. When uncertainty is involved, increments of B u d g e t allow the decision-maker to explore different sizes of the uncertainty set by constructing a trade-off curve to represent the robustness of the system in different performance points. The robustness ( α ^ ) is the maximum allowed deviation α that still ensures the level of performance in which the operation cost is below B u d g e t .

4. Formulation of the Proposed Regional WSS Model

This research presents an annual multi-year model for the management of water quantity allocation and storage management using optimization methods. We focus on the management of a regional WSS under hydrological uncertainty associated with the replenishment of natural resources. The objective that drives the decision is to minimize the present value of WSS operation cost, subject to physical, operational, and environmental constraints.
To apply the Info-Gap model, three elements are combined to assess the robustness to uncertainty: the system model, the uncertainty model, and the performance criteria [7]. As the drop in average annual precipitation is a major factor that is likely to influence the integrity of the SoG WSS, in this study we focus on the uncertainty in the surface water, represented as a vector R t o t = [   R t o t t = 1 ,   , R t o t t = T f ] T , to investigate its influence on the system’s robustness. Under these settings, the robustness   α is a measure of the maximum amount of uncertainty in the surface water that the system can tolerate while satisfying the constraint on the maximum budget (i.e., the performance criteria).

4.1. System Model—Physical and Operational Constraints

4.1.1. Water Conservation Law

The topology of the network is represented by the junction node connectivity matrix, where A N 2 × M has a row for each node and a column for each edge. The construction of matrix A is demonstrated by a small hypothetical water supply system (Figure 2) that was used for model development and testing. The system is represented as a directed graph consisting of N nodes connected by M edges. The M edges represent the links between nodes; links in which the direction of flow is not fixed can be represented by two edges, one in each direction. The nodes are composed of three subgroups: N 1 are source nodes (e.g., desalination plant and reservoir) with one outgoing link from each source node; N 2 are intermediate nodes, where two or more edges meet; and N 3 are demand nodes, with one or more incoming links to each demand node.
The nonzero elements in each row are +1 and −1 for incoming and outgoing edges, respectively. The first two columns in A correspond to the links, which leave source nodes, while the last two columns correspond to the links ingoing demand nodes.
For each year t, the following system of linear equation ensures water conservation at the network nodes:
A Q t = 0   t
where Q t is the vector of discharges in the links (MCM/Y).

4.1.2. Limits on Water Quantities

The amount of water in all the system components (i.e., desalinated water from each plant, the discharges in the links, and extracted amount from the natural resources) is limited by an upper bound which represents maximum capacity, and by a lower bound (which may be zero or above) that represents operational constraints:
Q t m i n     Q t     Q t m a x   t

4.1.3. Hydrological Water Balance in the Sources

The hydrological water mass balance ensures that the change in storage equals the difference between the total recharge and withdrawal over the operation horizon:
V t =   V 0 + i = 1 t R i P i = 1 t Q i   t
where V t is a vector of the water volume (MCM) in the reservoirs at the end of year t; V 0 is a vector of the initial water volume (MCM) in the reservoirs; P is a permutation matrix that extracts the outgoing flows of the reservoirs from the vector Q t ; and R t is a vector of the reservoirs recharge (MCM/Y). The recharge vector of the reservoirs is parametrized based on the total recharge in the watershed, R t t o t , based on predetermined spatial distribution coefficients:
R t = β s · R t t o t   t
where β s is a vector of spatial distribution coefficients for water sources.
The water volume is bounded between a minimum and maximum water volume that represent physical conditions and/or operational constraints on the storage state:
V m i n     V t     V m a x   t
It is noteworthy that some of the sources may not be capable of storing water (e.g., pumping directly from a river stream). In such cases, one can set V m i n =   V m a x = V 0 = 0 to represent a direct inflow to the system without the ability to store water between years.

4.2. Objective Function

The objective is to operate the system with a minimum total cost composed of production cost, conveyance cost, and unsatisfied water demand cost (i.e., damage functions) over the operation horizon of T f years. These costs are modeled as nonlinear convex functions w.r.t the flows in the network. That is, at each link we have a cost function ϕ l Q t l .
Additionally, we account for Environmental Cost ( E C ), which is a convex penalty function related to the final state of the water storage in the reservoirs ( V T f ) at the end of the operation horizon. This penalty represents the sustainability goal in which we require to return the reservoir storage to an adequate level at the end of the operation horizon. Henceforth, all values of water quantities ( Q t ) are in Million Cubic Meters per year (MCM/Y) and the costs are in Million New Israeli Shekels (MNIS). The total cost, T C , is given as:
T C = t = 1 T f l = 1 n l i n k s ϕ l Q t l + E C V T f
For the sample network in Figure 2, an illustrative production/conveyance cost on link Q a is given in Figure 3a; an illustrative damage function on link Q z 1 is given in Figure 3b; and an illustrative Environmental Cost function is given in Figure 3c.

4.3. Deterministic LP Optimization Model

The cost functions, damage functions, and the environmental cost are represented in the model as nonlinear convex functions, thus leading to a nonlinear optimization problem. By approximating these nonlinear functions using convex piecewise-linear functions, the overall problem can be formulated as a Linear Programming (LP) model which could be efficiently solved by available optimization solvers (in this study we used LINPROG solver available within the MATLAB optimization toolbox).
Piecewise linearization is performed by dividing the nonlinear convex functions into linear segments, then, owing to the convexity, the value of the piecewise function could be obtained by applying the maximum operator overall segments. Thus, for a given convex function ϕ Q , an approximated convex piecewise-linear function ϕ ^ Q could be formulated as follows:
ϕ l Q t l ϕ ^ l Q t l = max i = 1 m l ( a i l · Q t l + b i l )   l t
where m l is the number of linear segments in link l and a i l , b i l are the coefficient of the linear segments. Figure 3d shows a convex piecewise-linear function with two segments. The figure shows that the linear function of the first segment obtains the maximum values in the range of the first segment (i.e., −2:4), while the linear function of the second segment obtains the maximum values in the range of the second segment (i.e., 4:10). Thus, we can write the piece-wise linear function in terms of a maximum operator as shown in Equation (19). The maximum operator in Equation (19) could then be formulated as an LP by introducing an auxiliary scalar variable φ t l , which represents the value of the piece-wise linear function. The equivalent LP can be formulated as:
m i n   φ t l s u b j e c t   t o
a i l · Q t l + b i l   φ t l
Applying convex piecewise-linear approximation to all the cost functions guarantees that the objective function and the constraints will be formulated as an LP. The overall deterministic (i.e., assuming R t t o t is known) is given as:
min t = 1 T f l = 1 nlinks φ t l   + φ EC s u b j e c t   t o
AQ t = 0   t
V t =   V 0 +   i = 1 t β s · R t tot P   i = 1 t Q i t
a i l · Q t l + b i l   φ t l   t           l         i = 1 m l
a i EC · V T f + b i EC   φ EC   i = 1 m EC
Q t min     Q t     Q t max   t
V min     V t     V max   t
The formulation above is a deterministic formulation in which the recharge is an assumed known parameter. In the next section, we show the Info-Gap formulation of the optimization problem above.

4.4. Info-Gap Model

As explained previously, to apply the Info-Gap model, three elements are combined to assess the robustness to uncertainty: the system model, the uncertainty model, and the performance criteria [7]. The system model and the performance criteria (budget) were presented in the previous section. In the following section, we describe the uncertainty model formulation.
For the uncertainty Info-Gap model, there is a wide range of commonly used Info-Gap models. The formulation and choice of an uncertainty model depend on the type and the quantity of available information [1,10]. In our case, the total available surface water ( R t t o t ) is considered uncertain and may vary largely among the years of the operation horizon. Therefore, the deviation α from the nominal value R ^ t t o t is measured relative to the standard deviation ( σ R ). As such, we consider a fractional-error Info-Gap model [8,9]:
U α = R t t o t   :     R t t o t 0 , R t t o t R ^ t t o t σ R   α           t   α 0
This model creates an expending interval around the nominal value R ^ t t o t :
U α = R t t o t   :   max 0 , R ^ t t o t α × σ R R t t o t R ^ t t o t + α   × σ R         t   α 0
where U α is a nested subset around a point estimate R ^ t t o t , thus R t t o t deviates from the nominal value by an unknown amount of uncertainty α . The family of nested subsets (Figure 4) is bounding the variation of the recharge around the nominal value. The greater the value of α , the greater the possibilities of R t t o t . However, in our case, the amount of available surface water must have a non-negative value. This implies that:
0     α     R ^ tot σ R
This work uses the average annual operation budget as a performance criterion of the system. Namely, the budget is used in the Info-Gap framework as a threshold to define the maximum achievable robustness. Figure 4 shows that the larger the budget, the larger the deviation that can be tolerated. More specifically, B u d g e t n > B u d g e t 2 > B u d g e t 1 lead to robustness (i.e., allowed deviation within budget) α n > α 2 > α 1 .
The Info-Gap framework seeks the maximum robustness for a specified budget, thus the optimization problem can be formulated as detailed in Equations (32)–(34) below. For ease of notation, the decision variables Q t   , V t   ,   φ t l   ,   φ E C are collected in a decision vector x and the set of constraints in Equations (23)–(28) is defined as a set X . Note that the set of constraints depends on R t o t because of the constraint in Equation (24). That is, a feasible decision vector should be x X R t o t .
max α s u b j e c t   t o
max   R tot   U α min   x   X R tot t = 1 T f l = 1 nlinks φ t l + φ EC   Budget
0     α     R ^ tot σ R
The first constraint in the optimization problem above implies that the worst optimum solution of the cost inside the uncertainty set should be less than B u d g e t . This means that the system cost will be guaranteed to yield a cost less than B u d g e t for any realization R t o t inside the uncertainty set. Chinneck and Ramadan [22] study the worst optimum solution of linear programs with interval uncertainty sets. In Theorem 6 of their work, they prove that the worst optimum solution will be obtained on an extreme scenario of the interval uncertainty sets. In our context, it is either R t o t = R ^ t o t σ R α or R t o t = R ^ t o t + σ R α . Since the optimization model in Section 4.3 focuses on water shortages costs (e.g., damage cost for unsatisfied demand, environmental cost of low water levels), the relevant extreme scenario for our problem is R t o t = R ^ t o t σ R α . In light of the above, the optimization problem in Equations (32)–(34) could be written as the following LP:
max x , α α s u b j e c t   t o
t = 1 T f l = 1 nlinks φ t l + φ EC   Budget
AQ t = 0   t
V t =   V 0 +   i = 1 t β s · R ^ tot σ R α P   i = 1 t Q i   t
a i l · Q t l + b i l   φ t l   t             l             i = 1 m l
a i EC · V T f + b i EC   φ EC   i = 1 m EC
Q t min     Q t     Q t max   t
V min     V t     V max   t
0     α     R ^ tot σ R

5. Case Study

In this section, the Info-Gap model described in Section 4.4 is applied to the optimal multi-year management of the SoG WSS. The water system layout is shown in Figure 5. It has one reservoir (Sea of Galilee), 9 surface sources, 7 ground water sources, 3 water imports (representing water transfer from national water supply system), 16 demand zones, and 59 links.
In the current study, we consider the hydrological uncertainty of the surface water in the watershed. The total surface water is considered uncertain according to the uncertainty set defined in Equation (29). The surface water is distributed between nine surface water sources according to predetermined spatial distribution coefficients, which are obtained based on the analysis of historical records. Thus, the surface water capacity is uncertain, while the other sources (i.e., ground water production and import) are deterministic and can produce any amount between pre-specified bounds. In terms of cost, the surface water is free, the ground water costs range between 0 and 1 NIS per m 3 , while the water import costs range between 1.5 and 2.5 NIS per m 3 . As such, import sources can produce a large and reliable amount of water, but at a greater cost. This operational trade-off between the deterministic-expensive vs. uncertain-free sources will be examined in the results.
For the next 3 years, we consider 30 equally likely hydrological scenarios of the total surface water in the watershed, as shown in Figure 6. These forecasts aim to emulate different possible future scenarios, ranging from low precipitation series to high precipitation series. The nominal scenario (i.e., the scenario with average total inflow) is highlighted in green in Figure 6. This nominal scenario is used to construct the uncertainty set of the Info-Gap model as detailed in Equation (29), where σ R = 150 MCM.
The development and implementation of the mathematical model were programmed using MATLAB and the YALMIP toolbox [23]. YALMIP is used to model and solve the optimization model by interfacing external solvers via simple mathematical programming syntax. In this study, the optimization problem was solved using the LP solver, LINPROG, which is provided within the optimization toolbox of MATLAB.

5.1. Operation Policies

An important aspect of solving uncertain models is the sequence in which decisions alternate with uncertain observations [4]. That is, which decisions should be taken before uncertainty is revealed and which decisions are recourse decisions that adapt to the revealed uncertainty. In our context, we must distinguish between the operation policy that the decision-maker determines before uncertainty is revealed and recourse decisions that adapt to the revealed uncertainty. For example, if we consider one demand zone D Z which is fed by one ground water source G W and one uncertain surface water source S W (Figure 7), we cannot set both Q G W and Q D Z since the problem will be infeasible for all realizations of the uncertain source, except R = Q D Z Q G W . A better modeling approach is to set one of the decisions (e.g., Q G W ) as an operation policy that should be determined by the decision-maker before uncertainty is revealed, and the other one (e.g., Q D Z ) as a recourse decision that takes value depending on the observed R and the determined operation policy (i.e., Q D Z = R + Q G W ).
The operation policy represents the immediate commitment of the decision-maker that should be made in face of the uncertain future, while the recourse decisions are decisions that represent the reaction of the system to the revealed future scenario. In our case study, the WSS operation is solved for three years operation horizon, where the operation policy that represents the immediate commitment of the decision-maker is comprised of the ground water productions (seven decision variables) and the water imports (three decision variables) in the first year. In total, an operation policy is comprised of determining 10 decision variables in the first year of the operation horizon. These 10 variables could be determined using deterministic optimization or optimization under uncertainty methods. In this study, we examine three operational policies and show their robustness based on the Info-Gap framework.
First, we consider the Nominal Policy (NP), which is obtained by solving a deterministic optimization problem with the nominal scenario highlighted in Figure 6. In the NP, the decision-maker assumes that the nominal scenario will be realized and commits to ground water production and water import during the first year according to this assumption. Second, we consider the Scenario Based Policy (SBP), in which the decision-maker accounts for all 30 forecast scenarios in Figure 6. In this case, as the scenarios are equally likely, the decision-maker solves each of the scenarios and then commits to ground water production and water import during the first year based on the average amounts obtained from the solutions of the 30 scenarios. Third, we consider the Robust Policy (RP), which is obtained by solving the Info-Gap model Equations (35)–(43) with a budget of 400 MNIS that guarantees the robustness of 1.5, as will be explained in the next section. Figure 8 shows the obtained operation policies (i.e., ground water productions and water imports in the first year) using the three approaches discussed above.
As shown in Figure 8, the NP relies on ground water only, without the need to utilize expensive imported water. That is, if the nominal scenario is realized, the system can rely on local sources without the need for water transfer from the national water system. The SBP shows more ground water withdrawal as well as the utilization of water imports. This is because some of the scenarios considered (see Figure 6) are dry scenarios that require more ground water withdrawal and transferring water from the national water system. The different utilization of the sources (e.g., high amount in I m 2 compared to I m 1 and I m 3 ) is attributed to the different production and conveyance costs in the system. The RP exhibits the largest utilization of all ground water sources and water imports. This is because the RP, which is designed to guarantee a robustness of 1.5, accounts for a dry scenario that deviates from the nominal scenario by as high as 1.5 × σ R = 225 MCM.

5.2. Operational Policies Comparison

The robustness of the three operation policies discussed in the previous section is analyzed using the Info-Gap framework. The Info-Gap formulation in Equations (35)–(43) is applied using a two-stage recourse approach, where the operation policies are implemented at the beginning of the operation horizon before the uncertainty is revealed, and the other decisions are recourse decisions that are determined within the Info-Gap formulation. More specifically, each of the operation policies will impose equality constraints on the Info-Gap formulation by setting the 10 decision variables of the operation policy to the values shown in Figure 8. For example, to evaluate the robustness of the NP, the Info-Gap formulation in Equations (35)–(43) is solved with additional constraints that set the first-year ground water productions and water imports to the values shown in red in Figure 8.
Figure 9a,b shows the trade-off between robustness and annual budget for each of the three operation policies. These trade-offs are obtained by solving the Info-Gap formulation in Equations (35)–(43) for increasing series of annual budgets for each of the three operation policies. For each point in the robustness trade-off, we have the decision vector corresponding to that point. Namely, for any given budget the optimization problem provides the optimal decisions for the entire three years operation horizon that maximize the robustness of the system operation under the operation policy constraints.
For example, the trade-off of the NP shows robustness vs. budget in case the decision-maker relies on the nominal scenario to dictate the productions and water import in the first year. The different points on the trade-off show the maximum deviations (in standard deviations units) from the nominal surface water amount, which can be tolerated under different budgets. For example, if the NP is implemented and the annual budget is 238 MNIS, then the robustness is zero. This indicates that we cannot tolerate any uncertainty with a budget of 238 MNIS under the NP, while to tolerate a deviation of 1.5 × σ R = 225 MCM in the NP, an annual budget of 520 MNIS is required. The results also show that the NP can reach a maximum robustness of 2.895 (i.e., tolerate a deviation of 2.895 × σ R = 434.25 MCM) with an annual budget of 1500 MNIS.
Comparing the different trade-offs from the three operation policies shows that the RP is superior for budgets larger than 307 MNIS. This is demonstrated by higher robustness for the given budget obtained by the RP. For instance, when the budget is 800 MNIS the RP can reach a robustness of 2.451, while the SBP and the NP can only reach a robustness of 2.183 and 2.068, respectively. Another way to interpret the results is that the RP can tolerate a deviation of 1.5 × σ R = 225 MCM while guaranteeing a budget below 400 MNIS, while the SBP and NP will require budgets of 480 MNIS and 520 MNIS, respectively, to tolerate the same deviation.
The maximum robustness in all operation policies is 2.895. However, the RP reaches this robustness with a smaller budget of 1100 MNIS as opposed to 1400 MNIS and 1500 MNIS in the SBP and NP, respectively. To conclude, the results show that for a budget larger than 307 MNIS, the RP can cope better with uncertainty by maximizing the uncertainty horizon that can be tolerated for a given budget. Comparing the SBP with the NP shows that the SBP can tolerate larger uncertainty. This is because the SBP operation policy is derived from 30 scenarios that include dry scenarios, which result in a more robust operation policy compared to the NP.
Figure 9b shows the trade-off for low budgets. The results show that for a low budget (budget < 307 MNIS), the RP has less robustness, that is, it can tolerate less uncertainty compared to the SBP and NP. In other words, more budget is required to tolerate a small level of uncertainty in the RP. For example, to tolerate a deviation of 0.4 × σ R , the SBP and the NP can guarantee a budget below 262 MNIS and 268 MNIS, respectively, whilst the RP will have a higher budget of 285 MNIS. This result can be explained by the fact that RP is a worst-case oriented policy that is expected to perform optimally in extreme scenarios, while SBP and NP are expectation-oriented policies that are expected to perform optimally around the nominal scenarios. Nonetheless, for a robustness of 0.4, the additional budget in the RP compared to the NP is only 17 MNIS, while for a high robustness of 1.5 the NP will have a high additional budget of 100 MNIS compared to the RP. As such, for risk aversion, the decision-maker may prefer the RP, which behaves optimally in extreme scenarios.

5.3. Examining the Operation Policies under Monte-Carlo Simulations

For water resources management models, we are often able to generate a large number of scenarios that form an approximation of the underlying random process. These sets of scenarios may be obtained by simulation using stochastic models, historical data, and expert forecasting for extreme scenarios. One of the most common ways of generating individual scenarios is Monte-Carlo Sampling [2]. Using this approach, random variables are generated from their known (or assumed) distribution, thus replacing continuous probability space with a finite space consisting of the samples. The generated samples are then used to simulate the possible future values of the uncertain variables [5].
We would like to emphasize that the IGDT does not require Monte-Carlo simulation or any other probabilistic assumptions. In this section, we only use Monte-Carlo simulations to evaluate the different operation policies and compare the results to the budget guarantee provided by the IGDT. As indicated previously, in our case study we focus on dry scenarios and the resulting water shortages. Thus, we have generated 1000 uniformly random samples between R t t o t 1.5 × σ R and R ^ t t o t for the three years operation horizon. The SoG WSS model is solved for each of the 1000 samples under the three operation policies to obtain the budget histograms, as shown in Figure 10.
The results obtained by the Monte-Carlo simulations (Figure 10) show that the RP guarantees a budget below 400 MNIS for deviation up to 1.5 × σ R . This budget guarantee agrees with the IGDT analysis shown in Figure 9b for a robustness level of 1.5.
As discussed in the previous section, the NP and SBP are less robust compared to the RP under extreme scenarios. As such, both of these policies exceed the 400 MNIS budget guarantee of the RP, as shown in Figure 10b,c. That is, both require a budget greater than 400 MNIS to cope with a deviation of 1.5 × σ R . The budget violation in the NP and SBP relative to the RP budget guarantee is 25% and 10%, respectively. On the other hand, the average budgets in the histograms are 314, 303, and 308 for the NP, the SBP, and the RP, respectively, showing that the SBP performs better on average when compared to the RP. Nonetheless, the deviation in the average budget is (308 − 303)/303 = 1.7% compared to a budget guarantee reduction of (440 − 400)/400 = 10% for extreme scenarios. The empirical results obtained from the Monte-Carlo simulations analysis fortify the advantages of the RP. It is noteworthy that investigating the system robustness via the IGDT is more computationally efficient, since the budget required for each robustness value is obtained through the trade-off curve (e.g., Figure 9a) without the need to empirically rely on thousands of random simulations and assumptions on the probabilistic distribution of the random data.

5.4. Sensitivity Analysis

The IGDT robustness evaluation can also facilitate sensitivity analysis for important parameters in the model. For example, an important parameter in the SoG WSS is the minimum flow in link L30, which represents the minimum flow needed for recreational purposes and the natural conservation of the Upper Jordan River. In recent years, the annual average flow of the Upper Jordan River has dropped to roughly half of the annual average flow, which was typically over 100 MCM [24].
To examine the impact of minimum flow regulation, a sensitivity analysis was performed by running the Info-Gap model for different levels of Q l 30 m i n . Requiring minimum flow in link L30 implies that less surface water should be used in the upper part of the system in order to meet the minimum flow constraint. As such, the system is forced to rely on more expensive imported water and increased ground water production. This increases the operation cost, which leaves less budget to tolerate the uncertainty (i.e., reduced robustness).
Figure 11 shows the robustness vs. budget trade-off of the RP for different levels of minimum flow regulation in link L30. The comparison of the obtained trade-offs shows that the robustness vs. budget trade-off shrinks when the minimum flow regulation is stricter (i.e., Q l 30 m i n is higher). For example, the RP can tolerate a deviation of 2.895 × σ R = 434.25 MCM with a budget of 1100 MNIS when there is no minimum flow requirement, but it can only tolerate a deviation of 1.463 × σ R = 219.45 MCM for the same budget when a minimum flow regulation of 200 MCM is imposed on link L30. An important aspect of the results shown in Figure 11 is the overlapping trade-offs in the budget range of 279–350 MCM, which indicates that the minimum flow regulation in link L30 has little impact on the robustness within this budget range. Another important observation is that the length of the flat segments in the curves is increasing for increased minimum flow regulation in link L30. This indicates that for higher minimum flow regulation, the budget has less impact on the trade-off (i.e., the robustness is fixed regardless of the budget amount). In other words, the budget becomes a non-limiting factor in determining the robustness as the minimum flow requirement increases (i.e., other factors such as conveyance and production capacities become limiting factors).
We emphasize that the above is just a demonstration of how the IGDT could be used for sensitivity analysis. Many other sensitivity analyses might be performed to explore the impact on the system robustness. Future studies may use the IGDT for sensitivity analysis to explore other impacts on the system robustness, such as topological changes (e.g., adding new links or exploring planning alternatives).

6. Discussion and Conclusions

In this study, we developed a water management model based on the IGDT to investigate climate change impacts on WSS management decisions. We focused on the impact of climate change related to hydrological uncertainty, where only the replenishment into the natural resources is taken as an uncertain parameter. To incorporate explicit considerations of uncertainty, an appropriate decision-making framework under uncertainty should be selected such that the real nature of the problem will be reflected during the modeling process, and that the results are usable in the real world. The tools required to support decision making in the face of uncertainty must be flexible enough to allow the testing of trade-offs between various characteristics of future outcomes, efficient enough to run quickly for different inputs (i.e., computational efficiency), and transparent enough to convey to the decision-makers the full range of consequences of different possible decisions. We demonstrated the IGDT as a tool that preserves the properties above for the multi-year SoG WSS management problem.
Our results demonstrate the importance of introducing uncertainty into deterministic models when evaluating climate change impacts, as the solution can change dramatically when recharge uncertainty is considered. In cases of highly uncertain future conditions, the uncertainty could be represented in the Info-Gap model using the concept of uncertainty sets. We considered the sequence in which decisions alternate with observations and showed that in cases where all decisions are taken before the uncertainty is revealed, the problem may become infeasible for most of the realizations of the uncertain source. Thus, for policy selection, we implemented a two-stage recourse framework. In the first stage, some decisions (i.e., an operation policy) are implemented at the beginning of the operation horizon before the uncertainty is revealed, and this represents the immediate commitment made, while in the second stage recourse decisions are allowed to adapt to the revealed uncertainty.
The Info-Gap model enables the decision-maker to examine the trade-off between the amounts of robustness gained for every budget limit. That is, for predetermined operation policies, the performance of each policy is examined by searching for the robust optimal decision at each budget limit. In this study, we examined three operation policies and showed their robustness based on the Info-Gap framework. An important aspect related to the investigated SoG WSS was to examine the trade-off between the deterministic-expensive vs. uncertain-free sources. That is, to evaluate the impact of two available water sources with different characteristics on the system robustness, where the import sources (i.e., desalination water) can produce a large and reliable amount of water, but at a greater cost, compared to the local productions which are influenced by uncertainty. The results reveal that the operation policies show different preferences between imported water and relying on local productions. The NP relies mostly on the water drawn from the local production sources, whereas policies with higher robustness (i.e., SBP and RP) rely on increased amounts of imported water. In general, the NP and SBP are less robust compared to the RP under extreme scenarios.
From a budget perspective, a swap between the trade-off curves indicates that some operation policies may perform better only in part of the budget range. This is obtained by the RP, which demonstrated superior robustness starting from a given budget, namely, the RP can cope better with uncertainty by maximizing the uncertainty horizon that can be tolerated for a given budget. Likewise, as the SBP is a more conservative operation policy compared to the NP, the SBP can tolerate larger uncertainty, and thus it is more robust than NP. In conclusion, as the RP is a worst-case oriented policy, it is expected to perform optimally in extreme scenarios, while the SBP and NP are expectation-oriented policies that are expected to perform optimally in nominal scenarios. Therefore, a risk-averse decision-maker may prefer the RP to maximize the uncertainty that the WSS can tolerate in extreme hydrological scenarios.
A future extension to our Info-Gap model should include water quality modeling. The combined management of the quality and the quantity may change the optimal solution compared to quantity considerations alone. Moreover, we used a relatively simple interval uncertainty set to demonstrate the methodology. Choosing other shapes of the uncertainty set, such as ellipsoidal shape uncertainty sets, may be advantageous, as suggested in Housh et al. [25]. Thus, testing different uncertainty sets within the Info-Gap methodology is a matter for future research. Moreover, in our model, we considered the sustainability of the SoG WSS operation decisions by assigning a penalty value on the system state at the end of the operation horizon. More specifically, a penalty is incurred for leaving the system at the end of the planning horizon with a reservoir storage below a predetermined level. The future consideration of sustainability in the management plan should include a perspective with a relatively long time-horizon to assure that the SoG reservoir will not deplete, nor have high salinity at the end of the planning horizon.

Author Contributions

Conceptualization, M.H.; Methodology, T.A. and M.H.; Software, T.A. and M.H.; Validation, M.H.; Formal Analysis, T.A. and M.H.; Investigation, T.A. and M.H.; Data Curation, T.A.; Writing—original draft, T.A. and M.H.; Writing—review & editing, M.H.; Visualization, T.A.; Supervision, M.H.; Project Administration, M.H.; Funding Acquisition, M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Israeli Water Authority, grant number 4501683698 and The APC was funded by the Center for Cyber Law & Policy at the University of Haifa.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hipel, K.W.; Ben-Haim, Y. Decision making in an uncertain world: Information-gap modeling in water resources management. IEEE Trans. Syst. Man Cybern. Part C Appl. Rev. 1999, 29, 506–517. [Google Scholar] [CrossRef]
  2. Korteling, B.; Dessai, S.; Kapelan, Z. Using Information-Gap Decision Theory for Water Resources Planning under Severe Uncertainty. Water Resour. Manag. 2013, 27, 1149–1172. [Google Scholar] [CrossRef] [Green Version]
  3. Roach, T.; Kapelan, Z.; Ledbetter, R. Comparison of Info-gap and Robust Optimisation Methods for Integrated Water Resource Management under Severe Uncertainty. Procedia Eng. 2015, 119, 874–883. [Google Scholar] [CrossRef] [Green Version]
  4. Sen, S.; Higle, J. An Introductory Tutorial on Stochastic Linear Programming Models. Interfaces 1999, 29, 33–61. [Google Scholar] [CrossRef] [Green Version]
  5. Hine, D.; Hall, J.W. Information gap analysis of flood model uncertainties and regional frequency analysis. Water Resour. Res. 2010, 46, 1–18. [Google Scholar] [CrossRef]
  6. Hall, J.W.; Harvey, D. Decision making under severe uncertainties for flood risk management: A case study of info-gap robustness analysis. In Proceedings of the 8th International 25/25 Introduction Program Voodoostan Academistan Publicistan Conclusion Conference on Hydroinformatics, Concepcion, Chile, 12–16 January 2009. [Google Scholar]
  7. Ben-Haim, Y. Info-Gap Decision Theory (IG). In Decision Making under Deep Uncertainty; Marchau, V., Walker, W., Bloemen, P., Popper, S., Eds.; Springer: Cham, Switzerland, 2019. [Google Scholar]
  8. Hayes, K.R.; Barry, S.C.; Hosack, G.; Peters, G. Severe uncertainty and info-gap decision theory. Methods Ecol. Evol. 2013, 4, 601–611. [Google Scholar] [CrossRef]
  9. Matrosov, S.E.; Woods, M.A.; Harou, J.J. Robust Decision Making and Info-Gap Decision Theory for water resource system planning. J. Hydrol. 2013, 494, 43–58. [Google Scholar] [CrossRef]
  10. Ben-Haim, Y.; Demertzis, M. Decision making in times of Knightian uncertainty: An Info-Gap perspective. Econ. Open Access Open Assess. E J. 2016, 10, 1–29. [Google Scholar] [CrossRef] [Green Version]
  11. Woods, A.; Matrosov, E.; Harou, J.J. Applying Info-Gap decision theory to water supply system planning: Application to the Thames basin. In Proceedings of the Computer Control in the Water Industry (CCWI), Exeter, UK, 5–7 September 2011. [Google Scholar]
  12. Ben-Haim, Y. Info-Gap Decision Theory: Decisions under Severe Uncertainty, 2nd ed.; Academic Press: London, UK, 2006. [Google Scholar]
  13. Sade, R.; Rimmer, A.; Samuels, R.; Salingar, Y.; Denisyuk, M.; Alpert, P. Water Management in a Complex Hydrological Basin—Application of Water Evaluation and Planning Tool (WEAP) to the Lake Kinneret Watershed, Israel. In Integrated Water Resources Management: Concept, Research and Implementation; Borchardt, D., Bogardi, J., Ibisch, R., Eds.; Springer: Cham, Switzerland, 2016; pp. 35–57. [Google Scholar]
  14. Tal, A. The implications of climate change driven depletion of Lake Kinneret water levels: The compelling case for climate change-triggered precipitation impact on Lake Kinneret’s low water levels. Sci. Total Environ. 2019, 664, 1045–1051. [Google Scholar] [CrossRef] [PubMed]
  15. Gophen, M. Climate and Water Balance Changes in the Kinneret Watershed: A Review. Open J. Mod. Hydrol. 2020, 10, 21–29. [Google Scholar] [CrossRef] [Green Version]
  16. Gophen, M.; Meron, M.; Levin-Orlov, V.; Tsipris, Y.; Peres, M. Climate Change, Regional Water Balance and Land Use Policy, in the Watershed of Lake Kinneret (Israel). Open J. Ecol. 2020, 10, 200–224. [Google Scholar] [CrossRef]
  17. Hine, D.J.; Hall, J.W. Convex analysis of flood inundation model uncertainties and Info-Gap flood management decisions. In Proceedings of the ISSH—Stochastic Hydraulics, Nijmegen, The Netherlands, 23–24 May 2005. [Google Scholar]
  18. Markel, D. Monitoring and management Lake Kinneret (Sea of Galilee)–preserving Israel’s major surface water resource. EnviroGeoChimica Acta 2014, 1, 411–420. [Google Scholar]
  19. Givati, A.; Guillaume, T.; Rosenfeld, D.; Paz, D. Climate Change Impacts on Streamflow at the Upper Jordan River Based on an Ensemble of Regional Climate Models. J. Hydrol. Reg. Stud. 2019, 21, 92–109. [Google Scholar] [CrossRef]
  20. Givati, A.; Rosenfeld, D. Possible Impacts of Anthropogenic Aerosols on Water Resources of the Jordan River and the Sea of Galilee. Water Resour. Res. 2007, 43, W10419. [Google Scholar] [CrossRef]
  21. IWA. Determine Operation Policies for SoG WSS and Investigation of Planning Alternatives for Water Import, Using a Simplified Model—Final Report; Israel Water Authority, Operation Department: Tel-Aviv, Israel, 2019. (In Hebrew) [Google Scholar]
  22. Chinneck, J.; Ramadan, K. Linear programming with interval coefficients. J. Oper. Res. Soc. 2000, 51, 209–220. [Google Scholar] [CrossRef]
  23. Lofberg, J. YALMIP: A toolbox for modeling and optimization in MATLAB. In Proceedings of the 2004 IEEE International Conference on Robotics and Automation (IEEE Cat. No.04CH37508), Taipei, Taiwan, 2–4 September 2004; pp. 284–289. [Google Scholar]
  24. Tal, A. Climate change’s impact on Lake Kinneret: Letting the data tell the story. Sci. Total Environ. 2019, 685, 1272–1275. [Google Scholar] [CrossRef] [PubMed]
  25. Housh, M. Optimal Multi-Year Management of Regional Water Resources Systems under Uncertainty. Ph.D. Thesis, Technion, Haifa, Israel, 2011. [Google Scholar]
Figure 1. Climate change indicators: (a) Linear regression between availability of Sea of Galilee (SoG) waters and Years (1976–2017); (b) Linear regressions between years and annual means of daily air maximum and minimum temperature measured in KEFAR BLUM metrological station during 1976–2017; (c) Drop in SoG water levels over the past decade.
Figure 1. Climate change indicators: (a) Linear regression between availability of Sea of Galilee (SoG) waters and Years (1976–2017); (b) Linear regressions between years and annual means of daily air maximum and minimum temperature measured in KEFAR BLUM metrological station during 1976–2017; (c) Drop in SoG water levels over the past decade.
Sustainability 13 03152 g001aSustainability 13 03152 g001b
Figure 2. A small hypothetical water supply system (right) and the junction node connectivity matrix ( A ) of the network (left).
Figure 2. A small hypothetical water supply system (right) and the junction node connectivity matrix ( A ) of the network (left).
Sustainability 13 03152 g002
Figure 3. Cost functions for the sample network: (a) Production/conveyance cost on link Q a ; (b) Damage function on link Q z 1 ; (c) Environmental cost function in the reservoir; (d) Illustration of convex piecewise-linear function.
Figure 3. Cost functions for the sample network: (a) Production/conveyance cost on link Q a ; (b) Damage function on link Q z 1 ; (c) Environmental cost function in the reservoir; (d) Illustration of convex piecewise-linear function.
Sustainability 13 03152 g003aSustainability 13 03152 g003b
Figure 4. Nested subsets U α around a point estimate R ^ t t o t , where R t t o t deviates from R ^ t t o t by an unknown amount of uncertainty α.
Figure 4. Nested subsets U α around a point estimate R ^ t t o t , where R t t o t deviates from R ^ t t o t by an unknown amount of uncertainty α.
Sustainability 13 03152 g004
Figure 5. Sea of Galilee WSS layout.
Figure 5. Sea of Galilee WSS layout.
Sustainability 13 03152 g005
Figure 6. Thirty equally likely hydrological scenarios of the total surface water in the SoG watershed.
Figure 6. Thirty equally likely hydrological scenarios of the total surface water in the SoG watershed.
Sustainability 13 03152 g006
Figure 7. Example of one demand zone which is fed by one ground water source and one uncertain surface water source.
Figure 7. Example of one demand zone which is fed by one ground water source and one uncertain surface water source.
Sustainability 13 03152 g007
Figure 8. Flows of the examined operation polices for the SoG WSS. Note: NP = Nominal Policy; SBP = Scenario Based Policy; RP = Robust Policy.
Figure 8. Flows of the examined operation polices for the SoG WSS. Note: NP = Nominal Policy; SBP = Scenario Based Policy; RP = Robust Policy.
Sustainability 13 03152 g008
Figure 9. (a) Robustness vs. budgets trade-offs for the operation policies. (b) Zoom in for the operation policies’ robustness vs. budgets trade-offs.
Figure 9. (a) Robustness vs. budgets trade-offs for the operation policies. (b) Zoom in for the operation policies’ robustness vs. budgets trade-offs.
Sustainability 13 03152 g009
Figure 10. Distributions of operation budget obtained by Monte-Carlo Simulations for: (a) RP, (b) SBP and (c) NP.
Figure 10. Distributions of operation budget obtained by Monte-Carlo Simulations for: (a) RP, (b) SBP and (c) NP.
Sustainability 13 03152 g010aSustainability 13 03152 g010b
Figure 11. Robustness vs. budget trade-offs under RP, for different minimum flow regulations in link L30.
Figure 11. Robustness vs. budget trade-offs under RP, for different minimum flow regulations in link L30.
Sustainability 13 03152 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Housh, M.; Aharon, T. Info-Gap Models for Optimal Multi-Year Management of Regional Water Resources Systems under Uncertainty. Sustainability 2021, 13, 3152. https://doi.org/10.3390/su13063152

AMA Style

Housh M, Aharon T. Info-Gap Models for Optimal Multi-Year Management of Regional Water Resources Systems under Uncertainty. Sustainability. 2021; 13(6):3152. https://doi.org/10.3390/su13063152

Chicago/Turabian Style

Housh, Mashor, and Tomer Aharon. 2021. "Info-Gap Models for Optimal Multi-Year Management of Regional Water Resources Systems under Uncertainty" Sustainability 13, no. 6: 3152. https://doi.org/10.3390/su13063152

APA Style

Housh, M., & Aharon, T. (2021). Info-Gap Models for Optimal Multi-Year Management of Regional Water Resources Systems under Uncertainty. Sustainability, 13(6), 3152. https://doi.org/10.3390/su13063152

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