ABSTRACT
Drought is one of the most destructive natural disasters, becoming more extreme and less predictable due to climate change. Drought directly affects a region's water resources, leading to inadequate water availability and harming crops, animals, and humans. The objective of the present study is to assess the long-term trends of precipitation and drought characteristics for the Churu district of northeast Rajasthan, India. For this purpose, mean monthly precipitation data are collected through the India-WRIS from 1901 to 2022. The trend analysis uses statistical methods (M-K test & amp; SS estimator), and the results are compared with a graphical method (innovative trend analysis). Drought events are found using precipitation-based drought indices, i.e., standardized precipitation index, Z-score index, and percentage of normal precipitation index. Interestingly, the results of the statistical method are well matched with the results of the graphical method. The study concludes that during the annual precipitation a rising trend is observed, which can help to mitigate the severity of drought. Also, it is found that moderately dry droughts occur most frequently compared with other drought events. This study will help policy-makers and local administrators to take necessary action to mitigate the severity of drought.
HIGHLIGHTS
To assess the long-term trends of precipitation and drought characteristics.
Drought directly affects a region's water resources, leading to inadequate water availability and harming crops, animals, and humans.
INTRODUCTION
Globally, climate change receives more attention in relation to weather and intense incidents (Mehta & Yadav 2021a). In recent years, droughts have been occurring more frequently due to climate change and global warming (Surendran et al. 2019). Drought is one of the most destructive catastrophic events, because it directly influences the agricultural sector, environment, heath of the people, and the economy to a great extent (Fung et al. 2020; Mehta & Yadav 2021b). Drought is one of the most complicated water-related disasters because it is a phenomenon that can affect human activities more than other natural disasters (Mahmoudi et al. 2019). Drought is a temporary event which lasts for some duration, but it can destroy the health and wealth of the people living in that area. Generally, drought occurs when there is an absence of water for a long duration usually of months or years in a region (Jonathan & Suvarna Raju 2017). It can also be defined as a natural drop in precipitation over a long time-duration, which is generally a season or more, called meteorological drought (Mishra & Singh 2010; Verma et al. 2024). Generally, drought can take place in any region of the world, but its frequency and effects are more severe in arid and semi-arid regions (Wilhite 2000; Mehta & Yadav 2021c).
Meteorological drought usually occurrs when an area receives less than average precipitation during an extensive period of time (Quiring 2009). According to Burton et al. (2013), compared with other hydrological hazards like flood, it is very difficult to identify when droughts can occur. The development of drought policy to reduce its effects broadly depends on information about its severity and duration (Jain et al. 2015). Therefore, forecasting meteorological drought is necessary to decrease its adverse effects. This can be achieved through drought monitoring and forecasting by using drought indices (Chandrasekara et al. 2021; Verma et al. 2021; Verma et al. 2022a). A drought index is a number which is obtained by using hydrological and meteorological data as an input as well as drought indicators (Chandrasekara et al. 2021; Sahu et al. 2022, 2023). An index basically uses hydro-meteorological parameters like precipitation, temperature, runoff, streamflow, and evapotranspiration and gives a single numerical value as output (Patel et al. 2021; Mehta et al. 2022). The main benefit of employing drought indices is that they are easy to use as well as giving more accurate results to help policy-makers (Zargar et al. 2011; Halder et al. 2022; Verma et al. 2022b).
Also, it is mandatory to determine the past and present trend of precipitation for the study of droughts (Mehta & Yadav 2022a, 2022b). Analysis of precipitation trends helps in better forecasting of future droughts and sustainable water resources management (Perera et al. 2020; Verma et al. 2023a, 2023b). According to Patel & Mehta (2023), it is vital to recognize the trend of precipitation for drought and flood forecasting. Various researchers have analysed long-term precipitation trends and adopted different drought indices for the assessment of droughts based on precipitation data. Thomas & Prasannakumar (2016) applied different precipitation-based drought indices, namely, rainfall anomaly index (RAI), percentage departure of rainfall (D%), and standardized precipitation index (SPI) to determine the severity of droughts in Kerala. Mehta & Yadav (2022a, 2022b) used percentage departure of rainfall (D%), RAI, and SPI to study drought characteristics across Jalore district. Yadav et al. (2021) also utilized various drought indices, i.e., percentage departure (D%), RAI, and SPI to evaluate the severity of droughts (Shaikh et al. 2022; Mehta & Yadav 2023).
There are numerous drought indices that can be used to analyse the effect of droughts. However, in this study, precipitation-based drought indices (SPI, Z-score index (ZSI), and percentage of normal precipitation index (PNPI)) are used. To assess the trends in precipitation time-series, the widely known M-K test, Sen's slope (SS) estimator, and innovative trend analysis (ITA) method are used. The current research aims (i) to carry out the trend analysis of precipitation across the Churu district of Rajasthan and (ii) to examine the drought characteristics using precipitation-based drought indices.
STUDY AREA AND DATA COLLECTION
In this research, mean monthly precipitation data of the Churu district were gathered through the Indian Water Resource Information System (IWRIS) (indiawris.gov.in) for 1901 to 2022 (122 years).
METHODOLOGY
The current research work is carried out in two parts:
(1) trend analysis of precipitation time series;
(2) drought analysis based on precipitation.
Checking for effects of serial correlation
where implies the mean of sample size and n implies the size of sample.
If r1 is found between the given confidence interval, then the data are free from the effects of serial correlation, and there is no need to perform pre-whitening. Otherwise, pre-whitening should be conducted to terminate the effects of serial correlation (Yue et al. 2002).
M-K test and SS estimator method
In the present study, the M-K test was employed to detect the trend in precipitation time-series. The World Meteorological Organisation (WMO) suggested that the M-K test is best for identifying trends in time series (Aher & Yadav 2021). The benefit of using this test is that it is simple, it does not consider any assumptions, and it can handle the missing data (Porter et al. 2002). The M-K test was applied at a 5% significance level. To conduct the M-K test, XLSTAT 2022 software was used.
In the aforementioned equations, xj and xk denote precipitation data in years and n means the number of observations. If the value of S is positive, it means the data follow an increasing trend and vice-versa if it is negative, which shows a decreasing trend.
where q and t mean the tied group and number of observations, respectively.
If the value of ZMK is more than zero (ZMK > 0), then the data show an increasing trend. Furthermore, the data reveal a decreasing trend when the ZMK value is less than zero (ZMK < 0). On the other hand, the trend is considered as no trend if ZMK is zero (ZMK = 0).
If the value of βi > 0, it demonstrates upward trends in the data. Vice-versa, it reveals a downward trend if the value of βi > 0.
ITA method
To find trends in precipitation time-series, Şen (2012) discovered an ITA method. The main merit of the ITA method is that it is free from serial correlation, normality, or record length, compared with the frequently used trend analysis techniques, such as the M-K test and SS estimator (Pastagia & Mehta 2022). In this approach, the precipitation data are split into two equivalent parts from the start date to the last date (Caloiero 2020). Then both subseries are listed in either an increasing or a decreasing order. Further, the first part of subseries is plotted on the x-axis, while the second part of the subseries is plotted on the y-axis (Ali & Abubaker 2019). If every data point is obtained on a 1:1 line, it is interpreted as no trend existing in the data. On the other hand, data points that are acquired above the 1:1 line indicate a positive trend. However, a negative trend exists in the data when it is gathered below the 1:1 line (Besha et al. 2022).
Standardized precipitation index
SPI was found by McKee et al. (1993), which is one of the most popular drought indices to assess meteorological drought. Generally, it is used worldwide due to its simplicity as well as accuracy. The major benefit of using this index is that it requires only precipitation data as an input variable. Also, it can be used over multiple timescales to show the consequences of drought on water resources. In this study, three, six, and 12 months SPI are utilized to assess the intensity of drought over multiple time-scales. To calculate the SPI, DrinC software is used.
Table 1 illustrates the distribution of drought situations as per the SPI value.
SPI value . | Category . |
---|---|
More than 2.0 | Extremely wet |
1.5 to 1.99 | Severely wet |
1.0 to 1.49 | Moderately wet |
−0.99 to 0.99 | Near normal |
−1.0 to −1.49 | Moderately dry |
−1.5 to −1.99 | Severely dry |
Less than −2.0 | Extremely dry |
SPI value . | Category . |
---|---|
More than 2.0 | Extremely wet |
1.5 to 1.99 | Severely wet |
1.0 to 1.49 | Moderately wet |
−0.99 to 0.99 | Near normal |
−1.0 to −1.49 | Moderately dry |
−1.5 to −1.99 | Severely dry |
Less than −2.0 | Extremely dry |
Z-score index
Table 2 illustrates the distribution of drought situations as per the ZSI value.
ZSI value . | Category . |
---|---|
More than 2.0 | Extremely wet |
1.5 to 1.99 | Severely wet |
1.0 to 1.49 | Moderately wet |
−0.99 to 0.99 | Near normal |
−1.0 to −1.49 | Moderately dry |
−1.5 to −1.99 | Severely dry |
Less than −2.0 | Extremely dry |
ZSI value . | Category . |
---|---|
More than 2.0 | Extremely wet |
1.5 to 1.99 | Severely wet |
1.0 to 1.49 | Moderately wet |
−0.99 to 0.99 | Near normal |
−1.0 to −1.49 | Moderately dry |
−1.5 to −1.99 | Severely dry |
Less than −2.0 | Extremely dry |
Percentage of normal precipitation index
PNPI is the simple drought index found by Willeke et al. (1994). It is calculated by dividing the actual precipitation with its average value during a particular time period and expressed as percentage. The main benefit of using PNPI is that only precipitation data are needed for its computation, and it is used on the monthly timescale (Willeke et al. 1994).
Table 3 illustrates the distribution of drought situations as per the PNPI value.
PNPI value (%) . | Category . |
---|---|
70–80 | Normal |
55–70 | Moderately dry |
40–55 | Severely dry |
Less than 40 | Extremely dry |
PNPI value (%) . | Category . |
---|---|
70–80 | Normal |
55–70 | Moderately dry |
40–55 | Severely dry |
Less than 40 | Extremely dry |
RESULTS AND DISCUSSION
Effects of serial correlation in precipitation time-series
To investigate the effects of serial correlation the lag-1 autocorrelation coefficient (r1) was applied at a 5% significance level to the seasonal (winter, pre-monsoon, southwest (SW) monsoon, and post-monsoon) and annual precipitation time-series of the Churu district. The test revealed that the seasonal (winter, pre-monsoon, SW monsoon, and post-monsoon) time-series had the following values of the autocorrelation coefficient: −0.055, 0.235, 0.126, and 0.307, respectively. Furthermore, the annual time series was observed with r1 = 0.115. As the results of serial correlation were found within the permitted range (−1.178 to 1.178), it can be concluded that none of the time series (either seasonal or annual) has any significant serial correlation. Hence, in this research, the M-K test is conducted on the original precipitation data excluding the pre-whitening approach.
Results of precipitation trend analysis using M-K test and SS estimator
Precipitation duration . | Kendall's tau . | Trend Interpretation . | M-K test p-value . | Sen's slope . | Test interpretation . |
---|---|---|---|---|---|
January | −0.003 | Falling | 0.966 | 0 | No trend |
February | 0.005 | Rising | 0.931 | 0 | No trend |
March | 0.114 | Rising | 0.065 | 0.009 | Insignificant increasing trend |
April | 0.147 | Rising | 0.017 | 0.015 | Significant increasing trend |
May | 0.215 | Rising | 0 | 0.107 | Significant increasing trend |
June | 0.209 | Rising | 0.001 | 0.233 | Significant increasing trend |
July | 0.162 | Rising | 0.008 | 0.416 | Significant increasing trend |
August | 0.035 | Rising | 0.571 | 0.097 | Insignificant increasing trend |
September | 0.053 | Rising | 0.392 | 0.068 | Insignificant increasing trend |
October | 0.072 | Rising | 0.246 | 0.001 | Insignificant increasing trend |
November | 0.135 | Rising | 0.042 | 0 | No trend |
December | −0.104 | Falling | 0.107 | 0 | No trend |
Winter | 0.01 | Rising | 0.873 | 0.004 | Insignificant increasing trend |
Pre-monsoon | 0.201 | Rising | 0.001 | 0.152 | Significant increasing trend |
SW monsoon | 0.158 | Rising | 0.01 | 0.758 | Significant increasing trend |
Post-monsoon | −0.079 | Falling | 0.197 | −0.02 | Insignificant decreasing trend |
Annual | 0.199 | Rising | 0.001 | 1.041 | Significant increasing trend |
Precipitation duration . | Kendall's tau . | Trend Interpretation . | M-K test p-value . | Sen's slope . | Test interpretation . |
---|---|---|---|---|---|
January | −0.003 | Falling | 0.966 | 0 | No trend |
February | 0.005 | Rising | 0.931 | 0 | No trend |
March | 0.114 | Rising | 0.065 | 0.009 | Insignificant increasing trend |
April | 0.147 | Rising | 0.017 | 0.015 | Significant increasing trend |
May | 0.215 | Rising | 0 | 0.107 | Significant increasing trend |
June | 0.209 | Rising | 0.001 | 0.233 | Significant increasing trend |
July | 0.162 | Rising | 0.008 | 0.416 | Significant increasing trend |
August | 0.035 | Rising | 0.571 | 0.097 | Insignificant increasing trend |
September | 0.053 | Rising | 0.392 | 0.068 | Insignificant increasing trend |
October | 0.072 | Rising | 0.246 | 0.001 | Insignificant increasing trend |
November | 0.135 | Rising | 0.042 | 0 | No trend |
December | −0.104 | Falling | 0.107 | 0 | No trend |
Winter | 0.01 | Rising | 0.873 | 0.004 | Insignificant increasing trend |
Pre-monsoon | 0.201 | Rising | 0.001 | 0.152 | Significant increasing trend |
SW monsoon | 0.158 | Rising | 0.01 | 0.758 | Significant increasing trend |
Post-monsoon | −0.079 | Falling | 0.197 | −0.02 | Insignificant decreasing trend |
Annual | 0.199 | Rising | 0.001 | 1.041 | Significant increasing trend |
Results of precipitation trend analysis using ITA method
Season . | Low . | Medium . | High . |
---|---|---|---|
Winter | ↑ | ↓ | ↓ |
Pre-monsoon | ↑ | ↑ | ↑ |
SW monsoon | ↑ | ↑ | ↑ |
Post-monsoon | ↓ | ↑ | O |
Annual | ↑ | ↑ | ↑ |
Season . | Low . | Medium . | High . |
---|---|---|---|
Winter | ↑ | ↓ | ↓ |
Pre-monsoon | ↑ | ↑ | ↑ |
SW monsoon | ↑ | ↑ | ↑ |
Post-monsoon | ↓ | ↑ | O |
Annual | ↑ | ↑ | ↑ |
Comparative trend analysis of precipitation between M-K test, SS estimator, and ITA method
In this research, trends in monthly, seasonal, and annual time-series of precipitation are obtained employing the M-K test and SS estimator, and compared with the ITA method. Table 6 shows the comparative trend analysis of precipitation results between these three methods. By comparing results among these methods, it is found that, for the annual precipitation, all three methods (M-K test, SS estimator, and ITA) show an increasing trend. Also, pre-monsoon and SW monsoon seasons are observed with an increasing trend for all three methods (M-K test, SS estimator, and ITA). Moreover, the winter season reveals an increasing trend for statistical methods, while the graphical method represents the decreasing trend, which is highlighted in Table 6. Only the post-monsoon season is found with a decreasing trend in all three methods (M-K test, SS estimator, and ITA).
Season . | M-K test . | SS estimator . | ITA method . |
---|---|---|---|
Winter | Increasing | Increasing | Decreasing |
Pre-monsoon | Increasing | Increasing | Increasing |
SW monsoon | Increasing | Increasing | Increasing |
Post-monsoon | Decreasing | Decreasing | Decreasing |
Annual | Increasing | Increasing | Increasing |
Season . | M-K test . | SS estimator . | ITA method . |
---|---|---|---|
Winter | Increasing | Increasing | Decreasing |
Pre-monsoon | Increasing | Increasing | Increasing |
SW monsoon | Increasing | Increasing | Increasing |
Post-monsoon | Decreasing | Decreasing | Decreasing |
Annual | Increasing | Increasing | Increasing |
Drought analysis using SPI
In this research, SPI is computed for three different timescales (three, six and 12 months). Short-term drought events are analysed using a three-month SPI, and to analyse intermediate drought events, six-month SPI is employed. Twelve-month SPI is utilized to analyse long-term drought events. The calculation of SPI is done for the hydrological year October to September by using DrinC software.
Drought analysis using Z-score index
In this study, ZSI is also calculated for the three timescales, i.e., three, six and 12 months.
Drought analysis using percentage of normal precipitation index
In this study, PNPI is calculated for 12 months for the hydrological year from October to September.
CONCLUSION
In this study, the semi-arid district of northeast Rajasthan (Churu district) is selected and trend analysis is carried out for monthly, seasonal, and annual precipitation. To carry out the trend analysis, statistical methods (M-K test and SS estimator) for the time period of 1901–2022 (122 years) are used. However, the graphical method (ITA) is also utilized to detect a trend in the precipitation time-series. Also, to evaluate the severity of drought in the study area, various drought indices such as SPI, ZSI, and PNPI are utilized. The study concludes the following points:
The winter season shows a decreasing trend for the M-K test and SS estimator. On the other hand, the ITA method shows an increasing trend.
For pre-monsoon, SW monsoon, and annual seasons, an increasing trend is noticed by applying all three methods (M-K test, SS estimator, and ITA method).
Only the post-monsoon season revealed a decreasing trend under all three methods (M-K test, SS estimator, and ITA method).
Except only for the winter season, the results of M-K test and SS estimator match the results of the ITA method.
The graphical method (ITA method) is more accurate in detecting the trend because it does not consider any assumptions like the M-K test.
During the annual season, the rising trend of precipitation can help to mitigate the severity of drought.
From the results of drought analysis, it can be said that moderately dry events take place more frequently compared with other drought events.
The results of 12-month SPI and 12-month ZSI show that an extremely dry event occurred during 1918 with the severity of −3.21 and −2.11, respectively.
This study will help policy-makers and local administrators to take necessary action to mitigate the severity of drought as well as better management of water resources in the Churu district.
DECLARATIONS
All authors have read, understood, and have complied as applicable with the statement on ‘Ethical responsibilities of Authors’ as found in the instructions for authors.
FUNDING
This research received no external funding.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.