Recent Changes in Drought Events over South Asia and Their Possible Linkages with Climatic and Dynamic Factors

: South Asia is home to one of the fastest-growing populations in Asia, and human activities are leaving indelible marks on the land surface. Yet the likelihood of successive observed droughts in South Asia (SA) and its four subregions (R-1: semi-arid, R-2: arid, R-3: subtropical wet, and R-4: tropical wet and dry) remains poorly understood. Using the state-of-the-art self-calibrated Palmer Drought Severity Index (scPDSI), we examined the impact of different natural ocean variability modes on the evolution, severity, and magnitude of observed droughts across the four subregions that have distinct precipitation seasonality and cover key breadbaskets and highly vulnerable populations. The study revealed that dryness had signiﬁcantly increased in R-1, R-2, and R-4 during 1981–2020. Temporal analysis revealed an increase in drought intensity for R-1 and R-4 since the 2000s, while a mixed behavior was observed in R-2 and R-3. Moreover, most of the sub-regions witnessed a substantial upsurge in annual precipitation, but a signiﬁcant decrease in vapor pressure deﬁcit (VPD) during 1981–2020. The increase in precipitation and the decline in VPD partially contributed to a signiﬁcant rise in Normalized Difference Vegetation Index (NDVI) and a decrease in dryness. In contrast, a strong positive correlation was found between drought index and precipitation, and NDVI across R-1, R-2, and R-4, whereas temperature and VPD exhibited a negative correlation over these regions. No obvious link was detected with El-Niño Southern Oscillation (ENSO) events, or Indian Ocean Dipole (IOD) and drought evolution, as explored for certain regions of SA. The ﬁndings showed the possibility that the precipitation changes over these regions had an insigniﬁcant relationship with ENSO, IOD, and drought onset. Thus, the study results highlight the need for considering interactions within the longer climate system in describing observed drought risks rather than aiming at drivers from an individual perspective.

activities to drought occurrence and severity in a changing climate must, therefore, be recognized, both to prepare for their impacts and to guide a more sustainable management of water resources.
Here, we examined the long-term variations in the trend of drought events and their possible linkages with vegetation and dynamic drivers across SA and its subregions. The present study also investigated the influence of PET, vapor pressure deficit (VPD), and vegetation estimation approaches, on the changing features of droughts, particularly the effects of extreme droughts (dryness/wetness) due to the diverse climate of SA. The remainder of this manuscript is structured as follows: Section 2 highlights the study domain and datasets, whereas Section 3 provides the analytical methodology. Section 4 discusses the feasibility and implications of our results. Finally, the conclusions and limitations are presented in Section 5.

Study Domain
The study domain comprises eight countries: India, Pakistan, Afghanistan, Bangladesh, Nepal, Bhutan, Maldives, and Sri Lanka. SA is situated in the southern part of the Asian continent with geographical coordinates of 5.5 • -39.5 • N latitudes and 60.5 • -99.5 • E longitudes. It has a multifaceted topography covering an area of 5,135,752 km 2 (see Figure 1a). The region spans a diversity of climatic and geographical types, i.e., arid, drylands and deserts, tropical and subtropical, humid, alpine, and mountain climates and regions ( Figure 1b).
Climatologically, arid to semi-humid climate types characterize the region. Precipitation in SA is mostly dependent on monsoon during the summer season, where more than 75% of the annual precipitation is received [3]. Temperatures vary from around 5 • C in higher altitude regions during the cold season to around 50 • C on the lowland plain during the warm season [7]. SA and its subregions were declared as a hotspot for climate change, and play a vital role in neighboring countries' socioeconomic development. The climate of SA is heavily affected by both the SASM and the southern branch of the mid-latitude westerlies that bifurcate around the western edge of the Tibetan Plateau [24]. This multifaceted climate background sets the stage for complex compound interactions among temperature, precipitation, runoff, evapotranspiration, and glacial melt, all of which contribute to water resource availability and its variability in the region [25]. Extreme rainfall during the monsoon season can trigger floods and landslides, while declines in winter precipitation are likely enhance water scarcity and drought risk within the region [26,27].

