Syn-sedimentary compaction or consolidation is an important process in deltaic environments because it affects both the local morphodynamics and hydrodynamics as well as the delta-scale accommodation space. However, the impact of syn-depositional compaction on the sediment distribution and the interdependency between different delta areas related to the sediment budget are not fully understood. This paper simulates syn-depositional compaction using improved 1D grain-size compaction formulations, integrated into hydrodynamic and morphodynamic modelling software Delft3D. The updated code is used to model sedimentation in mud-rich deltas under various compaction rate scenarios, which represents the maximum compaction rate potential of sediment that experiences the highest overburden stress in the delta. The simulated deltas are analysed by first classifying their plan-view area development into depositional elements: distributary channel, underfilled channel, delta plain, mouth bar, delta front and pro delta depositional elements. Then, sedimentation by mass, accommodation space and depositional segment metrics are calculated using the interpreted depositional elements. The results for zero compaction rate scenarios (0 mm year−1) show that limited space-varying and temporal-varying accommodation is available to deposit sediment in the delta plain depositional element. Therefore, the sedimentation mainly occurs in the mouth bar depositional element. For low-mid compaction rate scenarios (0.01–1 mm year−1), the additional syn-depositional accommodation space in the delta plain depositional element increases sedimentation in this area, limiting sedimentation in the mouth bar depositional element. For high compaction rate scenarios (>1 mm year−1), a further increase in the accommodation space in the delta plain depositional element leads to lateral sedimentation attributed to channel relocation, where the sedimentation mainly occurs in the mouth bar depositional element. This study shows that, although considered a gradual process, syn-sedimentary compaction does impact long-term delta evolution by influencing the distribution of sedimentation in the delta.

In natural deltaic settings, mixed hydrodynamic forcings and sediment properties are known to influence the preserved delta deposits. One process that has not received much attention yet is syn-sedimentary compaction of clastic sediment on millennial-scale delta evolution. To study how compaction interacts with delta morphodynamics and preserved sediment, a modelling approach is proposed. A 1D grain-size dependent compaction model was implemented into Delft3D-FLOW, which provides an opportunity to understand the underexplored connection between grain sizes supplied to the deltas and sediment compaction. The compaction model allows deposited sediment to decrease in volume due to the accumulation of newly deposited sediments above or the elapsed time. Differences in morphological trends are presented for scenarios defined by the composition of sediment supply (mud rich and sand rich) and the maximum allowed compaction rate in the model (0–10 mm year⁻¹). The resultant deposits are classified into sub-environments: delta top, delta front and pro delta. The delta top geometry (e.g. area increase, rugosity and aspect ratio), sediment distribution alongshore and across sub-environments, and delta top accommodation (e.g. volume reduction and average water depth) are compared. The modelling results show that compaction of the underlying delta front and pro delta deposits increases the average water depth at the delta top, driving morphological variability observed in the mud-rich and sand-rich deltas. The morphological changes are more prominent in the mud-rich deltas, which experience larger compaction-induced volume reduction for the same scenario. Moreover, higher compaction rates further increase the delta top accommodation, resulting in more deposition and evenly distributed sediment at the delta top. This leads to a less significant area increase and a wider delta top with a smoother coastline. The presented modelling results bridge the knowledge gap on the influence of syn-sedimentary compaction on long-term delta morphodynamics and preserved sediment. These findings can be applied to unravel the controlling processes in ancient delta deposits and predict the evolution of modern systems under changing climates.

This paper presents a new model, using existing consolidation theory, suitable for long-term morphodynamic simulations; we refer to the dynamic equilibrium consolidation (DECON) model. This model is applicable for muddy systems at small suspended particulate matter (SPM) concentrations, where the sedimentation rates are smaller than the consolidation rates and small fractions of sand can be accounted for. Thus, the model assumes quasi-equilibrium of the consolidating bed. It is derived from the full consolidation (Gibson) equation and is implemented in a mixed Lagrangian-Eulerian bed model guaranteeing stable and non-negative solutions, while numeric diffusion remains small. The erosion and deposition of sand and mud is accounted for, whereas internal mixing (e.g., bioturbation) is modeled through diffusion. The parameter settings for the new consolidation model (the hydraulic conductivity, consolidation coefficient, and strength) can be obtained from consolidation experiments in the laboratory. The model reproduces one-dimensional consolidation experiments and the qualitative behavior of erosion and deposition in a tidal flume. The DECON model was also applied to more natural conditions, simulating fine sediment dynamics on a schematized mud flat and in a schematized tidal basin under tide and wave forcing. The computational results of the mudflat simulations compared well with the simulations with the full Gibson equation. For the tidal basin simulations, DECON predicted the expected landward tidal transport of fine sediment during tide-dominated conditions, while the tidal basin withstood erosion during the more energetic wave-dominated periods. Computational times for the morphodynamic simulations of the tidal basin example without waves increased by a factor of 5 when consolidation was included. For the simulations with waves, this increase in computational times was only a factor of 2, as simulations with waves are always expensive. Applying a complete consolidation model would be prohibitive. The DECON model therefore serves as a useful tool to simulate fine-sediment dynamics in complex wave- and tide-dominated conditions, as well as the effects of seasonal variations.