Abstract. Alluvial river estuaries consist largely of sand but are typically flanked by mudflats and salt marshes. The analogy with meandering rivers that are kept narrower than braided rivers by cohesive floodplain formation raises the question of how large-scale estuarine morphology and the late Holocene development of estuaries are affected by cohesive sediment. In this study we combine sand and mud transport processes and study their interaction effects on morphologically modelled estuaries on centennial to millennial timescales. The numerical modelling package Delft3D was applied in 2-DH starting from an idealised convergent estuary. The mixed sediment was modelled with an active layer and storage module with fluxes predicted by the Partheniades–Krone relations for mud and Engelund–Hansen for sand. The model was subjected to a range of idealised boundary conditions of tidal range, river discharge, waves and mud input. The model results show that mud is predominantly stored in mudflats on the side of the estuary. Marine mud supply only influences the mouth of the estuary, whereas fluvial mud is distributed along the whole estuary. Coastal waves stir up mud and remove the tendency to form muddy coastlines and the formation of mudflats in the downstream part of the estuary. Widening continues in estuaries with only sand, while mud supply leads to a narrower constant width and reduced channel and bar dynamics. This self-confinement eventually leads to a dynamic equilibrium in which lateral channel migration and mudflat expansion are balanced on average. However, for higher mud concentrations, higher discharge and low tidal amplitude, the estuary narrows and fills to become a tidal delta.
Abstract Alluvial fans develop their semi‐conical shape by quasi‐cyclic avulsions of their geomorphologically active sector from a fixed fan apex. On debris‐flow fans, these quasi‐cyclic avulsions are poorly understood, partly because physical scale experiments on the formation of fans have been limited largely to turbidite and fluvial fans and deltas. In this study, debris‐flow fans were experimentally created under constant extrinsic forcing, and autogenic sequences of backfilling, avulsion and channelization were observed. Backfilling, avulsion and channelization were gradual processes that required multiple successive debris‐flow events. Debris flows avulsed along preferential flow paths given by the balance between steepest descent and flow inertia. In the channelization phase, debris flows became progressively longer and narrower because momentum increasingly focused on the flow front as flow narrowed, resulting in longer run‐out and deeper channels. Backfilling commenced when debris flows reached their maximum possible length and channel depth, as defined by channel slope and debris‐flow volume and composition, after which they progressively shortened and widened until the entire channel was filled and avulsion was initiated. The terminus of deposition moved upstream because the frontal lobe deposits of previous debris flows created a low‐gradient zone forcing deposition. Consequently, the next debris flow was shorter which led to more in‐channel sedimentation, causing more overbank flow in the next debris flow and resulting in reduced momentum to the flow front and shorter runout. This topographic feedback is similar to the interaction between flow and mouth bars forcing backfilling and transitions from channelized to sheet flow in turbidite and fluvial fans and deltas. Debris‐flow avulsion cycles are governed by the same large‐scale topographic compensation that drives avulsion cycles on fluvial and turbidite fans, although the detailed processes are unique to debris‐flow fans. This novel result provides a basis for modelling of debris‐flow fans with applications in hazards and stratigraphy.
Abstract. In the list of challenges facing the world largest deltas, increased salt intrusion and its role in jeopardizing freshwater supply is often ranked very high. Yet, detailed process-based studies of salt intrusion at the whole delta scale are limited and the trends are regularly associated to global sea level rise. Here, using field measurements and a sophisticated 3D model that integrates the riverine, rural, estuarine, and coastal dynamics within one numerical domain, we study salt intrusion at the scale of the Mekong Delta in extensive detail. While many studies down-scale the salt intrusion problem to a topic within an estuary, we show that the continental shelf is an intrinsic component of the delta, and its physical processes, such as monsoon-driven ocean surge, directly influence salinity dynamics within the delta. Typical values of 20–40 cm surge over the continental shelf contribute to up to 10 km of further salt intrusion. The delta's estuarine system is also more sensitive than many other systems to upstream discharge variations. Furthermore, spring-neap variability plays a key role in salt intrusion in the delta. The estuarine variability from a stratified to a mixed system between neap and spring tides develops 3D processes such as estuarine circulation and tidal straining that become the main upstream salt transport mechanisms. The 3D nature of salinity dynamics, and the role of upstream and downstream processes, suggests that compromising on dimension or extent of the numerical domain, can limit the accuracy of predictions of salt intrusion in the delta. The study also showcases that riverbed incision in response to anthropogenic sediment starvation in the last two decades, has increased stratification, and activated or magnified 3D salt transport sub-processes that amplify upstream salt transport. With all the external forces on the delta namely climate change and altered hydrological regime by the upstream dams, due to deeper estuarine channels (driven by sand mining and upstream impoundments), the delta itself is far more vulnerable to even mild natural events. This exemplifies the fundamental importance of preserving the sediment budget and riverbed levels in protecting the world's deltas against salt intrusion.
The interaction of marine (tides and waves) and fluvial processes determines the sedimentary fill of coastal systems in the fluvial-to-marine (FTM) transition zone. Despite frequent recognition of tidal and wave influence in modern and ancient systems, our understanding of the relative importance of marine processes and their impact on mud deposition and reservoir architecture is limited. This study combined subsurface field observations and numerical simulations to investigate the relative importance of river flow, tides, waves, and mud input in governing the sedimentary fill in funnel-shaped basins along the FTM transition. Model simulations show a self-forming bar-built estuary with dynamic channels and sandy bars flanked by mud flats, which is in agreement with trends observed in nature. From three-dimensional virtual sedimentary successions, statistical tendencies for mud distribution and thickness were derived for the spectrum of marine and fluvial processes, and these values provide quantitative information on the net-to-gross ratio and mud architecture. The relative influence of marine and fluvial processes leads to a predictable facies organization and architecture, with muddier and more heterogeneous sediments toward the flanks. For the first time, our simulations allow the sedimentary fill in basins along the FTM transition to be related explicitly to hydrodynamic conditions, providing new insights into the morphosedimentary evolution of coastal systems, with implications for sequence stratigraphy.
Abstract. Vegetation enhances bank stability and sedimentation to such an extent that it can modify river patterns, but how these processes manifest themselves in full-scale estuarine settings is poorly understood. On the one hand, tidal flats accrete faster in the presence of vegetation, reducing the flood storage and ebb dominance over time. On the other hand flow-focusing effects of a tidal floodplain elevated by mud and vegetation could lead to channel concentration and incision. Here we study isolated and combined effects of mud and tidal marsh vegetation on estuary dimensions. A 2-D hydromorphodynamic estuary model was developed, which was coupled to a vegetation model and used to simulate 100 years of morphological development. Vegetation settlement, growth and mortality were determined by the hydromorphodynamics. Eco-engineering effects of vegetation on the physical system are here limited to hydraulic resistance, which affects erosion and sedimentation pattern through the flow field. We investigated how vegetation, combined with mud, affects the average elevation of tidal flats and controls the system-scale planform. Modelling with vegetation only results in a pattern with the largest vegetation extent in the mixed-energy zone of the estuary, which is generally shallower. Here vegetation can cover more than 50 % of the estuary width while it remains below 10 %–20 % in the outer, tide-dominated zone. This modelled distribution of vegetation along the estuary shows general agreement with trends in natural estuaries observed by aerial image analysis. Without mud, the modelled vegetation has a limited effect on morphology, again peaking in the mixed-energy zone. Numerical modelling with mud only shows that the presence of mud leads to stabilisation and accretion of the intertidal area and a slight infill of the mixed-energy zone. Combined modelling of mud and vegetation leads to mutual enhancement with mud causing new colonisation areas and vegetation stabilising the mud. This occurs in particular in a zone previously described as the bedload convergence zone. While vegetation focusses the flow into the channels such that mud sedimentation in intertidal side channels is prevented on a timescale of decades, the filling of intertidal area and the resulting reduction in tidal prism may cause the infilling of estuaries over centuries.
Abstract Predicting debris flow runout is of major importance for hazard mitigation. Apart from topography and volume, runout distance and area depends on debris flow composition and rheology, but how is poorly understood. We experimentally investigated effects of composition on debris flow runout, depositional mechanisms, and deposit geometry. The small‐scale experimental debris flows were largely similar to natural debris flows in terms of flow behavior, deposit morphology, grain size sorting, channel width‐depth ratio, and runout. Deposit geometry (lobe thickness and width) in our experimental debris flows is largely determined by composition, while the effects of initial conditions of topography (i.e., outflow plain slope and channel slope and width) and volume are negligible. We find a clear optimum in the relations of runout with coarse‐material fraction and clay fraction. Increasing coarse‐material concentration leads to larger runout. However, excess coarse material results in a large accumulation of coarse debris at the flow front and enhances diffusivity, increasing frontal friction and decreasing runout. Increasing clay content initially enhances runout, but too much clay leads to very viscous flows, reducing runout. Runout increases with channel slope and width, outflow plain slope, debris flow volume, and water fraction. These results imply that debris flow runout depends at least as much on composition as on topography. This study improves understanding of the effects of debris flow composition on runout and may aid future debris flow hazard assessments.
Cohesive floodplain sediment and vegetation are both thought to cause meandering river patterns. Our aims are to compare the isolated and combined effects of mud and vegetation on river planform and morphodynamics in the setting of intermediate-sized valley rivers. We use a numerical model for century-scale simulation of flow, sediment transport and morphology coupled with riparian vegetation settlement, growth and mortality as functions of species traits on which flow resistance depends. Mud fluxes were predicted by excess shear stress relations in combination with the active layer formulation. We found that valley-flooding water levels increase with vegetation density, causing a higher braiding intensity rather than meandering tendency. The shear stress during floods carves channels through the muddy floodplain surface. Higher mud concentration, on the other hand, increases floodplain aggradation, reduces the overbank flow frequency and ultimately causes formation of a single-thread channel. Vegetation causes mud to deposit closer to the river channel as a levee, showing that mud sedimentation and vegetation settling mutually enhance floodplain formation. However, mud and vegetation counteract in two ways. First, vegetation enhances floodplain accretion, which ultimately increases plant desiccation for high mud concentrations. Second, vegetation increases the tendency of periodic chute cutoffs in valleys. The chute cutoffs locally reset the landscape and create new windows of opportunity for the vegetation. Surprisingly, in systems with a high mud concentration this causes hysteretic loops of vegetation cover and delayed mud deposition. Ramifications for the interpretation of Palaeozoic fluvial facies are that even rootless vegetation, capturing cohesive mud closer to the river channel to form thicker floodplain on the point bar, can enhance the tendency to meander and, under high mud supply, form stable channels. However, meandering is more unlikely in narrower valley rivers with higher vegetation density.
Abstract Estuaries are dynamic landscapes with complex bar and channel patterns formed by interactions between tidal and fluvial currents. River discharge dampens the tidal wave, enhances the ebb flow, and supplies sediment to the estuary. However, it is largely unknown how river discharge influences overall estuary morphology. The objective of the current study is to quantify the control of river discharge on bar and channel dimensions and sediment transport throughout the estuary. To this end, a long‐term and large‐scale Delft3D‐2DH estuary model was designed with a suite of model runs undertaken where discharge systematically varied. Results show that tide‐dominated estuaries with significant river discharge can develop towards a dynamic equilibrium with a constant tidal prism through adjustment of channel dimensions to accommodate the supplied river discharge. It is essential to account for this morphodynamic adjustment when considering the transition from tide‐dominated estuaries to aggrading river‐dominated estuaries. After this transition, the estuary evolution depends on the discharge‐to‐width ratio. Tidal prism either decreases with higher river discharge as the tidal flow is dampened and the estuary aggrades, or increases when the estuary widens as it adjusts to the increase in total discharge. Additionally, results show that a higher river discharge increases the difference between the limit of flood‐dominant sediment transport and the limit of flow reversal, which has important implications for the preservation of the tidal signal in the stratigraphy. Estuary dimensions and channel patterns can be described as a function of river and tidal discharge. These findings indicate that the dynamic spatial component in numerical models is crucial in predicting trends in long‐term estuary morphology as well as in inverse predictions from stratigraphy.