Data Sets
The present study used a 30 m high-resolution digital elevation model (DEM) from the Shuttle Radar Topography Mission (SRTM) of the US National Aeronautics and Space Administration (NASA), to show the topography of the target region ( Figure 1a). The land use/cover change (LULC) map for 2018 at a 20 m high resolution was obtained from the European Space Agency Climate Change Initiative (ESACCI) [28], see Figure 1b. For climate variables, the high resolution (0.25 • ) daily precipitation was acquired from the Climate Hazards Group Infrared Precipitation with Stations (CHIRPS) version 2 product, from 1981 to the present [29]. The daily CHIRPS precipitation dataset has been used widely for different meteorological applications, for instance, the assessment of precipitation extremes and their features in many parts of the globe [30]. CHIRPS merges the satellite-based precipitation estimates with ground observations and models of terrain-based precipitation alteration to provide high resolution, spatially complete, and long-term continuous data (1981 to present). The dataset benefits from rain observation-based datasets characterized by sparser station density and remotely sensed rainfall products with limited temporal coverage. The observed minimum and maximum temperature products were collected from the Climate Prediction Center (CPC) at a spatial resolution of 0.5 • [31]. This study employed CHIRPS and CPC datasets to illustrate the observed climatology of the region; however, both datasets strictly follow the World Meteorological Organization (WMO) Figure 1. Spatial patterns of topography and elevation (a), and land use land cover (LULC) (b), across South Asia (SA) and its four subregions; (c) and (d) represent the annual precipitation and mean temperature distributions of spatial anomalies during 1980-2020; (e) and (f) demonstrate the annual and seasonal trends of precipitation (mm/year) and temperature (°C/year), and the regional averaged annual precipitation and mean temperature time series over SA. The blue and red solid lines in (e) and (f) characterize the precipitation and temperature changes over SA. The horizontal lines (black) indicate the average annual precipitation and mean temperature, while the dotted lines (blue and red) are the linear trends of annual precipitation and mean temperature, respectively. The South Asia (SA) and its four subregions; (c,d) represent the annual precipitation and mean temperature distributions of spatial anomalies during 1980-2020; (e,f) demonstrate the annual and seasonal trends of precipitation (mm/year) and temperature ( • C/year), and the regional averaged annual precipitation and mean temperature time series over SA. The blue and red solid lines in (e,f) characterize the precipitation and temperature changes over SA. The horizontal lines (black) indicate the average annual precipitation and mean temperature, while the dotted lines (blue and red) are the linear trends of annual precipitation and mean temperature, respectively. The asterisk *, **, *** represent trends passing the 95%, 99%, 99.9% significance levels, respectively.
To explore the spatiotemporal distributions of vegetation and their responses to dryness/wetness conditions, we used the U.S. NOAA Advanced Very High-Resolution Radiometer (AVHRR) Normalized Difference Vegetation Index (NDVI) daily product at 0.083 • spatial resolution [34]. The NDVI data employed in this study have a temporal coverage from 1981 to the present. This study used monthly precipitation, air temperature (maximum and minimum), relative humidity, wind components, latent and sensible heat fluxes, and surface air pressure from ERA5 to estimate the drought index from 1979 to 2020. The recently developed European Centre for Medium-Range Weather Forecasts (ECMWF) a global reanalysis product (ERA-5) was used in the current study [35]. Available water content (AWC) of soil was collected from the Harmonized World Soil Database (HWSD) [36].

Self-Calibrated Palmer Drought Severity Index (scPDSI)
We used the scPDSI as a state-of-the-art drought index to compute the intensity and frequency of extreme drought events over SA and its subregions. The scPDSI was developed using a simple two-layer water balance model [37], considering water supply (precipitation) and atmospheric water demand (PET). The precipitation, PET, and AWC of soil employed in scPDSI calculation need to be assessed under climatically suitable conditions [20]. The present study used the Penman-Monteith method for estimating the PET [19]. Long-term scPDSI series for the period 1981-2020 were considered using the tool provided by the PDSI [38], for which readers are referred to Jacobi et al. (2013) [38] and van der Schrier et al. (2013) [39]. According to the scPDSI dryness/wetness classification (scPDSI values range from -10 (dryness) to +10 (wetness)), values below -3 indicate severe to extreme droughts. However, in our analysis, minimum and maximum scPDSI values are limited to ±5 characterizing dryness and wetness conditions, respectively. A brief description of drought classification is presented in Table 1.

k-Means Algorithm
The present study employed a k-means algorithm for scPDSI data to better understand the drought evolution and frequency in the distinct subregions of SA. The k-means algorithm has various capabilities to minimize variability with the number of k clusters and maximize the variability of each centroid among all the data series [40,41]. The kmeans algorithm provides a better partition result when many grids, for instance, more than 5600 points, are considered. Various studies prefer to use the k-means algorithm while dividing the diverse climate into subregions [41]. Therefore, for a detailed assessment of droughts, we selected four distinct subregions, namely, semi-arid (R-1), arid (R-2), subtropical wet (R-3), and tropical wet and dry (R-4) subregions across SA.

Trends Analysis
To estimate trends, this study used the nonparametric Mann-Kendall (MK) trend test, which is widely used to identify trends/changes in auto-correlated SPEI data for a hydro-meteorological time-series [42,43]. The MK test does not require normal distribution and is effective against outliers in the time-series. The presence of auto-correlation in the dataset could probably affect the trend detection. Thus, a secular trend estimator was utilized to deduct from the original time-series while estimating auto-correlation [44]. The nonparametric modified-MK test was then employed in the scPDSI time series and the trend estimate, with p-values < 0.05 were considered statistically significant. The current study used the nonparametric Theil Sen s Slope (TS) test to assess the trend magnitude [45], while the Pettit test was also adopted to check the abrupt change points [46]. The TS estimator is a non-parametric statistical method for estimating the slope of a time series. Compared with other methods, the TS method shows strong robustness to the outliers over time series. Drought is a multifaceted natural hazard that changes with space and time. Thus, an inclusive understanding of drought requires studying its spatial variability. This study assessed the variation in magnitude for scPDSI data using the TS estimator for all the scPDSI grid points over SA. The black color dots were used to show the impact of the trend at a 0.05% significance level assessed using the MK test at each grid point in Figure 5.
To explore the interannual variability, this study also assessed precipitation, temperature, VPD, and NDVI anomalies, considering the spatially averaged annual values at each grid point. In our analysis, the negative values refer to dryness and vegetation stress, while positive values correspond to wetness and/or green vegetation. To enhance the robustness of our findings, we used a spatial correlation of various climate factors and drought index, i.e., precipitation, temperature, VPD, and NDVI, with associated dynamic drivers, i.e., ENSO and IOD, of the wetness and dryness coupling for the period of 1981-2020.

