*Article* **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**

**Huynh Vuong Thu Minh <sup>1</sup> , Van Pham Dang Tri <sup>2</sup> , Vu Ngoc Ut <sup>3</sup> , Ram Avtar <sup>4</sup> , Pankaj Kumar 5,\* , Trinh Trung Tri Dang <sup>6</sup> , Au Van Hoa <sup>3</sup> , Tran Van Ty <sup>7</sup> and Nigel K. Downes <sup>1</sup>**


**Abstract:** This study utilized MIKE 11 to quantify the spatio-temporal dynamics of water quality parameters (Biochemical Oxygen Demand (BOD<sup>5</sup> ), 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 significant agreement between the observed and simulated values of water quality parameters. Locations near to cage culture areas exhibited higher BOD<sup>5</sup> values than sites close to pond/lagoon culture areas due to the effects of numerous point sources of pollution, including upstream wastewater and out-fluxes 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 effluents have intensified and dispersed over time, water quality in the surrounding water bodies has degraded. The findings 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.

**Keywords:** EcoLab module; hydrodynamics modeling; surface water quality; one dimension; cage culture; pond/lagoon culture

### **1. 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–4]. In Southeast Asia, poorly managed aquacultural and agricultural activities are one of the most prominent sources of water pollution [2,5–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–12]. Due of the cumulative and synergistic impacts on water resources, water

**Citation:** Thu Minh, H.V.; Tri, V.P.D.; Ut, V.N.; Avtar, R.; Kumar, P.; Dang, T.T.T.; Hoa, A.V.; Ty, T.V.; Downes, N.K. 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. *Water* **2022**, *14*, 412. https:// doi.org/10.3390/w14030412

Academic Editor: Domenico Cicchella

Received: 19 November 2021 Accepted: 26 January 2022 Published: 29 January 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

quality management requires a fundamental understanding of the spatial and temporal variations in water characteristics, including the hydro-morphological, chemical and biological parameters [13]. Previously, several methods have been developed to predict, monitor and assess water quality. These include using hydrogeochemical analysis [14–17], the use of various quality indexes [2,18,19], using numerical modeling for scenario development [20,21] and using socio-hydrological approaches to assess the nexus between water and human well-being [22]. All these approaches seek to provide a better understanding of the drivers and interlinkages at work, and to help decision makers take evidence-based actions with regard to improved water resource management.

However, using traditional hydro-chemical analysis and statistical approaches for analyzing the dynamics of water quality have inherent limitations when considering the different environmental components in a holistic manner. Moreover, they tend to be resource (money and human power) intensive. As a result, numerical simulation models or tools, including environmental modeling, that are capable of detecting regional and temporal changes in current and future water quality or quantity parameters are currently gaining increased popularity among scientists and water management practitioners [23,24]. In addition, the application of numerical simulations can also save labor, time and money [24,25].

Various hydraulics models, such as the Hydrological Engineering Centre—River Analysis System (HEC-RAS), MIKE 11 and Vietnamese River System and Plain (VRSAP) have been previously applied in the VMD to assess the changes in the quantity and quality of water in rivers, as well as quantify the impact of land management practices on water quality [26–28]. These models are holistic in nature and attempt to take into account the all-important environmental processes [16,20,21]. Because of its robust nature, the MIKE 11 model has previously been extensively used to investigate water security issues, especially in Asia [29].

As a rapidly developing and lower riparian country, Vietnam is particularly susceptible to water resource changes [30,31]. Southeast Asia's longest river, the Mekong River, appears to have its source on the Tibetan plateau, and runs through China, Myanmar, Laos, Thailand and Cambodia before reaching Vietnam. At its end, the Vietnamese Mekong Delta (VMD) is one of the world's largest river deltas, with a dense network of rivers, canals and ditches, and covering over 4 million hectares, which is approximately 12 percent of Vietnam's natural land area [32]. The VMD is also the largest agricultural and aquacultural hub in Vietnam, accounting for 50% of rice, 65% of aquaculture and 70% of fruit production, as well as including 95% of exported rice and 60% of exported fish [33].

The Long Xuyen Quadrangle (LXQ), as shown in Figure 1, is the first area in Vietnam to collect and use water from the Mekong River via two main branches, the Bassac River and the Mekong River. However, in the last few decades, water pollution has been an increasing problem in the delta [1,2,8,34]. The principal sources of pollution are non-point sources, such as wastewater discharged without treatment from industrial zones and commercial activities along the banks of rivers or canals. Furthermore, within the LXQ, intensive freshwater farming activities in the provinces of An Giang, Kien Giang and Can Tho city affect directly the water quality. In spite of its high regional importance and socio-economic significance, very few studies have investigated the water resources of LXQ in an in-depth, comprehensive and holistic manner.

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.

**Figure 1.** Map of the Vietnamese Mekong Delta and the Long Xuyen Quadrangle (LXQ) with its upstream and downstream boundaries. Calibration and validation samples were collected in the provinces of An Giang, Kien Giang and Can Tho. Long Xuyen Canal = 1, Vinh Tre Canal = 2, Rac Gia-Long Xuyen Canal = 3 and Tam Ngan Canal = 4. **Figure 1.** Map of the Vietnamese Mekong Delta and the Long Xuyen Quadrangle (LXQ) with its upstream and downstream boundaries. Calibration and validation samples were collected in the provinces of An Giang, Kien Giang and Can Tho. Long Xuyen Canal = 1, Vinh Tre Canal = 2, Rac Gia-Long Xuyen Canal = 3 and Tam Ngan Canal = 4.

#### Against this background, our study aims to apply a hydrological simulation to **2. Materials and Methods**

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. **2. Materials and Methods**  The LXQ covers a large portion of the An Giang and Kien Giang provinces, and a 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].

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 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<sup>3</sup> ·s −1 . However, this can reach up to 24,000 m<sup>3</sup> ·s <sup>−</sup><sup>1</sup> during the rainy season, and recede to only 5000 m<sup>3</sup> ·s <sup>−</sup><sup>1</sup> 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

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<sup>5</sup> and the current land use map were collected from residential, industrial, aquacultural and agricultural areas to estimate the pollution load discharges.


