Analysis of the Effects of Water Management Strategies and Climate Change on the Environmental and Agricultural Sustainability of Urmia Lake Basin , Iran

In arid and semi-arid areas, unsustainable development of irrigated agriculture has reduced the water level of large lakes such as Aral, Urmia, Hamoon, and Bakhtegan. Urmia Lake, as a hyper saline and very shallow lake, located in the northwest of Iran, has water level reductions of about 40 cm each year over the past two decades. In this research, the indices of environmental and agricultural sustainability are evaluated using performance criteria influenced by climate change and water management strategies for the Zarrinehrud and Siminehrud River basins as the largest sub-basin of Urmia Lake basin. Modeling of hydrologic behavior of these basins is performed using WEAP21 model. The model is analyzed for three future emission scenarios (A2, A1B, and B1), for the period of 2015–2040 and five water management scenarios: (1) keeping the existing situation; (2) crop pattern change; (3) improving the conveyance and distribution efficiency; (4) combining the improvement of conveyance and distribution efficiency with improving the application efficiency using modern technology; and (5) the combination of crop pattern change with the improvement of total irrigation efficiency. The results show that the highest values of indices of environmental sustainability and agricultural sustainability are related to the scenario of combining the crop pattern change with improving the total irrigation efficiency under the B1 emission scenario (B1S4).


Introduction
In basins where water consumption is mostly for irrigation, sustainable management of water resources must consider two goals simultaneously: sustainable irrigated agriculture to ensure food security and protection of the environment.It is necessary to establish sustainable interaction between these two goals in the present and future, and at the same time decrease potential conflicts between these two goals with the help of methods such as the application of new methods of irrigation, preventing water loss in the conveyance paths, changing the crop pattern to less water use products, developing greenhouse cultivation, etc. [1].Increasing water demand often results in unsustainable water consumption and inadequate water remain for environmental protection.To save water for the environment, measures for saving agricultural water have been introduced as an agenda for environmental policy in many arid and semi-arid regions [2].Törnqvist and Jarsjö [3] investigated the hydrological effects of improving irrigation technologies on vast cotton fields in the Aral Sea basin located in Central Asia.The results of their research showed that the implementation of irrigation technologies will lead to water storage and will increase the discharge to the Aral Sea between 1 and 6 km 3 per year.
In arid and semi-arid areas, high-consumption agriculture has led to a decrease in the water level of large lakes such as Aral Sea, Urmia Lake, Hamoon Lake, and Bakhtegan Lake [4].Urmia Lake, the largest lake in Iran and the second most saline lake in the world, has been dying during the last two decades and this has resulted in similar or even greater socio-environmental impacts to the disaster of the Aral Sea [5].
Various studies have been conducted on Urmia Lake regarding the effective factors in reducing the inflow to the lake and they are referred to in the following: Jalili et al. [6] investigated the effects of climate change and human factors on fluctuations in water level of Urmia Lake.The results showed that the recent decrease in the water level of Urmia Lake is mostly due to human impacts rather than climate change.The results of some research studies [7][8][9][10] also confirmed that human factors have a greater impact on drying Urmia Lake than climatic factors.Shadkam et al. [11] examined the effects of climate change and the development of water resources (including constructing dams and increased irrigated agriculture) on the reduction of inflow to the Urmia Lake during the period of 1960 to 2010.The results showed that the annual inflow to Lake Urmia decreased by 48% during the studied period.About 60% of this decline was caused by climate change and about 40% by water resources development.Urmia Lake Basin is located in an area with the following problems: poverty, water stress, soil degradation, natural agro-constraints, and isolation [12].
Generally, hydrological models focus on understanding how water flows within a basin in response to hydrological events, while water resource planning models focus primarily on water allocation management (for example, deciding on supply and demand).The integration of these two components is possible using a single analytical framework.This framework is the Water Evaluation and Planning Model of version 21 (WEAP21), which is a bridge between the basin hydrology and water management by combining the physical hydrological process in the water management framework [13].
In connection with the application of the WEAP model in the analysis of water saving scenarios and the effects of climate change, several studies have been conducted and some of them are referred to in the following.
Joyce et al. [14] reviewed the agricultural water management strategies for adaptation to climate change using the WEAP model.Adaptation strategies include improving irrigation technologies and crop pattern change into more valuable or less water consuming products.Mehta et al. [15], using the WEAP model, examined the supply and demand of irrigation water, due to climate change and land use scenarios in the Cache Creek watershed in California.They tested three scenarios until 2099: (1) crop pattern change based on economic forecasts; (2) change to more varied cultivation patterns and more efficient water use; and (3) combining improvement of irrigation technology and crop pattern change.Their results showed that demand for irrigation water increases by 26% and 32%, respectively, under the B1 and A2 climatic scenarios at the end of the century.The greatest water saving is achieved in the scenario of combining varied cultivation patterns and more efficient water use, with improvements in irrigation technology, which reduces demand by about 12% below the mean time.Blanco-Gutiérrez [16] conducted a hydrologic-economic analysis of water management strategies for balancing water for nature (protecting wetlands) and water for food on the Guadiana River Basin in Spain.The results of their study showed that, in periods of long-term droughts (predicted increase in temperature between 1.5 and 3.6 • C and 10% to 20% reduction of precipitation by 2050), tensions will be intensified on agriculture, land resources and water resources by decreasing water availability and increasing demand for irrigation water.
Considering the lack of water resources and uncertainty in the climate and future demand, decision-makers need sustainability analysis techniques in water resource management [17].For the first time, Loucks developed an index for quantifying the sustainability of water resource systems using three criteria: reliability, resilience and vulnerability.He defined the sustainability index for assessing the effectiveness of various managerial options from the perspective of stakeholders of water and the environment.Sustainability index estimates system capacity to reduce its vulnerability.If a proposed policy leads to a more stable system, this index indicates that the system has a high capacity to reduce vulnerability in the future.Subsequently, the proposed index was used in numerous studies by researchers [18].Criteria of reliability, resilience and vulnerability measure various aspects of water resource system performance.These three criteria together provide one of the most comprehensive approaches for analyzing the probability of success or failure of a system, the rate of recovery (or return) of a system from undesirable conditions, and the quantification of expected results in satisfactory modes for long periods [19].Safavi et al. [20] analyzed the sustainability index of water resources and consumption in Zayandehrud River basin under the scenarios of integrated water resources management and planning.They defined the water resources sustainability index based on the performance criteria of time and volume reliability, resilience, vulnerability and maximum deficit.The results showed that reducing drinking, industry, and agriculture water, respectively, to 10%, 15%, and 30%, in line with the demand management scenario, and the implementation of Koohrang 3 and Behesht Abad plans in line with the scenario of supply management could improve or start improving the sustainability of water resources from weak mode to moderate to good.Yilmaz and Harmancioglu [21] used the WEAP model to evaluate water resource management in the Gediz River basin in Turkey.The initial input of the proposed model was the amount of surface water that is largely devoted to irrigation purposes.Therefore, the mutual relations of supply and demand in agricultural water consumption were the main focus of their study.The results of their study based on various indices such as sustainability index showed that Gediz River basin is completely sensitive to drought conditions and agricultural sector is significantly affected by irrigation deficit, especially in drought periods.
The main objective of this research was to evaluate the effects of future agricultural water management strategies under the influence of climate change on the inflow to the Urmia Lake and the agricultural condition during the years 2015-2040 using the WEAP21 model [22,23].Considering that the Zarrinehrud and Siminehrud River basins are the largest sub-basins in the Urmia Lake basin, this research was carried out on these two sub-basins.The evaluation of the effects of these strategies was done by indices of environmental sustainability, agricultural sustainability, and irrigation water deficit.