Vapor Pressure Deficit
In addition to precipitation, the other significant influencing factor of the scPDSI is the surface VPD, which was estimated based on the saturated vapor pressure (SVP) and the actual vapor pressure (AVP). Previous studies have highlighted that both drivers of surface VPD and soil moisture should be accounted for to explain and simulate recent observed changes in hydro-meteorological disasters (i.e., droughts) [47]. The VPD was calculated as the difference between SVP and AVP which were estimated based on the method given by Allen et al. (1994) [48]. The relative humidity used for the AVP calculation was acquired by employing Tetens's formula [49] with the factors constructed on saturation conditions over water [50].

Observed Annual and Seasonal Trends of Precipitation and Temperature
Figure 1c-f show the spatiotemporal characteristics of long-term precipitation and temperature (i.e., annual and seasonal) changes over the four distinct subregions of SA from 1981 to 2020. The precipitation anomaly exhibited great variability from region to region due to the diverse climate of SA (Figure 1c). The annual precipitation anomaly was obvious in R-1 with an amount of −350 mm, while relatively more precipitation was seen in R-3, followed by R-2 and R-4. In the case of seasonal variations, precipitation in the Spring season, Summer season, and Winter season exhibited decreasing trends, which were −0.2 mm/year, −0.5 mm/year, and −0.3 mm/year, respectively, whereas Autumn showed Remote Sens. 2022, 14, 3219 7 of 20 an increasing trend with the amount of 0.2 mm/year. The annual trend of precipitation in SA was −0.9 mm/year ( Figure 1e). These trends in seasonal precipitation show that the summer and winter seasons largely contributed to the change in annual precipitation. On the other hand, not all the trends in precipitation were statistically significant, except for that of the winter season.
In terms of climatological mean, R-4 had a relatively higher spatial temperature anomaly of 9.0 • C, followed by R-2 (6.0 • C), and R-3 (8.0 • C), while less temperature anomaly was observed in R-1 and the foothills of Himalayans, at a rate of 3.0 • C (Figure 1d). The annual and seasonal mean temperatures showed strong and significant warming in the recent past across SA and its subregions (Figure 1d,f). Seasonal mean temperature exhibited the largest increasing trend during the spring season at a rate of 0.03 • C/year, and summer saw the smallest warming trend. Furthermore, all the upward trends were statistically significant at the 99% level ( Figure 1f). Interestingly, the present study also witnessed increasing trends of warm temperature events in the southwestern and northeastern parts of the regions, which agreed with the findings of previous studies in the same regions [51,52]. The upsurge of frequency, intensity, and extent of daily warm extreme temperature and decline of extreme cold temperature may occur in parts of the world, resulting in the escalation of length and intensity of warm period or a heatwave in most land areas [3,5,53,54]. Figure 2 presents the spatial patterns of drought and its main influencing factors that include scPDSI (2a), PET (2b), VPD (2c), and NDVI (2d) across SA for the period 1981-2020, while regional percentage changes are shown in Figure 2e, respectively. For scPDSI (Figure 2a), the drying tendency was pronounced, with maximum intensity over R-2, R-3, and R-4, where a noticeable rise in mean temperature and decline in precipitation was observed in those regions (Figure 1c,d). Interestingly, scPDSI indicated a significant drying trend over large areas of all subregions, especially in southern Pakistan, eastern Afghanistan, northern and inner parts of India, as well as the eastern part of Bangladesh. The study region is vulnerable to climate change, and a sharp temperature rise will further accelerate the melting of glaciers in the region [52]. Recent studies have suggested that if global warming was measured below 1.5 • C, the temperature upsurge in the foothills of Himalayan and Thar deserts regions could be at least 0.8 • C higher than the global average [55,56]. For PET, an obvious increasing pattern was observed over R-3 and R-4, followed by R-1 and R-2, suggesting more drying trends in these regions. On the other hand, the dry spatial anomalies for PET (10 to 15 mm/month) were obvious in most SA subregions except R-1. The results showed that the sub-tropical westerly jets frequently appeared during dry anomalies, which are expected to bring dry and cold air into the northern parts of SA, resulting in lower precipitation and higher evaporation. Therefore, the regions observed a decline in precipitation and a rise in PET in the recent past, implying a strong effect of regional oceanic and atmospheric circulation patterns, i.e., monsoon onset and easterly/westerly systems [11,57].

