The Impacts of Climate Change on Water Resources and Crop Production in an Arid Region

: Climate change is one of the most pressing global issues of the twenty-ﬁrst century. This phenomenon has an increasingly severe impact on water resources and crop production. The main purpose of this study is to evaluate the impact of climate change on water resources, crop production, and agricultural sustainability in an arid environment in Iran. To this end, the study constructs a new integrated climate-hydrological-economic model to assess the impact of future climate change on water resources and crop production. Furthermore, the agricultural sustainability is evaluated using the multicriteria decision making (MCDM) technique in the context of climate change. The ﬁndings regarding the prediction of climate variables show that the minimum and maximum temperatures are expected to increase by about 5.88% and 6.05%, respectively, while precipitation would decrease by approximately 30.68%. The results of the research reveal that water availability will decrease by about 13.79–15.45% under different climate scenarios. Additionally, the ﬁndings show that in the majority of cases crop production will reduce in response to climate scenarios so that rainfed wheat will experience the greatest decline (approximately 59.95%). The results of the MCDM model show that climate change can have adverse effects on economic and environmental aspects and, consequently, on the sustainability of the agricultural system of the study area. Our ﬁndings can inform policymakers on effective strategies for mitigating the consequences of climate change on water resources and agricultural production in dry regions.


Introduction
Climate change's impact on agricultural production has raised serious global concerns in the twenty-first century [1,2]. Carbon dioxide (CO 2 ) emissions have been identified as the primary cause of climate change [3][4][5]. Human activities, such as the use of fossil fuels, environmental degradation, and land-use changes, have all contributed significantly to rising CO 2 levels in the atmosphere since the industrial revolution [6][7][8][9]. The global CO 2 concentration in the atmosphere increased from 288 ppm in 1750 to 415 ppm in 2021 [10]. This has resulted in higher global average temperatures and unpredictability of rainfall patterns [11]. Climate change and variability have far-reaching consequences for natural resources, human communities, and biodiversity [12]. Water resources and agriculture are most affected by climate change because it directly determines the availability of resources in terms of time and space [13]. Some researchers have concluded that climate change has a negative impact on groundwater table recharge, which affects irrigation [12,14]. Others believe that this phenomenon will increase agricultural water demand due to increased The literature review indicated that the effects of climate change on water resources and crop production have been studied all over the world. Xiong et al. [24] used climate scenarios of the regional climate model to investigate the consequences of climate change on water availability and cereal production in China in the 2020s and 2040s. The findings of this study revealed that there are insufficient water resources for cereal production, particularly in southern China, due to an increase in nonagricultural water demand and the occurrence of climate change. Palazzoli et al. [16] developed a soil and water assessment tool (SWAT) model to investigate the effects of future climate change on rainfed crop productivity and water resources in Nepal. Based on their results, they predicted significant potential changes in water resources availability (from −26 to +37%) and crop production (rice from −17 to +12%, wheat from −36 to +18% and maize from −17 to +4%). Sinnarong et al. [25] applied an econometric model to estimate the effect of climate change on rice production in Thailand. The results showed that temperature has a negative impact on rice production while precipitation has different regional effects on rice production. Additionally, the findings indicated that rice production under different climatic scenarios would decrease between 4.56% and 33.77%. Mostafa et al. [26] evaluated the impact of climate change on water resources and the agricultural sector of Egypt using climate and irrigation (CROPWAT) models and found the irrigation water requirement for wheat crop would rise by about 6.2% in 2050 and 11.8% in 2100. Furthermore, wheat production would decrease by approximately 8.6% and 11.1% in 2050 and 2100, respectively. Medellín-Azuara et al. [27] estimated the impact of climatic change on crop farming in California using the statewide agricultural production model (SWAP). They found that, by 2050, water supply, agricultural land use, and production of most crops will decrease in California due to climate change, such as rising temperature and declining precipitation. Shahvari et al. [28], using the SWAT model, assessed the impact of climate change on water resources and crop yield in Iran for the future. Their results revealed that future climate scenarios will lead to an increase in runoff in spring and autumn and a decrease in summer and winter. In addition, future climate change will reduce the yield of rainfed crops in the region. Lu et al. [29] constructed a new climate-economic model to analyze the effects of climate change on grain production and water resources. The findings of this study showed that irrigation water consumption has increased by about 100 billion m 3 because of climate change in China. This phenomenon also reduced the grain yield in this area by 1000 kg/hm 2 . The current gap in the existing literature is a comprehensive view of all meteorological, hydrological, economic, and sustainability aspects of climate change in the agricultural sector.
The study, therefore, assesses the effect of future climate change on water resources, crop production, and agricultural sustainability in the Mashhad plain under three climate scenarios (RCP 2.6, 4.5, and 8.5). Specifically, the study aims at (1) projecting climate variables using the Long Ashton Research Station Weather Generator (LARS-WG) model alongside HadGEM2-ES outputs and the three RCP scenarios; (2) assessing the impact of future climate change on water resources in the Mashhad plain using panel data model; (3) estimating the relationship between crop yield and climate variables, including minimum temperature, maximum temperature, and precipitation, using the GME (generalized maximum entropy) technique; (4) investigating the impact of climate change on cropping pattern, crop production, and water consumption in the selected area using PMP (positive mathematical programming) model; (5) evaluating the agricultural sustainability under climate scenarios using a MCDM (multicriteria decision making) method. The results of this study are expected to provide policymakers with insights into designing climate change mitigation policies to ensure food security and sustainable production in the region.
The contribution of this study to the literature is twofold. First, to the best of our knowledge, this is the first attempt to apply an integrated climate-hydrological-economic model to evaluate the effect of climate change on water resources, crop production, and cropping pattern in Iran. The second contribution of the study is found in the use of the MCDM approach to investigate the sustainability of agricultural activity at the regional level under different climate scenarios.
The study is structured as follows: Section 2 describes the study area, datasets, and methodology. Section 3 presents the results and discussions, and the last section concludes with the research and policy implications of these findings.

