Monthly, seasonal and annual trends of rainfall and temperature (both minimum and maximum) have been analyzed using the Mann–Kendall trend test (a non-parametric test) and Sen's slope estimator for Sagar division, India from 1988 to 2018. Sagar division is a drought-prone zone of Madhya Pradesh, India. The same analysis has been performed for two drought indices, the Standardized Precipitation Index (SPI) and Reconnaissance Drought Index (RDI). Both indices were calculated to see the trend in the drought for 35 rain-gauge stations belonging to the study area. The study revealed that the minimum temperature had increased more than the maximum temperature in the last 31 years. The strong similarity in the results of Sen's slope of SPI and RDI were seen for both significant and non-significant trends. Analysis of variance (ANOVA) testing validates the substantial similarity between SPI and RDI based on Sen's slope. It also indicated the suitability of RDI for future projection of drought using the general circulation models (GCMs) or regional climate models (RCMs) in meteorological drought as well as the agricultural drought category. In contrast, the SPI indicated the meteorological drought only. The distribution of trends of temperature and drought indices were presented using the kriging interpolation.

  • The non-parametric trend of drought: One rainfall based drought index, SPI and another rainfall and temperature-based index, RDI are used to see the behavior of drought in the study area, which is drought-prone.

  • The non-parametric trend of rainfall and temperature (maximum as well as minimum) to see the effect of both.

  • Spatial variation of drought as well as meteorological variables using kriging interpolation on ArcGIS: To see the variation by plotting Mann-Kendall Z and Sen's slope on ArcGIS.

  • One way ANOVA approach to validate the applicability of SPI and RDI indices for future drought projection.

  • The study compares both the indices based on Sen's slope values using the ANOVA approach and spatial representation by kriging interpolation.

Bundelkhand has been historically known as a drought-prone region of India, but the frequent and intense droughts have increased in recent decades. The scarcity of water in the semi-arid region, with poor soil and low productivity, further aggravates the problem of food security (Gupta et al. 2014). Hence, further study is expected to be alert to such drought events happening in the future. Analysis of trends is a primary tool to predict future drought scenarios. Yue et al. (2002) studied two non-parametric ranked based tests, namely the Mann–Kendall test and Spearman's rho test. Both tests were compared based on the power of the test by Monte Carlo simulation and gave very similar results of the trend. It is noted that the distribution of Mann–Kendall approaches more quickly than Spearman's rho. These tests become more powerful with an increase in the sample size and decrease when the number of variation grows. The Mann–Kendall trend test and Sen's slope estimator are helpful to find trends effectively for data of more than 30 years. The increasing and decreasing trend shows climate trend and helps to foresee future climate variables (Birara et al. 2018). The main advantage of the Mann–Kendall trend test is that it can be performed at any time scale. It is robust against distribution-free and has low sensitivity to sudden breaks in time series.

Jain & Kumar (2012) performed the Mann–Kendall trend test for the whole of India for rainfall and temperature. The magnitude of the trend in a time series is determined using Sen's estimator method. The test represents both positive and negative trends in the area, although they are not very significant. The test can be performed to determine the pattern for drought indices also. Two relevant drought indices; RDI (Reconnaissance Drought Index) and SPI (Standardized Precipitation Index) were used for 3, 6, 9, 12, 18 and 24-month time scales in 40 meteorological stations in Iran and it was seen that the correlation coefficient (R) decreases with an increase in the time scales between SPI and RDI and can be easily differentiated in the wetter climate for all time scales (Zarch et al. 2011). For the semi-arid region, SPI and RDI can both be performed for trend detection. This study is expected to organize region-specific drought monitoring by identifying the trend of precipitation and temperature along with the significant trend in drought. The kriging interpolation has been adopted to plot the change of temperatures (maximum and minimum), SPI and RDI on the ArcGIS tool. Kriging is one of the geospatial interpolation techniques based on statistics and it is a more advanced tool (Childs 2004). While plotting the Sen's slope values over the study area on the ArcGIS tool, there are many raster interpolation methods available such as inverse distance weightage (IDW), kriging and natural neighbor interpolation, etc. Of all the processes, kriging shows a better representation of spatial figures. The one way ANOVA method is chosen to discriminate both SPI and RDI indices. This method can compare more than two groups. ANOVA uses the F statistic to test if all groups have the same mean (Park 2009; Skidmore & Thompson 2013). It gives better comparative results than other methods such as the t-test.

