Intertidal mudflats are morphodynamic features present in many estuaries worldwide. Often located between vegetated shores and deep channels they comprise valuable ecosystems and serve to protect the hinterland by attenuating waves. Although mudflats are persistently present on yearly to decadal time scales, little is known on their morphodynamic adaptation to short-term variations in forcing such as storms, spring-neap tidal cycles, and sediment supply. This study aims to explore the morphodynamic resilience of mudflats to seasonal variations in forcing. First, we compare transects observed in South Bay, California, at 3- to 6-monthly intervals. Second, we present the results of a process-based, morphodynamic profile model (Mflat). Mflat is an open source, Matlab code that describes both cross-shore and alongshore tidal hydrodynamics as well as a stationary wave model. An advection-diffusion equation solves sediment transport while bed level changes occur by the divergence of the sediment transport field. Mflat reproduces the observed South San Francisco Bay profile in equilibrium with significant skill. Short-term variations in hydrodynamic forcing and sediment characteristics disturb the profile mainly at the channel-shoal edge. The modeled profile disturbance is consistent with observations. The modeled profile is remarkably resilient since it recovers to the equilibrium profile within weeks to months. The model results suggest that 3-monthly observation intervals are probably too long to discriminate processes responsible for the profile disturbance. These processes may include variations in sediment supply, mudflat erodibility, and wave action as well as the spring-neap tidal cycle.

@en

Anticipated sea level rise (SLR) threatens intertidal areas and associated ecosystems in estuaries worldwide. There is a need to develop validated modeling tools to assess the impact of SLR on estuarine morphodynamics. This study explores the morphological impact of SLR on a channel-shoal system in San Pablo Bay, a subembayment of San Francisco Bay, California, using a 3-D, process-based modeling approach (Delft3D) including density currents and wave action. The Bay underwent considerable morphologic development in response to variations in fluvial sediment load and discharge associated with a period of hydraulic mining for gold and later damming in the watershed. The availability of a unique 150-year, 30-year sequenced, bathymetric data set provided a rare opportunity for model validation. We investigate a 250-year period of morphodynamic evolution including a 150-year hindcast and a 100-year forecast with different SLR scenarios. The model shows significant skill in hindcasting volumes and patterns of bathymetric development during both net depositional (1856–1951) and erosional (1951–onward) periods. Forecasts show that SLR alters the Bay's erosional trend to a depositional trend again. Despite increased sediment trapping rates, the intertidal mudflats drown under all modeled SLR scenarios (42, 84, and 167 cm by end of the 21st century). Our work highlights the potential of using process-based models to assess the morphodynamic impact of SLR. The study also suggests that SLR can greatly increase the loss of intertidal area when landward migration is not possible. Sustainable management strategies are required to safeguard these valuable intertidal habitats.

@en

Suspended sediment concentration is an important estuarine health indicator. Estuarine ecosystems rely on the maintenance of habitat conditions, which are changing due to direct human impact and climate change. This study aims to evaluate the impact of climate change relative to engineering measures on estuarine fine sediment dynamics and sediment budgets. We use the highly engineered San Francisco Bay-Delta system as a case study. We apply a process-based modeling approach (Delft3D-FM) to assess the changes in hydrodynamics and sediment dynamics resulting from climate change and engineering scenarios. The scenarios consider a direct human impact (shift in water pumping location), climate change (sea level rise and suspended sediment concentration decrease), and abrupt disasters (island flooding, possibly as the results of an earthquake). Levee failure has the largest impact on the hydrodynamics of the system. Reduction in sediment input from the watershed has the greatest impact on turbidity levels, which are key to primary production and define habitat conditions for endemic species. Sea level rise leads to more sediment suspension and a net sediment export if little room for accommodation is left in the system due to continuous engineering works. Mitigation measures like levee reinforcement are effective for addressing direct human impacts, but less effective for a persistent, widespread, and increasing threat like sea level rise. Progressive adaptive mitigation measures to the changes in sediment and flow dynamics resulting from sea level rise may be a more effective strategy. Our approach shows that a validated process-based model is a useful tool to address long-term (decades to centuries) changes in sediment dynamics in highly engineered estuarine systems. In addition, our modeling approach provides a useful basis for long-term, process-based studies addressing ecosystem dynamics and health.

@en

