Next Article in Journal
Modelling Drought Risk Using Bivariate Spatial Extremes: Application to the Limpopo Lowveld Region of South Africa
Next Article in Special Issue
Analysis of Hydrometeorological Characteristics and Water Demand in Semi-Arid Mediterranean Catchments under Water Deficit Conditions
Previous Article in Journal
A Comparison of Wave Spectra during Pre-Monsoon and Post-Monsoon Tropical Cyclones under an Intense Positive IOD Year 2019
Previous Article in Special Issue
Temporal and Spatial Variability of Dryness Conditions in Kazakhstan during 1979–2021 Based on Reanalysis Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Long-Term Seasonal Drought Trends in the China-Pakistan Economic Corridor

Department of Earth and Atmospheric Sciences, University of Nebraska-Lincoln, Lincoln, NE 68588, USA
*
Author to whom correspondence should be addressed.
Climate 2023, 11(2), 45; https://doi.org/10.3390/cli11020045
Submission received: 9 December 2022 / Revised: 1 February 2023 / Accepted: 2 February 2023 / Published: 12 February 2023
(This article belongs to the Special Issue Coping with Flooding and Drought)

Abstract

:
In recent years, drought events have influenced agriculture, water-dependent industries, and energy supply in many parts of the world. The China–Pakistan Economic Corridor (CPEC) is particularly susceptible to drought events due to large-scale monsoon circulation anomalies. Using the 0.5 × 0.5 resolution rainfall and potential evapotranspiration data set from the Climatic Research Unit (CRU), we assessed the changes in seasonal drought variation and effects of climate variables on drought over the CPEC for the period of 1980 to 2018 using the Standardized Precipitation Evapotranspiration Index (SPEI). Our results show a statistically significant negative trend of SPEI over the hyper-arid region for two monsoons (December–February and June–September) and intra-monsoonal seasons (March–May and October–November), suggesting that the hyper-arid region (southern and southwestern part of the CPEC) is experiencing more frequent drought. A high probability for the occurrence of winter (30–35%) and summer (20–25%) droughts are observed in hyper-arid regions and gradually decreases from south to north of the CPEC. Decreasing seasonal rainfall and increasing potential evapotranspiration with increasing temperature in hyper-arid and arid regions resulted in frequent drought events during the winter monsoon season (from December to February). The findings from this study provide a theoretical basis for the drought management of the CPEC and a framework for understanding changes in drought in this region from climate projections.

1. Introduction

Drought, defined as when water resources are inadequate to meet the demand of people or the environment, is considered one of the most complex natural disasters due to its slow-creeping nature and nonstructural impacts [1]. Model projections of the future climate suggest an increase in drought and drought characteristics in land areas [2], especially over Central America, East Asia, the Amazon, Western Africa, and the Mediterranean [3]. The increasing frequency, severity, and intensity of drought has serious impacts on irrigation, agriculture, and hydropower, hindering the sustainable development of society [4]. For instance, the World Health Organization (WHO) estimated 55 million people who live in drought-vulnerable areas are affected by droughts every year. Climate projections also suggest increases in drought frequency and severity in Asia, Europe, the Americas, and Africa [1,5,6], and the WHO has also highlighted that 700 million people face a high risk of being displaced due to drought by 2030. Therefore, understanding the present and future drought risk in surface and groundwater, irrigation, agriculture, and socio-economic development over different parts of the world is essential for combating food and water insecurity.
South Asia has been among the perennially drought-prone regions of the world, where more than 1.5 billion people have experienced prolonged droughts [7] with high frequency and severity [8]. One of the most noticeable features of the southern Asian climate is typical of monsoon climates. The year-to-year variation of the monsoon rainfall associated with large-scale atmospheric and oceanic variability is a key driver for drought occurrence [9,10]. At the same time, high heat stress and increased air temperature can intensify the drought severity and intensity over the region [11,12] through high evaporation [13,14]. As a result of frequent drought events, drought-induced water scarcity has disrupted agricultural activities in South Asia, resulting in famines in the 19th and 20th centuries [15,16]. In the South Asian region, Pakistan is one of the most vulnerable countries to droughts strongly linked to water scarcity, food security, and the agro-economics of the country [17,18]. The climate risk index (CRI) ranked Pakistan as 5th most-affected country to climate change from 1999 to 2018 [19].
To understand the frequent drought in Pakistan, previous studies investigated the drought variability, trend, and associated mechanisms using different drought indices such as the Standard Precipitation Index (SPI) [17,20], the Reconnaissance Drought Index (RDI) [21], and the Standardized Precipitation Evapotranspiration Index (SPEI) [22], while most studies focus on drought assessment and monitoring in two cropping seasons (Rabi-November to April; Kharif-May to October). For instance, Ahmed et al. [18] found increasing drought intensity over dry and semi-arid regions for cropping seasons. Saadia et al. [17] investigated the precursor conditions that might be engaged for predicting droughts in Pakistan using an SPI. The occurrence of drought in Pakistan is closely related to anomalous total precipitation, near-surface air temperature, large-scale patterns of geopotential height, and prevailing soil moisture conditions [23]. These studies highlight the impact of hydroclimatic disasters (droughts and floods) on socio-economic, agricultural, and environmental conditions in Pakistan. Most of the climate studies in Pakistan focused on variability and trends of rainfall [24,25,26] and temperature [26,27,28].
To improve the social and economic well-being of people in underdeveloped regions of Pakistan, agriculture, health, education, energy, transport, and economic development programs are implemented [29,30] under the China–Pakistan Economic Corridor (CPEC). The CPEC is the pilot project of the Belt and Road Initiative (BRI) initiated in 2015 [31] as a collaboration between the Chinese and Pakistan governments [32]. The CPEC region is sensitive to climate change, desertification, salinization, drought, and human activity [33,34]. For instance, the CPEC is vulnerable to droughts as a result of seasonal changes in rainfall over the region, leading to food security [24]. Du et al. [35] found increasing frequency and strength of extreme precipitation, which caused floods in the CPEC region and projected increasing flood risk in the next 30 years. In addition, the CPEC has experienced dramatic changes in ecological vulnerability as a consequence of climate change and human activity [33,36]. In this context, the project shows some challenges to fully accomplishing the project objectives. Among these constraints, seasonal drought and flood hinder the project significantly.
According to the Ministry of Water Resources in Pakistan (Ministry of Water Resources, 2018), the scientific understanding of retrospective drought assessment, development of drought warning systems, and implementation of drought management plans for the CPEC is very limited [37]. The assessment of the spatio-temporal changes of the drought and flood in the CPEC is essential for adaptation to climate change, building resilient infrastructure, and for the reduction in the risks of disaster in the region. According to the author’s knowledge, few studies focus on the spatial distribution of climatological seasonal drought strongly connected with circulation regimes. Therefore, the present study aims to investigate the seasonal drought trends over the CPEC routes and their relationship with temperature, rainfall, and potential evapotranspiration during 1980–2018. In this study, we used a multi-scale SPEI index that is widely adopted for drought analysis [4,38,39].
The rest of the paper is organized as follows. Section 2 presents the study site, data, and analysis methods used for our study. Section 3 focuses on presenting seasonal drought probability, the spatial trend for drought, rainfall, temperature, and PET at a seasonal time scale. Within Section 4, we compare the results with previous studies and discuss the implications and extensions of our study. Finally, in Section 5, we summarize the study and draw our conclusions on the observed seasonal drought over the CPEC.

2. Materials and Methods

2.1. Study Area

The CPEC region is located in the northwest of South Asia and extends from Kashghar in Xingjiang province (China) to Gawadar port Baluchistan Pakistan [40,41]. Within the CPEC two routes are being constructed with a total length of 4918 km. As shown in Figure 1, the eastern and western routes are located in the East and West of Pakistan, respectively. The eastern route goes across the Punjab and Sindh provinces, while most parts of the Khyber Pakhtunkhwa and Baluchistan Provinces are linked with the western route. The elevation of CPEC is complex and varies from 0 to 8611 m above mean sea level. The mountainous northern part of the CPEC is mostly comprised of glaciers and glacial lakes, while plat terrain is observed in the east part of the CPEC [29,42,43]. The CPEC receives rainfall from the Indian Summer Monsoon (ISM) and westerly disturbances initiated from the Mediterranean Sea. During the summer season, monsoon depressions originating from the Bay of Bengal bring moisture to the east and northeast parts of the CPEC during the months of July to September with high intensity (~55% to 60% of the annual budget). The CPEC shows large spatio-temporal variation in climatic conditions. The coastal climate is observed in the southern part of the CPEC, while the tropical and humid climate is dominant in the central and northern parts of the CPEC [43]. The temperature in the winter is between 4–20 °C in northern Pakistan, while in some years, it is below zero. The average temperature over the southern part of the CPEC is between 20 and 25 °C during winter. The temperatures in the north (<15 °C) and south (35 °C) parts of the CPEC show significant temperature gradients during the summer months.
Based on the rainfall and temperature distribution, Rahman et al. [44] and Rahman et al. [45] have demarcated the four climate groups; hyper-arid (R-I), arid (R-II), humid (R-III), and glacial (R-IV) regions as shown in Figure 1. More than 60% of Pakistan’s population resides in rural areas [37], and most of the population is concentrated in the semi-arid to hyper-arid plains of the lower Indus basin, that experience frequent drought events. In this study, we also used previously defined four climate regions for further analysis. In irrigated areas of the CPEC, wheat, maize, and rice are the main food crops, while cotton and sugarcane are considered economic crops in the region.

