Spatial assessment of drought features over different climates and seasons across Iran

Drought is one of the most complex phenomena in the world; so, proper management is very important in monitoring and reducing its damage. For this purpose, Standard Precipitation Index (SPI), Standardized Precipitation Evapotranspiration Index (SPEI), and Reconnaissance Drought Index (RDI) indices were used to analyze the intensity and frequency of drought in the coastal wet, mountain, semi-mountain, semi-desert, desert, and coastal desert climates of Iran in four seasons, separately: autumn, winter, spring, and summer. Forty-three synoptic stations with a common statistical period of 50 years (1969–2019) were selected. The results showed that the trend of drought in winter and summer is increasing in all studied climates. The comparison of the results in the trend analysis of the drought showed the same trend, but the SPEI index compared to the other indicators showed a quicker response to changes in drier climates. The highest correlation (0.80–0.99) between SPI-RDI and SPEI-RDI indices in coastal desert, mountain, and semi-mountain climates and the lowest correlation (0.34) between SPI-SPEI and SPEI-RDI indices in semi-desert, desert, and coastal desert climates were obtained. SPI-RDI variations showed similar values in colder climates. The SPEI is based on precipitation and temperature data, and it has the advantage of combining multi-scalar character with the capacity to include the effects of temperature variability in the drought assessment. Thus, SPEI is recommended as a suitable index for studying and identifying the effect of climate change on drought conditions.


