Solution Approaches for the Management of the Water Resources in Irrigation Water Systems with Fuzzy Costs

Currently, the management of water networks is key to increase their sustainability. This fact implies that water managers have to develop tools that ease the decision-making process in order to improve the efficiency of irrigation networks, as well as their exploitation costs. The present research proposes a mathematical programming model to optimize the selection of the water sources and the volume over time in water networks, minimizing the operation costs as a function of the water demand and the reservoir capacity. The model, which is based on fuzzy methods, improves the evaluation performed by water managers when they have to decide about the acquisition of the water resources under uncertain costs. Different fuzzy solution approaches have been applied and assessed in terms of model complexity and computational efficiency, showing the solution accomplished for each one. A comparison between different methods was applied in a real water network, reaching a 20% total cost reduction for the best solution.


Sustainability Concept in Water Systems
The use of water resources is crucial in irrigation water systems, and currently, their management is also of the utmost importance since their sustainability must be improved by developing an Integrated Water Resources Management (IWRM). The work in [1] defined this as "a process, which promotes the coordinated development and management of water, land and related resources in order to maximize the resultant economic and social welfare in an equitable manner without compromising the sustainability of vital ecosystems". The IWRM requires the application of Adaptative Water Resources Management (AWRM), which is a process to improve the water policies and practices by applying the experience from previous years [2]. The need for increased sustainability and resilience, as well as the improvement of drought management results in water managers needing to develop or use mechanisms that support, make more resilient, and ease the decision-making process [3]. Different research works have analyzed the management of water sources in a basin, defining the optimized water level, the volume transfer between rivers, and/or the maximum farmer demand in the case of drought [4,5]. However, there are not many research works addressing the analysis and management of water networks when the system supplies an irrigation community. Such irrigation systems usually have different water sources. For this reason, the development of supporting tools is critical to apply optimization strategies that reduce the exploitation costs and ensure the required demand is met. Currently, the decision about the water source is usually done in a traditional manner, and it is not focused on the optimization of economic costs.
The modernization of irrigation systems supposes the use of new efficient techniques of cultivation with a positive margin where the climatic conditions are not favorable (e.g., long periods of drought, high temperatures, or other factors). However, the modernization of irrigation systems also presents some disadvantages. The main setback of this type of system is that the installed power should be increased by 2 kW/ha to guarantee the pressure in the irrigation system, and therefore, this increase directly causes the increase of energy costs and reduces the sustainability indicators [6].
Currently, the relationship between energy and water is evolving as we search for a new concept of sustainability [7]. Pressurized irrigation systems demand more energy, and this is currently one of the main uses of water linked to energy (49%) together with hydroelectric production (30.6%) and the highest (92%) in terms of water consumption. Furthermore, in recent years, the use of energy in this sector has grown rapidly [8].

