Using Remote Sensing and Multivariate Statistics in Analyzing the Relationship between Land Use Pattern and Water Quality in Tien Giang Province, Vietnam

: This study was carried out to understand how land use patterns influence surface water quality in Tien Giang Province using remote sensing and statistical approaches. Surface water qual ‐ ity data were collected at 34 locations with the frequency of four times (March, June, September, and November) in 2019. Water quality parameters were used in the analysis, including pH, temper ‐ ature, electrical conductivity (EC), total suspended solids (TSS), dissolved oxygen (DO), biological oxygen demand (BOD), chemical oxygen demand (COD), ammonium (N ‐ NH 4+ ), nitrite (N ‐ NO 2 − ), nitrate (N ‐ NO 3 − ), sulfate (SO 42 − ), orthophosphate (P ‐ PO 43 − ), chloride (Cl − ), total nitrogen (TN), total phosphorus (TP), and coliform. The relationship between land use patterns and water quality was analyzed using geographic information techniques (GIS), remote sensing (RS), statistical ap ‐ proaches (cluster analysis (CA), principal component analysis (PCA), and Krustal–Wallis), and weighted entropy. The results showed water quality was impaired by total suspended solids, nu ‐ trients (N ‐ NH 4+ , N ‐ NO 2 − , P ‐ PO 43 − ), organic matters (BOD, COD), and ions (Cl − and SO 42 − ). Kruskal– Wallis analysis results showed that all water quality parameters in the water bodies in Tien Giang Province were seasonally fluctuated, except for BOD and TN. The highest levels of water pollutants were found mostly in the dry season (March and June). The majority of the land in the study area was used for rice cultivation (40.64%) and residential (27.51%). Water quality in the study area was classified into nine groups corresponding to five combined land use patterns comprising residen ‐ tial–aquaculture, residential–rice cultivation, residential–perennials, residential–rice–perennial, and residential–rice–perennial crops–aquacultural. The concentrations of the water pollutants (TSS, DO, BOD, COD, N ‐ NH 4+ , N ‐ NO 2 − , Cl − , and coliform) in the locations with aquaculture land use patterns (Clusters 1 and 2) were significantly larger than those of the remaining land use patterns. PCA analysis presented that most of the current water quality monitoring parameters had a great impact on water quality in the water bodies. The entropy weight showed that TSS, N ‐ NO 2 − , and coliform are the most important water quality parameters due to residential–aquaculture and resi ‐ dential–rice cultivation; EC, DO, N ‐ NH 4+ , N ‐ NO 2 − , Cl − , and coliform were the significant variables for the land use type of residential–perennial crops; N ‐ NO 2 − , P ‐ PO 43 − , and coliform for the land use pattern of residential–rice cultivation–perennial crops) and N ‐ NH 4+ , N ‐ NO 2 − , Cl − , and coliform for the land use pattern of residential–rice cultivation–perennial crops–aquaculture. The current find ‐ ings showed that that surface water quality has been influenced by the complex land use patterns in which residential and rice cultivation may have major roles in causing water impairment. The results of the water quality assessment and the variation in water properties of the land use patterns found in this study provide scientific evidence for future water quality management. water quality was indicated by the higher concentrations of water pollutants in the dry season than the wet season, except for TN and BOD. Five major land uses including resi ‐ dential–aquaculture, residential–rice cultivation, residential–perennials, residential rice– perennial, and residential–rice–perennial–aquacultural were identified in Tien Giang Province, in which residential and rice cultivation were the majorities. CA analysis di ‐ vided the sampling sites in the study areas into nine groups distributed on five land use patterns. Water quality at Cluster 1 and Cluster 2 (the locations with aquaculture land use patterns) was lower than that in the other land use patterns. The PCA analysis results show that five potential water ‐ polluting sources resulting in water quality variation and all the current water quality monitoring variables were significant. The entropy weight calculation results show that TSS, N ‐ NO 2 − , and coliform were the important water param ‐ eters in the water bodies within the land use types of residential–aquaculture and resi ‐ dential–rice cultivation; EC, DO, N ‐ NH 4+ , N ‐ NO 2 − , Cl − , and coliform were significant in the residential–perennial land use type. Water quality in the land area of residential–rice cultivation–perennial crops was contaminated by N ‐ NO 2 − , P ‐ PO 43 − , and coliform. Finally, the parameters of N ‐ NH 4+ , N ‐ NO 2 − , Cl − , and coliform were the main concerns for water quality in mixed land use patterns of residential–rice cultivation–perennial crops–aqua ‐ culture. The current results reflect the temporal and spatial variation of surface water quality in various land use patterns in Tien Giang Province, and this information could be very useful for future water quality management.