Anomalous Changes in Different scPDSI Drivers
To document the reason beyond the scPDSI drying pattern, we further assessed the spatial distribution of VPD as a potential indicator for the period 1981 to 2020 (Figure 2c). We found that VPD had increased across most SA subregions, particularly in R-3 and R-4, due to a significant increase in temperature. The spatially anomalously high VPD contents were obvious over southern SA, including India, Sri Lanka, and Bangladesh, with enhanced VPD due to higher air temperature. In the case of NDVI anomalies (Figure 2d), the positive NDVI values were obvious over the arid, semi-arid, and subtropical regions in most drought-prone regions. This finding is consistent with spatially drying conditions observed in previous studies [3,58], which described that water scarcity was one of the main influencing factors of spatial variations in vegetation greenness over these regions. dry spatial anomalies for PET (10 to 15 mm/month) were obvious in most SA subre except R-1. The results showed that the sub-tropical westerly jets frequently appe during dry anomalies, which are expected to bring dry and cold air into the northern of SA, resulting in lower precipitation and higher evaporation. Therefore, the region served a decline in precipitation and a rise in PET in the recent past, implying a s effect of regional oceanic and atmospheric circulation patterns, i.e., monsoon onse easterly/westerly systems [11,57].  Panel (e) indicates the regional changes in all those variables over SA and its subregions. It should be noted that we estimate the regional changes as a proportion of the areas with (significant) reduction (increase) of the indices over the period.
Regional percentage changes in influencing factors across different subregions of SA were presented using boxplots from 1981 to 2020 (Figure 2e). PET exhibited the highest upsurge, followed by NDVI and VPD with a relatively slight change. The regional temporal changes further justified the spatial results for influencing factors to drought conditions in these regions. However, PET and NDVI observed the highest percentage changes, while slight variations were observed in VPD, suggesting having fewer effects than the rest of the variables. Conversely, the relative percentage changes in influencing factors demonstrated that the observed drought conditions were subjected to an upsurge by 45%, 65%, and 85% for PET and NDVI with extreme severe drought events in recent decades, respectively. Overall, it can be concluded from the results that the influencing factors (e.g., temperature, VPD, and NDVI) rapidly increased over major parts of SA, which could be more vulnerable to extreme droughts and temperature extreme events in the near future, due to rapid urbanization, population growth, and human-induced activities.

Observed Spatiotemporal Changes in Drought Events
The spatial distributions of moderate and extreme drought intensity and frequency over SA and its subregions are presented in Figure 3. Here, we considered the drought intensity and frequency over the total number of months, indicating the scPDSI values smaller than −2 (moderate) and/or greater than −4 (extreme). Thus, we only considered the months of scPDSI with the threshold <-2 and >-4 as moderate drought, and that of scPDSI with the threshold <-4 as extreme drought. In Figure 3a,b, R-2, R-3, and R-4 observed moderate drought intensity (−4.5), while reciprocal results over these subregions experienced extreme drought intensity with amounts of −6.0 to −4.5. To assess the temporal and spatial pattern of drought intensity (moderate and extreme) trends at each grid point, we investigated the regional trend in drought intensity for the four subregions of SA over the study period (see Figure 4). In R-1 (Figure 4a), extreme drought intensity detected a significant drying trend for both moderate and extreme drought intensity, suggesting that drought intensity had significantly increased over the recent past. However, the extreme drought intensity appeared to have largely increased since 2018. Similar patterns were seen in R-2 (Figure 4b), indicating an overall drying trend in moderate and extreme drought intensity during the study period, implying more intensified drought severity in a changing climate. The current study results align with the previous literature that detected allied patterns in these regions [57,61].  In terms of drought frequency, moderate drought frequency was observed about five to seven times per year in R-2 and R-3, followed by R-1 and R-4 with relatively less drought frequency. A similar pattern was seen for those subregions, but extreme drought frequency was observed with the amount of 14 to 20, respectively (see Figure 3c,d). The dominant pattern of extreme drought revealed a significant frequency in R-2, R-3, and R-4, which is in line with previous studies [25,59]. In addition, SA experienced a high frequency of severe droughts in R-2 and R-3. It is worth mentioning that extreme drought intensity and frequency in R-2 and R-3 were relatively sharper than in other subregions. Remarkably, the peak of drought frequency was observed in R-2, which is vulnerable to climate change, and a sharp temperature rise would further accelerate the extreme drought events in the region [60].
To assess the temporal and spatial pattern of drought intensity (moderate and extreme) trends at each grid point, we investigated the regional trend in drought intensity for the four subregions of SA over the study period (see Figure 4). In R-1 (Figure 4a), extreme drought intensity detected a significant drying trend for both moderate and extreme drought intensity, suggesting that drought intensity had significantly increased over the recent past. However, the extreme drought intensity appeared to have largely increased since 2018. Similar patterns were seen in R-2 (Figure 4b), indicating an overall drying trend in moderate and extreme drought intensity during the study period, implying more intensified drought severity in a changing climate. The current study results align with the previous literature that detected allied patterns in these regions [57,61]. During R-3 (Figure 4c), drying trends were observed for both the moderate and extreme drought intensity values during the first two decades (from 1981 to 2000). R-4 (Figure 4d) depicted the sharpest drying trend for moderate values, while the highest drying trends appeared for extreme drought since the 2000s. Recently, numerous researchers found pronounced decreasing trends in annual precipitation over the mountainous regions (i.e., Hindu Kush and foothills of the Himalayans) [62,63].

