Assessing Spatiotemporal Drought Dynamics and Its Related Environmental Issues in the Mekong River Delta

Drought is a major natural disaster that creates a negative impact on socio-economic development and environment. Drought indices are typically applied to characterize drought events in a meaningful way. This study aims at examining variations in agricultural drought severity based on the relationship between standardized ratio of actual and potential evapotranspiration (ET and PET), enhanced vegetation index (EVI), and land surface temperature (LST) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) platform. A new drought index, called the enhanced drought severity index (EDSI), was developed by applying spatiotemporal regression methods and time-series biophysical data derived from remote sensing. In addition, time-series trend analysis in the 2001–2018 period, along with the Mann–Kendal (MK) significance test and the Theil Sen (TS) slope, were used to examine the spatiotemporal dynamics of environmental parameters (i.e., LST, EVI, ET, and PET), and geographically weighted regression (GWR) was subsequently applied in order to analyze the local correlations among them. Results showed that a significant correlation was discovered among LST, EVI, ET, and PET, as well as their standardized ratios (|r| > 0.8, p < 0.01). Additionally, a high performance of the new developed drought index, showing a strong correlation between EDSI and meteorological drought indices (i.e., standardized precipitation index (SPI) or the reconnaissance drought index (RDI)), measured at meteorological stations, giving r > 0.7 and a statistical significance p < 0.01. Besides, it was found that the temporal tendency of this phenomenon was the increase in intensity of drought, and that coastal areas in the study area were more vulnerable to this phenomenon. This study demonstrates the effectiveness of EDSI and the potential application of integrating spatial regression and time-series data for assessing regional drought conditions.


