Understanding Rainfall Distribution Characteristics over the Vietnamese Mekong Delta: A Comparison between Coastal and Inland Localities

: This study examines the changing rainfall patterns in the Vietnamese Mekong Delta (VMD) utilizing observational data spanning from 1978 to 2022. We employ the Mann–Kendall test, the sequential Mann–Kendall test, and innovative trend analysis to investigate trends in annual, wet, and dry season rainfall, as well as daily rainfall events. Our results show significant spatial variations. Ca Mau, a coastal province, consistently showed higher mean annual and seasonal rainfall compared to the further inland stations of Can Tho and Moc Hoa. Interestingly, Ca Mau experienced a notable decrease in annual rainfall. Conversely, Can Tho, showed an overall decrease in some months of the wet season and an increase in dry season rainfall. Furthermore, Moc Hoa showed an increase in the number of rainy days, especially during the dry season. Principal component analysis (PCA) further revealed strong correlations between annual rainfall and extreme weather events, particularly for Ca Mau, emphasizing the complex interplay of geographic and climatic factors within the region. Our findings offer insights for policymakers and planners, thus aiding the development of targeted interventions to manage water resources and prepare for changing climate conditions.


Introduction
Southeast Asia consists of 11 countries, most of which are archipelagos [1,2].The complicated distribution of land, sea, and topography causes the climate to change dramatically across time and place [1,2].According to various studies [3,4], heatwaves have increased in recent decades, accompanied by weakened monsoon seasons, resulting in less rainfall and hotter weather.Vietnam is situated on the east coast of the Indochina Peninsula.It is divided from the rest of the peninsula by the Truong Son Mountain Range, also known as the Annamite Range, which runs parallel to the coast [5,6].Vietnam has a characteristic monsoon-influenced climate like that of mainland Southeast Asia.The wide range of topography, long span of latitude (Vietnam spans over 15 • of latitude), and influences from the East Sea result in climatic conditions that differ between regions.The climate is additionally impacted by the El Niño Southern Oscillation, which shapes monsoonal circulation and prompts intricate changes in rainfall and temperature patterns that differ spatially at a subnational level [7][8][9].Many studies also suggest that tropical cyclones affect rainfall patterns in Vietnam.According to Garcia, Vietnam is struck by an average of four to six typhoons per year [10].The central region of Vietnam experiences more tropical cyclones than other regions, peaking in October and November.The maximum number of tropical cyclones occurs in the northern region from July to September, while nontropical cyclones account for most of the total rainfall in the southern region [11].
Vietnam has been classified by the Intergovernmental Panel on Climate Change as one of the country's most likely to be impacted by climate change due to its extensive coastline, vast agriculturally significant river deltas and floodplains, location on the path of typhoons, and large impoverished populations that are dependent on climate-sensitive economies [12][13][14].Changes in rainfall patterns are significant for socioeconomic development planning and national security, in particular for forecasting the onset of droughts and flood events, as well as assessing surface and groundwater water resources and their quality [15].Many studies have highlighted the regional variations in rainfall patterns, notably decreases in the north of Vietnam and increases volumes in the south [16].The distribution and length of the rainy season and the number of rainy days has also changed.
The Vietnamese Mekong Delta (VMD) is a low-lying freshwater ecosystem, a hotspot of biodiversity and an agricultural production zone that contributes to the lives and livelihoods of more than 21 million people.The VMD offers the basis for inclusive socioeconomic development that can help eradicate poverty throughout the region.However, the VMD is additionally one of the regions of the world that is most affected by climate change.Those residing in this delta are particularly exposed to natural hazards that are caused or exacerbated by climate change.The number and intensity of climate-induced disasters, such as floods and droughts, are already increasing noticeably.In the VMD, 5% of the water comes from rainfall.The influence of geography and local natural conditions has affected the spatial distribution of rainfall, particularly the differences between rainfall variation in coastal areas and inland areas.Climate change is making weather events more extreme, not only by altering the mean and variability, making wet seasons wetter and dry seasons drier, but also through the spatial distribution of rainfall [17,18].Currently, little understood about how the characteristics of rainfall vary from coastal to inland areas or how these changes relate to certain largescale oscillations [18].
To evaluate recent climate change trends, non-parameter tests were conducted.Firstly, we applied the Mann-Kendall statistical test to detect monotonic trends in data.A large bevy of previous studies have used this statistical test to ascertain the presence of statistically significant trends in meteorological and climatic timeseries datasets [18][19][20][21][22][23].Furthermore, Sen's slope estimator can yield a trend slope estimate and is comparatively robust to outliers.Although this test is extensively used for trend detection, it has several problems related to serial correlation in timeseries data.To address this issue, modified Mann-Kendall (MMK) and the trend-free pre-whitening Mann-Kendall (TFPW-MK) tests have been created.These tests yield more accurate results than the MK test, but they are also heavily dependent on sample size and data distribution.As a result, innovative trend analysis (ITA) was established to address the shortcomings of the MK test [24,25].Prior studies have reported that ITA provides reliable results over the MK test to detect trends in climate parameters [19,[26][27][28].Furthermore, as of late, innovation trend analysis (ITA) has gained popularity [24], as it is capable of capturing nonlinear trends and offers additional insights into the data's structure, including trend identification across all value ranges [29].Serinaldi et al. recently reported that sample size and data dispersion also affect ITA [30].Thus, it is crucial to investigate rainfall trends using nonparametric methods other than the ITA method in order to better understand rainfall variability by comparing the results to those obtained using other methodologies [25].
We selected three distinct areas of the VMD to analyze changes in rainfall characteristics over time and under climate change conditions.(i) Ca Mau is a coastal plain province of the VMD.The climate in Ca Mau, in recent years, has not been as distinctly divided as before, with the rainy season is no longer concentrated from May to November and the dry season months no longer occurring from December to April of the following year.Now, rain is often observed scattered throughout all months of the year.Ca Mau has a very unique climate, with periods of wetness and dryness occurring in the same year [8].(ii) Can Tho city is located in the mid-lower region of the VMD, while (iii) Moc Hoa is located in the tropical monsoon and hot and humid climate zone.Because of the contiguity between the two regions of the Southeast and the Mekong Delta, this area not only has the typical characteristics of the Mekong Delta but also has the unique characteristics of the Southeast region.The issue of climate change has become a critical concern for the sustainable development of the VMD [25].As spatial and temporal rainfall characteristics are an important part of the hydrologic cycle, their examination will aid the improved understanding of flood, drought, saline water intrusion, and extreme weather [28,[31][32][33].The goal of this study is to evaluate how climate change affects the distribution of spatiotemporal rainfall in the three different areas of the VMD.The results can inform the scientific basis for local managers to draft and plan appropriate adaptation strategies for future changes in rainfall distributions.

Study Area
Datasets from the three meteorological stations were used to assess the changes in rainfall characteristics along the transect from the coastal to inland regions of the VMD (Figure 1).The first selected meteorological station is located in Ca Mau, Vietnam's and the VMD's most southern province.This coastal province is a peninsula encompassed by the sea on three sides and incised by a dense network of rivers and canals.In total, the province has a coastline of over 250 km.The second meteorological station we assessed is Can Tho station in Can Tho City.Can Tho is located in the middle and lower reaches of the VMD.The province stretches for over 55 km along the west bank of the Hau River, with a total natural land area of 1401.61 km 2 .The third station selected was Moc Hoa in Long An province.Long An is located in the transitional zone between Vietnam's southeast and southwest regions.Therefore, Moc Hoa not only exhibits the typical climatic characteristics of the VMD, but also the unique characteristics of the southeast region.

Rainfall Correlation Matrix and Principal Componant Analysis
A Pearson matrix correlation analysis was conducted to examine the rainfall characteristics that influence the annual wet season and dry season rainfall volumes and distributions.We used an accuracy index with the Pearson correlation coefficient as shown in Equation ( 1): Data collection: Daily rainfall data for the 45-year period from 1978 to 2022 was recorded at the three meteorological stations.

Rainfall Correlation Matrix and Principal Componant Analysis
A Pearson matrix correlation analysis was conducted to examine the rainfall characteristics that influence the annual wet season and dry season rainfall volumes and distributions.We used an accuracy index with the Pearson correlation coefficient as shown in Equation ( 1): The Evans standard (1996) [34] was adopted to determine the level of Pearson correlation matrix strengths (Table 1).Here, the R values vary from −1 to 1, with positive values indicating a positive correlation (two variables increase or drop together) and negative values indicating a negative correlation (one variable increases while the other decreases).The CA tool was applied to examine the data for spatial and temporal differences, which previous research has applied to hydrological and meteorological data [35][36][37].The PCA was used to transform the original variables into new principal components performed along directions of maximum variance.Also, reducing the contribution of less significant variables with minimal information loss was identified.In this study, the principal component was applied to identify which factors will be the most important parameters of water quality, as expressed in Equation (2) (see references [38][39][40]): where Z is the component score, a is component loading, x is component number, j is sample number, and m is a total number of variables.

Rainfall Trend Analysis
Because the rainfall data at the stations did not follow a normal distribution, nonparametric tests were employed in this study.We conducted the Mann-Kendall test, the sequential Mann-Kendall test, and innovative trend analysis, which have been used extensively in a number of past studies [23,[41][42][43][44].The Mann-Kendall statistic determines the presence of a trend and whether that trend is positive or negative based on the range of Mann-Kendall statistics (between −1 and 1).The positive values of the test statistics (S) indicate an increasing trend, whereas negative values suggest a falling trend.The Mann-Kendall equation, which is based on the S statistic (Equation ( 3)) [45,46], is as follows: where, x i and x j are sequential data values in timeseries i, j, n is the length of time series, and sign x i − x j is the sign function as shown in Equation (4): The variance is computed as shown in Equation ( 5): where n is the number of data points, g is the number of tied groups, and t p denotes the number of data values in the p th group.If there are no tied groups, this summary process can be ignored [44].A tied group is a set of sample data having the same value.In cases where the sample size n > 30, the standard normal test statistic Z s is computed using Equation ( 6) [23]: Positive values of Z s indicate increasing trends, while negative Z s values show decreasing trends.Testing trends is performed at the specific α significance level.When |Z s | > Z 1−α/2 , the null hypothesis is rejected, and a significant trend exists in the timeseries.Z 1−α/2 is obtained using the standard normal distribution table.In this study, significance levels α = 0.01 and α = 0.05 were used.At the 5% significance level, the null hypothesis of no trend is rejected if |Z s | > 1.96.It is also rejected if |Z s | > 2.576 at the 1% significance level.
Sen's slope estimator (Sen 1968) [47] can be used to estimate the slope or rate of change in hydrometeorological timeseries.The slopes of each pair of data are computed to provide an estimate of the slope Q, as shown in Equation (7): where X j and X k are the data values at times j and k (j > k), respectively.If there is only one datum in each time period, then N = n(n−1)

2
, where n is the number of time periods.If there are multiple observations in one or more time periods, then N < n(n−1)

2
. The n values of Q i are ranked from smallest to largest, and the median of slope or Sen's slope estimator is computed as shown in Equation ( 8): The Q med sign reflects a data trend, while its value indicates the steepness of the trend.To determine whether the median slope is statistically different than zero, one should obtain the confidence interval of Q med at a specific probability.The confidence interval of the time slope [48] can be computed as show in Equation (9): The S-MK test is used to determine if the series is increasing or decreasing with time.While the test displays the data graphically, it also determines the trend's starting point [49].The t value denotes the test statistic, which was calculated by summing the n i values obtained by counting the smaller ones from the previous ranks for each rank (Equation ( 10)).
The mean value, variance Var(t), and test statistic u(t) were calculated as shown in Equations ( 11), (12), and (13), respectively.The backward test statistic u'(t) is calculated similarly to u(t).11) Innovative trend analysis (ITA) was developed by Sen (2012) [24,50].ITA has shown some benefits over the Mann-Kendall test, which have been illustrated via a comparison of the two approaches [51].The timeseries is separated into two similar groups, each of which is classified in ascending order.Following that, both the first and second parts of the timeseries are plotted on the X and Y axes, respectively.If the data are collected on the 1:1 ideal line (45 • line), the timeseries exhibits no trend.When data are placed on the upward triangular edge of the ideal line, the timeseries continues to grow.When data are aggregated in the downward triangular region of the 1:1 line, the timeseries shows a falling trend [24,50].If the scatter points do not correspond to growing or decreasing trends, the data is then divided into "low", "medium", and "high" value groups [52].Thus, this method allows for the obvious identification of low-, medium-, and high-value trends in any climate timeseries [24,26,51].The "low", "medium", and "high" value groups were based on the percentiles [53] according to the low (X < X − S x ), medium ( X − S x < X < X + S x , and high levels (X > X − S x ).X. X denotes mean value, and S x denotes standard deviation.

Descriptive Statistics
In general, the results highlight different annual rainfall distribution volumes for all three assessed stations.Total annual rainfall decreased from the coastal to the inland stations along the transect (Table 2).The mean annual rainfall over the assessment period in Ca Mau, Can Tho, and Moc Hoa was 2360 mm, 1624 mm, and 1595 mm, respectively.The mean annual rainfall in Ca Mau was approximately 31% and 32% higher than that of Can Tho and Moc Hoa, respectively.Similarly, Ca Mau exhibited 30% and 32% more rainfall during the wet season than both Can Tho and Moc Hoa.It should be noted that, the mean rainfall in Ca Mau during the dry season was also approximately 43% and 33% higher than in Can Tho and Moc Hoa.Can Tho also recorded the lowest rainfall of all three stations during the dry season.The maximum daily rainfall and number of days with rainfall exceeding 50 mm were also higher in Moc Hoa than in Can Tho.The observed decrease in total annual rainfall from coastal to inland stations suggests a gradient influenced by geographical factors and, potentially, by monsoonal patterns [54].Hence, due to the direct exposure to monsoon winds, coastal areas like Ca Mau typically receive more rainfall, while inland areas, such as Moc Hoa, receive less [55].

Correlation and Principal Component Analysis (PCA) for Examining Rainfall Characteristics
Table 3 demonstrates that the annual rainfall volumes in Ca Mau had a very strong Spearman correlation with wet season rainfall, at 0.936, strong correlation with daily rainfall above 50 mm, at 0.64, and a moderate correlation with daily maximum rainfall, at 0.533.The wet season rainfall also showed a strong correlation with daily rainfall above 50 mm and maximum daily rainfall (r = 0.611) and moderate correlation with daily maximum rainfall (r = 0.594); however, the dry season rainfall has a strong correlation with dry season rainfall days and moderate correlation with the overall number of rainy days.In general, annual rainfall and wet season rainfall were heavily influenced by the number of daily rainfalls above 50 mm.Therefore, there are significant differences in rainfall between the wet and dry seasons across all stations, which is a characteristic feature of monsoon climates.More specifically, monsoon's influence is evident in the substantial increase in rainfall during the wet season, which is a direct result of the southwest monsoon that affects the region from May to October [55,56].However, rainfall during the dry season was most affected by the number of dry season days.These results indicate that the PCA results would provide more details related to the rainfall characteristics in Ca Mau.
Based on a diverse range of rainfall features in Ca Mau, Figure 2 displays the PCA results that separated the annual rainfall datasets into four major categories.Groups I and II experienced higher annual rainfall than the other categories.Group I often experienced substantial annual rainfall volumes in certain years as a result of La Niña.The La Niña years include 1989, 1995, 1996, 1999, 2000, 2006, 2008, 2011, and 2022 [57].Group II is especially notable because it includes years in which one year is afflicted by both El Niño during the dry season and La Niña during the wet season.This group is characterized by very high rainfall during the wet season and very low rainfall during the dry season.Consequently, Group II still exhibited high rainfall volumes, although they are less than those of Group I. Group III includes the years with the least amount of annual rainfall relative to the other groups.It is significant to note that the amount of rainfall in Group III was most impacted by El Niño in all years with low rainfall.The years when El Niño was pronounced are 1983, 1987, 1991, 2004, 2010, 2014, 2016, and 2019 [58].Consequently, the total annual rainfall volume is minimal because both wet and dry season rainfall volumes were low.In addition, other rainfall characteristics are low in Group III (Figure 2).Group IV shows higher rainfall figures than Group III because the years in this group had a high dry season rainfall, despite low wet season rainfall.Additionally, Figure 2 shows that, although the number of rainy days varied little, wet season rainfall, annual rainfall, daily rainfall exceeding 50 mm, and maximum daily rainfall all varied significantly in Ca Mau.These results highlight the monsoon's sensitivity to broader climatic oscillations [59].El Niño typically weakens the monsoon, leading to drier conditions, while La Niña strengthens it, resulting in more rainfall [60].Moreover, as certain years with higher rainfall volumes are often associated with La Niña events, this further emphasizes the monsoon's role in regional rainfall variability [61].IV shows higher rainfall figures than Group III because the years in this group had a high dry season rainfall, despite low wet season rainfall.Additionally, Figure 2 shows that, although the number of rainy days varied little, wet season rainfall, annual rainfall, daily rainfall exceeding 50 mm, and maximum daily rainfall all varied significantly in Ca Mau.These results highlight the monsoon's sensitivity to broader climatic oscillations [59].El Niño typically weakens the monsoon, leading to drier conditions, while La Niña strengthens it, resulting in more rainfall [60].Moreover, as certain years with higher rainfall volumes are often associated with La Niña events, this further emphasizes the monsoon's role in regional rainfall variability [61].Table 4 and Figure 3 demonstrate that the annual rainfall in Can Tho station (Var 1) was very strongly correlated to wet season rainfall, with a correlation coefficient of 0.923; strongly correlated to daily rainfall over 50 mm (r = 0.678); moderately correlated to the number of rainfall days and the number of wet season rainfall days (with correlation coefficient of 0.590 and 0.554, respectively).
The wet season rainfall (Var 2) showed a strong correlation with daily rainfall over 50 mm (r = 0.751) and moderately correlated with the number of wet season rainfall days, the daily maximum rainfall and number of rainfall days (correlation coefficient of 0,578,   3 demonstrate that the annual rainfall in Can Tho station (Var 1) was very strongly correlated to wet season rainfall, with a correlation coefficient of 0.923; strongly correlated to daily rainfall over 50 mm (r = 0.678); moderately correlated to the number of rainfall days and the number of wet season rainfall days (with correlation coefficient of 0.590 and 0.554, respectively).Values in bold are different from 0, with a significance level of alpha = 0.05.PCA classified the rainfall data series into four categories (Figure 3).Rainfall in Can Tho was less affected by the El Niño and La Niña phenomena than Ca Mau Station, except for Group I. Compared with the other groups, Group I exhibited the highest average annual rainfall (1905 mm), with nearly 80% of the high rainfall occurring during La Niña years.The high rainfall in Group I was attributed to high annual rainfall, high wet season rainfall, a high number of rainfall days, a high number of wet season rainfall days, and a high number of daily rainfall exceeding 50 mm.Whereas, within Group II, 27% of the years were normal, 36% were affected by both neutral and El Niño, and 27% were mixed (either affected by El Niño, La Niña) and 10% of the high rainfall occurred during La Niña years.However, Group II received the second highest annual rainfall after Group I (1758 mm).Groups III and IV received the same amount of annual rainfall and were smaller than Groups I and II, containing fewer years.Groups III (1464 mm) and IV (1421 mm) account for less than 50% of all El Niño years.Within the framework of this investigation, we can see that the El Niño-Southern Oscillation (ENSO) had a smaller impact on rainfall at Can Tho station than it did at Ca Mau.It should be noted that all variables of rainfall characteristics were found to fluctuate over the years.
Table 5 and Figure 4   PCA classified the rainfall data series into four categories (Figure 3).Rainfall in Can Tho was less affected by the El Niño and La Niña phenomena than Ca Mau Station, except for Group I. Compared with the other groups, Group I exhibited the highest average annual rainfall (1905 mm), with nearly 80% of the high rainfall occurring during La Niña years.The high rainfall in Group I was attributed to high annual rainfall, high wet season rainfall, a high number of rainfall days, a high number of wet season rainfall days, and a high number of daily rainfall exceeding 50 mm.Whereas, within Group II, 27% of the years were normal, 36% were affected by both neutral and El Niño, and 27% were mixed (either affected by El Niño, La Niña) and 10% of the high rainfall occurred during La Niña years.However, Group II received the second highest annual rainfall after Group I (1758 mm).Groups III and IV received the same amount of annual rainfall and were smaller than Groups I and II, containing fewer years.Groups III (1464 mm) and IV (1421 mm) account for less than 50% of all El Niño years.Within the framework of this investigation, we can see that the El Niño-Southern Oscillation (ENSO) had a smaller impact on rainfall at Can Tho station than it did at Ca Mau.It should be noted that all variables of rainfall characteristics were found to fluctuate over the years.
Table 5 and Figure 4 emphasize that the annual rainfall in Moc Hoa station (Var 1) has very strong correlation with wet season rainfall (r = 0.938) and daily rainfall exceeding 50 mm (r = 0.896) and a strong correlation with daily maximum rainfall (r = 0.638).Wet season rainfall (Var 2) has a correlation with daily rainfall exceeding 50 mm (r = 0.826) and the daily maximum rainfall (r = 0.624).As expected, dry season rainfall has a correlation with the number of dry season rainfall days at 0.631 (strong).Rainfall in the three regions with similar points depends primarily on the daily rainfall exceeding 50 mm.The second element influencing rainfall at Ca Mau and Moc Hoa is the daily maximum rainfall, whereas Can Tho station is affected by the number of rainfall days.Values in bold are different from 0, with a significance level of alpha = 0.05.Group I comprises all years with La Niña, and thus showed the most rainfall (average annual rainfall 2000 mm).Except for daily maximum rainfall, all rainfall variables were higher than those of the other groups.Group II showed the highest daily rainfall at 126 mm/day, whereas Group I had 115 mm/day.Overall, 92% of the years in Group II coincided with El Niño years, although the group exhibited the second highest average rainfall Rainfall in the three regions with similar points depends primarily on the daily rainfall exceeding 50 mm.The second element influencing rainfall at Ca Mau and Moc Hoa is the daily maximum rainfall, whereas Can Tho station is affected by the number of rainfall days.
Group I comprises all years with La Niña, and thus showed the most rainfall (average annual rainfall 2000 mm).Except for daily maximum rainfall, all rainfall variables were higher than those of the other groups.Group II showed the highest daily rainfall at 126 mm/day, whereas Group I had 115 mm/day.Overall, 92% of the years in Group II coincided with El Niño years, although the group exhibited the second highest average rainfall after Group I. Group III had the lowest rainfall, with 47%, 36%, and 18% for El Niño years, normal years, and La Niña years (1316 mm), respectively.Group IV is split 50/50 between El Niño and La Niña years.
Annual rainfall and rainfall during the wet season varied greatly over three locations.The number of wet season rainfall days at Ca Mau station, which is located in the coastal region, varied slightly.Furthermore, El Niño and La Niña have a greater influence on rainfall in coastal areas than inland localities.Furthermore, rainfall in Can Tho showed strong fluctuations in the number of rainy days.

Rainfall Trend Test
Figure 5 shows that the mean cumulative rainfall at all three study stations decreased during the period from 2000 to 2022 when compared with the period from 1978 to 2001.The mean cumulative rainfall in Ca Mau, Can Tho, and Moc Hoa fell by 14%, 4%, and 5%, respectively.Whereas, for Can Tho station, we observed that rainfall increased during the dry season.It should be noted that rainfall at Ca Mau station decreased in both the rainy and dry seasons.At Moc Hoa station, rainfall during the rainy season was seen to decrease, while it increased slightly during the dry season.
To investigate the rainfall characteristics, we conducted the Mann-Kendall trend test, the sequential Mann-Kendall test, and innovative trend analysis to better understand rainfall trends, including (i) the rainfall amount in each of the following seasons: annual, wet season, dry season, and the monthly (ii) daily maximum rainfall, daily rainfall exceeding 50 mm, and (iii) the number of days with rainfall (annually, wet and dry seasons).We also discovered that only the annual rainfall and wet season rainfall in Ca Mau had declined significantly at the 95% and 90% confidence level, respectively, using the Mann-Kendal test (Figure A1).August rainfall and Jun rainfall also had declined significantly at the 95% and 90% confidence level, respectively.The results of Sen's Slope calculation show the rate of decrease in annual rainfall, wet season rainfall, August rainfall, and Jun rainfall were found to be 8.25 mm/year, 6.8 mm/year, 3.2 mm/year, and 2.5 mm/year, respectively.However, the dry season rainfall and number of dry season rainfall day in Can Tho increased significantly at the 95% and 90% confidence level, respectively (Figure A2).The results of Sen's Slope calculation show the rate of increase in dry season rainfall and number of dry season rainfall days were 2.7 mm/year and 18 day/decade, respectively.In Can Tho station, we also found the increased significantly in August rainfall and Jun rainfall at the 95% and 90% confidence level, respectively.The results of Sen's Slope calculation show the rate of decrease in August rainfall and Jun rainfall were 1.8 mm/year and 1.77 mm/year.In Moc Hoa station, we found the number of rainfall days and the number of dry season rainfall days in Moc Hoa had increased significantly at the 95% and 90% confidence level, respectively (Figure A3).The results of Sen's Slope calculation show the rate of increase in the number of rainfall days and the number of dry season rainfall days were 3.2 day/decade and 2.5 day/decade, respectively.
Figure 5 shows that the mean cumulative rainfall at all three study stations decreased during the period from 2000 to 2022 when compared with the period from 1978 to 2001.The mean cumulative rainfall in Ca Mau, Can Tho, and Moc Hoa fell by 14%, 4%, and 5%, respectively.Whereas, for Can Tho station, we observed that rainfall increased during the dry season.It should be noted that rainfall at Ca Mau station decreased in both the rainy and dry seasons.At Moc Hoa station, rainfall during the rainy season was seen to decrease, while it increased slightly during the dry season.To investigate the rainfall characteristics, we conducted the Mann-Kendall trend test, the sequential Mann-Kendall test, and innovative trend analysis to better understand rainfall trends, including (i) the rainfall amount in each of the following seasons: annual, Figure 6 displays the trends of rainfall characteristics for the three assessed meteorological stations at a 95% confidence level using S-MK and ITA.
Between 1978 and 2022, the amount of annual rainfall in Ca Mau significantly decreased by 95% (Figure 6a).We discovered that the amount of annual rainfall overall has declined since 2007, with this decline being statistically significant since 2016.Innovative trend analysis classified this decrease in annual rainfall as "low" to "medium" level because most of the points plot is below the 1:1 line and based on the percentiles.The annual rainfall volumes in both the rainy and dry seasons decreased between 1978 and 2022.The annual rainfall volumes fell more during the wet season than during the dry season.This drop in rainfall during the wet season was statistically significant at 90%, whereas the drop during the dry season was not statistically significant.The annual rainfall at Can Tho station varied from year to year, with no discernible trend.We noticed that the overall wet season rainfall volume decreased during 1978-2022, but this is not statistically significant.Meanwhile, dry season rainfall volumes increased at the 95% confidence level over 45 years.The drop in wet season rainfall was because of a decrease in the number of rainy days during the wet season, and the increase in dry season rainfall was likewise because of an increase in the number of rainy days during the dry season.According to ITA results, the drop in wet season rainfall from 1978 to 2022 varies from "medium" to "high" levels, whereas the increase in dry season rainfall was at the "medium" level.
tion show the rate of decrease in August rainfall and Jun rainfall were 1.8 mm/year and 1.77 mm/year.In Moc Hoa station, we found the number of rainfall days and the number of dry season rainfall days in Moc Hoa had increased significantly at the 95% and 90% confidence level, respectively (Figure A3).The results of Sen's Slope calculation show the rate of increase in the number of rainfall days and the number of dry season rainfall days were 3.2 day/decade and 2.5 day/decade, respectively.Between 1978 and 2022, the amount of annual rainfall in Ca Mau significantly decreased by 95% (Figure 6a).We discovered that the amount of annual rainfall overall has declined since 2007, with this decline being statistically significant since 2016.Innovative trend analysis classified this decrease in annual rainfall as "low" to "medium" level because most of the points plot is below the 1:1 line and based on the percentiles.The annual rainfall volumes in both the rainy and dry seasons decreased between 1978 and 2022.The annual rainfall volumes fell more during the wet season than during the dry season.This Although annual rainfall in Moc Hoa decreased between 1978 and 2022, the change was not statistically significant.However, we found that the number of rain days increased at the 95% confidence level.This increase occurred mainly during the dry season, with a 95% confidence level.On the contrary, rainfall in the wet season decreased at a 90% significance level because of a decrease in rainfall intensity.Aside from that, we discovered a drop in the number of days with rainfall that exceeded 50 mm, while the number of rainy days in the wet season remained unchanged.To understand how the pattern of rain changes over the rainy season, we investigated the pattern of rainfall in both Ca Mau and Moc Hoa.We considered both seasons in relation to Can Tho station.Figure 7 shows a 95% statistically significant decrease in rainfall volume for the month of August for Ca Mau.Ca Mau station: In general, August rainfall decreased with a statistical significance of 95%.In particular, the rainfall in the month of August has decreased significantly since 1986 and has increased again thereafter.Since 1991, the rainfall in August has decreased, and the rainfall has decreased significantly again since 2007.Results from the IAT and S-MK show that the August rainfall decreased at "medium" levels during the period from 1978 to 1998 and at the "low" level after 1991.
Can Tho station: Rainfall in the dry season months increased, but this is insignificant.It should be noted that rainfall increased in January and February with a 95% confidence level, according to the Mann-Kendall test results.However, there was no trend in the S-MK test results.Consequently, additional monitoring time is needed over these two months to precisely estimate future changes in rainfall.Just for the month of April, we found evidence of an increasing trend at the 90% significance level, whereas June rainfall decreased at a 90% confidence level.There was no trend detected for the months of March, July, October, or November.Monthly rainfall tended to decrease in the majority of wet season months at Can Tho station (June to October).Only August showed a reduction in rainfall at the 95% confidence level (Figure 8).It should be mentioned that the early rainy season months (i.e., June) also declined, but only by 90%.This should be considered when observing fluctuations in rainfall patterns in Can Tho.Ca Mau station: In general, August rainfall decreased with a statistical significance of 95%.In particular, the rainfall in the month of August has decreased significantly since 1986 and has increased again thereafter.Since 1991, the rainfall in August has decreased, and the rainfall has decreased significantly again since 2007.Results from the IAT and S-MK show that the August rainfall decreased at "medium" levels during the period from 1978 to 1998 and at the "low" level after 1991.
Can Tho station: Rainfall in the dry season months increased, but this is insignificant.It should be noted that rainfall increased in January and February with a 95% confidence level, according to the Mann-Kendall test results.However, there was no trend in the S-MK test results.Consequently, additional monitoring time is needed over these two months to precisely estimate future changes in rainfall.Just for the month of April, we found evidence of an increasing trend at the 90% significance level, whereas June rainfall decreased at a 90% confidence level.There was no trend detected for the months of March, July, October, or November.Monthly rainfall tended to decrease in the majority of wet season months at Can Tho station (June to October).Only August showed a reduction in rainfall at the 95% confidence level (Figure 8).It should be mentioned that the early rainy season months (i.e., June) also declined, but only by 90%.This should be considered when observing fluctuations in rainfall patterns in Can Tho.
At Moc Hoa station, the number of dry season months declined, although not significantly at the 95% confidence level.Remarkably, the rainfall only declined in May at the 90% significance level.Comparing Moc Hoa to the Ca Mau and Can Tho stations, there has been minimal variation in rainfall for the assessment period from 1978 to 2022.
decreased at a 90% confidence level.There was no trend detected for the months of March, July, October, or November.Monthly rainfall tended to decrease in the majority of wet season months at Can Tho station (June to October).Only August showed a reduction in rainfall at the 95% confidence level (Figure 8).It should be mentioned that the early rainy season months (i.e., June) also declined, but only by 90%.This should be considered when observing fluctuations in rainfall patterns in Can Tho.

Discussion
By analyzing how local climatic factors influence rainfall distribution, the findings from this study highlight the significant spatiotemporal variability in rainfall across the three selected meteorological stations in the VMD.The distinct seasonal rainfall distribution [62], which is characterized by significant differences between the coastal and inland stations [63], highlights the spatial heterogeneity [64], which is significantly associated with monsoon's role in shaping local climate patterns [65].Specifically, the significant changes observed in seasonal rainfall (such as the decreased August rainfall at Ca Mau and the increased dry season rainfall at Can Tho) align with the known impacts of monsoonal shifts and changes in ENSO phases [66].Further, the maximum daily rainfall and the number of days with rainfall exceeding 50 mm describing the variations in the intensity and frequency of rainfall events are evidence of the role of monsoon dynamics in the region [67].The variability observed in our study, particularly the effects of El Niño and La Niña events, depicts broader global climate variability and its manifestation at local scales.The ENSO phenomenon plays a pivotal role in modulating global climate patterns, including the monsoon system in the VMD [68].Our results show a clear correlation between ENSO events and rainfall variability, with La Niña years associated with higher rainfall and El Niño years associated with reduced rainfall [69].By examining these correlations, our study provides insights on the local impacts of global climatic oscillations, as well as contributes to the understanding of how monsoon regions globally might respond to climate change.
Our study's findings on the annual and seasonal rainfall distributions across three meteorological stations in the VMD align with broader observations of hydroclimatic variability observed across the Mekong River Basin (MRB) [70].The observed rainfall patterns, which are characterized by significant spatial variability and a notable influence of monsoon dynamics, are consistent with the complex hydrological regime that governs the MRB [71,72].Comparative analyses with regional studies, such as those conducted by the Mekong River Commission (MRC) and various scientific and policy-oriented studies, reveal a consistent description of changing hydrological conditions across the MRB [73][74][75].These changes are attributed to both natural climatic fluctuations, including the impact of El Niño and La Niña phenomena [76], and anthropogenic factors, such as upstream dam construction [77,78] and water resource management practices [79].Hence, our study's observation of changed rainfall distributions and intensities can be interpreted as part of the broader hydrological shifts reported across the basin, which have significant implications for water security [80], agricultural productivity [81], and ecosystem health [82] in the VMD.
Furthermore, our analysis aligns with findings from studies highlighting the adverse impacts of climate change and human activities on the MRB's flow regimes.The joint efforts by the MRC and other regional bodies to address these challenges through data sharing [83] and collaborative research initiatives [84] provide a critical context for interpreting our study's results.Specifically, as discussed above, the correlation between rainfall patterns in the VMD and broader hydrological changes underscores the interconnectedness of local and regional climate systems.The hydrological shifts due to climate change have been responsible for the heightened frequency of floods and the corresponding economic damage in the region [85].Recent analyses have underscored the recurrent nature of floods in the VMD and the significant economic losses attributed to them [86].Specifically, floods, which are driven by sea level rise and increased precipitation during the wet season, pose substantial economic challenges, highlighting the necessity for reliable loss and risk assessments [87][88][89].In Can Tho City, for instance, floods have demonstrated that inundation duration significantly impacts the economy, with residential buildings and their contents being considerably affected [90].This highlights the need for effective flood risk management strategies that account for the unique flooding characteristics of the VMD, where households have adapted to long inundation durations and shallow water depths [91].Moreover, studies have emphasized the acute vulnerability of the VMD to coastal flooding, which is exacerbated by high tides and heavy rainfall [85,92].A noteworthy event detailed the inundation of streets in Can Tho City to a depth of 50 cm, damaging 91 hectares of crops.This event serves as a potent reminder of the region's perennial risk of flooding [93].Furthermore, the 2011 flood event in the Mekong Delta was among the most severe in recent history, resulting in 265 fatalities, damage to at least 449,000 houses, and an estimated economic loss of USD $600 million [93].These figures highlight the profound impact of flooding on the region, necessitating comprehensive strategies to mitigate future risks [94].Hence, these findings, which illustrate the variability in rainfall patterns across the VMD, offer crucial insights for policymakers, environmental planners, and local communities.By understanding the dynamics of rainfall distribution and its implications on flood frequency and severity, targeted interventions can be developed to enhance flood preparedness and response strategies [95][96][97][98][99][100][101][102].This is particularly relevant given the observed economic impacts of flooding events, which highlight the urgent need for improved flood loss modeling and risk assessments in the region [102].
Besides flooding, the variability in rainfall distribution and trends identified in our study also reflects the concerns raised by experts regarding water security in the VMD.The reported decrease in silt, sand, and nutrients [103,104], alongside the increasing incidents of saline intrusion [105] and coastal erosion [106], emphasizes the pressing need for a comprehensive strategy to ensure water security and sustainable development in the region.These findings call for enhanced cooperation among riparian countries, informed by accurate and timely hydrological data [107,108], to mitigate the adverse effects of climate variability and anthropogenic pressures on the MRB's hydrology.Hence, our study contributes to the growing body of evidence highlighting the significant impact of monsoon variability, climate change, and human activities on the hydrological dynamics of the VMD and the larger MRB.By providing detailed insights into rainfall patterns in the VMD, our research underscores the importance of region-specific studies in understanding the broader implications of hydroclimatic changes.Future research needs to focus on integrating local observations with regional climate models to improve predictions and inform water management strategies that can address the challenges posed by climate change and sustainable development in the MRB.

Conclusions
This study has provided a comprehensive analysis of rainfall distribution patterns across the Vietnamese Mekong Delta (VMD) over a 48-year period, uncovering significant spatial and temporal variability.Our investigation revealed a declining trend in total annual rainfall, with coastal regions, particularly Ca Mau, experiencing higher annual rainfall volumes compared with inland areas such as Can Tho and Moc Hoa.This spatial variability underscores the coastal environment's significant influence on rainfall distribution within the VMD.Moreover, our analysis detected a significant decrease in annual and wet season rainfall in Ca Mau, alongside a marked reduction in rainfall during specific months like August and June, highlighting the potential impact of regional climatic factors, including the monsoon system.This study's findings also point to the nuanced effects of global climate phenomena such as El Niño and La Niña on the region's hydrology, which is evidenced by contrasting trends in wet and dry season rainfall volumes.Notably, an increase in dry season rainfall volumes, particularly in inland areas, suggests a shifting climatic pattern, which is potentially influenced by the aforementioned global phenomena.This shift has significant implications for water resource management, agricultural planning, and flood risk mitigation in the VMD.While the analysis presented in this study offers valuable insights into rainfall trends and their implications for the VMD, it is acknowledged that the scope was limited to data from three meteorological stations.This limitation, coupled with the exclusion of other climatic factors such as temperature, humidity, and wind patterns, suggests the need for a more holistic approach in future research.Investigating the interplay between rainfall and other climatic variables, alongside their socioeconomic impacts, could provide a more comprehensive understanding of the region's climatic dynamics.In light of these findings, it is imperative for policymakers and stakeholders to consider the complex interplay of local, regional, and global climatic factors in their efforts to adapt to and mitigate the impacts of changing rainfall patterns in the VMD.The role of monsoons, the influence of El Niño and La Niña events, and the recurrent challenge of floods necessitate a strategic approach to water resource management and disaster preparedness.Future research should aim to incorporate predictive models to forecast climatic changes and their potential impacts, thereby enabling more effective planning and adaptation strategies for the region.Hence, this study not only highlights the critical need for enhanced climatic monitoring and analysis in the VMD but also calls for an integrated approach to understanding and managing the complex climatic and hydrological challenges facing the region.

Figure 1 .
Figure 1.Location of three studied meteorological stations in the VMD.

Figure 2 .
Figure 2. Contribution of the variables and observation of rainfall at Ca Mau station.

Figure 2 .
Figure 2. Contribution of the variables and observation of rainfall at Ca Mau station.

Figure 3 .
Figure 3. Contribution of the variables and observation of rainfall at Can Tho station.
emphasize that the annual rainfall in Moc Hoa station (Var 1) has very strong correlation with wet season rainfall (r = 0.938) and daily rainfall exceeding 50 mm (r = 0.896) and a strong correlation with daily maximum rainfall (r = 0.638).Wet season rainfall (Var 2) has a correlation with daily rainfall exceeding 50 mm (r = 0.826) and

Figure 3 .
Figure 3. Contribution of the variables and observation of rainfall at Can Tho station.The wet season rainfall (Var 2) showed a strong correlation with daily rainfall over 50 mm (r = 0.751) and moderately correlated with the number of wet season rainfall days, the daily maximum rainfall and number of rainfall days (correlation coefficient of 0.578, 0.46 and 0.424, respectively); however, the dry season rainfall (Var 3) had a very strong correlation with the number of dry season rainfall days and moderate correlation with the number of rainfall days.It can be seen that annual rainfall and wet season rainfall have a strong correlation with daily rainfall exceeding 50 mm while dry season rainfall has a very strong correlation with number of dry season rainfall days (r = 0.813).PCA classified the rainfall data series into four categories (Figure3).Rainfall in Can Tho was less affected by the El Niño and La Niña phenomena than Ca Mau Station, except for Group I. Compared with the other groups, Group I exhibited the highest average annual rainfall (1905 mm), with nearly 80% of the high rainfall occurring during La Niña years.The high rainfall in Group I was attributed to high annual rainfall, high wet season rainfall, a high number of rainfall days, a high number of wet season rainfall days, and a high number of daily rainfall exceeding 50 mm.Whereas, within Group II, 27% of the years were normal,

Figure 4 .
Figure 4. Contribution of the variables and observation of rainfall at Moc Hoa station.

Figure 4 .
Figure 4. Contribution of the variables and observation of rainfall at Moc Hoa station.

Figure 5 .
Figure 5.The comparison of daily rainfall and accumulated rainfall during the periods from 1978 to 2000 and 2011 to 2022 for Ca Mau (a,d), Can Tho (b,e), and Moc Hoa (c,f).

Figure 5 .
Figure 5.The comparison of daily rainfall and accumulated rainfall during the periods from 1978 to 2000 and 2011 to 2022 for Ca Mau (a,d), Can Tho (b,e), and Moc Hoa (c,f).

Figure 6 Figure 6 .
Figure 6 displays the trends of rainfall characteristics for the three assessed meteorological stations at a 95% confidence level using S-MK and ITA.

Figure 6 .
Figure 6.Rainfall trend patterns: (a) Annual rainfall decreasing trend in Ca Mau using S-MK.(b) Annual rainfall decreasing trend in Ca Mau using ITA.(c) Increased dry season rainfall in Can Tho using S-MK test.(d) Increased dry season rainfall in Can Tho using ITA.(e) Increased number of dry season rainfall days in Moc Hoa using S-MK test.(f) Increased number of dry season rainfall days in Moc Hoa using ITA.

Figure 7 .
Figure 7. Decreased August rainfall at Ca Mau station using S-MK test (a) and ITA (b).(a) Decreased August rainfall at Ca Mau station using S-MK test.(b) Decreased August rainfall at Ca Mau station using ITA.

Figure 7 .
Figure 7. Decreased August rainfall at Ca Mau station using S-MK test (a) and ITA (b).(a) Decreased August rainfall at Ca Mau station using S-MK test.(b) Decreased August rainfall at Ca Mau station using ITA.

Figure 8 .
Figure 8. Decreased August rainfall in Can Tho station using the S-MK test (a) and ITA (b).(a) Decreased August rainfall in Can Tho using S-MK test.(b) Decreased August rainfall in Can Tho using ITA.

Figure 8 .
Figure 8. Decreased August rainfall in Can Tho station using the S-MK test (a) and ITA (b).(a) Decreased August rainfall in Can Tho using S-MK test.(b) Decreased August rainfall in Can Tho using ITA.

Figure A2 .
Figure A2.Rainfall trend patterns in Can Tho station using Mann-Kendall test: (a) Dry season rainfall day increasing trend in Can Tho using Mann-Kendall test with a significance level of 95%.(b) No. dry season rainfall day increasing trend in Can Tho using Mann-Kendall test with a significance level of 90%.(c) August rainfall decreasing trend in Can Tho using Mann-Kendall test with a significance level of 95%.(d) Jun rainfall decreasing trend in Can Tho using Mann-Kendall test with a significance level of 90%.

Figure A2 .
Figure A2.Rainfall trend patterns in Can Tho station using Mann-Kendall test: (a) Dry season rainfall day increasing trend in Can Tho using Mann-Kendall test with a significance level of 95%.(b) No. dry season rainfall day increasing trend in Can Tho using Mann-Kendall test with a significance level of 90%.(c) August rainfall decreasing trend in Can Tho using Mann-Kendall test with a significance level of 95%.(d) Jun rainfall decreasing trend in Can Tho using Mann-Kendall test with a significance level of 90%.

Figure A3 .
Figure A3.Rainfall trend patterns in Moc Hoa station using Mann-Kendall test: (a) Number of rainfall days increasing trend in Moc Hoa using Mann-Kendall test with a significance level of 95%.(b) Number of dry season rainfall days increasing trend in Moc Hoa using Mann-Kendall test with a significance level of 95%.

Table 1 .
Strengths of the Pearson correlation matrix.

Table 2 .
Summary statistics of rainfall during the period from 1978 to 2022 for the three meteorological stations.

Table 3 .
Correlation matrix showing Spearman correlation analysis between rainfall characteristics in Ca Mau station.
Values in bold are different from 0, with a significance level of alpha = 0.05.

Table 4 and
Figure

Table 4 .
Correlation matrix showing Spearman correlation analysis between rainfall characteristics in Can Tho station.
Values in bold are different from 0, with a significance level of alpha = 0.05.

Table 5 .
Correlation matrix showing Spearman correlation analysis between rainfall characteristics in Moc Hoa station.
Values in bold are different from 0, with a significance level of alpha = 0.05.

Table 5 .
Correlation matrix showing Spearman correlation analysis between rainfall characteristics in Moc Hoa station.