Application of the Hidden Markov Bayesian Classiﬁer and Propagation Concept for Probabilistic Assessment of Meteorological and Hydrological Droughts in South Korea

: Drought is one of the most destructive natural hazards and results in negative e ﬀ ects on the environment, agriculture, economics, and society. A meteorological drought originates from atmospheric components, while a hydrological drought is inﬂuenced by properties of the hydrological cycle and generally induced by a continuous meteorological drought. Several studies have attempted to explain the cross dependencies between meteorological and hydrological droughts. However, these previous studies did not consider the propagation of drought classes. Therefore, in this study, to consider the drought propagation concept and to probabilistically assess the meteorological and hydrological drought classes, characterized by the Standardized Precipitation Index (SPI) and Standardized Runo ﬀ Index (SRI), respectively, we employed the Markov Bayesian Classiﬁer (MBC) model that combines the procedure of iteration of feature extraction, classiﬁcation, and application for assessment of drought classes for both SPI and SRI. The classiﬁcation results were compared using the observed SPI and SRI, as well as with previous ﬁndings, which demonstrated that the MBC was able to reasonably determine drought classes. The accuracy of the MBC model in predicting all the classes of meteorological drought varies from 36 to 76% and in predicting all the classes of hydrological drought varies from 33 to 70%. The advantage of the MBC-based classiﬁcation is that it considers drought propagation, which is very useful for planning, monitoring, and mitigation of hydrological drought in areas having problems related to hydrological data availability. propagation rate varied spatially all over South Korea. Our results are also consistent with these previous ﬁndings in terms of the spatial pattern of propagation.


Introduction
Drought is one of the most destructive natural hazards because it has negative impacts on the environment, agriculture, economics, and society, and occurs in most climatic zones over the world [1][2][3]. Many efforts have been made worldwide in the planning, monitoring, and mitigation of drought due to its impacts and complex nature. Droughts can be generally classified into four closely related categories: meteorological, agricultural, hydrological, and socioeconomic [4]. Among these drought types, meteorological drought is mainly dependent on deficient precipitation, while agricultural The objectives of this work are (1) to find the propagation probability of all the drought classes (i.e., meteorological to hydrological), which will give us an understanding of the propagation percentage of the different classes in the study area; (2) to probabilistically assess different classes of meteorological and hydrological drought using SPI and SRI, respectively; and (3) to identify the propagation of drought classes using a Markov Bayesian classifier.

Research Area and Data
South Korea ranks high among countries facing water scarcity problems and experiencing serious droughts since the late 1990s [40]. There are five large river basins in South Korea, namely, Han River, Nakdong River, Geum River, Seomjin River, and Yeongsan River. The Han River provides main sources of water for drinking, industry, irrigation, and hydropower generation to the capital and metropolitan areas, as shown in Figure 1. The Han River basin is further subdivided into 24 different sub-basins; the identification number, elevation, and stream network of these sub-basins are shown in Figure 1. For our study purpose, we selected four sub-basins marked with red colors in Figure 1 (enlarged figure). Climatic conditions of the selected four sub-basins are almost similar but hydrological conditions may vary depending upon factors like elevation, slope, and urbanization. These attributes of four selected sub-basins are different that may affect the hydrological cycle or hydrological conditions of the area [41]. Therefore, the selection of these sub-basins may provide a diverse spatial pattern of drought propagation. The drainage areas covered by sub-basins #1003, #1010, #1012, and #1018 are 2483.8 km 2 , 1587.4 km 2 , 1852 km 2 , and 1537.2 km 2 , respectively, and average annual precipitation varies from 1186 mm to 1291mm [42].
propagation in South Korea. Jehanzaib et al. [39] also investigated the effect of climate change on drought propagation using multi-model ensemble projections. These studies incorporated the drought propagation features in their work that revealed the importance of propagation in monitoring and mitigating drought.
Considering the limitations on the modeling of drought propagation, the novelty of this study is to employ a hidden Markov model as a dynamic Bayesian classifier for continuous assessment of meteorological and hydrological drought classes taking the propagation of drought classes. The objectives of this work are (1) to find the propagation probability of all the drought classes (i.e., meteorological to hydrological), which will give us an understanding of the propagation percentage of the different classes in the study area; (2) to probabilistically assess different classes of meteorological and hydrological drought using SPI and SRI, respectively; and (3) to identify the propagation of drought classes using a Markov Bayesian classifier.