Study Area
The Urmia Lake basin with an area of 51,876 km 2 is located in the mountainous area of northwest of Iran and among West Azarbaijan provinces with 42% of the basin area, East Azarbaijan with 37% of the basin area and Kurdestan with 11% of the basin area.Urmia Lake with an area of about 5100 km 2 covers 10% of the basin area and is shared between the provinces of West Azerbaijan and East Azerbaijan [24].Considering that Zarrinehrud and Siminehrud River basins provide about 52% of environmental flow requirements of Urmia Lake [25] and are among the main sub-basins of the Urmia Lake basin, this research is performed on these two sub-basins.Basin of Zarrinehrud-Siminehrud is located in the southern and southeastern parts of Lake Urmia and in terms of the size it is the largest sub-basin of Lake Urmia basin (Figure 1).The rivers of Zarrinehrud and Siminehrud with a maximum discharge of about 3 billion cubic meters per year in the period 1995-2014 are considered as water rich rivers in the country, which together with numerous natural and social potentials, have led to the development of agricultural activities in the watersheds of the basin.The length of Zarrinehrud River is about 240 km and the area of the basin is 11,840 km 2 .The length of Siminehrud River is roughly 200 km and the area of the basin is 3785 km 2 .The construction of irrigation networks such as Zarrinehrud (Miandoab), etc., and dams such as Shahid Kazemi Boukan dam with a reservoir volume of 810 MCM have created facilities at the basin level, as a result of which about 136,000 hectares of the lands of the basin are cultivated by irrigation either traditionally or by modern methods.

Introducing the WEAP Model
WEAP is a modeling platform that can provide integrated assessment of climate, hydrology, land use, irrigation facilities, water allocation and water management priorities of the watershed.The WEAP model uses a standard linear programming model to solve water allocation problems at any time step and its target function is to maximize the percentage of supplying demand centers' needs, with regard to supply and demand priority, mass balance and other constraints.All constraints are defined intermittently for each step of time, with regard to the priority of supply and demand.The WEAP model calculates the mass balance equilibrium of water for each node and branch at either daily or monthly time steps [22,23].Using the time series of the climate, the WEAP model calculates the hydrological cycle components by simulating the rainfall-runoff process at the surface of the catchment area [26].In this study, the soil moisture method of WEAP model was used for modeling hydrological reaction of basins and inter-basins [27].

Soil Moisture Method
This one-dimensional, two-bucket soil moisture method is based on empirical functions that describe evapotranspiration, surface runoff, interflow, baseflow, and deep percolation for a basin unit (Figure 2).For a basin subdivided into several subbasins with different fractional land use or soil- The rivers of Zarrinehrud and Siminehrud with a maximum discharge of about 3 billion cubic meters per year in the period 1995-2014 are considered as water rich rivers in the country, which together with numerous natural and social potentials, have led to the development of agricultural activities in the watersheds of the basin.The length of Zarrinehrud River is about 240 km and the area of the basin is 11,840 km 2 .The length of Siminehrud River is roughly 200 km and the area of the basin is 3785 km 2 .The construction of irrigation networks such as Zarrinehrud (Miandoab), etc., and dams such as Shahid Kazemi Boukan dam with a reservoir volume of 810 MCM have created facilities at the basin level, as a result of which about 136,000 hectares of the lands of the basin are cultivated by irrigation either traditionally or by modern methods.

Introducing the WEAP Model
WEAP is a modeling platform that can provide integrated assessment of climate, hydrology, land use, irrigation facilities, water allocation and water management priorities of the watershed.The WEAP model uses a standard linear programming model to solve water allocation problems at any time step and its target function is to maximize the percentage of supplying demand centers' needs, with regard to supply and demand priority, mass balance and other constraints.All constraints are defined intermittently for each step of time, with regard to the priority of supply and demand.The WEAP model calculates the mass balance equilibrium of water for each node and branch at either daily or monthly time steps [22,23].Using the time series of the climate, the WEAP model calculates the hydrological cycle components by simulating the rainfall-runoff process at the surface of the catchment area [26].In this study, the soil moisture method of WEAP model was used for modeling hydrological reaction of basins and inter-basins [27].

