Unsustainability Syndrome—From Meteorological to Agricultural Drought in Arid and Semi-Arid Regions

: Water is the most important resource for sustainable agriculture in arid and semi-arid regions, where agriculture is the mainstay for rural societies. By relating the water usage to renewable water resources, we deﬁne three stages from sustainable to unsustainable water resources: (1) sustainable, where water use is matched by renewable water capacity, ensuring sustainable water resources; (2) transitional, where water use occasionally exceeds renewable water capacity; and (3) unsustainable, with lack of water resources for agriculture, society, and the environment. Using available drought indicators (standardized precipitation index (SPI) and streamﬂow drought index (SDI)) and two new indices for agricultural drought (overall agricultural drought index (OADI) and agricultural drought index (ADI)), we evaluated these stages using the example of Fars province in southern Iran in the period 1977–2016. A hyper-arid climate prevailed for an average of 32% of the province’s spatio-temporal coverage during the study period. The area increased signiﬁcantly from 30.6% in the ﬁrst decade (1977–1986) to 44.4% in the last (2006–2015). The spatiotemporal distribution of meteorological drought showed no signiﬁcant negative trends in annual precipitation during 1977–2016, but the occurrence of hydrological droughts increased signiﬁcantly in the period 1997–2016. The expansion of irrigated area, with more than 60% of rainfed agriculture replaced by irrigated agriculture (especially between 1997 and 2006), exerted substantial pressure on surface water and groundwater resources. Together, climate change, reduced river ﬂow, and signiﬁcant declines in groundwater level in major aquifers led to unsustainable use of water resources, a considerable reduction in irrigated area, and unsustainability in agricultural production in the period 2006–2015. Analysis of causes and e ﬀ ects of meteorological, hydrological, and agricultural drought in the area identiﬁed three clear stages: before 1997 being sustainable, 1997–2006 being transitional, and after 2006 being unsustainable.


Introduction
Due to global climate change, the extent, intensity, and frequency of droughts in arid and semi-arid regions are increasing [1]. Droughts affect large areas and last a relatively long time compared with other natural disasters, such as floods and landslides [2]. Dryness characterizes arid and semi-arid To assess the consequences of hydrological and meteorological droughts and evaluate the possibility of these leading to agricultural drought, we analyzed changes in rainfed and irrigated farmland over the 40-year study period ). We also analyzed groundwater level variations from 1995 to 2015 in 12 major aquifers in the province's different regions (Figure 1d), as supplementary and supporting information for drought analysis. All hydro-climatological data, including river flow, groundwater levels, precipitation, and temperature, were provided by Fars Regional Water Authority. Agricultural data were obtained from the national and regional agricultural census 1972-2016 provided by the ministry of Agricultural-Jahad. Although different climates prevail in different parts of the province, the dominant climate is hyper-arid/temperate. This creates pressure on the agriculture sector, due to its high dependency on water resources. The expanding agriculture sector, river regulation, and increasing frequency of meteorological droughts have led to a reduction in available surface water in recent years [29], with major lakes in the province drying up [29,45] Thus, farmers are increasingly dependent on groundwater resources. As a result, a significant decrease in groundwater levels has been recorded [40,46], and river water levels are critically low [41,47]. For spatial analysis and to assess different types of droughts, we divided Fars province into six regions (southern, western, eastern, central, northwestern, and northeastern), based on the distribution of major cities (Figure 1b). To monitor variations in climate and meteorological drought, we obtained monthly precipitation and temperature data for the period 1977-2016 from 38 meteorological stations in Fars and a neighboring province (Figure 1b). A full list of stations is available in Appendix A, Table A1). Hydrological drought was assessed based on flow data from 24 gauging stations on 17 rivers in six main basins in Fars province (Maharloo-Bakhtegan, Sirjan desert, Mond, Heleh, Kol-o-Mehran, and Zohreh) ( Figure 1c). Geographically, the main rivers are located in the central, northwestern, and western regions of the province (Table A2 in Appendix A). Of these six basins, only Maharloo-Bakhtegan falls entirely within Fars province, and its main river, the Kor, discharges to Lake Bakhtegan (Figure 1c). The rivers in the Kol-o-Mehran, Mond, Heleh, and Zohreh basins discharge to the Persian Gulf, whereas the Bavanat River in the Sirjan desert basin discharges to Kavir-e Sirjan.