Introduction
In the Vietnamese Mekong Delta, rivers can be considered as important water resources serving a variety of purposes, including domestic, industrial, service, irrigation, aquaculture, etc. On the other hand, rivers play a crucial role in transporting numerous pollutants originated from human and natural activities [1]. Therefore, information on river water quality and influencing factors is essential to form an effective water management system. In recent years, factors affecting water quality such as land use/land cover, human activities, topographical characteristics, climate, and other natural processes [2] have become the concern of water resources managers [3][4][5].
In water quality monitoring in Vietnam, the water quality index (WQI) is frequently used to examine changes in water quality; however, this indicator cannot prove sources of water pollution because the index is calculated by normalizing the analytical values [6]. Moreover, most of the approaches to understanding surface-water-polluting sources in Vietnam are mainly based on the results of field surveys and interviews. Meanwhile, there have been several studies of the methods applied to assess and identify sources of water quality fluctuations [7][8][9][10][11][12]. Typically, land use and soil type can be particularly useful in explaining the relationship between water quality and environmental characteristics in the study area. This can be explained by the fact that land use and soil type have a direct impact on the permeation and dispersion rates of pollutants causing fluctuations in water quality [13]. Several recent studies have also shown that different land use patterns can influence the extent and causes of water pollution, the timing of water quality, and the balance of the ecosystems in water bodies [14][15][16]. However, the relationships between land use and water quality are not always consistent, depending on nature-and humanrelated activities in water bodies [17,18]. Therefore, prediction of potential sources and contaminants of surface water should be performed using multivariate statistical methods. Among the multivariate statistical methods, principal component analysis (PCA) and cluster analysis (CA) have been successfully used in the previous studies for the identification of key water contaminants, potential polluting sources, and classification of water quality [19][20][21].
Socioeconomic development in Tien Giang Province heavily depends on agriculture, which could lead to negative changes in water quality if water management policies are improperly incorporated into the development strategies. Therefore, reliable information on water quality characteristics is needed for the effective management of water resources. However, available information on the relationship between land use and surface water quality in the water bodies in Tien Giang Province is still limited. Therefore, the main objectives of this study are (1) to characterize surface water quality in the water bodies in Tien Giang Province; (2) to analyze the relationship between land use patterns and water quality; and (3) to identify important water quality parameters and potential sources of pollution influencing on surface water quality in the study area. The findings from the current study provide useful information for effective water quality management in Tien Giang Province and in Vietnam.

The Study Area
The study was conducted in Tien Giang, a coastal province in the Vietnamese Mekong Delta, with a coastline of 32 km long. The study area is a province located along the Tien River and is also the end of the flow of the Mekong River before it flows into the East Sea. The area is about 2510.5 km 2 . The province has a flat terrain with a slope <1% (West to East) and the elevation ranges from 0 to 1.6 m above sea level (0.8-1.1 m) [22]. In particular, the alluvial soil area accounts for 53% of the total area, accounting for most of the areas in Cai Be, Cai Lay, Chau Thanh, Cho Gao, My Tho, and a part of Go Cong Tay District belonging to the areas with freshwater sources. Regarding the soil, approximately 19.4% of the study area is the acid sulfate soil group distributed in the districts of Cai Be, Cai Lay, Tan Phuoc; 14.6% is alluvial soil that is salty. Meanwhile, sand dune soil group accounts for about 3%, scattered in the districts of Cai Lay, Chau Thanh, Go Cong Tay, and mostly concentrated in Go Cong Dong District; some other land groups account for about 10%. Annual rainfall varies from 1100 mm to 1400 mm, averaging about 1175 mm. The annual average temperature is about 27 °C. In addition, due to its proximity to the sea, this makes river basins in Tien Giang Province (coastal rivers) vulnerable to climate change hazards, particularly sea-level rise [22]. The speed of tidal transmission in estuaries is about 30 km h −1 (1.5 times higher than that of the Hau River). This may have caused salinity intrusion in some water bodies near the sea, with water backflow from the sea to the river about 1.2 m s −1 . Typically, some areas are likely to be salty in Tien Giang Province, including the districts of Tan Phu Dong, Go Cong Dong, Go Cong Tay, Cho Gao, the town of Go Cong, and the city of My Tho. Therefore, monitoring of water quality in the study area has been of interest to the locality and researchers.