Research Area and Data
South Korea ranks high among countries facing water scarcity problems and experiencing serious droughts since the late 1990s [40]. There are five large river basins in South Korea, namely, Han River, Nakdong River, Geum River, Seomjin River, and Yeongsan River. The Han River provides main sources of water for drinking, industry, irrigation, and hydropower generation to the capital and metropolitan areas, as shown in Figure 1. The Han River basin is further subdivided into 24 different sub-basins; the identification number, elevation, and stream network of these sub-basins are shown in Figure 1. For our study purpose, we selected four sub-basins marked with red colors in Figure 1 (enlarged figure). Climatic conditions of the selected four sub-basins are almost similar but hydrological conditions may vary depending upon factors like elevation, slope, and urbanization. These attributes of four selected sub-basins are different that may affect the hydrological cycle or hydrological conditions of the area [41]. Therefore, the selection of these sub-basins may provide a diverse spatial pattern of drought propagation. The drainage areas covered by sub-basins #1003, #1010, #1012, and #1018 are 2483.8 km 2 , 1587.4 km 2 , 1852 km 2 , and 1537.2 km 2 , respectively, and average annual precipitation varies from 1186 mm to 1291mm [42]. Precipitation and runoff data of the period 1967 to 2013 were used for investigating the meteorological and hydrological drought characteristics. Precipitation data were obtained from the Korea Meteorological Administration (https://web.kma.go.kr) for all meteorological stations within the selected sub-basins. For runoff data, the government authorities installed several gauges but firstly the data available was station-based and we needed basin-based data; secondly, the data available was discrete data with a limited period of time, which was not enough to compute the drought characteristics accurately and to train our model. Therefore, different simulated models were used for estimation of runoff for the sub-basins and the most famous and widely adopted models are the PRMS and TANK models. In our study, runoff data were generated for the same period as the precipitation data by using the TANK model, which is the most widely used model in South Korea for rainfall-runoff modelling [43]. Many studies verified the application of the TANK model for runoff simulation. For example, Song et al. [44] proved the accuracy of the TANK model in simulation of runoff at ungauged watersheds in South Korea, and Jaiswal et al. [45] evaluated the accuracy of the TANK model in simulating runoff in Australian watersheds.

Calculation of Standardized Drought Indices
The SPI is a very popular and commonly used drought index to quantify meteorological drought, which is recommended by the World Meteorological Organization (WMO). The fundamental advantage of SPI is its simplicity and ability to calculate drought severity at multiple time scales. Due to the multi-scale properties, the SPI can be used to monitor short-and long-term rainfall deficits, which can be applied for monitoring agricultural and hydrological drought. To calculate the SPI, daily precipitation data were first aggregated into a specific time scale and then a gamma distribution was fitted to the data, as is accepted in many previous studies [37,46,47]. After that, the data were transformed into a standard normal distribution with a mean of zero and standard deviation of one. The SPI represents the deviation of precipitation values from their mean, and positive values indicate surplus precipitation, while negative values represent deficit precipitation.
Following the development of the SPI, Shukla and Wood [7] proposed the SRI to characterize and monitor hydrological drought using runoff data. Similar to the SPI, the SRI is calculated by fitting a lognormal distribution [37,47], and is transformed into a standard normal variate with a mean of zero and standard deviation of one. Positive values of SRI indicate surplus runoff from the mean value, while negative values indicate deficit or drought conditions. Classification of wet and drought conditions based on the SPI and SRI is given in Table 1 [46].

