The morphodynamic stability of double-barred beaches is explored using a numerical model based on linear stability analysis. Surfzone hydrodynamics is described by coupling depth and wave averaged conservation of mass and momentum with the wave-energy and phase equations, including roller dynamics. A simplified sediment transport formulation relates flow to seabed changes. Under normal wave incidence an alongshore uniform coast with a cross-shore profile characterized by the presence of two sandbars, can be unstable, thereby resulting in the development of crescentic/rip channel patterns. Our 5 study demonstrates that sandbar coupling can be either in-phase (highs and lows of both sandbars are at the same alongshore position) or out-of-phase (highs and lows of one sandbar correspond to lows and highs of the other sandbar). In line with observations, results of numerical simulations show a large variability in the possible emerging bottom patterns. Our analysis indicates that the inner bar-modes are dominant for large height/depth differences between the two sandbars crests and small offshore wave heights, while patterns related to the outer sandbar dominate for small values of the difference in sandbar depth. 10 For intermediate differences between the two sandbars depths, patterns on both longshore bars appear to be fully coupled. For relatively larger waves and large depth over the outer sandbar, patterns develop close to the shoreline/inner surfzone. 1 Introduction Multiple sandbar systems have been observed in a variety of settings worldwide. We specifically focus on the dynamics of double sandbar systems in the surfzone where the sandbars almost constantly affect (and are affected by) wave transformation and 15 onshore/offshore exchanges of sediment. Alongshore changes in double sandbar configurations sometimes result in rhythmic patterns, usually called crescentic bars or rip channels (Figure 1).The development of alongshore patterns in multiple sandbar settings has been studied through both observations (e.g. Castelle et al., 2007Castelle et al., , 2015 and numerical studies (e.g. Klein and Schuttelaars, 2006;, and has also been considered in the wider framework of a conceptual model of sequential beach changes by Short and Aagaard (1993).
20This conceptual model as well as field data (e.g. Castelle et al., 2007) indicate that the inner and outer sandbars are likely to be characterized by different spatial and temporal scales. In a double-barred system like the one considered in this study, the inner sandbar usually displays crescentic features with an alongshore spacing (distance between consecutive sandbar horns) smaller than the one characterizing the outer sandbar. Moreover, in the case of accretionary conditions, Short and Aagaard (1993) assume that the inner sandbar responds faster than the outer one. The Short and Aagaard (1993) model also indicates that 25 1 https://doi.