**Table 1.** Summary of hydrology and water quality data collection and sources.

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 BOD5, 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).

### *2.1. Estimation of Pollution Load*

Previous studies have estimated pollution load in various Vietnamese river networks, and found the ratio of BOD5/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<sup>5</sup> concentration as an indicator. Pollution load was computed using the discharge estimate, using Equation (1) as shown below:

$$\mathbf{Q\_M} = \mathbf{Q} \times \mathbf{C\_i} \tag{1}$$

where Q<sup>M</sup> denotes the pollution load from the Mekong River (tons·year−<sup>1</sup> ), Q denotes the discharge (m<sup>3</sup> ·s −1 ) and C<sup>i</sup> 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:

$$\mathbf{Q\_D = P \times Q\_i} \tag{2}$$

where Q<sup>D</sup> is pollution load from the population (tons·year−<sup>1</sup> ), P is the population of the area (persons) and Q<sup>i</sup> is the domestic waste load of parameter i (kg·person−<sup>1</sup> ·year−<sup>1</sup> ) (Table 2).

**Table 2.** Pollution load estimation.


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:

$$\mathbf{Q\_{I}} = \sum\_{\mathbf{j=1}}^{n} \mathbf{V\_{j}} \times \mathbf{C\_{i,j}} \tag{3}$$

where Q<sup>I</sup> is the pollution load from industries, V<sup>j</sup> is the volume of annual wastewater discharged from industry j (m<sup>3</sup> ·year−<sup>1</sup> ), Ci,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:

$$\mathbf{Q}\_{\rm L} = \mathbf{n} \times \mathbf{Q}\_{\rm i} \tag{4}$$

where Q<sup>L</sup> is the pollution load from livestock (tons·year−<sup>1</sup> ), n is the number of livestock and poultry (unit) and Q*<sup>i</sup>* is the load of parameter i (kg·unit−<sup>1</sup> ·year−<sup>1</sup> ).

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:

$$\mathbf{Q}\_{\mathbf{A}} = \mathbf{Q}\_{\mathbf{i}} \times \mathbf{S} \times \mathbf{t} \tag{5}$$

where Q<sup>A</sup> is the pollution load from aquacultural activities (tons·year−<sup>1</sup> ), Q<sup>i</sup> is the load of the pollution source (kg·ha−<sup>1</sup> ·day−<sup>1</sup> ), 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].

### *2.2. 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<sup>5</sup> 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<sup>5</sup> 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. *Water* **2022**, *13*, x FOR PEER REVIEW 7 of 22

**Figure 2.** The hydrodynamics and water quality modeling procedure phases utilizing one-dimensional MIKE 11. **Figure 2.** The hydrodynamics and water quality modeling procedure phases utilizing onedimensional MIKE 11.

The HD module is built upon the Saint-Venant Abbott's continuity (Equation (6)) and

$$\frac{\partial \mathbf{Q}}{\partial \mathbf{x}} + \frac{\partial \mathbf{A}}{\partial \mathbf{t}} = \mathbf{q} \tag{6}$$

$$\frac{\partial \mathbf{Q}}{\partial \mathbf{t}} + \frac{\partial \left( \alpha \frac{\mathbf{Q}^2}{\mathbf{A}} \right)}{\partial \mathbf{x}} + \mathbf{g} \mathbf{A} \frac{\partial \mathbf{h}}{\partial \mathbf{x}} + \frac{n^2 \mathbf{g} \mathbf{Q} |\mathbf{Q}|}{A R^{4/3}} = \mathbf{0} \tag{7}$$

