Fluvial sediment supply towards the coast has been the subject of extensive research. Important aspects relate to the impact of sediment retaining hydropower dams, potential delta progradation, coastal sediment supply and delta vulnerability to sea level rise. Once validated, process-based models provide a valuable tool to address these aspects and offer detailed information on sediment pathways, distribution and budget in specific systems. This study aims to advance the understanding of the sediment dynamics and sediment budget in the Mekong Delta system. We developed a process-based model (Delft3D FM) that allows for coupling 2D area grids to 1D network grids. The flexible mesh describes both wide river sections and channel irrigation and drainage networks present in the Mekong Delta. We calibrated the model against observed discharge, salinity, suspended sediment concentration (SSC) and sediment flux. The model was able to skillfully describe seasonal variations of SSC and hysteresis of SSC and water discharge caused by the Tonle Sap Lake induced flow patterns and seasonally varying bed sediment availability in the channels. Model results suggest that the Mekong River delivers an amount of sediment, towards the delta which is much lower than the common estimate of 160 Mt/year. About 23% of the modeled total sediment load at Kratie reaches the sea. Our modeling approach is a useful tool to assess sediment dynamics under strategic anthropogenic interventions or climate change scenarios.

Building high dykes is a common measure of coping with floods and plays an important role in agricultural management in the Vietnamese Mekong Delta. However, the construction of high dykes causes considerable changes in hydrodynamics of the Mekong River. This paper aims to assess the impact of the high-dyke system on water level fluctuations and tidal propagation in the Mekong River branches. We developed a coupled 1-D to 2-D unstructured grid using Delft3D Flexible Mesh software. The model domain covered the Mekong Delta extending to the East (South China Sea) and West (Gulf of Thailand) seas, while the scenarios included the presence of high dykes in the Long Xuyen Quadrangle (LXQ), the Plain of Reeds (PoR) and the Trans-Bassac regions. The model was calibrated for the year 2000 high-flow season. Results show that the inclusion of high dykes changes the percentages of seaward outflow through the different Mekong branches and slightly redistributes flow over the low-flow and high-flow seasons. The LXQ and PoR high dykes result in an increase in the daily mean water levels and a decrease in the tidal amplitudes in their adjacent river branches. Moreover, the different high-dyke systems not only have an influence on the hydrodynamics in their own branch, but also influence other branches due to the Vam Nao connecting channel. These conclusions also hold for the extreme flood scenarios of 1981 and 1991 that had larger peak flows but smaller flood volumes. Peak flood water levels in the Mekong Delta in 1981 and 1991 are comparable to the 2000 flood as peak floods decrease and elongate due to upstream flooding in Cambodia. Future studies will focus on sediment pathways and distribution as well as climate change impact assessment.

Bathymetric data plays a major role in obtaining accurate results in hydrodynamic modeling of rivers, estuaries, and coasts. Bathymetries are commonly generated by spatial interpolation methods of data on a model grid. Sparse and limited data will impact the quality of the interpolated bathymetry. This study proposes an efficient spatial interpolation framework for producing a channel bathymetry from sparse, cross-sectional data. The proposed approach consists of three steps: (1) anisotropic bed topography data locations transformed to an orthogonal and smooth grid coordinate system that is aligned with its riverbanks and thalweg; (2) sample data are linearly interpolated to generate river bathymetry; and (3) the generated river bathymetry is converted into its original coordinates. The proposed approach was validated with a high spatial resolution topography of the Tieu estuarine branch. In addition, the proposed approach is compared with other spatial interpolation methods such as ordinary kriging, inverse distance weighting, and kriging with external drift. The proposed approach gives a nearly unbiased topography and a strongly reduced RMSE compared with the other methods. In addition, it accurately reproduces the thalweg. The proposed approach appears to be efficiently applicable for regions with sparse cross-sections. Moreover, river topography generated by the proposed approach is smooth including important morphologic features, making it suitable for two- and three-dimensional hydrodynamic modeling.

Fluvial sediment is the major source for the formation and development of the Mekong Delta. This paper aims to analyse the dynamics of suspended sediment and to investigate the roles of different processes in order to explore flux pattern changes. We applied modelling on two scales, comprising a large-scale model (the whole delta) to consider the upstream characteristics, particularly the Tonle Sap Lake's flood regulation, and a smaller-scale model (tidal rivers and shelf) to understand the sediment processes on the subaqueous delta. A comprehensive comparison to in-situ measurements and remote sensing data demonstrated that the model is capable of qualitatively simulating sediment dynamics on the subaqueous delta. It estimates that the Mekong River supplied an amount of 41.5 mil tons from April 2014 to April 2015. A substantial amount of sediment delivered by the Mekong River is deposited in front of the river mouths in the high flow season and resuspended in the low flow season. A sensitivity analysis shows that waves, baroclinic effects and bed composition strongly influence suspended sediment distribution and transport on the shelf. Waves in particular play an essential role in sediment resuspension. The development of this model is an important step towards an operational model for scientific and engineering applications, since the model is capable of predicting tidal propagation and discharge distribution through the main branches, and in predicting the seasonal SSC and erosion/deposition patterns on the shelf, while it is forced by readily available inputs: discharge at Kratie (Cambodia), GFS winds, ERA40 reanalysis waves, and TPXO 8v1 HR tidal forcing.