Introduction
Drought originates from precipitation shortage in a region over an extended period of time and it can be considered as an abnormal climatic event aligned with water resource deficits [1,2].Drought is a serious natural disaster that has adverse effects on the environment and human activities [3][4][5].Depending on the causes, droughts are often classified into four categories, including meteorological, agricultural, hydrological, and socioeconomic drought [6,7].Meteorological drought is mainly affected by the reduction of precipitation and high temperature, while agricultural drought is caused by soil moisture deficit [8,9].The hydrological drought occurs as a consequence of stream flow shortage and groundwater reduction and socioeconomic drought takes place when water demand for socioeconomic activities cannot be fulfilled by water supply [9].Historically, most drought studies were conducted using data from hydro-meteorological stations [10][11][12].The amount and quality of the meteorological data, however, vary between individual meteorological stations, creating challenges in monitoring and evaluating the drought events in the case of large-scale studies [13][14][15].It is widely accepted that drought is a regional phenomenon and occurs in virtually all climatic regimes [16].Hence, it is necessary to conduct a space and time analysis of this issue.
Drought indices are efficient measurements for examining the drought categories in a particular area [17].There is a difficulty in determining agricultural drought due to the direct consequences of meteorological drought on agriculture [7,9,18].Therefore, to understand and examine better agricultural drought variation, the computation of a drought monitoring index is performed prior to the first step [9].Numerous drought indices have been proposed to estimate drought variations [19][20][21].Progress in drought observation can be categorized into climatological and ecosystem groups [2,14].The climatological perspective focuses on several parameters, including rainfall, air temperature, soil moisture, potential evapotranspiration, and snow, while characteristics of plant and land cover from remotely sensed data are the main factors from an ecosystem point of view.Several studies have proved that an efficient drought measurement strategy should be based on composite and multi-index drought models since the range of vegetation responses to drought (and its variation) should be analyzed in order to take ecosystem diversity into account [22][23][24].The multi-index drought perspective, therefore, plays an important role in monitoring agricultural drought efforts in this study.
The need to explore more appropriate and efficient agricultural drought indices is increasing and currently being studied by scientists around the world [2].Compared to the commonly used drought indices, the Palmer drought severity index (PDSI) was the first effective index applied to explain agricultural drought [24,25].The parameters to compute the PDSI consist of precipitation, temperature, and soil moisture.However, the index is difficult to formulate a universal drought index due to the complexity of drought over space and time, as documented by Keyantash and Dracup [26] and Mu et al. [24].Additionally, the standardized precipitation index (SPI), which was developed by McKee et al. [19], has been used to examine meteorological drought.As it uses only precipitation for drought assessment, the calculation process for this index is simple and flexible based on the time scale, and hence, has been applied internationally [25,27].The index, however, is limited when using it for agricultural areas [8], because it does not consider other critical variables such as evapotranspiration and soil moisture, which can significantly influence drought conditions [28].
Besides, wide ranges of drought indices have been developed using remote sensing data [16,[29][30][31] for large-scale ecosystem drought monitoring due to its advantage for surface soil moisture detection [32,33].Among these, many indices are based on a single parameter [14], such as the normalized difference vegetation index (NDVI) [34,35], the enhanced vegetation index (EVI) [36], the vegetation condition index (VCI) [37], the standardized vegetation index (SVI) [38,39], the temperature condition index (TCI) [40], and the normalized difference temperature index [41,42].More recently, drought indices that use a combination of land surface temperature (LST) and NDVI [29,33], or NDVI and the ratio of evapotranspiration (ET) to potential evapotranspiration (PET) [20,24,25,43] have been developed.For instant, the drought severity index (DSI)-a worldwide applicable drought index-was developed and first published in 2013 by Mu et al. [24,25], integrating the MOD16 ET/PET and MOD13 NDVI products, based on a series of composite and multi-index models.It was efficiently used for drought monitoring and mitigation efforts in all vegetated areas from 2000 to 2011 [25].However, NDVI has lower sensitivity than EVI, which limits it when being applied to the detection of the vegetation greenness in high biomass regions [44,45].Several subsequent studies modified the DSI by replacing factors around the drought algorithm for EVI2 [30] or climate data derived from the advanced very high resolution radiometer (AVHRR) [25] and applied them for their typical regions.Hence, the analysis of the input factors that have a role in drought modeling is vital for the assessment of the drought status at the area of analysis.Moreover, drought variations are varied between geographic regions depending on the climatic conditions and characteristic of land cover and land use [25,27].Hence, this study aims at developing an efficient drought index using multiple parameters at a regional scale in order to estimate its drought status.Mekong River Delta (MRD) in Vietnam is selected to implement the study as it is a pivotal agricultural region of Vietnam and drought is a significant problem in the dry season along with salinity.The drought conditions in the MRD have a coherent pattern with total annual water supply of the Mekong River over the last 45 years [46], causing water shortages for agricultural land, especially severe during the dry season (from November to April).Extreme drought events were recorded in the study area in 1992, 1993, 1998, and 1999 [47].The recently extreme drought occurrences had taken place in 2003, 2006, and 2017 [48].
There are various types of remotely sensed data available to examine drought severity across various regions in the world, such as the advanced very high resolution radiometer (AVHRR) [44,49,50], the moderate resolution imaging spectroradiometer (MODIS) [24,33], Landsat [51,52], or Sentinel [53].Among these sensors, MODIS has great temporal resolution of 1-2 days with high radiometric resolution image (12 bit) and collects data for 36 spectral bands [54,55].Whereas, higher spatial resolution data such as Landsat (30 m) is challenged in providing an appropriate number of high-quality images [51,56] and Sentinel data is not sufficient to investigate long-term analysis as it has only been available since 2015 [57].In fact, several studies revealed the effectiveness of MODIS data for agricultural drought assessment in tropical delta, particularly in the Mekong River Delta (MRD) as adapting to challenges of persistent cloud cover due to its advantages of temporal resolution [33,58,59].Therefore, with a spatial resolution ranging from 0.25 to 1 km, MODIS data is appropriate for examining drought issues in the case study of Mekong River Delta [33,44,60,61].In this study, multi-temporal MODIS imagery was applied for analyzing spatiotemporal drought variations in the MRD.
Different types of spatiotemporal regression methods have been utilized to examine the spatiotemporal variations of environmental factors derived from remotely sensed data [62][63][64][65][66].The Mann-Kendal significance test and the Theil Sen slope estimator are powerful tools for determining the trend of variations in environmental factors because of their innovations, as compared to ordinary least-squares (OLS) regression [66,67].In addition, the Pearson correlation coefficient has been widely applied for investigating the relationship among variables in OLS regression [66].However, relationships among physical variables are applied at different spatial scales and in separate geographical locations.Hence, and in order to gain some insight into this process, local differences are thus needed to be considered [68].Brunsdon et al. [69] developed a local form of the standard linear regression method, called geographically weighted regression (GWR).It has been applied efficiently for evaluating the spatial non-stationarity and local correlation coefficient at the pixel level [70][71][72].Our study, therefore, proposed and applied the enhanced drought severity index (EDSI) algorithm-developed from MODIS satellite-derived LST, ET, PET, and EVI factors-and applied the Mann-Kendal and Theil Sen (MK-TS) as well as GWR, in order to conduct a spatiotemporal analysis of drought situations.
Overall, our research aimed at the (i) analysis of the spatiotemporal variations of a series of environmental parameters of interest (i.e.ET, PET, EVI, and LST), (ii) assessment of the local correlation among each other (considering also their standardized values), (iii) development of a new drought index based on the concept of the drought severity index (DSI) [24], and (iv) pattern dynamics and trend analysis in agricultural drought.The Mekong River Delta was selected as a case study because of the impact of drought along with salinity during the dry season on the region [57,73].Besides, the comparison of EDSI and other drought indices was performed to test our hypothesis and to enhance its capabilities.

Study Area
The Mekong River Delta (MRD) (located between 8 • 34 N-10 • 24 N and 104 • 27'E-106 • 48 E), which covers an area of approximately 40,000 km 2 , is a major agricultural region in the Lower Mekong River Delta (LMRD) (Figure 1).This delta attains several advantages, such as having a large area of fluvial soil, a criss-crossing river system, and the regularization of river runoff by Tonlé Sap Lake in Cambodia.The agricultural land accounts for 72.68% (approximately 31,283.14km 2 ) and more than 90% of this land use type is used for rice production [33,58].Open water, urban, and other land use types occupy 5.17%, 10.27%, and 11.88% of the total area, respectively.Besides, the major productive crops accounted for more than 50% of the rice-cultivated area, 71% of the aquaculture area, 30% of the value of agricultural production, and 54% of fisheries in Vietnam in 2015 [74].This delta is therefore, generally considered crucial for food security and the economic development of Vietnam.Most of the locations in the delta have similar climatic conditions (i.e., air temperature, precipitation, and evapotranspiration) as the major landform in the delta is flat plain.Besides, the air temperature remains stable through the whole year with a mean value of over 26 • C [12,75].This is an advantage condition for agricultural production.Regarding the landscape in the study area, it can be divided into upper (inner) delta and lower (outer) delta plain [76], and the land elevation between these parts is not significant.However, with elevations' difference in the range from 0.5 to 1.5 m above the sea level [77], the delta is identified as a low-lying area, which are extremely vulnerable to climate change impacts [58,78].In addition, located in a tropical monsoon and sub-equator region, climate in the MRD has two seasons.The rainy season lasts from May to October and the dry season begins in November and ends in April, with less than 20% of the total annual precipitation.Consequently, drought in MRD occurs significantly in the dry season and the most damaging effect of drought takes place in the last week of March and the first week of April.Therefore, it is recommended that examining drought events in the study site should focus on the dry season.The spatial variations of drought may be influenced by the character of the land surface rather than climatic factors.
In recent years, however, droughts, along with salinization in the dry season, have created a serious adverse impact on agriculture and water supply as a result of global warming and upstream hydropower dams [79][80][81].The study area, for instance, was affected by the most intense drought and salinity on record in 2016-2017, leading to a serious fresh water shortage and a decrease in cultivated area, especially severe in paddy land which is a form of rice-farming [10].Hence, it is necessary to examine the variations in drought intensity over time as a predictor of agricultural development.

Data
Earth observation data from the different satellite sensors used in this research is shown in Table 1.A range of MODIS products, including LST (MOD11A2), ET, PET (MOD16A2), and EVI (MOD13A1), obtained from the Land Processes Distributed Active Archive Center (LP DAAC), were the main sources when computing the EDSI, based on the spatiotemporal relationship between these factors.With a spatial resolution of 0.5-1 km, using MODIS products is appropriate as it enables us to map drought severity at a map scale of 1:2,000,000-1:1,000,000 [82].Whereas, there are only eleven weather stations in the study area, so using meteorological data is difficult for mapping spatial drought pattern.The MODIS products (i.e., MOD16A2, MOD13A1, and MOD11A2) were used to obtain ET, PET, EVI, and LST and their spatiotemporal trends, respectively.According to Running et al. [83], the MOD16 ET/PET algorithm used the well-known Penman-Monteith (PM) equation to calculate global remotely sensed ET [20,25], and integrate both PM and Priestley-Taylor methods to estimate PET [20,83] with a combination of remote sensing [84] and global meteorological data.The reflections of MOD16 ET/PET (i.e., evaporation from canopy interception, wet and moist soil surfaces, and transpiration from canopy stomata) and its relations to both surface energy partitioning and environmental constraints on ET are also described in several studies [20,24,25,83].Besides, the MOD13A1 provides a piece of information about vegetation indices (i.e., normalized difference vegetation index and enhanced vegetation index) and its potential associated with climate changes (e.g., increasing frequency increase and drought severity).The vegetation responses can be observed through vegetation greenness changes with land-atmosphere water, carbon and energy fluxes, and linked climate feedbacks [24,30,85].Additionally, the LST derived from MOD11A2 has been used as a key factor to examine drought severity [33,40,86] as a response of soil moisture and its texture [29].Based on the relationship among parameters and the subsequent computation of the drought severity index (DSI) [24], an integration of the operational MOD16 ET/PET, MOD13 EVI, and MOD11 LST products was employed to develop a remotely sensed agricultural drought index for the MRD as an enhancing piece of information of water stress.
A transformations are performed such that derived standardized factors of each environmental variable follow a normal distribution and reduce the deviations of long-term mean values because 8 day environmental variables are not normally distributed [4,87].In addition to satellite data, the meteorological data obtained from eleven meteorological stations on the ground (Figure 1) were used to calculate the reconnaissance drought index (RDI) [17] and the standardized precipitation index (SPI) [19].The RDI and SPI were then applied to examine the performance of remotely sensed drought indices.Then, either the locally spatial relationship among environmental parameters or its standardized values were explored to discover the significant probability of input data and drought model.Finally, a time-series of the EDSI was conducted to assess drought variations between space and time.

Enhanced Drought Severity Index
Drought variations, in fact, depend on various factors such as precipitation, soil moisture, and ET.Several studies [29,88,89] have reported that there is a significant relationship between drought and other relevant parameters (e.g., vegetation indices, LST, rainfall, ET).The ET reflects the status of ecosystem function and is directly related to water, carbon, and energy cycles of the land surface [45].Besides, LST is a major factor controlling the physical processes that are responsible for the land surface balance of water, energy, and CO 2 [33,40].Additionally, LST is often used to indicate drought intensity because the surface water loss by evapotranspiration increased under high thermal conditions [31,56].Besides, the behavior of LST in response to soil moisture provides information about its texture [90,91].Several studies have used LST as the main source of satellite-derived information in order to develop different types of drought indices [5,40,42].We therefore consider that these factors should be integrated into drought modeling, in order to increase their accuracy, because they enable us to provide more pieces of appropriate information on surface moisture and water deficits.
From the correlation among environmental factors, we have developed an integrated index, the enhanced drought severity index (EDSI), using standardized ratios of Z RT , Z LST , and Z EVI , as given in Equations ( 1) and (2).
The EDSI anomaly is dimensionless and might vary in (−∞, +∞) for rainy to dry period variations, respectively.The ratio (RT) between ET and PET in any given month, was developed by Yao et al. [92] and employs reanalysis of meteorological data and remotely sensed data, and the MODIS data MOD16 ET/PET products were used as a factor of the terrestrial water availability and the atmosphere [30,45].This indicator manifests a condition of wetness or dryness [24,25], as shown within Equation (3).Besides, both the vegetation and LST reflect the water-energy balance of the surface [29,33], and thus, are good indicators for drought assessment.Consequently, the additional environment variables (i.e., ET, PET, EVI, and LST) provide information about water content as an effective indicator of agricultural drought monitoring.Moreover, it is necessary to apply remotely sensed data for analyzing drought in a large area like the Mekong River Delta because only 11 meteorological stations have been used to monitor weather conditions of around a 40,000 km 2 area.
In addition, in this study, the performance of EDSI was subsequently assessed by performing regression analysis between EDSI and DSI developed from remotely sense data [24] versus RDI [27,93] and SPI [16,19,27] measured at the meteorological stations.
The drought severity index (DSI), which combines NDVI and the ET/PET ratio [24,25], is derived by the operational MOD 16 ET/PET and MOD13 NDVI products.The DSI is estimated using the standardized ratios of RT (Z RT ) and NDVI (Z NDVI ), as shown in Equations ( 4) and (5).
In Equations ( 4) and ( 5), σ Z and Z are the standard deviation (SD) and mean respectively, of Z.
Both of the meteorological drought indices (i.e., RDI and SPI) are computed by the following Equations, from ( 6) to (7): The RDI was developed to detect the water deficit in a more accurate way, representing a balance between input and output parameters in a water system [27,93].It is based on both cumulative rainfall (R) and potential evapotranspiration (PET).In Equation ( 6), the initial value α (i) k of RDI is calculated for the i-th year in a time basis of k (months), PET ij is the precipitation and potential evapotranspiration of the j-th month from January to December of the i-th year, R i is annual rainfall of i-th year respectively, and N is the available interval of years.Besides, the SPI was proposed by McKee et al. [19,27].It is based solely on precipitation measurements and accounted for how much rainfall for a given period of time has deviated from historically established norms [16], shown in Equation (7).R and δ are the average and standard deviation of rainfall in the whole period.This study used RDI and SPI to compare the performance of our newly developed drought index (EDSI) and the previous drought index (DSI).

Spatiotemporal Regression Analysis
A pixel-based time-series trend analysis of drought parameters was performed to examine the environmental factors and drought variations in the study area.In this model, the year sequence of input factors from 2001 to 2018 is treated as an independent variable and each image stack values of EVI, LST, ET, and PET and the drought index images (for the same period) are the corresponding dependent variables.
The space-time model can be used to examine the trend of environmental changes in the sense that the analysis generates slope coefficient values that represent the mean annual change of each environmental parameter [57].In addition, the model includes the Mann-Kendal (MK) significance test, a tool to evaluate the statistical significance of the spatiotemporal trend.The MK test [94] is an effective method for examining trends in time series analysis.Unlike the linear regression model, the MK test enables the detection of monotonic (not only linear) trends without any assumption about the data distribution [63,95].The MK test was first used to analyze trends in hydro-meteorological data [96,97].It statistically evaluates if the change in the values of a variable over a given period followed an upward or a downward trend.In our study, only the pixels which have statistically significant changes (p < 0.01) are further considered for analysis.This provided us with results to determine if each of these environmental measures has significantly changed (positively or negatively) at each pixel in the study area over the examined 18 year period.
Afterward, local correlation among ET, PET, EVI, and LST, as well as these values and drought indices, were conducted using geographically weighted regression (GWR).The GWR, developed by Brunsdon [69], is a spatial multivariate regression and has been widely employed to characterize spatially varying relationships [98].It is a nonparametric model of spatial drift that relies on a sequence of locally linear regressions to produce estimates for every point in space by using a subset of information from nearby observations [99].It is a relatively simple technique that extends the traditional regression framework of equation: where, y i is the i-th observation of the dependent variable, x ik is the i-th observation of the k-th independent variable, the ε i are independent normally distributed error terms with zero means, and each a k must be determined from a sample of n observations.The method can be used to spatialize discrete point data on the assumption that auxiliary, independent variables are known and continuous in space, or, technically, they can be provided as raster layers.The main results from this model are obtained considering a local r-value and slope coefficients within a p-value less than 0.01, which was then used to obtain the relationship among environmental parameters and drought indices.

Assessment of the Enhanced Drought Severity Index
The results of the regression analysis for the correlation among drought indices are presented in Table 2. Strong positive relationships between measured indices at ground observation and spectral indices were found.All regression models showed high r values that were above 0.6.However, EDSI illustrated a higher performance compared to DSI as it achieved better r and p values.We therefore analyzed the correlation among EDSI, SPI, and RDI to adjust the drought categories for EDSI (Table 3).Besides, with a majority of local r values around 0.9, strongly positive spatial correlation values across the Mekong River Delta between the EDSI and DSI, computed by remotely sensed data, are illustrated in Figure 2 using the geographically weighted regression.Areas with yellow to orange represent a certain level of differences between EDSI and DSI.In general, the EDSI classification scale is somewhat similar to DSI up to a certain extent.It also revealed that the performance of a single drought index may depend on specific circumstances.Although many drought indices have been developed in the world, none of them is comprehensive enough to be applicable to all regions of an enlarged nation [100], especially like Vietnam due to the differences in climatic conditions, topography, and local characteristic of land cover and land use.Vietnam is extended along the longitude (8 • -23.5 • N) leading to air temperature and precipitation patterns which are various from North to South.Most of the developed agricultural drought indices, which approached a single parameter or used the long-term deference between the rainfall and evapotranspiration, cannot fundamentally represent the water shortage processes [101,102].Besides, a limited number of meteorological stations in the MRD also creates a challenge to monitor drought at the pixel level.Although the climatic conditions and vegetation are quite similar over this delta, criteria such as land use and land cover (LULC) types and pattern play a crucial role in drought variations [57,73,80].Hence, it is recommended to define appropriate drought index for the region based on the multi-parameter perspective and through the application of remote sensing data and techniques.

Relationships among the Enhanced Vegetation Index (EVI), Land Surface Temperature (LST), Evapotranspiration (ET), and Potential Evapotranspiration (PET)
A general understanding of the correlation among environmental indicators is illustrated in Table 4. Overall, there was a high correlation among them, showing by the R 2 values from 0.64 to 0.93.However, it is observed that the correlation differs across space, as shown in Figure 3. Where, EVI, ET, PET, and LST are the enhanced vegetation index, the actual evapotranspiration, the potential evapotranspiration; and the land surface temperature, respectively.These maps show both positive and negative spatial correlation between EVI, LST, and evapotranspiration.The negative correlation between LST and evapotranspiration is unusual compared to conclusions obtained in previous studies [30,33,40].Several authors explored that PET or ET (which are measured from meteorological stations) are generally positively correlated to temperature (including water temperature and air temperature) [103][104][105][106][107].However, some observations suggest that the correlation between LST and evapotranspiration may be different in particular circumstances.Recently, a number of studies revealed that under limited solar energy, water conditions due to high moisture constraint and low, non-uniform vegetation cover in a dry season, the relationship between evapotranspiration and LST may show a negative correlation [66,108,109].We therefore consider, that the negative correlation performed in our study makes sense as it was implemented during the dry season, a similar context to those mentioned in the above referenced studies.In addition, the influence of LST on air temperature, and the strong relationship between them, has been shown in several studies [110][111][112][113][114]. Consequently, the changes in LST and PET in the study area showed a dependent behaviour in relation to EVI.
Moreover, it is obvious that there was a significant correlation among secondary factors when compared to the original parameters (i.e., ET/PET, LST, and EVI) and a difference in local coefficients (presented in Figure 4).This contributed to concluding that adding LST as an input of the arithmetic is of great importance to conduct in the study area.

The Spatio-temporal Variations of Environmental Factors in the Study Area
The spatio-temporal pattern of variation for each factor in the study area is seen in Figure 5. Overall, there was a significant change in ET, PET, EVI, and LST from 2001 to 2018, with an upward trend for coastal provinces and a downward trend for inland provinces.From 2001 to 2018, ET increased across most of the study area, the total land area with ET increase (6725.71km 2 ) was considerably larger than the total land area of ET decrease (2417.08 km 2 ).There was an increase in ET in Ben Tre, Tien Giang, and Dong Thap provinces with the green shown on the Figure 5a.Besides, a raising trend of PET pattern (shown in red and yellow) mostly distributed in Long An, Dong Thap, Tien Giang, and Ben Tre with a total area of PET increase of 12,537.47km 2 , while a declining area of 8897.18 km 2 concentrated in a part of coastal areas, Ca Mau and Kien Giang.Regarding the spatiotemporal pattern of EVI changes, it is highlighted that, across the delta, the total area of significant EVI increase was 4.2 times greater than the area of significant EVI decrease (3976.9km 2 versus 16,757.07km 2 ).This can be explained by the transformation of rice cultivation to dragon fruit, pineapples, and other kinds of annual plants that aims to adapt to natural hazards and contributes to agricultural development.Moreover, these plants grow over several years and their leaf canopy is larger than the paddy fields, leading to a significant EVI increase.
LST increased significantly across an even larger area of the delta-a total area of LST increase of 7784.54 km 2 versus LST decrease of 3849.71km 2 .The area of decrease was mainly concentrated on provinces of lower delta, including Ca Mau, Kien Giang, and Bac Lieu, while the areas of increase were provinces of upper delta.
The amount of decrease in PET and increase in LST was distributed across most of the provinces in the upper delta, while there was a significant variation in ET and EVI between the western and eastern parts of the study area.
Areas of significant changes for each measure at the provincial level are shown in Table 5.All provinces experienced an increasing trend in environmental factors.Several provinces (Dong Thap, Can Tho, and Vinh Long) did not have a decreasing trend in LST.The largest areas of ET and PET increase (which accounted for 67.29% and 66.40% of total increased evapotranspiration area, respectively) were found in Ca Mau, Bac Lieu, Kien Giang, and Soc Trang.The highest area of EVI increase was discovered in Long An, Kien Giang, and Soc Trang, with an enlarged percentage of 39.16%.However, Long An, Ben Tre, and Tien Giang showed the biggest area of LST increase (10.68%) in comparison with the rest of the provinces.The main reason may be due to the switch from paddy rice to tropical fruits, including dragon fruit growing and pineapples [115].These fruit leaves' area is small, and the cultivated distance is about 2.5 meters, so it has an arid area or water surface among them without vegetation.At the provincial level, Ca Mau and Kien Giang provinces had the highest significant increase in ET, PET, and EVI area because these regions were covered by the largest mangrove forest compared to other provinces in the delta [116,117].In contrast, the decreased value in ET, PET, and EVI in Long An and Tien Giang provinces can be found due to changes in land use from rice to dragon fruit and rapid urban sprawl.

The Spatiotemporal Variations of the EDSI in the Study Area
The annual drought risk from 2001 to 2018 was the period used to apply Equations ( 1)-( 6).An EDSI value below zero is treated as an annual change toward drought, while areas with EDSI values above zero were treated as annual change toward wetness.There were changes in the drought index along the study period according to the categories presented in Table 3, and annual drought appeared in the west and then expanded toward the coastal areas in later years (Figure 6).The annual drought increased in intensity and in coverage along some years and decreased along others during the 2001-2016 period.From 2016 to 2018, the level of drought became more severe, concentrated in the south coastal areas, especially in Ca Mau and Kien Giang.It is important to note that the increase in drought area in coastal areas reflects negative human influences [80,81].data can be acquired over a relatively short period of time by remote sensing.In addition, several weaknesses were discovered by Mu et al., in particular that SPI with the single parameter as rainfall is unreliable and may change with long duration and the RDI faces the challenge when used across separate basins and computed seasonally.Besides, the typical features of location are challenges for DSI as a universal drought index.

Conclusions
The spatiotemporal changes of ET, PET, EVI, LST, and EDSI during 2001-2018 were investigated using spatial regression in the Mekong River Delta (MRD).Positive and negative relationships between each parameter were explored across space and time.It was found that the changes in PET and LST are coherent with the EVI change in the study area.Additionally, the combination of the factors reflected better water stress as a response to agricultural drought, compared to other applied drought indices in the MRD.ET/PET is a key factor for surface water and the combination of EVI and LST reflected the soil moisture texture.Moreover, EDSI is composed of components with a direct physical meaning and improves water deficit assessment on ground surface, as compared to the DSI.Therefore, EDSI inherits the advantages of DSI, enhancing soil moisture information, by adding new factors, and was effectively applied in the MRD.
The study computed and mapped the spatiotemporal variation of drought based on the EDSI using MODIS products.This method assessed drought situations in the study area, more broadly than using a single parameter.The effectiveness of the index was shown by examining spatiotemporal correlation among the other drought indices (i.e., SPI, SDI, and DSI).The relationships of both drought indices (i.e., SPI and RDI) are significant, with r values ranging above 0.7.This contributes to consider that the EDSI plays pivotal importance in monitoring drought events, compared to the other applied drought indices in the study area.In addition, the analyzed results showed that drought is becoming more serious in coastal areas-especially in the category of above-moderate drought (D2).It was clearly demonstrated that the regions with larger areas with severe drought index are Ca Mau, Kien Giang, and Ben Tre.Besides, the agricultural drought events occurred mainly on the rice paddy land which accounts for most of the area and is recognized as a major resource of agricultural production in the study site.The government and residents, therefore, should take appropriate actions to minimize or mitigate the negative drought impacts in these areas.Moreover, the testing of the newly developed EDSI for assessing authentic drought conditions at each meteorological station and upon the types of land use was found to be accurate and operational.Consequently, it informs that there is an urgent need to invest an efficient and sustained drought warning system for the MRD.This will enable the local government to prepare better plans for drought mitigation.
These MODIS-derived indices, using a spatial-temporal approach, are suitable for the assessment and monitoring of environmental degradation in connection with drought on a broad scale.These results encourage applying the approach adopted in this study to monitoring spatial distribution and change over time in drought conditions around other tropical deltas in the world, which exhibit similar physical characteristics.The index is also aimed at improving drought index assessment, conducted daily in the world, from the authentic meaning of EDSI.A seasonal assessment, hence, is able to be implemented in the future for ESDI as a process of monitoring and forecasting drought events on a large scale in the MRD.Besides, replacing other study sites with application of the EDSI can be conducted to assess its effectiveness for various regions.

Figure 1 .
Figure 1.Location of the study area and list of meteorological stations in the Mekong River Delta (MRD).

Figure 2 .
Figure 2. The spatiotemporal correlation between EDSI and DSI in the MRD during the time period of analysis using the geographically weighted regression (GWR).

Figure 4 .
Figure 4. GWR analysis between (a) Z RT and Z EVI ; (b) Z RT and Z LST , and (c) Z LST and Z EVI local correlation coefficients (p < 0.01).

Figure 6 .
Figure 6.The spatiotemporal distribution of the annual drought index.Looking closely at the EDSI in 2017, a major drought level occurred across almost the entire study area, demonstrated by the fierce drought experienced in the dry season of 2016-2017 under the combined impact of upstream dams, climate change, and El Niño[80,118,119].Regarding drought impacts, it damaged 10,000 ha of rice in the winter season at a cost of 60 million USD during the period of 2003-2005[33] and this figure for the 2015-2017 period was estimated up to 360 million USD[75].In addition, it was reported that the influenced area, in the dry season of 2015-2016[120], can be detailed as: 238,276 hectares of rice crop, 29,277 hectares of annual plants, 6575 hectares of vegetables, and 79,000 hectares of shrimp farming.The total number of households facing water shortage was estimated by the Ministry of Agriculture and Rural Development at 211,261 (over approximately one million inhabitants) and the cost to implement the temporary solutions for water supply was more than 22.5 million USD[121].

Figure 7 .
Figure 7. Variations in drought levels in different provinces during the 2001-2018 period.The points show the real values in each province and the lines are the polynomial fitting of these points respectively, for each drought category.

Table 1 .
Earth observation data used in the study site.
* Data at different resolution were resampled to 1000 m resolution pixels.

Table 2 .
The relationship between the enhanced drought severity index (EDSI), the drought severity index (DSI) and the reconnaissance drought index (RDI), the standardized precipitation index (SPI) measured at the meteorological stations in the study area during the 2001-2018 period.

Table 3 .
The EDSI dynamic range and relative categories for dry conditions.

Table 4 .
Summary of results when applying GWR analysis.

Table 5 .
Concluding changes of environmental factors in the study area using Theil Sen slope analysis with p < 0.01.