Observed Drought Trends over Space and Time
To understand how drought events represent the spatial patterns of trends across SA and its subregions, we estimated pixel-wise spatial distributions of scPDSI trends using MK and TS estimators over the period 1981-2020 ( Figure 5). A pronounced drying trend was obvious in R-1, R-3, and R-4 at the rate of −4.5, −4.0, and −3.5/year (Figure 5a), while a relatively less wetting trend was observed in R-2 at the rate of 1.5/year. Interestingly, a sharp drying trend appeared with increasing drought intensity over the Hindu Kush and foothills of the Himalayans at a rate of −3.5/year. The results of the MK trends consistently matched the findings of the TS estimator (Figure 5b), indicating almost similar patterns over these regions. In SA and its subregions, therefore, the drying trends were obvious and significant for most areas, which could be due to increased temperature in recent decades over the region.

Observed Drought Trends over Space and Time
To understand how drought events represent the spatial patterns of trends across SA and its subregions, we estimated pixel-wise spatial distributions of scPDSI trends using MK and TS estimators over the period 1981-2020 ( Figure 5). A pronounced drying trend was obvious in R-1, R-3, and R-4 at the rate of −4.5, −4.0, and −3.5/year (Figure 5a), while a relatively less wetting trend was observed in R-2 at the rate of 1.5/year. Interestingly, a sharp drying trend appeared with increasing drought intensity over the Hindu Kush and foothills of the Himalayans at a rate of −3.5/year. The results of the MK trends consistently matched the findings of the TS estimator (Figure 5b), indicating almost similar patterns over these regions. In SA and its subregions, therefore, the drying trends were obvious and significant for most areas, which could be due to increased temperature in recent decades over the region.
The regional average drought time-series in the four subregions of SA further reaffirm the spatial patterns of trends shown above ( Figure 6). The time series revealed abrupt drying patterns with intensified droughts from 2000 to 2020 in R-1 and R-2, whereas R-3 observed a relatively less drying trend, and R-4 even showed a wetting trend over the 40-year period. Similar findings were reported in previous studies over the regions [3,5]. The previous works attributed the regional drying to the warming climate and/or fewer changes in precipitation. Other studies by Ehsan et al. (2020) [16] and   [20] suggested that R-2 (i.e., Pakistan, Afghanistan, and southwestern India), and R-4 (i.e., southern India) could face severe extreme temperatures by the mid-21st century. It is worth mentioning that R-1, R-2, and R-4 observed an abrupt increasing drying trend during the 2000s-2020s, which may have been caused by natural variability of the regional climate, and a rapid drought intensity/frequency shift in climate regime may occur in the future. MK and TS estimators over the period 1981-2020 ( Figure 5). A pronounced drying was obvious in R-1, R-3, and R-4 at the rate of −4.5, −4.0, and −3.5/year (Figure 5a) a relatively less wetting trend was observed in R-2 at the rate of 1.5/year. Interestin sharp drying trend appeared with increasing drought intensity over the Hindu Ku foothills of the Himalayans at a rate of −3.5/year. The results of the MK trends consi matched the findings of the TS estimator (Figure 5b), indicating almost similar pa over these regions. In SA and its subregions, therefore, the drying trends were o and significant for most areas, which could be due to increased temperature in rece ades over the region.  The regional average drought time-series in the four subregions of SA further reaffirm the spatial patterns of trends shown above ( Figure 6). The time series revealed abrupt drying patterns with intensified droughts from 2000 to 2020 in R-1 and R-2, whereas R-3 observed a relatively less drying trend, and R-4 even showed a wetting trend over the 40year period. Similar findings were reported in previous studies over the regions [3,5]. The previous works attributed the regional drying to the warming climate and/or fewer changes in precipitation. Other studies by Ehsan et al. (2020) [16] and Mishra et al. (2020) [20] suggested that R-2 (i.e., Pakistan, Afghanistan, and southwestern India), and R-4 (i.e., southern India) could face severe extreme temperatures by the mid-21st century. It is worth mentioning that R-1, R-2, and R-4 observed an abrupt increasing drying trend during the 2000s-2020s, which may have been caused by natural variability of the regional climate, and a rapid drought intensity/frequency shift in climate regime may occur in the future.

ScPDSI Possible Linkages with Different Climate Factors and Dynamic Drivers
To understand more deeply, we further assessed the temporal variability of VPD, NDVI, PET, precipitation, and temperature as dominant climatic factors influencing drought events over the diverse subregions of SA (see Figure 7). Figure 7 depicts the interannual variation in VPD, NDVI, PET, precipitation, and temperature anomalies in four subregions of SA from 1981 to 2020, based on the regional averages at each grid. During VPD (Figure 7a-d), the interannual changes in VPD were more distinct over R-1 and R-2,