Methods
To assess the consequences of hydrological and meteorological droughts and evaluate the possibility of these leading to agricultural drought, we analyzed changes in rainfed and irrigated farmland over the 40-year study period ). We also analyzed groundwater level variations from 1995 to 2015 in 12 major aquifers in the province's different regions (Figure 1d), as supplementary and supporting information for drought analysis. All hydro-climatological data, including river flow, groundwater levels, precipitation, and temperature, were provided by Fars Regional Water Authority. Agricultural data were obtained from the national and regional agricultural census 1972-2016 provided by the ministry of Agricultural-Jahad.

Methods
There are many methods and indices for identifying climate variability and its interrelation with water availability. Most of these methods measure drought and aridity, often linked to other related components, such as vegetation and soil moisture. The methodology for the present analysis was developed by combining Emberger climate classification, the Emberger aridity index, SPI, SDI, groundwater level fluctuation analysis, and agricultural development investigations. This helped identify agricultural development and changes in climate characteristics in the region. We applied the non-parametric Mann-Kendall test [48] for assessing trends in different types of drought in the study area.

Climate Variability Analysis (Emberger Method)
A method designed by Emberger [49] is widely used for climate classification in climatology, hydrology, ecology, and bioclimatic analysis [50][51][52]. The method considers the two important climate parameters, precipitation and temperature: where P is annual precipitation (mm), M is mean maximum temperature in the warmest month ( • K), and m is mean minimum temperature in the coldest month ( • K). Equation (1) considers the minimum temperature of the coldest month in the year, since vegetation growth is strictly related to thermal limits [49]. In addition, the rate of evaporation changes with temperature variations, so the term (M-m) represents evaporation and the continentality of a climate [50]. By varying m and Q in Equation (1), 36 different climate classes were defined and represented in the Emberger climatogram ( Figure 2).
To evaluate climate variability, the Emberger aridity index was calculated annually for each station. Then, by calculating three Emberger influence factors, i.e., humidity influence factor (HIF), thermal influence factor (TIF), and climate influence factor (CIF), the dominant humidity, thermal, and climate groups were defined in each region and the whole province for different decades (1977-1986, 1987-1996, 1997-2006, and 2007-2016) and for the whole study period .
where P is annual precipitation (mm), M is mean maximum temperature in the warmest month (°K), and m is mean minimum temperature in the coldest month (°K). Equation (1) considers the minimum temperature of the coldest month in the year, since vegetation growth is strictly related to thermal limits [49]. In addition, the rate of evaporation changes with temperature variations, so the term (M-m) represents evaporation and the continentality of a climate [50]. By varying m and Q in Equation (1), 36 different climate classes were defined and represented in the Emberger climatogram ( Figure 2).  There are six HIFs (Equation (2)), representing the influences of different humidity groups over an area or station. These HIFs are shown on the y-axis of the Emberger climatogram ( Figure 2). They comprise: HA (hyper-arid), AR (arid), SA (semi-arid), SU (sub-humid), HU (humid), and HH (hyper-humid), which are defined based on the value of Q in (Equation (1)): where HIF x is the humidity influence factor for group x (HA, AR, SA, SH, HU or SU) in a region or station, n is the number of meteorological stations in each region, NY x(i) is the number of years in which station i experienced humidity conditions of type x, and NY Hi is the number of years in which station i experienced humidity conditions of type H.
There are five TIFs (Equation (3)), representing the influences of different thermal groups over an area or station. These thermal groups are shown on the x-axis of the Emberger climatogram ( Figure 2). They comprise VC (very cold), CD (cold), CL (cool), TE (temperate), and HO (hot), which are defined based on the value of mean lowest temperature in the coldest month (m): where TIF y is thermal influence factor for group y (VC, CD, CL, TE, or HO) in a region or station, n is the number of meteorological stations in each region, NY (y)i is the number of years in which station i experienced thermal conditions of type y, and NY Ti is the number of years in which station i experienced thermal conditions of type T.
There are 30 CIFs (Equation (4)), representing the influence of each climate group over an area or station. Each climate group has two terms, where the first term indicates humidity group (HA, AR, SA, SH, HU, SU) and the second term indicates thermal group (VC, CD, CL, TE, HO); for example, HA-VC is hyper-arid-very cold: where CIF x-y is the climate influence factor of x and y (where x is a humidity group and y is a thermal group) in a region or station, H is a humidity group (1-6 are HA, AR, SA, SH, HU, and SU, respectively), T is a thermal group (1-5 are VC, CD, CL, TE, and HO, respectively), n is the number of meteorological stations in each region, N (x-y)i is the number of years in which station i experienced climate types of x and y in the selected region, and NY HTi is the number of years in which station i experienced humidity conditions of type H and thermal conditions of type T.