Land Use/Land Cover Classification
Land use/land cover data for the study area were extracted from remote sensing images using the QGIS 3.16 classification method (https://qgis.org/en/site/forusers/download.html, accessed on 16 December 2020) [26]. The study collected Landsat 8 images with a spatial resolution of 30 × 30 m from the available United States Geological Survey (USGS) database (https://earthexplorer.usgs.gov, accessed on 16 December 2020). In addition, the digitization of the map creates geospatial data areas in shapefile format for zoning the study area. After using Google Earth and the land use planning map of Tien Giang Province, the land use/land cover data have been further revised to reflect the current land use status for 2019. Land use/land cover is classified into six main categories, including water bodies, aquaculture land, rice paddy land, perennial cropland, forest land, and residential land.

Assessment of Surface Water Quality
For the purpose of assessing the statistically significant differences in surface water quality between the observed months, a test was used to check the validity of the dataset. In particular, the Shapiro-Wilk test was performed instead of the Kolmogorov-Smirnov test because the sample size (number of locations) in the study was less than 50 [27]. Then, a one-way ANOVA analysis was used to evaluate a statistically significant difference if the dataset had a normal distribution (significance level greater than 0.05); on the contrary, if the dataset did not have a normal distribution, then a Kruskal-Wallis (hierarchical oneway ANOVA) analysis was performed. Kruskal-Wallis is described as an alternative nonparametric test for ANOVA analysis. The results were presented in tabular form with the values including mean ± standard deviation (SD) and the notations of the difference with a confidence level of 95%.

Assessment of Surface Water Quality
In this study, cluster analysis (CA) was used to consider whether the water quality monitoring sites can be spatially grouped based on the similarity of water quality characteristics or not. This analysis was performed using Euclidean distance as a measure of similarity/dissimilarity [28]. Specifically, the clusters with the highest similarity would be grouped first, which means that the Euclidean distance is the lowest. Similarly, the similarity level of the clusters would decrease, and eventually, all clusters would be combined into a single cluster. The number of clusters in the study was determined by examining the results of the dendrogram, specifying the ratio of the quotient between the bonding distance of a particular cluster, and the maximum bonding distance ((Dlink/Dmax) × 100 < 60) was considered significant [1]. Each cluster shows its own water quality; in addition, CA analysis could be based on any number of variables, and these variables could be of any type.
Principal component analysis (PCA) was used to identify factors or sources of variation in surface water quality in the study area, in which eigenvalues was used as a measure of change and was explained by each PC [29]; the PCs were selected when eigenvalue coefficients were greater than 1 [30]. The contribution and importance of the initial parameters are shown in the PC through the weighting factors (loading). Furthermore, water samples affected by different sources and processes would have different water quality, which was indicated by a variety of water quality parameters. Therefore, sources of pollution were determined based on the weighting coefficient of water quality parameters in the main components [31]. The higher the weighting coefficient was, the greater the influence [32]. According to research by Liu et al. (2003) [33], loading factors were classified as "strong," "medium," and "weak," with absolute loading values of >0.75, 0.75-0.50, and 0.50-0.30, respectively. In addition, the entropy weight calculation of the surface water quality parameters in each cluster was carried out for a specific analysis in determining the main parameters of the areas affected by land use/land cover.
The weighted entropy was calculated as follows: Assume that m was the water quality evaluation parameter, and n was the monitoring position. Rij is the correlation of the ith parameter at position j.
where Cij is the concentration of the parameter i at the monitoring site j.
In most cases, the data size (number of replicates) for each contaminant is different. Therefore, the data should then be standardized before calculating the informational entropy value. The normalized values of the data and the information entropy were calculated using Equation (2) and Equation (3), respectively.
Standardized data: where pịj is the standardized data of the parameter i at site j.
The value of informational entropy (Hi): Finally, the weighted entropy (wi) were obtained from Equation (4) as follows: All calculations and statistical analysis were performed using Excel version 2016 software (Microsoft Crop., Washington, DC, USA), SPSS version 20 software (IBM Crop., Armonk, NY, USA) and Primer software version 5 (Primer-E Ltd., Plymouth, UK). Table 1 provides information on the mean values and standard deviations of the 16 water quality parameters over four observational periods. The pH value was seasonally changed with a tendency to shift gradually from March (7.2 ± 0.5) to September (7.8 ± 0.6). Climatic characteristics and seawater intrusion into 45-50 km (from the estuary) area of the region can be considered the main cause of pH high in September. On the other hand, the average temperature fluctuated from 29.0 ± 1.4 to 31.2 ± 1.0 °C, which tended to decrease gradually until September and increased again in November. The temperature difference may be partly due to the increase in air temperature in the dry season and decreasing in the rainy season. A previous study showed that any change in land cover, depth, and disturbance of runoff could cause temperature fluctuation [34]. In addition, temperature fluctuation depends on the time of the sampling. However, pH and temperature do not affect physical, chemical, and biological processes in water and are still in a suitable range for aquatic life [35,36].

Water Quality Characteristics of Tien Giang Province in 2019
TSS concentration was usually the lowest in the dry season and the highest in the rainy season [37]. This was also consistent with the fluctuation of TSS in surface water bodies in Tien Giang Province with TSS ranging from 77.7 ± 21.3 to 121.8 ± 50.1 mg L −1 . The reported TSS concentration in 2014-2015 ranged from 32.5 to 57.4 mg L −1 [38], which was lower than that found in the current study. The increase in TSS may have led to a number of obstacles in using water for domestic activities and for water treatment. DO and EC values over the survey were in the ranges of 3.2 ± 1.1-4.0 ± 1.2 mg L −1 , 432 ± 665.7-992.8 ± 1273.6 μS cm −1 , respectively, and DO value tended to be lower than that in the period of 2014-2015 (4.9-5.1 mg L −1 ) [38,39]. Consequently, water quality was impaired resulting in constraints for the activities of organisms. In contrast to TSS, EC and DO values in March and June were higher than those in September and November, and this was similar to the previous report observing EC and DO in water bodies in Ben Tre [40]. High EC values during March and June can be attributed to lower river water volume and higher temperature facilitating ionization [41]. Furthermore, the high amount of organic matter in water has led to a decline in DO concentration in the water [35,42]. This decrease was evidenced by the relatively high levels of BOD and COD with the values ranged from 8.0 ± 2.7 to 8.9 ± 2.6 mg L −1 and from 14.4 ± 4.8 to 17.3 ± 4.4 mg L −1 , respectively. The concentrations of BOD and COD were higher than that of 2014-2015 (7.3-11.0 mg L −1 and 11.6-18.0 mg L −1 for BOD and COD, respectively) [38,39]. There were statistically significant differences between seasons for COD (p < 0.05); however, no significant changes were observed in all four observations of BOD (p > 0.05). Since COD mainly consists of BOD, the variation of COD also has the same time variation as for BOD. COD increased in June and decreased in September, then tended to increase again in November. High levels of BOD and COD in this study indicated that the water bodies were organically polluted.  In this study, the concentrations of N-NH4 + , N-NO2 − , and N-NO3 − were in the ranges of 0.3 ± 0.4-0.5 ± 0.5 mg L −1 , 0.0 ± 0.0-0.2 ± 0.3 mg L −1 , 0.1 ± 0.1-0.4 ± 0.2 mg L −1 , respectively. There was seasonal variation, in which nitrogen species increased from March to June and then decreased in November. Variation in nitrogen-derived nutrients may be related to water biology, seasonal, and types of cultivated crops. N-NO3 − concentration remained at a safe level, i.e., not harmful to human health, at the pH value of 6.5-8.5 [43]. However, nutrient pollution caused by N-NO2 − has been found at a relatively high level [44]. When organic matter content is high, nitrate can be reduced to some degree to nitrite, which could account for the high concentrations of this pollutant during June and September [45]. On the other hand, river phosphate is not toxic to humans, animals, or fish and is a limiting factor in eutrophication [46]. The P-PO4 3− concentration ranged from 0.1 ± 0.2 to 0.1 ± 0.1 mg L −1 , indicating the risk of causing eutrophication in the water bodies in the study area. Meanwhile, TN concentration ranged from 3.1 ± 1.3 to 3.8 ± 1.5 mg L −1 , with the highest found in March and the lowest in November. N-NO3 − concentration is the main parameter that causes TN concentration in river water changes [47]. However, no significant difference was found for TN (p > 0.05), with the mean concentration ranging from 3.1 ± 1.3 to 3.8 ± 1.5 mg L −1 . Similar to the other parameters, TP concentration increased from March (0.1 ± 0.2 mg L −1 ) to June (0.3 ± 0.2 mg L −1 ) and then decreased gradually until November (0.14 ± 0.2 mg L −1 ). This was also reported in the previous study, TP concentrations in the months of the rainy season were higher than those in the dry season [48]; This is a different feature compared to the TN parameter. In general, an excess of nutrient concentration during the study period can affect the growth of aquatic plants, especially algae.
The concentrations of Cl − and SO4 2− ions in September and November (rainy season) were lower than that in the dry season, and there were statistically significant differences between the dry and wet months (p < 0.05). This observed result was found to be similar for electrical conductivity values. The concentrations of Cl − and SO4 2− ions ranged between 86.0 ± 235.7 and 366.4 ± 1104.5 mg L −1 , and 39.0 ± 15.4 and 67.6 ± 41.5 mg L −1 , respectively. In addition, based on the values of pH, EC, Cl − , and SO4 2− in March, it was shown that the water bodies in Tien Giang Province in 2019 were heavily affected by saline intrusionthis is considered a special case of water pollution and has been reported in many other coastal areas [49]. In addition, the discharge of wastewater could result in an increase in the concentration of these two ions in the water [43]. Cl − concentrations in the range below 70 mg L −1 are considered safe; however, Cl − can have harmful effects on agriculture if it is higher than 350 mg L −1 [50]. Therefore, Cl − concentrations in water bodies in June, September, and November may be unsafe and Cl − in March was considered unsuitable for irrigation.
Finally, the coliform density through the observational periods tended to increase gradually until the end of the year (November), ranging from 972.9 ± 718.3 to 2261.2 ± 1349.0 MPN 100 mL −1 . According to some previous reports of coliform in the water bodies in Dong Thap, Hau Giang, Can Tho, and Soc Trang, coliform density in the surface water in Tien Giang Province tended to be lower [37,[51][52][53]. It can therefore be shown that the water bodies were less likely to be directly affected by domestic wastewater and waste from livestock and humans. Therefore, microbiological pollution may not be the main pollution problem in the water bodies in Tien Giang Province.
In summary, the results of water quality analysis during the four observational periods revealed that the river water quality was polluted at an average level and that it tended to be more seriously polluted in the dry months of March and June. The issues of total suspended solids, nutrients, organic matters, and ions pollution should be monitored regularly. Among the water quality parameters, only BOD and TN did not show a clear difference between observed seasons. On the other hand, the seasonal variation of the remaining parameters has significantly influenced the pollution concentration of rivers.

Distribution of Land Use Patterns
Satellite imagery was used to analyze different land use patterns by the available knowledge about the study area. The spatial distribution of land use patterns in 2019 is shown in Figure 2. Agricultural cultivation was considered the main activity of Tien Giang Province at the time of the study. As can be seen, most of the area used for rice cultivation accounted for about 97,746.6 ha (approximately 40.6% of the total area), which was distributed mainly in the north and northwest of the province. This was consistent with the soil characteristics in the area; alluvial soil occupied most of the area in Cai Be, Cai Lay, Chau Thanh districts. Meanwhile, perennial crop cultivation was often concentrated along the Tien River basin from west to east with an area of about 8.9%, equivalent to 21,414.2 ha. Agricultural land can be polluted by transporting fine residues, fertilizers, pesticides, etc. [54]. In addition, in the coastal area land (east and southeast of the province) with reflectance factor used for aquaculture accounts for 21,209.8 ha (about 8.8% of the total area, distributed mainly in Tan Phu Dong and Go Cong Dong District). According to the socioeconomic report of the People's Committee of Tan Phu Dong District in 2019, this land was mainly used for shrimp farming and rice-shrimp cultivation. Residential land was the second highest with 66,151.5 ha (27.5%); residential areas were mainly concentrated around cities in the region and contiguous areas of Long An Province. Residential land indicated that this was an adverse factor for water quality due to the lack of wastewater treatment systems and impervious surfaces comprising the majority of the area [13,55]. The total area for forest land and water bodies was about 20,934.2 ha and 13,044.9 ha, respectively, accounting for 8.7% and 5.4% of the area. Forests were concentrated mainly in the Tan Phu Dong District. Forest land plays important role in the hydrological process by water permeation and reducing both surface runoff and soil erosion [13,56,57].  Figure 3 shows that 34 sampling sites in the study area were divided into nine clusters with Euclidean distances less than 2. In addition, the mean values of water quality parameters at each cluster are presented in Table 2 for comparing and evaluating water quality between clusters. The results from Figure 3 indicate that the Euclidean distance between Cluster 1 and Cluster 2, compared to the other clusters, also showed a significant difference in water quality between the locations of Cluster 1 and Cluster 2. This can be seen at the positions S29, S32, and S34, which are located near the sea and affected by the tidal regime and great seawater encroachment. In addition, these locations belong to the large concentration of aquaculture of the province. The results of the water quality assessment show that the locations in Cluster 1 and Cluster 2 had water quality considered the most polluted, in-Distance cluding pollution of organic matters, nutrients, microorganisms, ions). Electrical conductivity and Cl − were major parameters distinguishing the water quality of Cluster 1 and Cluster 2 from the other clusters. Specifically, values of EC and Cl − at the positions of Cluster 1 (EC: 1394.5 μS cm −1 , Cl − : 430.4 mg L −1 ) and Cluster 2 (EC: 2915.0 μS cm −1 , Cl − : 2332.1 mg L −1 ) ( Table 2) were much higher than those in the other positions. The average values per year of EC and Cl − at the sampling sites were calculated from four observations before the mean values of parameters were calculated in the cluster; therefore, the seasonal variation of these parameters has also partly affected the results of water quality classification in CA, especially for locations located near the sea such as in Cluster 1 and Cluster 2. Moreover, the difference of EC was not the exception for the possibility that there was a difference in the geological composition, different humidity of the water body [58], the types of mud and clay are abundant of soluble ions. As can be inferred from Table  2, the concentration of nutrients derived from nitrogen is an important factor in eutrophication occurring in coastal waters, while phosphorus has been found in freshwater bodies. This result has also been reported in previous research that nitrogen is the primary cause of eutrophication in many coastal ecosystems [59]. Moreover, coliform densities in the Cluster 1 and Cluster 2 were found at high density, with values of 2170 MPN 100 mL −1 , 3366.7 MPN 100 mL −1 , respectively.

Clustering Surface Water Quality in Tien Giang Province
Meanwhile, Cluster 3 covered only one site (S13), where there was a major surface cover of the residential and rice cultivation with a high concentration of organic substances and nutrients derived from nitrogen. For Cluster 4 and Cluster 5, positions S17, S21, S22, S25, and S27, respectively, were included. Both Cluster 4 and Cluster 5 were identified on perennial farmland and surrounded by residential areas, where fewer nitrogen fertilizers were applied. Therefore, in contrast to Cluster 3, water quality characteristics in these clusters were determined to have a lower concentration of nitrogen-containing wastes, mainly phosphorus-derived nutrients. This was consistent with the fact that in the Vietnamese Mekong Delta, compared with rice cultivation, the amount of fertilizer used for fruit/perennial crop cultivation was assessed at a higher level, and farmers overuse phosphorus but lack nitrogen and potassium. In addition, the tradition of perennial crop cultivation in the highlands of the Mekong Delta, where rainwater can easily wash away many nutrients. Therefore, weather problems are also a part of the reason for the water quality classification results of Clusters 4 and 5. Moreover, the locations in Clusters 4 and 5 were located far away from the main river system (Tien River), thus limiting drainage capacity in the rainy season.
Cluster 6 consisted of the sites S18 and S33, which were located on the edge or perimeter of the study area and were mainly surrounded by rice farming and residential area. Therefore, organic matter pollution has also been recorded at these sites. Cluster 7 included the positions S4, S11, S12, S14, S15, and S26, while Cluster 8 included the positions S7, S8, S9, and S23. Clusters 7 and 8 were nearly the same but differ in the pollution levels, and both were mainly represented by residential land and agricultural (rice and perennial) land. Finally, the remaining sites were grouped together into one cluster (Cluster 9), which was a cluster representing all different land uses. However, the water quality in this cluster had a relatively low level of pollution. These locations were located near and along the Tien River, which accounted for the high dilution and self-cleaning ability of the water source.

Identification of Critical Water Parameters and Potential Sources of Water Variation
PCA analysis was performed using 16 water quality parameters at 34 monitoring locations with the values of the average of four sampling periods; the analytical results are shown in Table 3. From Table 3, it can be seen that the values of EC, COD, BOD, N-NH4 + , and coliform negatively correlated with the component PC1, which has eigenvalue coefficient and the percentage of variance of 4.63 and 29.0%, respectively. These parameters may be closely related to the runoff from fertilized farmland and domestic wastes. Component PC2 had eigenvalues of 2.59, and the percentage of variance was 16.2%, which negatively correlated with the surface water parameters such as EC, DO, Cl − and positively correlated with TSS, N-NO2 − , and TP. This PC2 reflected the effects of seasonal variation, saline intrusion by electrical conductivity, and Cl − in water, which was the parameter directly related to total dissolved salts [61,62]. Meanwhile, TSS may be related to land use conversion, weather conditions (rain) and soil erosion, and sedimentation of basins [63]. Furthermore, the evaluation results in Section 3.1 and the PCA analysis have shown that EC, TSS, and Cl − parameters can be considered the three most markedly seasonal variations of water bodies in Tien Giang Province. Due to the important role of water regimes in the rainy season (rainfall accounting for 90% of the annual rainfall and water flow from upstream) and dry season (water flow from upstream), which affected the concentration of dissolved salts in the water. In addition, at lower DO levels, the denitrification reaction can promote the formation of N-NO2 − . This may be the reason why N-NO2 − and DO showed a negative correlation in PC2.
pH was found to have a negative correlation with P-PO4 3− and TP in the PC3. In contrast, pH, and TN were negatively correlated in PC6 with the Eigenvalues values, and the percentage of variance was 1.02 and 6.4%, respectively. In parallel, PC5 also correlated with nutrient-related sources of pollution (N-NO3 − and TP). Furthermore, the correlation parameters with PC4 included pH, N-NH4 + , N-NO2 − , N-NO3 − , and SO4 2− ; these PCs can be predicted to be affected by agriculture and aquatic wastewater with 38.4% variation of the physicochemical composition explained. Aquatic wastewater contains many organic substances such as organic carbon, uric acid, and nutrients such as nitrogen and phosphorus. These compounds were derived from fertilizers, antibiotics, uneaten food, and metabolism converted to dissolved organic matter and nutrients [64]. Additionally, agriculture has previously been identified as one of the main sources of nitrogen (N) and phosphorus (P) contributing to water pollution; about 82% of the articles were related [57,65]. PCA analysis showed that all parameters could cause pollution of river water quality; however, the importance of the parameters at different land use patterns was not specified. Therefore, a method to calculate the weights of each parameter for different land use/land cover types has been performed in the next section of this study. This calculation confirmed the specificity and scientificity of the PCA analysis in the study. This was because each type of land use/land cover could induce a different pollutant concentration and parameter [66]. According to the land use pattern distribution map and CA analysis (Figures 2 and 3), the sites in Cluster 1 and Cluster 2 were distributed in the water bodies with a high aquaculture land use pattern, while Cluster 4 and Cluster 5 were located in the water bodies with high population and perennial crop cultivation. Similarly, Cluster 3 and Cluster 6 were suggested to be the water bodies mainly affected by residential area-rice cultivation; Clusters 7 and 8 were the residential area-cultivating rice and perennial crops. Finally, the positions in Cluster 9 were affected by all of the mentioned land use patterns. Therefore, the computation of entropy weights at each land use type is divided into five water body sectors with different impact factors.
In this study, TSS, N-NO2 − , and coliforms had the highest levels of pollution in the use of residential land-aquaculture; this result was similarly noted in the previous study in aquacultural areas in An Giang [67]. However, for residential land use-rice cultivation, the number, and importance of polluting parameters tended to increase, including TSS, N-NO2 − , TP, and coliform. This has been reported in previous research [66]. Furthermore, Table 4 shows the great influence of residential land-perennial crop cultivation on water quality by the parameters EC, DO, N-NH4 + , N-NO2 − , Cl − , and coliform. The high degree of EC influence on this type of land use was consistent with the other studies reporting that perennial and residential (urban and rural) land were the two main land use types causing the increase in this parameter [16], which is considered to be the most polluting places for water quality [16]. However, the rating was in contrast to the average values of the parameters at Clusters 4 and 5; This difference may be due to the fact that the land used for perennial crops was not large enough. It was shown that the popularity of the usage pattern also greatly influenced the variation in water quality. The study has also found that the mixed land use pattern of residential, rice cultivation, and perennial crops can significantly affect N-NO2 − , P-PO4 3− , and TP. This type of land use has been identified as a major contributor to the increase in nutrients in the water and has been reported in several studies [16,68]. The highest levels of coliform importance were found in the mixed land use pattern that combined all activities of residential, rice farming, perennial, aquaculture. N-NH4 + , N-NO2 − , and Cl − were the next most important parameters for quality variation in the study area with the complex socioeconomic activities. Temperature, pH, SO4 2− , N-NO3 − , and TN did not show a high degree of importance in the water bodies for different soil covers. This result was consistent with the water quality characteristics evaluated in Section 3.1, in which it has been suggested that temperature, pH, SO4 2− , N-NO3 − , and TN were still at safe levels. In addition, according to the previous study of Ly and Giao (2018) [69] and Giao (2020) [70], there was also no significant change in temperature and pH between water bodies. In general, the concentration of nitrogen (N-NH4 + , N-NO2 − , and TN) and orthophosphate (P-PO4 3− , and TP) tended to increase in the areas with a higher density of agricultural and residential land use, suggesting that these high levels may be the cause of wastewater inputs in residential areas and the fertilization of nitrogen and phosphate fertilizers in cultivated areas. In addition, coliform densities were found to be of high importance in most of the residential affected water bodies. Therefore, in order to reduce the impacts of economic activities on water quality, it is necessary to establish a suitable soil cover for each water body that both brings economic and environmental efficiency. Table 4. Entropy weights of water quality parameters for different land use types.

Land Use
Residential-Aquaculture

Conclusions
The results of the current study show that water quality in the study area was polluted by suspended solids, nutrients, organic matter, and ions. The seasonal change of water quality was indicated by the higher concentrations of water pollutants in the dry season than the wet season, except for TN and BOD. Five major land uses including residential-aquaculture, residential-rice cultivation, residential-perennials, residential riceperennial, and residential-rice-perennial-aquacultural were identified in Tien Giang Province, in which residential and rice cultivation were the majorities. CA analysis divided the sampling sites in the study areas into nine groups distributed on five land use patterns. Water quality at Cluster 1 and Cluster 2 (the locations with aquaculture land use patterns) was lower than that in the other land use patterns. The PCA analysis results show that five potential water-polluting sources resulting in water quality variation and all the current water quality monitoring variables were significant. The entropy weight calculation results show that TSS, N-NO2 − , and coliform were the important water parameters in the water bodies within the land use types of residential-aquaculture and residential-rice cultivation; EC, DO, N-NH4 + , N-NO2 − , Cl − , and coliform were significant in the residential-perennial land use type. Water quality in the land area of residential-rice cultivation-perennial crops was contaminated by N-NO2 − , P-PO4 3− , and coliform. Finally, the parameters of N-NH4 + , N-NO2 − , Cl − , and coliform were the main concerns for water quality in mixed land use patterns of residential-rice cultivation-perennial crops-aquaculture. The current results reflect the temporal and spatial variation of surface water quality in various land use patterns in Tien Giang Province, and this information could be very useful for future water quality management.