Assessing the Effects of Drought on Rice Yields in the Mekong Delta

: In contrast to other natural disasters, droughts may develop gradually and last for extended periods of time. The World Meteorological Organization advises using the Standardized Precipitation Index (SPI) for the early identiﬁcation of drought and understanding of its characteristics over various geographical areas. In this study, we use long-term rainfall data from 14 rain gauge stations in the Vietnamese Mekong Delta (1979–2020) to examine correlations with changes in rice yields. Results indicate that in the winter–spring rice cropping season in both 2016 and 2017, yields declined, corresponding with high humidity levels. Excessive rainfall during these years may have contributed to waterlogging, which in turn adversely affected yields. The results highlight that not only drought, but also humidity has the potential to adversely affect rice yield.


Introduction
Droughts occur in virtually all climatic zones, and cause adverse impacts to natural habitats, ecosystems, society and the economy [1].Droughts have been considered as a silent hazard [2][3][4][5].Several recent studies have indicated that agricultural production has been found to decline significantly due to increased drought frequency and severity [4,5].For example, the severe Australian drought in 2006 reduced the national winter cereal crop by 36%, leaving many farmers in a financial crisis [6].Drought in the Kingdom of Eswatini in 2015-2016 caused a 30% reduction in income [6].Cunha et al. (2019) notes that most regions of Brazil have faced extreme droughts over the last 60 years [7].Furthermore, 26 drought events in the USA from 1980 to 2019 have resulted in estimated losses of at least USD 249 billion [8].
Droughts also severely influence the socioeconomic development of the Lower Mekong Basin (LMB) counties of Thailand, Cambodia, Laos and Vietnam, and negatively affect the livelihoods of about 85% to 90% of people living in the rural areas of the region.The LMB experienced extreme drought events in 1992, 1999, 2003, and 2015-2016, which led to massive economic losses due to crop damage and adverse effects on the environment and farmers' livelihoods.The most extreme drought, in 2016, caused 1.7 billion USD of damages in Thailand and caused water shortages in 18 out of 25 provinces, affecting 2.5 million people, in Cambodia.In Vietnam, the cost of the drought in 2016 was estimated at 669 million USD, and the estimated recovery costs amounted to approximately 1.5 billion USD [9].
In the Vietnamese Mekong Delta (VMD), drought has significantly reduced rice yields.Unlike other natural disasters, drought develops slowly, develops in a ladder-like pattern, and is difficult to predict [10].As a result, a drought index is required to assess the effects of drought.Many drought indicators are currently utilized both domestically and internationally to assess the level of drought and the level of moisture [11].For varied topographic and hydrological conditions, individual indices may be appropriate.As a result, to improve the accuracy of drought detection in areas with numerous water sources, various indicators can be used.The standardized precipitation index (SPI), which is based on rainfall over several time periods, was chosen by the World Meteorological Organization for the assessment of meteorological drought in 2009 [1,12].Additionally, different levels of constraint can be assessed using the SPI indicator [10].Furthermore, the SPI can be used to assess different levels of drought [10].Rainfall data are frequently recorded in a variety of locations, making it easy to calculate the SPI.Meteorological drought was assessed by using SPI [13,14] and the impacts on the development and growth of rice were evaluated [11].Additionally, SPI can be predicted using a climate change dataset.Climate change data were utilized in a study by Ty et al. to calculate the SPI for the VMD [15].Their findings indicated that while the frequency of drought will not increase until the year 2047, it will become increasingly more severe in the VMD, where the probability of a drought occurring across the entire region ranged from 23% to 34% [15].Their study also revealed that dry and wet conditions could be correlated with the average annual rice yield, as well as seasonal rice yields [2].Extreme wet conditions are common in the coastal areas of the VMD under the influence of tropical storms, and can occur for several months.In particular, Ca Mau province exhibits this characteristic [16].Meanwhile, the inland provinces of the VMD have been found to mostly be affected by prolonged drought [17,18].
Droughts also strongly influenced saline intrusion in the VMD; the storage in Tonlé Sap and flow to Kratié are two main upstream factors that have impacts on the VMD's downstream water supply and saline intrusion.Major salinity events, as well as droughts, occurred in the VMD during years of climatic extremes [17,20].In addition, El Niño phenomena impacted the dry season of 2015-2016, causing the rainy season to arrive late and end early, with a decreased volume of rainfall.Additionally, low discharge volumes in the Mekong River have allowed saline waters further upstream.The timing, intensity, and frequency of climate and weather extremes in Vietnam specifically, and in Southeast Asia generally, including the Mekong River basin, have noticeably increased according to the examination of hydrometeorological data from 1975 to 2020.As a result, increased frequency and magnitude of salinity and drought events in the VMD are causing increased economic and environmental repercussions [21,22].
To identify a timely pathway for drought adaptation, it is therefore vital to recognize and understand the spatial and temporal degrees of drought.It is established that the SPI is a somewhat appropriate index for the topographical and climatic features of the VMD.The region's drought state is therefore evaluated using the SPI.The main objectives of this study were to estimate SPI for rainy periods of three, six, and nine months, as well as for the summer-autumn rice cropping season (April to July) and autumn-winter season (December to March), to assess agricultural dryness and create isodrought maps for each season.