Materials and Methods
This study used a variety of methods to achieve its research goals. The LARS-WG model was used to downscale the climate variables (maximum and minimum temperatures and precipitation), and a regression model was used to forecast the groundwater availability in the Mashhad plain. The sensitivity of yield to climate change was estimated using the GME technique. The cropping pattern was then evaluated under climate change using a PMP model. Finally, economic, social, and environmental indicators were ranked using an integrated MCDM method. The complete structure of the framework is presented in Figure 2.

Study Area
The present case study is in Northeast Khorasan Razavi Province, Iran, between the latitude 35 • 59 to 37 • 03 N and longitude 58 • 22 to 60 • 06 E, covering an area of approximately 9957 km 2 . The plain is bounded on the north by Hezar Masjid heights, on the northwest by Atrak river basin, on the south by Binaloud Mountain, and on the southeast by Jamroud river basin [30]. It has a semi-arid to arid climate, with an average annual rainfall from 1991 to 2015 is about 262 mm [31]. The average monthly temperature in this plain is reported to be between 11.6 • C and 26.7 • C. Furthermore, the average annual evapotranspiration ranges from 236 to 310 (mm). The location of the study area is shown in Figure 3. Around 3 million people live in the study basin and rely mostly on groundwater resources for drinking and agricultural cultivation. Climate change and a lack of water resource management in this area have resulted in a 12-m water drop in the water table over a 20-year period [32].

Data Collection
The observed daily time series data of minimum temperature, maximum temperature, and precipitation for Mashhad synoptic station during the period 1979-2016 were obtained from Iran's Meteorological Organization. To estimate the model (Equation (6) Figure 4 shows the growing seasons of crops (wheat, barley, alfalfa, potato, corn, tomato, melon, onion, sugar beet, and cucumber) in the Mashhad plain. Data and information on outputs prices, inputs costs, technical coefficients, crop yield, and resources availability were gathered through face-to-face interviews with farmers in 2016-2017 cropping season.

Meteorological Model
The LARS-WG model is a random generator that uses statistical downscaling techniques to generate meteorological data [33]. Because of the repeated calculations, it requires less input data and is simpler and more efficient than other models [34,35]. Racsko et al. [36] proposed this model, which Semenov et al. [37] later revised and developed. This model's sixth version (LARS-WG6) was updated in 2018 for downscaling the coupled model intercomparison project phase (CIMP5) [34]. The HadGEM2-ES universal model data were used in this study to project climate variables during the horizon in 2045 under three climate scenarios, namely RCP 2.6, RCP 4.5, and RCP 8.5. The LARS-WG model was implemented using daily data of maximum and minimum temperature and precipitation from 1979 to 2016.
The performance of the LARS-WG statistical model was evaluated by comparing the simulated and observed maximum and minimum temperatures, as well as precipitation, using the following statistics [2,38,39]: coefficient of determination (R 2 ), normalized root mean square error (NRMSE), root mean square error (RMSE), mean absolute deviation (MAD), and mean square error (MSE) (see Table 1).
the lower values show a better model [40] MAD Note: N is the number of data points, S n is the simulated values, O n is the observed values, O is the mean of the observed values, and S is the mean of the simulated values.