Sagar division comes under the Bundelkhand region in Madhya Pradesh, India. It consists of five districts (Sagar, Chhatarpur, Tikamgarh, Panna and Damoh). It is situated at 23°49′48″ N latitude and 78°42′36″ E longitude where the nearest city is Sagar, which is the headquarters of the division. The area lies in the heart of India, found below the Indo-Gangetic Plain to the north, with the undulating Vindhyan mountain range spreads from the northwest to the south (Gupta et al. 2014). The location of the Sagar division can be seen in Figure 1.

Figure 1

Location of Sagar division in Madhya Pradesh, India.

Figure 1

Location of Sagar division in Madhya Pradesh, India.

Close modal

Daily rainfall data of 35 rain gauge stations in the Sagar division were collected from the Department of Water Resources (BODHI), Bhopal (M.P.), India. Daily gridded data (0.5 × 0.5° resolution) of maximum and minimum temperature for 31 years (1988–2018) were collected from NASA (https://power.larc.nasa.gov/data-access-viewer/) for 35 stations in the Sagar division. The description of rainfall data along with the rain gauge stations in the Sagar division is shown in Table 1.

Table 1

Rain gauge stations in Sagar division and rainfall data availability

S. no.DistrictRain gauge stations
1. Sagar Sagar (1977–2018), Shahgarh (1993–2018), Jaisinagar (1988–2018), Kesli (1992–2018), Malthon (1988–2018), Rehli (1977–2018), Rahatgarh (1988–2018), Bina (1988–2018), Khurai(1977–2018), Banda (1977–2018), Deori (1977–2018) and Garhakota (1993–2018) 
2. Chhatarpur Chhatarpur (1988–2018), Bijawar (1988–2018), Buxwaha (1988–2018), Laundi (1988–2018), Nowgaon (1988–2018), Rajnagar (1994–2018), Gaurihar (1994–2018) and Badamalhera (1988–2018) 
3. Tikamgarh Tikamgarh (1988–2018), Baldevgarh (1988–2018), Jatara (1988–2018), Palera (1991–2018), Niwari (1988–2018), Prithvipur (1988–2018) and Orcha (1993–2018) 
4. Panna Panna (1988–2018), Pawai (1988–2018), Shahnagar (1988–2018), Gonour (1988–2018) and Ajaigarh (1988–2018) 
5. Damoh Damoh (1998–2018), Hatta (1998–2018), Jabera (1998–2018) 
S. no.DistrictRain gauge stations
1. Sagar Sagar (1977–2018), Shahgarh (1993–2018), Jaisinagar (1988–2018), Kesli (1992–2018), Malthon (1988–2018), Rehli (1977–2018), Rahatgarh (1988–2018), Bina (1988–2018), Khurai(1977–2018), Banda (1977–2018), Deori (1977–2018) and Garhakota (1993–2018) 
2. Chhatarpur Chhatarpur (1988–2018), Bijawar (1988–2018), Buxwaha (1988–2018), Laundi (1988–2018), Nowgaon (1988–2018), Rajnagar (1994–2018), Gaurihar (1994–2018) and Badamalhera (1988–2018) 
3. Tikamgarh Tikamgarh (1988–2018), Baldevgarh (1988–2018), Jatara (1988–2018), Palera (1991–2018), Niwari (1988–2018), Prithvipur (1988–2018) and Orcha (1993–2018) 
4. Panna Panna (1988–2018), Pawai (1988–2018), Shahnagar (1988–2018), Gonour (1988–2018) and Ajaigarh (1988–2018) 
5. Damoh Damoh (1998–2018), Hatta (1998–2018), Jabera (1998–2018) 

Statistical test for trend and variability analysis

Trend analysis is the prediction of future outcomes by using historical results based on their trend. The non-parametric Mann–Kendall test is used to obtain patterns using historical data and its magnitude is evaluated using Sen's slope estimator (Mondal et al. 2012; Rahmat et al. 2015). The increase or decrease in trend is checked, which is based on the normalized test statistics (Z) value. When Z is positive, the trend is said to be increasing and when Z is negative, it is said to be decreasing. The trend's slope gives the annual rate and direction of change.

Mann–Kendall test (non-parametric)

The Mann–Kendall Statistic (S) is given by:
(1)
Here, i = 2,3,….,n; j= 1,2,….,i–1 and:
(2)
For a sample size >10, a normal approximation to the Mann–Kendall test may be used. For this, the variance of S is obtained as:
(3)
where tp is the number of ties for the pth value and q is the number of tied values.
Then the standardized statistical test is computed by:
(4)

Theil–Sen's slope estimator

The magnitude of the trend is estimated by Sen's slope method which is conducted by calculating the slope (β) as a change in measurement per change in time:
(5)
where Q’ is the slope between data points Xt’ and Xt, Xt’ is the data measurement at time t’ and Xt is the data measurement at time t.
The median slope gives Sen's slope estimator as:
(6)
where N is the number of calculated slopes. A positive value of β indicates an increase in the trend, and a negative value indicates a decrease in trend in the time series.

Assessment of drought indices

Model description

The DrinC (Drought Indices Calculator) model is used to evaluate two drought indices, SPI and RDI (Tigkas et al. 2015). The model was developed at the National Technical University of Athens for providing an adaptable and straightforward interface to calculate drought indices. The model is applicable for several locations, especially for arid and semi-arid regions. Sagar division comes under the semi-arid region. SPI is calculated only by rainfall data (McKee et al. 1993), while RDI uses both rainfall and temperature data. It makes it more effective as two meteorological elements are involved there (Tsakiris & Vangelis 2005; Tigkas et al. 2015). SPI is based on the probability of precipitation at any time scale (Tian et al. 2018), while RDI is based on both cumulative rainfall (P), and potential evapotranspiration (PET), in which P is measured and PET is calculated using the Hargreaves method (Tsakiris et al. 2007).

Standardized Precipitation Index (SPI)

SPI, based on probability of precipitation for any time scale, is computed as:
(7)
where X = precipitation for the station; Xm = mean precipitation; = standardized deviation.

Reconnaissance Drought Index (RDI)

RDI is based on both cumulative precipitation (P) and potential evapotranspiration (PET), in which P is measured and PET is calculated. The initial value (αk) of RDI is calculated for the ith year of k (months) as follows:
(8)
in which Pij and PETij are the precipitation and potential evapotranspiration of the jth month of the ith year and N is the total number of years of available data. The values of αk follow both the log-normal and the gamma distributions satisfactorily. By assuming that the log-normal distribution is applied, RDIst is calculated using Equation (9):
(9)
in which y(i) is the ln (k(i)), is its arithmetic mean and is its standard deviation.

Analysis of variance (ANOVA)

ANOVA is a statistical technique that is used to check if the means of two or more groups are remarkably different from each other (Skidmore & Thompson 2013). Statistician and evolutionary biologist Ronald Fisher developed it. ANOVA verifies the effect of one or more factors by comparing the means of different samples at a significance level α = 0.05.

The null hypothesis in ANOVA is valid when all the sample means are equal, or they do not have any significant difference. ANOVA follows F-distribution which does not have any negative values because between and within-group, variability is always positive due to each deviation's square.

The non-parametric Mann–Kendall trend tests of rainfall and temperature were analyzed on monthly, seasonal and annual time series. If the normalized test statistics Z > 1.96 or Z < –1.96, null hypothesis (HO) is rejected at the 95% significance level.

Rainfall trend

Monthly, annual and seasonal non-parametric trends of rainfall evaluation for 35 rain gauge stations (Table 1) of Sagar division are mentioned in Table 1. Out of five districts, Sagar has the maximum number of stations with a significant increasing and decreasing trend (nine stations). Sagar station showed a negative rainfall trend in August (β = –3.783 mm/year) and in the spring season (β = –1.053 mm/year). Jaisinagar and Rahatgarh stations showed a higher decreasing trend in August at β = –7.263 mm/year and β = –10.62 mm/year, respectively, as shown in Figure 2.

Figure 2

Monthly, seasonal and annual rainfall trend of stations showing significant trends of Sagar district.

Figure 2

Monthly, seasonal and annual rainfall trend of stations showing significant trends of Sagar district.

Close modal

Baldevgarh station in Tikamgarh district indicated a higher negative trend in rainfall in August (β = –6.379 mm/year), annual (β = –11.737 mm/year) and also in the rainy season (β = –13.044 mm/year). That means a significant deficiency is seen in future rainfall in the month of August, annually and in the rainy season. Similarly, Nowgaon and Gaurihar stations saw a massive decreasing trend in August.

The less increasing rainfall rate has been seen in most of the rain gauge stations (<1 mm/year). Deori station was the only station that showed β as 4.75 mm/year in July.

It is noticed that the decreasing trend is more dominant compared to the increasing trend of rainfall for most of the stations of the Sagar division, as shown in Figure 3.

Figure 3

Monthly, seasonal and annual rainfall trend of stations showing significant trends of Tikamgarh, Chhatarpur and Damoh districts.

Figure 3

Monthly, seasonal and annual rainfall trend of stations showing significant trends of Tikamgarh, Chhatarpur and Damoh districts.

Close modal

Temperature trend

Temperature is one of the prime keys responsible for climate change. Daily maximum and minimum temperature data of five districts of Sagar division (Sagar, Chhatarpur, Tikamgarh, Panna, and Damoh) over 31 years (1988–2018) were collected and arranged in a monthly, seasonal and annual basis. It is evident that temperature has been increasing day by day which means that the increasing trend in the results was expected. The non-parametric test displayed no significant trend for maximum temperature but a significant increasing trend for minimum temperature for most of the stations on the monthly time scale.

April became warmer in terms of the lower temperature range for almost all of the districts except Panna station. Sagar district mostly suffered due to the increasing minimum temperature in April (β = 0.045 °C/year) and May (β = 0.03 °C/year) and also in summer (β = 0.031 °C/year) and spring (β = 0.016 °C) seasons. Panna district showed a negative trend in the month of May (β = 0.03 °C/year) and annually (β = 0.012 °C/year), as shown in Figure 4. The results of the significant trend of temperature at seasonal and annual scales are presented in Figure 5.

Figure 4

Stations showing a significant trend of temperature at monthly, seasonal and annual scales.

Figure 4

Stations showing a significant trend of temperature at monthly, seasonal and annual scales.

Close modal
Figure 5

Seasonal and annual minimum temperature trend in Sagar Division, India (Mann Kendall–Z).

Figure 5

Seasonal and annual minimum temperature trend in Sagar Division, India (Mann Kendall–Z).

Close modal

The results revealed that the minimum temperature had increased more as compared to the maximum temperature over the last 31 years. Increasing temperature affects the capability of the plant to produce grain efficiently which means a rise in minimum temperature has a bigger effect on crop yield than a maximum temperature rise. Hence, it becomes necessary to carry out a broad range study of temperature behavior and its impact, which will be seen in the future climate by predicting them.

Drought indices trend

Monthly and annual trends of SPI and RDI indices performed for 35 stations were both evaluated using the DrinC model (Tigkas et al. 2015). The results of the significant trend of SPI and RDI are presented in Table 2. Sen's slope of rainy months JJAS (June, July, August and September) over the study area were plotted using kriging interpolation on ArcGIS, as shown in Figure 6.

Table 2

Stations showing significant trends either of SPI or RDI or of both

DistrictStationTime scales showing significant resultsSPI
RDI
Z-ValueSen's slope (β)Z-ValueSen's slope (β)
Tikamgarh Baldevgarh Annual –1.96 –0.046 –1.75 –0.041 
Palera August –0.97 –0.063 –2.06 –0.061 
Orcha August –1.94 –0.064 –1.99 –0.065 
Prithvipur Annual –1.98 –0.031 –1.59 –0.025 
Sagar Jaisinagar August –2.62 –0.056 –2.34 –0.052 
Rehli August –2.08 –0.069 –2.19 –0.068 
September –2.05 –0.041 –2.09 –0.041 
Annual –2.14 –0.035 –1.87 –0.035 
Rahatgarh August –1.89 –0.065 –2.68 –0.059 
Bina August –2.12 –0.044 –2.28 –0.041 
Khurai August –2.71 –0.06 –2.43 –0.052 
Banda August –1.63 –0.062 –2.5 –0.06 
Chhatarpur Nowgaon August –2.11 –0.041 –1.8 –0.044 
Baxwaha August –2.28 –0.048 –2 –0.045 
Gaurihar August –2.21 –0.077 –2.13 –0.08 
Laundi August –2.07 –0.047 –1.89 –0.046 
Panna No significant trend 
Damoh No significant trend 
DistrictStationTime scales showing significant resultsSPI
RDI
Z-ValueSen's slope (β)Z-ValueSen's slope (β)
Tikamgarh Baldevgarh Annual –1.96 –0.046 –1.75 –0.041 
Palera August –0.97 –0.063 –2.06 –0.061 
Orcha August –1.94 –0.064 –1.99 –0.065 
Prithvipur Annual –1.98 –0.031 –1.59 –0.025 
Sagar Jaisinagar August –2.62 –0.056 –2.34 –0.052 
Rehli August –2.08 –0.069 –2.19 –0.068 
September –2.05 –0.041 –2.09 –0.041 
Annual –2.14 –0.035 –1.87 –0.035 
Rahatgarh August –1.89 –0.065 –2.68 –0.059 
Bina August –2.12 –0.044 –2.28 –0.041 
Khurai August –2.71 –0.06 –2.43 –0.052 
Banda August –1.63 –0.062 –2.5 –0.06 
Chhatarpur Nowgaon August –2.11 –0.041 –1.8 –0.044 
Baxwaha August –2.28 –0.048 –2 –0.045 
Gaurihar August –2.21 –0.077 –2.13 –0.08 
Laundi August –2.07 –0.047 –1.89 –0.046 
Panna No significant trend 
Damoh No significant trend 
Figure 6

Sen's slope of rainy months JJAS (June, July, August and September) for (a) SPI and (b) RDI.

Figure 6

Sen's slope of rainy months JJAS (June, July, August and September) for (a) SPI and (b) RDI.

Close modal

The four stations of Tikamgarh had a significant trend (out of seven) for either SPI or RDI or both. Baldevgarh station presented an annual decreasing trend (β = –0.046) of SPI, while Prithvipur also had an annual negative trend (β = –0.031) of SPI. For Palera and Orcha, both stations had a diminishing trend in RDI for the month of August.

Six stations in Sagar (out of 12) displayed a significant trend of either SPI or RDI or for both. Jaisinagar showed a decreasing trend for both SPI (β = –0.056) and RDI (β = –0.052) in August while Khurai and Bina stations showed a negative trend of both SPI and RDI in August. Rehli was the only station which showed a significant trend of both SPI and RDI in two months (August and September) simultaneously. Four stations in Chhatarpur (out of eight) had a significant trend of either SPI or RDI or both. Baxwaha and Gaurihar stations showed a negative trend for both SPI and RDI in August.

It is noticed that no significant trend is seen for any of the stations of Damoh and Panna for both SPI and RDI. There are many similarities found in the results of Sen's slope of SPI and RDI, as shown in Table 2. SPI depends only on rainfall while RDI depends on rainfall as well as temperature. Hence, RDI can give better results for the trend with two variables.

The performance of RDI with two meteorological variables (rainfall and temperature) was found to be quite good, along with the SPI with one meteorological variable (rainfall). Hence, RDI can replace SPI if better correlation is shown between both of them. Sen's slope of both indices showed nearly similar results on an 85% significance level while Z < 1.44 or Z > –1.44. The ANOVA test was performed to see the level of discrimination of both SPI and RDI.

One way ANOVA

For all the stations, values of SPI and RDI are nearly equal, but the discrimination of SPI and RDI is necessary based on their increasing and decreasing behavior. Thus, Sen's slope values were taking for the ANOVA test on an 85% significance level. There are two groups with a single variable, and that is why one way ANOVA is applied, shown in Table 3.

Table 3

Variance results

GroupsCountSumAverageVariance
β of SPI 15 –0.786 –0.0524 0.0001798 
β of RDI 15 –0.754 –0.0503 0.0002005 
GroupsCountSumAverageVariance
β of SPI 15 –0.786 –0.0524 0.0001798 
β of RDI 15 –0.754 –0.0503 0.0002005 

There was a deficient level of discrimination found between Sen's slope of SPI and RDI (Hypothesis is accepted, F ≪ Fcrit), as shown in Table 4. That means RDI shows substantial accuracy towards projections of future drought with two climate variables.

Table 4

ANOVA results

Source of variationSum of squares (SS)Degrees of freedom (df)Mean square (MS)FF crit
Between groups 0.00003 0.000034 0.179 4.196 
Within groups 0.00532 28 0.000190   
Total 0.00535 29    
Source of variationSum of squares (SS)Degrees of freedom (df)Mean square (MS)FF crit
Between groups 0.00003 0.000034 0.179 4.196 
Within groups 0.00532 28 0.000190   
Total 0.00535 29    

It is concluded that the decreasing trend of rainfall is more dominant as compared to the increasing trend for most of the stations of the Sagar division, especially Baldevgarh station, which showed the higher negative trend of rainfall in August (β=–6.379 mm/year), annual (β=–11.737 mm/year) and also in the rainy season (β = –13.044 mm/year). That means significant deficiency was seen in future rainfall in August, annually and in the rainy season. Similarly, Jaisinagar, Rahatgarh, Nowgaon and Gaurihar stations showed a massive decreasing rainfall trend in August.

The minimum temperature had increased more as compared to the maximum temperature during the last 31 years of the study region. No significant trend was seen for maximum temperature in the study region, but for minimum temperature an increasing trend was observed for all stations, in particular the month of April gets warmer in the lower temperature range for almost all of the districts. Sagar district mostly suffered in April (β = 0.045 °C/year) and May (β = 0.03 °C/year) as well as in summer (β = 0.031 °C/year) and spring (β = 0.016 °C/year) seasons.

For most of the stations of the Sagar division, the trend of drought was found to be significant, especially in August. It was also noticed that no significant trend was seen for any of the stations of Damoh and Panna districts. There are many similarities found in SPI, and RDI values at the level of severity of drought and also similar results are seen in Sen's slope of SPI and RDI.

Sen's slope values of SPI and RDI for August presented a very similar magnitude. The significant trend was mostly seen in this month. Figure 6 represented the Sen's slope of rainy months JJAS (June, July, August and September) over the study area.

The ANOVA result of SPI and RDI also displayed substantial similarity to Sen's slope. It means RDI can give better results of the trend with two meteorological variables along with SPI. The study also helps to choose a suitable drought index to project future drought using the general circulation models (GCMs) or regional climate models (RCMs) in meteorological drought as well as the agricultural drought category. At the same time, the SPI indicated the meteorological drought only.

The focus of the current study was to check the applicability of drought indices for the Sagar division so that further drought could be forecast using specific climate models such as GCMs and RCMs. Specific bias correction methods could be applied such as linear scaling (LS), general quantile mapping (GEQM), gamma quantile mapping (GAQM), and power transmission (PT) to downscale the model data (Homsi et al. 2020). Drought indices could also be modeled by the support vector regression (SVR), SVR model coupled with firefly algorithm (SVR-FA), Gene Expression Programming (GEP), and M5 model trees (MT) to forecast drought (Moazenzadeh et al. 2018; Shamshirband et al. 2020).

All relevant data are included in the paper or its Supplementary Information.

Birara
H.
Pandey
R. P.
Mishra
S. K.
2018
Trend and variability analysis of rainfall and temperature in the Tana basin region, Ethiopia
.
J. Water Clim. Change
9
(
3
),
555
569
.
doi:10.2166/wcc.2018.080
.
Childs
C.
2004
Interpolating surfaces in ArcGIS spatial analyst
.
ArcUser
3235
,
569
.
Gupta
A. K.
Nair
S. S.
Ghosh
O.
Singh
A.
Dey
S.
2014
Bundelkhand Drought: Retrospective Analysis and Way Ahead
.
National Institute of Disaster Management
,
New Delhi
, pp.
1
148
.
Homsi
R.
Shiru
M. S.
Shahid
S.
Ismail
T.
Harun
S. B.
Al-Ansari
N.
Chau
K. W.
Yaseen
Z. M.
2020
Precipitation projection using a CMIP5 GCM ensemble model: a regional investigation of Syria
.
Eng. Appl. Comput. Fluid Mech.
14
(
1
),
90
106
.
Jain
S. K.
Kumar
V.
2012
Trend analysis of rainfall and temperature data for India
.
Curr. Sci.
102
(
1
),
37
49
.
McKee
T. B.
Doesken
N. J.
Kleist
J.
1993
The relationship of drought frequency and duration to time scales
.
Proc. 8th Conf. Appl. Climatol.
17
(
22
),
179
183
.
Moazenzadeh
R.
Mohammadi
B.
Shamshirband
S.
Chau
K. W.
2018
Coupling a firefly algorithm with support vector regression to predict evaporation in northern Iran
.
Eng. Appl. Comput. Fluid Mech.
12
(
1
),
584
597
.
Mondal
A.
Kundu
S.
Mukhopadhyay
A.
2012
Rainfall trend analysis by Mann–Kendall test: a case study of north-eastern part of Cuttack district, Orissa
.
Int. J. Geol. Earth Environ. Sci.
2
(
1
),
70
78
.
Park
H. M.
2009
Comparing group means: t-tests and one-way ANOVA using Stata, SAS, R, and SPSS. Working Paper. The University Information Technology Services (UITS) Center for Statistical and Mathematical Computing, Indiana University, USA
.
Rahmat
S. N.
Jayasuriya
N.
Bhuiyan
M. A.
2015
Precipitation trends in Victoria, Australia
.
J. Water Clim. Change
6
(
2
),
278
287
.
doi:10.2166/wcc.2014.007
.
Shamshirband
S.
Hashemi
S.
Salimi
H.
Samadianfard
S.
Asadi
E.
Shadkani
S.
Kargar
K.
Mosavi
A.
Nabipour
N.
Chau
K. W.
2020
Predicting Standardized Streamflow index for hydrological drought using machine learning models
.
Eng. Appl. Comput. Fluid Mech.
14
(
1
),
339
350
.
Skidmore
S. T.
Thompson
B.
2013
Bias and precision of some classical ANOVA effect sizes when assumptions are violated
.
Behav. Res. Methods
45
(
2
),
536
546
.
doi:10.3758/s13428-012-0257-2
.
Tian
L.
Yuan
S.
Quiring
S. M.
2018
Evaluation of six indices for monitoring agricultural drought in the south-central United States
.
Agric. For. Meteorol.
249
,
107
119
.
doi:10.1016/j.agrformet.2017.11.024
.
Tigkas
D.
Vangelis
H.
Tsakiris
G.
2015
Drinc: a software for drought analysis based on drought indices
.
Earth Sci. Inform.
8
(
3
),
697
709
.
doi:10.1007/s12145-014-0178-y
.
Tsakiris
G.
Vangelis
H.
2005
Establishing a drought index incorporating evapotranspiration
.
Eur. Water
9
(
10
),
3
11
.
Tsakiris
G.
Pangalou
D.
Vangelis
H.
2007
Regional drought assessment based on the Reconnaissance Drought Index (RDI)
.
Water Resour. Manage.
21
(
5
),
821
833
.
doi:10.1007/s11269-006-9105-4
.
Yue
S.
Pilon
P.
Cavadias
G.
2002
Power of the Mann–Kendall and Spearman's rho tests for detecting monotonic trends in hydrological series
.
J. Hydrol.
259
(
1–4
),
254
271
.
doi:10.1016/S0022-1694(01)00594-7
.
Zarch
M. A. A.
Malekinezhad
H.
Mobin
M. H.
Dastorani
M. T.
Kousari
M. R.
2011
Drought monitoring by reconnaissance drought index (RDI) in Iran
.
Water Resour. Manage.
25
(
13
),
3485
.
doi:10.1007/s11269-011-9867-1
.