Soil Moisture Method
This one-dimensional, two-bucket soil moisture method is based on empirical functions that describe evapotranspiration, surface runoff, interflow, baseflow, and deep percolation for a basin unit (Figure 2).For a basin subdivided into several subbasins with different fractional land use or soil-type areas, the mathematical formulation to compute the storage change in the first layers is expressed in terms of a water balance as follows [28,29]: where z 1,j is the relative soil water storage, a fraction of the total effective water storage in the root zone layer in area j (dimensionless); Rd j is the soil water holding capacity of the area j (mm); Pe is the effective precipitation (mm); PET(t) is the reference potential evapotranspiration (mm/day, potential evapotranspiration is derived using a Penman-Monteith formulation); k c,j is the crop coefficient for area j; and RRF j is the Runoff Resistance Factor of the land cover.Higher values of RRF j lead to less surface runoff.P e (t)z is the surface runoff; f j k s,j z 2 1,j is the interflow from the first soil layer for area j; f j is the partition coefficient related to the land cover type, soil, and topography for area j, which divides flow into horizontal f j and vertical (1 − f j ) flows; and k s,j is the saturated hydraulic conductivity of the root zone layer for area j [mm/time].
Water 2018, 10, x FOR PEER REVIEW 5 of 25 type areas, the mathematical formulation to compute the storage change in the first layers is expressed in terms of a water balance as follows [28,29]: where z1,j is the relative soil water storage, a fraction of the total effective water storage in the root zone layer in area j (dimensionless); Rdj is the soil water holding capacity of the area j (mm); Pe is the effective precipitation (mm); PET(t) is the reference potential evapotranspiration (mm/day, potential evapotranspiration is derived using a Penman-Monteith formulation); kc,j is the crop coefficient for area j; and RRFj is the Runoff Resistance Factor of the land cover.Higher values of RRFj lead to less surface runoff.The change of storage in the second layer (dz2/dt) is computed by: where Smax is the deep percolation from the upper layer storage; and ks2 is the saturated hydraulic conductivity of the lower storage (mm/time).

Model Setup
The model was developed for Zarrinehrud and Siminehrud River basins using WEAP21 software.The schematic of the model is shown in Figure 3.The change of storage in the second layer (dz 2 /dt) is computed by: where S max is the deep percolation from the upper layer storage; and k s2 is the saturated hydraulic conductivity of the lower storage (mm/time).

Model Setup
The model was developed for Zarrinehrud and Siminehrud River basins using WEAP21 software.The schematic of the model is shown in Figure 3. Two points of the environmental flow requirements of Lake Urmia (from Zarrinehrud (ZFR) and Siminehrud (SFR) Rivers).

Sources of Data and Information
In this research, a 20-year period (1995-2014) was considered for calibration and validation of the model.Meteorological stations inside and near the basins studied were used for six synoptic stations of Zarrineh Obato, Tekab, Saghez, Mahabad, Miandoab, and Maragheh for meteorology parameters of precipitation, temperature, relative humidity, wind speed and cloudiness fraction.The mean of each of the meteorological parameters for each sub-basin was estimated using the Thiessen Polygons method.In the WEAP21 model, the basins of Zarrinehrud and Siminehrud were divided into upstream and downstream sub-basins according to the hydrometric gauge stations used.Zarrinehrud basin includes Tekab, Saghez, Shahin Dezh, and right bank and left bank of Zarrinehrud irrigation network and Siminehrud basin includes Boukan 1 and Boukan 2 agriculture catchments (Figure 3).Two points of the environmental flow requirements of Lake Urmia (from Zarrinehrud (ZFR) and Siminehrud (SFR) Rivers).

Sources of Data and Information
In this research, a 20-year period (1995-2014) was considered for calibration and validation of the model.Meteorological stations inside and near the basins studied were used for six synoptic stations of Zarrineh Obato, Tekab, Saghez, Mahabad, Miandoab, and Maragheh for meteorology parameters of precipitation, temperature, relative humidity, wind speed and cloudiness fraction.The mean of each of the meteorological parameters for each sub-basin was estimated using the Thiessen Polygons method.In the WEAP21 model, the basins of Zarrinehrud and Siminehrud were divided into upstream and downstream sub-basins according to the hydrometric gauge stations used.Zarrinehrud basin includes Tekab, Saghez, Shahin Dezh, and right bank and left bank of Zarrinehrud irrigation network and Siminehrud basin includes Boukan 1 and Boukan 2 agriculture catchments (Figure 3).

Calibration and Validation of the WEAP Model
Calibration and validation of a model is necessary to minimize errors and to ensure the validity and reliability of the model results [17].In this research, the model was calibrated by comparing the discharge data simulated in a monthly scale with observation data at selected hydrometric gauge stations.An 11-year period (from water year 1994-1995 to 2004-2005) was used to calibrate the model and a period of nine years (from water year 2005-2006 to 2013-2014) was used to validate the model.
In this study, to simulate flow discharge data, the parameters of non-agricultural crop coefficient, soil water capacity, root zone conductivity and preferred flow direction were calibrated using PEST (Parameter Estimation) [28].The values of non-calibrated parameters were selected based on spatial and temporal conditions and the default values of the model.During this process, the flow simulated by the model was compared to the flow observed in the selected hydrometric gauge stations.
For quantitative evaluation of the results of calibration and validation of the model, Pearson's correlation coefficient, Nash-Sutcliffe Efficiency (NSE) and Index of Agreement (d) were used.The relations of the Nash-Sutcliffe Efficiency and Index of Agreement are as follows [31,32].
where Q obs i is the observed discharge for the ith period, Q mean obs is the mean of the observed discharge for the entire period, Q sim i is the simulated discharge for the ith period, and n is the total number of data.

