*Article* **Shifting of Meteorological to Hydrological Drought Risk at Regional Scale**

**Awais Naeem Sarwar <sup>1</sup> , Muhammad Waseem 1,\* , Muhammad Azam <sup>2</sup> , Adnan Abbas <sup>3</sup> , Ijaz Ahmad <sup>1</sup> , Jae Eun Lee 4,\* and Faraz ul Haq <sup>1</sup>**


**Abstract:** The drought along with climate variation has become a serious issue for human society and the ecosystem in the arid region like the Soan basin (the main source of water resources for the capital of Pakistan and the Pothohar arid region). The increasing concerns about drought in the study area have brought about the necessity of spatiotemporal analysis and assessment of the linkage between different drought types for an early warning system. Hence, the streamflow drought index (SDI) and standard precipitation index (SPI) were used for the analysis of the spatiotemporal variations in hydrological and meteorological drought, respectively. Furthermore, statistical approaches, including regression analysis, trend analysis using Mann Kendall, and moving average, have been used for investigation of the linkage between these drought types, the significance of the variations, and lag time identification, respectively. The overall analysis indicated an increase in the frequency of both hydrological and meteorological droughts during the last three decades. Moreover, a strong linkage between hydrological and meteorological droughts was found; and this relationship varied on the spatiotemporal scale. Significant variations between hydrological and meteorological droughts also resulted during the past three (3) decades. These discrepancies would be because of different onset and termination times and specific anthropogenic activities in the selected basin for the minimization of hydrological drought. Conclusively, the present study contributes to comprehending the linkage between hydrological and meteorological droughts and, thus, could have a practical use for local water resource management practices at the basin scale.

**Keywords:** drought SDI; SPI; linkage; propagation

### **1. Introduction**

Drought is a condition considered as the deficit of water, including surface, ground, or atmospheric water, for a long period [1]. Irrespective of climatic sites, the drought could occur worldwide even in humid and wet environments [2]. Drought is generally classified into meteorological, hydrological, socioeconomic, and agricultural drought [3]. Among these four types of droughts, hydrological is the most important form as sustainable water resources management is heavily dependent on hydro-meteorological information. Various numbers of factors are associated with the onset of hydrological drought, and meteorological drought is one of the main influencing factors. The root cause of drought onset is the deficiency of rainfall across a large area for a very long period and is notated as a meteorological drought [4]. A meteorological drought can develop quickly because it is primarily caused by a lack of precipitation [5], and if the lack of precipitation spreads

**Citation:** Sarwar, A.N.; Waseem, M.; Azam, M.; Abbas, A.; Ahmad, I.; Lee, J.E.; Haq, F.u. Shifting of Meteorological to Hydrological Drought Risk at Regional Scale. *Appl. Sci.* **2022**, *12*, 5560. https://doi.org/ 10.3390/app12115560

Academic Editors: Andrea Chiozzi, Željana Nikoli´c and Elena Benvenuti

Received: 24 March 2022 Accepted: 27 May 2022 Published: 30 May 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

to specific regions, the meteorological drought can be transformed into a hydrological drought, and then into an agricultural drought. The drought propagation process refers to the transition of different types of droughts.

The relationships between different types of droughts are studied by several authors, for example, [6] resulted that hydrological drought events occurred approximately seven months after meteorological drought events. The streamflow drought can be investigated using meteorological drought information at the annual time scale [7]. The interrelationship between runoff and meteorological drought has also been investigated and it was concluded that the significant dependency of hydrological drought on meteorological drought [8]. The impact of meteorological forcing on hydrological droughts has been computed [9]. A copula-based joint meteorological-hydrological drought index has been used to model the relationship between meteorological and hydrological droughts upstream and downstream of the Kasilian basin [10]. The transmission of meteorological droughts to hydrological droughts and the influence factors have been explored and investigated [11]. The entropy theory was used to create a hybrid drought index that combines hydrological, meteorological, and agricultural data and was used to investigate the drought condition in Northwest China [12]. In literature, the relationships between hydro-meteorological droughts, climatic variables, and human activities were also explored [13,14]. For example, [15] examined the relationship between various hydrologic and meteorological drought indices considering natural and human factors for different basins in the contiguous United States. Similarly, [16–19] studied the dependence structure between meteorological drought and hydrological drought using different approaches and resulted that the occurring time lag between meteorological and hydrological droughts is critical for sustainable drought management.