Estimation of Propagation Probability
The concept of propagation was used in this study to probabilistically assess the different classes of drought. Since hydrological drought is a slower developing phenomenon, a propagation probability matrix provides a clear visualization of the relationships of different classes of different drought types at a given time step. To estimate the propagation probability of different classes of droughts in this study, we employed a conditional probability approach to estimate the propagation probability of each class defined in Table 1, i.e., {1, 2, 3, 4, 5} from the SPI to SRI. Equation (1) estimated the conditional probability for propagation of each class.
Atmosphere 2020, 11, 1000 5 of 15 where C h and C m are the classes to which the hydrological and meteorological drought belong, respectively, as mentioned in Table 1.

Markov Bayesian Classifier (MBC)
The Markov Bayesian Classifier is a probability classifier based on Bayes' rule, the dynamics of which are governed by a hidden Markov chain. We considered five classes according to Table 1, and a vector Y R n of the data was allocated to either of these classes. Formally, there is a pair Z = (Y, J), where Y is a random vector, and J is a random variable that assigns class information to Y, i.e., J: Ω → {1, 2, 3, 4, 5}, where Ω is the sample space. We observed only Y, where J is hidden. Thus, a rule (function) is needed to assign an observed vector Y as accurately as possible. We assumed at this stage that the mean and covariance matrices, u i and v i of the classes, are known.
The next task is to derive a classifier that assigns a probability to the event that an observation is from class i. Denoting the prior P(J = i) by p(i), the class (posterior) probabilities for the observed y can be calculated by Bayes' rule.
where P(y|J = i) is the conditional probability of class i. Further, when P(J = i|Y = y) is known, the Bayesian classifier allocates y to r, g(y) = r. Class probabilities express uncertainty: the closer the probability estimate is to one or to zero, the less uncertainty there is in a decision. The classifier defined in Equation (2) could be static. However, our aim was to form a dynamic probability classifier that also describes the dependence between observations when dealing with time series data, Z t = (Y t , J t ).
The normality assumption of the classes now takes the form of Y t~N (u jt , v jt ).
To model time dependence in a mathematically tractable way, we postulated that J t could be described as a homogenous Markov chain, where the data generating process has five hidden classes. For each class, the likelihood of various observations is one of the five multi-normal densities. A Markov chain generates switching between classes. When in class i, the process is said to be working in regime i. These Markov probabilities are collected into the Markov matrix P = p ij . The model is now defined by the regime distribution and the Markov matrix. When class priors at t−1, P i (t−1) are known, the regime (posterior) probabilities for a given y t are Equations (2) and (3) define the Markov Bayesian Classifier (MBC). The introduced model is a special case of a hidden Markov model (HMM), and Y t is an observed time series depending on an unobserved Markov chain J t . We collected the data from 1967-2013; the data from 1967-2000 were used for training the model, and the data from 2001-2013 were used for application of the model. The parameters of the classes were estimated using the attributes of the SPI and SRI, which were further used for classification. By following the propagation concept, the probabilistic assessment of each class of drought based on SPI and SRI of current month t depends on the probability value of the SPI drought classes at previous month t−1. For example, the probability of occurrence of different classes of drought based on SPI and SRI in March depends on the probability of the SPI classes in February. The probability value for each class can be obtained by iterating Equation (4). 15 . . . . . . . . .
Atmosphere 2020, 11, 1000 6 of 15 where X = [J, K] represents both the SPI and SRI drought classes. P(X t+1 = i|Y t ) is the probability of both the SPI and SRI drought classes in the next time step given the probability of the SPI drought classes in the current time step J t . P ij is the transition probability matrix for both the SPI and SRI classes in the next time step conditioned on the SPI classes in the current time step.