Meteorological, Hydrological, and Agricultural Drought Analysis
The standardized precipitation index is used to indicate meteorological drought, which is caused by a deficiency of precipitation [28,53]. To calculate SPI, the long-term precipitation data were fitted to a probability distribution, which was then transformed into a standard distribution (see Appendix B). Zero means normal conditions, while positive and negative values indicate wet and dry conditions, respectively (Table 1). Table 1. Standardized precipitation index (SPI) classification [29].
The streamflow drought index uses standardized annual streamflow volumes to define values lower than the mean streamflow by at least one standard deviation [15]. Based on the SDI values, the state of hydrological drought is then defined ( Table 2). SDI is calculated as: where F y is the annual stream flow for a given year (y), and F and S are the mean and standard deviation of cumulative streamflow for the whole study period. Table 2. Definition of states of hydrological drought with the streamflow drought index (SDI) [15].

Criterion Description of State
We developed two additional indices, the overall agricultural drought index (OADI) and agricultural drought index (ADI), to evaluate changes in agriculture and drought in the study area: where A y , A i is the annual farmed area for a given year (y), A and S are the mean and standard deviation of the farmed area for the whole study period, and A 1−y (1−y) and S (1−y) are the mean and standard deviation of farmed area between the first and yth year in the study period. We calculated these indices for irrigated, rainfed, and total farmed areas in Fars province. In addition to the indices for drought analysis, we used the non-parametric Mann-Kendall test [48] to determine the significance of a trend in rainfall, discharge, groundwater level, and farmed area. This test is widely used for assessing trends in hydrology and climatology data [29,54].  Figure A1 in Appendix A). The most variety in climate (17 types) was observed in the central region of the province, while the least was observed in the western and northeastern regions, eight types each. The dominant climate type in the province was HA-TE, hyper-arid and temperate, with a CIF value of 14.9% (Figures 2 and 3). Temporal analysis showed that the CIF for HA-TE climate increased from 12.7% to 20.8% between the first (1977-1986) and last (2007-2016) decades of the study period ( Figure 3, Figure A1 in Appendix A). The least common types of climate, CIF less than 0.3% ( Figure A1 in Appendix A), were HU-VC (humid and very cold) and HU-CD (humid and cold), which were only observed in the northwestern region, with CIF 0.4% and 0.9% respectively ( Figure A3 in Appendix A). The diversity of climate types in the northwestern, western, central, and southern regions decreased, while in eastern and northeastern regions it increased Based on the HIF values obtained, the dominant climate was HA, with 32.3% spatiotemporal coverage during the study period. Its dominance significantly increased from 30.6% in the first decade (1977)(1978)(1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)) to 44.4% in the last (2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015). The next most important humidity group was AR, HIF 31.4%, followed by SA, HIF 24.9%; meanwhile SH, HIF 11.1%, and HU, HIF 0.3%, were less well represented. All these groups showed a reduction in HIF over time ( Figure A1 in Appendix A). The HH type was not observed in any region. The HU type was only observed in the northwestern region, and the SH type in northwestern (HIF 32.6%) and central (16%) regions ( Figures A2 and A3 in Appendix A), whereas other humidity types, HA, AR, SA, and SU, were observed in all regions.

Spatiotemporal Climate Variability
Based on the TIF values, VC, very cold, was the dominant thermal group, 28.7% and then TE 26.5%, CL 21.3%, CD 19%, and HT 4.5%. The TIF for VC and CD showed an increase over time, whereas it decreased for the remaining thermal groups (Figures 3 and A2-A5 in Appendix A).

Meteorological Drought
The lowest annual value of SPI in Fars province and its regions was observed in 2008 ( Figure  4a). The last decade of the study period (2007-2016) experienced the most extreme droughts in different periods (1, 2, 3, 5, and 10 years) in the different regions. The longest drought period started

