Abstract

Drought is a creeping phenomenon that slowly holds an area over time and can be continued for many years. The impacts of drought occurrences can affect communities and environments worldwide in several ways. Thus, assessment and monitoring of drought occurrences in a region are crucial for reducing its vulnerability to the negative impacts of drought. Therefore, comprehensive drought assessment techniques and methods are required to develop adaptive strategies that a region can undertake to reduce its vulnerability to drought substantially. For this purpose, this study proposes a new method known as a regional comprehensive assessment of meteorological drought (RCAMD). The Standardized Precipitation Index (SPI), Standardized Precipitation Evapotranspiration Index (SPEI), and Standardized Precipitation and Temperature Index (SPTI) are jointly used for the development of the RCAMD. Further, the RCAMD employs Monte Carlo feature selection (MCFS) and steady-state probabilities (SSPs) to comprehensively collect information from various stations and drought indices. Moreover, the RCAMD is validated on the six selected stations in the northern areas of Pakistan. The outcomes associated with the RCAMD provide a comprehensive regional assessment of meteorological drought and become the initial source for bringing more considerations to drought monitoring and early warning systems.

1. Introduction

Drought is a multifaceted phenomenon triggered by a deficiency of precipitation, and its related impacts have severe effects on weather-related events, natural ecosystems, forestry, economy, agriculture, and environment [15]. It progressively holds an area over time, can be persisted for a long time, and distressed agricultural [68], environmental [911], and socioeconomic conditions [1214]. Furthermore, it exhibits substantial spatial and temporal variability in various climates and regions. Several authors have proposed various procedures and frameworks to address the spatial and temporal variability of drought events [1520]. However, it is considered a highly variable complicated phenomenon, and it is challenging to discover its onset and termination periods [2124]. The complication in drought assessment and monitoring underpins the need for new drought assessment and monitoring methods and procedures [2528].

Wilhite and Glantz [1] categorized the drought into several categories, i.e., “meteorological, agricultural, hydrological, and socioeconomic.” Yihdego et al. [29] have defined meteorological drought as a prolonged precipitation deficit over time. The precipitation data have been used as a single input variable to mark meteorological drought occurrences and onsets [25, 26, 28, 3033]. The continuous shortfall in precipitation interlinks the meteorological drought to the agricultural drought. The agricultural drought manifests itself as a deficiency in precipitation, a deficit in soil moisture condition, crop failure, etc. [34, 35]. Further, the prolonged period without rainfall becomes the root of the hydrological drought [36, 37]. Hydrological drought manifests itself as decreased streamflow eviction and falling water level in lakes, groundwater, or reservoirs [38]. The hydrological drought can be damaging and cause severe societal impacts if not alleviated timely. The drought of socioeconomic concerns the supply and demands of the economic goods and is associated with the other three types of drought [39]. An extended period with a deficit in precipitation leads to crop failure issues, a shortage of water supply, and industrial and economic productivity [40]. Increasing demand for goods can lead to exploitation, resulting in vast socioeconomic influences and conflicts. In the recent past, drought has become one of the most dangerous natural hazards and disturbed economic and environmental sectors worldwide [4144].

Distinctively, drought has been assessed under meteorological, agricultural, hydrological, and socioeconomic aspects by developing various indices that have been discussed and employed in various publications [4558]. The indices are essential components for assessing and monitoring drought since they simply quantify the complicated interrelationships between varying climate and climate-related parameters [5963]. Wilhite et al. (2000) have defined that indices are developed to communicate information related to climate anomalies to diverse users and allow researchers to evaluate climate anomalies quantitatively in terms of spatial extent. Several drought indices are established and employed to quantitatively assess the impacts of several kinds of droughts to provide helpful information for planning, organizing, and various management applications of water resources associated with several users and the environment [4551, 6467].