Climate Change
To investigate the effect of climate change on water resources of the basins, the HadCM3 climatic model was used under the emission scenarios A2, A1B and B1 based on the IPCC (Intergovernmental Panel on Climate Change, Geneva, Switzerland) [33] for the period 2015-2040.SRES A2, A1B and B1 have a similar trajectory to RCP8.5, RCP6.0 and RCP4.5, respectively [34].The output of climatologic models does not have the spatial and temporal analysis accuracy required for hydrologic models, therefore, it is necessary to downscale the output of climatic models for the region in question.In this research, the LARS-WG model was used to downscale the data.The inputs of this model are minimum and maximum daily temperature, daily precipitation and daily sunshine hours [35].The baseline period in this study was 1986-2010 for six synoptic stations of Zarrineh Obato, Tekab, Saghez, Mahabad, Miandoab, and Maragheh.The model was evaluated by comparing baseline data and data generated by the model using RMSE, MSE, MAE, R 2 , Bias and NSE.After evaluating the performance of the LARS-WG model using the relevant indices and ensuring the suitability of the model, data for the period 2015-2040 were generated.After data generation for the upcoming period from rainfall, minimum temperature and maximum temperature, the monthly means of these parameters were calculated for the upcoming period (2015-2040).Then, their changes were calculated relative to the mean values of the baseline period and the temperature and precipitation diagrams were plotted for each of the stations.
Because LARS-WG is a stochastic weather generator, it was run 100 times for each emission scenario.These 100 realizations per emission scenario were used in the WEAP model for computing water fluxes which were averaged afterwards.Although Guo et al. [36] concluded that 25 realizations are sufficient and there is no significant difference between 25 and 100 realizations, we followed the recommendation of the LARS-WG manual.

Evaluation of the Impact of Climate Change to the Hydrological Component Using the WEAP Model
To analyze the future impact of climate change on streamflow and ET, the WEAP model is used for the hydrological processes using the A2, A1B and B1 climate projection from the HadCM3 and LARS-WG.
To quantify the change in hydrology of the basin as an impact of climate change, the comparison of streamflow and ET in baseline period and future emission scenarios is performed [37].To determine the baseline conditions, the validated WEAP model was used to generate the streamflow and ET in the Zarrinehrud and Siminehrud River basins during the period 1986-2010.Similar steps were applied to generate the streamflow and ET during future period (2015-2040) using downscaled temperature and precipitation corresponding to A2, A1B and B1 climate projections.

Developing Scenarios
In this research, scenarios have been developed based on management strategies and future climate change.The main goal of the management strategies is to increase water supply, increase irrigation efficiency or reduce agricultural water demand, so that the saved water can be allocated to Urmia Lake and increase its sustainability.To develop future climate scenarios, outputs of rainfall and temperature downscaled by the LARS-WG model for the HadCM3 general circulation model (GCM) for the three emission scenarios of A2, A1B and B1 for the period 2015-2040 were used.
Specifications of agriculture catchments for the current account ( 2015) is described in Table 1.Selected scenarios are combination of management strategies and climate change as follows: (1) Base scenario (S0) The base scenario is based on a future climate prediction (A2, A1B, and B1).The cropped area and future crop pattern are considered constant for the base scenario and are equal to 2014 data.However, domestic water demand varied according to the predicted population growth rate and will be applied to all other scenarios.
Water 2018, 10, 160 S1 is a combination of climate change and a change in the present crop pattern with crops which have lower water requirements.In this scenario, high-consumption crops of alfalfa and sugar beets will be replaced by low-consumption wheat, barley and rapeseed.In each agriculture catchment, a small amount of alfalfa cultivation will not change to provide local requirement.
(3) Scenario 2 (S2) S2 combines climate change and improvement of conveyance and distribution efficiency by 15% of present status to reduce water losses.In this scenario, improvements in conveyance and distribution efficiency will be achieved through rehabilitation of the Zarrinehrud irrigation network, dredging and canal lining, water conveyance through pipelines, and education and technical support for farmers.
(4) Scenario 3 (S3) S3 is a combination of S2 and improving application efficiency by changing the systems of surface irrigation to sprinkler and drip irrigation.In this scenario, the application efficiency of sprinkler irrigation system is considered as 75% and the application efficiency of drip irrigation system is considered as 95%.
(5) Scenario 4 (S4) S4 is a combination of S1 and S3.In this scenario, in addition to changing the crop pattern to lower-consumption crops, the conveyance, distribution, and application efficiency (total irrigation efficiency) also increased.

Evaluation Indicators
To assess the considered scenarios, Environmental Sustainability Index (ESI), Agricultural Sustainability Index (ASI) and Irrigation Water Deficit (IWD) were considered (Table 2).

Indicator Description
Environmental Sustainability Index (ESI) The temporal aggregation of supply/demand ratio time series (only for environmental flow requirements) based on performance criteria when satisfactory value is equal to 1 (full coverage) Agricultural Sustainability Index (ASI) The temporal aggregation of supply/demand ratio time series (for irrigation only) based on performance criteria when satisfactory value is considered between 0.8 and 1 Irrigation Water Deficit (IWD) Indicates annual unmet demand for irrigation (annual mean value is used in evaluations); 10 6 m 3 ESI is used for assessing the status of providing the environmental flow requirements (EFRs) of Urmia Lake from the basins of Zarrinehrud and Siminehrud and indicators of agricultural sustainability and irrigation water deficit are used to evaluate the irrigation status of the Zarrinehrud and Siminehrud basins.

Performance Criteria
The sustainability index (SI) is calculated by integrating the reliability, resilience and vulnerability performance criteria.

