A Model-Based Approach for Improving Surface Water Quality Management in Aquaculture using MIKE 11: A Case of the Long Xuyen Quadangle, Mekong Delta, Vietnam

: This study utilized MIKE 11 to quantify the spatio-temporal dynamics of water quality parameters (Biochemical Oxygen Demand (BOD 5 ), Dissolved Oxygen (DO) and temperature) in the Long Xuyen Quadrangle area of the Vietnamese Mekong Delta. Calibrated for the year of 2019 and validated for the year of 2020, the developed model showed a signiﬁcant agreement between the observed and simulated values of water quality parameters. Locations near to cage culture areas exhibited higher BOD 5 values than sites close to pond/lagoon culture areas due to the effects of numerous point sources of pollution, including upstream wastewater and out-ﬂuxes from residential and tourism activities in the surrounding areas, all of which had a direct impact on the quality of the surface water used for aquaculture. Moreover, as aquacultural efﬂuents have intensiﬁed and dispersed over time, water quality in the surrounding water bodies has degraded. The ﬁndings suggest that the effective planning, assessment and management of rapidly expanding aquaculture sites should be improved, including more rigorous water quality monitoring, to ensure the long-term sustainable expansion and development of the aquacultural sector in the Long Xuyen Quadrangle in particular, and the Vietnamese Mekong Delta as a whole.