Along with the numerous indices proposed for assessing the meteorological drought, some specific indices are extensively used. In particular, Palmer’s Drought Severity Index (PDSI) was presented and used [68, 69]. The index was created to “ measure the cumulative departure of moisture supply.” The PDSI is commonly used by the United States (USA). Further, instead of precipitation variability, the PSDI expands its assurance of drought on water supply and demand. The PDSI comprises important determinants, including data on soil temperature and precipitation. By incorporating these determinants as inputs, the PDSI analyzes four terms in the water balance equation (“evapotranspiration, moisture, soil recharge, and runoff”). Another extensively used index for the characterization of the meteorological drought is the Standardized Precipitation Index (SPI) [46, 7073]. The SPI comprises only a single determinant, which is precipitation, and thus, SPI uses precipitation as an input to describe the water deficit. SPI is a renowned index, extensively used to assess and monitor meteorological drought. The SPI is less complicated than the PDSI. Therefore, it can be applied in any place by transforming the precipitation data from a skewed distribution to a normal distribution. Moreover, SPI with longer time scales can indicate the agricultural and hydrological drought [71, 74, 75]. For instance, the SPI for a nine-month time scale with a value less than −1.5 is an alert for the agricultural drought [59]. The streamflow, reservoir level, etc., can be reflected by positioning SPI at a twelve-month time scale. Therefore, SPI is famous and operational in numerous papers and publications [51, 76]. Further, the Standardized Precipitation Evapotranspiration Index (SPEI) is also a well-known index proposed by [50] that triggers the effect of temperature variability on drought estimation. Numerous analyses have employed SPEI for the drought evaluation [7783], and Standardized Precipitation and Temperature Index (SPTI) by [84] is also considered in multiple studies for the assessment of meteorological drought [76, 8486].

Considerable research has been done to quantify and understand the complex and meteoric nature of the drought [77, 78, 80, 81, 83, 8790]. However, the manifestation of the drought nature is very complex [91]. The complexity of determining its pattern reinforces the development of new techniques and methods [92, 93]. The appropriate methods and procedures can help to minimize its meteoric influence in various parts of the world [87] [94] [90, 95]. However, the applications of the new methods may be better described by investigating drought at the regional level. Recently, numerous studies have been done to timely examine the drought occurrences in various regions. Therefore, the study of the particular region has significant importance; thus, current research is applied to the specific region. The selected region has a homogeneous pattern of drought occurrences concerning specific drought indices and a time scale (one-month time scale) [76, 9699]. Ali et al. [96] examined meteorological drought based on three indices (SPI, SPEI, and SPTI). The study found that the three indices provide similar information about the selected region for the particular time scale. Hence, investigating meteorological drought from the selected homogeneous locations using several meteorological indices (SPI, SPEI, and SPTI) becomes counterproductive. This issue underpins the use of some new drought assessment methods that provide comprehensive information based on these indices. Therefore, this study proposes a new method, known as regional comprehensive assessment for meteorological drought (RCAMD). The RCAMD comprehensively collects information from several stations and drought indices using Monte Carlo feature selection (MCFS) and steady-state probabilities (SSPs). Further, the RCAMD mainly helps to overcome two issues. For instance, the first phase of the RCAMD chooses important stations more comprehensively for three indices from six homogeneous stations. In the presence of influential climatic factors in estimating the drought indices, the second phase of RCAMD characterizes several drought classes more comprehensively and accurately among the three indices (SPI, SPEI, and SPTI). Moreover, the six stations in the northern areas of Pakistan are selected to validate RCAMD. The findings associated with the RCAMD propose a comprehensive regional assessment of meteorological drought and create the initial basis for taking more considerations for assessing and monitoring drought at the regional level.

2. Material and Methods

2.1. Description of the Study Area