Relationship between SPI and SRI
In this study, we selected a three-month time scale for both SPI and SRI to reflect the seasonal characteristics of both drought types. To demonstrate that the time scale selection was correct, a Pearson correlation analysis was performed [35]. The estimation of the correlation coefficient between the two drought indices was very important and beneficial for the setup of MBC. Because we are going to analyze the propagation of drought classes from SPI to SRI, if the SPI and SRI were poorly correlated with each other, the model results may not be accurate [48]. Sattar and Kim [37] evaluated the relationship between SPI/SRI and SPEI/SRI for the Han River basin and found high values of the correlation coefficient between them. Similarly, Sattar et al. [47] also found high values of the correlation coefficient between the SPI and the SRI for historical and future periods under the Representative Concentration Pathway (RCP) 8.5 climate change scenario. It is clear from Figure 2 that, in our study, the correlation coefficient between the SPI and SRI is significant in all the selected sub-basins, varying from 0.77 to 0.88.
Atmosphere 2020, 11, x FOR PEER REVIEW 6 of 16 where X = [J, K] represents both the SPI and SRI drought classes. P(Xt+1 = i|Yt) is the probability of both the SPI and SRI drought classes in the next time step given the probability of the SPI drought classes in the current time step Jt. Pij is the transition probability matrix for both the SPI and SRI classes in the next time step conditioned on the SPI classes in the current time step.

Relationship between SPI and SRI
In this study, we selected a three-month time scale for both SPI and SRI to reflect the seasonal characteristics of both drought types. To demonstrate that the time scale selection was correct, a Pearson correlation analysis was performed [35]. The estimation of the correlation coefficient between the two drought indices was very important and beneficial for the setup of MBC. Because we are going to analyze the propagation of drought classes from SPI to SRI, if the SPI and SRI were poorly correlated with each other, the model results may not be accurate [48]. Sattar and Kim [37] evaluated the relationship between SPI/SRI and SPEI/SRI for the Han River basin and found high values of the correlation coefficient between them. Similarly, Sattar et al. [47] also found high values of the correlation coefficient between the SPI and the SRI for historical and future periods under the Representative Concentration Pathway (RCP) 8.5 climate change scenario. It is clear from Figure 2 that, in our study, the correlation coefficient between the SPI and SRI is significant in all the selected sub-basins, varying from 0.77 to 0.88. Based on the criterion used to define drought classes in Table 1, we compared the SPI and SRI to investigate the propagation relationship between them. Figure 3 clearly illustrates temporal evolution of SPI and SRI and provides a symbolic example of propagation of different classes from the SPI to SRI in the sub-basin #1003. The different combination of colors represents different classes. This comparative analysis between the SPI and SRI showed a very good explanation of the drought classes propagation features. For example, some marked areas of Figure 3 illustrate that during 1968, 1971, 1973-1975, Atmosphere 2020, 11, 1000 7 of 15 1977, 1981-1983, 1988, and 1996, moderate, severe, and extreme classes of SPI were propagated into moderate, severe, and extreme classes of SRI. In comparison to the previous work on propagation in South Korea, Shin et al. [20] also found propagation of meteorological to hydrological drought at different time scales and concluded the propagated events in the year 1977, 1980-1982, 1988, 1992, 1994, 1996, 2001, and 2014-2015.
Based on the criterion used to define drought classes in Table 1, we compared the SPI and SRI to investigate the propagation relationship between them. Figure 3 clearly illustrates temporal evolution of SPI and SRI and provides a symbolic example of propagation of different classes from the SPI to SRI in the sub-basin #1003. The different combination of colors represents different classes. This comparative analysis between the SPI and SRI showed a very good explanation of the drought classes propagation features. For example, some marked areas of Figure 3 illustrate that during 1968, 1971, 1973-1975, 1977, 1981-1983, 1988, and 1996, moderate, severe, and extreme classes of SPI were propagated into moderate, severe, and extreme classes of SRI. In comparison to the previous work on propagation in South Korea, Shin et al. [20] also found propagation of meteorological to hydrological drought at different time scales and concluded the propagated events in the year 1977, 1980-1982, 1988, 1992, 1994, 1996, 2001, and 2014-2015.