Hydrological Model
This section investigates the effects of climate change on groundwater resources in the Mashhad plain. The panel data model was used to forecast the groundwater depth in the plain. To this end, the following groundwater conceptual model was estimated using panel data regression [41,42] related to piezometers as cross units from 2000 to 2016. lnH t = α 0 + α 1 lnH t−1 + α 2 lnP t−1 + α 3 lnTmin t−1 + α 4 lnTmax t−1 (1) where t indicates time, H t is predicted groundwater depth, H t−1 is the initial groundwater depth, P is total monthly precipitation (mm), Tmin and Tmax are monthly minimum and maximum temperature ( • C), respectively, and α 0 to α 4 are the model coefficients that should be estimated. After estimating the amount of changes in groundwater depth because of climate change, Equation (2) was used to calculate the amount of changes in groundwater resources [43][44][45][46].
where ∆V is the groundwater storage change (m 3 ), A is geographical area (m 2 ), S y is specific yield (dimensionless), and ∆H is average depth change (m).

Generalized Maximum Entropy (GME) Model
In this study, the Cobb-Douglas Yield Response (CDYR) model was used to assess the sensitivity of crops yield to climatic variables, such as minimum and maximum temperatures, as well as precipitation [47,48]. The CDYR model, after taking the logarithm of both sides of the equation, can be presented as follows: where t is the set of years; Y is the yield of crops (wheat, barley, alfalfa, potato, corn, tomato, melon, onion, sugar beet, and cucumber); Tmin and Tmax are average minimum and maximum temperature; and P is total growing season precipitation. It is worth noting that the impact of additional factors affecting crop production, which were not incorporated in Equation (3), was covered in the residual terms [49]. Due to the lack of access to crop yields data in the research area, we only analyzed data of 34 years (1983-2016). In this case, because of the limited sample size, traditional estimation approaches, such as ordinary least squares (OLS), may result in parameter estimation with excessive variance [49,50]. To address this issue, following Moreno et al. [51], we used the GME estimator. The GME approaches are founded on Shannon's entropy information measure and the generalized maximum entropy theory [50,52]. Instead of calculating the mean and variance of coefficients directly, GME estimator considers a probability distribution for the coefficients and error terms [52]. Assume that y is dependent on K independent variables x k (k = 1, . . . , k): where y is a (T × 1) vector of observations for y, X is a (T × K) matrix of observations for the x k variables, β is the (K × 1) vector of explanatory variable coefficients, and ε is a (T × 1) vector of residual terms. To estimate explanatory variable coefficients using GME, firstly we reparametrize the regression model, and then we recast the coefficients and residual in terms of discrete probability distributions.
In this method, each β k is assumed to be as a discrete point with M dimension (M ≥ 2). Let Z k = [Z k1 , . . . , Z kM ] be support points for parameter β k , which are symmetrical around zero. Additionally, the probability mass function of z k is defined as P k = [P k1 , . . . , P kM ] such that: Then, β can be represented as follows: where Z is a (K × KM) matrix of support values, and P is a (KM × 1) of vector of unknown weights. The unknown error is defined as follows: where w g = w g1 , . . . , w gJ is a vector of weights, andv g = v g1 , . . . , v gJ (J ≥ 2) is a set of support points. The error vector is presented as follows: Now, Equation (9) can be extended as follows: In Equation (9), y, Z, and V are known vectors, and P and w are unknown vectors that are estimated using GME, which is defined as follows: w gj ln w gj (10) and is subject to: The GME techniques are created by solving the optimization problem (Equation (10)) while taking constraint into account (Equations (11)(12)(13)). Equation (11) is a condition for the compatibility of the probability of the posterior distribution of the coefficients and the residual terms with the observations. Equations (12) and (13) are normalization constraints for the probabilities.