Reliability (Rel)
Reliability is defined as the probability of any specific C t value, whose C t is time series of any selected indicator C, within the range of satisfactory values as follows: Number of satisfactory values Total number of simulation period (5) Resilience (Res) Resilience is a criterion for describing the rate of recovery from an unsatisfactory condition.It is the probability of occurrence of a satisfactory C t+1 value following an unsatisfactory C t value and is defined as follows: Number of times a satisfactory value occurs after an unsatisfactory value Total number unsatisfactory values ( 6) Vulnerability is a statistical measure of size (magnitude) or duration of failures in a given time series.The size (magnitude) of a failure is a value where a C t exceeds the upper limit of satisfactory values (UL(C t )), or below the lower limit of satisfactory values (LL(C t )).In this research, vulnerability is defined as the expected extent-vulnerability, in which the duration of failures is eliminated and defined as follows: Total individual extents of failure Total number of individual extents of failures (7) In the above equations, a satisfactory value is the "supply/demand ratio" that for environmental flow requirements is equal to 1 (full coverage) and for irrigation is considered between 0.8 and 1.Other values are considered as unsatisfactory values.

Sustainability Index (SI)
Sustainability is a relatively new concept for measuring the performance of water resource systems over long periods [43].This index is proposed to facilitate the assessment of a water management policy and to compare it to other policies [44].The sustainability index changes from 0 for the lowest and the worst possible value up to 1 for the highest and best possible value [21].In this research, the sustainability index is calculated as Equation ( 8).This index assures the properties of the SI defined by Loucks [39], but also has the following improvements: Content, Scaling, and Flexibility [45].ESI and ASI are calculated based on Equation (8) with Equations ( 9) and (10).
In Equation ( 8), Rel, Res and Vul are Reliability, Resilience and Vulnerability, respectively.
The index used for ESI is supply/demand ratio (Se/De) of Urmia Lake, in which demand is Urmia Lake's environmental flow requirements (EFRs) and the index used for ASI is supply/demand ratio (Sa/Da) of agriculture catchments for reliability, resilience and vulnerability performance criteria.
For a group comparison of water users, sustainability based on the group (SG) was defined as the weighted average of sustainability indices.SG is used to calculate the sustainability for a group k with water users from i to j belonging to this group [45].
Water 2018, 10, 160 In Equation ( 11), W i is relative weight for the ith water user, which varies between 0 and 1 and the sum is equal to 1.In this study, the water users are equal to the agriculture catchments (i) and the k group is Zarrinehrud and Siminehrud basins.For Zarrinehrud basin, j = 5 (number of agriculture catchments) and for Siminehrud basin, j = 2 (number of agriculture catchments).
If SI for each water user is weighted by its annual water demand, SG for the kth group is defined as follows: in which: Water demand i (13)

Environmental Flow Requirements (EFRs) of Urmia Lake
Abbaspour and Nazaridoost [46] determined the environmental flow requirements of Urmia Lake using an ecological approach.The three variables of ecology, water quality and water quantity were considered as environmental indicators to calculate the lake environmental flow requirements using the ecological approach.The salinity threshold of 240 ppt is considered as the water quality requirement for Artemia Urmiana.The aggregation of these variables showed that Urmia Lake requires about 3100 MCM of water flow per year for a sustainable ecosystem, which is equivalent to the ecological level of the lake water level (1274.1 m.a.s.l.).
In addition, based on long-term historical data, the Zarrinehrud basin provides about 41% and the Siminehrud basin provides about 11% of the Urmia Lake's environmental flow requirement [25].Therefore, the share of supplying the Urmia Lake's environmental flow requirements from the Zarrinehrud and Siminehrud basins is estimated as 1271 and 341 MCM per year, respectively.

WEAP Model of the Study Area
The monthly simulated discharge by the WEAP model is in good agreement with the observed discharge in the hydrometric gauge stations of the entrance of Boukan Dam, Nezam Abad, Dashband and Miandoab Bridge (Figure 4).The values of the correlation coefficient indices, Nash-Sutcliffe Efficiency (NSE) and Index of Agreement (d) are also shown in Figure 4.If SI for each water user is weighted by its annual water demand, SG for the kth group is defined as follows: in which:

Environmental Flow Requirements (EFRs) of Urmia Lake
Abbaspour and Nazaridoost [46] determined the environmental flow requirements of Urmia Lake using an ecological approach.The three variables of ecology, water quality and water quantity were considered as environmental indicators to calculate the lake environmental flow requirements using the ecological approach.The salinity threshold of 240 ppt is considered as the water quality requirement for Artemia Urmiana.The aggregation of these variables showed that Urmia Lake requires about 3100 MCM of water flow per year for a sustainable ecosystem, which is equivalent to the ecological level of the lake water level (1274.1 m.a.s.l.).
In addition, based on long-term historical data, the Zarrinehrud basin provides about 41% and the Siminehrud basin provides about 11% of the Urmia Lake's environmental flow requirement [25].Therefore, the share of supplying the Urmia Lake's environmental flow requirements from the Zarrinehrud and Siminehrud basins is estimated as 1271 and 341 MCM per year, respectively.

WEAP Model of the Study Area
The monthly simulated discharge by the WEAP model is in good agreement with the observed discharge in the hydrometric gauge stations of the entrance of Boukan Dam, Nezam Abad, Dashband and Miandoab Bridge (Figure 4).The values of the correlation coefficient indices, Nash-Sutcliffe Efficiency (NSE) and Index of Agreement (d) are also shown in Figure 4.
In general, it can be said that the agreement between observed and simulated values both in the calibration period and in the validation period is good with respect to the indices and using the WEAP21 model is recommended for water resource planners to simulate and evaluate the different scenarios.

Precipitation and Temperature
The results of the comparison of the average monthly precipitation and temperature observations and those simulated using the LARS-WG model based on the relevant indices for the Zarrineh Obato station are shown in Table 3.For Tekab, Saghez, Mahabad, Miandoab and Maragheh stations, the indices of RMSE, MSE, MAE, R 2 , BIAS and NSE for Precipitation were between 2.53 and 4.73, 6.38 and 22.41, 2.07 and 3.66, 0.96 and 0.99, −2.08 and 0.38, and 0.94 and 0.99, respectively; for Tmin, they were between 0.16 and 0.26, 0.03 and 0.07, 0.13 and 0.22, 0.99 and −0.06, and 0.04 and 0.99, respectively; and, for Tmax, they were between 0.24 and 0.38, 0.06 and 0.14, 0.21 and 0.27, 0.99 and 0.04, and 0.1 and 0.99, respectively.The results show that the model is well able to make possible changes in the future according to the future climate scenarios.In general, it can be said that the agreement between observed and simulated values both in the calibration period and in the validation period is good with respect to the indices and using the WEAP21 model is recommended for water resource planners to simulate and evaluate the different scenarios.

Precipitation and Temperature
The results of the comparison of the average monthly precipitation and temperature observations and those simulated using the LARS-WG model based on the relevant indices for the Zarrineh Obato station are shown in Table 3.For Tekab, Saghez, Mahabad, Miandoab and Maragheh stations, the indices of RMSE, MSE, MAE, R 2 , BIAS and NSE for Precipitation were between 2.53 and 4.73, 6.38 and 22.41, 2.07 and 3.66, 0.96 and 0.99, −2.08 and 0.38, and 0.94 and 0.99, respectively; for Tmin, they were between 0.16 and 0.26, 0.03 and 0.07, 0.13 and 0.22, 0.99 and −0.06, and 0.04 and 0.99, respectively; and, for Tmax, they were between 0.24 and 0.38, 0.06 and 0.14, 0.21 and 0.27, 0.99 and 0.04, and 0.1 and 0.99, respectively.The results show that the model is well able to make possible changes in the future according to the future climate scenarios.

Future Climate
The average annual precipitation, and minimum and maximum temperature for the studied stations for the base period (1986-2010) and the future period (2015-2040) are shown in Figures 5 and 6.The precipitation of the studied stations in the future period is increasing for all emissions scenarios as compared to the base period.Most precipitation of emission scenarios for all stations is related to the B1 emission scenario and the lowest amount to the A2 emission scenario.In addition, the minimum and maximum temperature of the studied stations in the future period will increase for all emission scenarios as compared to the base period.

Future Climate
The average annual precipitation, and minimum and maximum temperature for the studied stations for the base period (1986-2010) and the future period (2015-2040) are shown in Figures 5 and  6.The precipitation of the studied stations in the future period is increasing for all emissions scenarios as compared to the base period.Most precipitation of emission scenarios for all stations is related to the B1 emission scenario and the lowest amount to the A2 emission scenario.In addition, the minimum and maximum temperature of the studied stations in the future period will increase for all emission scenarios as compared to the base period.

Future Climate
The average annual precipitation, and minimum and maximum temperature for the studied stations for the base period (1986-2010) and the future period (2015-2040) are shown in Figures 5 and  6.The precipitation of the studied stations in the future period is increasing for all emissions scenarios as compared to the base period.Most precipitation of emission scenarios for all stations is related to the B1 emission scenario and the lowest amount to the A2 emission scenario.In addition, the minimum and maximum temperature of the studied stations in the future period will increase for all emission scenarios as compared to the base period.

Change in Streamflow
Figure 7 shows the projected changes of the average annual streamflow in the future period (2015-2040) under A2, A1B and B1 emission scenarios compared with the baseline period (1986-2010) for Zarrinehrud and Siminehrud River basins.As shown in Figure 7, streamflow is projected to decrease in the future for all emission scenarios in comparison with the baseline period.
Water 2018, 10, x FOR PEER REVIEW 14 of 25

Change in Streamflow
Figure 7 shows the projected changes of the average annual streamflow in the future period (2015-2040) under A2, A1B and B1 emission scenarios compared with the baseline period (1986-2010) for Zarrinehrud and Siminehrud River basins.As shown in Figure 7, streamflow is projected to decrease in the future for all emission scenarios in comparison with the baseline period.

Change in Actual Evapotranspiration (ETc)
The result of the ETc analysis shows that the average annual ETc under the all emission scenarios is projected to increase compared to the baseline period (Figure 8).In this study, precipitation and temperature are projected to increase in the future.Potential ET increase due to increase in temperature which results in an increase in actual ET (ETc).As the ETc increases in the future, the streamflow decreases since the increase in ETc will be more than the increase in precipitation.

Change in Actual Evapotranspiration (ETc)
The result of the ETc analysis shows that the average annual ETc the all emission scenarios is projected to increase compared to the baseline period (Figure 8).
Water 2018, 10, x FOR PEER REVIEW 14 of 25

Change in Streamflow
Figure 7 shows the projected changes of the average annual streamflow in the future period (2015-2040) under A2, A1B and B1 emission scenarios compared with the baseline period (1986-2010) for Zarrinehrud and Siminehrud River basins.As shown in Figure 7, streamflow is projected to decrease in the future for all emission scenarios in comparison with the baseline period.

Change in Actual Evapotranspiration (ETc)
The result of the ETc analysis shows that the average annual ETc under the all emission scenarios is projected to increase compared to the baseline period (Figure 8).In this study, precipitation and temperature are projected to increase in the future.Potential ET increase due to increase in temperature which results in an increase in actual ET (ETc).As the ETc increases in the future, the streamflow decreases since the increase in ETc will be more than the increase in precipitation.In this study, precipitation and temperature are projected to increase in the future.Potential ET increase due to increase in temperature which results in an increase in actual ET (ETc).As the ETc increases in the future, the streamflow decreases since the increase in ETc will be more than the increase in precipitation.4. The highest increase for reliability and resilience and highest decrease for vulnerability for all agriculture catchments is related to S4, or the crop pattern change with increasing irrigation efficiency.According to Table 4, the implementation of water management scenarios has increased reliability and resilience and reduced vulnerability.Performance criteria of reliability, resilience and vulnerability for Zarrinehrud and Siminehrud River basins in connection with the supply of Urmia Lake's environmental flow requirements are presented in Table 5.The highest increase for reliability and resilience and highest decrease for vulnerability is again related to S4.According to Table 5, the implementation of water management scenarios has increased reliability and resilience and reduced vulnerability.For Zarrinehrud and Siminehrud River basins, there is a combination of three indicators (ASI, IWD and ESI) against five scenarios.Three combinations are obtained for the three emission scenarios of A2, A1B and B1, because these combinations change with climatic conditions.Table 6 presents sustainability indices of agriculture (according to Table 4) and irrigation water deficit for agriculture catchments.ASI varies between 0 and 1. ASI values close to 1 indicate good performance.As shown in Table 6, in S0, ASI values for all agriculture catchments are less than 1 for A2, A1B, and B1 emission scenarios.This means that all agriculture catchments are under water stress and significant improvements for ASI can be achieved under different policies.The ASI value for all other considered scenarios has increased, but the highest amount of ASI for all agriculture catchments is related to B1S4, or the combination of crop pattern change with increasing total irrigation efficiency under B1 emission scenario.In addition, because the agriculture catchments of Saghez, Tekab and Shahin Dezh have higher potential to improve sustainability due to lower irrigation efficiency, in these areas, the increase in ASI is higher.
According to Table 6, water management strategies (S1-S4) lead to reduced water demand (by changing the crop pattern) and increased water supply (by increasing the irrigation efficiency), thereby mitigating the effects of climate change.
Irrigation water deficit (IWD) or unmet demand is a valuable index for evaluation.As shown in Table 6, in S0, IWD values for all agriculture catchments are significant.This deficit in many areas can reduce yield and reduce farmer's income.However, by the proposed solutions, IWD can be greatly improved.According to Table 6, the IWD value in B1S4 shows the highest improvement for all agriculture catchments.
Water 2018, 10,160 In Figure 9, the average amount of water is released from the Zarrinehrud and Siminehrud River basins to Urmia Lake in the future (2015-2040) under various emission scenarios is shown.The share of the environmental flow requirements (EFRs) of Lake Urmia from the Zarrinehrud and Siminehrud River basins is estimated as 1271 and 341 MCM per year, respectively.As shown in Figure 9, the maximum average amount of water entering the lake from these basins is related to the scenario of mixing the crop pattern change with the improvement of total irrigation efficiency under the B1 emission scenario (B1S4).Under this scenario, the average annual flow of water entering to Urmia Lake from the basins of Zarrinehrud and Siminehrud is 1292 and 351 MCM per year, and, as a result, the lake's environmental flow requirements will be supplied from these basins.
greatly improved.According to Table 6, the IWD value in B1S4 shows the highest improvement for all agriculture catchments.
In Figure 9, the average amount of water is released from the Zarrinehrud and Siminehrud River basins to Urmia Lake in the future (2015-2040) under various emission scenarios is shown.The share of the environmental flow requirements (EFRs) of Lake Urmia from the Zarrinehrud and Siminehrud River basins is estimated as 1271 and 341 MCM per year, respectively.As shown in Figure 9, the maximum average amount of water entering the lake from these basins is related to the scenario of mixing the crop pattern change with the improvement of total irrigation efficiency under the B1 emission scenario (B1S4).Under this scenario, the average annual flow of water entering to Urmia Lake from the basins of Zarrinehrud and Siminehrud is 1292 and 351 MCM per year, and, as a result, the lake's environmental flow requirements will be supplied from these basins.Table 7 shows the indices of environmental sustainability (according to Table 5), agricultural sustainability and irrigation water deficit (according to Table 6) for the Zerrinehrud and Siminehrud River basins.ASI is achieved based on weighted average of agriculture catchments (group sustainability) for each basin.The irrigation water deficit index is also based on the total deficit of agriculture catchments for each basin.According to Table 7, the highest ESI and ASI for both basins are related to B1S4.In addition, the greatest improvement in IWD is related to B1S4.In B1S4, both agriculture and the environment (Urmia Lake's environmental flow requirements) are stable.Notes: 1 Environmental Sustainability Index (ESI). 2 Agricultural Sustainability Index (ASI). 3Irrigation Water Deficit (IWD).
The aforementioned results show that, by adopting suitable management strategies, it is possible to improve agricultural sustainability and irrigation deficit, save water consumption, use saved water for environmental use (Urmia Lake) and improve its sustainability.The results of this research are in line with the results of other researchers [1,3,15,21].The researchers concluded that improving irrigation efficiency would increase agricultural sustainability and save water consumption, and the saved water can be used for other intentions, including the environment.

Conclusions
In this research, environmental and agricultural sustainability indices are evaluated using performance criteria influenced by water management strategies and climate change.The irrigation water deficit index for agriculture catchments is also analyzed.The main results are summarized as follows: 1.
The current unsustainable conditions in the Urmia Lake basin are due to climate change, water stress, poverty, agricultural development, cultivation of high-consumption crops, soil degradation, and low irrigation efficiency.

2.
The climate scenarios of this study show an increase in the average temperature, precipitation, and actual evapotranspiration, and a decrease in streamflow for the future period (2015-2040) compared to the baseline period (1986-2010) for the Zarrinehrud and Siminehrud River basins.

3.
In this study, the sustainability index (SI) was successfully used to assess water management strategies and climate change.