Most studies in the literature focused on determining hydro-meteorological drought and drought propagation; however, it remains unexplored. Besides that, the lack of understanding considering hydrological drought response to meteorological drought in different regions raises an unanswered question for basin-scale drought risk management [20]. To find answers to these questions, a comprehensive study of droughts that span multiple geographic areas and last for extended periods is required. Strategic water resources can be implemented more efficiently with accurate water-based information based on regional drought characteristics. Furthermore, given Pakistan's erratic, scarce, unstable climate and current drought situation [20–29], drought propagation information is critical at the regional and national levels, as it can provide the appropriate and consistent information required for efficient water management and drought early warning systems. Drought propagation information can aid in the prevention of significant economic losses as well as decision-making. Early warning systems are estimated to save hundreds of lives, save 2.7 billion Euros in natural disaster losses, and generate billions of supplementary benefits per year in Europe based on the optimization of economic production in the energy and agriculture sectors. In the near past, Pakistan has also experienced recurrent droughts and more severe droughts are anticipated in near future due to climate change [22].

The Soan River is a tributary of Pakistan's mighty Indus River Basin. It flows from the Murree Mountains into the Indus via the Dhoke Pathan hydrological station. It is the main hydrological unit for the Pothohar arid region of Pakistan. The importance of the Pothohar region cannot be repudiated as the region of Pothohar is rich in agriculture and agriculture is the backbone of the country's economy [30]. However, there are canal irrigation resources in this area and the crop production depends entirely on rainfall at suitable timings. Pothohar region gets rainfall during both the winter and summer seasons, which is the main source of crop germination, flowering, and maturity. If any of the seasons fail to bring rain, then soil moisture depletes resulting in drought and ultimately significant crop damage. More specifically, the crops in these areas depend on monsoon rains in summer and on western rains in winter. If the summer season fails to bring rains in the area, then it will cause huge crop yield loss in that specific year. Similarly, if the westerly system also fails to bring rains to the area, crops are affected badly, resulting in severe

droughts. Historically, this Pothohar region was severely affected by substantial variation in the rainfall and prolonged rain shortage especially from 1999 to 2002 [31]. Hence, the objectives of the present study include the exploration of spatiotemporal evolutions of meteorological and hydrological droughts in the study basin, evaluation of the link between meteorological and hydrological droughts, identifying the lag time, and the investigation of the differences between hydrological and meteorological droughts. variation in the rainfall and prolonged rain shortage especially from 1999 to 2002 [31]. Hence, the objectives of the present study include the exploration of spatiotemporal evolutions of meteorological and hydrological droughts in the study basin, evaluation of the link between meteorological and hydrological droughts, identifying the lag time, and the investigation of the differences between hydrological and meteorological droughts. **2. Materials and Methods**

in the area, then it will cause huge crop yield loss in that specific year. Similarly, if the westerly system also fails to bring rains to the area, crops are affected badly, resulting in severe droughts. Historically, this Pothohar region was severely affected by substantial

*Appl. Sci.* **2022**, *12*, x FOR PEER REVIEW 3 of 14

#### **2. Materials and Methods** *2.1. Study Area*

### *2.1. Study Area* The Soan River, initiating from the Murree highlands, is a major tributary of the In-

The Soan River, initiating from the Murree highlands, is a major tributary of the Indus River and a key hydrological entity of the Potohar region of Pakistan. The Soan River flows through Chirah and Dhoke Pathan hydrological gauging stations (Figure 1) before being a part of the Indus River. The total drainage area of the Soan River basin is 6842 km<sup>2</sup> with an elevation ranging from 265 to 2274 m. The Soan basin is characterized by gentle to the steep slope and monsoon fed streams which generate almost all of the basin flow with mean annual precipitation of 1465 mm, and the mean annual temperature ranges from 8 to 22 ◦C. The Simly Dam, which spans the Soan River, provides drinking water to Islamabad, Pakistan's capital, as well as water for irrigation activities in the Pothohar region [32]. Furthermore, for a better understanding and evaluation of the spatial variation of drought, the study area was divided into two sub-basins: Chirah (hereafter sub-basin 1) and Dhoke Pathan (hereafter sub-basin 2). In the current study, monthly precipitation data of five meteorological stations within the basin (Murree, Zero Point, airport, Dhamial, and SAWCRI) from 1983 to 2015 was collected from Pakistan Meteorological Department (PMD). The average value of four stations lying in sub-basin 2 was used to estimate the meteorological drought in sub-basin 2, whereas the data of the remaining one station, i.e., Murree was used for meteorological drought analysis at sub-basin 1. Similarly, the monthly streamflow data of two hydrological gauging stations, i.e., Chirah and Dhoke Pathan was collected from Water and Power Development Authority (WAPDA) for the same period. dus River and a key hydrological entity of the Potohar region of Pakistan. The Soan River flows through Chirah and Dhoke Pathan hydrological gauging stations (Figure 1) before being a part of the Indus River. The total drainage area of the Soan River basin is 6842 km<sup>2</sup> with an elevation ranging from 265 to 2274 m. The Soan basin is characterized by gentle to the steep slope and monsoon fed streams which generate almost all of the basin flow with mean annual precipitation of 1465 mm, and the mean annual temperature ranges from 8 to 22 °C. The Simly Dam, which spans the Soan River, provides drinking water to Islamabad, Pakistan's capital, as well as water for irrigation activities in the Pothohar region [32]. Furthermore, for a better understanding and evaluation of the spatial variation of drought, the study area was divided into two sub-basins: Chirah (hereafter sub-basin 1) and Dhoke Pathan (hereafter sub-basin 2). In the current study, monthly precipitation data of five meteorological stations within the basin (Murree, Zero Point, airport, Dhamial, and SAWCRI) from 1983 to 2015 was collected from Pakistan Meteorological Department (PMD). The average value of four stations lying in sub-basin 2 was used to estimate the meteorological drought in sub-basin 2, whereas the data of the remaining one station, i.e., Murree was used for meteorological drought analysis at sub-basin 1. Similarly, the monthly streamflow data of two hydrological gauging stations, i.e., Chirah and Dhoke Pathan was collected from Water and Power Development Authority (WAPDA) for the same period.

**Figure 1.** Location map of the study area. **Figure 1.** Location map of the study area.

### *2.2. Spatiotemporal Analysis of Drought 2.2. Spatiotemporal Analysis of Drought* 2.2.1. Meteorological Drought Assessment

2.2.1. Meteorological Drought Assessment For the assessment of meteorological drought, the Standardized Precipitation Index (SPI) was used due to its effectiveness, applicability, and suitability at different time scales for many case studies, e.g., [28,33,34]. Considering the hydrological year (October to September), SPI was calculated at four different time scales, i.e., SPI-3 (October–December), For the assessment of meteorological drought, the Standardized Precipitation Index (SPI) was used due to its effectiveness, applicability, and suitability at different time scales for many case studies, e.g., [28,33,34]. Considering the hydrological year (October to September), SPI was calculated at four different time scales, i.e., SPI-3 (October–December), SPI-6 (October–March), SPI-9 (October–June), and SPI-12 (October–September). Furthermore, for the calculation of SPI at different time scales, the following expression was used for cumulative precipitation (*Ri*,*K*).

$$R\_{i,K} = \sum\_{j=1}^{3k} P\_{i,j} \text{ i } = 1, 2, 3, \dots, n \text{ j } = 1, 2, \dots, \dots, 12 \text{ and } k = 1, 2, 3, 4 \tag{1}$$

where *Pi*,*<sup>j</sup>* corresponds to monthly precipitation, *j* denotes the particular month of the hydrological year, *I* indicates the hydrological year, and k reveals the timescale (e.g., SPI-3, *k* = 1 indicating 3-month time scale October-December, and similarly *k* = 4 is for October to September).

SPI computation includes fitting distributions to precipitation data (*Ri*,*K*) and estimating a probability density function (PDF) and cumulative distribution function (CDF). These functions are further transformed into standardized distributions with unit standard deviations and zero means yielding the SPI value using the expression below (Equation (2)). The SPI value could range from <−2 to >+2 [35], whereas the positive value of SPI indicates a wet condition while the negative value represents a drought condition. Moreover, in this study, the drought thresholds were as follows: weak drought (−1.0 ≤ index < 0), moderate drought (−1.49 ≤ index < −1.0), severe drought (−1.99 ≤ index < −1.50), and extreme drought (index ≤ −2).

$$\text{SPI}\_{i,k} = \frac{R\_{i,k} - R\_{i,k}}{S\_k} \tag{2}$$

where *S<sup>k</sup>* and *R<sup>k</sup>* are the standard deviation and average value of the precipitation, respectively.

### 2.2.2. Hydrological Drought Assessment

The Streamflow Drought Index (SDI) was developed by [36] and is widely used to characterize hydrological drought events. SDI computation also includes fitting distributions to runoff data and estimation of PDF and CDF and then transforming into a standardized distribution, yielding the SDI value [13]. Similarly, the positive SDI value represents wet conditions, while the negative value represents drought conditions. Furthermore, SPI thresholds were adopted for hydrological drought analysis.

$$V\_{i,k} = \sum\_{j=1}^{3k} Q\_{i,j} \text{ i } = 1, 2, 3, \dots, n \text{ j } = 1, 2, 3, \dots, 12 \text{ and } k = 1, 2, 3, 4 \tag{3}$$

$$\text{SDI}\_{i,k} = \frac{V\_{i,k} - \overline{V}\_{i,k}}{S\_k} \tag{4}$$

where *Vi*,*<sup>k</sup>* the cumulative runoff for reference period *k* and the *i* is the hydrological year; *S<sup>k</sup>* is the standard deviation and *V<sup>k</sup>* is the mean value of runoff.

The SDI can also be calculated at various time scales, such as one, three, six, and twelve months. For the estimation of SPI-3, SDI-3, SPI-6, SDI-6, SPI-9, SDI-9, SPI-12, and SDI-12, the cumulative sums of precipitation and runoff for 3, 6, 9, and 12-month timescales were used. When calculating SPI-3 and SDI-3, for example, the cumulative sum for October was obtained by adding the following two-month data (November and December) to the October data [20]. SPI-6 and SDI-6 cumulative sums were calculated by adding six months from October to March.

Furthermore, the Mann-Kendall method [37] was used to estimate the trend in SPI-3, SPI-6, SPI-9, and SPI-12 at a significance level of 5% and similarly for SDI values at different time scales.

### *2.3. Identification of Linkage between Meteorological and Hydrological Drought and Lag Time*

A regression analysis is a statistical technique for estimating the magnitude of the effect of any independent time series (hereafter the meteorological drought index) on any dependent time series (hereafter the hydrological drought index) [38]. In particular,

regression analysis is used to better understand the relationships between independent and dependent time series. As a result, a linear regression was used to evaluate the relationship and temporal changes in drought. Coefficient of determination (R<sup>2</sup> ), slope, and *p*-value were calculated for each relationship using the two-tailed t-test, and R2 was used to indicate the best relationship between SDI and SPI [38]. Equation (5) depicts the mathematical expression of regression between any independent and any dependent time. In the current study, several simple linear regression models were developed between SPI and SDI of the same time scale, e.g., SPI-3 versus SDI-3; SPI-6 versus SDI-6; and similarly for other time scales.

$$Y = \beta\_0 + aX + \varepsilon \tag{5}$$

where *Y* is the dependent time series, *a* is the regression coefficient, *X* is the independent time series, *β*<sup>0</sup> is the intercept, and *ε* is an error term.

Furthermore, in this study, moving average (running average) analysis was performed at multiple time scales to better assess the correlation between SPI and SDI. It is a quite simple statistical analysis tool used as the fluctuation, trend, or lagging indicator. For a given time series and a fixed size subset, the first value of the moving average is calculated by averaging the fixed subset of the given time series. After this, the subset is shifting forward by including the next number and excluding the first number in the subset. For example, Let (*X*1, *t*1), (*X*2, *t*2), . . . . . . , (*Xn*, *tn*) represents SPI or SDI time series, *X*1, *X*2, . . . , *X<sup>n</sup>* are the values of SPI or SDI; against time periods *t*1, *t*2, . . . , *tn*, respectively. The first two moving averages of order *k* (here 2) can be calculated as given in Equations (6) and (7). A similar procedure was adopted up to the nth value.

$$\begin{array}{c} \text{For 1st MA}: \\ \text{MA}\_K = \frac{\text{X}\_1 + \text{X}\_2 + ... \dots + \text{X}\_k}{K} \end{array} \tag{6}$$

$$\begin{array}{c} \text{For 2nd MA}: \\ \text{MA}\_K = \frac{\text{X}\_2 + \text{X}\_3 + ... \dots + \text{X}\_{k+1}}{K} \end{array} \tag{7}$$

### **3. Results**

### *3.1. Spatiotemporal Analysis of Droughts*

### 3.1.1. Analysis of Meteorological Drought

