The increasing rate of occurrence of extreme events (droughts and floods) and their rapid transition magnify the associated socio-economic impacts with respect to those caused by the individual event. Understanding of spatio-temporal evolution of wet–dry events collectively, their characteristics, and the transition (wet to dry and dry to wet) is therefore significant to identify and locate most vulnerable hotspots, providing the basis for the adaptation and mitigation measures. The Upper Jhelum Basin (UJB) in South Asia was selected as a case study, where the relevance of wet–dry events and their transition has not been assessed yet, despite clear evidence of climate change in the region. The standardized precipitation evapotranspiration index (SPEI) at the monthly timescale was applied to detect and characterize wet and dry events for the period 1981–2014. The results of temporal variations in SPEI showed a strong change in basin climatic features associated with El Niño–Southern Oscillation (ENSO) at the end of 1997, with the prevalence of wet and dry events before and after 1997 respectively. The results of spatial analysis show a higher susceptibility of the monsoon-dominated region towards wet events, with more intense events occurring in the eastern part, whereas a higher severity and duration are featured in the southwestern part of the basin. In contrast, the westerlies-dominated region was found to be the hotspot of dry events with higher duration, severity, and intensity. Moreover, the surrounding region of the Himalaya divide line and the monsoon-dominated part of the basin were found to be the hotspots of rapid wet–dry transition events.
There is growing evidence that recent warming is leading to significant alteration in the hydrological cycle, exacerbating extreme weather events in general (Peterson et al., 2012) in many regions of the world. Extreme weather events such as floods and droughts and their rapid successions (recurrent spells) during the past few decades have taken a heavy toll on both life and property. Moreover, such events can have large impacts on water availability, agriculture and food security, power production, and natural ecosystems (He et al., 2019; Sheffield and Wood, 2012). These events are projected to regionally intensify and be more frequent within the context of global warming, underscoring the importance of research on wet–dry extreme weather events collectively. The climate change projections for the Asian continent in the sixth Assessment Report (AR6) of the Intergovernmental Panel on Climate Change (IPCC) reported that during the 21st century South Asia is likely to face more intense and frequent heatwaves and humid heat stress, whereas both annual and summer monsoon precipitation will increase, with enhanced inter-annual variability (medium confidence) (Arias et al., 2022). Various studies at local, basin, national, and regional scales already documented and acknowledged the vulnerability to climate change of that region (He and Sheffield, 2020; Zhao et al., 2020; Visser-Quinn et al., 2019; He et al., 2017).
Typically, wet and dry events are generally considered independently in water resource management and planning. However, these events are inherently interconnected and governed by the same underlying hydrological processes and atmospheric dynamics, which may augment hydro-climatic variability under the influence of climate change (He and Sheffield, 2020). A number of rapid wet–dry events in the last decade acknowledged the relevance of sequences of wet and dry events. For example, California's large-scale flood event in 2017 occurred at the offset of prolonged drought (2011–2016) (He et al., 2017; NOAA National Centers for Environmental Information, 2018). South Carolina observed an abrupt transition (within a week) from drought to flood in September 2015 (He and Sheffield, 2020). Other examples include the successive drought and flood events of 2010–2012 and 2015–2016 in the UK (Parry et al., 2013) and Tasmania, Australia respectively (CSIRO, 2018). Such abrupt flood–drought transitions pose a substantial risk for water management practices, especially for reservoir operation, as a trade-off should be set between short-term flood control and long-term water storage imperatives to satisfy water demand (He and Sheffield, 2020). This has aroused widespread concern in the scientific community to understand the wet–dry interplay under a changing environment.
During the past few decades, significant effort was put forward towards the adoption of a multi-hazard approach (consideration of both types of extreme hydrological conditions at the same time) in developing resilience to climate change. Kourgialas (2021) analysed floods and droughts collectively in the Mediterranean agricultural region and proposed water-saving and flood protection measures for adapting to the inevitable adverse effects of climate change. Visser-Quinn et al. (2019) identified hotspot regions in the UK where a spatio-temporally concurrent increase in the number of flood and drought events was projected. Zhao et al. (2020) investigated the rapid transition of flood and drought events under present and future climate change in the Hanjiang Basin and found more frequent drought-to-flood rapid-transition events of higher intensity in the 21st century. Other examples include the analysis of rapid drought-to-flood transitions in river basins in China (Yan et al., 2013) and in England and Wales (Parry et al., 2013). These studies employed the peak-over-threshold (POT) method and various indices recommended by the World Meteorological Organization (WMO) for the detection and characterization of extreme wet and dry events (floods and droughts).
Some commonly used indices are the standardized precipitation index (SPI) (McKee et al., 1993), standardized precipitation evapotranspiration index (SPEI) (Vicente-Serrano et al., 2010), Palmer drought severity index (PDSI) (Palmer, 1965), normalized difference vegetation index (NDVI) (Tucker, 1979), standardized drought indices (SDI) (Svoboda and Fuchs, 2016), and standardized anomaly index (SAI) (Katz and Glantz, 1986). Among these indices, SPI and SPEI are more widely accepted for the following reasons: (a) simple to calculate; (b) require few input data (precipitation and temperature), that are easily accessible in most cases; (c) standardized indices, which facilitate the comparison of different climatic zones; and (d) can be calculated at multiple timescales, depending on the objective. For instance, SPI and SPEI at short timescales (1, 2, 3, or 6 months) better reflect the meteorological and agricultural drought, while longer timescales (12, 24, or 48 months) are usually considered in hydrology (Kourgialas, 2021). The calculation of SPI and SPEI is mathematically similar, but it differs in the input parameters. The SPI only uses precipitation, whereas SPEI is based on the climatic water balance. Many studies advocate the use of SPEI, rather than SPI, due to its link to potential evapotranspiration (PET), which makes it more sensitive in the context of global warming (Himayoun and Roshni, 2019; Yao et al., 2018; Huang et al., 2017; Vicente-Serrano et al., 2010).
In this study, attempts were made to understand the regional evolution of wet–dry events collectively, their characteristics, and their transition (wet to dry and dry to wet) for different severity levels ranging from moderate to extreme. Here, the term “wet and dry events” does not necessarily imply observed flood and drought events, unless explicitly mentioned. There exists a basic difference between a flood and a wet event. The former has a short duration effect (e.g. a few hours or days) while the latter is regarded as a long period without precipitation shortage (e.g. several months or years) (Wu and Chen, 2019).
The proposed framework was implemented with reference to the Upper Jhelum
Basin (UJB), where the relevance of wet–dry events and their transition have
not been assessed yet, despite clear evidence of climate change in the
region. The UJB is located in the western Himalaya and shared by Pakistan and
India. The region already witnessed an increase in extreme
hydro-meteorological events in the last few decades, but these events are
expected to become even more pronounced in the coming future (Pachauri et al., 2014). A study conducted over the northern highlands of Pakistan
investigated the trends in time distribution patterns (TDPs) and return
periods for event-based extreme precipitation for the period of 1961 to 2014
and found maximum values of 20- and 50-year return levels of TDP for the UJB
(Zaman et al., 2020). Another study conducted on a portion of the UJB located in
Kashmir, India, uses the SPEI for spatio-temporal characterization of
drought events only (Himayoun and Roshni, 2019). Akhtar et al.
(2020) investigated the correlation of meteorological and hydrological
drought using the SPEI and the standardized streamflow index (SSI) over the Upper
Indus Basin (UIB), including UJB. They validated the results with a
historically prolonged drought event observed in Pakistan
(1999–2002). Another
study employed the locally weighted SDI and compared it with SPI and SPEI
on 10 meteorological stations within Pakistan (Ali et al., 2019).
Ullah et al. (2021a) evaluated four reanalysis products for drought
assessment in Pakistan using SPI and SPEI at multiple timescales. All above-mentioned studies put a focus on drought event characteristics only,
whereas the wet events and transition of wet–dry events were overlooked.
This study attempts to fill this gap by addressing the following specific
points.
How does climate change influence the evolution of the regional wet–dry
events? How comparatively frequent were wet or dry events in the past? What is the average transition time of wet-to-dry and dry-to-wet events? Which parts of the basin are hosting hotspots for rapid wet–dry transition
events?
The most widely used index, SPEI, is here adopted to detect and characterize
wet and dry events of different severity levels (moderate, severe, and
extreme). The analysis was carried out both at each grid cell and averaged
over the basin, using corrected ERA5 precipitation and observed temperature
data for a period of 35 years (1981–2014).
The Upper Jhelum Basin (UJB) has a latitudinal extent stretching from
Location of the UJB and spatial distribution of climatic stations.
The climate of the UJB is influenced by dynamic local and regional weather
systems, and the topography of the high mountains causes a huge variability
in the spatial and seasonal distribution of precipitation (Dolk et al.,
2020). Two distinct precipitation patterns (i.e. western disturbances and
monsoon) exist in the basin. The western disturbances bring precipitation in
the form of snow during the winter season. The monsoon pattern brings liquid
rainfall during summer seasons. The monsoon precipitation pattern dominates
in the two lower sub-basins, i.e Poonch and Kanshi, and progressively loses
strength northward towards the foothills of the western Himalaya, where the
influence of western disturbances is predominant (Neelum and Kunhar sub-basins). The basin average annual precipitation and temperature is about
1150
The daily observed precipitation and temperature data of 15 climatic
stations located within the political boundary of Pakistan were collected
from the Pakistan Meteorological Department (PMD) and Water and Power
Development Authority (WAPDA). For the Indian region, Indian
Meteorological Department (IMD) daily gridded precipitation and temperature
datasets, derived from a dense network of meteorological stations for the
Indian mainland (Pai et al., 2014), were extracted at five stations and used
for that region. The analysis was carried out for a period of 34 years
(1981–2014), due to the availability of observed data. In fact there are
only a few climatic stations where data are available starting from 1971,
but the number of stations would not be enough for the spatial analysis. The
observed temperature data were used to calculate potential evapotranspiration
(PET) using the Thornthwaite equation (Thornthwaite, 1948) due to data
limitation. A study conducted by Beguería et al. (2014) compared the
SPEI values calculated with three different methods (Penman–Monteith,
Hargreaves, and Thornthwaite) and found small differences in humid regions.
Mavromatis (2007) also reported similar outcomes of PET methods for drought
index calculation. Afterwards PET values were interpolated at 0.25
SPEI, the most widely used index, was adopted to detect and characterize wet
and dry events of different severity levels (moderate, severe, and extreme).
The SPEI supports comparisons over time and space, as proxies of wet
and dry conditions from both the meteorological and agricultural
perspectives. Although the SPEI was originally proposed for drought
monitoring, it can also be used as a tool to detect flood risk. The
calculation procedure of SPEI involves two steps: fitting a log-logistic
distribution to the monthly climatic water balance (P-PET) time series and
then transforming the cumulative probability of the fitted distribution into a
standard normal distribution (with mean zero and variance 1). According to
this distribution method, the probability distribution function of a
variable
SPEI classification of dry and wet events (from Chen et al., 2020).
In this study, three characteristics (severity, duration, and intensity) of wet and dry events were calculated for each pixel. Following Spinoni et al. (2014), the duration (D) of a wet–dry event is the length of time (months) that the index is consecutively above or below a truncation value; the severity (S) refers to the cumulative value of the index from the first month to the last month of the wet–dry event, and it represents the water surplus and deficit respectively; and the intensity (I) of an event is the ratio of severity (S) to duration (D). These characteristics were computed for each event and then further the total wet and dry event duration (TWD and TDD), total wet and dry severity (TWS and TDS), total wet and dry intensity (TWI and TDI), average wet and dry event duration (AWD and ADD), average wet and dry severity (AWS and ADS), average wet and dry intensity (AWI and ADI), maximum wet and dry event duration (MWD and MDD), maximum wet and dry severity (MWS and MDS), and maximum wet and dry intensity (MWI and MDI) were calculated for a period of 34 years (1981–2014).
The wet–dry (WD) ratio is defined as the natural logarithm of the ratio of the
total number of wet months (
The total number of transitions and their average transition time (
The wet–dry rapid-transition event is defined as the consecutive occurrence
of wet and dry months/events. For instance, a dry (or wet) event occurring
in the
The basin average SPEI time series at 1-month (SPEI-1), 3-month (SPEI-3),
6-month (SPEI-6), and 12-month (SPEI-12) timescales is presented in Fig. 2.
It can be seen that the study domain mostly experienced moderate-to-severe
wet–dry events, whereas the extreme wet–dry events (
Temporal variations in SPEI at 1-, 3-, 6-, and 12-month timescales over UJB for the period 1981–2014.
This study focuses on the short-timescale conditions to analyse frequent variations in climatic conditions and their interplay; therefore, more detailed analysis was carried out at the monthly timescale. Moreover, the floods and flash droughts are not clearly associated with long-term SPEI, because the averaging effect of long-term accumulated precipitation and temperature surpasses the signal of extreme precipitation and temperature over a short period. Flash drought is a relatively new type of drought. Currently, there is not a universally accepted definition or criteria for flash drought, though there is general consensus on the principle of rapid onset or intensification characterized by moisture deficits and abnormally high temperatures for a period lasting at least 3 weeks (Lisonbee et al., 2021; Otkin et al., 2018; Hunt et al., 2009). This highlights the usefulness of SPEI at the monthly scale in representing flood and flash drought events. It is noted that the terms “wet–dry events” or “wet–dry months” present similar meaning for our study, as the analysis was made at the monthly time step. A clearer picture of the monthly evolution of wet–dry events of different severity levels and their variability can be seen in Table 2. The SPEI-1 values fluctuate remarkably from one month to another. For example, an extremely wet October in 1987 was followed by a severely dry November, and a severely wet June occurred at the tail of the longest drought spell in May 2001. Such rapid transition from wet to dry and from dry to wet events was more prominent during the first half of the study period (before the year 1997). Another interesting observation concerns the strong change in the basin climatic features which can be noticed around the years 1997–1998. During the first half of the study period (1981–1997), the dominancy of wet events of different categories prevails whereas the basin conditions lean towards dryer conditions during the second half of the period (1998–2014).
Temporal variations in monthly SPEI over UJB from 1981–2014. The blank cells show normal months, and the different severity levels are presented as EW – extremely wet; ED – extremely dry; SW – severely wet; SD –severely dry; MW – moderately wet; and MD – moderately dry). The line between 1997 and 1998 indicates the strong change in the basin climatic features.
Annual variations in the number of months affected by dry–wet events (SPEI
Annual variations in the number of months affected by wet–dry conditions during 1984–2014. The brown and blue colours present dry and wet months respectively. Different shades of the colours define the different severity levels (EW – extreme wet; ED – extreme dry; SW – severe wet; SD – severe dry; MW – moderate wet; MD – moderate dry).
The wet–dry event characteristics (duration, severity, and intensity) were computed for each pixel to analyse their spatial distribution. Pixel-based analysis shows the location of the most vulnerable parts of the basin, providing the basis for future decisions on adaptation and mitigation measures. In this study, the total, average, and maximum values of duration, severity, and intensity were computed for the study period (1981–2014). The maps of wet and dry duration are displayed in Fig. 4. Overall, the study area encountered relatively more wet months than dry months during the whole study period. The total wet duration (TWD) and the total dry duration (TDD) vary from 66 to 80 and from 61 to 65 months respectively for most of the basin. The low-elevation parts in the south of the basin show the highest value of TWD whereas the TDD is higher across the Himalaya divide line than in other parts of the basin. The Himalaya divide line is a line in the middle of the UJB at the Pir Panjal mountainous range, separating the dominance of the two precipitation patterns: westerlies in the north-facing slopes and monsoon in the south-facing slopes of the line (Archer and Fowler, 2008).
Spatial distribution of total wet duration (TWD), total dry duration (TDD), average wet duration (AWD), average dry duration (ADD), maximum wet duration (MWD), and maximum dry duration (MDD) for the period 1981–2014.
The average wet and dry event durations (AWD and ADD) were found to be similar throughout the basin with a slight difference in the range of 1–2 weeks. However, their spatial patterns were found to be mostly complementary. Maximum wet and dry event durations (MWD and MDD) exhibit high values in two distinct parts of the basin. The MWD is about 6–7 months in the east of the basin, which is located in Kashmir, India, whereas it varies between about 4–5 months and 2–3 months in the northwest and southwest parts of the basin. For the MDD, the northwest and central parts of the basin show higher values (4–5 months) than the remaining parts (2–3 months).
The spatial distribution of total, average, and maximum severity of wet–dry events is presented in Fig. 5. All wet–dry severity maps show similar spatial patterns as wet–dry duration maps. In terms of total wet severity (TWS) and total dry severity (TDS), the wet and dry hotspots are located in the south and middle (across Himalaya divide line) of the basin respectively. Unlike the spatial patterns of TDD, the TDS is relatively higher in the north of the basin above the Himalaya divide line. This shows more intense dry events in this part of the basin. The underlying reason for higher TDS could be the higher warming rates in western Himalaya, hosted in the north of the basin. The average severity of wet and dry events is categorized from moderate to severe levels. The average wet severity (AWS) exhibits random spatial patterns, whereas the average dry severity (ADS) is relatively higher in the north of the basin. Observed spatial patterns of maximum wet severity (MWS) and maximum dry severity (MDS) were similar to those of MWD and MDD. The eastern part of the basin experienced wet events of higher severity than the western one, whereas the most severe dry events affected the northwest and central parts of the basin.
Spatial distribution of total wet severity (TWS), total dry severity (TDS), average wet severity (AWS), average dry severity (ADS), maximum wet severity (MWS), and maximum dry severity (MDS) for the period 1981–2014.
Figures 6 illustrates the spatial distribution of intensities of wet–dry
events, calculated as the ratio of severity to duration. The total wet
intensity (TWI) and total dry intensity (TDI) vary from moderate to severe
with a noted range of 1.44 to 1.55 and
Spatial distribution of total wet intensity (TWI), total dry intensity (TDI), average wet intensity (AWI), average dry intensity (ADI), maximum wet intensity (MWI), and maximum dry intensity (MDI) for the period 1981–2014.
The WD ratio features the dominance of wet or dry events for the period of 34 years (1981–2014). The WD ratio for the three severity levels (moderate, severe, and extreme) at pixel basis is presented in Fig. 7. The positive and negative values of WD ratio depict the prevalence of wet and dry events respectively. As the figure shows, higher frequencies of moderate dry events with respect to moderate wet events were found throughout the basin except for a few pixels in the south. By contrast, severe to extreme wet events are more frequent for most parts of the basin. The highest positive values of WD ratio for extreme level of hazard were found in the southwest of the basin, which shows the higher susceptibility of the area towards extreme wet events. Moreover, the analysis of wet–dry event characteristics also revealed the prevalence of wet events with higher duration and severity over monsoon-dominated regions.
Spatial distribution of wet–dry (WD) ratio derived for three
levels of severity (moderate, severe, and extreme) during the period
1981–2014. Blue (
The number of transitions and their average transition time for wet-to-dry and dry-to-wet events for the period 1981–2014 are presented in Figs. 8 and 9. As expected, the number of transitions for wet-to-dry and dry-to-wet events was the highest for the moderate level of events, followed by severe and extreme levels of events. Consequently the average transition time from wet-to-dry and dry-to-wet events was found to be the highest for the extreme level of event followed by severe and moderate levels of events. The number of transitions for moderate, severe, and extreme levels of events varies from 15 to 26, from 6 to 16, and from 1 to 5 respectively. Overall, the number of transitions for dry-to-wet events is larger than the wet-to-dry events for severe and extreme levels of events, whereas the opposite was found for the moderate level of events. The transition time for moderate, severe, and extreme levels of events varies from 1.8 to 6.5, from 1.8 to 16.75, and from 3.5 to 187.0 months respectively. Overall, 53.57 % and 17.86 % of pixels in the UJB showed longer transition time from wet to dry than from dry to wet for moderate and extreme levels, whereas the opposite was seen for severe events.
Number of transitions from wet-to-dry (left) and dry-to-wet (right) events for three levels of severity (moderate, severe, extreme) for the period 1981–2014.
Average transition time (
The wet–dry rapid transition is the consecutive occurrence of wet and dry months of any severity level. The frequency of wet-to-dry (wet month followed by dry month) and dry-to-wet (dry month followed by wet month) rapid-transition events was computed for each grid cell and is shown in Fig. 10. The frequency of wet–dry transition events varies/ranges from 5 to 20 events during the 34 years of the study period. About 50 % of pixels in the UJB encountered a higher number of wet events terminated at dry months. The spatial distribution of frequency of wet–dry rapid-transition events revealed that the wet-to-dry events are less frequent over the westerlies-dominated region of the basin, whereas the southwestern part of the basin was more affected by the abrupt wet-to-dry events. By contrast, abrupt dry-to-wet events are found to be more frequent over pixels surrounding the Himalaya divide line, whereas the remaining part of the basin depicts less incidence of dry-to-wet events.
Frequency of occurrence of abrupt events, wet to dry (left) and dry to wet (right), during the period 1981–2014.
This study attempts to investigate the spatiotemporal variations in wet–dry events collectively, their characteristics (duration, severity, intensity), and transition from wet-to-dry and dry-to-wet events during the period 1981–2014 in the Upper Jhelum Basin (UJB) in South Asia. The SPEI, which incorporates both precipitation and potential evapotranspiration, was used to extract and analyse the wet–dry events. The whole analysis was carried out at the monthly timescale, but the temporal evolution of the basin-averaged index was also simulated at multiple timescales (1, 3, 6, and 12 months). The reason for selecting the monthly timescale for this study is that it is expected to provide the best performance in detecting floods and flash droughts, as longer time steps are more appropriate for long-term droughts only and not for floods.
The results of temporal variations in SPEI showed that the study domain mostly encountered moderate to severe wet–dry events, whereas the extreme wet–dry events rarely occurred during the study period. The results of basin-average SPEI at multiple timescales revealed that the response of SPEI to the deviations in climatic features varies with the accumulation time. Therefore, shorter timescales are more appropriate for detecting frequent seasonal and inter-annual variations, whereas longer timescales provide useful information regarding the signature of the events over the region (Ayugi et al., 2020; Du et al., 2013). Furthermore, the SPEI time-series plots capture the observed extreme floods and drought events that occurred in the basin during the study period well: for instance, the longest drought event occurred from the late 1990s to early 2000s, as evident in Fig. 2 and Table 2. The drought started in 1998 and was considered to be the worst in the history of Pakistan. The drought spell in 2001–2002 resulted in water shortage of up to 51 % of normal supplies (Ahmad et al., 2004). Likewise, the notable flooding events, usually flash floods ranging from moderate to severe, occurred in the years 1988, 1992, 1994, 1997, 2007, and 2014 (Bhat et al., 2019) and were well captured by the SPEI, confirming its valuable contribution to this type of analysis.
An interesting clue to the changing climate is the strong change that occurred in the basin at the end of 1997 (Table 2). Before this change (1981–1997), wet events of different severity levels predominated in the basin, whereas dryer conditions prevailed after 1997. However, it still needs to be investigated whether dryer conditions are expected to continue in the future or whether a large multi-decadal variation is taking place. This strong change in the basin climate coincides with the strongest El Niño–Southern Oscillation (ENSO) event in the winter season of 1997–1998, where the Oceanic Niño Index (ONI) peaked at 2.3 and influenced the climate conditions all over the world (MRCC, 2021). The 1998–2002 drought in southwestern Asia, accompanied by the most severe drought conditions in the last 50 years, was also a result of this strong ENSO event (Ain et al., 2020; Ahmed et al., 2018). ENSO is the primary mode of inter-annual variability, having great influence on global weather and climate via atmospheric circulations (Ullah et al., 2021a). Many researchers reported the close association between variations in atmospheric circulation patterns and climatic variables, extreme weather phenomena like drought and flood (Luca et al., 2020; Omidvar et al., 2016; Sun et al., 2015). Kenyon and Hegerl (2010) examined the response patterns of hydroclimate extremes to ENSO over global land areas and stated a significant decrease in precipitation extremes over Southeast Asia, Indonesia, Australia, and the northernmost region of South America during El Niño phases, whereas in the southern tier of the United States and the region from Argentina to southern Brazil heavy precipitation increased during El Niño phases, and vice versa during La Niña phases. The strength of such connections for Pakistan was also demonstrated in several studies. El Niño suppresses monsoon rainfall activity over Pakistan, while La Nina has a negative impact on winter precipitation over Pakistan (Farooqi et al., 2005; Khan, 2004). Ullah et al. (2021a) found significant impacts of three large-scale climate indices, i.e sea surface temperature (SST) and multivariate El Niño–Southern Oscillation (ENSO4.0), on seasonal droughts across Pakistan.
The results of wet–dry event characteristics (duration, severity, intensity) at pixel basis outline the greater susceptibility of the westerlies-dominated region to dry events with higher duration, severity, and intensity. The dryer conditions in this region could be explained with the increasing rates of global warming over the mountainous region of the basin, also reported by many researchers (Rashid et al., 2020; Shafiq et al., 2020; Zaz et al., 2019). Studies by Negi et al. (2018) and Dimri and Dash (2012) also confirm that most of the western Himalayan region recorded a significant warming trend from 1975 onwards in particular. This is also supported by the tree-ring chronologies of the region, which indicate a rapid growth of the tree rings in recent decades, especially at higher altitudes (Borgaonkar et al., 2009). The impact of global warming on short-term dry events (soil moisture drought) is not straightforward as rising temperature did not necessarily cause increase in actual ET, especially in arid and semiarid regions (Trenberth et al., 2014; Sheffield et al., 2012). In fact the rate and amount of ET results from a complex interaction of temperature, radiation balance, precipitation rates and vegetation physiological control, rather than being exclusively limited by one of these factors. For flash drought, the rapid soil moisture decline should be a result of the intensification of ET driven by higher temperature, which is very common in humid and semi-humid regions, where soil moisture can sustain higher ET amounts up to a few weeks (Yuan et al., 2019). Further decrease in winter and spring precipitation leads to water deficit conditions in this part of the basin. The worst drought event period (2000–2001), partially induced by a stronger ENSO in winter, was also due to the low winter and spring precipitation, as shown in Table 2. During 2000–2001, winter and spring seasons were moderately to severely dry, whereas the monsoon and autumn seasons observed normal months. By contrast, the higher duration and severity of wet events were detected in the monsoon-dominated region, implying that floods mainly occurred during monsoon season with heavy rainfall along with snowmelt. However, the eastern part of the basin was the hotspot of more intense wet events. The above discussion is also supported by the historic database of observed flood events, as most of these events occurred during monsoon season.
The results of the WD ratio showed the prevalence of severe to extreme wet events for most of the basin, while the dry events of moderate severity level were more frequent in the study domain. The southwestern part of the basin, located in the monsoon-dominated region, was found to be the hotspot for the extreme wet events. Moreover, the analysis of wet–dry event characteristics also revealed the prevalence of wet events with higher duration and severity over the same monsoon-dominated region. The spatial patterns of average transition time from one extreme type to the other type was found to be heterogeneous and different for the three severity levels. Overall, a greater number of pixels took a shorter time to switch from dry to wet events than from wet to dry events. Apart from the average transition period, the study domain also experienced rapid transition of wet–dry events. In general, the surrounding region of the Himalaya divide line and the monsoon-dominated part of the basin were found to be the hotspots of rapid wet–dry transition. The rapid wet–dry swings could be explained in the context of global warming. In a warmer climate, increased evapotranspiration rates in response to increased temperature could elevate the drought risk and frequency. At the same time, the prospect of localized heavy precipitation causing floods is expected to increase in response to increased atmospheric moisture content due to increased evapotranspiration rates (He and Sheffield, 2020; Krishnan et al., 2020). Further warming-induced changes in global climate variability, such as El Niño and La Niña, can cause more inter-annual variability or persistence in global weather and climate, significantly affecting regional precipitation and temperature distribution in space and time (Ullah et al., 2021b). Further compelling scientific evidence of human interventions, such as boosted human water intake and land use changes, exacerbates the extreme flood and drought risk hazard.
To conclude, knowledge of wet–dry event characteristics and their rapid transition provides meaningful insight into the geographical hotspots of compound extreme events, which could be of practical value to inform a group of stakeholders (researchers, local authorities, policy makers, relief agencies, non-governmental organizations (NGOs), and (re)insurance companies) on the potential risk. In general, results contribute to hydrological predictability and risk assessment and therefore effectively support disaster preparedness and risk management, ensuring the regional water, food, and socio-economic security and stability against the background of a changing environment. Future work should explore to what extent future wet–dry event frequency will respond to anthropogenic forcing, internal atmospheric processes, and human interventions.
All calculations and plots were produced using ArcMap (version 10.8) and R (version3.3.2) by making use of open-source R packages. The potential evapotranspiration (PET) and standardized precipitation evapotranspiration index (SPEI) were calculated with the R package “SPEI” (version 1.7). The interpolation of PET values was done using the R package “hydroTSM” (Version 0.6-0). Data and code written to generate the figures are available from the corresponding author upon request.
The ERA5 precipitation data can be accessed online (DOI:
The supplement related to this article is available online at:
This paper was conceptualized by RA and GG. RA performed the data analysis and visualization. The original draft was written by RA and revised by GG.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This article is part of the special issue “Recent advances in drought and water scarcity monitoring, modelling, and forecasting (EGU2019, session HS4.1.1/NH1.31)”. It is a result of the European Geosciences Union General Assembly 2019, Vienna, Austria, 7–12 April 2019.
The authors are grateful to the Pakistan Meteorological Department (PMD) and Water and Power Development Authority (WAPDA) for sharing the data.
The authors received funding from the Cooperation Agreement PFK PhD programme 2019–2022 “Partnership for Knowledge-Platform 2: Health and WASH (WAter Sanitation and good Hygiene)” of the AICS-Italian Agency for Development Cooperation to attend higher education programmes in Italy in favour of non-Italian citizens.
This paper was edited by Brunella Bonaccorso and reviewed by two anonymous referees.