Meteorological Drought
The lowest annual value of SPI in Fars province and its regions was observed in 2008 ( Figure 4a). The last decade of the study period (2007-2016) experienced the most extreme droughts in different periods (1, 2, 3, 5, and 10 years) in the different regions. The longest drought period started in 2006 in the eastern and southern regions, and in 2007 in the other regions and the province as a whole. Although the SPI value was mostly negative during latter decades in all regions (Figure 4c), there was no significant negative trend in precipitation in any region ( Figure A2 in Appendix A). The wettest decade in the eastern region began in 1996, while for the whole province and other regions it began in 1990. The year 1993 was recognized as extremely wet, with the highest value of SPI in the whole province and in the southern region. In the central, northwestern, western, eastern, and northeastern regions, the highest SPI values were in 2005, 2004, 1997, 1999, and 2004, respectively ( Figure 4). Before 1990, the SPI value showed more random fluctuation, but after 1990, two wet (1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006) and dry (2006-2016) periods were identified in all regions (Figure 4c).
The SPI distribution showed that normal conditions, with SPI values between −0.50 and 0.5, were the most common, occurring in 18% of years ( Figure 4b). Based on the mean SPI values for Fars province, normal to moderately wet conditions occurred in 25% of years during the study period. This distribution shows that meteorological drought years were less common in the province than wet and normal years. The SPI distribution showed that dry years were more common in southern and eastern regions of Fars.
Water 2020, 12, x FOR PEER REVIEW 10 of 29 The SPI distribution showed that normal conditions, with SPI values between -0.50 and 0.5, were the most common, occurring in 18% of years ( Figure 4b). Based on the mean SPI values for Fars province, normal to moderately wet conditions occurred in 25% of years during the study period. This distribution shows that meteorological drought years were less common in the province than wet and normal years. The SPI distribution showed that dry years were more common in southern and eastern regions of Fars.

Hydrological Drought
The flow at all 24 gauging stations across Fars province (Figure 1c and Figure A3 in Appendix A) showed a significant negative trend in 1977-2016 ( Figure 5), while before 1997 there was no significant trend for any gauging station ( Figure A4 in Appendix A).
In the period 1977-1996, for more than 70% of stations, 17 out of 24, a positive tendency as with

Hydrological Drought
The flow at all 24 gauging stations across Fars province (Figures 1c and A3 in Appendix A) showed a significant negative trend in 1977-2016 ( Figure 5), while before 1997 there was no significant trend for any gauging station ( Figure A4 in Appendix A).  Table A2 in Appendix A.  Table A2 in Appendix A.

Agricultural Drought
During the study period (1977-2016), irrigated area showed a significant positive trend and rainfed area a significant negative trend, whereas total irrigation showed a non-significant negative tendency (Figure 7a1,7a3). Based on the annual irrigated and rainfed cultivated area, two main  Table A2 in Appendix A.