Morphodynamic equilibrium is a widely adopted yet elusive concept in the field of geomorphology of coasts, rivers and estuaries. Based on the Exner equation, an expression of mass conservation of sediment, we distinguish three types of equilibrium defined as static and dynamic, of which two different types exist. Other expressions such as statistical and quasi-equilibrium which do not strictly satisfy the Exner conditions are also acknowledged for their practical use. The choice of a temporal scale is imperative to analyse the type of equilibrium. We discuss the difference between morphodynamic equilibrium in the "real world" (nature) and the "virtual world" (model). Modelling studies rely on simplifications of the real world and lead to understanding of process interactions. A variety of factors affect the use of virtual-world predictions in the real world (e.g., variability in environmental drivers and variability in the setting) so that the concept of morphodynamic equilibrium should be mathematically unequivocal in the virtual world and interpreted over the appropriate spatial and temporal scale in the real world. We draw examples from estuarine settings which are subject to various governing factors which broadly include hydrodynamics, sedimentology and landscape setting. Following the traditional "tide-wave-river" ternary diagram, we summarize studies todate that explore the "virtual world", discuss the type of equilibrium reached and how it relates to the real world.@en

Model results show different adaptation time scales of the channel mudflat system, but no equilibrium in the strict sense. Our numerical, process-based modeling approach allows for a detailed analysis of the underlying processes, which is part of ongoing efforts.@en

Peak river flows transport fine sediment, nutrients, and contaminants that may deposit in the estuary. This study explores the importance of peak river flows on sediment dynamics with special emphasis on channel network configurations. The Sacramento-San Joaquin Delta, which is connected to San Francisco Bay (California, USA), motivates this study and is used as a validation case. Besides data analysis of observations, we applied a calibrated process-based model (D-Flow FM) to explore and analyze high-resolution (∼100 m, ∼1 h) dynamics. Peak river flows supply the vast majority of sediment into the system. Data analysis of six peak flows (between 2012 and 2014) shows that on average, 40 % of the input sediment in the system is trapped and that trapping efficiency depends on timing and magnitude of river flows. The model has 90 % accuracy reproducing these trapping efficiencies. Modeled deposition patterns develop as the result of peak river flows after which, during low river flow conditions, tidal currents are not able to significantly redistribute deposited sediment. Deposition is quite local and mainly takes place at a deep junction. Tidal movement is important for sediment resuspension, but river induced, tide residual currents are responsible for redistributing the sediment towards the river banks and to the bay. We applied the same forcing for four different channel configurations ranging from a full delta network to a schematization of the main river. A higher degree of network schematization leads to higher peak-sediment export downstream to the bay. However, the area of sedimentation is similar for all the configurations because it is mostly driven by geometry and bathymetry.@en

Estuarine tidal mudflats form unique habitats and maintain valuable ecosystems. Historic measurements of a mudflat in San Fancsico Bay over the past 150 years suggest the development of a rather stable mudflat profile. This raises questions on its origin and governing processes as well as on the mudflats’ faith under scenarios of sea level rise and decreasing sediment supply. We developed a 1D morphodynamic profile model (Delft3D) that is able to reproduce the 2011 measured mudflat profile. The main, schematised, forcings of the model are a constant tidal cycle and constant wave action. The model shows that wave action suspends sediment that is transported landward during flood. A depositional front moves landward until landward bed levels are high enough to carry an equal amount of sediment back during ebb. This implies that, similar to observations, the critical shear stress for erosion is regularly exceeded during the tidal cycle and that modelled equilibrium conditions include high suspended sediment concentrations at the mudflat. Shear stresses are highest during low water, while shear stresses are lower than critical (and highest at the landward end) along the mudflat during high water. Scenarios of sea level rise and decreasing sediment supply drown the mudflat. In addition, the mudflat becomes more prone to channel incision because landward accumulation is hampered. This research suggests that sea level rise is a serious threat to the presence of many estuarine intertidal mudflats, adjacent salt marshes and their associated ecological values.

@en

In estuaries suspended sediment concentration (SSC) is one of the most important contributors to turbidity, which influences habitat conditions and ecological functions of the system. Sediment dynamics differs depending on sediment supply and hydrodynamic forcing conditions that vary over space and over time. A robust sediment transport model is a first step in developing a chain of models enabling simulations of contaminants, phytoplankton and habitat conditions. This works aims to determine turbidity levels in the complex-geometry delta of the San Francisco estuary using a process-based approach (Delft3D Flexible Mesh software). Our approach includes a detailed calibration against measured SSC levels, a sensitivity analysis on model parameters and the determination of a yearly sediment budget as well as an assessment of model results in terms of turbidity levels for a single year, water year (WY) 2011. Model results show that our process-based approach is a valuable tool in assessing sediment dynamics and their related ecological parameters over a range of spatial and temporal scales. The model may act as the base model for a chain of ecological models assessing the impact of climate change and management scenarios. Here we present a modeling approach that, with limited data, produces reliable predictions and can be useful for estuaries without a large amount of processes data.

@en