Modeling long-term delta dynamics reveals persistent geometric river avulsion locations

Abstract

River deltas grow through repeated stacking of sedimentary lobes, the location and size of which are determined by channel avulsions (relatively sudden changes in river course). We use a model coupling fluvial and coastal processes to explore avulsion dynamics under a range of wave energies and sea-level-rise rates and find that the primary control on avulsion location and delta lobe size in our model is the critical superelevation ratio (SER), the amount of channel aggradation relative to the surrounding floodplain that is required to trigger an avulsion. The preferred avulsion location arises because of geometric constraints – a preferential avulsion node occurs at the break in floodplain slope that develops as the river progrades and/or sea level rises. This concavity develops in our model because the river profile aggrades and erodes via linear diffusion, whereas the diffusion of the floodplain topography is limited to episodic crevasse splays. These results are in contrast to recent modeling work, which was motivated by laboratory experiments and assumes a union between river channel and floodplain aggradation rates, and where avulsion nodes are driven by backwater hydrodynamics. The preferred avulsion length in our model scales well with laboratory, field, and model results without including hydrodynamic backwater effects. This work suggests an alternative mechanism to explain avulsion locations on deltas where floodplain topography aggrades and/or diffuses more slowly than the river channel profile, and it points to the need to elucidate river channel and floodplain connectivity over large space and time scales, and how the connectivity varies from one type of delta to another.