4.
The results of this research showed that water management strategies including crop pattern change with increasing irrigation efficiency is necessary for achieving environmental and agricultural sustainability.By applying these strategies, the performance criteria of reliability and resilience are increased and the vulnerability decreases, thus the sustainability index increases.

5.
The highest amount of environmental sustainability index (ESI) is related to the scenario of combining the crop pattern change with improving the total irrigation efficiency under the B1 emission scenario (B1S4).In this scenario, the average annual flow of water entering to Urmia Lake from the basins of Zarrinehrud and Siminehrud is 1292 and 351 MCM per year, respectively, and, as a result, the lake's environmental flow requirements (EFRs) will be supplied from these basins.6.
The highest amount of agricultural sustainability index (ASI) is related to the scenario of combining the crop pattern change with improving the total irrigation efficiency under the B1 emission scenario (B1S4).In this scenario, the average amount of irrigation water deficit (IWD) for Zarrinehrud and Siminehrud river basins were 33.1 and 5.2 MCM per year, respectively.7.
Water management strategies (S1-S4) lead to reduced water demand (by changing the crop pattern) and increased water supply (by increasing the irrigation efficiency), thereby mitigating the effects of climate change.8.
In the context of B1S4, the balance between agricultural sustainability (water challenge for food) and environmental sustainability (water challenge for the environment) will be achieved.9.
For more sustainability of Urmia Lake and reduction of its decline under other scenarios, other solutions can be considered, such as reducing the area of cultivated lands, conveying water between basins, cutting off the domestic water of Tabriz city and allocating it from another basin, and taking into account the consequences of the solutions.10.The water management strategies discussed in this study are preliminary, as they did not consider the socio-economic consequences.Changing the crop pattern might reduce farmer's income, while increasing the irrigation efficiency using new technologies might increase the costs.It is therefore suggested that the economic aspects of the proposed water management strategies and the benefits of Urmia Lake restoration be investigated in subsequent studies.

Figure 1 .
Figure 1.The map of the studied area.

Figure 1 .
Figure 1.The map of the studied area.
is the interflow from the first soil layer for area j; fj is the partition coefficient related to the land cover type, soil, and topography for area j, which divides flow into horizontal fj and vertical (1 − fj) flows; and ks,j is the saturated hydraulic conductivity of the root zone layer for area j [mm/time].

Figure 2 .
Figure 2. Conceptual diagram and equations incorporated in the soil moisture method [30].

Figure 2 .
Figure 2. Conceptual diagram and equations incorporated in the soil moisture method [30].

Figure 3 .
Figure 3. Schematic of WEAP model for Zarrinehrud and Siminehrud River basins.This modeling includes: • Two main rivers (Zarrinehrud and Siminehrud), two canals of water diversion (the right bank canal (RPC) and the Left bank canal (MC) of the Zarrinehrud irrigation network); • One reservoir dam (Shahid Kazemi Boukan dam); • Five aquifers (Saghez, ShahinDezh, Miandoab, Boukan, and Hajiabad); • Six domestic demand points (Tekab, Saghez, Boukan, Shahin Dezh, Miandoab, and Tabriz (in the form of inter-basin conveyance from Zarrinehrud basin to Tabriz city with the capacity of about 150 MCM per year)); • Eleven catchments (four non-agriculture catchments including upstream and downstream of Zarrinehrud basin and upstream and downstream of Siminehrud basin, and seven agriculture catchments including Tekab, Saghez, Shahin Dezh, right bank and left bank of Zarrinehrud irrigation network, Boukan 1, and Boukan 2); • Twenty-three lines of runoff/infiltration; • Nineteen transmission links; • Five return flows; • Four hydrometric gauge stations (including the entrance of Boukan Dam and Nezam Abad located on the Zarrinehrud River and Dashband and Miandoab Bridge on the Siminehrud River); and • Two points of the environmental flow requirements of Lake Urmia (from Zarrinehrud (ZFR) and Siminehrud (SFR) Rivers).

Figure 3 .
Figure 3. Schematic of WEAP model for Zarrinehrud and Siminehrud River basins.

Water 2018 ,
10, x FOR PEER REVIEW 11 of 25

Figure 5 .
Figure 5.Comparison of mean precipitation of the studied stations in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 6 .
Figure 6.Comparison of mean minimum and maximum temperature of the studied stations in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 5 .
Figure 5.Comparison of mean precipitation of the studied stations in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 5 .
Figure 5.Comparison of mean precipitation of the studied stations in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 6 .
Figure 6.Comparison of mean minimum and maximum temperature of the studied stations in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 6 .
Figure 6.Comparison of mean minimum and maximum temperature of the studied stations in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

3 .
Impact of Climate Change on the Hydrology in the Study Area

Figure 7 .
Figure 7.Comparison of mean streamflow in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 8 .
Figure 8.Comparison of mean actual evapotranspiration (ETc) in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 7 .
Figure 7.Comparison of mean streamflow in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 7 .
Figure 7.Comparison of mean streamflow in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 8 .
Figure 8.Comparison of mean actual evapotranspiration (ETc) in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

Figure 8 .
Figure 8.Comparison of mean actual evapotranspiration (ETc) in the future period (2015-2040) under different emission scenarios as compared to the baseline period (1986-2010).

4 .
Evaluation of Scenarios 3.4.1.Performance Criteria Performance criteria of reliability, resilience and vulnerability for agriculture catchments are presented in Table

Table 3 .
Performance evaluation of the LARS-WG model in the basic period using different indices.

Table 3 .
Performance evaluation of the LARS-WG model in the basic period using different indices.

Table 3 .
Performance evaluation of the LARS-WG model in the basic period using different indices.

Table 4 .
Performance criteria of reliability, resilience and vulnerability for agriculture catchments.

Table 5 .
Performance criteria of reliability, resilience and vulnerability for environment.

Table 6 .
Evaluation of agricultural sustainability indices and irrigation water deficit for agriculture catchments.