Positive Mathematical Programming (PMP) Model
The present study used an economic modeling system composed of the PMP model to analyze and assess the effects of climate change (the decline in available water and the changes in crop yields) on the cropping pattern of selected crops and farmers' gross revenues in Mashhad plain.
In recent decades, PMP has been widely used to evaluate the effects of climate change on the agricultural sectors [53][54][55][56][57][58]. The main objective of the PMP is to improve the accuracy of modeling farmers' behavior in the context of an optimization model utilizing observed values from the baseline year [59]. There are two primary reasons for interest in this approach: firstly, in the presence of incomplete and insufficient data, alternative approaches, such as traditional econometrics, are unable to model farmer behavior; secondly, optimization models cannot properly calibrate farm-level models [60]. In the current study, the impacts of changes in climate variables, such as temperature and rainfall on crop pattern, were simulated in the framework of a developed PMP. The study's empirical model comprises a nonlinear objective function as well as constraints, such as water, labor force, fertilizer, and land. Following Röhm and Dabbert [61] and Radmehr and Shayanmehr [59], PMP is constructed in three stages: (1) solve a linear optimization programming model and obtain shadow prices, (2) use a generalized maximum entropy (GME) approach to calibrate crop yield parameters, (3) solve a nonlinear optimization programming model that includes the objective function and constraints (from step one), as well as calibrated yield functions (obtained in the second step). In order to simplify, the nonlinear optimization model developed in the third step is presented in this section as follows: Max TMG = ∑ e ∑ r P e α e,r X e,r − β e,r X 2 e,r − ∑ e ∑ r C e,r X e,r (14) and is subject to: ∑ r=land ∑ e w e X e,r ef ≤ b r water (15) ∑ e a e,r X e,r ≤ b r land, labor, fertilizer, and machinery (16) X e,r ≥ 0 (17) In this expression, e is the set of different crops; r is the set of production inputs (land, water, labor, machinery, and fertilizers); TMG shows the total gross margin in the region; X e,r is a decision variable that represents the amount of input r used for crop e; α e,r and β e,r are coefficients of yield function that calibrated using the GME approach (for more information about the details of the technique, see Paris and Howitt [62]); C e,r is the unit cost of input r for crop e; w e is the water requirement of crop e; ef is technical efficiency of irrigation water use; a e,r is the technical coefficient of input r for crop e, which shows the amount of input r to produce a unit of crop e; b r is the total available input r. Equation (14) indicates the objective functions that maximize the total gross margin of production in the irrigated area. Equation (15) is the constraint of water that represented the amount of water allocated for agricultural production and should be less than total water availability for crop production in the region. Equation (16) is the constraint of inputs that shows the amount of input allocated for crop production to be less than total input availability for crop production in the region. Finally, the non-negative constraint (Equation (17)) states that the decision variable (X e,r ) must be greater than or equal to zero.