Figure 2 illustrates the meteorological drought index calculated at selected time scales for both sub-basins of the study area. Based on the analysis, it was observed that there was an increase in the frequency of drought after 1998 specifically in the case of sub-basin-1 when SPIs mostly remained negative. These results are in good agreement with the results previously published on droughts in Pakistan [19,20,23,29,39]. The variations in drought between both sub-basins are due to the more relative decrease in precipitation in subbasin-1, compared to sub-basin-2 [30,32]. Moreover, the numbers of drought events were different at different computed time scales. Significant inconsistencies were found when the 3-month time scale was compared with the 6-month time scale, and a comparison between the 6-month and 12-month time scales gave slight variations in results. Results reveal the frequent occurrence of severe drought events at each sub-basin of the study area, and these were more sensitive to short-term drought compared to long-term droughts. According to the results, the recurrence of drought is more frequent and consecutive at 1 to 6 months' timescale in sub-basin 1 and during the 12-month timescale, there exists a sharp difference between wet and dry episodes. while the sub-basin 2 is found to be more sensitive to drought (mild to severe) at a 1-month timescale as compared to 3, 6, and 12-month timescales. The results also specified that the drought events state varied with the increase in SPI time scales, mostly a declining trend in sub-basin 2. This is due to the reason that computation of 12-month time scale SPI involves the aggregation of total precipitation from October to September and includes both wet and dry seasons. However, the 3-month SPI time scale considers only the sum of three-month precipitation, and the 3-month time scale (March-May) of the study basin was normally a dry season.

However, the 3-month SPI time scale considers only the sum of three-month precipitation, and the 3-month time scale (March-May) of the study basin was normally a dry season.

However, the 3-month SPI time scale considers only the sum of three-month precipitation, and the 3-month time scale (March-May) of the study basin was normally a dry season.

*Appl. Sci.* **2022**, *12*, x FOR PEER REVIEW 6 of 14

**Figure 2.** Temporal variations of meteorological drought periods during 1983–2015 in sub–basins 1 and 2. **Figure 2.** Temporal variations of meteorological drought periods during 1983–2015 in sub–basins 1 and 2. **Figure 2.** Temporal variations of meteorological drought periods during 1983–2015 in sub–basins 1 and 2.

Figure 3 illustrates the severity of meteorological drought based on the 12-month time scale for both hydrological sub-basins and the difference in their SPI-12 values during the study period (i.e., 1983–2015). Based on the analysis, it resulted that the years from 1998 to 2004 were the driest years of the selected time series. Moreover, the difference in meteorological drought (SPI-12) severity at both sub-basins represents the notable spatial variations across the Soan river basins. The Mann-Kendall test was also used on SPI-12 to figure out the trend in meteorological drought at sub-basin 1 and sub-basin 2. The analysis showed an increase in meteorological drought severity during the past 33 years, and a significant decreasing trend in SPI–12 values was observed downstream, i.e., sub-basin 2. The Z-value of the Mann-Kendall test was 1.78 at upstream sub-basin 1, and 1.71 at downstream sub-basin 2. Figure 3 illustrates the severity of meteorological drought based on the 12-month time scale for both hydrological sub-basins and the difference in their SPI-12 values during the study period (i.e., 1983–2015). Based on the analysis, it resulted that the years from 1998 to 2004 were the driest years of the selected time series. Moreover, the difference in meteorological drought (SPI-12) severity at both sub-basins represents the notable spatial variations across the Soan river basins. The Mann-Kendall test was also used on SPI-12 to figure out the trend in meteorological drought at sub-basin 1 and sub-basin 2. The analysis showed an increase in meteorological drought severity during the past 33 years, and a significant decreasing trend in SPI–12 values was observed downstream, i.e., sub-basin 2. The Z-value of the Mann-Kendall test was 1.78 at upstream sub-basin 1, and 1.71 at downstream sub-basin 2. Figure 3 illustrates the severity of meteorological drought based on the 12-month time scale for both hydrological sub-basins and the difference in their SPI-12 values during the study period (i.e., 1983–2015). Based on the analysis, it resulted that the years from 1998 to 2004 were the driest years of the selected time series. Moreover, the difference in meteorological drought (SPI-12) severity at both sub-basins represents the notable spatial variations across the Soan river basins. The Mann-Kendall test was also used on SPI-12 to figure out the trend in meteorological drought at sub-basin 1 and sub-basin 2. The analysis showed an increase in meteorological drought severity during the past 33 years, and a significant decreasing trend in SPI–12 values was observed downstream, i.e., sub-basin 2. The Z-value of the Mann-Kendall test was 1.78 at upstream sub-basin 1, and 1.71 at downstream sub-basin 2.

**Figure 3.** Difference in meteorological drought (SPI–12) severity at both sub-basins representing the spatial variations across the Soan river basins. **Figure 3.** Difference in meteorological drought (SPI–12) severity at both sub-basins representing the spatial variations across the Soan river basins. **Figure 3.** Difference in meteorological drought (SPI–12) severity at both sub-basins representing the spatial variations across the Soan river basins.

