An Optimization-Evaluation Agricultural Water Planning Approach Based on Interval Linear Fractional Bi-Level Programming and IAHP-TOPSIS

: In this study, an interval linear fractional bi-level programming (ILFBP) model was developed for managing irrigation-water resources under uncertainty. The ILFBP can fully address system fairness, uncertainties, and the leader–follower relationship of decision makers in the optimization process, which can better reflect the complexity of real decision-making process and help formulate reasonable water policies. An interactive fuzzy coordination algorithm based on satisfaction degree was introduced to solve the ILFBP model. In order to evaluate the applicability of optimization schemes, the interval analytic hierarchy process (IAHP) and the interval technique for order preference by similarity to an ideal solution (TOPSIS) method were integrated as IAHP-TOPSIS. To verify its validity, the developed optimization-evaluation framework was applied to an irrigation water management case study in the middle reaches of the Shiyang River Basin, located in the northwest China. The ILFBP model results show that the total water allocation is [6.73, 7.37] × 10 8 m 3 , saving nearly 0.9 × 10 8 m 3 more than the current situation. The benefit per unit of water is [2.38, 2.95] yuan / m 3 , nearly 0.4 yuan / m 3 more than the status quo, and the Gini coefficient is within a reasonable range of [0.35, 0.38]. The ILFBP model can well balance economic benefits and system fairness. Through the evaluation bases on IAHP-TOPSIS, the results of ILFBP show better water allocation effects and applicability than the other two models in this study area. Furthermore, due to various characteristics such as geographical location, population and area, there are three irrigation districts, Xiying, Donghe, and Qinghe, showing higher importance than others when considering regional water allocation. These findings can provide useful information for limited water resource managers and help decision makers determine effective alternatives of water resource planning under uncertainty.