2.2. Data and Methodology

In this study, we use monthly rainfall, temperature, and potential evapotranspiration data from the Climatic Research Unit gridded Time Series (CRU TS) available at 0.5 × 0.5 degrees. This data set is developed by the National Centre for Atmospheric Science (NCAS) at the University of East Anglia (https://crudata.uea.ac.uk/cru/data/hrg/ accessed on 12 August 2022) and is widely used for drought analysis [46,47,48]. The spatial distribution of soil moisture in different seasons is plotted using CPC Soil Moisture V2 (https://psl.noaa.gov/data/gridded/data.cpcsoil.html accessed on 16 January 2023).

Calculation of Standardized Precipitation Evapotranspiration Index (SPEI)

Standardized Precipitation Index (SPI) and Standardized Precipitation Evapotranspiration Index (SPEI) [49,50,51,52,53,54] are most widely adopted for drought studies over the world because these multi-scaler indexes can be used to monitor and assess drought characteristics. Compared to the SPI, the influence of PET on drought occurrence frequency, intensity, and severity [55,56,57] is represented by SPEI [58]. In this study, we calculate the SPEI index for the investigation of seasonal drought variations in the CPEC. In addition, the SPI index is also used to support the validity of the SPEI results.
According to Vicente-Serrano et al. [58], the balance between precipitation (P) and potential evapotranspiration (PET) is the basic form of SPEI (Equation (1)) which can be calculated at different timescales using Equation (2) or Equation (3)
D = P P E T
D i , j   k = l = j k + 1 j ( P i , l P E T i , l )   if   j k
D i , j   k = l = 13 k + j 12 ( P i 1 , l P E T i 1 , l ) + l = 1 j P i , l P E T i , l   if   j < k
where monthly potential evapotranspiration and precipitation are denoted by PET and P, respectively, while the aggregated value of monthly precipitation minus potential evapotranspiration for k timescale (months) is represented by D. The indexes i (year) and j (month) depend on the selected timescale (k). The log-logistic probability distribution is used to transform D i , j k , into standardized units to get multi-scalar SPEI [58].
The gridded rainfall and PET are then used to calculate the SPEI at 1-month (SPEI-1), to 12-month (SPEI-12) timescales. Previous studies have set “−1” SPEI as the threshold value to identify drought [59,60,61], and the negative SPEI values is used to determine moderate, severe, and extreme drought. In this study, SPEI-3 at February and SPEI-3 at May are selected to represent droughts in DJF and MAM seasons, respectively. For instance, as the DJF season spans from December to February, in order to consider the drought over this season, SPEI-3 February is calculated using rainfall and evapotranspiration from December to February. For the JJAS season, SPEI-4 at September is extracted from SPEI Index accumulated over four months (SPEI-4).
The nonparametric Mann-Kendall test (MK test) and Sen’s slope estimator based on Kendall’s tau (τ) have been frequently used to analyze trends in drought [62,63,64]. We also used MK test and Sen’s slope estimator to identify the spatial drought trend for the different seasons. We looked for a significant increase or decrease of seasonal rainfall, temperature, and PET at the 90% significance level.

3. Results

3.1. Rainfall Distribution of the CPEC

Figure 2 shows the annual cycle of rainfall distribution of different sub-regions in Pakistan. Across the entire region of Pakistan (Figure 2a) and three different climate regions (Figure 2c–e) produce bimodal rainfall distributions where the dominant peak occurs in JJA, and another peak is observed in March. The long-term average monthly rainfall in R-I is less than 20 mm for all the months except July and August. In R-II, the average monthly rainfall is less than 35 mm, while July and August receive 50–55 mm rainfall (Figure 2b,c). R-III gets more rainfall in all the months compared to the other regions, where the maximum rainfall (140 mm) was recorded from July to August (Figure 2d). R-IV produces a different distribution than the other regions (Figure 2f), with monthly rainfall from June to November of less than 30 mm and maximum rain (50 mm) observed in March and April. The lowest rainfall for all the regions is recorded in October and November, where R-I and R-II get less than 5 mm per month while R-III gets 20 mm rainfall for this period.
The findings highlight that rainfall distribution over the CPEC is closely associated with the monsoon circulation over the South Asian region. Ali et al. [25] revealed that the ISM and winter monsoon bring a significant amount of rainfall to Pakistan. July and August are the active periods of the ISM for Pakistan compared to the four months of the summer monsoon rainfall season. Wang et al. [65] found that heavy rainfall occurs in July and August in this region, and Bhatti et al. [66] observed many extreme rainfall events in August. During the December to February period, extratropical storms initiated in the Mediterranean are responsible for the rain in the northwestern parts of Pakistan [67].
Figure 3 shows the spatial distribution of season rainfall, maximum temperature, potential evapotranspiration, and soil moisture over Pakistan for the 1980–2018 periods. During the DJF and MAM seasons, R-I and R-II received less than 50 mm per month of rainfall, while R-III gets relatively more rainfall than other regions (Figure 3a,b). The JJAS rainfall is more concentrated in the R-III, while other regions get less rainfall, especially the western part of R-I (Figure 3c). Compared to the other season, JJAS rainfall contributes to annual total rainfall in R-I, R-II, and R-III at 57.0%, 53.6%, and 52.1%, respectively; however, for R-IV, it only contributes 24% to the annual total rainfall. During the ON season, Pakistan experiences dry conditions with less rainfall (Figure 3d), contributing less than 5% in R-I, R-II, and R-III and 9.3% in R-IV. Rain from the ISM and winter monsoon seasons contributes more than 70% of the annual rainfall in all regions except R-IV.
As shown in Figure 3i–l, the lowest Tmax is observed in R-IV during the DJF season compared to other regions. R-I recorded the highest Tmax (>40 °C) in the JJAS season, followed by MAM and ON seasons. In general, the temperature in the eastern part of R-II and R-III is higher than in the western part of these two regions (Punjab and Sindh provinces). Similar to our findings, station-based observations also found that the eastern part of the country [29] is warmer (26 to 30 °C) than the western and southwestern hilly areas (21 to 25 °C).
To determine the regional dryness and wetness and variations in meteorological conditions, PET can be used as an indicator [68,69]. Therefore, the spatial distribution of long-term PET climatology for Pakistan in different seasons is investigated (Figure 3). The lowest and highest PET is recorded in DJF and JJAS seasons, respectively. Among the regions, the lowest PET is recorded in R-IV for all seasons. Notably, the high PET is observed in R-I, and R-II occurs when Tmax is largest, suggesting that PET has a strong correlation with temperature over the region. Adnan et al. [69] used station-based PET observations and revealed that evapotranspiration in Pakistan was positively correlated with temperature, solar radiation, and wind speed. Importantly, Smakhtin and Schipper [70] found an increase in evapotranspiration under high temperatures and low humidity, which worsens the drought condition in both arid and humid regions. Based on these findings, the high temperature and low humidity may exacerbate the drought intensity and severity over R-I and R-II. Soil moisture plays a critical role in agriculture and can be used to the monitoring of the intensity of droughts [71]. As shown in Figure 3m–p, the lowest soil moisture is observed in R-I and R-II, while R-III shows the highest soil moisture irrespective of the seasons. Compared to other panels, soil moisture in the four regions is primarily associated with seasonal precipitation, while excesses in evapotranspiration (ET) due to high temperature has an additional impact, especially over the hyper-arid and arid regions.

3.2. Drought Variation in the CPEC

To identify the temporal variation in drought periods, a diagram showing the SPEI at different timescales (1 to 12 months) between the period 1980 and 2018 was generated for different climate regions in Pakistan (Figure 4). The cumulative effects of the temporal changes in wet and dry periods become more apparent with longer periods as the time scale of SPEI increases. As shown in Figure 4a,b, strong wet events are observed in R-I and R-II during the 1980–1998 period, which is reversed after 1999. For instance, the 1999–2003 and 2014–2018 periods exhibit continuous droughts over R-I (Figure 4a). Interestingly, we found that R-II and R-III also experienced prolonged drought events from 1999 to 2003, while the drought magnitude (−2.5 to −1.5) is much higher than R-I (−1.5 to −1.0). Similar to our findings, Ahmad et al. [72] also identified this period as the most severe drought condition in the last 50 years in Pakistan. Ahmad et al. [73] revealed that the prolonged drought from 1999 to 2003 reduced the irrigation water supply, cropped area, and crop production of Sindhi (−4% to −40%) and Baluchistan (−8% to −20%) provinces in Pakistan.
Compared to the other regions, R-IV experiences fewer drought events (Figure 3d). The findings highlight that R-I and R-II have been more susceptible to droughts in recent decades. Similar to our drought evaluation in R-I and R-II, the Indo-Gangetic plains also experienced many drought events during 2002–2004, 2009–2014, and 2015–2017 [74]. We also found that extreme wet events are dominated in the active monsoon region in Pakistan (R-III) might influence frequently occurring flood events reported over the region. Meanwhile, we were able to identify seasons with strong precipitation using these diagrams. For example, in 2010 and 2015, record-breaking monsoon rainfall was well reproduced in R-III and R-IV (Figure 4c,d), which was associated with massive flooding events (citation for flooding events). Shah et al. [75] reported severe floods in Punjab province (R-III) in 2010 with extreme rainfall in the northern parts of Pakistan. We also compared the SPEI and SPI time series in each climate region to validate the SPEI (Figure S1). The correlation between SPEI-3 and SPI-3 in all the regions is statistically significant at a 95% confidence level suggesting that the SPEI well captures the drought events in different climate regions.

3.3. Seasonal Drought Frequency in the CPEC

In addition to seeing how the indexes for SPEI changed for the regions, we also calculated the spatial distribution of seasonal drought frequency in Pakistan for the 1980–2018 period (Figure 5). During the DJF season, almost all of the country experiences droughts more than 15% of the time, with the strongest drought occurrence found in the southern part of the domain, especially in R-I (Figure 5a). R-I and the south part of the R-II region experience droughts more than 20% of the time during the MAM period as compared to the rest of the area (Figure 5b). During the JJAS season, drought occurrence frequency in the southwest part of the R-I region is much higher than in other regions (Figure 5c). The southern and southeast parts of R-I and R-II also experienced drought events more than 20% of the time during the ON season for the 1980–2018 period (Figure 5d). R-I experiences frequent drought events (more than 25% of the time), while drought is less frequently observed in humid and glacial regions. Similar to our findings, Jamro et al. [76] also found frequent seasonal droughts in the arid region of Balochistan province, Pakistan. Naz et al. [77] used station data sets to calculate the SPI index and observed that the western/southwestern part of R-I and R-II has frequent droughts, which is consistent with our findings.

3.4. Spatial Trend of Standardized Precipitation Evapotranspiration Index (SPEI)

While understanding the temporal variability in each region is illuminating, it is also important to consider how SPEI changes with time across the domain without being constrained by our four regions. Figure 6 shows the spatial distribution of the trend (per decade) of SPEI for the seasonal drought in Pakistan. The negative SPEI trend indicates increasing drought conditions, while the positive SPEI represents increased wetness. During the DJF and MAM seasons, a statistically significant decreasing SPEI trend is observed over R-I, suggesting that droughts are increasing from 1998 to 2018 (Figure 6a,b). The south and southwest parts of R-I, R-II, and northern parts of R-IV recorded statistically significant decreasing SPEI (Figure 6c). As shown in Figure 6d, most parts of R-I and R-II exhibit statistically significant negative SPEI, while the rest of the country has a nonsignificant negative trend for the 1980–2018 period (Figure 6d). Considering the spatial distribution of drought trends, this figure demonstrates that, for the most part, Pakistan is more vulnerable to seasonal droughts under a warming climate. This may ultimately affect crop production since Ashraf et al. [78] found that the agricultural vulnerability to drought risks increased significantly in Baluchistan (the southwest part of R-I and R-II). The calculated SPI based on station data also shows that seasonal drought is increasing in Baluchistan province [77].

3.5. Spatial Trend of Rainfall, Evapotranspiration, and Temperature

Drought occurs when there is below-average rainfall over for a protracted period, such as a few months or a season resulting in water scarcity over the region. Recent studies have found that sparse rainfall is a major reason for drought in India [74], and less rainfall coupled with increasing temperature is likely to cause seasonal drought in Balochistan, Pakistan [78]. Therefore, we investigate the seasonal rainfall trend for 1980–2018 to explain the observed drought trend over Pakistan (Figure 7). The rainfall in the DJF season shows a significant decrease over the western and southwestern parts (>3 mm per decade) of R-I and R-II, while the rest of the region (including part of R-III) depicts a statistically nonsignificant decreasing trend (Figure 7a). Ullah et al. [24] used station-based rainfall data and also found a decreasing trend of DJF rainfall in these regions.
During the MAM season, the spatial pattern of the rainfall trend is quite similar to the DJF season, but the trend magnitude is comparatively less and a small part of R-I shows a statistically significant negative rainfall trend (Figure 7b). We also identify increasing rainfall in R-III and decreasing rainfall in R-I and part of R-II during the JJAS season (Figure 7c), with notable decreases over the southwest/west part of the regions. While these are prominent trends, they are not statistically significant due to the large amount of variability from year to year in these regions (Figure S2).
Rainfall in the ON season (Figure 7d) declined by around 1 mm per decade from 1980 to 2018. We found that very little rainfall receives during the ON season, thus a declining trend in the rainfall. Hanif et al. [79] also observed a decreasing precipitation trend in the southern part of Pakistan as well as along the coastal belt and parts of Sindh province.
The role of PET in drought is more important under global warming and changing atmospheric water demands [80,81]. Therefore, the spatial distribution of the long-term trend of PET over Pakistan is analyzed for the 1980–2018 period (Figure 8). During the two monsoon seasons (DJF and JJAS), a statistically significant increasing trend of PET is observed over R-I, while the trend magnitude in the JJAS season is higher than in the DJF season (Figure 8a,c). The increasing trend magnitude of PET for the MAM season in R-I and R-II is larger than in other seasons (Figure 8b). In the humid region, we observed statistically significant decreasing PET for all the seasons, while the largest decreasing trend is observed for the JJAS season (Figure 8c). The PET over the south and southeast part of the R-I is increasing during the ON season (Figure 8d).
In contrast, the increasing trend of PET over R-III and the upper part of the R-II are observed during the JJAS season. Ahmed et al. [82] also observed statistically significant PET trends in annual and two-crop growing seasons in recent years in Pakistan. Similar to our findings Goroshi et al. [83] found decreasing PET over monsoon regions in India.
Increasing temperature generally leads to increased evaporative demand [84,85]. Sharma and Mujumdar [86] revealed that the combined effects of decreasing rainfall and abnormally high temperatures (heatwaves) encourage drought conditions to build and intensify than their occurrence. This study also analyzed the long-term trend of maximum and minimum temperatures in different seasons (Figure 9). The maximum and minimum temperature over Pakistan in MAM and ON seasons showed a statistically significant increasing trend where the magnitude of the trend for both Tmax and Tmin (>0.4 °C per decade) in the MAM season is higher than in other seasons (Figure 9b,f). Iqbal et al. [87] also observed the sharpest increases in Tmax and Tmin in the MAM season based on station-based observation.
During the JJAS season, Tmax and Tmin in R-I and south, southwest part of R-II showed a statistically significant increasing trend, while Tmax in R-IV also depicts a statistically significant increasing Tmax (Figure 9c). Previous studies found that Tmax and Tmin in temperature in Punjab, Pakistan increased significantly at the annual time scale [88]. Station-based observation also indicates warming in most temperature extreme indices over Pakistan [29]. In R-I and R-IV, we observed that Tmin increases faster than Tmax during all the seasons except the DJF season in R-I (Figure 9a,e). Similarly, other countries in the South Asian region have depicted an increasing trend of Tmax and Tmin where the magnitude of increasing Tmin is larger than Tmax [89,90].

3.6. Interactions between Precipitation, Evapotranspiration, and SPEI

Precipitation and potential evapotranspiration are considered to calculate the SPEI in determining drought [91]. Therefore, we investigate the spatial correlation between precipitation, evapotranspiration, and SPEI, as shown in Figure 10. During the MAM, JJAS, and ON seasons, the statistically significant negative correlation between SPEI and rainfall was observed (Figure 10a–c) with a 0.9 correlation coefficient; however, the relationship between SPEI and rainfall in the winter season is comparatively less than in other seasons (Figure 10d). Yao et al. [51] revealed that rainfall was the main factor causing droughts in Xinjiang, China. The correlation between SPEI and PET over all regions is positive and statistically significant for MAM, JJAS, and ON seasons (Figure 10 a–c). The correlation between SPEI and PET is weaker (<0.4) during the DJF season (Figure 10h) and not significant for most of the domains. The findings of this study are consistent with the results of other sensitivity analyses of rainfall PET and SPEI [51,92]. However, this study shows that the relationship between SPEI, Rainfall, and PET depends on the regions and seasons.

4. Discussion

Previous studies have suggested that the CPEC is highly sensitive to climate change [24,41] since most parts of the CPEC are in arid and semi-arid areas, where water resources are already scarce [41]. Frequent drought events, increasing agricultural production, and increasing population have further worsened the demand for water resources in this region [93]. To better understand the historical change in a drought, we used a CRU gridded data set and multi-scalar drought indices (SPEI) to investigate seasonal droughts in both time and space over the CPEC routes during 1980–2018. We found that across most of the southern half of Pakistan (R-I and R-II), a statistically significant decrease in SPEI was observed during the period.
When considering the water balance equation, rainfall is the supply side of the water balance that represents one aspect of drought. Actual evapotranspiration (ET) is the demand side of the water balance which depends on the moisture availability and evaporative demand (PET). Evapotranspiration is a function of temperature, humidity, solar radiation, and wind speed. With increasing temperatures, PET will increase the demand side of the water balance resulting in more severe droughts. Hyper-arid and arid regions receive less rainfall compared to other regions, and high temperatures over the region are associated with more evapotranspiration, suggesting that R-I and R-II are more susceptible to short-term and long-term droughts. Mishra and Liu [94] revealed that the risk of drought increases with fewer light precipitation days, and prolonged dry spells increase with increasing temperature. In addition, increasing the temperature intensifies the drought conditions and supports drought propagation through land-atmospheric interaction [95,96], suggesting that increasing temperature contributes to intensified drought conditions in Pakistan through an increase in evaporative demand. According to Konapala et al. [97], variations of seasonal and annual mean precipitation and evaporation influence patterns of water availability [98], crop yield [99], and ecology [100], ultimately impacting society and ecosystems. Increased PET over the R-I and R-II, in turn, enhances evapotranspiration on land and speeds up the drying of soils. Consequently, drought conditions are exacerbated and intensify droughts over the region.
In Pakistan, two cropping seasons are based on the monsoon rainfall: Rabi crop depends on the winter monsoon (DJF) rainfall, while Kharif crops are cultivated during the summer monsoon (JJAS) [101,102]. Under the warming climate, the decreasing rainfall and increasing evapotranspiration intensify plant drought stress and reduce/degrade both Rabi and Kharif vegetation in water-limited hyper-arid and arid regions. In parallel with our findings, Kharif and Rabi aridity over the southern part of the country is increased [101]. It was just anticipation based on the assumption that increased drought due to decreasing rainfall and rising temperature intensifies potential evapotranspiration (PET) over the region. As a result of water scarcity induced by frequent drought events, agricultural production will decrease in arid and semi-arid regions and end with food insecurity for 70–80% of the population in Pakistan who are engaged with the environment and weather-vulnerable farming [78].
As a response to increasing temperature, glacier retreat ultimately leads to changes in water discharge with gradually decreasing available ice reserves [103]. Kääb et al. [104] found that glaciers in the Jammu–Kashmir region (located in R-IV) are melting at the rate of 0.66 m/year under global warming. With this perspective, increasing Tmax and Tmin in the glacial region in Pakistan results in glaciers melting faster than normal, which can trigger flash flooding in the northern part of Pakistan. For instance, a glacial lake outburst flood in 2022 destroyed Hassanabad Bridge in Pakistan’s northern Gilgit-Baltistan region. We found that over the humid region (R-III) an increase in summer rainfall was associated with decreased drought. Safdar et al. [105] also found an increase in monsoon-hit areas located in the northern part of Pakistan with rainfall ≥2.5 mm/day. According to the climate projection studies using CMIP6 models [106], the ISM is increasing by 0.33 mm/day and 5.3% per Kelvin of global warming. In Pakistan, monsoon rainfall is one of the major reasons for floods, while the snowmelt intensifies flooding in rivers during the monsoon period [75].
Monsoonal variability and inconsistency can be partially attributed to large-scale circulations such as IOD and ENSO [105]. Previous studies in Pakistan revealed that rainfall variability is strongly associated with large-scale teleconnections [17,107]. Therefore, understanding the influence of the tropical sea surface temperature (SST) anomalies on observed seasonal drought variability over Pakistan may help explain the increase in intensity and frequency of droughts in recent decades. To identify the dominant precursory oceanic signals on seasonal drought characteristics, numerical simulations with a climate system model and observation-based statistical analysis will be needed. For instance, statistical analysis on reanalysis data found that circulations associated with Pacific and Indian Ocean SSTs are highly associated with drought episodes across the Pakistan region [17].
Except for meteorological forcings and their variability, droughts are initiated and maintained by several factors, such as soil and vegetation feedback, agricultural practices, and management choices, including irrigation and grazing density. Among these factors, soil moisture variability and trends are crucial to quantify the intensity of droughts [71]. In this study, drought analysis is based only on temperature and precipitation and does not make use of soil moisture content and moisture deficits in Pakistan. Therefore, future studies may wish to investigate the influence of drought on soil moisture in both the surface and root zones, as well as identify how changes in soil moisture impact the changing drought using in situ observations, integrating satellite soil moisture data or soil moisture data from hydrological or land surface models.
Though it is challenging to set up drought policies and preparedness plans, the outcomes of the process can significantly increase societal resilience to drought and drought-related disasters. Based on the findings, drought in the different climate regions depends on the climatological seasons. In addition, we identified hyper-arid and arid regions as the most drought-prone regions in the CPEC. These facts encourage the improvement and application of seasonal and shorter-term drought forecasts. Furthermore, we highlighted the priority of developing integrated monitoring and drought early warning systems and associated information delivery at regional scales, especially over the hyper-arid and arid regions. The findings of this study are also helpful in raising awareness at the regional level and developing a host of preventive and mitigation measures against droughts. In the drought policy and preparedness process, inventory data and financial resources available and identifying groups at risk are crucial. We observe that drought in the hyper-arid and arid regions is intensifying. Therefore, government and non-government organizations can prioritize identifying the people at high risks, such as agricultural and livestock communities. Establishing a drought task force is recommended to understand drought mitigation techniques, risk analyses (economic, environmental, and social aspects), and drought-related decision-making processes.

5. Conclusions

The CPEC is very sensitive to natural disasters such as floods and droughts under global warming. Therefore, monitoring drought and floods at seasonal timescales is a vital and timely task to mitigate the impacts of recurrent and unexpected drought and flood events. For this reason, this study investigates the seasonal drought variation and its associated trends from 1980 to 2018 over the region using a multi-scalar drought index (SPEI) which represents the effect of precipitation and PET on drought characteristics. For calculating SPEI at a seasonal time scale, long-term (1980–2018) 0.5° × 0.5° resolution gridded data of monthly PET and rainfall was used. The results suggest that hyper-arid and western and southwestern parts of arid regions experience frequent seasonal droughts compared to the rest of the region in Pakistan. We also found that an increase in drought occurrence in hyper-arid and southern and southwestern parts of the arid regions was attributed to decreasing seasonal rainfall and increasing potential evapotranspiration with increasing temperature.
This study provides a scientific understanding of seasonal droughts for drought monitoring and implementing drought early warning systems at local and national levels. However, effective communication of the results at local, state, and national levels, with environmental, social, and economic perspectives, is crucial. In addition, predicting the probability of occurrence of seasonal droughts and establishing a skillful drought prediction system is essential and useful to society concerning water resource management and agricultural planning. Niranjan Kumar et al. [108] found that predicting seasonal droughts, especially during the monsoon season is challenging. Therefore, we have to improve our regional and global models’ ability to simulate the drivers of hydrological change across this region and increase model skill in seasonal to decadal prediction of precipitation change. The findings of our study can be used for implementing drought adaptation strategies for future droughts expected in the coming decades. Therefore, identifying regional drought risk and improving the prediction of regional drought risk under future scenarios is crucial to provide stakeholders with useful information in order to prepare for changes in extreme regional events.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cli11020045/s1, Figure S1 Time series of SPEI-3 and SPI-3 for (a) hyper-arid (b) arid (c) humid and (d) glacial regions in CPEC from 1980 to 2018. The positive and negative SPEI-3 values are shown by blue and pink colored shadings, respectively. The black line indicates an SPI-3 time series. Figure S2: Standard deviation of seasonal rainfall during (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) season over the China Pakistan Economic Corridor for 1980–2018; Figure S3: Spatial distribution of self-calibrating Palmer Drought Severity Index using Pen-man-Monteith (sc_PDSI_pm) during the (a) DJF (December to February; SPEI-3 at February), (b) MAM (March to May; SPEI-3 at May), (c) JJAS (June to September SPEI-4 at September), and (d) ON (October to November; SPEI-2 at November) Season over Pakistan for 1980–2018. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.

Author Contributions

Conceptualization, S.S. and R.D.D.; methodology, S.S.; formal analysis, S.S.; writing—original draft preparation, S.S.; writing—review and editing, R.D.D.; visualization, S.S.; supervision, R.D.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The gridded rainfall, temperature, and potential evapotranspiration data are freely available from the National Centre for Atmospheric Science (NCAS) at the University of East Anglia (https://crudata.uea.ac.uk/cru/data/hrg/ accessed on 12 August 2022).

Acknowledgments

We thank the National Centre for Atmospheric Science (NCAS) and the Climatic Research Unit (CRU) at the University of East Anglia for providing the CRU TS dataset. We thank the Editor and three anonymous reviewers for their comments which helped improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vicente-Serrano, S.M.; Peña-Angulo, D.; Beguería, S.; Domínguez-Castro, F.; Tomás-Burguera, M.; Noguera, I.; Gimeno-Sotelo, L.; El Kenawy, A. Global drought trends and future projections. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2022, 380, 20210285. [Google Scholar] [CrossRef] [PubMed]
  2. Nasrollahi, N.; AghaKouchak, A.; Cheng, L.; Damberg, L.; Phillips, T.J.; Miao, C.; Hsu, K.; Sorooshian, S. How well do CMIP5 climate simulations replicate historical trends and patterns of meteorological droughts? Water Resour. Res. 2015, 51, 2847–2864. [Google Scholar] [CrossRef]
  3. Spinoni, J.; Naumann, G.; Carrao, H.; Barbosa, P.; Vogt, J. World drought frequency, duration, and severity for 1951–2010. Int. J. Climatol. 2014, 34, 2792–2804. [Google Scholar] [CrossRef]
  4. Liu, C.; Yang, C.; Yang, Q.; Wang, J. Spatiotemporal drought analysis by the standardized precipitation index (SPI) and standardized precipitation evapotranspiration index (SPEI) in Sichuan Province, China. Sci. Rep. 2021, 11, 1280. [Google Scholar] [CrossRef] [PubMed]
  5. Vicente-Serrano, S.M.; Domínguez-Castro, F.; McVicar, T.R.; Tomas-Burguera, M.; Peña-Gallardo, M.; Noguera, I.; López-Moreno, J.I.; Peña, D.; El Kenawy, A. Global characterization of hydrological and meteorological droughts under future climate change: The importance of timescales, vegetation-CO2 feedbacks and changes to distribution functions. Int. J. Climatol. 2020, 40, 2557–2567. [Google Scholar] [CrossRef]
  6. Touma, D.; Ashfaq, M.; Nayak, M.A.; Kao, S.-C.; Diffenbaugh, N.S. A multi-model and multi-index evaluation of drought characteristics in the 21st century. J. Hydrol. 2015, 526, 196–207. [Google Scholar] [CrossRef]
  7. Aadhar, S.; Mishra, V. Increased Drought Risk in South Asia under Warming Climate: Implications of Uncertainty in Potential Evapotranspiration Estimates. J. Hydrometeorol. 2020, 21, 2979–2996. [Google Scholar] [CrossRef]
  8. Mishra, V.; Tiwari, A.D.; Aadhar, S.; Shah, R.; Xiao, M.; Pai, D.S.; Lettenmaier, D. Drought and Famine in India, 1870–2016. Geophys. Res. Lett. 2019, 46, 2075–2083. [Google Scholar] [CrossRef]
  9. Shelton, S.; Pushpawela, B. Observed southwest monsoon rainfall changes in Sri Lanka and possible mechanisms. Model. Earth Syst. Environ. 2022, 8, 4165–4175. [Google Scholar] [CrossRef]
  10. Pathak, A.; Ghosh, S.; Martinez, J.A.; Dominguez, F.; Kumar, P. Role of Oceanic and Land Moisture Sources and Transport in the Seasonal and Interannual Variability of Summer Monsoon in India. J. Clim. 2016, 30, 1839–1859. [Google Scholar] [CrossRef]
  11. Davis, K.F.; Chhatre, A.; Rao, N.D.; Singh, D.; DeFries, R. Sensitivity of grain yields to historical climate variability in India. Environ. Res. Lett. 2019, 14, 064013. [Google Scholar] [CrossRef]
  12. Mishra, V.; Thirumalai, K.; Singh, D.; Aadhar, S. Future exacerbation of hot and dry summer monsoon extremes in India. NPJ Clim. Atmos. Sci. 2020, 3, 10. [Google Scholar] [CrossRef]
  13. Iyakaremye, V.; Zeng, G.; Ullah, I.; Gahigi, A.; Mumo, R.; Ayugi, B. Recent Observed Changes in Extreme High-Temperature Events and Associated Meteorological Conditions over Africa. Int. J. Climatol. 2022, 42, 4522–4537. [Google Scholar] [CrossRef]
  14. Sun, C.; Zhu, L.; Liu, Y.; Hao, Z.; Zhang, J. Changes in the drought condition over northern East Asia and the connections with extreme temperature and precipitation indices. Glob. Planet. Chang. 2021, 207, 103645. [Google Scholar] [CrossRef]
  15. Aadhar, S.; Mishra, V. A substantial rise in the area and population affected by dryness in South Asia under 1.5 °C, 2.0 °C and 2.5 °C warmer worlds. Environ. Res. Lett. 2019, 14, 114021. [Google Scholar] [CrossRef]
  16. Chandrasekara, S.S.K.; Kwon, H.-H.; Vithanage, M.; Obeysekera, J.; Kim, T.-W. Drought in South Asia: A Review of Drought Assessment and Prediction in South Asian Countries. Atmosphere 2021, 12, 369. [Google Scholar] [CrossRef]
  17. Saadia, H.; Farhan, S.; Arfan, A.; Alina, H.; Irfan, U. Droughts over Pakistan: Possible cycles, precursors and associated mechanisms. Geomat. Nat. Hazards Risk 2021, 12, 1638–1668. [Google Scholar] [CrossRef]
  18. Ahmed, K.; Shahid, S.; Nawaz, N. Impacts of climate variability and change on seasonal drought characteristics of Pakistan. Atmos. Res. 2018, 214, 364–374. [Google Scholar] [CrossRef]
  19. Eckstein, D.; Künzel, V.; Schäfer, L.; Winges, M. Global Climate Risk Index 2020; Germanwatch: Bonn, Germany, 2019. [Google Scholar]
  20. Qaisrani, Z.N.; Nuthammachot, N.; Techato, K.; Asadullah; Jatoi, G.H.; Mahmood, B.; Ahmed, R. Drought variability assessment using standardized precipitation index, reconnaissance drought index and precipitation deciles across Balochistan, Pakistan. Braz. J. Biol. Rev. Brasleira Biol. 2022, 84, e261001. [Google Scholar] [CrossRef]
  21. Ashraf, M.; Ullah, K.; Adnan, S. Satellite based impact assessment of temperature and rainfall variability on drought indices in Southern Pakistan. Int. J. Appl. Earth Obs. Geoinf. 2022, 108, 102726. [Google Scholar] [CrossRef]
  22. Akhtar, T.; Mushtaq, H.; Hashmi, M.Z.R. Drought monitoring and prediction in climate vulnerable Pakistan: Integrating hydrologic and meteorologic perspectives. Hydrol. Earth Syst. Sci. Discuss. 2020, 2020, 1–29. [Google Scholar] [CrossRef]
  23. Saleem, F.; Arshad, A.; Mirchi, A.; Khaliq, T.; Zeng, X.; Rahman, M.M.; Dilawar, A.; Pham, Q.B.; Mahmood, K. Observed Changes in Crop Yield Associated with Droughts Propagation via Natural and Human-Disturbed Agro-Ecological Zones of Pakistan. Remote Sens. 2022, 14, 2152. [Google Scholar] [CrossRef]
  24. Ullah, S.; You, Q.; Ullah, W.; Ali, A. Observed changes in precipitation in China-Pakistan economic corridor during 1980–2016. Atmos. Res. 2018, 210, 1–14. [Google Scholar] [CrossRef]
  25. Ali, G.; Sajjad, M.; Kanwal, S.; Xiao, T.; Khalid, S.; Shoaib, F.; Gul, H.N. Spatial–temporal characterization of rainfall in Pakistan during the past half-century (1961–2020). Sci. Rep. 2021, 11, 6935. [Google Scholar] [CrossRef] [PubMed]
  26. Khan, F.; Ali, S.; Mayer, C.; Ullah, H.; Muhammad, S. Climate change and spatio-temporal trend analysis of climate extremes in the homogeneous climatic zones of Pakistan during 1962–2019. PLoS ONE 2022, 17, e0271626. [Google Scholar] [CrossRef] [PubMed]
  27. Hussain, A.; Cao, J.; Hussain, I.; Begum, S.; Akhtar, M.; Wu, X.; Guan, Y.; Zhou, J. Observed Trends and Variability of Temperature and Precipitation and Their Global Teleconnections in the Upper Indus Basin, Hindukush-Karakoram-Himalaya. Atmosphere 2021, 12, 973. [Google Scholar] [CrossRef]
  28. Latif, Y.; Yaoming, M.; Yaseen, M.; Muhammad, S.; Wazir, M.A. Spatial analysis of temperature time series over the Upper Indus Basin (UIB) Pakistan. Theor. Appl. Climatol. 2020, 139, 741–758. [Google Scholar] [CrossRef]
  29. Ullah, S.; You, Q.; Ullah, W.; Ali, A.; Xie, W.; Xie, X. Observed changes in temperature extremes over China–Pakistan Economic Corridor during 1980–2016. Int. J. Climatol. 2019, 39, 1457–1475. [Google Scholar] [CrossRef]
  30. Kouser, S.; Subhan, A.; Abedullah. Uncovering Pakistan’s Environmental Risks and Remedies under the China-Pakistan Economic Corridor. Environ. Sci. Pollut. Res. 2020, 27, 4661–4663. [Google Scholar] [CrossRef]
  31. Ahmed, Z.S. Impact of the China–Pakistan Economic Corridor on Nation-Building in Pakistan. J. Contemp. China 2019, 28, 400–414. [Google Scholar] [CrossRef]
  32. Saad, A.; Ijaz, M.; Asghar, M.U.; Yamin, L. China-Pakistan economic corridor and its impact on rural development and human life sustainability. Observations from rural women. PLoS ONE 2020, 15, e0239546. [Google Scholar] [CrossRef]
  33. Wu, H.; Guo, B.; Fan, J.; Yang, F.; Han, B.; Wei, C.; Lu, Y.; Zang, W.; Zhen, X.; Meng, C. A novel remote sensing ecological vulnerability index on large scale: A case study of the China-Pakistan Economic Corridor region. Ecol. Indic. 2021, 129, 107955. [Google Scholar] [CrossRef]
  34. Ali, M.; Sajjad, W.; Haleem, A. Climate engineering: A strategic approach to combat environmental potential risks associated with Pak-China Economic corridor (CPEC) Development. Rev. Environ. Health 2021, 36, 143–144. [Google Scholar] [CrossRef] [PubMed]
  35. Du, S.; Wu, R.; Sun, H.; Yan, D.; Xue, J.; Liao, W.; Tuo, Y.; Zhang, W. Projection of Precipitation Extremes and Flood Risk in the China–Pakistan Economic Corridor. Front. Environ. Sci. 2022, 10, 887323. [Google Scholar] [CrossRef]
  36. Kamran, M.; Bian, J.; Li, A.; Lei, G.; Nan, X.; Jin, Y. Investigating Eco-Environmental Vulnerability for China–Pakistan Economic Corridor Key Sector Punjab Using Multi-Sources Geo-Information. ISPRS Int. J. Geo-Inf. 2021, 10, 625. [Google Scholar] [CrossRef]
  37. Young, W.J.; Arif, A.; Tousif, B.; Edoardo, B.; Stephen, D.; Garthwaite, W.R., III; Michael, G.E.; Leb, C.; Lucy, L.; Ian, M.; et al. Pakistan: Getting More from Water; World Bank: Washington, DC, USA, 2019. [Google Scholar]
  38. Zhu, X.; Liu, T.; Xu, K.; Chen, C. The impact of high temperature and drought stress on the yield of major staple crops in northern China. J. Environ. Manag. 2022, 314, 115092. [Google Scholar] [CrossRef]
  39. Bera, B.; Shit, P.K.; Sengupta, N.; Saha, S.; Bhattacharjee, S. Trends and variability of drought in the extended part of Chhota Nagpur plateau (Singbhum Protocontinent), India applying SPI and SPEI indices. Environ. Chall. 2021, 5, 100310. [Google Scholar] [CrossRef]
  40. Makhdoom, A.S.; Shah, A.B.; Sami, K. Pakistan on the roadway to socio-economic development: A comprehensive study of China Pakistan Economic Corridor (CPEC). Gov.-Annu. Res. J. Political Science. 2018, 6, 37–42. [Google Scholar]
  41. Li, Y.; Chen, Y.; Duan, W.; Cao, M.; Qin, J. Spatiotemporal variation in irrigation water requirements in the China–Pakistan Economic Corridor. Sci. Rep. 2022, 12, 17258. [Google Scholar] [CrossRef]
  42. Muhammad, S.; Shiyin, L.; Muhammad, A.; Muhammad, A.; Muhammad, Z.; Sarfraz, H.; Sher, M. Risks of Glaciers Lakes Outburst Flood along China Pakistan Economic Corridor. In Glaciers and the Polar Environment; Masaki, K., Danilo, G., Niccolò, D., Eds.; IntechOpen: Rijeka, Croatia, 2020; Chapter 6. [Google Scholar] [CrossRef]
  43. Asmat, U.; Athar, H.; Nabeel, A.; Latif, M. An AOGCM based assessment of interseasonal variability in Pakistan. Clim. Dyn. 2018, 50, 349–373. [Google Scholar] [CrossRef]
  44. Rahman, K.U.; Shang, S.; Shahid, M.; Wen, Y. Performance Assessment of SM2RAIN-CCI and SM2RAIN-ASCAT Precipitation Products over Pakistan. Remote Sens. 2019, 11, 2040. [Google Scholar] [CrossRef]
  45. Rahman, U.K.; Shang, S.; Shahid, M.; Wen, Y. An Appraisal of Dynamic Bayesian Model Averaging-based Merged Multi-Satellite Precipitation Datasets over Complex Topography and the Diverse Climate of Pakistan. Remote Sens. 2020, 12, 10. [Google Scholar] [CrossRef]
  46. Morsy, M.; Moursy, F.I.; Sayad, T.; Shaban, S. Climatological Study of SPEI Drought Index Using Observed and CRU Gridded Dataset over Ethiopia. Pure Appl. Geophys. 2022, 179, 3055–3073. [Google Scholar] [CrossRef]
  47. Naimark, J.G.; Fiore, A.M.; Jin, X.; Wang, Y.; Klovenski, E.; Braneon, C. Evaluating Drought Responses of Surface Ozone Precursor Proxies: Variations with Land Cover Type, Precipitation, and Temperature. Geophys. Res. Lett. 2021, 48, e2020GL091520. [Google Scholar] [CrossRef]
  48. Chiang, F.; Mazdiyasni, O.; AghaKouchak, A. Amplified warming of droughts in southern United States in observations and model simulations. Sci. Adv. 2018, 4, eaat2380. [Google Scholar] [CrossRef] [PubMed]
  49. Potopová, V.; Štěpánek, P.; Možný, M.; Türkott, L.; Soukup, J. Performance of the standardised precipitation evapotranspiration index at various lags for agricultural drought risk assessment in the Czech Republic. Agric. For. Meteorol. 2015, 202, 26–38. [Google Scholar] [CrossRef]
  50. Chen, H.; Sun, J. Changes in climate extreme events in China associated with warming. Int. J. Climatol. 2015, 35, 2735–2751. [Google Scholar] [CrossRef]
  51. Yao, J.; Zhao, Y.; Chen, Y.; Yu, X.; Zhang, R. Multi-scale assessments of droughts: A case study in Xinjiang, China. Sci. Total Environ. 2018, 630, 444–452. [Google Scholar] [CrossRef]
  52. Mallya, G.; Mishra, V.; Niyogi, D.; Tripathi, S.; Govindaraju, R.S. Trends and variability of droughts over the Indian monsoon region. Weather Clim. Extrem. 2016, 12, 43–68. [Google Scholar] [CrossRef]
  53. Spinoni, J.; Naumann, G.; Vogt, J.V. Pan-European seasonal trends and recent changes of drought frequency and severity. Glob. Planet. Chang. 2017, 148, 113–130. [Google Scholar] [CrossRef]
  54. Shelton, S.; Ogou, F.K.; Pushpawela, B. Spatial-Temporal Variability of Droughts during Two Cropping Seasons in Sri Lanka and Its Possible Mechanisms. Asia-Pac. J. Atmos. Sci. 2021, 58, 127–144. [Google Scholar] [CrossRef]
  55. Dai, A. Drought under global warming: A review. Wiley Interdiscip. Rev. Clim. Chang. 2011, 2, 45–65. [Google Scholar] [CrossRef]
  56. Sheffield, J.; Wood, E.F. Projected changes in drought occurrence under future global warming from multi-model, multi-scenario, IPCC AR4 simulations. Clim. Dyn. 2008, 31, 79–105. [Google Scholar] [CrossRef]
  57. Sheffield, J.; Wood, E.F.; Roderick, M.L. Little change in global drought over the past 60 years. Nature 2012, 491, 435–438. [Google Scholar] [CrossRef]
  58. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef]
  59. Guo, H.; Bao, A.; Liu, T.; Jiapaer, G.; Ndayisaba, F.; Jiang, L.; Kurban, A.; De Maeyer, P. Spatial and temporal characteristics of droughts in Central Asia during 1966–2015. Sci. Total Environ. 2018, 624, 1523–1538. [Google Scholar] [CrossRef]
  60. Bae, S.; Lee, S.-H.; Yoo, S.-H.; Kim, T. Analysis of Drought Intensity and Trends Using the Modified SPEI in South Korea from 1981 to 2010. Water 2018, 10, 327. [Google Scholar] [CrossRef]
  61. Wu, J.; Chen, X.; Yao, H.; Gao, L.; Chen, Y.; Liu, M. Non-linear relationship of hydrological drought responding to meteorological drought and impact of a large reservoir. J. Hydrol. 2017, 551, 495–507. [Google Scholar] [CrossRef]
  62. Yue, Y.; Liu, H.; Mu, X.; Qin, M.; Wang, T.; Wang, Q.; Yan, Y. Spatial and temporal characteristics of drought and its correlation with climate indices in Northeast China. PLoS ONE 2021, 16, e0259774. [Google Scholar] [CrossRef] [PubMed]
  63. Yacoub, E.; Tayfur, G. Spatial and temporal of variation of meteorological drought and precipitation trend analysis over whole Mauritania. J. Afr. Earth Sci. 2020, 163, 103761. [Google Scholar] [CrossRef]
  64. Mekonen, A.A.; Berlie, A.B.; Ferede, M.B. Spatial and temporal drought incidence analysis in the northeastern highlands of Ethiopia. Geoenvironmental Disasters 2020, 7, 10. [Google Scholar] [CrossRef]
  65. Wang, S.-Y.; Davies, R.E.; Huang, W.-R.; Gillies, R.R. Pakistan’s two-stage monsoon and links with the recent climate change. J. Geophys. Res. Atmos. 2011, 116, 15760. [Google Scholar] [CrossRef]
  66. Bhatti, A.S.; Wang, G.; Ullah, W.; Ullah, S.; Fiifi Tawia Hagan, D.; Kwesi Nooni, I.; Lou, D.; Ullah, I. Trend in Extreme Precipitation Indices Based on Long Term In Situ Precipitation Records over Pakistan. Water 2020, 12, 797. [Google Scholar] [CrossRef]
  67. Hunt, K.M.R.; Turner, A.G.; Shaffrey, L.C. Falling Trend of Western Disturbances in Future Climate Simulations. J. Clim. 2019, 32, 5037–5051. [Google Scholar] [CrossRef]
  68. Zhang, H.; Wang, L. Analysis of the variation in potential evapotranspiration and surface wet conditions in the Hancang River Basin, China. Sci. Rep. 2021, 11, 8607. [Google Scholar] [CrossRef]
  69. Adnan, S.; Ullah, K.; Khan, A.H.; Gao, S. Meteorological impacts on evapotranspiration in different climatic zones of Pakistan. J. Arid Land 2017, 9, 938–952. [Google Scholar] [CrossRef]
  70. Smakhtin, V.U.; Schipper, E.L.F. Droughts: The impact of semantics and perceptions. Water Policy 2008, 10, 131–143. [Google Scholar] [CrossRef]
  71. Zeri, M.; Williams, K.; Cunha, A.P.M.A.; Cunha-Zeri, G.; Vianna, M.S.; Blyth, E.M.; Marthews, T.R.; Hayman, G.D.; Costa, J.M.; Marengo, J.A.; et al. Importance of including soil moisture in drought monitoring over the Brazilian semiarid region: An evaluation using the JULES model, in situ observations, and remote sensing. Clim. Resil. Sustain. 2022, 1, e7. [Google Scholar] [CrossRef]
  72. Ahmad, I.; Ambreen, R.; Sultan, S.; Sun, Z.; Deng, W. Spatial-Temporal Variations in January Precipitation over the Period of 1950–2000 in Pakistan and Possible Links with Global Teleconnections: Geographical Perspective. Am. J. Clim. Chang. 2014, 3, 378–387. [Google Scholar] [CrossRef]
  73. Ahmad, S.; Hussain, Z.; Qureshi, A.S.; Majeed, R.; Saleem, M. Drought Mitigation in Pakistan: Current Status and Options for Future Strategies; Working Paper 85; International Water Management Institute (IWMI): Colombo, Sri Lanka, 2004. [Google Scholar]
  74. Rawat, S.; Ganapathy, A.; Agarwal, A. Drought characterization over Indian sub-continent using GRACE-based indices. Sci. Rep. 2022, 12, 15432. [Google Scholar] [CrossRef]
  75. Shah, S.M.H.; Mustaffa, Z.; Teo, F.Y.; Imam, M.A.H.; Yusof, K.W.; Al-Qadami, E.H.H. A review of the flood hazard and risk management in the South Asian Region, particularly Pakistan. Sci. Afr. 2020, 10, e00651. [Google Scholar] [CrossRef]
  76. Jamro, S.; Channa, F.N.; Dars, G.H.; Ansari, K.; Krakauer, N.Y. Exploring the Evolution of Drought Characteristics in Balochistan, Pakistan. Appl. Sci. 2020, 10, 913. [Google Scholar] [CrossRef]
  77. Naz, F.; Dars, G.H.; Ansari, K.; Jamro, S.; Krakauer, N.Y. Drought Trends in Balochistan. Water 2020, 12, 470. [Google Scholar] [CrossRef]
  78. Ashraf, M.; Arshad, A.; Patel, P.M.; Khan, A.; Qamar, H.; Siti-Sundari, R.; Ghani, M.U.; Amin, A.; Babar, J.R. Quantifying climate-induced drought risk to livelihood and mitigation actions in Balochistan. Nat. Hazards 2021, 109, 2127–2151. [Google Scholar] [CrossRef]
  79. Hanif, M.; Khan, A.; Adnan, S. Latitudinal precipitation characteristics and trends in Pakistan. J. Hydrol. 2013, 492, 266–272. [Google Scholar] [CrossRef]
  80. Stahle, D.W. Anthropogenic megadrought. Science 2020, 368, 238–239. [Google Scholar] [CrossRef]
  81. Wang, Y.; Wang, S.; Zhao, W.; Liu, Y. The increasing contribution of potential evapotranspiration to severe droughts in the Yellow River basin. J. Hydrol. 2022, 605, 127310. [Google Scholar] [CrossRef]
  82. Ahmed, K.; Shahid, S.; Chung, E.-S.; Nawaz, N.; Khan, N.; Rasheed, B. Divergence of potential evapotranspiration trends over Pakistan during 1967–2016. Theor. Appl. Climatol. 2020, 141, 215–227. [Google Scholar] [CrossRef]
  83. Goroshi, S.; Pradhan, R.; Singh, R.P.; Singh, K.K.; Parihar, J.S. Trend analysis of evapotranspiration over India: Observed from long-term satellite measurements. J. Earth Syst. Sci. 2017, 126, 113. [Google Scholar] [CrossRef]
  84. Basso, B.; Martinez-Feria, R.A.; Rill, L.; Ritchie, J.T. Contrasting long-term temperature trends reveal minor changes in projected potential evapotranspiration in the US Midwest. Nat. Commun. 2021, 12, 1476. [Google Scholar] [CrossRef]
  85. Dong, W.; Li, C.; Hu, Q.; Pan, F.; Bhandari, J.; Sun, Z. Potential Evapotranspiration Reduction and Its Influence on Crop Yield in the North China Plain in 1961–2014. Adv. Meteorol. 2020, 2020, 3691421. [Google Scholar] [CrossRef]
  86. Sharma, S.; Mujumdar, P. Increasing frequency and spatial extent of concurrent meteorological droughts and heatwaves in India. Sci. Rep. 2017, 7, 15582. [Google Scholar] [CrossRef] [PubMed]
  87. Iqbal, M.A.; Penas, A.; Cano-Ortiz, A.; Kersebaum, K.C.; Herrero, L.; del Río, S. Analysis of recent changes in maximum and minimum temperatures in Pakistan. Atmos. Res. 2016, 168, 234–249. [Google Scholar] [CrossRef]
  88. Nawaz, Z.; Li, X.; Chen, Y.; Guo, Y.; Wang, X.; Nawaz, N. Temporal and Spatial Characteristics of Precipitation and Temperature in Punjab, Pakistan. Water 2019, 11, 1916. [Google Scholar] [CrossRef]
  89. Shelton, S.; Pushpawela, B.; Liyanage, G. The long-term trend in the diurnal temperature range over Sri Lanka from 1985 to 2017 and its association with total cloud cover and rainfall. J. Atmos. Sol.-Terr. Phys. 2022, 227, 105810. [Google Scholar] [CrossRef]
  90. Sheikh, M.M.; Manzoor, N.; Ashraf, J.; Adnan, M.; Collins, D.; Hameed, S.; Manton, M.J.; Ahmed, A.U.; Baidya, S.K.; Borgaonkar, H.P.; et al. Trends in extreme daily rainfall and temperature indices over South Asia. Int. J. Climatol. 2015, 35, 1625–1637. [Google Scholar] [CrossRef]
  91. Kim, D.; Rhee, J. A drought index based on actual evapotranspiration from the Bouchet hypothesis. Geophys. Res. Lett. 2016, 43, 10–277. [Google Scholar] [CrossRef]
  92. Vicente-Serrano, S.M.; Lopez-Moreno, J.-I.; Beguería, S.; Lorenzo-Lacruz, J.; Sanchez-Lorenzo, A.; García-Ruiz, J.M.; Azorin-Molina, C.; Morán-Tejeda, E.; Revuelto, J.; Trigo, R. Evidence of increasing drought severity caused by temperature rise in southern Europe. Environ. Res. Lett. 2014, 9, 044001. [Google Scholar] [CrossRef]
  93. Howells, M.; Hermann, S.; Welsch, M.; Bazilian, M.; Segerström, R.; Alfstad, T.; Gielen, D.; Rogner, H.; Fischer, G.; Van Velthuizen, H. Integrated analysis of climate change, land-use, energy and water strategies. Nat. Clim. Chang. 2013, 3, 621–626. [Google Scholar] [CrossRef]
  94. Mishra, A.; Liu, S.C. Changes in precipitation pattern and risk of drought over India in the context of global warming. J. Geophys. Res. Atmos. 2014, 119, 7833–7841. [Google Scholar] [CrossRef]
  95. Brunner, M.I.; Swain, D.L.; Gilleland, E.; Wood, A.W. Increasing importance of temperature as a contributor to the spatial extent of streamflow drought. Environ. Res. Lett. 2021, 16, 024038. [Google Scholar] [CrossRef]
  96. Miralles, D.G.; Gentine, P.; Seneviratne, S.I.; Teuling, A.J. Land–atmospheric feedbacks during droughts and heatwaves: State of the science and current challenges. Ann. N. Y. Acad. Sci. 2019, 1436, 19–35. [Google Scholar] [CrossRef]
  97. Konapala, G.; Mishra, A.K.; Wada, Y.; Mann, M.E. Climate change will affect global water availability through compounding changes in seasonal precipitation and evaporation. Nat. Commun. 2020, 11, 3044. [Google Scholar] [CrossRef] [PubMed]
  98. Arnell, N.W.; Gosling, S.N. The impacts of climate change on river flood risk at the global scale. Clim. Chang. 2016, 134, 387–401. [Google Scholar] [CrossRef]
  99. Challinor, A.J.; Watson, J.; Lobell, D.B.; Howden, S.M.; Smith, D.R.; Chhetri, N. A meta-analysis of crop yield under climate change and adaptation. Nat. Clim. Chang. 2014, 4, 287–291. [Google Scholar] [CrossRef]
  100. Chen, Z.; Wang, W.; Fu, J. Vegetation response to precipitation anomalies under different climatic and biogeographical conditions in China. Sci. Rep. 2020, 10, 830. [Google Scholar] [CrossRef]
  101. Ahmed, K.; Shahid, S.; Wang, X.; Nawaz, N.; Khan, N. Spatiotemporal changes in aridity of Pakistan during 1901–2016. Hydrol. Earth Syst. Sci. 2019, 23, 3081–3096. [Google Scholar] [CrossRef]
  102. Syed, A.; Raza, T.; Bhatti, T.T.; Eash, N.S. Climate Impacts on the agricultural sector of Pakistan: Risks and solutions. Environ. Chall. 2022, 6, 100433. [Google Scholar] [CrossRef]
  103. Zekollari, H.; Huss, M.; Farinotti, D. Modelling the future evolution of glaciers in the European Alps under the EURO-CORDEX RCM ensemble. Cryosphere 2019, 13, 1125–1146. [Google Scholar] [CrossRef]
  104. Kääb, A.; Berthier, E.; Nuth, C.; Gardelle, J.; Arnaud, Y. Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas. Nature 2012, 488, 495–498. [Google Scholar] [CrossRef]
  105. Safdar, F.; Khokhar, M.F.; Arshad, M.; Adil, I.H. Climate Change Indicators and Spatiotemporal Shift in Monsoon Patterns in Pakistan. Adv. Meteorol. 2019, 2019, 8281201. [Google Scholar] [CrossRef]
  106. Katzenberger, A.; Schewe, J.; Pongratz, J.; Levermann, A. Robust increase of Indian monsoon rainfall and its variability under future warming in CMIP6 models. Earth Syst. Dynam. 2021, 12, 367–386. [Google Scholar] [CrossRef]
  107. Aamir, E.; Hassan, I. The impact of climate indices on precipitation variability in Baluchistan, Pakistan. Tellus A Dyn. Meteorol. Oceanogr. 2020, 72, 1833584. [Google Scholar] [CrossRef]
  108. Niranjan Kumar, K.; Rajeevan, M.; Pai, D.S.; Srivastava, A.K.; Preethi, B. On the observed variability of monsoon droughts over India. Weather Clim. Extrem. 2013, 1, 42–50. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of the China–Pakistan Economic corridor routes and different climate regions. The black color lines demarcate hyper-arid, arid, humid, and glacial regions named R-I, R-II, R-III, and R-IV, respectively.
Figure 1. Location of the China–Pakistan Economic corridor routes and different climate regions. The black color lines demarcate hyper-arid, arid, humid, and glacial regions named R-I, R-II, R-III, and R-IV, respectively.
Climate 11 00045 g001
Figure 2. Annual cycle of the monthly total rainfall (unit: mm) in (a) Pakistan (c) hyper-arid (R-I), (d) arid (R-II), (e) humid (R-III), and (f) glacial (R-IV) regions for the 1980–2018 period. The blue boxes demarcate the months of the ISM season (from June to September). The long-term mean is shown by black dots, while the lines of the upper and lower whisk represent maximum and minimum rainfall, respectively. 5th (lower margin of the box), and 95th (upper margin of the box) percentiles of the rainfall are shown. The x-axis represents the month of the year starting from December. The four climate regions are shown in (b).
Figure 2. Annual cycle of the monthly total rainfall (unit: mm) in (a) Pakistan (c) hyper-arid (R-I), (d) arid (R-II), (e) humid (R-III), and (f) glacial (R-IV) regions for the 1980–2018 period. The blue boxes demarcate the months of the ISM season (from June to September). The long-term mean is shown by black dots, while the lines of the upper and lower whisk represent maximum and minimum rainfall, respectively. 5th (lower margin of the box), and 95th (upper margin of the box) percentiles of the rainfall are shown. The x-axis represents the month of the year starting from December. The four climate regions are shown in (b).
Climate 11 00045 g002
Figure 3. Spatial distribution of rainfall (mm/month) during the (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) Season over Pakistan for 1980–2018. Temperature (°C), potential evapotranspiration (PET; mm/month), and soil moisture (Soilm; mm/month) in four climate seasons are shown in (eh), (il) and (mp), respectively. Hyper-arid (R-I), arid (R-II), humid (R-III), and glacial (R-IV) regions are demarcated and named in (a).
Figure 3. Spatial distribution of rainfall (mm/month) during the (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) Season over Pakistan for 1980–2018. Temperature (°C), potential evapotranspiration (PET; mm/month), and soil moisture (Soilm; mm/month) in four climate seasons are shown in (eh), (il) and (mp), respectively. Hyper-arid (R-I), arid (R-II), humid (R-III), and glacial (R-IV) regions are demarcated and named in (a).
Climate 11 00045 g003
Figure 4. Temporal variation of the SPEI at different timescales (1 to 12 months) for (a) R-I, (b) R-II, (c) R-II, and (d) R-IV from 1985 to 2015. The R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively. The SPEI-n below −1 (or greater than 1) is selected as a threshold level for the drought (or wet) event.
Figure 4. Temporal variation of the SPEI at different timescales (1 to 12 months) for (a) R-I, (b) R-II, (c) R-II, and (d) R-IV from 1985 to 2015. The R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively. The SPEI-n below −1 (or greater than 1) is selected as a threshold level for the drought (or wet) event.
Climate 11 00045 g004
Figure 5. Spatial distribution of seasonal drought frequency (%) for (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) seasons. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively, for the 1980–2018 periods.
Figure 5. Spatial distribution of seasonal drought frequency (%) for (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) seasons. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively, for the 1980–2018 periods.
Climate 11 00045 g005
Figure 6. Spatial distribution of Mann–Kendal trend (per decade) of Standardized Precipitation Evapotranspiration Index (SPEI) during the (a) DJF (December to February; SPEI-3 at February), (b) MAM (March to May; SPEI-3 at May), (c) JJAS (June to September SPEI-4 at September), and (d) ON (October to November; SPEI-2 at November) Season over Pakistan for 1980–2018. The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.
Figure 6. Spatial distribution of Mann–Kendal trend (per decade) of Standardized Precipitation Evapotranspiration Index (SPEI) during the (a) DJF (December to February; SPEI-3 at February), (b) MAM (March to May; SPEI-3 at May), (c) JJAS (June to September SPEI-4 at September), and (d) ON (October to November; SPEI-2 at November) Season over Pakistan for 1980–2018. The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.
Climate 11 00045 g006
Figure 7. Spatial distribution of Mann–Kendal trend (mm/decade) of rainfall during the (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) Season over the China Pakistan Economic Corridor for 1980–2018. The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.
Figure 7. Spatial distribution of Mann–Kendal trend (mm/decade) of rainfall during the (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) Season over the China Pakistan Economic Corridor for 1980–2018. The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.
Climate 11 00045 g007
Figure 8. Spatial distribution of Mann–Kendal trend (mm/decade) of Potential Evapotranspiration (PET) during the (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) Season over the China Pakistan Economic Corridor for 1980–2018. The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid and glacial regions, respectively.
Figure 8. Spatial distribution of Mann–Kendal trend (mm/decade) of Potential Evapotranspiration (PET) during the (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) Season over the China Pakistan Economic Corridor for 1980–2018. The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid and glacial regions, respectively.
Climate 11 00045 g008
Figure 9. Spatial distribution of Mann–Kendal trend (°C/decade) of maximum temperature (Tmax) during the (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) Season over the China Pakistan Economic Corridor for 1980–2018. Lower panels (eh) are the same as top panels but for the minimum temperature (Tmin). The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.
Figure 9. Spatial distribution of Mann–Kendal trend (°C/decade) of maximum temperature (Tmax) during the (a) DJF (December to February), (b) MAM (March to May), (c) JJAS (June to September), and (d) ON (October to November) Season over the China Pakistan Economic Corridor for 1980–2018. Lower panels (eh) are the same as top panels but for the minimum temperature (Tmin). The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.
Climate 11 00045 g009
Figure 10. Spatial distribution of correlation between rainfall and SPEI for (a) MAM (March to May), (b) JJAS (June to September), (c) ON (October to November), and (d) DJF (December to February) seasons over the China Pakistan Economic Corridor for 1980–2018. The lower panel (eh) is the same as the top panels but for PET. The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.
Figure 10. Spatial distribution of correlation between rainfall and SPEI for (a) MAM (March to May), (b) JJAS (June to September), (c) ON (October to November), and (d) DJF (December to February) seasons over the China Pakistan Economic Corridor for 1980–2018. The lower panel (eh) is the same as the top panels but for PET. The hatched area shows a statistically significant trend at a 90% confidence level. R-I, R-II, R-III, and R-IV represent the hyper-arid, arid, humid, and glacial regions, respectively.
Climate 11 00045 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shelton, S.; Dixon, R.D. Long-Term Seasonal Drought Trends in the China-Pakistan Economic Corridor. Climate 2023, 11, 45. https://doi.org/10.3390/cli11020045

AMA Style

Shelton S, Dixon RD. Long-Term Seasonal Drought Trends in the China-Pakistan Economic Corridor. Climate. 2023; 11(2):45. https://doi.org/10.3390/cli11020045

Chicago/Turabian Style

Shelton, Sherly, and Ross D. Dixon. 2023. "Long-Term Seasonal Drought Trends in the China-Pakistan Economic Corridor" Climate 11, no. 2: 45. https://doi.org/10.3390/cli11020045

APA Style

Shelton, S., & Dixon, R. D. (2023). Long-Term Seasonal Drought Trends in the China-Pakistan Economic Corridor. Climate, 11(2), 45. https://doi.org/10.3390/cli11020045

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop