The accurate simulation of wetting–drying processes in floodplains and coastal zones is a challenge for hydrodynamic modelling, especially for long time simulations. Indeed, dedicated numerical procedures are generally time-consuming, instabilities can occur at the wet/dry front, rapid transition of wet/dry interface and mass conservation are not always ensured.We present the extension of an existing wetting–drying algorithm in two space dimensions and its application to a real case. The wetting–drying algorithm is implemented in Second-generation Louvain-la-Neuve Ice-ocean Model (www.slim-ocean.be), a discontinuous Galerkin finite element model solving the shallow water equations in a fully implicit way. This algorithm consists in applying a threshold value of fluid depth for a thin layer and a blending parameter in order to guarantee positive values of the water depth, while preserving local mass conservation and the well balanced property at wet/dry interfaces.The technique is first validated against standard analytical test cases (Balzano 1, Balzano 3 and Thacker test cases) and is subsquently applied in a realistic domain, the Tonle Sap Lake in the Mekong River Basin, where the water level can vary by about 10 m between the dry and the wet season.

The Mekong river is one of the largest rivers in the world, which flows through six countries of Southeast Asia (China, Myanmar, Laos, Thailand, Cambodia and Vietnam). Its hydro-sedimentary regime is changing rapidly, as a consequence of a regional shift of land use (agriculture, road, etc.), damming, sand mining and climate changes, among others. This study assesses the behavior of particles transported in suspension in the Lower Mekong River (LMR), along approximately 1700 km from fluvial to estuarine environments. Suspended sediment properties were estimated, simultaneously with hydrodynamic conditions, during three field campaigns. In addition, further investigations were performed in the laboratory to assess the structures of particles (flocculated or not), their capacity to flocculate (and the impacts on siltation), under a wide range of sediment concentration (20–30,000 mg.L⁻¹). This study confirms that suspended sediment transported in the LMR are predominantly (75% by volume) flocculi (or freshly eroded soils aggregates), with median aggregated particle size in the range 10–20 μm and median settling velocity of the order of 0.01–0.1 mm s⁻¹. These flocculi are robust under the hydrodynamic conditions (turbulence and suspended sediment concentration – SSC) existing in the LMR. Laboratory investigations reveal the existence of a threshold sediment concentration (400 mg.L⁻¹), beyond which flocculation and sedimentation increase of orders of magnitudes. Thus, concentration that exceeds this threshold might promote the formation of so-called fluid mud layers. Because of the nonlinear response of flocculation and sedimentation with SSC and considering the ongoing changes at a regional scale in the LMR, higher occurrence of fluid mud layers in the fluvial upstream waterbodies might be anticipated, and a lower occurrence in estuaries and alongshore where the concentration decrease. The geomorphology could be impacted, with an over-siltation in dams and an exacerbated erosion of the muddy-mangrove coast.

The Mekong river, located in Southeast Asia, is one of the largest rivers in the world. It is facing serious problems related to sediment transport, e.g. the reduction of sediment volume supply to the delta and estuary. This study focuses on the physical properties of suspended particles in the upstream fluvial section (Luang Prabang, Laos) and in the region of tidal influence (Song Hau estuary, Vietnam). These sites are respectively the upper and lower limits of the Lower Mekong River (LMR), which are experiencing contrasted transport modes. The main objective of this study is to quantify the proportion of inert sand and of cohesive particles in the water through comprehensive measurements of particle size, settling velocity and flocculation in the water column. Two field campaigns were undertaken in the summer monsoon season. The results point out the predominance of flocculi, with sizes of 18 ± 5 μm in the fluvial environment and 15 ± 5 μm in the estuary and the similarities in patterns, mixture of sand (12 ± 5 %) and mud (88 ± 5 %). They also provide a statistically representative set of suspended particle populations in the upstream and downstream regions that could be used for improving numerical models.

Large rivers often present a river–lake–delta system, with a wide range of temporal and spatial scales of the flow due to the combined effects of human activities and various natural factors, e.g., river discharge, tides, climatic variability, droughts, floods. Numerical models that allow for simulating the flow in these river–lake–delta systems are essential to study them and predict their evolution under the impact of various forcings. This is because they provide information that cannot be easily measured with sufficient temporal and spatial detail. In this study, we combine one-dimensional sectional-averaged (1D) and two-dimensional depth-averaged (2D) models, in the framework of the finite element model SLIM, to simulate the flow in the Mahakam river–lake–delta system (Indonesia). The 1D model representing the Mahakam River and four tributaries is coupled to the 2D unstructured mesh model implemented on the Mahakam Delta, the adjacent Makassar Strait, and three lakes in the central part of the river catchment. Using observations of water elevation at five stations, the bottom friction for river and tributaries, lakes, delta, and adjacent coastal zone is calibrated. Next, the model is validated using another period of observations of water elevation, flow velocity, and water discharge at various stations. Several criteria are implemented to assess the quality of the simulations, and a good agreement between simulations and observations is achieved in both calibration and validation stages. Different aspects of the flow, i.e., the division of water at two bifurcations in the delta, the effects of the lakes on the flow in the lower part of the system, the area of tidal propagation, are also quantified and discussed.