The morphological development of fluvial and tidal systems is forecast more and more frequently by models in scientific and engineering studies for decision making regarding climate change mitigation, flood control, navigation and engineering works. However, many existing morphodynamic models predict unrealistically high channel incision, which is often dampened by increased gravity-driven sediment transport on side-slopes by up to two orders of magnitude too high. Here we show that such arbitrary calibrations dramatically bias sediment dynamics, channel patterns, and rate of morphological change. For five different models bracketing a range of scales and environments, we found that it is impossible to calibrate a model on both sediment transport magnitude and morphology. Consequently, present calibration practice may cause an order magnitude error in either morphology or morphological change. We show how model design can be optimized for different applications. We discuss the major implications for model interpretation and a critical knowledge gap.

Natural levees are common features in river, delta and tidal landscapes. They are elevated near‐channel morphological features that determine the connection between channel and floodbasin, and consequently affect long‐term evolution up to delta‐scales. Despite their relevance in shaping fluvial‐tidal systems, research on levees is sparse and often limited to fluvial or non‐tidal case studies. There is also a general lack of understanding of the role of vegetation in shaping these geomorphic units, and how levee morphology and dimensions vary in the transition from fluvial to coastal environments, where tides are increasingly important. Our goal is to unravel the effects of fluvial‐tidal boundary conditions, sediment supply and vegetation on levee characteristics and floodbasin evolution. These conditions were systematically explored by 60 large‐scale idealized morphodynamic simulations in Delft3D which self‐developed levees over the course of one century. We compared our results to a global levee dataset compilation of natural levee dimensions. We found that levee height is determined by the maximum water level, provided sufficient levee building sediments are available. Discharge fluctuations increased levee width and triggered more levee breaches, i.e. crevasses, that effectively filled the fluvio‐tidal floodbasin. The presence of wood‐type (sparse) vegetation further increased the number of crevasses in comparison with the non‐vegetated scenarios. Conversely, reed‐type (dense) vegetation strongly dampened tidal amplitude and reduced the accommodation space and sedimentation further into the floodbasin, resulting in narrower levees, no crevasses and limited floodbasin accretion. However, dense vegetation reduced tidal forces which allowed levee growth further downstream. Ultimately, the levees merged with the coastal barrier, eliminating the floodbasin tides entirely. Our results elucidate the mechanisms by which levee and crevasse formation, and vegetation may fill fluvio‐tidal wetlands and affect estuary evolution. This brings new insights for geological reconstructions as well as for the future management of deltas and estuaries under sea‐level rise. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd

Summary: Sediment‐stabilizing and ‐destabilizing organisms, i.e. microphytobenthos (biofilms) and macrozoobenthos (bioturbators), affect the erodibility of muddy sediments, potentially altering large‐scale estuarine morphology. Using a novel eco‐morphodynamic model of an idealized estuary, we investigate eco‐engineering effects of microphytobenthos and two macrozoobenthic bioturbators. Local mud erodibility is based on species pattern predicted through hydrodynamics, soil mud content, competition and grazing.Mud resuspension and export is enhanced under bioturbation and prevented under biostabilization through respective exposure and protection of the supra‐ and intertidal. Bioturbation decreases mud thickness and bed elevations, which increases net mud fluxes. Microphytobenthos reduces erosion, leading to a local mud increase of intertidal sediments.In multi‐species scenarios, an effective mud‐prone bioturbator strongly alters morphology, exceeding that of a more abundant sand‐prone moderate species, showing that morphological change depends on species traits as opposed to abundance. Altering their habitat, the effective mud‐prone bioturbator facilitates expansion of the sand‐prone moderate bioturbator. Grazing and species competition favor species distributions of dominant bioturbators. Consequently, eco‐engineering affects habitat conditions while species interactions determine species dominance.Our results show that eco‐engineering species determine the mud content of the estuary, which suggests large effects on the morphology of estuaries with aggravating habitat degradation.

Nearshore morphological modelling is challenging due to complex feedback between hydrodynamics, sediment transport and morphology bridging scales from seconds to years. Such modelling is, however, needed to assess long-term effects of changing climates on coastal environments, for example. Due to computational efficiency, the sediment transport driven by currents and waves often requires a parameterization of wave orbital velocities. A frequently used parameterization of skewness-only was found to overfeed the coast unrealistically on a timescale of years—decades. To improve this, we implemented a recently developed parameterization accounting for skewness and asymmetry in a morphodynamic model (Delft3D). The objective was to compare the effects of parameterizations on long-term coastal morphodynamics. We performed simulations with default and calibrated sediment transport settings, for idealized coastlines, and compared the results with measured data from analogue natural systems. The skewness-asymmetry parameterization was found to predict overall stable coastlines within the measured envelope with wave-related calibration factors within a factor of 2. In contrast, the original parameterization required stronger calibration, which further affected the alongshore transport rates, and yet predicted erosion in deeper areas and unrealistic accretion near the shoreline. The skewness-asymmetry parameterization opens up the possibility of more realistic long-term morphological modelling of complex coastal systems.