scPDSI Possible Linkages with Different Climate Factors and Dynamic Drivers
To understand more deeply, we further assessed the temporal variability of VPD, NDVI, PET, precipitation, and temperature as dominant climatic factors influencing drought events over the diverse subregions of SA (see Figure 7). Figure 7 depicts the interannual variation in VPD, NDVI, PET, precipitation, and temperature anomalies in four subregions of SA from 1981 to 2020, based on the regional averages at each grid. During VPD (Figure 7a-d), the interannual changes in VPD were more distinct over R-1 and R-2, followed by R-3 and R-4, respectively. For instance, the interannual VPD trend had declined over R-3 and R-4 in the initial epoch (i. e., during 1990-2000). The considerable decline in interannual VPD over these regions could be partly attributed to higher surface temperature and/or lesser changes in precipitation. However, highly fluctuating temporal variability was seen from region to region. It is worth stating that the interannual VPD pattern was leading to vapor pressure variability due to the regional precipitation dependency on soil moisture for precipitation. Moreover, we also found that the frequency of extreme surface VPD had decreased substantially during earlier epochs in R-4, whereas an increased frequency was seen during later epochs. Overall, the NDVI exhibited a declining trend in all of the climatic regions, especially in earlier decades. However, an extreme declining trend appeared in recent decades, suggesting that vegetation stress continued even after drought ended (Figure 7e-h). This finding is consistent with drying conditions observed in previous studies [64,65]; they found that in drought conditions in R-1, R-3, and R-4, water scarcity was one of the main drivers of variations in vegetation stress. It should be noted that these regions had witnessed frequent and intense climate extremes in the recent past, which can directly be accredited to the increasing trend of aridity and temperature [3,25].
In addition, PET, temperature, and precipitation are the major climatic factors of dryness/droughts that regulate the climate system from regional to global [25,64]. In Figure 7i-l, the regional precipitation pattern had exhibited vigorous changes with a consistent declining trend in the early period (1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000), while an inclining trend was observed in the latter period (2001-2020) of the study period. A similar pattern appeared for surface air temperature (Figure 7m-p), depicting sharp and accelerated regional temperature during the 1990s and onwards, suggesting the aggression of climate change in SA and its sub-climatic regions. Interestingly, in the latter decades (2001-2020), a significant increasing trend of surface temperature was revealed in R-1, R-3, and R-4 at the rate of 0.6 • C, 0.8 • C, and 0.9 • C/year, respectively. Results of this study further support the findings of previous studies [52,[66][67][68], which described a significant decline in precipitation (10 to 15 mm/year), while the incline in air temperature (0.6 to 0.9 • C) detected drought/dry conditions in utmost parts of the study region. The temporal trend of PET observed an increasing pattern over many regions of SA except R-1, with higher intensity during 1981-2000 and severe extreme intensity during the period 2000-2020 (Figure 7q-t). However, the linear trend revealed that SA regional precipitation (PET) had enhanced at the rate of 5 (0.5) mm/year (R-1), 6 (0.4) mm/year (R-2), 8 (0.7) mm/year (R-3), and 10 (0.8) mm/year (R-4), respectively. It is worth mentioning that the simultaneous rise in surface air temperature and PET emphasized that surface air temperature had played a more escalating role than precipitation in the development of PET in the recent past. Though decadal change varied from region to region, the SA precipitation abruptly declined, which could possibly have an important role in intensifying the development of PET. To summarize, it can be concluded that both precipitation and temperature had largely influenced PET in diverse regions of the target region, which could have ultimately affected the inclination of SA droughts in the study regions. Remote Sens. 2022, 14, 3219 14 of 21  We further explored the linkages between droughts (scPDSI) and hydroclimatic and vegetation variables; we examined the spatial variations in climate forcing and determined its relationship with regional droughts. We assessed the nexus between drought and temperature, precipitation, PET, VPD, and NDVI for the period 1981-2020 ( Figure 8) to better understand their controlling roles. The results of the study indicated that the scPDSI spatial correlation with precipitation was expected, as precipitation participated in scPDSI computation, and the contribution of precipitation is accounted for in the model (Figure 8a). Nevertheless, using different aspects of precipitation to explore drought's relationships serves as a further verification step for the drought. With this aim, we computed the sensitivity of drought to precipitation changes for each grid point. However, an inverse correlation pattern appeared for temperature, and PET exhibited a negative correlation (−0.6 to −0.8) over those regions, where a positive correlation (0.6) appeared in R-4 for both temperature and PET (Figure 8b,c). The regions that revealed a strong negative association between drought and temperature and PET, indicated that the extent of correlation fluctuated from region to region due to a significant increase in temperature, particularly over SA. Interestingly, many regions of SA had revealed a strong positive correlation (0.6 to 1) between precipitation and drought index, including R-1, R-2, and R-3 (Figure 8a). Similarly, in NDVI (Figure 8d), correlation patterns changed as expected with the same significance and order of magnitude as in precipitation and PET. Our results support the findings of Hussain et al. (2017) [69], who reported that precipitation favored vegetation growth, and drought events may disturb vegetation growth over the target regions. For VPD (Figure 8e), the study also found a decline in VPD over R-2, R-3, and R-4 during the same period; however, a substantial positive association between scPDSI and VPD was observed in R-1 at the rate of 0.8, respectively. vegetation variables; we examined the spatial variations in climate forcing and determined its relationship with regional droughts. We assessed the nexus between drought and temperature, precipitation, PET, VPD, and NDVI for the period 1981-2020 ( Figure 8) to better understand their controlling roles. The results of the study indicated that the scPDSI spatial correlation with precipitation was expected, as precipitation participated in scPDSI computation, and the contribution of precipitation is accounted for in the model (Figure 8a). Nevertheless, using different aspects of precipitation to explore drought's relationships serves as a further verification step for the drought. With this aim, we computed the sensitivity of drought to precipitation changes for each grid point. However, an inverse correlation pattern appeared for temperature, and PET exhibited a negative correlation (−0.6 to −0.8) over those regions, where a positive correlation (0.6) appeared in R-4 for both temperature and PET (Figure 8b,c). The regions that revealed a strong negative association between drought and temperature and PET, indicated that the extent of correlation fluctuated from region to region due to a significant increase in temperature, particularly over SA. Interestingly, many regions of SA had revealed a strong positive correlation (0.6 to 1) between precipitation and drought index, including R-1, R-2, and R-3 (Figure 8a). Similarly, in NDVI (Figure 8d), correlation patterns changed as expected with the same significance and order of magnitude as in precipitation and PET. Our results support the findings of Hussain et al. (2017) [69], who reported that precipitation favored vegetation growth, and drought events may disturb vegetation growth over the target regions. For VPD (Figure 8e), the study also found a decline in VPD over R-2, R-3, and R-4 during the same period; however, a substantial positive association between scPDSI and VPD was observed in R-1 at the rate of 0.8, respectively.
The results showed that the depletion in precipitation and soil moisture led to an escalation in VPD, decreasing vegetation growth. Moreover, decoupling of precipitation and surface VPD could separately intensify the energy [65,70]. Atmospheric VPD and surface air temperature were negatively associated in the majority of subregions of SA, mostly determined by land-atmospheric trajectories [70]. Thus, in regions where surface VPD and air temperature were negatively associated (e.g., R-2, R-3, and R-4), the landatmospheric trajectories influenced dry and hot extremes [71,72]. For example, low precipitation and high temperature led to concurrent extreme events [20].  The results showed that the depletion in precipitation and soil moisture led to an escalation in VPD, decreasing vegetation growth. Moreover, decoupling of precipitation and surface VPD could separately intensify the energy [65,70]. Atmospheric VPD and surface air temperature were negatively associated in the majority of subregions of SA, mostly determined by land-atmospheric trajectories [70]. Thus, in regions where surface VPD and air temperature were negatively associated (e.g., R-2, R-3, and R-4), the land-atmospheric trajectories influenced dry and hot extremes [71,72]. For example, low precipitation and high temperature led to concurrent extreme events [20].
Next, we studied changes in the effect of ENSO and IOD on the scPDSI variability in the study regions (see Figure 9a-c). The extent and spatial distribution of relationships between the multi-variate ENSO index and the NDVI were largely similar in R-2, R-3, and R-4, with a few exceptions. Regions with a substantial correlation between ENSO and NDVI showed considerable decrease for R-2, R-3, and R-4 (from −0.6 to −1) and increase for R-1 (from 0.6 to 0.8), respectively. The ENSO spatial correlations showed stronger dry conditions over R-2, R-3, and R-4 during NDVI and scPDSI. However, wet conditions also exhibited strengthening over R-1 during ENSO (Figure 9a-c). In contrast, the VPD revealed a relatively higher spatial correlation with the ENSO index in most of the target regions (Figure 9c). Results indicated that ENSO patterns affected a similar area within most regions with enhanced magnitudes in all climate indicators, which inferred that the variability of ENSO interconnection remained largely static with warming [62,73]. Given the dominant role of multivariate ENSO events in determining drought events in a warmer climate, these findings emphasized the significance of the current ENSO-drought association and the related mechanisms [69,74,75]. For instance, the dynamic variabilities in SASM precipitation usually resulted in extreme wet and dry episodes, leading to flooding and drought events, respectively, with devastating socio-economic impacts [10]. The magnitude and onset timing of monsoon precipitation are uncertain, which often results in extreme events such as droughts and floods in SA regions [69,74]. Next, we studied changes in the effect of ENSO and IOD on the scPDSI variability in the study regions (see Figure 9a-c). The extent and spatial distribution of relationships between the multi-variate ENSO index and the NDVI were largely similar in R-2, R-3, and R-4, with a few exceptions. Regions with a substantial correlation between ENSO and NDVI showed considerable decrease for R-2, R-3, and R-4 (from −0.6 to −1) and increase for R-1 (from 0.6 to 0.8), respectively. The ENSO spatial correlations showed stronger dry conditions over R-2, R-3, and R-4 during NDVI and scPDSI. However, wet conditions also exhibited strengthening over R-1 during ENSO (Figure 9a-c). In contrast, the VPD revealed a relatively higher spatial correlation with the ENSO index in most of the target regions (Figure 9c). Results indicated that ENSO patterns affected a similar area within most regions with enhanced magnitudes in all climate indicators, which inferred that the variability of ENSO interconnection remained largely static with warming [62,73]. Given the dominant role of multivariate ENSO events in determining drought events in a warmer climate, these findings emphasized the significance of the current ENSO-drought association and the related mechanisms [69,74,75]. For instance, the dynamic variabilities in SASM precipitation usually resulted in extreme wet and dry episodes, leading to flooding and drought events, respectively, with devastating socio-economic impacts [10]. The magnitude and onset timing of monsoon precipitation are uncertain, which often results in extreme events such as droughts and floods in SA regions [69,74].  During IOD (Figure 9d-f), our results revealed an extensive and consistent negative correlation of the IOD anomalies (e.g., R-2, R-3, and R-4) and a positive correlation of the IOD anomalies (e.g., R-1) with NDVI, scPDSI, and VPD, respectively, implying that IOD patterns were contributing to droughts in these regions. On the contrary, we also found that the IOD and drought index had a diverse impact across these regions. For example, the IOD positively correlated with the drought index in R-1, but a negative correlation was observed in R-2, R-3, and R-4. This indicated that the positive phase of IOD conditions promoted droughts over R-1. For example, El Niño conditions may co-occur with the positive phase of IOD as they tend to determine warmer SSTs over the southeastern Indian Ocean and cooler SSTs over the southern Indian ocean through the atmospheric bridge and oceanic Indonesia through-flow [75][76][77][78][79][80][81]. The study did not negate the potential impact of large-scale oceanic indices but rather postulated the hypothesis that these controls need to be reviewed with clear links toward drought variability over SA and its subregions. Another possible reason could be the quality of gridded datasets, since it is recorded for synoptic-scale use, and climatological processes usually occur on longer-time scales, thus it may lack such large-scale signals. The present study has attempted only long-term statistical composites and is subjected to limitations. Hence, a comprehensive numerical experimental validation is required to further investigate the multifaceted connection of these parameters at play. Furthermore, though the oceanic variability may have little influence, the effect of Tibetan Plateau thermodynamics and Himalayan thermal controls can be deeply examined. Thus, a detailed investigation is required to understand how the positive phase of IOD and oceanic modes interrelate with El Niño to affect drought features over diverse regions such as SA and, as a result of compound heatwaves and droughts.

Conclusions
Droughts are linked mainly to ecological, environmental, and social factors. Given the cumulative regional and global interconnections of socioeconomic systems, understanding the recent observed features of droughts and their potential dynamic drivers in a warmer climate is key for a wide range of interrelated climate sensitive factors. The present study examined the spatiotemporal variability of long-term droughts and their characteristics over SA and its diverse climatic regions during 1981-2020. To explore different climate influencing factors and vegetation changes, we assessed the drought variability in response to temperature, precipitation, PET, NDVI, and VPD, at each grid point over the same period. To better understand the effect of teleconnections on drought evolution, this study used ENSO and IOD as potential dynamic drivers during 1981-2020.
Results indicated that drought severity and frequency fluctuated from region to region. However, significant severe droughts were observed in R-4, followed by R-2 and R-3, with more frequent droughts during 1981-2020. In terms of temporal trend, we also found significant prolonged drying trends in R-2, R-3, and R-4, which further clarified the spatial drying patterns over those regions; however, a relatively wetting trend was observed in R-1, for both, over space and time. Generally, our results showed a significant effect of surface air temperature, PET, NDVI, and VPD on drought events in several regions of SA. Precipitation had substantially declined, while NDVI and PET had increased extensively over R-2, R-3, and R-4, from 1981 to 2020. Additionally, VPD had increased in many subregions of SA due to a substantial upsurge in surface air temperature, particularly in R-2 and R-4. In addition, the correlations between drought index and climate influencing factors were relatively stronger than 0.8 for nearly 80% of the grid points, which included most of the regions where the variations in the means of monthly accumulated PET and VPD values depicted large changes. Generally, a strong positive correlation was found between drought index, precipitation, and NDVI across R-1, R-2, and R-4, whereas temperature, PET, and VPD exhibited a negative correlation over these regions except R-4, respectively. However, no obvious connection was detected with ENSO events, or IOD and drought onset timing, as explored for certain regions of SA, especially R-3 and R-4. The findings deduced the probability that the precipitation variability over these regions had an insignificant relationship with ENSO, IOD, and drought evolution. The trend in temperature and precipitation could be due to extreme climatic events being the driver, namely, their particular magnitude, frequency, and intensity during each month. Thus, future work can consider prospective lead-lag interactions between multivariate ENSO and precipitation in space and time compound drought and heatwave analyses.