### 3.1.2. Analysis of Hydrological Droughts 3.1.2. Analysis of Hydrological Droughts

3.1.2. Analysis of Hydrological Droughts Figure 4 depicts the temporal variations of hydrological drought periods in sub-basins 1 and 2 at different time scales, indicating the occurrence of hydrological drought over the last three decades. Hydrological drought, like meteorological drought, was more Figure 4 depicts the temporal variations of hydrological drought periods in sub-basins 1 and 2 at different time scales, indicating the occurrence of hydrological drought over the last three decades. Hydrological drought, like meteorological drought, was more Figure 4 depicts the temporal variations of hydrological drought periods in sub-basins 1 and 2 at different time scales, indicating the occurrence of hydrological drought over the last three decades. Hydrological drought, like meteorological drought, was more sensitive to short-term droughts than long-term droughts. Drought events occurred continuously from 1983 to 2015, with the driest period for the study area being 1998–2004. Drought patches were also randomly observed in 1993, 1998–2004, and 2009. Sub-basin 1 was stressed by hydrological drought, and the number of drought events was high from 1985 to 2004, with a slight wet condition from 2005-to 2015. However, in the case of sub-basin 2,

the results showed that the frequency of hydrological drought events was increasing in the basin from 1983 to 2015. Figure 5 shows the illustration of the difference in hydrological drought severity (SDI-12) for two sub-basins and indicates notable variations across the two sub-basin. Furthermore, based on the MK test, a statistically significant trend was observed for downstream sub-basins 2 with a Z-value of 1.83, whereas, upstream sub-basin 1 showed opposite behavior with an insignificant statistical trend. However, there still exists an increasing trend of drought at upstream sub-basin 1. basin 2, the results showed that the frequency of hydrological drought events was increasing in the basin from 1983 to 2015. Figure 5 shows the illustration of the difference in hydrological drought severity (SDI-12) for two sub-basins and indicates notable variations across the two sub-basin. Furthermore, based on the MK test, a statistically significant trend was observed for downstream sub-basins 2 with a Z-value of 1.83, whereas, upstream sub-basin 1 showed opposite behavior with an insignificant statistical trend. However, there still exists an increasing trend of drought at upstream sub-basin 1. basin 2, the results showed that the frequency of hydrological drought events was increasing in the basin from 1983 to 2015. Figure 5 shows the illustration of the difference in hydrological drought severity (SDI-12) for two sub-basins and indicates notable variations across the two sub-basin. Furthermore, based on the MK test, a statistically significant trend was observed for downstream sub-basins 2 with a Z-value of 1.83, whereas, upstream sub-basin 1 showed opposite behavior with an insignificant statistical trend. However, there still exists an increasing trend of drought at upstream sub-basin 1.

sensitive to short-term droughts than long-term droughts. Drought events occurred continuously from 1983 to 2015, with the driest period for the study area being 1998–2004. Drought patches were also randomly observed in 1993, 1998–2004, and 2009. Sub-basin 1 was stressed by hydrological drought, and the number of drought events was high from 1985 to 2004, with a slight wet condition from 2005-to 2015. However, in the case of sub-

sensitive to short-term droughts than long-term droughts. Drought events occurred continuously from 1983 to 2015, with the driest period for the study area being 1998–2004. Drought patches were also randomly observed in 1993, 1998–2004, and 2009. Sub-basin 1 was stressed by hydrological drought, and the number of drought events was high from 1985 to 2004, with a slight wet condition from 2005-to 2015. However, in the case of sub-

*Appl. Sci.* **2022**, *12*, x FOR PEER REVIEW 7 of 14

*Appl. Sci.* **2022**, *12*, x FOR PEER REVIEW 7 of 14

**Figure 4.** Temporal variation of hydrological drought periods in sub-basins 1 and 2 at different time scales. **Figure 4.** Temporal variation of hydrological drought periods in sub-basins 1 and 2 at different time scales. **Figure 4.** Temporal variation of hydrological drought periods in sub-basins 1 and 2 at different time scales.

**Figure 5.** Difference in Hydrological drought (SDI–12) severity at both sub-basins representing the spatial variations across the Soan river basin. **Figure 5.** Difference in Hydrological drought (SDI–12) severity at both sub-basins representing the spatial variations across the Soan river basin. **Figure 5.** Difference in Hydrological drought (SDI–12) severity at both sub-basins representing the spatial variations across the Soan river basin.