Political Context
The growing demand for water joined with increasing water scarcity [9] has led to increasing interest in global modeling of water resources [10]. There are several factors that influence the supply and industrial and agricultural water use. These factors are characterized by an uncertain context (e.g., climatic conditions and energy prices, among others). This requires that decision-makers be resilient enough to deal with uncertainty in an efficient manner [11].
Based on the review [12], the United Nations Conference on Water held in Argentina in 1977 defined a set of objectives to work on before the end of the 20 th Century. One of these main objectives was to ensure an adequate supply of quality water to meet the planet's socioeconomic needs [12]. Afterwards, during the International Conference on Water and Environment, held in Dublin in 1992, it was agreed that water should be considered as an economic good, a principle that has been greatly discussed during and after this event. At the International Conference on Freshwater held in 2001 in Bonn, a set of actions necessary to mobilize financial resources was identified, and among these actions, the improvement of economic water-use efficiency was defined as a priority. Based on the recent historical milestones of water management, the definition of efficient water resource management strategies that administer water in a sustainable and resilient direction seems necessary. These strategies should not only focus on basin management; they must consider the distribution in the water systems.
In particular, the Mediterranean is a region characterized by water scarcity and increasing energy demand [13]. The irrigation systems need energy to be: (i) originated from them, (ii) distributed to the users, and (iii) regenerated during waste-water treatment [14]. These interdependencies across water, energy, and food define the Water-Energy-Food (WEF) nexus that highlights the importance of the interactions among these sectors [15]. This study is focused on the WEF nexus, as it is addressed to the optimization of the provision of water to be used in irrigation systems to produce food, taking into account the energy costs, in particular the extraction costs, which depend on the water's source.
The literature shows different water pricing approaches to obtain the appropriate prices to allocate water efficiently [16,17]. Every day, water becomes increasingly vital as demand for food and water increases and water scarcity becomes a reality [16]. In Spain, the relationship between water and energy has reached great significance as a consequence of the increase in energy prices since 2007. Moreover, the modernization of irrigation systems has led to a reduction in water consumption per area, but the increase in hydraulic efficiency has also led to an increase of energy consumption [9]. Water prices are determined by the EU Water Framework Directive [18], which provides a range of water-pricing tools that have been adopted by law in the EU Member States and have been applied in Southeast Europe on a voluntary basis [13]. Therefore, based on the recovery cost principle of water, which was defined by the Directive, an improvement of the water management is necessary to define the new sale prices for users, as well as the increase in the energy efficiency of these water systems. Decreasing the recovery cost principle and minimizing the water price that is paid by farmers in order to improve the efficiency are the main objectives of this study, as well as the resilience and sustainability of the water supply systems.
Moreover, it is worth mentioning that there is also a great cost variability of irrigation water depending on the water source and the location of the water source. In particular, the water scarcity and the low profit margins are the main problems for the farmers in the Mediterranean area. Therefore, these problems represent a high percentage of the production costs. In Spain, the sale price of water varies from 0.08 €/m 3 in some areas of Valencia (located in the mid-eastern part of Spain in the Mediterranean region) to 0.60 €/m 3 in Alicante (the southeast of Spain) (desalinated water) or Andalusian (the south of Spain) (groundwater source) [9]. It is common that the same irrigation community has different water sources, which usually have different extraction costs (e.g., an irrigation water system can receive water from groundwater that has a different piezometric head, water transferred between basins, or reused water that was treated in a wastewater treatment plant). For this reason, the decision about the most adequate water source from an economic viewpoint is critical. This selection significantly affects the final price that farmers have to pay [19], showing that there is a vital need for resilient and effective resource management strategies. This management has to be focused on societal, environmental, and economic aspects to promote sustainable development. Moreover, they also pointed to the need for decision-making mechanisms in order to address the question of resource utilization completely. To deal with this imperative need, different solution approaches have been developed.

Mathematical Programming Modelling Applied to Water Systems: Initial Overview
The aforementioned approaches were mainly focused on mathematical programming modelling and optimization aspects. To have more detailed knowledge about these solution approaches, a systematic literature review was carried out using the Web of Science database, covering all years and domain categories. The first approximation was focused on the definition of sets of general keywords to quantify the number of publications related to the mathematical approaches. Table 1 shows a summary of the results obtained. As an example, to illustrate the importance of mathematical approaches, the search of the keyword "business model" generated 10,465 results compared to the more than 200,000 results provided by the keyword "mathematical model". In light of this, it seems suitable to address the research on solution approaches related to mathematical programming modelling.   The first sets of keywords were then extended with the word "water". This combination of keywords provided a smaller number of results. As the search focused on the main topic of this research, the number of results decreased considerably. This tendency showed that mathematical solution approaches seemed to be universally accepted as reasonable mechanisms to solve miscellaneous problems. Therefore, they are also suitable for the management of the water sources. When the search was narrowed down further, adding "irrigation" to the set of originals keywords, and after the analysis of such results, it could be evinced that the proposal of solution approaches related to irrigation systems has been under-researched. Moreover, with the objective of delimiting the search further and excluding results that were not directly related to the main topic of the present research, only the following areas were chosen as the most connected with this study: water resources, agronomy, environmental sciences, agricultural engineering, agriculture multidisciplinary, engineering environmental, and green sustainable science technology. The results obtained in terms of the number of publications per year are also shown in Table 1, from 2009 to 2019. These findings were in concordance with the ones found by [20], who performed a review of mathematical programming applications in water resource management under uncertainty. Based on this review, it is worth highlighting the growing importance of developing pricing strategies to improve the water management in water irrigation networks.
In the literature, several and different mathematical programming methods were used to support water management in water irrigation networks under uncertainty. These methods, together with some developments and applications, are described below: • Multistage stochastic programming: Some studies use this method. For example, the work in [21] developed an interval multistage water allocation model to optimize water allocation between different growth stages to obtain the maximum food production in reservoir irrigation systems characterized by inputs' uncertainties. The study developed by [22] considered a fuzzy probability distribution based multistage stochastic robust programming method. This model supported regional water supply management. The developed model was applied to a water resources management system with three water users. • Stochastic dynamic programming: Among the references analyzed and in order to show their applications, the work in [23] used stochastic dynamic programming to model a farmer's choice whether to invest in a sprinkler irrigation system or in a more water efficient drip irrigation system under uncertainty. The work in [24] developed a stochastic dynamic programming model, but in this case, in the context of hydro-economic models to maximize irrigation benefits while minimizing the costs of power generation within a power market. The work in [25] also developed a stochastic dynamic programming model with fuzzy state variables for irrigation of multiple crops. This model, in which the reservoir storage and soil moisture of the crops are considered as fuzzy numbers and the reservoir inflow is considered as a stochastic variable, has the main objective of minimizing crop yield deficits, resulting in optimal water allocations to the crops by maintaining storage continuity and soil moisture balance.