Agricultural Drought
During the study period , irrigated area showed a significant positive trend and rainfed area a significant negative trend, whereas total irrigation showed a non-significant negative tendency (Figure 7a1,a3) Figure 7a3).
The expansion in cultivated area in the first period  was only in irrigated area, 14,079 ha y −1 , while the rainfed area decreased significantly, -2430 ha y −1 (Figure 7a1,a2). In the second period, both irrigated and rainfed areas showed significant decreases, −20,300 and −7260 ha y −1 , respectively (Figure 7a1,a2).
Based on the OADI values, the irrigated area was below the mean, and the rainfed area was above the mean in 85% of years in 1977-1996, but 80% of years in 1997-2016 (Figure 7b2,b3). The main change in agriculture occurred from 1977 to 2006, when rainfed farmland, low water requirement, was widely replaced by irrigated area, high water requirement (Figure 7b2,b3). The irrigated area reached a maximum 870,000 ha in 2006 and then decreased to 397,000 ha in 2016 (Figure 7a). The greatest rainfed area, 396,000 ha, was in 1983, a normal year according to SPI value (Figure 4), and the least, 60,000 ha, was in 2000, a severe meteorological drought according to SPI value (Figures 4 and 7b). A considerable increase in irrigated area, of around +21,000 ha y −1 , occurred between 1997 and 2006, resulting in an average area of 733,000 ha (Figure 7b3). In the decades 1987-1996 and 2007-2016, the average areas were 580,000 and 610,000 ha, respectively. The contribution of rainfed irrigation to total cultivated area was 41 % during the period 1977-1986, but declined to less than 20% in the period 2007-2016, while irrigated area increased from 59 to 80% in the same period.
The ADI values revealed agricultural drought in Fars province (Figure 7c). Total farmland in most years before 2006 showed positive values of ADI, so no agricultural drought, except in years when the province was affected by extensive metrological drought (Figure 7c3). Generally, rainfall deficit is a major limiting factor for rainfed agriculture. In Fars provenience, in addition to rainfall deficit, the change in agriculture type from rainfed to irrigated farming was another reason for the decreasing ADI values for rainfed farmland in recent years (Figure 7c2). In contrast, the ADI for irrigated area was positive for most years before 2006 (Figure 7c1). This shows that metrological and hydrological droughts did not have any influence on irrigated agriculture and its development.
It can be concluded that the increase in irrigated area was mainly achieved based on groundwater exploitation, as there was a lack of available surface water according to hydrological drought analysis (SDI in Figure 6; trend analyses in Figure A4 in Appendix A). This lowered the groundwater level in the 12 major aquifers in different regions of Fars province (Figure 8). The greatest lowering was recorded in the Arsenjan aquifer in the eastern region of the province during 1995-2015 ( Figure 8). The least groundwater lowering was observed in the Shiraz and Neyriz aquifers, around 10 m over 20 years (Figure 8). The groundwater lowering impact was affected by aquifer area, with small aquifers showing more depletion. For example, the Arsenjan, area 101.6 km 2 , showed the greatest depletion in volume, about 0.212 km 3 , while in Shiraz, area 207 km 2 , the depletion in volume was about 0.123 km 3 . Total, rainfed, and irrigated areas had negative values of ADI after 2006. In that period, it can be concluded that meteorological drought, SPI analysis, hydrological drought, SDI analysis, and the considerable depletion in groundwater resources display evidence of significant agricultural drought (Figure 7c). Based on the OADI values, the irrigated area was below the mean, and the rainfed area was above the mean in 85% of years in 1977-1996, but 80% of years in 1997-2016 (Figure 7b2,b3). The main change in agriculture occurred from 1977 to 2006, when rainfed farmland, low water requirement, was widely replaced by irrigated area, high water requirement (Figure 7b2,b3). The irrigated area reached a maximum 870,000 ha in 2006 and then decreased to 397,000 ha in 2016 (Figure 7a). The greatest rainfed area, 396,000 ha, was in 1983, a normal year according to SPI value (Figure 4), and the least, 60,000 ha, was in 2000, a severe meteorological drought according to SPI value (Figures 4  and 7b). A considerable increase in irrigated area, of around +21,000 ha y −1 , occurred between 1997 and 2006, resulting in an average area of 733,000 ha (Figure 7b3). In the decades 1987-1996 and 2007-2016, the average areas were 580,000 and 610,000 ha, respectively. The contribution of rainfed irrigation to total cultivated area was 41 % during the period 1977-1986, but declined to less than 20% in the period 2007-2016, while irrigated area increased from 59 to 80% in the same period.
The ADI values revealed agricultural drought in Fars province (Figure 7c). Total farmland in most years before 2006 showed positive values of ADI, so no agricultural drought, except in years when the province was affected by extensive metrological drought (Figure 7c3). Generally, rainfall deficit is a major limiting factor for rainfed agriculture. In Fars provenience, in addition to rainfall deficit, the change in agriculture type from rainfed to irrigated farming was another reason for the decreasing ADI values for rainfed farmland in recent years (Figure 7c2). In contrast, the ADI for Although soil moisture is a well-known indicator of agricultural drought [55][56][57], we developed two new indices, OADI and ADI, for agricultural development and drought, based on the change in cultivated area. These two indices complement each other in interpreting drought and changes in irrigated area, as shown in this study. The main advantage these two new indices is that they evaluate the situation based on real agricultural functions in the area. Needed data can be obtained from overall agricultural reports, and in comparison, the soil moisture data are sparsely measured and the available remote sensing data have a low resolution.
two new indices, OADI and ADI, for agricultural development and drought, based on the change in cultivated area. These two indices complement each other in interpreting drought and changes in irrigated area, as shown in this study. The main advantage these two new indices is that they evaluate the situation based on real agricultural functions in the area. Needed data can be obtained from overall agricultural reports, and in comparison, the soil moisture data are sparsely measured and the available remote sensing data have a low resolution.