Introduction
Today, many countries are experiencing the negative effects of global warming and climate change (Bachmair et al. 2016;Mazdiyasni and AghaKouchak 2015;Huang et al. 2016;Ficklin et al. 2015). Drought and water scarcity can cause serious crises in human social life, such as the competition over access to water resources, endangering public health, and adverse environmental and ecological impacts (Hosseini-Moghari and Araghinejad 2015; Haile et al. 2019). Therefore, considering the climatic conditions, especially in arid and semi-arid regions, the study of drought status can be an important role in providing forecasting to achieve the stability of systems (Sharafi et al. 2016). Since drought is a natural and recurring phenomenon that occurs due to reduced rainfall over a specified and relatively long period of time (Sharma and Mujumdar 2017), one of the important components of drought planning is the provision of timely and reliable climate information, including seasonal forecasts, which assist decision makers at all levels in managerial decisions. This information, if used properly, is crucial to understanding the effects of drought on various components of the hydrological cycle.
However, with rising temperatures and, consequently, increased reference evapotranspiration (ET ref ), declining precipitation and severe climate change are putting more pressure on water resources and agriculture; therefore, decision makers must have accurate information about the amount of evapotranspiration and drought risk assessment for each region (Byakatonda et al. 2018;Vicente-Serrano et al. 2011;Wang et al. 2012). Assessing drought caused by climate change using long-term meteorological data can be effective in predicting climate change (Guo et al. 2017), but with the complexity of the drought, it is much more difficult to quantify. Lack of direct methods for drought assessment has led to the use of drought indicators that proxy the effects of drought on hydrological systems (Svoboda and Fuchs 2016;Byakatonda et al. 2018;Yihdego et al. 2019). The first step is the selection of climatic parameters affecting the accurate assessment of ET ref and its effect on the accurate assessment of drought indicators. For example, it is not yet clear whether an increase in temperature automatically translates into an increase in ET ref demand, as one theory holds that other meteorological variables such as wind speed, relative humidity, and sunshine hours may not vary with temperature (McEvoy et al., 2012;Vicente-Serrano et al. 2015). However, so far, more than 50 methods of estimating ET ref have been presented in the form of combined methods, aerodynamic and empirical, which often have different results according to meteorological data (Sharafi and Mohammadi Ghaleni 2021a, b). In addition, the accuracy and sensitivity of the ET ref equation is very important, which can even provide important evidence for identifying models that have a more accurate estimate of drought under climate change conditions. For example, complex models such as the Penman-Monteith equation based on FAO 56 (PM-FAO 56 ) is often recommended because of its ability to better represent climatic variables that affect ET ref (Lu et al. 2016). To calculate accurately the PM-FAO 56 equation, required climatic variables include maximum, minimum, and average temperatures, precipitation, relative humidity, wind speed, and solar radiation (Kisi et al. 2015).
Among the various drought indices, the SPI has been studied worldwide (Güner Bacanli 2017;Pathak and Dodamani 2019). However, since SPI does not consider the effect of meteorological variables (maximum and minimum temperature, relative humidity and ET ref ) in the accurate assessment of drought. For this purpose, Vicente-Serrano et al., 2010a, b) proposed SPEI and RDI, respectively (Vicente-Serrano et al. 2015;Mohammed and Scholz 2017). In addition to precipitation, the effect of ET ref is also measured in SPEI and RDI indices (Vicente-Serrano et al. 2010b). Therefore, it is necessary to use various drought indicators to identify climatic regions with similar drought behavior (Dai 2011;Pathak and Dodamani 2019).
Numerous studies have been conducted on the evaluation of SPI, SPEI, and RDI drought indices in different climates (Cloppet 2011;Yihdego et al. 2019;Pathak and Dodamani 2019;Mehr and Vaheddoost 2020). In addition to climatic parameters, other factors such as land use, vegetation, human factors, biodiversity, etc. also have a significant impact (Mohammed and Scholz 2019;Yihdego et al. 2019;Wable et al. 2019;Jehanzaib et al. 2020;Wang et al. 2020;Salimi et al. 2021). Mahmoudi et al. (2019) compared the performance of seven precipitation-based drought indices and found SPI the best index for drought monitoring in Iran. Liu et al. (2017) found that estimating trends in SPI for various timescales is inadequate for accepting the variations in droughts or determining the impact of climate change on the frequency and intensity of droughts. Therefore, it is essential to evaluate the variations in spatial levels of different drought features. Tsakiris et al. (2007) in their study compared SPI and RDI indices in the Mediterranean climate. To evaluate the performance of the RDI index, they used its correlation with the SPI index and stated that RDI has a high behavioral similarity with SPI. Potop and Možný (2011) examined the application of the SPEI index as a new drought index in the Czech Republic and concluded that regardless of the time scale, this index has the ability to detect the severity of drought under the influence of increasing temperature conditions. Their results also showed that the assessment of drought over 10-year periods was due to the increase in ET ref . According to the mentioned cases, it seems that different indicators of meteorological drought monitoring in different climates do not always show the same behavior in monitoring the meteorological drought situation. Also, in two separate studies on the comparative evaluation of SPI with RDI that identified different climate conditions, RDI was recommended in the drought monitoring systems of Iran (Jamshidi et al. 2011;Zarch et al. 2011). Also, Vicente-Serrano et al. (2015 and Salimi et al. (2021) reported that the SPEI, in caparison with SPI, had better capability in representing drought situation.
Iran is located in the arid and semiarid zones of the world. Also, according to the present authors' best knowledge, no studies have focused on examining the ET ref -based drought indices exclusively at the climate and season levels. The current study aims to compare the multivariate meteorological drought indices for Iran's climate by considering the four seasons that consist of six major climatic regions of the country. Therefore, the study includes the following objectives: (1) Investigation of trends associated with the meteorological variables (Tmin. Tmax, precipitation and ET ref ); (2) Assessment of drought characteristics, estimated from SPI, RDI, and SPEI and its spatiotemporal variation at various time scales, and (3) Comparison of these drought indices in different climatic regions and various seasons within the study area.

Time and location scales
Iran, with area 1.65 million km 2 , is located at longitude 25-40°N and latitude 44-65°E and over of the country, involved with desert and semi-desert climates. The highest elevations are found in the Zagros and Alborz mountain ranges, which block moisture from the mountain, semimountain (northwest and west), and coastal wet (north) climates, resulting in uneven distributions of temperature and precipitation across the country (Fig. 1). Therefore, in order to evaluate the hydrological drought indices, 50 years of long-term meteorological data were gathered from 43 synoptic stations . All the important geographical areas of the country were covered by selected stations. Climatic data included values of maximum temperature (T max ), minimum temperature (T min ), average temperature (T mean ), relative humidity, dew point temperature, precipitation, sunshine hours, cloudy days, and solar radiation. Observations that were suspicious or incorrect were removed from the list. The month and day formats were replaced by the day format of the year (DOY) (Sharafi and Karim 2020). The stations studied based on climatic classification of Alijani et al. (2008) in 6 climatic zones included wet coastal, mountain, semi-mountain, semi-desert, desert, and coastal desert. Figure 1 shows the geographical location and climatic classification based on the values of T min , T max , precipitation and ET ref for the four seasons at the studied stations.

Data analysis
First, initial statistical tests were performed on meteorological data. For this purpose, the Mann-Kendall test was used to detect the trend, and Grabs-Back test was used to find data. In the meantime, missing data such as temperature, sunshine, and relative humidity were restored through FAO guidelines, and precipitation and wind speed parameters through average ratios. The PM-FAO 56 empirical equation was used for ET ref calculations (Sharafi and Ghaleni 2021a, b).
The T mean , T min , and T max over Iran and the whole study period 1969-2019 are 18.3, 11.5, and 25.1 °C, respectively. Mean annual precipitation and ET ref in the same T a +273 u 2 e s − e a Δ + 1 + 0.34u 2 In 67% of the country, mean annual precipitation is less than 156 mm (Fig. 1). The spatial differences of 16.65 °C for mean temperature and 1107 mm per year for precipitation show the climate diversity in Iran, as also reported in previous researches (Madani, 2014;Panahi et al. 2020).
Annual average T min varied from − 1.5 to 21.8 (≈7.1 °C) in autumn season, − 6.4 to 18.4 (≈3.1 °C) in winter season, 6.4 to 26.8 (≈15.6 °C) in spring season, and 10.9 to 30 (≈20.2 °C) in summer season. Annual average T max varied from 12.3 to 33.3 (≈20.3 °C) in autumn season 6.3 to 27.2 (≈14.8 °C) in winter season, 21.6 to 44.5 (≈30 °C) in spring season, and 28.1 to 48.7 (≈35.5 °C) in summer season. Furthermore, in desert and coastal desert climates, there is no significant difference between the T min and T max in spring and summer seasons (Fig. 1). These recordings have about the same gradient. The valleys in the mountain are normally colder in winter compared to the foothills. However, summers in these foothills are normally warmer.
Annual precipitation declines from > 1162, 385, and 340 mm (coastal wet, semi-mountain, and mountain climates, respectively) to < 216, 170 and 83 mm (semidesert, coastal desert, and desert climates, respectively). These recordings have about the same gradient in autumn, winter, and spring seasons. Also, annual ET ref increased from > 3469, 3383, and 2327 mm (coastal desert, desert, and semi-desert climates, respectively) to < 1545, 1415, and 1048 mm (semi-mountain, mountain coastal wet climates, respectively). Table 1 shows the statistical values in the form of heating map and significance at the level of 1 and 5%. According to the Mann-Kendall test results, many stations show an increase in T min and T max during the fall, winter, spring, and summer seasons. The slope of the warming trend was much steeper in winter and summer. In coastal wet and coastal desert climates, the trend of increasing T min was observed during the four seasons.

Mann-Kendall test
In general, during the last 50 years, in coastal wet, mountain, semi-mountain, semi-desert, desert, and coastal desert climates, about 70, 65, 80, 77, 80, and 90% of precipitation, respectively, occurs in autumn and winter, respectively. Therefore, the study of this climatic parameter has a very important role in better assessment and understanding of drought indicators. Accordingly, in all studied climates, a trend of reduced precipitation was observed, especially in winter; however, in mountain and semi-mountain climates, this declining trend was more severe. Also, a decrease in precipitation in spring season was observed for stations in coastal wet and mountain climates. On the other hand, in most semi-mountain, desert, and coastal desert climates, the amounts of increase in precipitation were reported in autumn season, although these values were not significant.
According to the results of preliminary studies, ET ref values in most of the stations studied in different climates have an increasing trend, which had an increasing and significant trend in winter and summer. The difference in ET ref values in coastal wet and coastal desert climates is about 2400 mm per year. This amount of ET ref in coastal desert climate is even higher than in desert climate, which is due to the high RH in this climate.

Determinant of change point (CP) using Pettit's test
Time of significant change in time series (changing point) is of particular importance in describing their behavior. In order to detect the CP in time series, several methods are used. In this paper, the non-parametric test proposed by Pettit (1979) was used. This test determines the time of occurrence of significant changes in the mean of the time series. The test statistic is μ t , N , which confirms whether the two sample sets X 1 ,…., X t and X t + 1 …… X N are from the same statistical population, and is defined as Eq.
(2) (Gao et al., 2012): where N; total number of observations, X t ; the first set and X j ; is the second set. The K N test statistics and the associated probability of occurrence are calculated by Eqs. (3) and (4): In the next stage, the drought situation in the monthly period was determined based on SPI, RDI, and SPEI indicators in each station and based on each of the indicators. In this study, the drought situation was examined on an annual scale in 12 consecutive months ending in September of each year, because in terms of agriculture and water resources, each crop year or water year is in this period. On the other hand, on an annual scale, the effect of seasonal fluctuations in precipitation and relative humidity is somewhat combined, so the drought (humidity) in each season was evaluated independently of other seasons in this study, in which this is also a special case on a quarterly scale. In terms of agriculture and water resources, the study of seasonal fluctuations in rainfall and drought is also very important (Sharafi et al. 2016).

SPI index
The SPI index was developed by McKee et al. (1993) to determine periods and severity of drought. SPI calculates precipitation deficit at different time scales (Yihdego et al. 2019). This index is calculated for each climate based on the recording of long-term precipitation.
In this method, first, the appropriate statistical distribution is fitted to the long-term precipitation statistics, then, the cumulative distribution function is converted to a normal distribution using equal probabilities, so that in the period under review, the average standardized precipitation values equal zero and deviate from their standard by one (Edwards 1997). To fit the precipitation Table 1 The results of Mann-Kendall test for the studied stations ;significant at 1% with increasing decreasing trend, Δ; significant at 5% with increasing trend, ; significant at 5% with decreasing trend , ; significant at 1% with trend, ; not significant with increasing trend and o; not significant with decreasing trend τ ≤ -0.2 -0.2 < τ < -0.0.5 -0.05≤ τ ≤ 0.05 0.05 < τ < 0.2 τ ≥ 0.2 data, the gamma distribution is the most accurate distribution. In some climates, however, the gamma distribution may be fitted as the most appropriate distribution for precipitation data (Logan et al. 2010). In fact, SPI is a variable of the standardized normal distribution function whose cumulative probability value is equal to the cumulative probability value of the variable obtained from the obtained gamma distribution (Svoboda and Fuchs 2016).

RDI index
This index was first introduced by Tsakiris and Vangelis (2005). In order to calculate the normal log distribution function on the cumulative values, the ratio of precipitation to ET ref is fitted and then the obtained values are standardized (Tsakiris et al. 2007 where a k is the arithmetic mean of the values in the years under review. To calculate the standardized RDI (RDI st ), it is assumed that the a k values follow the normal log distribution. a k values for different years are calculated using Eq. (7): where Y k is equal to lna k , Y , and ̂ are the arithmetic mean and standard deviation of the values of Y k , respectively. Thus, the RDI is calculated in three steps, which include initial values, normalized values, and standardized values. Standardized values can be compared with SPI (Hatefi et al. 2016;Yihdego et al. 2019). Calculating this index, like the method presented for calculating the SPI index, requires estimating the cumulative probability values of D i values by fitting a probability density function. Since the values of D i lead to negative values from the lower bound, the twoparameter probability functions cannot be selected for this value. Vicente-Serrano et al. (2010a), by examining different three-parameter functions, identified the threeparameter log-logistic probability density function as having the best fit for D i values. The general form of the probability density function of this function is shown in Eq. (9): where the parameters , , and are the scale, shape, and principal parameters for the values of D i in the domain ≥ D ≤ ∞ , respectively. The form of the three-parameter log-logistic cumulative probability function is also calculated according to Eq. (10): Finally, the cumulative probability values for the cumulative probability values for the D i values are converted to a standardized normal function with a mean of zero and a standard deviation of one equal to the SPI values. The classical Abramowitz and Stegun function (1965) is an estimate of the SPEI value using the values of the function F(x) (Eq. (11)): The values of W for different time periods are calculated from Eq. (12): where P is the probability value in excess of the set value D (P = 1-F(x)). If the value of P is greater than 0.5, the P value is replaced by 1-P and the SPEI result symbol is replaced. The value of constant coefficients in Eq. (8) is equal to:

Drought category
SPI concepts have been used in the development of RDI and SPEI indices. Thus, the values of the different classes of drought severity in them are similar to those presented by SPI, as shown in Table 2 (Edwards 1997). According to this method, the drought period occurs when the SPI is continuously negative and reaches a value of − 1 or less and ends when the SPI value becomes positive.

Correlation analysis
Based on the results of the Pearson correlation coefficient, the highest correlation was observed between SPI and RDI during the four seasons: autumn, winter, spring, and summer (R = 0.95-1). On the other hand, the lowest correlation coefficient between SPI and SPEI was observed in desert climates and semi-desert in quadruple climates (R ≤ 0.79). By moving from a wetter climate (such as coastal wet and mountain) to a drier climate (such as desert and coastal desert), the correlation between the SPI and SPEI, SPEI, and RDI indices has been greatly reduced. Due to the direct effect of T min and T max in the evaluation of SPEI and RDI indices, the difference between T min and T max should be evaluated. Accordingly, in humid (such as coastal wet) and very dry climates (such as semi-desert, desert, and coastal desert), the difference between T min and T max in winter (cold season) and summer (warm season) was reported 1 and 2-3° Celsius, respectively. However, the difference between T min and T max between cold and hot seasons in mountain and semi-mountain climates was between 11 to 18 and 13 to 20 °C (Fig. 2). The Pearson correlation matrix shows that the SPI drought index is better than the SPEI drought index in wet climates (Salimi et al. 2021). Also, they reported that wet climates have less ET ref . It means that the drought indices in the humid climates affect the hydrological drought at drier seasons and the maximum correlation is observed for wet seasons. Furthermore, there are seasonal variations in crosscorrelation coefficients, which can be explained by the fact that metrological drought indices in response to the deficit in precipitation may differ significantly, depending on seasonal variations. Additionally, it can be concluded that a stronger relationship begins when precipitation starts after a long, dry summer season in Iran. In this stage, the water moves both through the soil and over the surface via a range of hydrologic processes such as infiltration and runoff throughout the autumn and winter. Simultaneously, it appears that during a dry season (summer) there is a more robust relationship between the SPI and RDI. It can be concluded that the SPEI has a delayed response to the SPI during the dry season, whereas in coastal wet climate, the prompt response of the SPEI and RDI to the SPI is dominant. In this perspective, it is acknowledged that the characteristics of the meteorological drought, too, are significantly dependent on the seasonal variations.

CP analysis
Drought characteristics were extracted from SPI, SPEI, and RDI time series from different climates and seasons. This provides a comprehensive view of meteorological droughts  Table 3 shows the trend of significance, the CP of drought, the first mutation, and the second mutation in six climates and in four seasons. Based on P-value results, SPI, SPEI, and RDI indices were significant in all climates and seasons studied (60% ≤ coverage < 80% among all stations in each climate). This confirms that according to the results of the studied indices, during the past half century, drought was found in all climates and seasons of Iran. Also, according to the results of Table 3, the CP was reported to be decreasing in almost all climates and seasons. For example, the CP in the coastal wet climate began in November 1988. CP has been observed in mountain, semi-mountain, and coastal desert climates since November 1998 and was the same in almost all seasons. This shows that the changes in these climates were due to the human effects of previous years. Factors such as the destruction of forests and pastures (conversion to agricultural fields, especially in mountain and semi-mountain climates) and the uncontrolled construction of dams in this period and in the years before the occurrence of successive droughts have had a direct impact on the occurrence and registration of droughts. Furthermore, based on the research of Tabari et al. (2013), almost all climates experienced extreme droughts in Iran.  1998-1999, 1999-2000, and 2000-2001 were the driest years during the studied period.

Drought severity
For each season in different climates, drought severity and the total number of events were identified using thresholds of − 1 to − 1.5, − 1.5 to − 2 and more than − 2 (moderate, severe and extreme drought, respectively). For SPI and SPEI, more drought severity and the total number of events were identified at thresholds of more than − 2. Spatial patterns for the maximum severity characteristic were similar for all three thresholds, and the results not be shown. The results for the thresholds of − 1 to − 1.5, − 1.5 to − 2 and more than − 2 for different climates and seasons can be found in Figs. 3, 4, 5, and 6. In general, the results showed that, for drought severity more than − 2, there was comparatively little difference among different seasons (as shown by the clusters in Fig. 3), or around the country in the meteorological drought frequency extracted from SPI, SPEI, and RDI time series (Fig. 4). While Iran has an order of magnitude precipitation gradient across the country, there is little difference in the severity and frequency of the meteorological drought indices. For SPI with increasing in drought severity, there is little variation between the four seasons' catchments, especially for the threshold of more than − 2 in different climates (Fig. 3). This confirms that SPI index varies only modestly across the country over the threshold of more than − 2 once the precipitation has been standardized. However, in SPEI and RDI index, these changes were much greater at low drought intensities (− 1 to − 1.5) and with increasing drought severity, the changes between seasons, especially at thresholds greater than − 2, were more significant. In the SPEI index, changes in drought severity were greater in winter and spring for coastal wet, mountain and coastal desert climates. However, drought severity changes in autumn, spring, and summer in semi-mountain and semi-desert climates were greater than in winter. Also, according to the results of the RDI index, the highest drought severity was observed in the thresholds higher than − 2 in the four seasons. However, most of the changes in drought severity in  Class Autumn Winter Spring Summer -1 to -1.5 -1.5 to -2 <= -2 winter were observed in coastal wet and mountain climates. Similarly, Barker et al. (2016) found little spatial variation in the number and severity of meteorological events among catchments in the UK. Figure 4 shows the distribution of drought intensities at the threshold of − 1 for the surveyed indices in the four seasons. Table 1 shows that Iran has an order of negative precipitation gradient across the country (especially in winter season). Therefore, there is little difference in the median of drought severity. However, there may be differences in the results of drought indices in stations with similar climates. These differences in drought severity have also been found elsewhere, e.g. in coastal wet or coastal desert, where spatial variation was found to be the result of climatic variation across the mentioned climates. This result is confirmed by Panahi et al. (2020), who noted that the coastal wet and mountain climates of Iran have more variability in drought severity and that the desert and coastal desert climates are subject to longer dry spells, so in practice, these zones experience droughts in opposition.

Number of events
Compared to the SPI, the RDI shows fewer total events for both thresholds of − 1.5 to − 2 and more than − 2, which can be attributed to the stronger persistence (i.e., a smaller fluctuation) that is more likely to be characterized by the SPI (Farahmand and AghaKouchak 2015;Kwon et al. 2019). Furthermore, more frequent drought events appear to occur in thresholds of more than − 2 for the SPI, SPEI, and RDI indices (Fig. 5). It was also clearly seen that the total number of events for the SPEI is significantly longer than that of the SPI and RDI (especially in − 1 to − 1.5 threshold). On the other hand, SPI and RDI display negative values more frequently than does SPEI, and SPI recovers from wet states more quickly than does SPEI. More total events with the threshold − 1 to − 1.5 particularly stand out in coastal wet, mountain and semi-mountain climates for SPI, while an increased frequency of moderate droughts is identified in all climates for SPEI and RDI indices.
In general, significant severe droughts in terms of both the total number of events and drought severity primarily occurred in all climates of Iran especially in wet seasons (autumn and winter). Also, the spatial distribution of drought characteristics associated with SPI is more dependent on the threshold, while SPEI and RDI are less sensitive to the threshold (Fig. 5). Figure 6 shows the total number of events detected by each of the drought indices in different climates and seasons in Iran. Accordingly, the highest number of events recorded by the SPI index was more than all other drought indicators in all climates and seasons. The number of drought events reported by the SPI index in all climates and seasons was greater than 30. This occurred while the highest record was detected by the SPEI in the autumn and winter seasons in coastal desert climate (e.g. Bandar Abbas station with 40 seasons).
In general, all three indices confirmed that in all climatic classes, the highest total number of recorded events occurred in the autumn season. However, on average, the number of drought events in the autumn for the SPI index was 19 and 32% higher than the SPEI and RDI indices, respectively. The correlation analysis between different drought indices indicates that the minimum and maximum correlation was in desert and coastal wet climates, respectively. So, the average correlations were obtained in desert (0.74) and coastal wet (0.96) climates. The relative frequency of different drought classes in SPEI is more than RDI and SPI, especially in desert and coastal desert climates. Results of Pettit's test demonstrated significant change points in all three indices SPI, SPEI, and RDI. The change points occurred mainly from 1998 to 2009 and the hydrological years of 1998-1999, 1999-2000, and 2000-2001 were the driest years during the studied period. In general, the main change point of drought occurred in the 1998-1999 period, so that after this year, the longest duration and the most intensity appeared in 1998-2002 periods over Iran.
Also, the total number of drought events recorded in desert and coastal desert climates did not differ significantly in any of the seasons. However, these conditions were different for semi-desert, semi-mountain, mountain and coastal wet climates, so the total number of detected events in autumn and summer were more than in other seasons. This issue confirmed that the role of declining precipitation gradient in autumn is very disturbing (Fig. 6).

Conclusion
The phenomenon of climate change and its effects is one of the major challenges of water resources management in this century. It is important to determine how drought conditions change based on the diversity of different seasons and climates from various aspects, because many long-term plans require a drought forecasting system to monitor the future status of precipitation, drought and wet periods for different climates. Based on the results, in all studied climates, there is a high correlation SPI-RDI index pair (R = 0.99). The lack of difference between SPI and RDI can be considered as not determining the ET ref in the occurrence of drought. As a result, the SPI Fig. 3 Boxplots showing meteorological drought severity based on SPI, SPEI, and RDI using thresholds of − 1 to − 1.5, − 1.5 to − 2 and ≤ − 2 for each season ◂ index can monitor drought alone by considering precipitation. In addition, in mountain, semi-mountain, semidesert, desert and coastal desert climates, the correlation between SPI-SPEI and SPI-RDI varied (0.40 to 0.89), in which case the SPEI index is a suitable measure for these climates, because considering the two criteria of precipitation and ET ref , the fitting of values is done with the help of a three-parameter distribution (three-parameter Class Autumn Winter Spring Summer -1.5 to -1 -2 to -1.5 <= -2

Fig. 5
Boxplots showing the total number of events based on SPI, SPEI, and RDI using thresholds of − 1 to − 1.5, − 1.5 to − 2 and ≤ − 2 for each season logistics) and as a result will have better control over the form of distribution. The comparison of the results of all three indices in analyzing the trend of the number of drought events at different thresholds and in different seasons showed the trend of changes of three indices in the same zone, but the SPEI index compared to the other two indices in semi-desert, desert, and coastal desert climates reveals a higher frequency of drought and shows a faster response to moisture changes. Spatial analysis of drought trends using Pettit's test Autumn Winter Spring Summer Fig. 6 Maps showing the total number of events based on SPI, SPEI, and RDI using a threshold of − 1 for each season of Iran showed that drought trends are decreasing in most mountain and semi-mountain climates and increasing in semi-desert, desert, and coastal desert climates. This severity was shown by the SPEI drought index first and then by the SPI and RDI indices with more severity.
The results of this study showed that the SPI index is a very simple index and does not require complex data, but, due to its high sensitivity to decrease and increase of precipitation, it can be said that it is an indicator that monitors more accurately. Therefore, it is suggested that after examining the efficiency of indices in the climates and in fact its localization, an important phenomenon such as drought must be investigated. Drought's increase or decrease has economic, social, and political consequences. It is important to properly manage water resources in different climates and in different seasons, especially for food security that depends on climate diversity and water resources. It is generally suggested that in semi-desert, desert, and coastal desert climates, the use of SPI is insufficient to monitor drought, and due to the role of temperature as a limiting factor for humidity and high ET ref , the SPEI index should be used for effective and accurate monitoring of drought. SPEI: Standardized Precipitation Evapotranspiration Index; SPI: Standard Precipitation Index; Spr.: Spring; Sum.: Summer; T a : Average air temperature (°C); T max : Maximum temperature; T mean : Average temperature; T min : Minimum temperature; u 2 : Average wind speed at 2 m height (m s −1 ); VW: Severely wet; Win.: Winter; γ: Psychometric constant (kPa °C ); Δ: The slope of the saturation vapor pressure function (kPa °C ); μ 1 : Mutation 1; μ 2 : Mutation 2 Author contribution S.S. and M.M.G. prepared data and performed model runs, designed the study, interpreted the results, and wrote the manuscript.
Data availability All data used in this paper for Iran were extracted from metrological stations' data and can also be requested from the corresponding author.