•
Inexact programming including fuzzy and interval based programming: The work in [26] formulated a fuzzy mathematical programming model for a multi-reservoir system applied to a three reservoir system in the Upper Cauvery River basin, South India. The study tried to minimize the sum of deviations of the irrigation withdrawals from their target demands, on a monthly basis, over a year. Another study developed an interval-fuzzy two stage stochastic quadratic programming model. The goal was to allocate the limited irrigation water to different crops, maximizing the net benefit under uncertainty and to analyze how water allocation schemes change under different climate change scenarios [27].

•
Nonlinear programming: The work in [28] used a non-linear programming model to estimate farmers' willingness-to-pay for irrigation water that maximizes revenue from crop production under different shortage levels. In this case, Monte Carlo simulation was implemented considering model parameters' uncertainty to assess the variation of farmers' willingness-to-pay and avoid water shortage [28]. Other authors used nonlinear programing for the optimization of profitability and productivity in an irrigation command area with conjunctive water use options [29].

•
Multiobjective fuzzy linear programming: The work in [30] proposed a multiobjective fuzzy linear programming irrigation planning model for the evaluation of the management strategy in the case study of the Jayakwadi irrigation project, Maharashtra, India. Three conflicting objectives, net benefits, agricultural production, and labor employment, were considered in the irrigation planning scenario. However, the objectives pursued in this study are far from the main aim of the present research. In the same line, the work in [31] proposed a model of multiobjective fuzzy linear programming based on fuzzy parametric programming to solve the problem of optimal cropping pattern in an irrigation system. The objective of the irrigation planning model is to find out an optimal cropping pattern that maximizes simultaneously the net benefits, crop production, employment generation, and manure utilization.
In some cases, a combination of the previous methods has also been used to deal with aspects related to water irrigation networks under uncertainty.
Water management in water irrigation networks to satisfy the growing demand of water from different water sources is a subject that involves various uncertain impact factors. However, to our best knowledge, no evidence in the literature about the uncertainty in the different water source costs has been identified. For this reason, the present research is focused on developing strategies to improve the water management in water irrigation networks, studies related to water management resources in basins being outside the scope of this paper.
The different mathematical programming methods offer valuable information for water managers who will be able to choose more resiliently the water source depending on extraction costs. However, the aforementioned real cases were characterized through uncertainty scenarios in which statistical data were not reliable or even available. In these cases, the determination based models of the probability distributions for uncertain data may not be the best option, and therefore, stochastic approaches are not a feasible solution method [32,33]. In this context, fuzzy mathematical programming can prove to be an alternative approach, in which water managers can model the different types of uncertainties inherent to water source management processes. In this sense, the fuzzy approach can be employed to incorporate epistemic uncertainty or a lack of knowledge in the input parameters with analytical models or fuzziness in their objectives and also as a resolution technique of multi-objective programming models. Two main fuzzy mathematical applications can be highlighted. On the one hand, fuzzy mathematical programming can be employed to incorporate epistemic uncertainty or lack of knowledge in the input parameters to analytical models or fuzziness in their objectives [32][33][34][35][36][37]. On the other hand, fuzzy optimization can be used as a resolution technique of multi-objective programming models [36][37][38].
Based on all these findings, mathematical approaches are vital to make water management decisions more resilient and efficient. Moreover, due to the uncertainty that characterizes the operation of some irrigation systems, fuzzy techniques are fundamental to deal with these difficulties.

Research Goals
The previous section showed the needs to go in depth into the development of management techniques that improve the exploitation costs of water systems. In this line, the present research develops a mathematical programming model that considers the inherent uncertainty regarding operation costs. To do so, three fuzzy approaches were considered: (i) the first index of Yager, (ii) the third index of Yager, and (iii) Lai and Hwang. In the last case, the proposed approach generates multiple objectives considering the three different objectives. In order to transform this multi-objective model into an equivalent mixed integer linear programming model able to be solved by a commercial solver, several fuzzy multi-objective approaches were applied: Zimmerman, Werners, Selim, and Ozkarahan, and Torabi and Hassini.
The use of fuzzy methods is crucial to deal with the uncertainty of the results, and as previously mentioned, the use of fuzzy modelling is scarce with respect to the enhancement of irrigation systems.
The proposed model is able to make several and interrelated decisions with many variables and inaccurate or imprecise data in order to choose the most adequate origins of the water sources, minimizing the sale price and, therefore, improving the evaluation of water managers. The main objective of the present study is to develop a comparison, using different approaches, of the total acquisition cost of water sources. The comparison will be used to establish future research lines, which will support the decision-making process about the water sources as a function of different variables (e.g., quality, aquifer levels, and crops, among others).
This manuscript is structured as follows: The Introduction is given in this section. Section 2 defines the methodology of the formulation model, as well as the solution approaches. The third section presents a case study and the results obtained when the mathematical model is applied, to summarize finally the main conclusions in Section 4.

Methodology
The methodology of this research was focused on getting the minimum exploitation costs in a water irrigation system considering both maximum inlet flows, the maximum available volume of each water source, as well as the cost of each water resource, which can be fixed or variable. This methodology could be applied to any water pressurized systems and for their validation; it was applied in a real case study. The procedure has greater applicability and serves as a supporting mechanism for water irrigation managers to improve the sustainability of their systems.

Formulation Model
This section proposes the mathematical programming formulation for the management of the water sources problem [39]. This programming incorporates fuzzy costs in order to analyze the variation of source acquisition as a function of water costs. The hydraulic balance is proposed in Equation (2). The hydraulic balance is subject to different restrictions of the maximum flow and consumed volume. These constraints are defined in Equations (3)- (10). Moreover, the mathematical model considers the costs, which can be fixed or variable as they depend on the electric power price. The consideration of the costs is given in Equation (1), as well as the constrain (11). The sets of indexes, parameters, and decision variables in the fuzzy model are defined as follows: Indexes: i ∈ I Procurement water sources of the water network m ∈ M Procurement methods t ∈ T Time periods; in this case, 8760 periods were considered (one year) k ∈ K Months in the year Sets: Set of time periods in month k (720 periods for months that have 30 days) Constants of the model: Objective Function (1) aims to minimize the total costs related to the procurements of water from the different sources, including storage costs and also fixed and variable procurement costs. . This type of possibility distribution is represented by a triangular fuzzy number determined by the average or the most frequent value, the most optimistic, and the most pessimistic. For example, Figure 3 represents the possibility distribution for the fuzzy storage costs. This is subject to: Constraint (2) corresponds to the storage balance equation in the reservoirs.
Constraints (3) and (4) establish the safety stock and maximum storage capacity constraint, respectively.
Constraint (5) limits the required amounts from each source in each time period with respect to its maximum supply capacity.
Constraint (6) ensures that only one method is used for each source in each period.
Constraint (7) generates the activation of decision variable Yimt when decision variable Qimt is higher than 0.
Constraints (8) and (9) correspond to the limitation of the monthly volume for each source and the time limitation for each source and method, respectively.

Solution Approaches
In this section, different approaches for considering fuzzy costs in mathematical programming models are applied to the previously presented model. The considered approaches are the first index of Yager, the third index of Yager, and Lai and Hwang's approach.
A mathematical programming model with fuzzy costs may be formulated as follows: s. t. ≤ ∈ ∈ ∈ ℝ ∈ aij, bi ∈ M, j ∈ N, and ̃ ∈ ℝ where ℝ is the set of fuzzy numbers whose membership function represents the lack of precision of the values in objective function costs. The membership function is defined as follows: If triangular fuzzy numbers are considered in the form ̃ = , , , where rj, cj, and Rj correspond to the left (optimistic), center, and right (pessimistic) values, the membership function is given in the following way: According to Herrera and Verdegay [40,41], if the linear expression = ∑̃ is considered, in which ̃ are the fuzzy numbers with memebership functions similar to (13) and xj ≥ 0, then the membership function of can be formulated as follows: In order to obtain the optimal solution for (12), a ranking function f could be applied as Herrera and Verdegay [40,41] demonstrated. According to them, the first and third indices of Yager correspond to the following formulations.

First Index of Yager
The first function defuzzifies the uncertain costs by applying the concept of the center of gravity considering the following expression: Where the measure of the importance of the value z is expressed by g(z). If linear weights and triangular fuzzy numbers are considered, the problem (12) corresponds to the following equivalent formulation [40,41]: s. t. ≤ ∈ ℝ ∈

Third Index of Yager
The expression of the third index of Yager is as follows: where is the α-level set of and is the mean value of the elements of . Therefore, the optimal solution of (12) could be obtained by solving the following problem [40,41]: s. t. ≤ ∈ ℝ ∈

Lai and Hwang's Approach
Considering the problem defined by (12), each fuzzy cost is formulated as a triangular fuzzy number, which defines a fuzzy objective function with a triangular shape defined by the pessimistic, average or central, and optimistic values, according to the considered uncertain costs. In this sense, Lai and Hwang proposed to transform a mathematical programming model with fuzzy costs into a multiobjective mathematical programming model. This transformation was composed by three objective functions, which had the aim of minimizing the average costs (z1), as well as the difference between pessimistic and average costs (z3). Besides, simultaneously, the model maximized the difference between normal and optimistic costs (z2). Thus, minimization of the original fuzzy objective can be obtained by pushing these three critical points in the direction of the left-hand side. Therefore, the following auxiliary problem is obtained taking into account the previous nomenclature: In order to solve the obtained multiobjective mathematical programming model, several approaches can be applied, for example based on the use of fuzzy sets theory and goal programming. Readers are referred to the studies by Alavidoost et al. [42] and Torabi and Hassini [43] .

First Index of Yager
According to the first index of Yager [40,41] the following equivalent model is obtained: Objective function: and Constraints (2) to (11).

Third Index of Yager
An auxiliary crisp model is obtained by applying the third index of Yager [44,45] for ranking fuzzy numbers as follows: Objective function: and Constraints (2) to (11).

Lai and Hwang's Approach
The work in [46] introduced a solution method using uncertain costs, in which a fuzzy objective function is modelled with triangular possibility distributions into three different objective functions. According to [46], the original model for the management of water sources is converted into the following multiple objective model: This model minimizes the most possible value of the uncertain costs (14), whereas it maximizes the possibility of obtaining the lowest costs (15) and minimizes the possibility of obtaining the highest costs (16). The rest of the model is composed of Constraints (2) to (11).
In order to solve the previous multiple objective model, applying the fuzzy goal programming approaches of Zimmermann [47], Werners [48], Selim and Ozkarahan [49], and Torabi and Hassini [43], the membership functions (μ1, μ2, μ3) for each objective function (z1, z2, z3) are formulated according to [50] as follows: where , , and , , are the lower and upper bounds of each objective function. These values can be determined by the decision-maker according to his/her experience and personal criteria or by minimizing and maximizing separately each objective function as suggested by several authors such as [45] or [49]. Figure 4 shows the boundaries of the membership curves denoted by Equations (25)