The substantial climate changes have increasingly become a primary global task that endangers ecological, human, and natural systems [100102]. Pakistan is extremely in danger of the undesirable influences of climate change, specifically extreme hydrometeorological activities [103108]. The selected region is located in the northeastern part of Pakistan, spread over 72,971 square km, almost half of which covers peaks of mountains, glaciers, highlands, and lakes. The selected region has structural significance for other parts of the country. It has a key role in the agricultural sectors and the reservoir system of the country [109, 110]. However, it is highly at risk of climate change due to its geological composition, fragile mountain, topography, ecosystem, geographic locations, socioeconomic conditions, and scattered population [111]. Thus, the selected region requires more consideration for assessing the drought manifestations by developing comprehensive and proficient methods and procedures. Hence, the RCAMD method is designed for the selected region, enhancing the ability to assess drought events and facilitating drought monitoring and water resource management in the selected area (Figure 1).

2.2. Data and Methods

The data ranging from January 1971 to December 2017 are processed in the current analysis. The six stations in the northern areas are selected to calculate the indices (SPI, SPEI, and SPTI). These indices use information from the indicators (precipitation and temperature) to classify drought classes in the selected stations. The data of these indicators have been used in several publications [86, 97, 98, 112114]. The various drought classes of the selected stations and indices are used to propose RCAMD. The new proposed RCAMD uses MCFS and SSP to assess the information more intensively for the drought classes from the selected stations and indices. Moreover, the outcomes associated with RCAMD comprehensively assess drought classes for the homogeneous region. The RCAMD develops a new way of taking more consideration for evaluating and monitoring drought at the regional level.

2.2.1. Monte Carlo Feature Selection (MCFS)

Niaz et al. [85] used MCFS for selecting informative stations for their analysis. They applied MCFS in the Punjab region of Pakistan and selected important stations based on SPTI. However, in this study, MCFS uses three drought indices (SPI, SPEI, and SPTI) to select important stations. Thus, the MCFS selects an important meteorological station for each drought index. For example, for SPI the MCFS selects Astore as an important station; for SPEI, the MCFS selects Gilgit as an important station. Further, for SPTI, the MCFS selects Astore as an important station for the preliminary analysis. The MCFS input enables the RCAMD to collect information from various stations comprehensively. The suitable stations are chosen based on relative importance (RI) values. The mathematical detail about the RI is given in [85]. For the current analysis, the Astore Station with RI value of 0.1385 is higher than other selected stations for SPI. For SPTI and SPEI, the Astore and Gilgit are selected, respectively. In Astore, the RI value for SPTI is 0.1920, while SPEI for Gilgit has RI value of 0.7617.

2.2.2. Steady-State Probabilities (SSPs)

A Markov process can be expressed as the probabilities come up to the SSP when certain periods have been passed. The comprehensive mathematical details associated with the SSP of the Markov chain were described in Stewart (2009). The application of SSP is provided in several publications [76, 97]. Niaz et al. [85] used SSP as a weighting scheme from the long-run time-series data for different drought classes in the northern region of Pakistan. The proposed weighting scheme was used to accumulate information from the selected homogeneous stations. Further, Niaz et al. [97] used SSP to substantiate the prevalence of drought intensities in the northern region of Pakistan. Moreover, Niaz et al. [98] proposed a new technique based on SSP to accumulate information from various indices. Recently, Niaz et al. [99] incorporated SSP in their study to assess the probability of drought severity in the selected region. The SSP is used broadly in several publications to develop new methods and procedures [9799, 113]. Therefore, in RCAMD, SSP is used to propagate weights for various drought categories over several stations and indices to achieve a particular aspect. In the current analysis, SSP mainly helps to characterize the new vector of drought classes. The inclusion of MCFS and SSP in RCAMD makes the study innovative. This innovation provides a comprehensive procedure to collect information from several stations and indices.

2.2.3. Regional Comprehensive Assessment of Meteorological Drought (RCAMD)