Introduction
Globally, demand for freshwater resources continues to increase. Freshwater sources are increasingly required to meet growing domestic, agriculture, aquaculture and industrial uses, while at the same time they suffer from increased pollution and natural and anthropogenic interventions and changes to the environment driven by strong population and economic growth [1][2][3][4]. In Southeast Asia, poorly managed aquacultural and agricultural activities are one of the most prominent sources of water pollution [2,[5][6][7][8][9]. Here, countries find it challenging to manage surface water quality due to both point and non-point sources of pollutants, as well as the increased widespread use of chemicals and drugs in the context of limited land and water resources and more intensive culture methods [10][11][12]. Due of the cumulative and synergistic impacts on water resources, water Against this background, our study aims to apply a hydrological simulation to investigate the spatio-temporal dynamics of key water quality parameters using MIKE 11. The hydrodynamic module, the structural operation module, the advection/dispersion module and the EcoLab module are all inclusive of the MIKE 11 model system. The outcomes of this study will firstly help decision makers to understand the current situation, and secondly help the design of future management options for improved water resource management.

Materials and Methods
The LXQ covers a large portion of the An Giang and Kien Giang provinces, and a small portion of Can Tho city, with a total area of approximately 0.5 million hectares. It is bordered to the north by the Bassac River and the Vietnamese-Cambodian border, to the south by the Cai San Canal, and to the west by the West Sea of Vietnam (Gulf of Thailand) [23]. Like the VMD as a whole, the topography of the LXQ, is relatively low and flat, with ground elevations of 0-1.0 m above mean sea level accounting for over 80% of the area [35,36]. The tropical monsoon climate of LXQ, has two primary seasons, dry and wet, and

Materials and Methods
The LXQ covers a large portion of the An Giang and Kien Giang provinces, and a small portion of Can Tho city, with a total area of approximately 0.5 million hectares. It is bordered to the north by the Bassac River and the Vietnamese-Cambodian border, to the south by the Cai San Canal, and to the west by the West Sea of Vietnam (Gulf of Thailand) [23]. Like the VMD as a whole, the topography of the LXQ, is relatively low and flat, with ground elevations of 0-1.0 m above mean sea level accounting for over 80% of the area [35,36]. The tropical monsoon climate of LXQ, has two primary seasons, dry and wet, and it is hot and humid all year round. As a result, LXQ's average annual temperature, rainfall and humidity are approximately 27 • C, 1200 mm and 80%, respectively [36].
The LXQ is significantly affected by its geographical location, its monsoon climate and the upstream river network and the tidal regime [35,37]. Within the LXQ, the yearly average flow of the river systems is around 14,000 m 3 ·s −1 . However, this can reach up to 24,000 m 3 ·s −1 during the rainy season, and recede to only 5000 m 3 ·s −1 during the dry season [36]. During the rainy season, annual floods inundate roughly 70% of the total LXQ area with water levels of 1.0-2.5 m for 3-5 months a year. This episodic flooding has both Water 2022, 14, 412 4 of 21 positive and negative effects to the socio-economic signature of this region. On the positive side, floodwaters bring large volumes of water for agriculture, aquaculture, domestic uses and industrial operations, and provide the region with nutrient-rich sediments, as well as helping to wash out pollutants and salinity from the soil. Flood disasters, on the other hand, severely damage infrastructure, interrupt both community and livelihood activities and jeopardize agricultural and fishery production [38,39]. As a result, both the provincial and national governments have made significant investments in local infrastructure to regulate water levels and protect the region through a series of full-dyke and semi-dyke systems [36]. However, this has resulted in changes to the hydrometeorological regime and fluxes, and lessened the potential for water pollutant dispersion.
To simulate water quality in the complex and dense river network of LXQ, this study utilized the hydrodynamics and EcoLab modules, which are the foundation of the onedimensional (1-D) MIKE 11 model.
Various hydro-meteorological data used as input for the modeling were kindly provided from a variety of sources, including the Southern Institute of Water Resources Research (SIWRR), the Southern Region Hydro-Meteorological Centre (SRHMC), the Department of Natural Resources and Environment (DoNRE), and the Department of Agriculture and Rural Development (DARD) ( Table 1). Additional data regarding the population, wastewater discharge, pollution load of BOD 5 and the current land use map were collected from residential, industrial, aquacultural and agricultural areas to estimate the pollution load discharges. Within the study area, hourly observations of the discharge were conducted at two stations at Chau Doc and Vam Nao, and these were used as the upstream boundary conditions, and hourly observations of the water level were conducted at two stations at Rach Gia and Can Tho, and these were used as the downstream boundary conditions. The hourly discharge and water level at Long Xuyen were obtained for calibration and verification for the years 2019 and 2020 for the period from 00:00 a.m. on 1st January to 23:59 p.m. on 31st May. Average monthly water quality data, such as BOD 5 , DO and temperature, were also obtained at eight stations from January to May in 2019 and 2020 for calibration and validation, respectively.
The total estimated fishery production for the year 2019 was 532.6 thousand tons, with pangasius production accounting for 412 tons in the study area. In 2020, due to the impact of the COVID-19 pandemic, the export of pangasius encountered many difficulties with the selling price becoming low, so farmers and businesses cut back on feed to prolong farming time, waiting for the price to increase, leading to a decrease in the harvest, with 211 thousand tons produced from the month of January to May 2020 [40]. Fish are mostly cultured in ponds and lagoons on both banks of the Bassac River, as well as along the tributaries Bay Tre, Xa Doi Canal, Cai Sao Canal, Don Dong Canal and Moi Canal ( Figure 1 and Table A1).

Estimation of Pollution Load
Previous studies have estimated pollution load in various Vietnamese river networks, and found the ratio of BOD 5 /Chemical Oxygen Demand (COD) to be around 0.65 [41,42]. This suggests that the majority of organic pollutants are soluble and easily decomposable. However, because of a lack of observed data in the study area, the pollutant load was calculated using the BOD 5 concentration as an indicator. Pollution load was computed using the discharge estimate, using Equation (1) as shown below: where Q M denotes the pollution load from the Mekong River (tons·year −1 ), Q denotes the discharge (m 3 ·s −1 ) and C i is the concentration of parameter i (mg·L −1 ). For domestic sources, pollution load was calculated based on the population statistics in the study area. The pollution emission coefficient per capita was calculated using Equation (2), as shown below: where Q D is pollution load from the population (tons·year −1 ), P is the population of the area (persons) and Q i is the domestic waste load of parameter i (kg·person −1 ·year −1 ) ( Table 2). For industrial areas, the pollution load was estimated by multiplying the industrial discharge (Q) by the pollutant emission coefficient of the industrial type, using Equation (3) as shown below: where Q I is the pollution load from industries, V j is the volume of annual wastewater discharged from industry j (m 3 ·year −1 ), C i,j is the concentration of substance i in the wastewater of industry j (mg·L −1 ) and n is the number of industries in the region. The pollution load from livestock production activities was calculated using the total annual livestock herd and the unit of discharge load for livestock and poultry, using Equation (4) as shown below: where Q L is the pollution load from livestock (tons·year −1 ), n is the number of livestock and poultry (unit) and Q i is the load of parameter i (kg·unit −1 ·year −1 ). The pollution load arising from aquaculture sources was calculated based on the aquacultural area and the coefficient of each waste generation for each different form of aquaculture, using Equation (5) as shown below: where Q A is the pollution load from aquacultural activities (tons·year −1 ), Q i is the load of the pollution source (kg·ha −1 ·day −1 ), S is the area of land used for farming (ha), and t is the time of farming in the year (day). Table 2 shows the calculation of the waste load generated on the basis of system emissions according to UNEP (1984) [43], San Diego-McGlone (2000) [41]. This estimation method has been successfully applied to many studies [44,45].

Model Setup
MIKE is a suite of software applications developed by the Danish Hydraulic Institute (DHI), consisting of different models (MIKE 11, MIKE 21, MIKE 3, MIKE SHE, Mouse and MIKE Basin), to accurately analyze, model and simulate rivers, lakes, estuaries and coastal environments. MIKE 11 includes the hydrodynamic module (HD) and EcoLab. The HD module permits the simulation of water levels, discharge and the discharge of wastewaters, while the EcoLab module (WQ module) describes how pollutants travel and disperse along rivers or channels over time. The HD module, applied on open-channel flows, solves finite differences from Saint-Venant equations consisting of the mass conservation and fluid momentum conservation, based on the following assumptions: (i) the flow is a dimension, with depth and velocity varying in the longitudinal direction of the channel; (ii) the bottom slope is small, and scour and deposition are negligible and the channel bed is fixed; (iii) flow everywhere is parallel to the bottom (i.e., wavelengths are large compared with water depths); (iv) the flow is sub-critical; (v) the water is incompressible and homogeneous, i.e., without significant variation in density; and (vi) the lateral inflow does not affect velocity in the channel.
The EcoLab module is based on the conservation of mass, which is a basic principle of the water quality model. It involves performing a mass balance for a defined control volume over a specified period of time. The EcoLab module is coupled to the AD module; while the EcoLab module deals with the transforming processes of compounds in the river, the AD module is used to simulate the simultaneous transport process. The AD equation is based on the following main assumptions: (i) the considered substance is completely mixed over the cross-section, (ii) the substance is conservative or subject to a first-order reaction (linear decay); and (iii) Fick's diffusion law is applied. Fick's law assumes that the mass flux is proportional to the gradient of the mean concentration and that the flux is in the direction of decreasing concentration. Dependent on the nature of the water quality problem under consideration, the model can be adjusted to different levels of detail. The complexity of the model ranges from the most simple version, which includes only BOD 5 and DO, through the introduction of sediment/water interactions and the inclusion of inorganic nitrogen (ammonia and nitrate), to the most complex level, where the BOD 5 is divided into three forms: dissolved, suspended and deposited. The MIKE 11 model was selected due to its robustness to represent a complex system, its flexibility to include possible future changes and its ability to build different plausible scenarios at different spatial scales. This study used a modified version of a MIKE 11 model to simulate hydrological dynamics over the entire LXQ area in consideration of different components, such as dykes, drainage, sluice gate operations, tidal influences and flood waters [23,46].
The following diagram depicts the modeling process used to replicate the area hydrodynamics and water quality ( Figure 2). In the MIKE 11 model, the hydrodynamic (HD), advection-dispersion (AD) and the ecological (EcoLab) modules are the three basic components.
The HD module is built upon the Saint-Venant Abbott's continuity (Equation (6)) and momentum (Equation (7)) equations. The LXQ stations in both the upstream and downstream sections are strongly affected by a diurnal tide in the West Sea and a semilunar diurnal tide regime in the East Sea. In the East Sea, the maximum tidal range is quite high (3.0-3.5 m) and the average tide range is approximately 2.5 m, while the tidal range in the West Sea is approximately 1 m [47]. The high average tidal amplitude in the East Sea varies from 2.2 m to 3.8 m, while the low average tidal amplitude in the West Sea is around 0.5 m [48]. As a result, the high tide period was maintained for a short duration, but the low tide period was maintained for a longer duration. The model includes 1581 river or canal segments, 1193 water storage sluice gate structures, over 7500 simulated water level nodes and 4700 simulated discharge nodes (flow). The majority of the river and canal network is seen in Figure 3. The HD module is built upon the Saint-Venant Abbott's continuity (Equation (6)) and momentum (Equation (7)) equations. The LXQ stations in both the upstream and downstream sections are strongly affected by a diurnal tide in the West Sea and a semi-lunar diurnal tide regime in the East Sea. In the East Sea, the maximum tidal range is quite high (3.0-3.5 m) and the average tide range is approximately 2.5 m, while the tidal range in the West Sea is approximately 1 m [47]. The high average tidal amplitude in the East Sea varies from 2.2 m to 3.8 m, while the low average tidal amplitude in the West Sea is around 0.5 m [48]. As a result, the high tide period was maintained for a short duration, but the low tide period was maintained for a longer duration. The model includes 1581 river or canal segments, 1193 water storage sluice gate structures, over 7500 simulated water level nodes and 4700 simulated discharge nodes (flow). The majority of the river and canal network is seen in Figure 3.
where Q is the discharge (m 3 . s −1 ), t is time (sec), A is the flow cross-sectional area (m 2 ), q represents lateral inflow per unit length (m 3 .s −1 .m −1 ), g is the gravitational acceleration (m.sec −2 ), h represents the height of the water level above sea level, n is the resistance coefficient (s.m −1/3 ), x is the direction, R is the hydraulic or resistance radius (m) and α is the momentum distribution coefficient (e).

∂Q ∂x
where Q is the discharge (m 3 ·s −1 ), t is time (sec), A is the flow cross-sectional area (m 2 ), q represents lateral inflow per unit length (m 3 ·s −1 ·m −1 ), g is the gravitational acceleration (m·s −2 ), h represents the height of the water level above sea level, n is the resistance x is the direction, R is the hydraulic or resistance radius (m) and α is the momentum distribution coefficient (e). The AD module of the MIKE 11 model simulates transportation based on the onedimensional equation of mass conservation for dissolved or suspended material (Equation (8)). As a result, this module requires the HD module's outputs, such as discharge and water level, cross-section area and hydraulic radius.

∂AC ∂t
where C is concentration (mg·L −1 ), D is the dispersion coefficient (m 2 ·s −1 ), A is the crosssectional area (m 2 ), K is the linear decay coefficient, C 2 is the source or sink concentration, q is the lateral discharge (m 3 ·s −1 ·m −1 ), x is the space coordinate (m) and t is the time coordinate (sec).  The AD module of the MIKE 11 model simulates transportation based on the onedimensional equation of mass conservation for dissolved or suspended material (Equation (8)). As a result, this module requires the HD module's outputs, such as discharge and water level, cross-section area and hydraulic radius.
where C is concentration (mg·L −1 ), D is the dispersion coefficient (m 2 ·s −1 ), A is the crosssectional area (m 2 ), K is the linear decay coefficient, C2 is the source or sink concentration, q is the lateral discharge (m 3 ·s −1 .m −1 ), x is the space coordinate (m) and t is the time coordinate (sec). The EcoLab module, which is based on a traditional water quality model, can simulate six levels of natural processes, ranging from the simplest BOD5-COD relationship to complex water quality processes, such as nitrification, denitrification, sediment precipitation and resuspension and sediment oxidation reduction, etc. Based on data availability, this study considered DO, BOD5 and temperature as water quality indicators. The release of organic waste into rivers causes a drop in DO. As a result, the model explains the relationship between BOD5 and DO using Equation (9): where Kd is the liner decay coefficient, BODd is the BOD decay, θ is the temperature coefficient for BOD decay and T is the temperature. The above equation describes the degradation of dissolved organic materials. The biological activities of aquatic habitats are determined by the amount of DO in the water. The DO regime indicators reveal the degree of organic load and the intensity of the breakdown and mineralization events.

Calibration and Validation
Changes in the model's parameters, such as Manning's hydraulic roughness coefficient (n) for the HD module, the diffusion coefficient for the AD module, and other parameters in the EcoLab module, such as the degradation of the BOD coefficient (level 1), were used to calibrate and validate the model. Furthermore, both were carried out to assure dependable performance by trial and error until the computed data matched the  The EcoLab module, which is based on a traditional water quality model, can simulate six levels of natural processes, ranging from the simplest BOD 5 -COD relationship to complex water quality processes, such as nitrification, denitrification, sediment precipitation and resuspension and sediment oxidation reduction, etc. Based on data availability, this study considered DO, BOD 5 and temperature as water quality indicators. The release of organic waste into rivers causes a drop in DO. As a result, the model explains the relationship between BOD 5 and DO using Equation (9): where K d is the liner decay coefficient, BOD d is the BOD decay, θ is the temperature coefficient for BOD decay and T is the temperature. The above equation describes the degradation of dissolved organic materials. The biological activities of aquatic habitats are determined by the amount of DO in the water. The DO regime indicators reveal the degree of organic load and the intensity of the breakdown and mineralization events.

Calibration and Validation
Changes in the model's parameters, such as Manning's hydraulic roughness coefficient (n) for the HD module, the diffusion coefficient for the AD module, and other parameters in the EcoLab module, such as the degradation of the BOD coefficient (level 1), were used to calibrate and validate the model. Furthermore, both were carried out to assure dependable performance by trial and error until the computed data matched the observed data. Hourly water level data from the Long Xuyen, Rach Gia and Chau Doc stations, as well as hourly discharge data from the Vam Nao station, were used for calibration. The study used water level data from 2020 to validate the model. Because of frequent use in the scientific works, the Nash-Sutcliffe efficiency (NSE), correlation coefficients (R) and root mean square error (RMSE) were employed to verify the model's performance for calibration and validation. The correlation coefficient (R) measures how strongly two variables are related to each other.
Formulas to calculate these three parameters are shown in Equations (10)- (12), starting with the correlation coefficient (Equation (10)): The Nash-Sutcliffe efficiency (NSE) measure determines the magnitude of residual variation in comparison to recorded data variance: The root mean square error (RMSE) is a measure of how different two datasets are, comparing one predicted value to a known or observed value: where X i is the observed data at time i, Y i is the simulated data at time i, X is the mean value of the observed data X = 1 n ∑ n i=1 X i and Y is the mean value of the simulated data Y = 1 n ∑ n i=1 Y i . Following that, using the best output data from the HD module, the AD and Eco-Lab modules were calibrated and verified. The calibration step continued until the best modeling output was obtained, while the validation step was used to test the calibrated parameters of the model. The rating of the selected efficiency criteria for R, NSE and RMSE follows the past work of Moriasi (2007) [49].

Calibration and Validation Results of HD Modeling
During the calibration and validation procedure, the Manning's hydraulic roughness coefficient (n) was found to be 0.03 for the global value (and varied between 0.015 and 0.075 for the local values). The hydrographs of simulated and observed tidal amplitude and water levels in 2019 and 2020 at Long Xuyen station were plotted, and the results are shown in Figures 4 and A1, respectively, whilst, the accuracy of the hydrodynamic model is reported in Table 3 [50], which revealed that the tidal amplitude in the East Sea of Vietnam fluctuates approximately less than 1.5 m., and that the tidal amplitude in the river tends to decrease compared that of the coast, this is also consistent with the study of Gagliano (1968) [51]. In addition, the correlation between the observed and simulated water levels and the tidal amplitudes showed strong agreement, as seen in Figure 5. This suggests that the developed HD module has a good performance and can be used for further simulation activities.
Water 2022, 13, x FOR PEER REVIEW 10 of 22 East Sea of Vietnam fluctuates approximately less than 1.5 m., and that the tidal amplitude in the river tends to decrease compared that of the coast, this is also consistent with the study of Gagliano (1968) [51]. In addition, the correlation between the observed and simulated water levels and the tidal amplitudes showed strong agreement, as seen in Figure  5. This suggests that the developed HD module has a good performance and can be used for further simulation activities.

Calibration and Validation Results of Water Quality Modeling
The EcoLab module is linked to the advection-dispersion (AD) module, which describes both the transformation and transport processes of pollutants. Therefore, the relationship between BOD5 and DO in different conditions can be used to describe the fluvial water quality. In the calibration step, the dispersion coefficients were seen to vary between 50 and 700, and BOD decay between 0.1 and 1.5. Figure 6 illustrates that the BOD5 concentration difference between simulated and observed data varies by about 12 percent, 22 percent, 24 percent and 37 percent for first, second, third and fourth sites, respectively. The observed BOD5 data of the first site were highest, while the accuracy of the fourth site was lowest. The lack of statistics on pollution load at site 4 is due to its proximately to Vinh Te Canal, which forms the boundary between Vietnam and Cambodia.

Calibration and Validation Results of Water Quality Modeling
The EcoLab module is linked to the advection-dispersion (AD) module, which describes both the transformation and transport processes of pollutants. Therefore, the relationship between BOD 5 and DO in different conditions can be used to describe the fluvial water quality. In the calibration step, the dispersion coefficients were seen to vary between 50 and 700, and BOD decay between 0.1 and 1.5. Figure 6 illustrates that the BOD 5 concentration difference between simulated and observed data varies by about 12 percent, 22 percent, 24 percent and 37 percent for first, second, third and fourth sites, respectively. The observed BOD 5 data of the first site were highest, while the accuracy of the fourth site was lowest. The lack of statistics on pollution load at site 4 is due to its proximately to Vinh Te Canal, which forms the boundary between Vietnam and Cambodia.
The simulation findings of water quality metrics demonstrate that the trend for temporal variation for the pollutants were similar between the main river channels and the smaller rivers and canals. The validation results of BOD 5 concentration are shown in Figure 7.   The simulation findings of water quality metrics demonstrate that the trend for temporal variation for the pollutants were similar between the main river channels and the smaller rivers and canals. The validation results of BOD5 concentration are shown in Figure 7. For the first, second, third and fourth sites, the errors between the simulated and observed data were 30 percent, 62 percent, 32 percent and 44 percent, respectively. The accuracy of the first location was the highest, while the accuracy of the second locality was the lowest. Rapid riverbank and riverbed erosion and sedimentation in the research area For the first, second, third and fourth sites, the errors between the simulated and observed data were 30 percent, 62 percent, 32 percent and 44 percent, respectively. The accuracy of the first location was the highest, while the accuracy of the second locality was the lowest. Rapid riverbank and riverbed erosion and sedimentation in the research area has led to recent changes in the slope and cross-section of rivers and canals. Moreover, at site no. 2, as well as sites 3 and 4, located in the canal with low discharge volumes, the immediate pollution concentration is often strongly affected directly by the discharge source. In site no. 1, on the main river channel, a much higher discharge volume resulted in quicker pollution dilution, while pollutant concentrations fluctuated less in small canals. The water quality model performed DO concentrations are in good agreement with the observed DO concentrations in 2019 and 2020 ( Figure A2). For the first, second, third and fourth sites, the difference between simulated and observed DO values, respectively, varied by about 14 percent, 11 percent, 14 percent and 13 percent in 2019, and 7 percent, 13 percent, 12 percent and 10 percent in 2020.
An ANOVA was used to further analyze water quality data, i.e., BOD 5 , DO and temperature, at stations along the Bassac River, infield canals and aquacultural site canals for the years 2019 and 2020. The aquacultural regions had the highest average BOD 5 (12.88 mg/L), followed by the infield canals (11.55 mg/L) and the Bassac River (9.07 mg/L) in 2019. The difference between BOD 5 in aquaculture and BOD 5 in the Bassac River, on the other hand, was only detected at a 5% significance level in the study. In the aquacultural area, the typical DO and temperature were around 4 mg/L and 30 • C, respectively. At a 5% significance level, the aquacultural area's average DO (4.47 mg/L) and temperature (29.78 • C) in 2020 differed from the river's average DO (5.04 mg/L) and temperature (29.08 • C). The aquatic region had the lowest BOD 5 (13.65 mg/L), while the Bassac River (15.47 mg/L) and the infield canals (15.17 mg/L) displayed the highest. We found that water pollution, such as BOD 5 concentrations in the year of 2020, was lower than in the year of 2019. We suspect the reason is due to the impact of the COVID-19 pandemic, leading to decreased fish production and resulting in a decreased BOD 5 load.
However, according to QCVN 08-MT: 2015/BTNMT, column A2 of the National Technical Regulation on surface water quality, the surface water quality deterioration by aquaculture in 2019 and 2020 found in this study was far in exceedance of the desired water quality threshold. The DO concentration of water, for instance, was lower than the QCVN 08 threshold. Hence, when utilizing this water for domestic purposes, it is advised that people pretreat it carefully before use to ensure no long-term adverse health effects. This water quality model may be used by key stakeholders to forecast water quality in the aquacultural area so that appropriate in situ water treatment measures can be located and implemented. Furthermore, numerous studies have revealed transboundary environmental degradation caused by the Mekong River's flow through two major tributaries, the Mekong River and the Bassac River. The result highlights the need for surface water quality control checks in border areas. Figure 8 shows the spatial distribution of BOD 5 max, BOD 5 min and BOD 5 max-BOD 5 min. The parameters are highest in the northwest and decrease gradually towards the southeast. High BOD 5 levels can be found in both Chau Doc and Tri Ton. Chau Doc is an urban area bordering Cambodia with moderate tourism activities. The main sources of wastewater in Chau Doc are from urban areas and tourism (average BOD 5 max and BOD 5 min concentrations are 25 mg/L and 17 mg/L, respectively). The population density in Chau Doc and the larger city of Long Xuyen was 963 and 2368 people per square kilometer in the year 2020, respectively. The remaining urban areas, particularly in Kien Giang province, have vastly lower population densities, whilst the impacts of urbanization on fluvial water quality remain minimal in coastal areas. Water quality in and around Long Xuyen is heavily affected not only by production activities, but also by waste from urban residential areas. Tri Ton, a semi-mountainous area in An Giang, is a popular tourist destination that is also undergoing rapid urbanization. Here, moderate levels of water quality were recorded, with the main source of pollution load being from surrounding aquacultural activities. Here, averages of BOD 5 max and BOD 5 min were 16 mg/L and 12 mg/L, respectively. However, on the main river (Bassac River) and inland canals, BOD 5 max was 17 mg/L and 13 mg/L, respectively. Aquacultural practices in the area to the north of Bassac River are typically in the form of in-river cage culture farming, while to the south, aquacultural production is typically practiced in ponds and lagoons. Overall, our results showed that the areas on or directly adjacent to the main river channels showed more variability due to the influence of higher flow volumes and exhibit, therefore, a better self-cleaning capacity [2,3]. Furthermore, according to Minh et al. [2,42], high pollution levels in the Vinh Te Canal have a degrading impact on the wider study area.
in the year 2020, respectively. The remaining urban areas, particularly in Kien Giang province, have vastly lower population densities, whilst the impacts of urbanization on fluvial water quality remain minimal in coastal areas. Water quality in and around Long Xuyen is heavily affected not only by production activities, but also by waste from urban residential areas. Tri Ton, a semi-mountainous area in An Giang, is a popular tourist destination that is also undergoing rapid urbanization. Here, moderate levels of water quality were recorded, with the main source of pollution load being from surrounding aquacultural activities. Here, averages of BOD5 max and BOD5 min were 16 mg/L and 12 mg/L, respectively. However, on the main river (Bassac River) and inland canals, BOD5 max was 17 mg/L and 13 mg/L, respectively. Aquacultural practices in the area to the north of Bassac River are typically in the form of in-river cage culture farming, while to the south, aquacultural production is typically practiced in ponds and lagoons. Overall, our results showed that the areas on or directly adjacent to the main river channels showed more variability due to the influence of higher flow volumes and exhibit, therefore, a better selfcleaning capacity [2,3]. Furthermore, according to Minh et al. [2,42], high pollution levels in the Vinh Te Canal have a degrading impact on the wider study area. High DO levels were found in the main river (Bassac) and main canals (Vinh Te Canal), while DO5 max-DO5 min, on the other hand, tends to follow the northwest direction ( Figure 9). High DO levels were found in the main river (Bassac) and main canals (Vinh Te Canal), while DO 5 max-DO 5 min, on the other hand, tends to follow the northwest direction ( Figure 9). tively. Moreover, (a, b) also show the location of urban and tourism areas, and provincial boundaries. These arrows show the downtrends of BOD5 concentration (a, b) and the downtrends of BOD5 max-BOD5 min.
High DO levels were found in the main river (Bassac) and main canals (Vinh Te Canal), while DO5 max-DO5 min, on the other hand, tends to follow the northwest direction ( Figure 9).

Discussion
Since the Streeter and Phelps model [52], more than a hundred further water quality models have been developed. The surface water quality models have progressed through many major stages, from assessing a single factor of water quality to multiple factors of water quality, from steady-to non-steady-state hydraulic models, from a point source model to a coupling-model of point and nonpoint sources and from one-dimensional to two-dimensional and three-dimensional models [53,54].
MIKE 11 is an advanced model with water quality management capabilities [55]. According to several studies, MIKE 11 necessitates a large amount of data and, in the case of using a small time step to ensure the stability of the model, simulation step takes a long time [56]. However, MIKE 11 also allows the user to select a simpler model that is appropriate for the target and dataset available. The EcoLab module, in particular, provides six different levels of conceptual water quality modeling complexity.
Moreover, some challenges remain, as models typically require a large amount of accurate data, an absence of which may result in a significant difference between the reality and the simulation results. Further, models do not include pollutant mechanisms, and are ambiguous about contaminant migration, meaning that we cannot predict how contaminants will migrate. Combination models are becoming popular for water quality modeling due to the complexity of water quality issues. In the case of water quality modeling, these models include combination models, artificial intelligence models and system integration [55].
In this study, MIKE 11 provided good results for the hydrodynamics simulation be-

Discussion
Since the Streeter and Phelps model [52], more than a hundred further water quality models have been developed. The surface water quality models have progressed through many major stages, from assessing a single factor of water quality to multiple factors of water quality, from steady-to non-steady-state hydraulic models, from a point source model to a coupling-model of point and nonpoint sources and from one-dimensional to two-dimensional and three-dimensional models [53,54].
MIKE 11 is an advanced model with water quality management capabilities [55]. According to several studies, MIKE 11 necessitates a large amount of data and, in the case of using a small time step to ensure the stability of the model, simulation step takes a long time [56]. However, MIKE 11 also allows the user to select a simpler model that is appropriate for the target and dataset available. The EcoLab module, in particular, provides six different levels of conceptual water quality modeling complexity.
Moreover, some challenges remain, as models typically require a large amount of accurate data, an absence of which may result in a significant difference between the reality and the simulation results. Further, models do not include pollutant mechanisms, and are ambiguous about contaminant migration, meaning that we cannot predict how contaminants will migrate. Combination models are becoming popular for water quality modeling due to the complexity of water quality issues. In the case of water quality modeling, these models include combination models, artificial intelligence models and system integration [55].
In this study, MIKE 11 provided good results for the hydrodynamics simulation because the input data for the study area was of a high temporal resolution (hourly). However, the water quality monitoring data consisted of only monthly averages. As such, changes in water quality were not clearly visible. The use of an Artificial Neural Network (ANN) to overcome this limitation produced significantly better results than hydrodynamics and water quality modeling, as evidenced by goodness-of-fit indices. The ANN model, developed by McCulloch [57], consists of three layers, and the number of layer increases with the complexity. The ANN model has been used to predict some key water quality parameters in recent years, with results indicating that its accuracy is adequate for practical purposes.
Previous research also suggests that the ANN performs better in predicting water level, as well as nitrate and phosphate, compared to Support Vector Machine models [58,59]. Moreover, Rabindra et al. reveal that ANN does not require other physical parameters in the modeling process, which can reduce the complexities of modeling the system [60]. Most numerical models, on the other hand, have not been fully validated against field experimental data, which often requires large investments of time and capital to obtain. Data-based models have been used to replace the numerical model because water quality predictions can be made using only accumulated data. The ANN is one of many datadriven techniques that has been widely applied due to its efficiency in predicting and forecasting water quantity and quality variables in river systems [61].

Conclusions
The hydrodynamics module developed in this study was calibrated and validated to be in good agreement with the observed data. After successful flow simulation, the AD and EcoLab modules of MIKE 11 were used to simulate the surface water quality of the LXQ using discharge and water level data from existing the HD module in the dry season. The results of the BOD 5 transmission simulation meet the practical needs of the study area by assisting the appropriate management and planning of aquacultural areas in order to reduce pollution and the impact of harmful chemicals from other sources High BOD 5 pollution was seen in most pond and lagoon aquacultural sites during the dry season in both 2019 and 2020. This is because ponds and lagoons have frequently slower growth stages, are less well ventilated than cage freshwater rearing areas, with limited water exchange and self-cleaning. However, the areas close to cage culture farming areas were found to have higher BOD 5 than areas surrounding pond/lagoon culture. The reason is due to the influences from other point pollution sources, such as upstream wastewater urban and tourism areas. Therefore, before locating and promoting aquacultural production, it is necessary to test site water quality. Furthermore, aquacultural areas should be located and zoned far from urban and tourist areas. Overall, the management and effective sustainable use of water resources needs to be strengthened to ensure the overall suitability of aquacultural locations.
Additionally, considering the data scarcity, several automatic measurement stations are proposed to aid the understanding of the spatio-temporal dynamics of water quality parameters, as well as data for improved model calibration and verification. Furthermore, diligent monitoring of all important hydrogeochemical parameters, such as salinity, isotopes/tracer elements, etc., which help in deciphering the seawater-freshwater mixing mechanisms in this dynamic hydrological system, should be considered as a future course of work.   Appendix B Table A1. Location of aquaculture ponds and lagoons in the study area.