The Zimmerman Solution Method
The previous multiple objective linear programming model can be transformed into an equivalent model with a single objective, according to the approach of [46], maximizing an auxiliary variable ϵ 0,1 , which represents the minimum degree of fulfilment of all the objectives: (28) subject to: and Constraints (2) to (11).

The Werners Solution Method
According to [47], the multiple objective model can be transformed into an equivalent single objective models as follows: subject to: and Constraints (2) to (11). Besides, λ0 has an equivalent meaning to Zimmermann's solution method. In contrast, γ corresponds to the compensation coefficient among the objectives.

Selim and Ozkarahan's Solution Method
Based on the solution method proposed by [47,48], a new fuzzy goal programming approach was proposed based on a weighted sum to incorporate the relative importance degree of each objective function in order to represent the decision-maker's preferences. The equivalent model obtained applying the Selim-Ozkarahan solution method is presented as follows: and Constraints (2) to (11) and (33) to (35).

Torabi and Hassini's Solution Method
The fuzzy goal programming approach proposed by [49] considers a convex combination of the minimum satisfaction degree and the weighted sum of the membership function values for each objective function. In this sense, the obtained single objective equivalent crisp model applying the Torabi-Hassini approach is: and Constraints (2) to (11) and (29) to (31). Among the solution approaches considered, the first and third indices of Yager are characterized by their simplicity of application since they perform the defuzzification of uncertain values through average operations. The model considers the central, the more optimistic, and the more pessimistic values of the triangular fuzzy numbers. However, the approach proposed by Lai and Hwang is based on the resolution of a multi-objective mathematical programming model. If this proposal is compared with Yager's indices, this presents a higher level of difficulty in order to obtain a solution by using additional proper solution methods, such as weighting objectives, the ξ-constraint, metaheuristics, or goal programming, among others. In this study, different fuzzy goal programming approaches have been considered. The obtained solution by applying Zimmerman's [47] solution method might be neither unique nor efficient because it is focused on the optimization of the objective with a lower degree of satisfaction according to the decision-maker. Similar to Zimmerman's approach, the Werners [48] solution method does not take into account the preferences of the decision-maker regarding their relative importance, although it introduces a compensation coefficient among them. In order to overcome the previous deficiencies, Selim and Ozkarahan [49] and Torabi and Hassini [43] presented solution methods in which the coefficient of compensation allowed accounting for the effect of the change in objective weights. According to [42], in these two approaches, a small change in the coefficient of compensation may yield a large amount of variability in the obtained solution (especially for lower values), and therefore, it is difficult to set this value in a precise way. Moreover, the formulations of Selim and Ozkarahan [49] and Torabi and Hassini [43] are more complex than those of Zimmerman or Werners and need more decision variables and computational time, but nevertheless, they can obtain better results.