Introduction
Water is an indispensable and essential resource for human survival and development.
With the development of social economy, the conflicts between supplies and demands of water resources have become increasingly prominent [1]. As the large freshwater consumer, agriculture accounts for more than 70% of the total freshwater consumption in the world [2,3]. Especially in the semi-arid and arid regions, water deficits limit sustainable development of the local agriculture, and thus more efforts should be made to improve the usage efficiency of limited irrigation water [4].
The agricultural water management systems are complicated with multiple uncertainties, including inexactness, fuzziness, and randomness existed in precipitation and river runoff, crop market prices variations, and irrigation target adjustments. The necessity of dealing with these uncertainties brings a lot of difficulties in the planning processes of irrigation water management [5,6]. To address these uncertainties in practical problems, many methods have been proposed for helping decision makers formulate better water policies [7][8][9][10][11]. Generally, these methods can be categorized into fuzzy programming (FP), stochastic programming (SP), and interval programming (IP). Among these categories, IP is the most widely used planning method under uncertainty due to its low data requirements [12,13]. Although IP can deal with the uncertainties expressed as interval numbers, it cannot fully reflect multiple decision-making hierarchy in the practical problems. To obtain the satisfied irrigation water allocation schemes, the various requirements from decision makers of different levels should be taken into account within one framework when formulating allocation models [14]. To reflect the tradeoff between two decision-making levels with different concerns, the bi-level programing (BP) approach with a two-level structure was developed for supporting regional water policies formulation [15]. A leader-follower decision-making strategy was incorporated into the optimization process of BP to generate satisfied decision-making plan for both the upper-level and the lower-level decision makers. However, more challenges exist when attempting to tackle interval uncertainty and coordinate the preferences of different-level stakeholders simultaneously.
In recent years, some researches were conducted to integrate interval linear programming (ILP) and BP into one framework to address such problems in irrigation water management [16]. An interval linear bi-level programming (ILBP) approach was thus developed and applied to real-world problems of water resources planning to address different-level decision-making preferences and uncertainties expressed as intervals [17]. These researches proved that ILBP would be an effective tool for decision makers to formulate water policies. However, the concerns of decision makers in these researches mainly focused on the economic and environmental factors such as benefits and water volumes water resource allocation systems. It was neglected that water resources are not only accounted as economic resources but also represent a kind of social resource. Unfair allocations among water users sometimes may lead to serious social conflicts. Therefore, the fairness issue needs to be considered in the water resource allocation process [18,19]. Many indicators have been proposed for measuring system fairness, such as the Gini coefficient, weighted proportional fairness, and so on [20][21][22]. Among them, the Gini coefficient is the most widely used measurement of inequality [23,24]. The Gini coefficient was proposed by the Italian economist Corrado Gini in 1912 to measure the balance of income in a country or region. It is usually calculated via the Lorenz curve [25], and has been applied to the fairness studies in many fields of education, economy, and resources management [22,[26][27][28]. Therefore, the Gini coefficient is considered as a proper configuration goal of water resources planning to guarantee the fairness of water allocation. However, because the Gini coefficient is always presented as fractional, it is necessary to introduce the interval linear fractional programming (ILFP) model into the ILBP framework to solve such a problem. This results in the development of an interval linear fractional bi-level programming (ILFBP) model, which can simultaneously tackle interval uncertainty, coordinate the preferences of different-level decision makers and reflect system fairness in one framework to obtain satisfied water allocation schemes.
Moreover, although the mathematical optimization models can generate water allocation schemes under uncertainty, few of them are incorporated with the results evaluation methods to measure their applicability and effectiveness in specific cases. This limitation may reduce the guiding significance of these optimization planning results in real-world practice [29][30][31]. Therefore, the evaluation methods based on actual regional conditions are necessary for measuring impacts of implementing the optimization results in various agricultural regions [32]. However, compared to general evaluation methods proposed for assessing the deterministic optimization results, the evaluation of optimization results under interval uncertainty would be more difficult. In order to address such a challenge, an interval analytic hierarchy process (IAHP) method has been proposed to address system uncertainties associated with systematic errors, incomplete information and fuzzy subjective judgment, and shown good performance in evaluating multiple schemes under interval uncertainty [33,34]. Generally, the calculation of IAHP often involves the ordering of interval numbers and the typical method is to directly use the mean of the interval numbers (i.e., the average of the upper and lower bounds) [35]. The limitation of this approach is that it ignores the effects of interval lengths and might lead to unreasonable comparison results. Besides, the possibility matrix has been used in many researches as a possible way to sort interval numbers [36,37], but it becomes complicated and inefficient when dealing with more than two indicators. More effective methods are thus needed to sort interval numbers in the results evaluation process. As the interval-parameter TOPSIS (ITOPSIS) could sort interval numbers via calculating the distances of alternative schemes to the ideal scheme, the integrated method of IAHP and ITOPSIS (IAHP-TOPSIS) is considered as an effective way to deal with interval numbers' sorting in the evaluation of optimization results [38][39][40]. Through the IAHP-TOPSIS method, the effects evaluation of the optimization schemes can be carried out and would further provide decision-making support for local water policies development. There have been few studies that can address both optimization and evaluation processes within one framework and simultaneously consider interval characteristics associated with agriculture water management systems.
The main objective of this study is to develop an optimization-evaluation approach for agricultural water planning based on the ILFBP and IAHP-TOPSIS methods to help decision makers formulate reasonable regional water policies. The optimization model (ILFBP) proposed for optimal irrigation water allocation and the evaluation model (IAHP-TOPSIS) used for evaluating the effectiveness of optimization schemes will be integrated within one framework under interval uncertainty to address the complexities of practical problems. Such an approach can: (1) Tackle interval numbers existed in both optimization and evaluation processes; (2) consider the fairness issue in water allocation; (3) make trade-offs between different concerns of two-level decision makers; (4) generate multiple water resource allocation schemes under different runoff levels; and (5) provide the corresponding evaluation results of the obtained allocation schemes. The proposed approach will also be applied to a real-world case study to demonstrate its applicability and effectiveness, wherein the regional water managers need to identify more effective irrigation water allocation policies to support the sustainable development of local agriculture.

Methodology
This section presents the development of an optimization-evaluation framework for an agricultural water allocation system, which entails the following three major components: (1) The parameter preprocessing and development of interval linear fractional bi-level programming (ILFBP); (2) the ILFBP model for irrigation water resources optimal allocation; and (3) the IAHP-TOPSIS method for evaluating the effectiveness of the optimization schemes. These components and their connections are shown in Figure 1 in detail.

Interval Linear Fractional Bi-Level Programming (ILFBP)
The general ILFBP problem can be described follows [41,42]: Subject to: where the upper-level variables are x 1 ∈ R n1 and the lower-level variables are x 2 ∈ R n2 ; c 11 , c 21 , c 12 , c 22 , d 11 , d 21 , d 12 , d 22 , α 1 , β 1 , α 2 , β 2 are interval parameters; G is the set of bi-level constraints; f 1 and f 2 are the upper-level and lower-level objective functions, respectively; the upper-level decision maker (ULDM) and the lower-level decision maker (LLDM) accept the leader-follower Stackelberg game [43].
To solve model (1), an interactive fuzzy coordination algorithm is proposed through analyzing the interrelationship between the objective functions and the constraints [16,41]. The upperand lower-bound submodels of the ILFBP can be decomposed [44]. In particular, when the number of interval parameters in the ILFBP model is relatively small, the upper-or lower-bound parameters corresponding to two ILFBP submodels can be determined respectively by comparing all the combinations of the intervals. If the objective functions of the ULDM and LLDM problems are formulated as the maximization goals, they need to be converted as the equivalent minimization goals, i.e., max f = min(− f ), and then the general ILFBP model can be solved according to following steps.
(1) The ILFBP submodel corresponding to the lower-bound objectives can be formulated as: Subject to: (2) The ULDM problem and the LLDM problem can be solved independently with the solutions of . The ranges of the decision variables (x − 1 ) should be around within the maximum tolerances (r − 1 ) of their most preferred values. The membership functions that specify the satisfaction degrees of x − 1 can be presented as follows: . The LLDM problem can be solved with the optimal solutions of (x L− 1 , x L− 2 ), which in turn provides the ULDM problem with the objective value of f − 1 by considering f − 1 ≥ f − 1 unacceptable in practice. The membership function that specifies the satisfaction degree of the ULDM objective value can be stated as: In contrast, the LLDM objective value has the following membership function: . (4) Consequently, the solution of submodel (2) with maximum overall satisfaction for both the ULDM and LLDM can be obtained through the following model: where λ − is the overall satisfactory degree of the lower-bound ILFBP submodel; I is a column vector with the same dimension as µ x − 1 or x − 1 and all elements of I are equal to 1. Thus, the solutions of the lower-bound ILFBP submodel and the corresponding objective values (i.e., x − 1opt , x − 2opt , f − 1opt , and f − 2opt ) can be obtained through model (6). (5) Accordingly, the upper-bound ILFBP submodel can be formulated as follows: (Upper-level objective) Subject to: Similar to the solution process of submodel (2), the solution of submodel (7) with maximum overall satisfaction for both the ULDM and LLDM can be obtained through the following model: Subject to: where: The solutions of the upper-bound ILFBP submodel and the objective values (i.e., x + 1opt , x + 2opt , f + 1opt , and f + 2opt ) can be obtained from model (8). (7) Finally, the solutions for the ILFBP model can be generated as follows:

The ILFBP Model for Agricultural Water Resources Management
To achieve effective management of limited irrigation water resources, the decision makers often desire optimal allocation schemes to avoid high costs and large water consumptions, especially for the arid and semi-arid areas. Moreover, the decision makers need to consider the system fairness, balance the different-level preferences, and address interval uncertainties in practical problems [45][46][47].
To deal with these challenges of agricultural water resources management, an ILFBP model can be formulated as follows: (Upper-level objective) where f 1 is the upper-level objective of maximizing the system's net economic benefits (yuan); i is the index for irrigation districts (i = 1, 2, . . .
where f 2 is the lower-level objective function to reflect the fairness of water allocation among irrigation districts via the Gini coefficient; and R i are populations of the irrigation district i. The Gini coefficient varies within the range of [0, 1]. The smaller Gini coefficient represents the higher fairness of water allocation among irrigation districts. Conversely, the larger Gini coefficient represents the lower fairness of water allocation. In regard to income distribution levels, a Gini coefficient less than 0.2 stands for absolute evenness, that between 0.2 and 0.3 represents relative evenness, that between 0.3 and 0.4 indicates basically reasonable distribution, and that between 0.4 and 0.5 indicates unreasonable distribution. When a Gini coefficient reaches 0.5 or more, it means a significant gap in income levels. Generally, a Gini coefficient higher than 0.4 is used as a sign of unfairness.
Subject to: (1) Water availability constraint where TW h are the total water supplies (m 3 ) under the flow level h; and TW ih denote the water supplies (m 3 ) for the irrigation district i under flow level h. (2) Irrigation water demand constraint where IR ± ijh,min and IR ± ijh,max are respectively the minimum and maximum crop water requirement (m 3 ) of crop j in the irrigation district i under flow level h.

The IAHP-TOPSIS Evaluation Approach Coupling the IAHP and ITOPSIS Methods
The IAHP method is widely used for multi-criteria decision making (MCDM) under uncertainty [48]. It can build an evaluation framework based on the analysis of interval parameters, and obtain the interval weights of the indicators. However, more effective methods are needed to sort interval numbers when using the IAHP. The interval-parameter TOPSIS method (ITOPSIS) can analyze the uncertain information of the scheme results and sort interval numbers through measuring the distances between solution alternatives and the ideal solution. Therefore, for improving the assessment accuracy under interval uncertainty, the IAHP and ITOPSIS methods are integrated to evaluate the optimal allocation results of water resources. Particularly, the IAHP method is used for calculating the indicators' interval weights and the ITOPSIS method is used for ranking the uncertain water resources allocation schemes. The specific steps can be shown as follows [33,40,49]: (1) Suppose a MCDM problem with m alternatives X m (m = 1, 2, · · · , M), and n criteria Y n (n = 1, 2, · · · , N). All indicators in MCDM are interval numbers, and thus an interval-number decision matrix A ± = (a ± mn ) M × N can be established. (2) Normalize the decision matrix to obtain a standardized interval-number decision matrix R ± = (r ± mn ) M × N as follows: where A n is the norm of the nth column vector and we have: )), n = 1, 2, · · · , N. Please note the reverse indicators need to be converted to positive indicators as: Therefore, all the metrics are positive metrics. (3) Invite experts to compare the importance of indicators. According to the experts' opinions, a judgment matrix can be built as follows with the interval numbers [31]: n , x + n , which corresponding to the largest eigenvalues of the judgment matrices B − and B + , respectively. Then calculate the interval weights using the interval number eigenvalue method as follows [50].
(5) Judgment matrix consistency test. If 0 ≤ k ≤ 1 ≤ l, the consistency of the interval judgment matrix B is reasonable [36]. If k > 1 or l < 1, the consistency of the judgment matrix is poor. The results should be sent back to the experts, and the judgment matrix should be corrected or re-judged until the consistency requirements are met. (6) Multiply the columns of normalized decision matrix by the associated weights: It can be calculated based on the calculation principle of interval number multiplication [49]. (7) Respectively determine the positive ideal and negative ideal solutions: The distance from the scheme alternative X m to the negative ideal solution is: . (10) Calculate the proximity of each scenario to a positive ideal.
(11) Rank the alternatives according to the relative closeness to the ideal alternatives, the larger value of C m means better effects of the corresponding scheme.

Study Area
Shiyang River Basin (101 , as shown in Figure 2, is one of three inland rivers in the Hexi Corridor of Gansu Province in northwest China. The Shiyang River Basin has a typical arid continental climate, which is featured by low and irregular rainfalls, high evaporation rates, and long drought periods [51][52][53]. Since the 1960s, the population of the Shiyang River Basin has increased significantly, the scale of agricultural production has expanded, and the underground water resources in the basin have been overexploited. These changes caused the deterioration of the ecological environment and water shortage problems. In order to restore the ecological environment of the Shiyang River Basin, the local government formulated a series of policies such as the "Shiyang River Basin Key Management Plan" [54] and a series of supportive policies. However, in recent years, due to the decreased river runoff, the water shortage problems in the middle reaches of the Shiyang River Basin are becoming more and more serious. Irrigation is the largest water consumer in the middle reaches of the Shiyang River Basin, accounting for 88.28% of the total water consumptions [29]. The low efficiency of irrigation water usage further deepens the gap between water demands and supplies. Therefore, the practical allocation methods of irrigation water resources are highly desired to improve water usage efficiency in the middle reaches of the Shiyang River Basin. The study area was a mountain plain located in the middle reaches of the Shiyang River Basin, which is called the Wuwei Basin. It includes two districts: The Liangzhou District of Wuwei City and the Yongchang County of Jinchang City. The dominated industry of the Wuwei Basin is agriculture and the local economy is relatively undeveloped. There are six rivers in the Wuwei Basin, namely the Dongda River, the Xiying River, the Jinta River, the Zamu River, the Huangyang River, and the Gulang River. These rivers can provide agricultural irrigation water for six irrigation districts located in the south of the study area, including Gulang (GL), Huangyang (HY), Zamu (ZM), Jinta (JT), Xiying (XY), and Donghe (DH). In addition, there are four groundwater irrigation districts in the north of the study area, including Qingyuan (QY), Jingyang (JY), Yongchang (YC), and Qinghe (QH) [55]. The surface water flows and groundwater withdrawals can be represented as stochastic distributions. Aiming to allocate the limited irrigation water resources in the Wuwei Basin, the decision makers from different hierarchical levels are responsible for making decision preferences from different viewpoints. The upper-level decision makers always rank the economy development as the primary goal, while the lower-level decision makers always focus on the fairness of water allocation among the irrigation districts to avoid social conflicts. Therefore, a comprehensive planning approach based on the ILFBP model and the IAHP-TOPSIS method would be applied to the study area to help water managers formulate well-rounded water allocation schemes under uncertainty. from different hierarchical levels are responsible for making decision preferences from different viewpoints. The upper-level decision makers always rank the economy development as the primary goal, while the lower-level decision makers always focus on the fairness of water allocation among the irrigation districts to avoid social conflicts. Therefore, a comprehensive planning approach based on the ILFBP model and the IAHP-TOPSIS method would be applied to the study area to help water managers formulate well-rounded water allocation schemes under uncertainty.

Data Collection and Processing
The inputs of the developed model include the estimated water demands for irrigation districts, the available supplies of surface water and groundwater, the relevant crop production data, the regional meteorological data, and the economic data. In detail, the available supplies of groundwater are based on the collected data from the Shiyang River Basin Key Governance Projects [56]. The available supplies of surface water are based on historical data of the six rivers' annual runoffs from 1955 to 2009. By fitting P-III curve based on the original runoff data (Figure 3) , m represents the number of years with annual runoffs not less than a certain flow value, and n represents the total number of years

Data Collection and Processing
The inputs of the developed model include the estimated water demands for irrigation districts, the available supplies of surface water and groundwater, the relevant crop production data, the regional meteorological data, and the economic data. In detail, the available supplies of groundwater are based on the collected data from the Shiyang River Basin Key Governance Projects [56]. The available supplies of surface water are based on historical data of the six rivers' annual runoffs from 1955 to 2009. By fitting P-III curve based on the original runoff data (Figure 3), the runoffs of six rivers can be classified as different flow levels as: Specifically, P < 12.5% corresponds to the very high flow level, 12.5% < P < 37.5% corresponds to the high flow level, 37.5% < P < 62.5% corresponds to the middle flow level, 62.5% < P < 87.5% corresponds to the low flow level, and P > 87.5% corresponds to the very low flow level (where P = m n + 1 × 100%, m represents the number of years with annual runoffs not less than a certain flow value, and n represents the total number of years with the flow data). Therefore, the probability levels (i.e., 0.125, 0.25, 0.25, 0.25, and 0.125) and the corresponding water amounts under the five flow levels can be obtained. In this study, a confidence interval with a significance level of 0.05 is used to obtain the interval ranges of water supplies. Table 1 shows the available surface water supplies of six rivers under different flow levels. with the flow data). Therefore, the probability levels (i.e., 0.125, 0.25, 0.25, 0.25, and 0.125) and the corresponding water amounts under the five flow levels can be obtained. In this study, a confidence interval with a significance level of 0.05 is used to obtain the interval ranges of water supplies. Table  1 shows the available surface water supplies of six rivers under different flow levels.    [55]. Table 2 shows the other parameters of the irrigation districts needed to solve the developed optimization model, such as the population data, irrigation water usage efficiencies, water productivities of different crops, areas of irrigation districts, and water prices. According to the water price adjustment policy of Wuwei City, the water price consists of two parts: The metering water price and the basic water price. For surface water supplies, the basic and metering water prices are respectively 30 yuan/ha and 0.157 yuan/m 3 ; for groundwater supplies, the basic and metering water price are respectively 120 yuan/ha and 0.01 yuan/m 3 .  The irrigation water demands can be obtained according to ET c = K c × ET 0 , where ET c is the actual crop evapotranspiration (mm/day), K c is the crop coefficient, and ET 0 is the reference crop evapotranspiration (mm/day). The ET 0 value can be calculated based on the FAO 56 Penman-Monteith method [57]. The meteorological data of the Wuwei Sub-basin from 1956 to 2017 are obtained from the China Meteorological Data Sharing Service System. Annual rainfall data from 1956 to 2017 are collected from the observation stations in Wuwei and the effective rainfall coefficient is 0.8 [58]. The effective rainfall is estimated by the product of actual rainfall and the effective precipitation coefficient. The main crops in the Wuwei Basin include wheat, corn, and economic crops, and their market prices come from the 2013-2016 China Agricultural Products Price Survey Yearbook.  [55]. Table 2 shows the other parameters of the irrigation districts needed to solve the developed optimization model, such as the population data, irrigation water usage efficiencies, water productivities of different crops, areas of irrigation districts, and water prices. According to the water price adjustment policy of Wuwei City, the water price consists of two parts: The metering water price and the basic water price. For surface water supplies, the basic and metering water prices are respectively 30 yuan/ha and 0.157 yuan/m 3 ; for groundwater supplies, the basic and metering water price are respectively 120 yuan/ha and 0.01 yuan/m 3 .  Figure 4 shows the optimization allocation results of each irrigation district under different flow levels. Obviously, with the decreased incoming water, the configuration target was reduced, and the water allocation in each irrigation district was also gradually reduced.   Figure 4 shows the optimization allocation results of each irrigation district under different flow levels. Obviously, with the decreased incoming water, the configuration target was reduced, and the water allocation in each irrigation district was also gradually reduced. Except for the irrigation districts of Xiying, Donghe, and Qinghe, the water allocation results for other irrigation districts could reach the allocation target. The Xiying irrigation district has the largest planting area, thus faces the most serious water shortage. The main water source for irrigation of the Xiying irrigation district is the Xiying River. The Xiying River is the largest tributary of the Shiyang River system, and the Xiying Reservoir is located at the exit of the mountain. The Xiying Reservoir not only provides irrigation water for the Xiying irrigation district, but also has the functions of preventing drought and flood, supplying safe drinking water for humans and animals, and transferring water to downstream. The current irrigation schemes cannot meet the water demands, resulting in serious waste of water resources and low irrigation water usage efficiency. These reasons basically answer the key question why water shortage often occurs in the Xiying irrigation district. Figure 5 shows the water allocation results for different crops in different irrigation districts under the medium flow level. The regional water allocation for wheat, maize, and the economic crops were respectively [5.25, 5.   Except for the irrigation districts of Xiying, Donghe, and Qinghe, the water allocation results for other irrigation districts could reach the allocation target. The Xiying irrigation district has the largest planting area, thus faces the most serious water shortage. The main water source for irrigation of the Xiying irrigation district is the Xiying River. The Xiying River is the largest tributary of the Shiyang River system, and the Xiying Reservoir is located at the exit of the mountain. The Xiying Reservoir not only provides irrigation water for the Xiying irrigation district, but also has the functions of preventing drought and flood, supplying safe drinking water for humans and animals, and transferring water to downstream. The current irrigation schemes cannot meet the water demands, resulting in serious waste of water resources and low irrigation water usage efficiency. These reasons basically answer the key question why water shortage often occurs in the Xiying irrigation district. Figure 5 shows the water allocation results for different crops in different irrigation districts under the medium flow level. The regional water allocation for wheat, maize, and the economic crops were respectively [5.25, 5.69] × 10 8 m 3 , [1.23, 1.34] × 10 8 m 3 , and [0.23, 0.25] × 10 8 m 3 . Aggressive decision makers may choose the upper bounds of water allocations to obtain higher food production and economic benefits, but they also need to bear relatively higher risks; conservative decision makers may choose the lower bounds of water allocations with relatively lower levels of benefits and risk. In general, the water allocation of food crops production is far greater than that economic crops, which has the similar trends with other previous studies in this study area [47,55]. For modern agricultural development, the water managers of irrigation districts should be able to adjust the regional planting structure and the irrigation modes in time. The water management department should guide local farmers to plant more economic crops to improve water usage efficiency. Additionally, the water supplies of economic crops should be increased to improve the agricultural economic benefits in irrigation districts.

The ILFBP Model Results of the Water Allocation
Water 2019, 11, x FOR PEER REVIEW 15 of 24 risk. In general, the water allocation of food crops production is far greater than that economic crops, which has the similar trends with other previous studies in this study area [47,55]. For modern agricultural development, the water managers of irrigation districts should be able to adjust the regional planting structure and the irrigation modes in time. The water management department should guide local farmers to plant more economic crops to improve water usage efficiency. Additionally, the water supplies of economic crops should be increased to improve the agricultural economic benefits in irrigation districts.  Figure 6 shows the comparison of the total regional water allocation among the current situation, the results of the ILFBP model, and the results of the economic model without considering the system fairness. The total water allocation of the ILFBP model in the study area was [6.73, 7.37] × 10 8 m 3 , while that of the economic model is [7.04, 7.97] × 10 8 m 3 and that of the status quo was [7.57, 8.37] × 10 8 m 3 . The total water allocation of the ILFBP model was the least one among the three schemes. Except for Xiying and Donghe irrigation districts, the water allocation of the economic model and ILFBP model was almost the same. The crop planting area and irrigation water allocation in the Xiying and Donghe irrigation districts were the largest among the study areas. However, the population of these two irrigation districts was not the largest. Therefore, in order to improve the water allocation fairness between irrigation districts, the ILFBP model appropriately reduced the water allocation of the Xiying and Donghe irrigation district.  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3   GL  HY  ZM  JT  XY  DH  QH Figure 6 shows the comparison of the total regional water allocation among the current situation, the results of the ILFBP model, and the results of the economic model without considering the system fairness. The total water allocation of the ILFBP model in the study area was [6.73, 7.37] × 10 8 m 3 , while that of the economic model is [7.04, 7.97] × 10 8 m 3 and that of the status quo was [7.57, 8.37] × 10 8 m 3 . The total water allocation of the ILFBP model was the least one among the three schemes. Except for Xiying and Donghe irrigation districts, the water allocation of the economic model and ILFBP model was almost the same. The crop planting area and irrigation water allocation in the Xiying and Donghe irrigation districts were the largest among the study areas. However, the population of these two irrigation districts was not the largest. Therefore, in order to improve the water allocation fairness between irrigation districts, the ILFBP model appropriately reduced the water allocation of the Xiying and Donghe irrigation district. risk. In general, the water allocation of food crops production is far greater than that economic crops, which has the similar trends with other previous studies in this study area [47,55]. For modern agricultural development, the water managers of irrigation districts should be able to adjust the regional planting structure and the irrigation modes in time. The water management department should guide local farmers to plant more economic crops to improve water usage efficiency. Additionally, the water supplies of economic crops should be increased to improve the agricultural economic benefits in irrigation districts.  Figure 6 shows the comparison of the total regional water allocation among the current situation, the results of the ILFBP model, and the results of the economic model without considering the system fairness. The total water allocation of the ILFBP model in the study area was [6.73, 7.37] × 10 8 m 3 , while that of the economic model is [7.04, 7.97] × 10 8 m 3 and that of the status quo was [7.57, 8.37] × 10 8 m 3 . The total water allocation of the ILFBP model was the least one among the three schemes. Except for Xiying and Donghe irrigation districts, the water allocation of the economic model and ILFBP model was almost the same. The crop planting area and irrigation water allocation in the Xiying and Donghe irrigation districts were the largest among the study areas. However, the population of these two irrigation districts was not the largest. Therefore, in order to improve the water allocation fairness between irrigation districts, the ILFBP model appropriately reduced the water allocation of the Xiying and Donghe irrigation district.  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3  m=1  m=2  m=3   GL  HY  ZM  JT  XY  DH  QH Figure 7 shows the economic benefits of the three water allocation schemes. The economic benefits of the economic model, ILFBP model, and the status quo were respectively [16.49, 23.06] × 10 8 yuan, [16.01, 21.73] × 10 8 yuan, and [15.66, 22.04] × 10 8 yuan. The economic model only maximized the regional overall economic benefits brought by the water allocation, so more water was allocated to the irrigation areas with more economic crops to improve the overall economic benefits. Therefore, the economic model could bring higher economic benefits. The benefits of the actual situation and the ILFBP model were similar, but the interval ranges of the system benefits were different. The upper-bound benefits of the status quo were larger than those of the ILFBP model, indicating that the current situation could obtain greater benefits. However, the larger interval range of the benefits under the status quo was larger indicates greater level of uncertainties existed in the benefits of the status quo. It seems that more irrigation water could bring more economic benefits, but the marginal benefits existed in the relationship between economic benefits and water allocation amounts, which may lead to the decreased water usage efficiency when using more irrigation water. Therefore, the water allocation schemes could not be evaluated only through comparing the amount of irrigation water and economic benefits. An appropriate evaluation method is still needed to compare the effects of water allocation schemes and provide more useful decision-making information.
Water 2019, 11, x FOR PEER REVIEW 16 of 24 Figure 6. The total water allocation of the three water allocation schemes in each irrigation district. Note: m = 1, m = 2, and m = 3 represent the economic model, ILFBP model, and status quo, respectively. Figure 7 shows the economic benefits of the three water allocation schemes. The economic benefits of the economic model, ILFBP model, and the status quo were respectively [16.49, 23.06] × 10 8 yuan, [16.01, 21.73] × 10 8 yuan, and [15.66, 22.04] × 10 8 yuan. The economic model only maximized the regional overall economic benefits brought by the water allocation, so more water was allocated to the irrigation areas with more economic crops to improve the overall economic benefits. Therefore, the economic model could bring higher economic benefits. The benefits of the actual situation and the ILFBP model were similar, but the interval ranges of the system benefits were different. The upper-bound benefits of the status quo were larger than those of the ILFBP model, indicating that the current situation could obtain greater benefits. However, the larger interval range of the benefits under the status quo was larger indicates greater level of uncertainties existed in the benefits of the status quo. It seems that more irrigation water could bring more economic benefits, but the marginal benefits existed in the relationship between economic benefits and water allocation amounts, which may lead to the decreased water usage efficiency when using more irrigation water. Therefore, the water allocation schemes could not be evaluated only through comparing the amount of irrigation water and economic benefits. An appropriate evaluation method is still needed to compare the effects of water allocation schemes and provide more useful decision-making information.

Evaluation Results
In order to reflect the economic, social, and resource effects of water-allocation schemes, the evaluation model indicators such as the Gini coefficient, the irrigation water productivity, the benefits per unit of water usage, the irrigation loss per unit of area, and the water shortage rate were calculated based on the results obtained from three optimization models. The values of the indicators of the three schemes are shown in Table 3.

Evaluation Results
In order to reflect the economic, social, and resource effects of water-allocation schemes, the evaluation model indicators such as the Gini coefficient, the irrigation water productivity, the benefits per unit of water usage, the irrigation loss per unit of area, and the water shortage rate were calculated based on the results obtained from three optimization models. The values of the indicators of the three schemes are shown in Table 3. Y 1 , · · · Y 5 represent the above-mentioned five indicators, respectively, and the weights of the five indicators are represented by the weight vector w. The IAHP was used for comparing the estimation indicators to build an interval judgment matrix and calculate the results of the index weights as shown in Table 4. The calculation results of k = 0.92 < 1 and l = 1.06 > 1 could pass the consistency test. The weighted normalization decision matrix is shown in Table 5.
The relative closeness of each scheme is C 1 = 5.36 × 10 −1 , C 2 = 5.49 × 10 −1 , C 3 = 5.12 × 10 −1 . From the above results, the ranking result of the three schemes can be obtained as The evaluation results show that the applicability of the ILFBP model is better than other two models in this study area. The regional economic, resource, and social indicators can be improved by adjusting the water allocation schemes for different irrigation districts. Among the three schemes, the state quo showed the highest level of fairness among the three schemes with the Gini coefficient of only [0.29, 0.31], but the regional economic growth and local irrigation water productivity were also limited by strong fairness concerns. In the economic model, in order to obtain greater economic benefits, more water was distributed to the irrigation districts with more economic crops, such as Xiying, Donghe, and Qinghe. However, this leads to unfair water distribution between irrigation districts. The Gini coefficient of the economic model was [0.39, 0.41], which was close to triggering the warning sign. The ILFBP model coordinates the decision-making preferences of two-level stakeholders with the considerations of system economic benefits and system fairness. These efforts contribute to the social contradiction caused by the extremely unfair water allocation among irrigation districts. The Gini coefficient of the ILFBP scheme was lower than 0.4. The superiority of the ILFBP model could also be reflected through the indicator comparison, such as irrigation loss, irrigation water productivity, and benefit per unit water. These results indicate that the water allocation scheme obtained by ILFBP could effectively reduce unnecessary water waste and make most of the indicators closer to the positive ideal solution.

Discussion
This study used the ILFBP model as a planning management tool to improve irrigation-water use efficiency and further increase available ecological water in the Shiyang River Basin. The results obtained from ILFBP model provided specific adjustment and allocation recommendations under different flow levels, which contributed to the sustainable development of local agriculture.
The effects evaluation of the water allocation schemes could test the applicability of the ILFBP model to the specific regions, and further improve the implementation significance of the optimization schemes for practical problems. The IAHP-TOPSIS evaluation method could deal with uncertainties and provide more practical strategies than existing deterministic evaluation models. The uncertainties associated with the water allocation variables and the subjective preferences could be fully reflected in the expert evaluation process.
In order to evaluate the effects of the ILFBP model comprehensively, two models, respectively the status quo situation and the single economic goal without the fairness concern, were built and solved. The obtained results showed that the ILFBP model could coordinate the economic benefits and fairness of system under uncertainty. Moreover, it could improve the water allocation effects from multiple perspectives, and obtained an irrigation-water allocation scheme with a high satisfaction degree for different stakeholders. The water allocation schemes obtained from the ILFBP model were obviously more reasonable than the status quo, and could meet the policy making requirements of the Shiyang River Basin. These efforts could help local water managers of the Shiyang River Basin formulate more suitable water policies.
Compared with other previous studies on water resources allocation in the Shiyang River Basin [29,47,51], the ILFBP model could obtain a satisfied irrigation-water allocation scheme by considering the system fairness and making tradeoff between the upper-and lower-level decision makers under uncertainty. Moreover, the optimization-evaluation framework developed in this study could not only provide specific water allocation schemes, but also evaluate the effectiveness and the applicability of the optimization schemes, making the results more convincing and more realistic.
Furthermore, the optimization-evaluation framework developed in this study could be applied to other regions for improving water resource management at a government, regional, local, and basin level. The specific planning benefits and the evaluation results are associated with irrigation water allocation, ecological water allocation, different flow levels, and interval range of social and economic parameters. These data can be obtained through field research or model programming. Therefore, the proposed optimization-evaluation framework under uncertainty is universal and has significant implications for supporting regional water allocation practices, especially in arid and semi-arid regions.

Conclusions
In this study, an interval linear fractional bi-level programming (ILFBP) model was developed for irrigation-water allocation planning under uncertainty. The bi-level planning (BP), interval planning (IP) and linear fractional programming (LFP) techniques are integrated to deal with the problems in reality. The interactive fuzzy coordination algorithm was used for solving the developed model, and the interval solutions under different flow levels could be obtained. In order to reasonably evaluate the effects of the optimization schemes under uncertainty, the interval analytic hierarchy process (IAHP) method and the interval TOPSIS method were simultaneously applied for evaluating the obtained allocation schemes. The developed optimization-evaluation approach was applied to a real-world case study in the middle reaches of the Shiyang River Basin. Compared with the previous studies, the developed framework had the following advantages. (1) It could consider the decision-making preferences of multi-level stakeholders and reflect different priorities existed in the decision-making process. The goals of the upper-level decision makers could be prioritized, and then be weighed among the two levels of decision makers, and finally the decision-making plan with higher overall satisfaction could be obtained; (2) the fairness goal could be added to the water resources optimization model to ensure system fairness of allocation schemes under uncertainty; (3) it could quantify the subjective intention of the experts and reflect the uncertainty in this process; and (4) it could provide the corresponding effects evaluation for optimization schemes under interval uncertainty. The proposed approach could help decision makers develop and identify more appropriate water-irrigation schemes. This framework could also be applied to other practical problems of resources and environmental management for developing well-rounded strategies under various system complexities.

Parameters The ILFBP Methods Parameters and Variables
Meaning and Description U Upper-level L Lower-level opt The optimal values x 1 , x 2 The upper-and lower-level decision variables f 1 , f 2 The upper-and lower-level objective functions c 11 , c 21 , c 12 , c 22 , d 11 , d 21 , d 12 , d 22 Interval parameters G The set of bi-level constraints r 1 The maximum tolerances of the upper-level decision variables f 1 , f 2 The worst values of the upper-and lower-level objective functions µ x 1 (x 1 ) The satisfaction degrees of the upper-level decision variables The satisfaction degrees of the upper-and lower-level objective functions λ The overall satisfactory degree of the ILFBP model The prices of crop j in the irrigation district i (yuan/kg) Y ij The water productivities of crop j in the irrigation district i (kg/m 3 )

The application of the ILFBP model Parameters and Variables
Meaning and Description A i Irrigated areas of subarea i (ha) k i The irrigation water use efficiency P ijh Effective precipitation (m 3 ) JB i Basic water price (yuan/ha) JL i The metering water price (yuan/m 3 ) R i Populations of the irrigation district i TW h The total water supplies under the flow level h (m 3 )

TW ih
Water supplies for the irrigation district i under flow level h (m 3 ) IR ± ijh,min , IR ± i jh,max The minimum and maximum crop water requirement of crop j in the irrigation district i under flow level h (m 3 )

The IAHP-TOPSIS evaluation methods Parameters and Variables
Meaning and Description X m Alternatives (m = 1, 2, · · · , M) Y n Indicators (n = 1, 2, · · · , N) A ± = (a ± mn ) M × N Interval-number decision matrix R ± = (r ± mn ) M × N Standardized interval-number decision matrix A n The norm of the jth column vector The interval weights x ± n Eigenvectors k, l Parameters required to calculate interval weights S ± = (s ± mn ) M × N Weighted normalized decision matrix The distance from the scheme alternative X m to the positive and the Negative ideal solution, respectively X + , X − The positive ideal and negative ideal solutions C m The relative closeness