Propagation Probability of Drought Classes
We understood from Figure 3 that there was a strong linkage between the SPI and SRI in terms of propagation of different classes of meteorological to hydrological drought, so quantification of this linkage in terms of probability of each class from the SPI to SRI is necessary for further analysis. The propagation probability of each class for the selected sub-basins based on Equation (1) is shown in Figure 4. From Figure 4, it can be seen that when all the sub-basins have extremely wet classes of SPI, they have high probabilities that the SRI also has extremely wet classes and lower probabilities to have moderate wet classes. If the SPI has a moderate wet class, then there is a high probability that the SRI also has a moderate wet class and lower probability to have a near normal class in all subbasins. Likewise, when there is near normal class for SPI, there is high probability that the SRI may have a near normal class and lower probability to have moderate wet and moderate drought class for all sub-basins. In addition, when the SPI has a moderate drought class, there is a high probability that the SRI will have a near normal class and a lower probability to transform into a moderate as well as severe and extreme class of drought for all sub-basins. Lastly, for the SPI to have a severe and extreme drought class, then there is a high probability that the SRI will have a severe and extreme drought class in sub-basin #1003; in sub-basin #1010 it has an equal probability to have a near normal class as well as severe and extreme class of SRI; and in sub-basin #1012 and #1018 it has a high probability to propagate into a moderate drought class of SRI. From Figure 4, it is clear that the propagation probability values of all classes from SPI to SRI have different values in all sub-basins.

Propagation Probability of Drought Classes
We understood from Figure 3 that there was a strong linkage between the SPI and SRI in terms of propagation of different classes of meteorological to hydrological drought, so quantification of this linkage in terms of probability of each class from the SPI to SRI is necessary for further analysis. The propagation probability of each class for the selected sub-basins based on Equation (1) is shown in Figure 4. From Figure 4, it can be seen that when all the sub-basins have extremely wet classes of SPI, they have high probabilities that the SRI also has extremely wet classes and lower probabilities to have moderate wet classes. If the SPI has a moderate wet class, then there is a high probability that the SRI also has a moderate wet class and lower probability to have a near normal class in all sub-basins. Likewise, when there is near normal class for SPI, there is high probability that the SRI may have a near normal class and lower probability to have moderate wet and moderate drought class for all sub-basins. In addition, when the SPI has a moderate drought class, there is a high probability that the SRI will have a near normal class and a lower probability to transform into a moderate as well as severe and extreme class of drought for all sub-basins. Lastly, for the SPI to have a severe and extreme drought class, then there is a high probability that the SRI will have a severe and extreme drought class in sub-basin #1003; in sub-basin #1010 it has an equal probability to have a near normal class as well as severe and extreme class of SRI; and in sub-basin #1012 and #1018 it has a high probability to propagate into a moderate drought class of SRI. From Figure 4, it is clear that the propagation probability values of all classes from SPI to SRI have different values in all sub-basins.
The variation in these propagation probability values and drought classes was due to the complex nature of hydrological cycle. From Figure 1, it can be seen that the elevation and stream network is different for all selected sub-basins and, furthermore, the vegetation, land-use, and land-cover characteristics of these sub-basins are also different. It is confirmed from previous studies that soil moisture conditions, topography, vegetation cover, and groundwater level of the catchment area affect the transformation of meteorological conditions [34]. These factors can cause a different propagation probability and different condition of drought classes. Sattar et al. [38] found that the propagation rate varied spatially in the Han River basin for historical and future periods. Sattar et al. [47] also concluded that the propagation rate varied spatially all over South Korea. Our results are also consistent with these previous findings in terms of the spatial pattern of propagation. soil moisture conditions, topography, vegetation cover, and groundwater level of the catchment area affect the transformation of meteorological conditions [34]. These factors can cause a different propagation probability and different condition of drought classes. Sattar et al. [38] found that the propagation rate varied spatially in the Han River basin for historical and future periods. Sattar et al. [47] also concluded that the propagation rate varied spatially all over South Korea. Our results are also consistent with these previous findings in terms of the spatial pattern of propagation.