Multi-Criteria Decision-Making Approach
The current study examines the effects of climate change on agricultural sustainability using three types of indicators: social, environmental, and economic indicators. The social index is based on the sub-index of farm employment (FE), while the environmental index includes the sub-indices of nitrogen balance (NB), phosphorus balance (PB), and water consumption (WC). The economic index is introduced using the sub-indices of total gross margin (GM) and profit-to-water consumption ratio (PW) [59,[63][64][65]. Many complex decision-making issues employ MCDM models [66]. Analytic hierarchy process (AHP) and technique for order preferred preference by similarity to ideal solution (TOPSIS) are two approaches of this model that have been used in numerous studies [63,[67][68][69]. The main advantages of AHP method are the ability to (i) depict the rationale of human choice; (ii) evaluate the relative performance of alternatives based on the simple algorithm; and (iii) define flexibly the selection set [70][71][72].
To select and rank indicators, an integrated AHP and TOPSIS method is used. This method consists of the eight steps listed below: Step 1. Build a decision matrix.
First, a decision matrix is created, which can be expressed as follows: where A i is the alternative; F j is the evaluation indicators; and f ij is the performance value of A i with respect of F j .
Step 2. Construct the normalized decision matrix (r ij ) using following formula: Step 3. Compute the weight (w j0 ) of the indicators.
The relative importance of various indicators is determined with respect to the objective, and weights of indicators are given based on their importance.
K indicates an n × n pair-wise comparison matrix: In an arbitrary random reciprocal matrix K, each criterion k ij is the relative importance of i th alternatives compared to the j th indicators [73]. Therefore, it expresses that the higher values of k ij show stronger preference of k i over k j . In the matrix k, k ij = 1, when i = j and k ji = 1/k ij .
Geometric mean method is employed for normalization and determines the importance degree of the indicators [74]. If W i indicates the importance degree for the i th attribute, then: E indicates an n-dimensional column vector, which defines the sum of the weighted values of the importance degree of indicators. Then:  (24): λ max is maximum Eigen value that can be obtained as follow [74]: The consistency of evaluation in AHP is measured by consistency ratio (CR). Consistency ratio is defined as Equation (26): where RI indicates the inconsistency index of a random matrix. If the value of consistency ratio is less than 0.10, the evaluation of the importance of degrees of attributes is acceptable.
where J and J are the indicators with positive and negative polarity, respectively.
Step 6. Compute the relative distance of each A i from A + and A − [63].
Step 7. Determine the relative closeness (C i ) to the best alternative [74].
Step 8. Rank the alternatives.
The alternative that has the highest value of C i is selected as the best option.

Projecting Climate Variables
Using the LARS-WG model, climatic parameters of maximum and minimum air temperature and precipitation were predicted in the Mashhad plain based on data of 1979-2016. To calibrate and ensure the accuracy of the LARS-WG model, the simulated data are compared with the observed data on a monthly scale, as shown in Figure 5. Assessment of the monthly average of maximum air temperature, minimum air temperature, and precipitation shows a good agreement among all three parameters. LARS-WG performance also was investigated using R 2 , NRMSE, RMSE, MAD, and MSE indicators, as presented in Table 2. The model successfully downscaled the minimum and maximum temperature as well as precipitation, according to the evaluation of these indicators. The higher values of R 2 (>0.98) and the lower values of RMSE (0.21-2.09), MAD (0.17-1.69), MSE (0.04-4.39), and NRMSE (0.95-9.96) for this period reveal that the simulated precipitation and temperature data are acceptable. After evaluating the accuracy of the model, the climate scenarios were generated by downscaling HadGEM2 outputs under climate scenarios of RCP 2.6, 4.5, and 8.5 on the horizons of 2045. (We considered a 30-year period (2016-2045) to investigate the effects of climate change on crop production. Because the term climate is a long-term shift in the weather pattern, it is an average of at least 30 years of weather condition of a particular place [75][76][77].) The percentage change in climate variables was then calculated and compared to the benchmark year (2016). The obtained findings are shown in Table 3. According to the results, the minimum and maximum temperatures are expected to increase by about 5.88% and 6.05%, respectively, while precipitation would decrease by approximately 30.68%.

Evaluating the Impacts of Climate Change on Water Resources
The groundwater depth was forecasted because of climatic change using a panel data model under three scenarios in the Mashhad plain. The first step in the analysis was to determine whether the variables were stationary. This was performed using the Im-Pesaran-Shin (IPS) and ADF-Fisher-type panel unit root tests. Table 4 shows, for all variables in the model, the null hypothesis of unit roots is rejected with a 99% confidence level. As a result, the model's variables were all stationary. For panel data, the random-effects and fixedeffects models were used. To identify which should be used, the Hausman test was used. As shown in Table 4, the null hypothesis of no correlation between regional effects and independent variables is rejected. As a result, a fixed-effect model with regional-specific effects was used.
The results of the panel data model for sensitivity of the groundwater depth to climate variables are presented in Table 5. The findings indicated that maximum and minimum temperatures had a positive impact on the groundwater depth in the Mashhad plain from 2000 to 2016. Furthermore, precipitation was negatively and significantly related to groundwater depth in the plain. Therefore, as precipitation increased, the groundwater depth decreased, resulting in more water in the well, according to many previous studies, including Shahvari et al. [28] and Izady et al. [41].  The results regarding the percentage changes in groundwater depth and water availability under climatic scenarios compared to the baseline are presented in Table 6. According to the findings, the depth of groundwater in the Mashhad plain is expected to rise about 13.79% in RCP 2.6, 13.50% in RCP 4.5, and 15.45% in RCP 8.5. Furthermore, water availability will increase by approximately 13.25%, 13.26%, and 14.84% in response to RCP 2.6, RCP 4.5, and RCP 8.5, respectively.

Assessing the Impacts of Climate Change on Crop Yield
To assess crop yield sensitivity to temperature and precipitation, the yield response function was estimated using the GME model. The estimated results are displayed in Table 7. The Cobb-Douglas functional form was used in the estimation of yield response functions. Therefore, the estimated parameters in Table 5 show the elasticity values. Based on the results obtained from this table, increasing maximum and minimum temperatures reduce crop yield in many crops. In the case of irrigated wheat, for example, the results show that a 1% increase in maximum temperature results in a 1.17% decrease in yield. In addition, a 1% increase in the minimum temperature increases irrigated wheat yield by approximately 0.76%. Precipitation has a negative impact on yield of some crops due to increase humidity or the potential spread of diseases and pests [78]. Irrigated wheat, rainfed wheat, irrigated barley, rainfed barley, alfalfa, corn, sugar beet, melon, cucumber, and tomato yield are positively influenced by precipitation. Table 7. Estimates of the impact of climatic variables on crop yield using the GME model.  Table 8. The findings imply that the yield of irrigated wheat, rainfed wheat, irrigated barley, rainfed barley, alfalfa, sugar beet, onion, and cucumber will decrease in response to all three climate scenarios. Alfalfa crop is prone to experience a decrease in yield between 19.07% and 25.20% and is expected to emerge as a highly vulnerable crop in 2045. In addition, climate change will increase the yield of corn, potato, melon, and tomato crops under all scenarios. With changing climate, corn yield will rise more than other crops. This result is in line with the results of Almaraz et al. [79] and Zhang et al. [80].

Evaluating the Impacts of Climate Change on Crop Production and Cropping Pattern
In this section, changes in crop yields and water availability under climate change were incorporated into the PMP model to assess the impact of climate change on cropping pattern, crop production, and water consumption. It is worth noting that the cropping pattern is defined as a combination of agricultural crops that are grown in a particular geographical area [81]. The percentage changes in cropland under three climate scenarios in 2045, as compared to the base year, are shown in Figure 6. The results imply that RCP 2.6, RCP 4.5, and RCP 8.5 scenarios will decrease the cultivated area of irrigated wheat, rainfed wheat, irrigated barley, rainfed barley, alfalfa, sugar beet, onion, cucumber, and tomato in 2045 in the Mashhad plain. On the contrary, the cultivated area of corn, potato, and melon will increase under two scenarios of RCP 4.5 and RCP 8.5 in 2045 relative to the base year. As shown in the figure, the most significant reduction in cultivated land is related to rainfed wheat and RCP 4.5, where the area under cultivation will reduce by approximately 51.16%. The percentage changes in crop production under climate scenarios for 2045, as compared to the baseline, are presented in Table 9. The results of this table show that production of crops, such as irrigated wheat, rainfed wheat, irrigated barley, rainfed barley, alfalfa, sugar beet, onion, and cucumber, will decrease in response to all three climate scenarios. Therefore, the biggest decline (59.95%) will occur for rainfed wheat in RCP 4.5 scenario. In addition, corn, potato, and melon production will increase by about 17.92%, 5.33%, and 7.14%, respectively, under climate scenarios. This increase in production is due to an increase in yield or area under cultivation, or both, as discussed in the previous section. Additionally, in the presence of climate change in 2045, the Mashhad plain will experience the biggest increase in crop production in corn and RCP 4.5.
Given the decline in production of most crops that account for a large share of the region's production, it is reasonable to conclude that the occurrence of climate change poses a serious threat to crop production and food security in the region. Figure 7 shows the water consumption for each crop under current conditions and different climatic scenarios. As expected, under climate scenarios, due to the reduction in water availability and crop yield and, consequently, reduced cropland, water consumption will decrease for crops, such as irrigated wheat (21.40%), irrigated barley (21.88%), alfalfa (27.70%), sugar beet (7.90%), onion (8.71%), cucumber (4.52%), and tomato (4.33%), while there is an increase in water consumption for corn, potato, and melon because of improved yield. Overall, it can be stated that water consumption in climatic scenarios will be reduced by about 13.30%, 13.31%, and 14.90% under RCP 2.6, RCP 4.5, and RCP 8.5, respectively, compared to baseline conditions.

Assessing the Impacts of Climate Change on Agricultural Sustainability
In the present study, an attempt has been made to evaluate the effects of climate change on agricultural sustainability in the study area using the TOPSIS approach. In the first step, we employed the AHP method to determine the weight of indicators and sub-indicators through interviews with 15 agricultural experts and specialists. Table 10 shows the weights assigned to each indicator and sub-indicator. According to the results of this table, economic, environmental, and social indicators have the greatest importance with 52%, 33%, and 14%, respectively. Due to the low income of farmers in Iran, improving the profitability of agricultural activities is critical to achieving sustainable development in the agricultural sector [60]. In addition, over the past few decades, the excessive use of chemical fertilizers and the uncontrolled extraction of groundwater resources in the process of food production in agriculture has caused irreversible environmental damage [16]. This highlights the importance of environmental issues in assessing the sustainability of Iran's agricultural sector. In the economic indicator, the importance of the "GM" and "PW" sub-indicators are 58% and 42%, respectively. Because of Iran's high unemployment rate, the FE index was deemed as the mere social indicator, with a weight of 14%. Among the environmental sub-indicators, WC is the most important environmental sub-indicator, followed by the NI and PB. In the second step, using the results obtained from the PMP model, the value of each of the sustainability indicators in the base conditions and climatic scenarios was calculated, which forms the decision matrix (see Table 11). Table 11 shows under the climate scenarios that the values of economic indicators are lower and environmental indicators are higher than the baseline conditions. Table 12 indicates the normalized decision matrix. Figure 8 depicts the distribution of positive and negative ideal options across the various sustainability indicators. In the final step, current and climatic conditions were ranked based on sustainability indicators (see Table 13). Based on the results of Table 13, the base condition has the highest Ci values (0.77), indicating that the phenomenon of climate change can be considered a serious threat to the agricultural sustainability in this region. The findings in this section of the research are in line with the study of Karandish et al. [82] and Ghanian et al. [83].    To ensure consistency in the results of Table 13, a sensitivity analysis of the economic, social, and environmental criteria weights used in the TOPSIS model is performed. To achieve this goal, a total of seven experiments has been conducted to compare the impact of potential changes in the weights of economic, social, and environmental criteria (see Table 14). In Experiment 1, all sub-indicators have the same weight (16.66%); in Experiment 2-7, the weight of one criterion is higher than the weight of the remaining criteria. The results of the sensitivity analysis described indicate that in all experiments (except Experiment 4), the base condition is selected as the best alternative. In other words, given the importance of economic and environmental indicators in Iran's agricultural sector, it is reasonable to expect that if the harmful effects of climate change are not properly managed, the phenomenon will have a negative impact on agricultural sustainability.

Conclusions
The climate simulation model predicted a 6% increase in minimum and maximum temperatures, as well as a 30% decrease in precipitation. Additionally, the results showed that water availability will decrease between 13% and 15% under different climate scenarios. Crop yields were found to be negatively affected by increasing maximum and minimum temperatures. Precipitation affects crop yield in different ways, with positive effects on irrigated wheat, rainfed wheat, irrigated barley, rainfed barley, alfalfa, corn, sugar beet, melon, cucumber, and tomato yields and negative effects on potato and onion yields. Overall, future climate change is expected to reduce the yield of irrigated wheat, rainfed wheat, irrigated barley, rainfed barley, alfalfa, sugar beet, onion, and cucumber, while the effects will be reversed for corn, potato, melon, and tomato. The results of the PMP model showed that changes in crops yield and water availability will lead to a reduction in the cultivated area of most crops in 2045, among which dryland wheat will experience the greatest decrease (51%). The results of the evaluation of the effects of climate change on agricultural sustainability show that this phenomenon can have adverse economic and environmental effects on the agricultural system of the study area. As a result, it can have a negative impact on agricultural sustainability if not properly managed.
In entirety, the findings of this study reflect the fact that water and food security in the region will be severely adversely affected by climate change in the future. Nevertheless, by continuing to support population growth policies, uncontrolled extraction of groundwater, and expansion of urbanization in the presence of climate change, more severe irreversible effects on water and food resources in the control area are expected. These results underscore the necessity of implementing adaptation policies, such as reforming the cropping pattern and production technologies, as well as the introduction of drought-tolerant varieties to reduce the detrimental effects of climate in the region.