Intertidal shoals are key features of estuarine environments worldwide. Climate change poses questions regarding the sustainability of intertidal areas under sea-level rise (SLR). Our work investigates the SLR impact on the long-term morphological evolution of unvegetated intertidal sandy shoals in a constrained channel-shoal system. Utilizing a process-based model (Delft3D), we schematize a short tidal system in a rectangular (2.5 × 20 km) basin with a high-resolution grid. An initial, mildly sloping, bathymetry is subjected to constant semidiurnal tidal forcing, sediment supply, and small wind-generated waves modeled by SWAN. A positive morphodynamic feedback between hydrodynamics, sediment transport, and morphology causes the emergence of large-scale channel-shoal patterns. Over centuries, tide-residual sediment transport gradually decreases leading to a state of low morphological activity balanced by tides, waves, and sediment supply. Tidal currents are the main driver of the SLR morphodynamic adaptation. Wave action leads to wider and lower shoals but does not fundamentally change the long-term morphological evolution. SLR causes increased flood dominance which triggers sediment import into the system. Shoals accrete in response to SLR with a lag that increases as SLR accelerates, eventually causing intertidal shoals to drown. Seaward shoals near the open boundary sediment source have higher accretion rates compared to landward shoals. Similarly, on a shoal-scale, the highest accretion rates occur at the shoal edges bounding the sediment suppling channels. A larger sediment supply enhances the SLR adaptation. Waves help distribute sediment supplied from channels across shoals. Adding mud fractions leads to faster, more uniform, accretion and muddier shoals under SLR.

Intertidal shoals are pronounced morphological features found in many estuaries worldwide. Apart from maintaining an ecologically unique intertidal environment, shoals also protect adjacent dyke systems by attenuating waves. The fate of sandy shoals under anticipated sea level rise (SLR) scenarios is underexplored. The current research investigates the long-term morphodynamic evolution of estuarine sandy shoals under forcing by short fetch, locally generated wind-waves, tides, and SLR by means of a numerical, process-based model (Delft3D). The focus lies on a sheltered shoal complex in the Western Scheldt, the Netherlands. Starting from the initial, 1963 bathymetry, we model 50-year morphodynamic development with schematized wind-wave forcing. We analyze in detail the impact of locally generated wind-waves on shoal formation. Finally, we impose regional SLR of 1.10 m and 1.95 m for 100 years. Model results show that, on the spatial scale of intertidal flats, small, locally generated wind-waves lower and widen the shoals while the adjacent channels deepen. However, on the estuarine system scale, wind-waves do not lead to fundamentally different channel–shoal patterns and morphodynamic evolution trends. This suggests that channel–shoal formation is mainly due to tide residual sediment transports, with wind-waves playing a secondary role. SLR leads to a notable intertidal area loss, despite a continuous heightening of the shoals, implying that morphodynamic adaptation lags behind SLR. The inclusion of wind-waves does not fundamentally change the reaction of the estuarine shoal to SLR. Future research may focus on exploring the impact of including multiple sediment classes.

Intertidal shoals are vital components of estuaries. Tides, waves, and sediment supply shape the profile of estuarine shoals. Ensuring their sustainability requires an understanding of how such systems will react to sea level rise (SLR). In contrast to mudflats, sandy shoals have drawn limited attention in research. Inspired by a channel-shoal system in the Western Scheldt Estuary (Netherlands), this research investigates governing processes of the long-term morphodynamic evolution of intertidal estuarine sandy shoals across different timescales. We apply a high-resolution process-based numerical model (Delft3D) to generate a channel-shoal system in equilibrium and expose the equilibrium profile to variations in wave forcing and SLR. Combined tidal action and wave forcing initiate ridge formation at the seaward shoal edge, which slowly propagates landward until a linear equilibrium profile develops within 200 years. Model simulations in which forcing conditions have been varied to reproduce observations show that the bed is most dynamic near the channel-shoal interface. A decrease/increase in wave forcing causes the formation/erosion of small tidal levees at the shoal edge, which shows good resemblance to observed features. The profile recovers when regular wave forcing applies again. Sandy shoals accrete in response to SLR with a long (decades) bed-level adaptation lag eventually leading to intertidal area loss. This lag depends on the forcing conditions and is lowest near the channel and gradually increases landward. Adding mud makes the shoal more resilient to SLR. Our study suggests that processes near the channel-shoal interface are crucial to understanding the long-term morphodynamic development of sandy shoals.

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.