Case Study
A real irrigation network was analyzed in order to show the innovation developed by this research. The water supply was used by [39] in order to analyze the necessary storage to supply the irrigation demand. The water system under study was located in Alicante province, supplying 260 hectares where there were vineyard crops. The water system was supplied from a small reservoir (its useful volume was 527,500 m 3 ). The storage is defined as a function of water level using Equation (38): S = 1259 WL 2 + 43,556 WL (38) where S is the storage in m 3 and WL is the water level of the reservoir in m. The consumption nodes (there were 110 irrigation consumption nodes) were located between 380 and 590 m above sea level. The pipelines that transport the water to the consumption nodes were manufactured from, iron and their range of diameters was from 80 to 550 mm.
The model was implemented with the MPL modelling language v5.0 [51], and the resolution was carried out with the optimization solver Gurobi v7.5.2 [52] in which a time limit of 180 s was set. Finally, it is worth mentioning that a Microsoft Access 2010 database managed the input and output data of the model.
The data used in this model were the same as those used by [39] to make a comparison between both mathematical programming methods applied on the same water system. However, in this case, the total costs were determined considering the cumulative storage cost to make possible the comparison among the different proposed methods. Figure 2 shows the water demand of the case study for each month, as well as the annual cumulative required water. This water demand was satisfied from different water sources. Source 1 came from water transfer of another basin, and its price was fixed over time. The water of Source 2 was obtained from wastewater treatment, and its price was also fixed over time. Finally, the water from Source 3 came from desalination. The prices of Sources 4 and 5 were variable since they came from wells, and therefore; their prices depended on the energy price.  Table 2 shows water prices for each of the five water sources. When the price was variable (Sources 4 and 5), the water price was also variable, and therefore, Table 2 shows the price for each period depending on the type of agreement ( Table 3). The possibility distributions associated with uncertain costs were modelled using triangular fuzzy numbers, which were constructed setting a maximum deviation of 15% for the most optimistic costs and 25% for the most pessimistic costs with respect to the most frequent ones. Table 2. Water prices as a function of the different sources according to [39]. With regard to the variable costs, they were related to the energy price. In Spain, there are different agreements that regulate this aspect. One of the agreements considers six periods from P1 to P6. Each period depends on the hour and date according to Table 3.  P2  P4  P5  P5  P4  P2  P2  P6  P4  P5  P4  P2  9  P2  P2  P4  P5  P5  P3  P2  P2  P6  P4  P5  P4  P2  10  P1  P1  P4  P5  P5  P3  P2  P2  P6  P4  P5  P4  P1  11  P1  P1  P4  P5  P5  P3  P1  P1  P6  P4  P5  P4  P1  12  P1  P1  P4  P5  P5  P3  P1  P1  P6  P4  P5  P4  P1  13  P2  P2  P4  P5  P5  P3  P1  P1  P6  P4  P5  P4  P2  14  P2  P2  P4  P5  P5  P3  P1  P1  P6  P4  P5  P4  P2  15  P2  P2  P4  P5  P5  P4  P1  P1  P6  P4  P5  P4  P2  16  P2  P2  P3  P5  P5  P4  P1  P1  P6  P3  P5  P3  P2  17  P2  P2  P3  P5  P5  P4  P1  P1  P6  P3  P5  P3  P2  18  P1  P1  P3  P5  P5  P4  P1  P1  P6  P3  P5  P3  P1  19  P1  P1  P3  P5  P5  P4  P2  P2  P6  P3  P5  P3  P1  20  P1  P1  P3  P5  P5  P4  P2  P2  P6  P3  P5  P3  P1  21  P2  P2  P3  P5  P5  P4  P2  P2  P6  P3  P5  P3  P2  22  P2  P2  P4  P5  P5  P4  P2  P2  P6  P4  P5  P4  P2  23  P2  P2  P4  P5  P5  P4  P2  P2  P6  P4  P5  P4  P2 Moreover, the maximum supply capacity of each source is shown in Table 4. The data are expressed in hourly and monthly capacity per source. Additionally, Sources 4 and 5 had a limited number of hours for supplying water depending on the month and the selected method (Table 5).    Table 6 summarizes the solutions yielded by the proposed models. From Table 6, it can be determined that the two models, based on the first and third indices of Yager, provided the worst results in terms of total costs compared to Lai and Hwang's approach. In this sense, Lai and Hwang's approach (1992) generated four different fuzzy solutions, depending on the solution method applied. Among them, Selim and Ozkarahan's approach obtained the best solution in terms of most possible (z1), most optimistic (z2), and most pessimistic (z3) total costs. Therefore, the decision-maker could develop a water management plan with a triangular possibility distribution for total storage and procurements costs (z1, z2, z3) = (278,450, 236,683, 348,063).  Table 6 also shows the computational time used by each model as the percentage of the one used by the first index of Yager to highlight the order of magnitude of additional computation. As can be seen, although there were differences among the different methods, these were not significant as the solutions tended to be rather fast.