The RCAMD employs MCFS and SSP to mainly determine drought events that are likely to occur in the region from numerous stations and drought indices. The MCFS technique is used to accumulate comprehensive information on several time-series data of meteorological stations. The mathematical detail of MCFS is given in Niaz et al. [85]. In the first phase of the RCAMD, the MCFS allows the selection of more important stations based on several selected indices. Three drought indices (SPI, SPEI, and SPTI) are used in RCAMD to determine important stations. Hence, the MCFS selects an important meteorological station for each drought index separately. The criteria for selecting an important station are based on relative importance (RI). The higher values corresponding to any stations show that the stations are important for the preliminary investigation. For example, based on the higher value of RI using SPEI, the MCFS chooses Gilgit as an important station, and for SPI, the MCFS takes Astore as an important station. Moreover, for SPTI the MCFS picks Astore as an important station for the computation of RCAMD. In the second phase of the RCAMD, the SSP is applied to characterize several drought classes among the three indices (SPI, SPEI, and SPI). The complete mathematical detail related to the SSP of the Markov chain is given in Stewart (2009). The SSP is used in several publications to develop new procedures and methodologies [76, 97]. The SSP characterizes various drought categories among selected stations and indices in this study. The SSP for each drought category ((“Extremely Dry (ED),” Extremely Wet (EW),” “Severely Dry (SD),” “Severely Wet (SW),” Median Dry (MD),” “Median Wet (MW),” and “Normal Dry (ND)”) for each index (SPI, SPEI, and SPTI) in the particular station () can be expressed in a vector as . The obtained SSP for the varying drought categories can be described as the visit of the drought category in the long run. These long-run SSP of several drought categories is counted as weights. These weights are further utilized for the computation of RCAMD. The calculation of RCAMD is based on the vector of the stationary drought categories propagating on different drought indices, which can be identified as follows:

The obtained limiting probabilities , , can be referred to as the proportion or average of long-run probabilities of the drought states or categories for the varying indices (SPI, SPTI, and SPEI) on selected stations. These probabilities are used as weights for the computation of the RCAMD, which assigns the comprehensive weights to the varying drought categories from the selected stations. The flowchart of the RCAMD is given in Figure 2. Moreover, the drought states among selected drought indices (SPI, SPTI, and SPEI) that take maximum weights are chosen for the RCAMD. Hence, in the current research, the RCAMD selects the appropriate vector of drought classes from the time-series dataset for January 1971 to December 2017. The RCAMD enables a clearer, though complicated, representation of how interconnected indices are further associated and linkable to a distinctive set of comprehensive outcomes. Further, RCAMD can be utilized to locate the proper vector of drought classes for any long time-series data in a homogeneous environment. The outcomes associated with the RCAMD provide a comprehensive regional assessment of meteorological drought and become the initial source for bringing more considerations to drought monitoring and early warning systems.

3. Results

The time-series data are collected for six meteorological stations from the northern areas of Pakistan. The varying features, including mean, 1st quartile, median, 3rd quartile, kurtosis, and standard deviation (St.Dev) of the precipitation, are given in Table 1. Table 2 contains the varying characteristics of the minimum temperature. The various features of the maximum temperature are given in Table 3. Further, these climatological features are presented in various figures. For example, the climatological features of the monthly precipitation observed in varying stations are presented on various maps in Figure 3. The climatological characteristics of the observed minimum temperature in various stations are presented in Figure 4. The climatological features of the observed maximum temperature in various stations are presented in Figure 5. Further, the drought categories are classified according to Li et al. [115]. The varying behavior of the classes can be observed in the selected time-series data. However, for simplicity, the results for the specific year (2017) based on SPI, SPEI, and SPTI are provided. In Table 4, the results for the year 2017 based on SPI are given. The varying drought classes can be observed in varying months of the selected year. Further, the index values corresponding to each drought category are provided. Table 5 contains the classified values based on SPEI, and the classified values observed in varying months and their corresponding index values based on SPTI are given in Table 6. The temporal behavior in the selected period, January 1971 to December 2017, in varying stations for the SPI at a one-month time scale (SPI-1) is presented in Figure 6. Figure 7 presents the temporal behavior of the SPEI on a one-month time scale (SPEI-1) at selected stations. Figure 8 shows the temporal behavior of the SPTI at a one-month time scale (SPTI-1) at various stations. Moreover, the varying drought categories observed in various stations for SPI at a one-month time scale are presented in Figure 9. Figure 10 contains various maps for the varying drought categories observed in various stations for SPEI at a one-month time scale. The drought categories classified based on SPTI at a one-month time scale and observed in varying stations are presented in Figure 11. The northern zones of Pakistan (i.e., Astore, Bunji, Chilas, Gupis, Skardu, and Gilgit) have a homogeneous pattern of the drought classes for the specific drought indices and time scale [76, 96, 97, 99] and therefore selected for the current analysis. Three indices, SPI, SPEI, and SPTI, have shown a significant correlation and provided similar information in varying stations at a one-month time scale [76, 96, 116]. However, this study found a gap in the above research and proposed a new method that provides more comprehensive results. The mentioned research considered all stations for their analysis. Thus, considering all stations for drought analysis in a region with a similar pattern of drought occurrences seems counterproductive. It underpins a new gap that should be tackled by comprehensively accumulating information. Based on this gap, this study proposed to provide a more comprehensive drought assessment procedure for the region. For this purpose, the current research comprehensively offers an RCAMD method to accumulate information from numerous stations and drought indices. The RCAMD is based on two phases. In the first phase of the RCAMD, the MCFS technique is applied. The MCFS was utilized by Niaz et al. [85] for selecting more illustrative stations in the region of Punjab in Pakistan. The mathematical detail of the MCFS is available in [85].

Further, the selected stations for the current analysis have shown a similar pattern for all stations. Therefore, it underpins a rationale to apply MCFS for selecting only important stations for the analysis. Thus, the MCFS is applied in the current analysis to select more important stations among the selected stations for various drought indices. This selection of the important stations is based on the relative importance (RI) values (Figure 12). The corresponding higher values of RIs in any station show that the station is to consider for the drought assessment. For example, the Astore Station with RI value of 0.1385 for SPI is higher than other selected stations. For SPTI and SPEI, the Astore and Gilgit are selected, respectively. In Astore, the RI value for SPTI is 0.1920, while SPEI for Bunji has RI value of 0.7617 (Table 7).

Moreover, in the presence of influential climatic factors in estimating the drought indices, the second phase of RCAMD comprehensively characterizes numerous drought categories among the selected indices (SPI, SPEI, and SPI) (Figure 13). Niaz et al. [76] proposed a method based on a steady-state weighting scheme. They selected the classes from various stations based on the maximum weights; hence, the classes that received maximum weights among different stations were selected for the analysis. The weights from three indices (SPI, SPEI, and SPTI) for the varying drought categories for a specific year, 2017, are provided in Tables 810, respectively. Recently, Niaz et al. [97] proposed a weighting scheme based on steady-state probabilities for selecting classes among the three indices (SPI, SPEI, and SPI). These indices are correlated for a one-month time scale and present similar information for the six stations in the northern areas [85, 96, 117]. The mathematical detail of the weighting scheme is available in [76].

Similarly, based on the mentioned studies, this study uses the SSP as a weighting scheme in the second phase of the RCAMD for selecting varying drought classes. Conclusively, to accomplish a specific task (i.e., characterize drought classes more comprehensively), therefore, in RCAMD, SSP is utilized to disseminate weights for several drought categories over various stations and indices. The use of SSP mainly characterizes the new vector of drought classes. The RCAMD suggests a comprehensive regional method for assessing meteorological drought and developing the base for taking more considerations for evaluating and monitoring drought at the regional level.

4. Discussion

The data with varying features (precipitation, maximum and minimum temperature) are processed for the current analysis. The six stations in the northern areas are designated for data processing. The observed data are sufficient to calculate the varying SDI (SPI, SPEI, and SPTI). These SDIs are used to assess the drought severity in the selected region. The classification criteria are adopted from Li et al. [115] to characterize drought severity for the selected stations. The characterization and monitoring of the drought occurrences are vital components for the management and planning of water resources, mitigation strategies, and the creation of a climate-resilient society [2527, 118120]. Therefore, this study proposes an RCAMD to comprehensively and accurately characterize drought occurrences. The RCAMD employs MCFS and SSP to collect information from several stations and drought indices. The selected stations have a homogeneous pattern of drought occurrences among each other for specific indices. Ali et al. [96] cited these stations as homogeneous regions in their study. They used three SDIs (SPI, SPEI, and SPTI) and found that these stations are more and less similar in a specific time scale for three indices. Recently, Niaz et al. [76] considered these stations as homogeneous and developed a comprehensive index procedure to assimilate information. Further, Niaz et al. [98] considered these stations as homogeneous and proposed a regional-level propagation framework that is used to collect information from various indices. However, this study found a gap in their research; the mentioned research had considered all stations for the analysis, given that those stations were homogeneous. Hence, it was counterproductive to study all stations in a homogeneous environment. It underpins a gap addressed in this current research by accumulating more comprehensive information. The present research proposes a new method, RCAMD, which provides more comprehensive results. In the first phase, the RCAMD employed MCFS to separately provide important stations for each index. For example, there are three indices (SPI, SPEI, and SPTI) and six stations in the current analysis. The MCFS uses SPI for selecting important stations among six selected stations. Then, MCFS uses SPEI to select the important station from six selected stations, and similarly, it employs SPTI for selecting an important station from the selected stations. Hence, three vectors of the observations are computed by MCFS for each index separately in the first phase. In the second phase, using SSP the RCAMD provides comprehensive information about various drought classes among selected indices and stations. Hence, the results related to the RCAMD provide a comprehensive assessment of meteorological drought at the regional level and bring a new method to consider more on drought assessment and monitoring. The RCAMD can efficiently work for early warning and mitigation policies. It can be used to make better management and planning to enhance the capabilities of forecasting procedures to decrease the vulnerability of society to drought and its forgoing impacts.

5. Conclusion

Drought is one of the multifaceted natural hazards that has adverse impacts on the economy, water resources, and other environmental structures worldwide. However, the assessment and analysis of drought are crucial, specifically to sound water resource planning and management at the regional level. Therefore, the assessment and monitoring of drought in a region are thus vital to decrease its vulnerability to negative impacts. Therefore, this study proposes an RCAMD. The RCAMD employs MCFS and SSP to accumulate information from several stations and drought indices comprehensively. The three commonly used SDIs are jointly analyzed for the computation of RCAMD. The RCAMD is performed at the six designated stations in the northern areas of Pakistan. The results related to the RCAMD provide a comprehensive assessment of meteorological drought at the regional level and bring a new method to take more consideration on drought assessment and monitoring. Moreover, the RCAMD considers the initial state, and the transition probabilities are constant by assuming time homogeneous progression; however, it can be considered temporal characteristics to improve drought monitoring efficiency for the selected stations. Further, the results of RCAMD would be entertained for the current scenario and application site; however, it cannot be generalized for other climatic conditions. The climatology conditions of the selected stations will change the outcomes and influence the extrapolations. Moreover, the categorization of the given data from other indices in the selected stations can implicitly be useful to increase the capabilities for drought monitoring.

Data Availability

The data and codes used to prepare the manuscript are available from the corresponding author and can be provided upon request.

Ethical Approval

All procedures followed were in accordance with the ethical standards of the Helsinki Declaration of 1975, as revised in 2000.

All authors voluntarily agreed to participate in this research study and agreed to publication; there is no legal constraint in publishing the data used in the manuscript.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

All authors contributed equally.

Acknowledgments

The author extends his appreciation to the Deanship of Scientific Research at King Khalid University for funding this work under grant number (RGP.2/34/43), received by Mohammed M. Almazah (http://www.kku.edu.sa/).