Assessment of Drought Classes
In this section, we demonstrated the MBC results for assessing the drought classes of SPI and SRI based on propagation. Equations (2) and (3) provide a detailed description of the parameter values and the probability of each class for the training period of 1967-2000. The nature of a dynamic Markov chain provides a reasonable description of the transition probability matrix for each month of the year, yielding twelve transition probability evaluations. We assumed here that the current time is t and class values are known at every point in the past up to t. When we now try to look from the past into the future, the situation turns into a genuine HMM problem. Moreover, the MBC assumes that the dynamic process follows the first-order Markov chain process, i.e., each class of the month is

Assessment of Drought Classes
In this section, we demonstrated the MBC results for assessing the drought classes of SPI and SRI based on propagation. Equations (2) and (3) provide a detailed description of the parameter values and the probability of each class for the training period of 1967-2000. The nature of a dynamic Markov chain provides a reasonable description of the transition probability matrix for each month of the year, yielding twelve transition probability evaluations. We assumed here that the current time is t and class values are known at every point in the past up to t. When we now try to look from the past into the future, the situation turns into a genuine HMM problem. Moreover, the MBC assumes that the dynamic process follows the first-order Markov chain process, i.e., each class of the month is dependent on the previous class of the month, and the transition probability of each class is independent of time in a dynamic process.
The results of the MBC model for classification of the SPI at X t+1 based on the probability of classes in each previous month J t for the application period 2001-2013 are shown in Figure 5 We further employed the MBC to assess classes of SRI at Xt+1 based on the propagation probability of the SPI classes in a previous month Jt for the application period 2001-2013. The results are shown in Figure 6. The MBC for a given time step was set as the most likely drought class with a high probability value. Similar to the above evaluation, we compared the classification results of the MBC with the observed SRI, specifically for dry classes of drought.  We further employed the MBC to assess classes of SRI at X t+1 based on the propagation probability of the SPI classes in a previous month J t for the application period 2001-2013. The results are shown in Figure 6. The MBC for a given time step was set as the most likely drought class with a high probability value. Similar to the above evaluation, we compared the classification results of the MBC with the observed SRI, specifically for dry classes of drought. In sub-basin #1003, dry classes were observed

Performance Evaluation of the MBC
We employed confusion matrices as shown in Table 2 to evaluate the probabilistic classification results. A confusion matrix gives a better idea of the accuracy of the MBC in predicting the frequency of each class in terms of matches to the frequency of actual classes determined by the observed SPI From previous studies, the drought events observed in the Han River basin were in different months of the year 2004, 2006, 2007, 2008, 2009, 2012, and 2013 [21,41,49]. The comparison results of the MBC with the observed SPI, SRI, and previous studies, revealed that classification results of the MBC match most of the time with the observed SPI, SRI, and the previous findings in predicting the classes. However, in some years, the MBC failed to predict the expected classes of drought, and another shortcoming is that the start and end time of the classes also vary in comparison to SPI and SRI. However, the SPI and SRI provided a probability of one for an optimal class and probability of zero for all other classes while the MBC was classified with the associated probability for each class, which is indicated by the height of the bars in different colors. We also investigated spatial variation in probability values of each class in all sub-basins. For example, Figure 6a shows a Class 5 observed with a high probability value in sub-basin #1003 in the years 2011-2012, while the Class 5 values were rare or a low probability value in the other sub-basins.