Results
If Table 6 is analyzed in order to compare the results for each model used, considering the first index of Yager (FY) as the basis of the comparison, the total cost of the third index of Yager was 1% lower. The most possible total costs (z1) when Lai and Hwang's approach was used was 97.4% comparing it to FY. The different approaches varied between 278,450 (Selim and Ozkarahan's approach) and 278,579 € (Zimmerman's approach). However, the variation between the maximum and minimum costs was 0.05%. Therefore, the use of the different approaches was indifferent with respect to obtaining the optimal solution. When the most optimistic and pessimistic costs were analyzed (Table 6), a reduction of 17.8% was obtained, while the pessimistic costs showed an increase of 21.86 %. The annual volume used of each source as a function of the considered solution approach is shown in Table 7. Besides, Table 7 shows the percentage used of each considered water source for each solution approach.  Figure 3 shows an example of the water system management by using linear programming with Selim and Ozkarahan's approach. When all results over time were examined, Source 1 had a frequency of 19.7%. The frequency related to Source 2 was 1.13%, and the rest of the time, the water was obtained from Source 4 considering the energy period P6. Therefore, Source 4 was the main one used to minimize the exploitation cost of the water system. This trend was similar for the different statistics models developed in this research, where Source 4 was the main source to supply the water demand, considering P6. Regarding the operation of the irrigation system, Figure 6 shows the third week of June, for which the cumulative volume in the reservoir (i.e., storage) oscillated between 46,308 and 44,815 m 3 . This trend was stable over the planning horizon once the water system reached continuous time. Therefore, the planning demonstrated that it was not necessary to reach the maximum water level in the reservoir in the winter season since the water sources were enough to meet the annual water demand, minimizing the exploitation costs. Accordingly, the proposed model showed a significant improvement in comparison with the results given by the water manager, enhancing the current management [26]. In this sense, the proposed model allowed distributing the energy, storage, and management costs over the planning horizon simultaneously with incomes provided from selling water to farmers. Otherwise, the current planning method used by water manager concentrates most of the amounts of required water in winter months, and consequently, it is necessary to have large flows of cash available or, in some cases, to apply for loans from credit institutions. Table 8 presents the computational efficiency of the evaluated models. The results, from a computational viewpoint, were obtained by setting an upper bound of CPU time limited at 180 s and a stopping criterion for the gap of 0.5% so that the computational process stopped if either criterion was met. Moreover, data related to the number of constraints, decision variables, integer elements, and the nonzero and density of the coefficient matrix (defined as the percentage of nonzero coefficients in the total number of coefficients) are shown. The number of decision variables should be obtained by multiplying the size dimensions of each variable. However, in order to obtain a better computational performance, some filters related to impossible combinations of decision variables related to the case study were applied. As the evaluated models had a similar structure and the same input data, number of variables, integers, and constraints, the comparison presented very small or null variations. For example, the number of integer elements was the same for all the considered approaches, while the number of constraints and variables was slightly higher for Lai and Hwang's approach, and particularly so for Werners' and Selim and Ozkarahan's approaches. On the other hand, the number of nonzero elements was considerable higher for all the solution methods associated with Lai and Hwang's approach in comparison with Yager's indices.  Table 8 also shows how Lai and Hwang's approach required a considerably higher amount of CPU time compared to Yager's approaches. Moreover, the different approaches to solve Lai and Hwang's model also needed a higher number of iterations to obtain an optimal mixed integer linear programming solution, especially Zimmerman's solution method. Additionally, although Lai and Hwang's approach had a larger number of constraints and integer variables, this does not imply additional information storage requirements, but an increased modelling complexity. Finally, Werners' model and mainly Selim-Ozkarahan's and Torabi-Hassini's approaches required more parameter settings for the decision-maker.

Conclusions
Different solution approaches were applied and evaluated in terms of model complexity and computational efficiency, and the solution accomplished was presented. Lai and Hwang's approach, which was solved by Werners' method, implied more modelling complexity and CPU time, but provided more flexibility for decision-makers to achieve a fuzzy solution according to their preferences.
When the comparison was applied to a real water network, Lai and Hwang's approach reached the best solution in which the reduction of total costs was near 20% compared to the first index of Yager. The decrease of operational costs was very significant since this reduction was directly related to the reduction of energy consumption (e.g., pumped water or wastewater treatment, which established the water price). Therefore, a sustainability increase was reached when the decrease of costs was achieved.
Finally, the manuscript enables developing further lines of research. These lines can be focused on: (1) testing the models with real-world problems in a rolling horizon; (2) comparing the models with other solution approaches; (3) using other membership functions patterns; (4) including carbon footprint and the selection of sources of energy; (5) incorporating financial issues in order to optimize the generated cash flows and therefore improve the cost effectiveness.