ies 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 netwhere Q is the discharge (m<sup>3</sup> ·s −1 ), t is time (sec), A is the flow cross-sectional area (m<sup>2</sup> ), q represents lateral inflow per unit length (m<sup>3</sup> ·s −1 ·m−<sup>1</sup> ), g is the gravitational acceleration (m·s −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 <sup>α</sup> is the momentum distribution coefficient (e).

work is seen in Figure 3. ∂Q <sup>∂</sup>x<sup>+</sup> ∂A <sup>∂</sup>t = q (6) 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.

$$\frac{\partial \mathbf{A} \mathbf{C}}{\partial \mathbf{t}} + \frac{\partial \mathbf{Q} \mathbf{C}}{\partial \mathbf{x}} - \frac{\partial}{\partial \mathbf{x}} \left( \mathbf{A} \mathbf{D} \frac{\partial \mathbf{A}}{\partial \mathbf{x}} \right) = -\mathbf{A} \mathbf{K} \mathbf{C} + \mathbf{C}\_2 \mathbf{q} \tag{8}$$

where Q is the discharge (m3. s−1), t is time (sec), A is the flow cross-sectional area (m2), q represents lateral inflow per unit length (m3.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 where C is concentration (mg·L −1 ), D is the dispersion coefficient (m<sup>2</sup> ·s −1 ), A is the crosssectional area (m<sup>2</sup> ), K is the linear decay coefficient, C<sup>2</sup> is the source or sink concentration, q is the lateral discharge (m<sup>3</sup> ·s −1 ·m−<sup>1</sup> ), x is the space coordinate (m) and t is the time coordinate (sec).

momentum distribution coefficient (e).

**Figure 3.** The river and canal network in the LXQ. Blue squares denote the HD (discharge in the upstream and water levels in the downstream) and EcoLab (BOD5, DO and temperature) boundaries. **Figure 3.** The river and canal network in the LXQ. Blue squares denote the HD (discharge in the upstream and water levels in the downstream) and EcoLab (BOD<sup>5</sup> , DO and temperature) boundaries.

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 + ∂QC <sup>∂</sup><sup>x</sup> <sup>−</sup> <sup>∂</sup> <sup>∂</sup><sup>x</sup> <sup>൬</sup>AD ∂A ∂x ൰ = -AKC + C2q (8) 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, BOD<sup>5</sup> 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<sup>5</sup> and DO using Equation (9):

$$\frac{\text{dBOD}\_{\text{d}}}{\text{d}\_{\text{l}}} = \text{K}\_{\text{d}} \cdot \text{BOD}\_{\text{d}} \cdot \theta^{(\text{T}-\text{20})} \tag{9}$$

q is the lateral discharge (m3·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 simuwhere K<sup>d</sup> is the liner decay coefficient, BOD<sup>d</sup> is the BOD decay, θ is the temperature coefficient for BOD decay and T is the temperature.

late 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 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.

#### of organic waste into rivers causes a drop in DO. As a result, the model explains the rela-*2.3. Calibration and Validation*

tionship between BOD5 and DO using Equation (9): dBODୢ d୲ = Kୢ · BODୢ · θ(ିଶ) (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. *2.3. 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.

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

Formulas to calculate these three parameters are shown in Equations (10)–(12), starting with the correlation coefficient (Equation (10)):

$$\mathbf{R} = \frac{\frac{\sum\_{i=1}^{n} \left(\mathbf{x}\_{i} - \bar{\mathbf{X}}\right) \left(\mathbf{y}\_{i} - \bar{\mathbf{Y}}\right)}{\left(\mathbf{n} - \mathbf{1}\right)}}{\sqrt{\frac{1}{\mathbf{n} - \mathbf{1}} \sum\_{i=1}^{n} \left(\mathbf{X}\_{i} - \bar{\mathbf{X}}\right)^{2}} \sqrt{\frac{1}{\mathbf{n} - \mathbf{1}} \sum\_{i=1}^{n} \left(\mathbf{Y}\_{i} - \bar{\mathbf{Y}}\right)^{2}}}\tag{10}$$

The Nash–Sutcliffe efficiency (NSE) measure determines the magnitude of residual variation in comparison to recorded data variance:

$$\text{NSE} = \frac{\sum\_{i=1}^{n} \left(\chi\_{\text{i}} - \bar{\boldsymbol{\chi}}\right)^{2} - \sum\_{i=1}^{n} \left(\chi\_{\text{i}} - \chi\_{\text{i}}\right)^{2}}{\sum\_{i=1}^{n} \left(\chi\_{\text{i}} - \bar{\boldsymbol{\chi}}\right)^{2}} \tag{11}$$

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:

$$\text{RMSE} = \sqrt{\frac{1}{n} \sum\_{i=1}^{n} (\chi\_i - \chi\_i)^2} \tag{12}$$

where X<sup>i</sup> is the observed data at time i, Y<sup>i</sup> is the simulated data at time i, X is the mean value of the observed data X = <sup>1</sup> <sup>n</sup> ∑ n <sup>i</sup>=<sup>1</sup> X<sup>i</sup> and Y is the mean value of the simulated data Y = <sup>1</sup> <sup>n</sup> ∑ n <sup>i</sup>=<sup>1</sup> Y<sup>i</sup> .

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].

### **3. Results**

### *3.1. 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. Overall, the model's NSE, R and RMSE values were in the range of 0.84–0.90, 0.94–0.96 and 0.15–0.22 for water level performance, respectively. This indicates the model's overall good performance. While the 2019 calibration performed better than the 2020 validation, the NSE value for 2019 is 0.89, compared to 0.85 in 2020, and the RMSE value for 2019 (0.17) is likewise lower than that of 2020 (0.21). During this time, the river morphology, river network, the number of dykes and the sluice systems all changed slightly. In 2019, the comparison of tidal amplitudes performed quite well, with 0.89, 0.73 and 0.04 for NSE, R and RMSE, respectively. The results for tidal amplitude tide verification in 2020 were 0.71, 0.79 and 0.05 for NSE, R and RMSE, respectively. The average tidal amplitude at the Long Xuyen station was 0.8, and varied from 0.52 m to 1 m in 2019. In 2020, the tidal average amplitude was 0.85, and varied from 0.49 m to 1 m. According to research by Phan (2019) [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

*Water* **2022**, *13*, x FOR PEER REVIEW 10 of 22

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. 5. This suggests that the developed HD module has a good performance and can be used for further simulation activities.

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

**Figure 4.** Time series of daily simulated and observed amplitudes in both 2019 (**a**) and 2020 (**b**) at Long Xuyen station.

Long Xuyen station. **Table 3.** Summary of correlation coefficient (R), Nash–Sutcliffe efficiency (NSE) and root mean square error (RMSE) of water levels in the dry season (Jan–May) in 2019 and 2020 at Long Xuyen station.

**Figure 4.** Time series of daily simulated and observed amplitudes in both 2019 (**a**) and 2020 (**b**) at


amplitude Jan–May 0.89 0.04 0.73 0.71 0.05 0.79

March 0.88 0.18 0.95 0.84 0.21 0.96 April 0.87 0.18 0.95 0.84 0.22 0.96 May 0.89 0.17 0.95 0.86 0.21 0.95 Jan–May 0.89 0.17 0.95 0.85 0.21 0.96

> **Figure 5.** The correlation plots between observed and simulated water level data at an hourly scale for the dry seasons in (**a**) 2019 and (**b**) 2020. The correlation plots between observed and simulated tidal amplitude at a daily scale for the dry seasons in (**c**) 2019 and (**d**) 2020. **Figure 5.** The correlation plots between observed and simulated water level data at an hourly scale for the dry seasons in (**a**) 2019 and (**b**) 2020. The correlation plots between observed and simulated tidal amplitude at a daily scale for the dry seasons in (**c**) 2019 and (**d**) 2020.

#### *3.2. Calibration and Validation Results of Water Quality Modeling 3.2. 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. 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<sup>5</sup> 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<sup>5</sup> 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<sup>5</sup> 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<sup>5</sup> concentration are shown in Figure 7.

12

18

Simulated Observed

QCVN08\_MT:2015 BTNMT

*Water* **2022**, *13*, x FOR PEER REVIEW 12 of 22

12

18

Simulated Observed

QCVN08\_MT:2015 BTNMT

Site 2

Site 1

**Figure 6.** Calibrated results of BOD5 concentration at the four sites within the study area for the year 2019. The National Technical Regulation on surface water quality (QCVN 08MT:2015) was approved by the Ministry of Environment and Natural Resources (MoNRE) in 2015. **Figure 6.** Calibrated results of BOD<sup>5</sup> concentration at the four sites within the study area for the year 2019. The National Technical Regulation on surface water quality (QCVN 08MT:2015) was approved by the Ministry of Environment and Natural Resources (MoNRE) in 2015. poral 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.

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 ca-

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

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,

An ANOVA was used to further analyze water quality data, i.e.*,* BOD5, 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 BOD5 (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 BOD5 in aquaculture and BOD5 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 BOD5 (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 BOD5 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

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

to decreased fish production and resulting in a decreased BOD5 load.

13 percent, 12 percent and 10 percent in 2020.

nals.

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., BOD5, 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<sup>5</sup> (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<sup>5</sup> in aquaculture and BOD<sup>5</sup> 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<sup>5</sup> (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<sup>5</sup> 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<sup>5</sup> 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<sup>5</sup> max, BOD<sup>5</sup> min and BOD<sup>5</sup> max–BOD<sup>5</sup> min. The parameters are highest in the northwest and decrease gradually towards the southeast. High BOD<sup>5</sup> 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<sup>5</sup> max and BOD<sup>5</sup> 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<sup>5</sup> max and BOD<sup>5</sup> min were 16 mg/L and 12 mg/L, respectively. However, on the main river (Bassac River) and inland canals, BOD<sup>5</sup> 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. tural 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.

Mekong River and the Bassac River. The result highlights the need for surface water qual-

Figure 8 shows the spatial distribution of BOD5 max, BOD5 min and BOD5 max–BOD5 min. The parameters are highest in the northwest and decrease gradually towards the southeast. High BOD5 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 BOD5 max and BOD5 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 aquacul-

*Water* **2022**, *13*, x FOR PEER REVIEW 14 of 22

ity control checks in border areas.

**Figure 8.** BOD5 concentration in the dry season of 2020. Here, (**a**–**c**) represent the maximum, minimum and the difference between the maximum and minimum BOD5 concentration values, respectively. 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. **Figure 8.** BOD<sup>5</sup> concentration in the dry season of 2020. Here, (**a**–**c**) represent the maximum, minimum and the difference between the maximum and minimum BOD<sup>5</sup> concentration values, respectively. Moreover, (**a**, **b**) also show the location of urban and tourism areas, and provincial boundaries. These arrows show the downtrends of BOD<sup>5</sup> concentration (**a**, **b**) and the downtrends of BOD<sup>5</sup> max–BOD<sup>5</sup> 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). High DO levels were found in the main river (Bassac) and main canals (Vinh Te Canal), while DO<sup>5</sup> max–DO<sup>5</sup> min, on the other hand, tends to follow the northwest direction (Figure 9).

(**a**) (**b**)

(**c**)

**Legend**

< 5 5.1 – 10 11 – 15 16 – 20 21 – 25 26 – 30 31 – 40 41 – 50 51 – 60 > 60 mountains

BOD5 Max/Min (mg/L)

**Legend**

< 5 5.1 - 10 11 - 15 16 - 20 21 - 25 26 - 30 31 - 40 41 - 50 51 - 60 > 60 mountains

BOD5 Max/Min (mg/L)

Fish pond/cage location

**Figure 8.** BOD5 concentration in the dry season of 2020. Here, (**a**–**c**) represent the maximum, minimum and the difference between the maximum and minimum BOD5 concentration values, respectively. 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

Fish pond/cage location

mTourism area oUrban 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.** DO concentration in the dry season of 2020. Here, (**a**–**c**) represent the maximum, minimum and difference between the maximum and minimum DO concentration values, respectively. These arrows show the downward trend of DO concentration (9a and 9b) and the downward trend of DO max–DO min. **Figure 9.** DO concentration in the dry season of 2020. Here, (**a**–**c**) represent the maximum, minimum and difference between the maximum and minimum DO concentration values, respectively. These arrows show the downward trend of DO concentration (9a,b) and the downward trend of DO max–DO min.

#### **4. Discussion 4. Discussion**

integration [55].

practical purposes.

max–BOD5 min.

(Figure 9).

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]. 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. 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 con-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

taminants 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

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

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 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].

### **5. 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<sup>5</sup> 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<sup>5</sup> 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<sup>5</sup> 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.

**Author Contributions:** Conceptualization, V.N.U., V.P.D.T. and H.V.T.M.; methodology, V.N.U., V.P.D.T., T.T.T.D. and H.V.T.M.; software, T.T.T.D., A.V.H., V.P.D.T. and H.V.T.M.; writing—original draft preparation, V.N.U., V.P.D.T., R.A., P.K., T.V.T., T.T.T.D., N.K.D. and H.V.T.M.; writing—review

and editing, V.N.U., V.P.D.T., R.A., P.K., T.V.T., N.K.D. and H.V.T.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study is funded in part by the Can Tho University Improvement Project VN14-P6, supported by a Japanese ODA loan.

*Water* **2022**, *13*, x FOR PEER REVIEW 18 of 22

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

### **Appendix A**

**Figure A1.** *Cont*.

*Water* **2022**, *13*, x FOR PEER REVIEW 19 of 22

*Water* **2022**, *13*, x FOR PEER REVIEW 19 of 22

**Figure A1.** Comparison of simulated and observed water levels at the Long Xuyen station in the dry season 2019*.* **Figure A1.** Comparison of simulated and observed water levels at the Long Xuyen station in the dry season 2019. season 2019*.*

 **Figure A2.** Calibrated and verified results of DO concentration at the four locations within the study **Figure A2.** Calibrated and verified results of DO concentration at the four locations within the study area for the years 2019 and 2020.

area for the years 2019 and 2020.

### **Appendix B**

**Table A1.** Location of aquaculture ponds and lagoons in the study area.


### **References**