Discussion
There were clear signs of climate change in Fars province during the study period, with a reduction in the number of different climate types across the province. According to the Emberger index, the land area with a hyper-arid climate increased in the province over the years. This meant that availability of renewable water declined in recent decades, although, based on climate analysis, there was no significant negative trend in precipitation ( Figure A2 in Appendix A). This reduction in water resources could be related to climate variability and/or meteorological drought in the latter part of the study period.
Among the different categories of drought, meteorological drought based on SPI was found to be a non-stationary stochastic phenomenon in trend analysis for the province, while hydrological drought based on flow data clearly showed a negative trend, with significant lowering of groundwater levels. After 1997, hydrological drought conditions dominated, based on SDI values. Irrigated area continued to expand until 2006, when it reached its absolute maximum in the history of the province. The lack of sufficient surface water, the dominant negative phase of SDI ( Figure 6) in main rivers after 1997, could confirm that the expanding by about 400,000 ha of irrigated agriculture between 1997 and 2006 was mainly established based on groundwater exploration. It exerted more pressure on groundwater resources, as seen in groundwater level in major aquifer. A metrological drought period from 2006 onwards led to a further reduction in rainfed area, accelerated hydrological drought, reduced groundwater recharge, and increased groundwater overexploitation. Consequently, the agricultural drought period that started in 2006 resulted in a significant decrease in cultivated area. Overall, before 1997, hydrological drought was a stochastic phenomenon and almost followed the meteorological drought, whereas after 1997 it was mainly due to increasing water use in headwaters and metrological drought.
While drought is viewed locally as a force majeure event in most cases, detailed research has shown that drought is often related to human activities and could be avoided with proper management [58]. In Fars province, overexploitation of water resources, both groundwater and surface water, led to more frequent hydrological droughts after 1997. In this study, hydrological drought was observed over more than 72% of the period 1997-2016, whereas meteorological drought occurred during only 40% of that period. In contrast, in the period 1977-1996, hydrological drought was only observed for 42% of it, and more than 75% of gauging stations showed positive tendency in flow. Thus, the use of water resources and the irrigated area before 1997 were matched by the potential renewable water resources in the province, but after 1997, water consumption exceeded the renewable water available. After 1997, groundwater was the main source of water for irrigation of new cultivated areas and the groundwater resource became depleted. This led to a significant reduction in cultivated area after 2006, due to unavailability of water resources to maintain existing cultivated area in the whole province. The decreasing contribution of rainfed agriculture, from 40% to 20%, was another driver of unsustainability in the use of water resources and of increasing hydrological droughts in the period 1997-2016. Thus the 10 years after 2006 with hydrological drought resulted in agricultural drought, as reflected in decreasing cultivated area. Furthermore, Fars experienced several meteorological droughts in later years, and these undeniably contributed to water shortage. In the past, for example, before 1990, when the irrigated area was smaller and rainfed area made up a greater proportion of total cultivated area, groundwater resources were able to compensate for occasional rainfall deficits and the overall stress on water resources was minor.
The dry period particularly after 2006 and changing from rainfed to irrigated cultivated area led to a tipping point where water demand exceeded the sustainable level for water resources and potential for renewable water. The reduction in river flow and significant depletion in groundwater level confirms this unsustainability in use of water resources. In arid and semi-arid regions, groundwater is a secure source of water that should be saved for drought periods and water shortages. Basing agriculture on groundwater resources has led to agricultural droughts and has also had many socioeconomic and environmental impacts; e.g., land subsidence [59], abandonment in rural areas [60], depression, suicide of rural farmers [61,62], and household food consumption [63]. Overall, we identified three clear stages in Fars province in the 40-year study period.
The second stage (1996-2006) saw a transition from sustainability to unsustainability. It was a golden economic period for farmers, with total cultivated area reaching a maximum and the proportion of irrigated area increasing to 80%, mainly based on groundwater. The new economic conditions increased the expected water demand of local inhabitants beyond a sustainable state. This led to rapid loss of considerable groundwater resources, which altered the interaction between groundwater and surface water, reducing surface water recharge from groundwater [29]. Reduced flow in rivers and increasing frequency of hydrological droughts increased the pressure on the water resources system towards unsustainability.
The third stage from 2006 onwards was unsustainability, with a significant reduction in agricultural area being due to unavailable water resources and the overexploitation of groundwater with deeper wells and greater energy demand for pumping. Continued low flow in the rivers and recurring hydrological droughts, led to exacerbating problems in the operation of installed hydraulic structures, such as dams, on rivers [29]. Insecurity in groundwater resources and unsustainability in agriculture led to increasing internal migration and decreasing cultivated area. There were also environmental impacts, such as desiccation of lakes and wetlands, land subsidence, and desertification [64][65][66].
Overall, we cannot conclude that the whole of Fars province is in stage three, but much of it appears to be. As in other arid and semi-arid regions, the province must seek to move water usage in a sustainable and efficient direction. They must change crop patterns by using low consumption crops; and revise the water allocation policy and prioritize water for consumption. Optimize the reservoir operation and managed aquifers. Drainage control could be effective for irrigated farmland; farmers should select crops that are more suitable for the arid climate and need less water per hectare; and smallholder farms should be supplied with information on seeding and harvesting times for each crop type, for better water usage efficiency. Groundwater must be protected as a vital resource for diminishing the impact of climate variability and climate change.