Performance Evaluation of the MBC
We employed confusion matrices as shown in Table 2 to evaluate the probabilistic classification results. A confusion matrix gives a better idea of the accuracy of the MBC in predicting the frequency of each class in terms of matches to the frequency of actual classes determined by the observed SPI and SRI. The frequency of occurrence of each class for the MBC was compared with the actual classes of the observed SPI and SRI. Table 2 shows that the MBC presented the most matches with the observed SPI and SRI in determining all classes. Further, to quantify this confusion matrix result, we estimated the accuracy and precision from the confusion matrix for each class. The accuracy and precision of the classification can be determined by taking the ratio of the frequency of the predicted class to the frequency of the actual class, and by taking the ratio of the frequency of the predicted class to the sum of the frequency of the predicted classes, respectively. Table 2. The frequency of each drought class classified by the MBC and standardized drought indices.  Tables 3 and 4 shows that the accuracy and precision results are quite satisfactory for classification of both the SPI and SRI drought classes. In the case of SPI, the accuracy estimated for Class 1 varied from 52% to 62%; for Class 2 from 60% to 68%; for Class 3 from 63 % to 76%; for Class 4 from 50% to 58%; and for Class 5 from 36% to 50%. Similarly, in comparison with the SRI, the accuracy in estimating classes by the MBC varied for Class 1 from 50 % to 63%; for Class 2 from 46% to 70%; for Class 3 from 59% to 67%; for Class 4 from 45% to 67%; and for Class 5 from 33% to 50%. Generally, from Tables 3  and 4, the accuracy and precision results indicated that the model shows fair performance in classifying all classes in all the sub-basins. In comparison to previous studies, which only focused in characterizing individual drought types and to those that utilized the propagation concept to compare the features of two drought types, our results are a step forward from those findings, as we utilized these features to predict the classes of hydrological drought based on the meteorological drought classes. Therefore, the overall results of our study specify that the MBC was reasonable in predicting the onset and end of drought events, and can be utilized for monitoring and management of short-term drought risk.  Table 4. Accuracy and precision of the MBC in the classification of hydrological drought classes. (%)   1  2  3  4  5  1  2  3  4  5   1003  63  52  59  45  40  58  85  46  26  50   1010  61  70  62  67  50  57  56  78  50  50   1012  55  46  67  57  33  53  64  62  27  33   1018  50  51  59  58  42  45  55  73 24 71

Conclusions
In this study, we utilized the concept of propagation of drought classes to probabilistically assess the meteorological and hydrological drought classes represented by the standardized precipitation index (SPI) and standardized runoff index (SRI), respectively. We proposed an iteration of feature extraction, classification, and application for assessment of both the SPI and SRI drought classes. The observed SPI and SRI were used to extract the different features of the drought classes, and classification was performed for each month using a Bayesian model. The dynamic model Markov chain model provided reasonable descriptions of the meteorological and hydrological drought classes. We referred to the proposed model as a Markov Bayesian Classifier (MBC). Overall, the results revealed high probability values of the propagation percentage of the different classes from SPI to SRI, as are further employed in the MBC. The classification results of the MBC were compared with the observed SPI, SRI, and previous studies, demonstrating that the results of the MBC match the observed SPI, SRI, and previous findings most of the time when predicting the classes. However, in some years, MBC failed to predict the expected classes of drought, and another shortcoming is that the start and end time of the classes also vary in comparison to SPI and SRI. The accuracy of the MBC model in predicting all the classes of meteorological drought varies from 36 to 76% and in predicting all the classes of hydrological drought varies from 33 to 70%. Similarly, the precision of the MBC model in predicting all the classes of meteorological drought varies from 42 to 86% and in predicting all the classes of hydrological drought varies from 24 to 85%. However, the SPI and SRI provided a probability of one for an optimal class and probability of zero for all other classes, while the MBC was classified with the associated probability for each class. The MBC-based drought assessment model can easily be generalized to other study areas. The benefit of this model is that it considers the drought propagation, which is beneficial for