Study Area
With a total area of up to 40.6 thousand square kilometers (12% of the country's total area), the VMD is home to 13 provinces and cities with a combined population of 17,273,630 (18% of the nation's total population) in the year 2019 Figure 1.The population density was 1.45 times greater than the national average at 423 people per square kilometer.The VMD contributes roughly 50% of the nation's rice production, 95% of its export rice production, about 65% of its aquaculture production, 60% of its export fish, and 70% of its types of fruits [23].Furthermore, rice is grown on more than 3.2 million hectares, accounting for more than 50% of global rice production and more than 90% of global rice exports each year [24].The predicted area for autumn-winter rice production in 2021 was 719.7 thousand ha, a loss of 4.4 thousand ha from the autumn-winter crop in the previous year, due to the conversion to fruit crops and the impact of insufficient seeds being sowed for the summer-autumn crop in time for rice to be harvested [25].

Estimated SPI for Each Period and Rice Cropping Seasons
Numerous drought indices, such as SPI, SPEI, and ARI, have been developed to describe drought and help with the quantitative evaluation of drought severity.Some of these have been selectively applied in some areas of the VMD [16,[28][29][30][31].Among the indices, SPI has been widely used worldwide [2,[4][5][6][7]11,14,32,33].Therefore, the SPI was selected to calculate the impact of drought on agriculture and establish its frequency using time steps of 3, 6, and 9 months, as well as for the seasons (summer-autumn and winterspring) for the entire VMD.According to SPI 3, SPI 6, and SPI 9 months, summer-autumn and winter-spring crops, the current climate is arid.For each period, the SPI was estimated and the hierarchy of humidity and drought based on the SPI computed and graded in accordance with Formula (1) and Table 1, respectively.The VMD contains abundant water resources (the Tien and Hau rivers), but these are unevenly distributed according to the seasons, and influenced by tides in three directions (East Sea, West Sea and the area bordering the East and West seas) through the river system.In the dry season, the average flow is about 2500 m 3 /s, yet there are times when it is lower than 1500 m 3 /s.The impacts of hydropower, economic activities, excessive groundwater exploitation and climate change have led to a decrease in upstream flows reaching the VMD, while a lack of water in the dry season increases the risk of drought and saltwater intrusion.The discharge volumes of the Mekong River water flowing into the VMD decrease significantly in March or April, leading to saline intrusion [26].The dry season of 2015-2016 was considered to be a high salinity season, but in the dry season of 2019-2020, a salinity drought began to appear from mid-December 2019, nearly 1 month earlier than in the dry season of 2015-2016 and more than 3 months earlier compared with the long-term average [26].
Both drought and saltwater intrusion affected production in the VMD during the 2015-2016 dry season.According to [27], the salinity drought in 2019-2020 lasted for nearly two-thirds of the dry season and remained at a high level throughout the year [27].At its peak, in 2016, areas 25 km inland were affected; but in 2019 the area affected extended more than 40 km inland.Notably, water with a salinity of 4‰ penetrated inland to distances of 100 km along the Hau River, 75 km along the Ham Luong River, 70 km along the Co Chien River, and almost 60 km along the Cua Tieu River, the Cua Dai River and the Cai Long River, respectively.

Estimated SPI for Each Period and Rice Cropping Seasons
Numerous drought indices, such as SPI, SPEI, and ARI, have been developed to describe drought and help with the quantitative evaluation of drought severity.Some of these have been selectively applied in some areas of the VMD [16,[28][29][30][31].Among the indices, SPI has been widely used worldwide [2,[4][5][6][7]11,14,32,33].Therefore, the SPI was selected to calculate the impact of drought on agriculture and establish its frequency using time steps of 3, 6, and 9 months, as well as for the seasons (summer -autumn and winterspring) for the entire VMD.According to SPI 3, SPI 6, and SPI 9 months, summer-autumn and winter-spring crops, the current climate is arid.For each period, the SPI was estimated and the hierarchy of humidity and drought based on the SPI computed and graded in accordance with Formula (1) and Table 1, respectively.Assuming a precipitation series at a given timescale X = x 1 , x 2 , x 3 , . . ., x n , n = 1, 2, 3, . . ., n represents the length of the series, and the distribution function f (x, θ) can best describe the distribution of the precipitation series X.
where θ is the parameter set of the function f (x) that can be estimated based on the sample series X.
The cumulative distribution function can be expressed as: Equation ( 2) can be used to determine the cumulative probability of each observed precipitation event, xi.The SPI then undergoes an equiprobability transformation from the cumulative probability to the unit variance, zero mean, standard normal random variable Z.
where Φ −1 (•) is the inverse of a standard normal cumulative distribution function Φ (•) with zero mean and unit variance.Considering the fact that some kinds of distribution functions are not defined for x = 0, such as the Gamma function, an adapted statistic H (x) can be expressed as Equation (4).
Climate 2023, 11, 13 5 of 17 where q is the probability of no precipitation, which can be computed by simply dividing the number of months (N) by the number of months (n) that have no precipitation (Vicente-Serrano, 2006) [34].Equation ( 5) can then be used to describe the SPI computation process.
In this study, we used the SPI-based drought classification for Vietnam, in which SPI = −0.25 signifies the onset of drought conditions [17,35].Drought begins when the SPI first falls below the negative value of 0.25 or less, and ends when the SPI returns to a positive value of 0.25 or higher (Table 1).

Temporal Analysis
The application of the seasonal Mann-Kendall test was used for the piezometric time series data to enable further data interpretation for drought conditions.The magnitude of the slope was computed using the Sen's slope test as shown in Equation ( 6) [36].
Kendall's S statistic S i for each season is summed to form the overall statistic, S k .Mann-Kendall statistics provide an indication of whether a trend exists, and if so, whether the trend is positive or negative, based on a range between −1 and 1, for each pair of data points depending on the sign of the difference between the values.The Mann-Kendall equation, which is based on the value of the S statistic, is calculated using the equation shown as Equation (7).
where x i and x j are sequential data values, n is the length of time series, and: Mann (1945) and Kendall (1975) [37,38] stated that when n ≥ 8, S is almost normally distributed with the following mean and variance: Value of S at near zero, positive, and negative imply no significant, upward, and downward trends, respectively.The magnitude of S measures the strength of the trend.Note that if the absolute value of S is greater than the critical value of S, this indicates that trends are statistically significant.
where t i is the number of ties of the extent i, and n is the number of observations in the set.The standard Z sk statistic is calculated as follows: Climate 2023, 11, 13 6 of 17 where Z sk pursues the standard normal distribution with µ = 0 and δ = 1.The Sen's slope was used to determine the magnitude of the trend.The Z statistic is formed and compared with a critical point from the standard normal distribution.The magnitude of a time series trend was evaluated by a simple nonparametric procedure developed by Sen; T i is between any two values of a time series (which gives a total of n(n−1)/2 pairs of data could be calculated by: where x j and x i are values at time j and i (1 < i < j < n), n is the number of data.The median of these n values of T i is represented as Sen's estimator of slope.Positive and negative values of T i indicate upward and downward trends in the time series, respectively.

Spatial Analysis of Drought Status
Agglomerative Hierarchical Clustering (AHC) is a common approach applied in many previous studies, using the squared Euclidean distance and Ward's method [39][40][41].This technique can group single objects into clusters with similarity based on certain criteria, making data simple and easy to comment on.The fourteen stations were grouped using hierarchical cluster analysis.As stations of similarity, Ward's approach, employing squared Euclidean distances, was used [39].The findings of the cluster analysis were interpreted using the dendrogram.The linkage distance, expressed as the dissimilarity level, was used to calculate the geographic variability of water quality using AHC.The connectivity distance shown on the x axis was then standardized by multiplying the quotient by 100 [39,40].

Temporal Analysis Using SPI 12−Month Drought
The results of the temporal analysis using the SPI 12−month drought data over the VMD are shown in Figure 2. The Mann-Kendall test showed that although the severity of the drought tended to decrease at coastal stations, with the exception of Soc Trang station, the trend was not statistically significant.Meanwhile, interior meteorological stations showed reduced rainfall, suggesting that the degree of drought is increasing.Among the stations (Tri Ton, Moc Hoa, Can Tho, Cang Long, and Chau Doc), only Tri Ton station showed an increased drought trend at the computed p-value, lower than the significant level alpha = 0.05.

Cluster Analysis Using AHC Technique
Compared with SPI 6 and SPI 9, SPI 3 highlighted more variation and susceptibility to extremes (too dry or too humid) (Figure A1).Figures 3-5 display the findings of a cluster analysis for the three consecutive periods of 3 months, 6 months, and 9 months.According to the cluster analysis results, there were the same number of clusters and stations in each cluster for the 3, 6, and 9-month drought periods (SPI 3, SPI 6, and SPI 9).These were Cluster 1 (Ca Mau), Cluster 2 (Bac Lieu, Ba Tri, Cao Lanh, and Vinh Long), and Cluster 3 (Cang Long, Can Tho, Chau Doc, Moc Hoa, My Tho, Rach Gia, Soc Trang, Vi Thanh, Tri Ton).Cluster 1 has only one station, Ca Mau, which is the most unique station in the Mekong Delta, experiencing both extreme drought and extreme flooding.The drought severity in Cluster 2 tends to decrease, although this decrease is not statistically

Cluster Analysis Using AHC Technique
Compared with SPI 6 and SPI 9, SPI 3 highlighted more variation and susceptibility to extremes (too dry or too humid) (Figure A1).Figures 3-5 display the findings of a cluster analysis for the three consecutive periods of 3 months, 6 months, and 9 months.According to the cluster analysis results, there were the same number of clusters and stations in each cluster for the 3, 6, and 9-month drought periods (SPI 3, SPI 6, and SPI 9).These were Cluster 1 (Ca Mau), Cluster 2 (Bac Lieu, Ba Tri, Cao Lanh, and Vinh Long), and Cluster 3 (Cang Long, Can Tho, Chau Doc, Moc Hoa, My Tho, Rach Gia, Soc Trang, Vi Thanh, Tri Ton).Cluster 1 has only one station, Ca Mau, which is the most unique station in the Mekong Delta, experiencing both extreme drought and extreme flooding.The drought severity in Cluster 2 tends to decrease, although this decrease is not statistically significant, according to the Mann-Kendall test results.Cluster 3 contains stations where the severity of drought fluctuates over time.

SPI 3-Month, SPI 6-Month, SPI 9-Month, and SPI 12-Month Drought Patterns
AHC analysis for SPI 3, SPI 6 and SPI 9 gave similar results with three clusters: Cluster 1 (Ca Mau), Cluster 2 (Bac Lieu, Ba Tri, Cao Lanh, Vinh Long), Cluster 3 (Cang Long, Can Tho, Chau Doc, Moc Hoa, My Tho, Rach Gia, Soc Trang, Vi Thanh, Tri Ton). Figure A2 shows the seasonal SPI for each meteorological station in the VMD, and Figures 6 and  7 show the mean seasonal SPI and mean seasonal rice yields in the VMD.AHC analysis for SPI 3, SPI 6 and SPI 9 gave similar results with three clusters: Cluster 1 (Ca Mau), Cluster 2 (Bac Lieu, Ba Tri, Cao Lanh, Vinh Long), Cluster 3 (Cang Long, Can Tho, Chau Doc, Moc Hoa, My Tho, Rach Gia, Soc Trang, Vi Thanh, Tri Ton). Figure A2 shows the seasonal SPI for each meteorological station in the VMD, and Figures 6 and 7 show the mean seasonal SPI and mean seasonal rice yields in the VMD.
The seasonal mean SPI in the VMD shows that a slight to moderate drought occurred in both the summer-autumn and winter-spring rice cropping seasons during both 1978 and 1995.During the period from 1996-2008, our data show that most of the summer-autumn and winter-spring seasons were moderately wet periods.However, there was a partial drought in 2002 for both summer-autumn and winter-spring seasons.The period from 2009-2019 displays a clear difference between the summer-autumn and autumn-winter seasons; specifically, the winter-spring season was wet, while the summer-autumn season exhibited drought.These results also support Lien's similar research [19].AHC analysis for SPI 3, SPI 6 and SPI 9 gave similar results with three clusters: Cluster 1 (Ca Mau), Cluster 2 (Bac Lieu, Ba Tri, Cao Lanh, Vinh Long), Cluster 3 (Cang Long Can Tho, Chau Doc, Moc Hoa, My Tho, Rach Gia, Soc Trang, Vi Thanh, Tri Ton). Figure A2 shows the seasonal SPI for each meteorological station in the VMD, and Figures 6 and  7 show the mean seasonal SPI and mean seasonal rice yields in the VMD.The seasonal mean SPI in the VMD shows that a slight to moderate drough in both the summer-autumn and winter-spring rice cropping seasons during and 1995.During the period from 1996 -2008, our data show that most of the autumn and winter-spring seasons were moderately wet periods.However, t partial drought in 2002 for both summer-autumn and winter-spring seasons.T from 2009-2019 displays a clear difference between the summer-autumn and winter seasons; specifically, the winter-spring season was wet, while the su tumn season exhibited drought.These results also support Lien's similar resea Winter-spring rice yields declined in 2016 and 2017 at the same time as hi ity levels (Figure 5).The excessive rainfall during this time may have contribu Winter-spring rice yields declined in 2016 and 2017 at the same time as high humidity levels (Figure 5).The excessive rainfall during this time may have contributed to waterlogging, which in turn may have affected rice yield.Years of poor rice yield in both the summer-autumn and the autumn-winter crops first began in 2001 and 2002, corresponding with the presence of mild drought in both seasons (Figure 5).In 2009, the summer-autumn rice yield corresponded with slight drought.In contrast to the overall trend, rice production was low in the research area for the two consecutive years (2016 and 2017) which coincided with a slight drought in the summer-autumn crop, but high humidity in the winter-spring crop.The area affected by the drought was not taken into account in this study.However, Lien et al. [19] demonstrated that the occurrence of drought corresponds to times when the VMD experienced recent losses in agricultural production due to both drought and saltwater intrusion, such as in the years 2002, 2004-2005, 2006-2007, and 2009-2010 [19].Another study also revealed that rice crops were highly affected by drought and saline intrusion in the VMD [42].These results are also consistent with the timing of yield reduction during drought.Therefore, it can be seen that not only drought, but also humidity potentially reduces rice yields.
The SPI, SPI3, SPI6, and SPI9 analysis of the summer-autumn and winter-spring seasons reveals that the overall pattern of drought in the VMD ranged from low to medium drought in the period 1978-2015, while wet conditions occurred in the period 2016-2008, and mild drought from 2008-2019.In order to examine drought variance through time and space, the study mapped 3 periods using SPI 12-month drought patterns (Figure 8).The West coastal region showed moisture levels varying from slightly moist to severely wet throughout each of the three periods during 1980 and 2019, whether experiencing drought or wet periods.Regularly occurring slight to moderate droughts affected the northwest hinterlands of the VMD.Thereafter, the extent of the drought slowly subsided and conditions became wet towards the East and West seas.Between 1978 and 1995, a large drought spread from the interior to the East Coast and then progressively subsided toward the West Coast.The spatial distribution of drought during the period 1996-2008 is comparable to that of 1978-1995; however, during this time, the distribution of drought in both the interior and coastal East Sea areas was decreasing, while the West Coast areas were becoming more humid and with increased humidity levels.SPI 12 provides a general overview of how wetness and drought are distributed in the VMD.It has been pointed out that seasonal variation led to a change in rice seeding time to ensure productivity [43], and the reduction of precipitation caused a drop in rice yield [44].Therefore, differences in the level of seasonal variation and in 3-month, 6-month, or 9-month droughts, as well as the crop's growth stages, need to be carefully investigated.

Conclusions
The fluctuation of SPI values in consecutive periods of 3-, 6-, and 9-month drought differs in the three groups: Cluster 1 (Ca Mau), Cluster 2 (Bac Lieu, Ba Tri, Cao Lanh, and Vinh Long), and Cluster 3 (Can Tho, Can Tho, Chau Doc, Moc Hoa, My Tho, Rach Gia, Soc Trang, Vi Thanh, and Tri Ton).These results demonstrate that not only physical, but also geographical factors influence the distribution of drought in the VMD.Furthermore, drought in the VMD varies seasonally.Linking drought events to seasonal SPIs reveals that the seasonal SPI is responsible for the majority of average yield reduction events in the VMD.Rice productivity was shown to be affected not only by drought, but also by high humidity.
The VMD's pattern of drought generally demonstrates that both rice cropping seasons experienced drought from 1978-1995.High humidity levels were present in both the summer-autumn and winter-spring harvests between 1996 and 2005.It should be noted that there is a distinction between the two seasons from 2006-2019: summer-autumn (c) The West coastal region showed moisture levels varying from slightly moist to severely wet throughout each of the three periods during 1980 and 2019, whether experiencing drought or wet periods.Regularly occurring slight to moderate droughts affected the northwest hinterlands of the VMD.Thereafter, the extent of the drought slowly subsided and conditions became wet towards the East and West seas.Between 1978 and 1995, a large drought spread from the interior to the East Coast and then progressively subsided toward the West Coast.The spatial distribution of drought during the period 1996-2008 is comparable to that of 1978-1995; however, during this time, the distribution of drought in both the interior and coastal East Sea areas was decreasing, while the West Coast areas were becoming more humid and with increased humidity levels.SPI 12 provides a general overview of how wetness and drought are distributed in the VMD.It has been pointed out that seasonal variation led to a change in rice seeding time to ensure productivity [43], and the reduction of precipitation caused a drop in rice yield [44].Therefore, differences in the level of seasonal variation and in 3-month, 6-month, or 9-month droughts, as well as the crop's growth stages, need to be carefully investigated.

Conclusions
The fluctuation of SPI values in consecutive periods of 3-, 6-, and 9-month drought differs in the three groups: Cluster 1 (Ca Mau), Cluster 2 (Bac Lieu, Ba Tri, Cao Lanh, and Vinh Long), and Cluster 3 (Can Tho, Can Tho, Chau Doc, Moc Hoa, My Tho, Rach Gia, Soc Trang, Vi Thanh, and Tri Ton).These results demonstrate that not only physical, but also geographical factors influence the distribution of drought in the VMD.Furthermore, drought in the VMD varies seasonally.Linking drought events to seasonal SPIs reveals that the seasonal SPI is responsible for the majority of average yield reduction events in the VMD.Rice productivity was shown to be affected not only by drought, but also by high humidity.
The VMD's pattern of drought generally demonstrates that both rice cropping seasons experienced drought from 1978-1995.High humidity levels were present in both the summer-autumn and winter-spring harvests between 1996 and 2005.It should be noted that there is a distinction between the two seasons from 2006-2019: summer-autumn crops tend to be dry, but winter-spring crops are typically wet.However, considering the geographical characteristics and comparing these results taking different spatial scales into consideration will give these research findings more credibility and relevance; this is a limitation for this work and can be considered for future study.Moreover, there is a strong positive association between drought and salinity.This emphasizes the importance of holistic management plans looking into the system dynamics.

Climate 2023 , 19 Figure 1 .
Figure 1.The study area, with locations of the 14 meteorological stations in the 13 provinces of the VMD.

Figure 1 .
Figure 1.The study area, with locations of the 14 meteorological stations in the 13 provinces of the VMD.

Figure 2 .
Figure 2. SPI 12−month drought trends over the 14 meteorological stations of the VMD.The x axis denotes SPI 12−month drought.

Figure 2 .
Figure 2. SPI 12−month drought trends over the 14 meteorological stations of the VMD.The x axis denotes SPI 12−month drought.

Figure 5 .
Figure 5. Continuous term of 9 months, SPI 9-month drought divided into 3 clusters by AHC analysis.

Figure 6 .Figure 5 .
Figure 6.The arithmetic mean SPI of summer-autumn and winter-spring rice cropping seasons

Figure 6 .Figure 6 .Figure 7 .
Figure 6.The arithmetic mean SPI of summer-autumn and winter-spring rice cropping seasons across 14 meteorological stations in the VMD for the period 1978-2019.

Figure 7 .
Figure 7.The arithmetic means of summer-autumn and winter-spring rice yields in the 13 provinces in the VMD.

Table 1 .
The classification of drought and wet periods according to the SPI.

Table 1 .
The classification of drought and wet periods according to the SPI.
Continuous term of 3 months, SPI 3-month drought divided into 3 clusters by AHC analysis.Continuous term of 6 months, SPI 6-month drought divided into 3 clusters by AHC analysis.Continuous term of 6 months, SPI 6-month drought divided into 3 clusters by AHC analysis.
Figure 3. Continuous term of 3 months, SPI 3-month drought divided into 3 clusters by AHC analysis.Figure 3. Continuous term of 3 months, SPI 3-month drought divided into 3 clusters by AHC analysis.