Conclusions
In this study, we examined the transition from sustainability to unsustainability in water resource use and agriculture in an arid and semi-arid region and the temporal causes and effects of meteorological, hydrological, and agricultural droughts. We developed a new way to evaluate climate change and variability, based on two main elements of the Emberger aridity index, humidity influence factor (HIF) and climate influence factor (CIF), and considering the change in variety of climate types. We also developed two new indices for agricultural drought analysis, overall agricultural drought index (OADI) and agricultural drought index (ADI), which are based on changes in irrigated and rainfed farmed area. We demonstrated the outcomes using the example of Fars province, one of the most important agricultural regions in southern Iran. We showed that the number of climate types has decreased over time in different regions of the province, from a maximum of 20 in the period 1987-1996 to 14 in the period 2007-2016. A hyper-arid climate dominated, with on average 32% spatio-temporal coverage for the 40-year study period and a significant increase from 32.3% in 1977-1986 to 44.4 in 2006-2015. Around 75% of years between 1977 and 2016 were normal and wet in terms of meteorological conditions and uniformly distributed across the whole period, while the frequency of hydrological and agricultural droughts increased particularly in later years (2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015). Regarding water resources' sustainability and resilience, we identified three clear stages during the 40-year period: (1) Sustainable use of water resources before 1997, when water demand was matched by renewable water capacity and irrigated farmland area, was small relative to rainfed area. (2) Transitional use of water resources in 1997-2006, with increasing irrigated area accompanied by hydrological drought, leading to significant groundwater depletion and reduced resilience of the water resources. (3) Unsustainable use of water resources after 2006, with a lack of water resources, particularly groundwater, together with climatological and hydrological droughts, resulting in agricultural droughts and a considerable decrease in farmland area in Fars province.    Table A2.  Table A2.   Table A2.

Appendix B
To calculate the SPI in a specific location, the precipitation data for a long-term period are fitted to a probability distribution and then converted into a normal distribution. Estimating the SPI involves describing the frequency distribution of precipitation, using a gamma probability density function (Equation (A1)): where x is precipitation amount, α is a shape parameter, β is a scale parameter, and the gamma function is expressed as: The maximum likelihood method is used to estimate the optimal values of α and β.
where n is the number of precipitation records and x is average precipitation amount. The cumulative probability for a given month then can be obtained from: Setting t = X β , the incomplete gamma function can be obtained from: It is possible to have several zero values in a sample set. In order to account for zero values, the cumulative probability function for gamma distribution is modified as: where q and (1−q) denote the probability of zero and non-zero precipitation, respectively. The SPI is then derived from the cumulative probability: for 0.5 < H < 1: SPI = +(k − c 0 + c 1 × k + c 2 × k 2 for 0 < H < 0.5: SPI = −(k − c 0 + c 1 × k + c 2 × k 2 where c 0 = 2.515517, c 1 = 0.802853, c 2 = 0.010328, d 1 = 1.432788, d 2 = 0.189269, and d 3 = 0.001308 [67].