Reinterpreting the Bruun Rule in the Context of Equilibrium Shoreline Models

: Long-term (>decades) coastal recession due to sea-level rise (SLR) has been estimated using the Bruun Rule for nearly six decades. Equilibrium-based shoreline models have been shown to skillfully predict short-term wave-driven shoreline change on time scales of hours to decades. Both the Bruun Rule and equilibrium shoreline models rely on the equilibrium beach theory, which states that the beach proﬁle shape equilibrates with its local wave and sea-level conditions. Integrating these two models into a uniﬁed framework can improve our understanding and predictive skill of future shoreline behavior. However, given that both models account for wave action, but over different time scales, a critical re-examination of the SLR-driven recession process is needed. We present a novel physical interpretation of the beach response to sea-level rise, identifying two main contributing processes: passive ﬂooding and increased wave-driven erosion efﬁciency. Using this new concept, we analyze the integration of SLR-driven recession into equilibrium shoreline models and, with an idealized test case, show that the physical mechanisms underpinning the Bruun Rule are explicitly described within our integrated model. Finally, we discuss the possible advantages of integrating SLR-driven recession models within equilibrium-based models with dynamic feedbacks and the broader implications for coupling with hybrid shoreline models. process as a key component of this mechanism. The derivation of the


Introduction
Sandy beaches are dynamic environments, responding to a variety of complex processes interacting on different temporal and spatial scales [1]. As sea-level rise (SLR) is accelerating due to climate change [2], reliable projections of shoreline change on long (decadal to centennial) time scales are critical for coastal managers and decision makers [3,4]. SLR-driven shoreline recession occurs on time scales from decades to centuries as a result of the interaction between short-and long-term processes (e.g., wave action and SLR). Accurately integrating the SLR-driven erosion into more comprehensive shoreline change models is therefore necessary to improve our understanding and the predictability of long-term shoreline evolution in the context of climate change.
For the past 60 years, the Bruun [5] model (known as "Bruun Rule") has been the most widely used method to estimate long-term beach recession due to SLR, and its use in contemporary applications keeps growing ( Figure 1). The core assumption of the Bruun Rule is that a beach maintains a constant elevation profile and migrates upwards and landwards under the influence of SLR, based on the concept of equilibrium beach profile adjustment, under a number of assumptions (outlined in Section 2). The validity of the Bruun Rule has been recently demonstrated in laboratory settings under certain conditions [6,7].  [5] from 1962 to 2020, divided by: papers in which the Bruun Rule is employed (bright orange); non-accessible papers where the use of the Bruun Rule is uncertain (dark shaded orange); and papers citing Bruun [5] but not employing the Bruun Rule (black). Source: Google Scholar ® search engine. Bruun's model calculates shoreline recession explicitly as a function of SLR, and implicitly from the long-term effect of combined hydrodynamic processes (e.g., waves, tides, currents, etc.). To the authors' knowledge, the physical mechanisms of shoreline recession driven by the combination of waves and SLR has never been explicitly described.
In recent years, the development of equilibrium-based shoreline models (ESMs) enabled a skillful simulation of wave-driven shoreline change on cross-shore-transportdominated coasts [8][9][10][11][12]. ESMs parameterize complex sediment transport processes and rely on shoreline data to calibrate free model parameters. Such parametrization allows for the modeling of shoreline behavior on time scales ranging from hours to decades at a low computational cost. These models are conceptually similar to the Bruun Rule in that the equilibrium beach profile adjusts to the forcing conditions. ESMs model the short-term dynamic beach response to incident wave conditions rather than adjustments due to sea-level change. During storms, the upper portions of the beach profile erode, and the sediment is deposited offshore in response to increased incident wave energy (see, e.g., Yates et al. [9]). When wave height declines, the sediment is moved back to the shore, driving the profile to its original shape. ESMs show skill in predicting erosion/accretion cycles on hourly (single storm) to yearly/decadal time scales (e.g., Splinter et al. [11], Yates et al. [9], Castelle et al. [13], Lemos et al. [14]).
Applying ESMs over long time scales, where sea-level rise becomes important, is a critical step toward more comprehensive long-term shoreline projections. The most widely used solution is to couple ESMs with the Bruun model [15][16][17][18][19]. However, so far, the effects of hydrodynamic processes implicitly modeled by the Bruun Rule have never been explored in detail, and applications coupling these processes disregarded the conceptual relationship between the two models. Therefore, coupling ESMs with SLRdriven recession models requires a re-examination of the basic hypotheses and underlying physical processes of these models to verify their compatibility, modify them if necessary and lay the foundations for further developments. Such a re-examination is particularly warranted given the accelerating SLR trajectories and the increasing reliance on the Bruun Rule (see Figure 1).
In this paper, we present a novel interpretation of the physical mechanism of SLRdriven erosion in the context of equilibrium beach theory, focusing on the interaction between short-and long-term drivers (i.e., wave action and SLR). We use the proposed concept (i) to provide a new perspective on the Bruun Rule's underlying physics and (ii) to analyze the integration of SLR-driven recession into ESMs. Finally, we identify and discuss potential future research avenues to extend the validity of integrated ESMs beyond the limited set of conditions where the stand-alone Bruun Rule is applicable. The remainder of the paper includes a brief review of the Bruun model and the proposed interpretation of the generalized physical mechanism of SLR-driven shoreline recession, with an overview on ESMs and their potential role in quantifying SLR-driven shoreline change (Section 2); an analysis of the integrated SLR-driven recession and ESMs in a generalized case and under Bruun's assumptions (Section 3); and a model application to an idealized test case at the cross-shore-transport-dominated Truc Vert beach (Section 4). The results and implications for future improvements in shoreline modeling are discussed in Section 5. Conclusions are drawn in Section 6.
A list of symbols and notations (Table S1) and a list of acronyms (Table S2) used in this paper are included in the Supplementary Material. Appendix A provides mathematical insights on how the Bruun Rule can be expressed using the proposed conceptualization. Appendix B provides details on integrating ESMs and SLR-driven recession models with the dynamic feedbacks described in Section 3.2.

The Bruun Model
The Bruun Rule [5], which estimates a chronic beach recession due to SLR, is a onedimensional representation of two-dimensional processes on a cross-shore beach profile. The model assumes that: (i) a theoretical equilibrium beach profile exists and is maintained over time; (ii) sediment volume is conserved cross-shore between the limits of the active beach profile (down to a so-called "depth of closure"); (iii) SLR always induces long-term shoreline retreat (as the beach profile migrates upward and landward with SLR in order to conserve sediment volume); and (iv) there is unlimited available sediment from the subaerial beach [20,21].
The theoretical equilibrium profile shape is determined by the local sediment and wave climate characteristics [22]. Therefore, the preservation of the equilibrium profile shape implies that the local wave climate characteristics remain roughly constant.
According to the Bruun model, under the assumptions listed above, the long-term beach profile responds to SLR with a landward and upward translation (Figure 2), producing a shoreline retreat (S Bruun ) given by: where tan(α) is the mean slope of the active beach profile between the offshore limit of sediment exchange (i.e., the depth of closure) and the stable emerged beach location (e.g., emerged berm, dune toe or dune crest). The depth of closure is determined by the local sediment size and the frequency of extreme waves over a given period of time, and thus depends on the local wave climate [23]. This last point highlights another strong link between the Bruun model and the wave climate. As the depth of closure depends on the time window observed and the respective distribution of extreme wave events [24], the assumption of a constant depth of closure implies a long-term stationarity of the wave climate. The derivation of the original Bruun model was essentially based on an intuitive unification of physical beach characteristics and nearshore processes [21]. An implicit assumption of the Bruun model is that, in response to SLR, the sediment is redistributed over the profile by the long-term integrated effects of wave action [21]. This mechanism can be interpreted as the sea level setting a new non-equilibrium beach profile state, and the incident wave action mobilizing or relaxing the beach back toward equilibrium. Intuitively, if SLR occurred in the absence of sediment suspension and redistribution processes (e.g., waves), the shoreline regression would only result from the inundation of the landward profile (i.e., passive flooding [25]). This implies that the local wave climate must provide sufficient energy to mobilize sediment and reshape the profile quickly enough for the profile relaxation to keep pace with SLR [26]. Therefore, the shoreline recession modeled by the Bruun Rule (Equation (1)) could be decomposed in two contributions: an instantaneous geometric shoreline retreat due to the flooding of the beach, and a long-term profile relaxation resulting from the cumulated action of nearshore hydrodynamic processes such as wave breaking. While the first is a well understood process (known as passive flooding), there is no explicit description of the physical mechanism driving the second contribution of the Bruunian response to SLR. Figure 2. Schematic of the original Bruun model, where: SLR is the sea-level rise and S Bruun the respective shoreline recession over the considered time lap; h and L are, respectively, the vertical and cross-shore extent of the active beach, delimited seaward by the depth of closure (d c ) and landward by the berm, dune toe or dune crest; and tan(α) is the average slope of the active beach profile.
Over the past two decades, several alternative interpretations of the model's underlying assumptions were provided to critically evaluate its validity as an independent predictive tool [27][28][29][30]. However, they did not investigate this second contribution. Instead, they highlighted processes omitted by the Bruun Rule, such as the influence of local geology and sediment supply [27], landward sediment transport, dune erosion and aeolian sediment transport, and proposed modified versions of the model to include these processes with additional terms [16,26,28]. The Bruun model was developed to account only for the contribution of SLR to long-term coastal recession and does not consider other possible contributions to erosion, such as the long-term variability of the wave climate [21]. Therefore, the applicability of the Bruun model alone is restricted to a limited range of coasts [21].

A physical Interpretation of the Beach Response to Sea-Level Rise
Consider an idealized equilibrium beach profile, where the shoreline (S) is defined as the instantaneous intersection of the beach profile and mean sea level ( Figure 3a). Here, the foreshore is defined as the intertidal portion of beach between the mean high water (MHW) level and the mean low water (MLW) level ( Figure 3). With a stable mean sea level, the local wave climate produces seasonal fluctuations in the cross-shore beach profile around its theoretical equilibrium position. A rise in sea level floods the foreshore (Figure 3b), leading to a shoreline recession (Passive Flooding, S PF ) proportional to the average foreshore slope, as follows: where tan(β) is the mean foreshore slope. The shoreline recession due to passive flooding is purely geometric: it is the simple result of the change of the reference shoreline position (e.g., MSL), and occurs in the absence of sediment transport [25]. However, with respect to the new MSL (MSL' in Figure 3b), the beach slope across the profile is generally higher when the profile is concave (Figure 3b), as is the case for the majority of real beaches. A generally steeper (more accreted) profile is, on average, more balanced with low energy waves and less balanced with high energy waves [31]. Hence, the elevated mean sea-level position (MSL') brings the "beach-wave climate" system into a disequilibrium state where the current wave conditions are larger than in the equilibrium state. Consequently, continuous wave action tends to reshape the beach profile toward its equilibrium shape (Figure 3c), driving a further landward displacement of the shoreline position (Wave Reshaping, S WR ).
Recall that the beach profile in equilibrium with MSL' (Figure 3c) represents the mean position of the seasonal profile fluctuations. Therefore, depending on the energy provided by the incident wave conditions after a rise in sea level, the resulting, time-dependent effect of wave reshaping (here called d S WR ) may not immediately reach S WR (Figure 3c), but would be gradually achieved with the cumulative wave action. In this mechanism, S PF occurs simultaneously with SLR, while the associated S WR develops on the time scale of several wave events. When the Bruun Rule assumptions (Section 2.1) are satisfied, namely the assumption of stationary long-term wave climate and constant SLR rate, over periods associating no trends in wave climate (e.g., dT ≥ O(year)), the mean rates of passive flooding and wave reshaping effects are constant (see Appendix A for more details). Therefore, under these specific conditions, on the time scale of dT, the Bruun Rule can be expressed in terms of cumulated passive flooding and wave reshaping effects, as follows: Differently from the passive-flooding component (S PF ), the wave reshaping term (S WR ) results from sediment transport processes. In fact, S WR is the result of the change in wave-energy efficiency on sediment transport over a fixed amount of time, dt, due to the disequilibrium introduced by the sea-level change. Thus, quantifying S WR requires the evaluation of the disequilibrium increase induced by SLR, and the functional correspondence between shoreline response and disequilibrium state. This is discussed in the next two subsections.

Disequilibrium and Beach Profile Relaxation in Equilibrium Shoreline Models
Equilibrium shoreline models (ESMs) assume the existence of an equilibrium beach profile and rely on the concept of disequilibrium between the current beach state and the current wave conditions [31]. A general ESM schematization is given by: where F is the incident wave thrust function, ∆D the disequilibrium state, and k +/− a model free response rate parameter that has different values for accretion (∆D > 0) and erosion (∆D < 0) events, and can assume different physical meanings based on the formulation of F and ∆D (e.g., equilibrium time scale, or wave power efficiency). The time-varying disequilibrium condition ∆D can be defined in reference to the current beach profile position [8,9,12,15,17,32,33] or to the wave history [10,11,34]. Typical expressions of disequilibrium conditions based on shoreline position (∆D 1 ) and wave history (∆D 2 ) are shown in Equations (5) and (7), respectively: where E is the incident wave energy, S is the current shoreline position and E eq is the wave energy in equilibrium with S and the beach profile, i.e., the wave energy that would cause no change to the current shoreline position S, defined, for instance, by a linear relationship: with a and b being two empirical parameters, as in Yates et al. [9]; where H s is the incident significant wave height, T p is the peak wave period and (H s / T p ) eq is an equilibrium wave condition defined as function of the past wave conditions (e.g., a weighted average of the past H s / T p over a period Φ, calibrated for the site of application; see, e.g., Splinter et al. [11]). Depending on the approach, ∆D can thus be quantified in different ways, using either equilibrium wave energy or wave steepness, as in Equations (5) and (6), respectively. Each approach to evaluate ∆D (and F) is associated to a different physical interpretation to the k parameter [35]. However, regardless of the way it is defined, the disequilibrium state ∆D at a given time expresses the direction of shoreline response as well as the efficiency of the incoming wave in shaping the profile toward its theoretical equilibrium state.

Sea-Level-Rise-Driven Disequilibrium within ESMs
In the context of SLR, the explicit relationship between the disequilibrium state of the beach (∆D) and the wave-driven shoreline change provided by ESMs can be used to quantify the wave reshaping effect (dS WR ) responding to the sea-level-induced disequilibrium.
For instance, consider an ESM based on the current shoreline position (Equation (5)), where the equilibrium wave energy (E eq ) is expressed by Equation (6). An increase in MSL (dMSL) induces passive flooding (dS PF ), corresponding to a landward shift of the reference shoreline (Figure 3b), i.e., a smaller S value. Recall that this shift represents only a translation of the reference system, with no morphologic change to the beach profile. The second consequence of SLR, observed in Figure 3b, is a general increase in the bottom slope relative to the local wave climate propagating over the new MSL. In morphologic terms, this steeper beach profile is no longer in equilibrium for the mean wave climate at the site of interest. In particular, the steeper beach profile is closer to the equilibrium profile for low energy wave conditions, and farther to the equilibrium profile for high energy wave conditions [31]. If the mean wave climate does not change, wave events are more likely to result in an overall (erosive) disequilibrium of the "beach-wave climate" system ( Figure 3c). As the equilibrium condition is defined in terms of wave energy (E), the SLR-induced disequilibrium can be expressed as a decrease in equilibrium wave energy (dE eq,SLR < 0), and can be injected in Equation (6): This indicates that, for a given shoreline position (S), less wave energy is required to erode the beach and restore the equilibrium. Or, alternatively, for a given wave condition, the mean erosive efficiency of waves is increased, which likely results from an increased bottom slope across the beach profile due to SLR and wave breaking tending to occur closer to the shore.
Further insights on dE eq,SLR are provided in Section 3.2.

Integrating Sea-Level Rise in Equilibrium Shoreline Models
SLR-driven erosion can be integrated with ESMs using different methods. In turn, the ESMs' response to SLR can be different depending on how sea-level change is implemented in the model and whether the applied disequilibrium condition is influenced by the sea level. Below, we analyze approaches where wave-driven shoreline change and the effects of SLR are evaluated separately and combined linearly (Section 3.1), and approaches where the ESM's disequilibrium condition interacts dynamically with the effects of SLR (Section 3.2).

Linear Combinations of SLR Effects and ESMs
If the ESM's equilibrium formulation is based on wave history (Equation (7)), the impact of SLR can be integrated by linearly combining the ESM with passive flooding and wave reshaping models at each time step, and the differential equation for shoreline change reads: ESMs with disequilibrium conditions that are a function of shoreline position (e.g., based on Yates et al. [9], as in Equation (5)) can be sensitive to sea-level change or to the permanent flooding dS PF that follows. For this type of ESMs, the shoreline changes caused by short-term (ST) and long-term (LT) processes can be treated as independent components to avoid spurious feedbacks between the ESM and the SLR impact model: so that the ESM and SLR effects (S PF + S WR ) are computed independently over the simulated period and then simply added. In such approaches, the equilibrium energy term in Equation (5) is a function of the sole current shoreline position driven by the wave action in the absence of SLR (E eq = E eq (S)). For applications where Bruun's assumptions are satisfied, the SLR-driven shoreline retreat in Equation (9) and Equations (10)- (12) can be given by the Bruun model (Equation (3)). Such linear combination between ESMs and the Bruun model has been done in several previous applications (see, e.g., D'Anna et al. [18,19], Robinet et al. [36] and Vitousek et al. [17]). However, these studies did not explicitly acknowledge the physical mechanisms underlying the Bruun model (i.e., S PF + S WR ).

SLR Effects and ESMs with Dynamic Interactions
Below, we illustrate an alternative to the approach described by Equations (10)-(12), in which short-and long-term shoreline change are resolved in the same equation system at each time step (dt) for ESMs with disequilibrium conditions that are a function of shoreline position (e.g., Equation (5)), thereby enabling feedbacks between the modeled processes. Here, the wave reshaping (dS WR ) is explicitly computed at each time step within the ESM as a response to passive flooding (and thus to SLR). If SLR is applied at a given time step, the consequent dS PF is added on the short-term to the ESM (Equation (4)), accounting for the shoreline retreat induced by passive flooding. If the equilibrium condition is expressed as E eq (S) = aS + b (as per Yates et al. [9]), then the first (geometric) effect of passive flooding results in a disequilibrium change dE eq,SLR , which affects the wavedriven shoreline response and results in the wave reshaping (dS WR ) effect, as explained in Section 2.4. Quantifying dE eq,SLR , i.e., the change in equilibrium wave energy due to SLR, is not straightforward. To our knowledge, the literature does not provide a physical expression of the increase in wave-driven erosion efficiency as a function of SLR. Therefore, to go one step further in the quantification of this complex interaction, we make the reasonable (though approximate) assumption that the increase in wave efficiency (dE eq,SLR ) has the same magnitude as the apparent increase in wave energy associated with the passive flooding shift (i.e., dE eq,SLR = −a dS PF , which is further explored in Appendix B). Such assumption is realistic, as it translates that the steepening of the beach profile relative to the risen MSL is equivalent to the effect of an apparent accretion event that wave reshaping tends to reverse. Hence, for a given wave energy, additional erosion must occur (resulting from an elevated SLR), and, consequently, the magnitude of this wave-efficiency-driven erosion (dS WR ) is the same as dS PF . We further discuss the development of this assumption and its consistency with the Bruun Rule in Appendix B. The assumption above is based on intuitive physical considerations and is used in the present paper for demonstrational purposes in the following sections. We discuss alternative methods to estimate dE eq,SLR that do not rely on this assumption in Section 5. Here, where the disequilibrium condition is based on the shoreline position, the disequilibrium change caused by SLR (dE eq,SLR ), which produces dS WR , can be expressed adopting the assumption above and injected into Equation (5), obtaining: where the (*) symbol indicates that the disequilibrium condition responds to short-term changes in the shoreline position accounting for the geometric changes in the reference shoreline position due to passive flooding. Recall that although the additional SLR-induced disequilibrium is expressed as a function of dS PF , the resulting dS WR is of a different nature (Section 2.2). In response to the altered equilibrium owing to SLR, the model captures a timedependent, wave-driven erosion process at each modeled time step to profile reshaping (dS WR ). In this approach, the instantaneous adjustment to SLR is included by adding the sole dS PF term to the ESM and using the equilibrium condition in Equation (13): Approaches that modified the equilibrium condition of ESMs to integrate additional physical processes have been previously proposed in the literature. For example, Jaramillo et al. [37] introduced shoreline trends associated with long-term sources and sinks of sediment (e.g., gradients in longshore sediment transport or beach nourishments), using an empirical linear term. Other approaches that modified the ESM's disequilibrium condition were applied under the validity of Bruun's assumptions, although they were based on Bruun-like [15,38] or empirical terms [8,33] and did not identify separate contributions to SLR-driven erosion, such as S PF and S WR . Figure 4 illustrates the behavior of a model using the approach proposed in Equation (15) in response to a perturbation of MSL, in a simple idealized setup. Equation (15) is forced with constant wave forcing F (where F = E 0.5 and E = H s 2 /16) and a stepwise increase in sea level of 30 cm (Figure 4a). In this example, we apply realistic values of the model parameters a, b, k + and k − (see Table 1), and a mean wave height H s of 2.5 m, which are derived from the test case described in Section 4.1. Then, we repeated the application increasing and decreasing the mean wave height by 50% (i.e., 3.75 m and 1.25 m, respectively).  After a first period of stable shoreline position, in equilibrium with the applied constant wave energy, Figure 4b shows the two separate effects of passive flooding and wave reshaping, i.e., an instantaneous retreat corresponding to the immediate passive flooding induced by the stepwise SLR, followed by a gradual erosion produced by the wave action to relax the beach profile into a new equilibrium condition. Such modeled behavior shows the same characteristics of shoreline response to step-wise SLR observed in laboratory settings under stationary wave conditions [6]. In response to the change in sea level, the relaxation time of the beach profile toward the equilibrium position associated with different constant wave conditions would have been shorter (longer) for more (less) energetic waves (Figure 4b). This means that for a fast-changing sea level and persistent low-energy wave conditions, the instantaneous passive flooding and the slow development of the wave reshaping effect can be more easily isolated. Moreover, such conditions may not allow the full relaxation of the beach profile (and therefore dS WR ) until the occurrence of a high-energy wave event (e.g., storm).

Test Case
The following application aims at showing that when the assumptions underpinning the Bruun model (Section 2.1) are satisfied, the proposed mechanism explains the physics of the Bruun model shoreline retreat. We apply the variants of the combined SLR-and wave-driven shoreline change developed above, in Equations (9)-(15), to model 88 years of shoreline change at Truc Vert beach, in southwest France (Figure 5a,b), in an idealized setup that satisfies the main Bruun model's assumptions (i.e., cross-shore transport only, conservation of a long-term trend in wave energy and constant SLR rate). The following test case is intended as a proof-of-concept study rather than a site-specific hazard assessment. Hence, the results are generalizable to cross-shore-dominated beaches around the world where Bruun's conditions apply. Shoreline change at Truc Vert is essentially driven by cross-shore processes (see Castelle et al. [13], Robinet et al. [39] and D'Anna et al. [19] for more details). ESMs were found to reproduce the shoreline evolution over a decadal time scale at Truc Vert with very good results, using equilibrium based on wave history [11,13,18], but also based on the shoreline position [13,19], i.e., using Equations (6) and (5), respectively. The shoreline variability at Truc Vert beach is characterized by seasonal and interannual cycles with an amplitude of tens of meters, making this site a good benchmark for ESMs.
Historical wave conditions (significant wave height H s , peak period T p , mean direction D m ) over the period between 01 January 2012 and 31 December 2019, associated with no long-term shoreline trend on the same period, are extracted from the NORGAS-UG hindcast (Ifremer-LOPS-SHOM) at the Cap Ferret buoy located offshore, in a~50-m depth (Figure 5a). A synthetic wave time series, constructed by replicating the historical (2012-2019) time series end-to-end over 80 years, is then used to force the model from 01 January 2012 to 31 December 2099. Such a simplified approach is sufficient to compare different models under a periodic wave climate, which is required for the application of the Bruun model. Breaking wave conditions are computed using the direct formula of Larson et al. [40], which is derived from a simplified solution of the combined wave energy flux conservation equation and Snell's law in the assumption of shore-parallel bathymetry contours. SLR rates were evaluated at this site using SROCC projections [2] for the high emission RCP 8.5 scenario, considering the regional fingerprints of sea-level change contributors, over the simulation period ( Figure 6). This corresponds to a SLR rate of~1 cm/year by the mid-21st century [2,41]. The main characteristics of the active beach profile (Figure 5c) are estimated using a 2-m resolution topo-bathymetry and a high-resolution digital elevation model described in D'Anna et al. [18]. The mean profile slope (tan(α)) required by the Bruun model is calculated using the beach profile ( Figure  5c) between the dune toe and the depth of closure obtained from the Hallermeier formula (see Equation (3) of Hallermeier [23]). The latter computes the yearly depth of closure as a function of the wave height exceeded for 12 h in a year (H 0.137% ) and the corresponding wave period. Given that this calculation depends on the period of the wave time series over which extreme wave conditions are extracted, we applied the Hallermeier formula to each year of the available wave data and adopted the median of the resulting values (~14 m). The resulting mean profile slope is tan(α) = 0.014. The time average foreshore slope tan(β) is estimated from the portion of beach profile between the MHW (+1.  ESMs are run at a 1-h time step using the two different disequilibrium approaches described in Section 3. In the first approach (hereafter Y09, referring to the Yates et al. [9]-like model) the disequilibrium condition is based on a prior shoreline position (Equation (5), with E eq = aS + b, while the square root of wave energy at breaking is used as forcing (F), in line with models based on this type of disequilibrium condition [9,12,14,42,43]. The second approach (hereafter S14, referring to the Splinter et al. [11]-like model) adopts a disequilibrium condition based on the offshore wave history (Equation (7)), with (H o,s /T p ) eq defined as the weighted average of H o,s /T p over the previous Φ days, and the square root of wave power at breaking as forcing, similarly to models using this type of approach [10,11,18,36]. For each approach, the ESM k +/− and respective parameters (a and b for Y09, and Φ for S14) are calibrated using the Simulated Annealing optimization technique [44] described in D'Anna et al. [18] on Equation (4). Compared to D'Anna et al. [18], who used the Truc Vert shoreline data up to 2017, the time series of mean shoreline position estimated here with topographic surveys is extended to the end of 2019 [45]. The calibration is performed over the period 01 January 2012-31 December 2019, which shows no long-term shoreline trend (Figure 7a). For the purpose of the present idealized case application, we assumed that no SLR occurs over the calibration period, so that the models reproduce the shoreline evolution as the only response to the wave climate, with no long-term trend. The optimized parameters are reported in Table 1. The resulting Y09 and S14 models produce a Root-Mean-Square-Error of 4.68m and 6.99m, respectively, and a coefficient of determination R 2 of 0.83 and 0.62, respectively. Over the calibration period, Y09 captures the interannual shoreline variability better, producing a higher R 2 . However, both models visually demonstrate good skill and capture the overall shoreline behavior (Figure 7a). A preliminary simulation over the period 2012-2099 with no SLR is performed for each model to verify that the ESM does not produce any long-term shoreline trend in the absence of SLR ('Y09 NoSLR' and 'S14 NoSLR' in Figure 7b and Table 1).

Integration of SLR Impact Model into ESMs
Three simulations were run, coupling first S14 and Y09 with the Bruun model (B) and no feedback (Equations (9) and (12), respectively), and then the Y09 with a "passive flooding" (PF) model feeding back into the equilibrium condition (Equation (15). As the configuration of the applications satisfies the main Bruun model assumptions (no long-term variability of the wave climate, absence of processes other than SLR inducing shoreline trends, no short-term sea level fluctuations), the simulation results are expected to be consistent with the Bruun model and reproduce the cumulated effects of S PF and S WR . The combined (Y09 + PF)* models reproduce the effects of SLR on the shoreline consistently with the individual Bruun model (Equation (1)). This suggests that introducing dS PF (PF model) into the ESM, the consequent disequilibrium increase (dE eq,SLR ), which reproduces the dS WR , explains the Bruun model's results ( Figure 8b). As SLR is applied to (Y09 + PF)* at each time step, the relative passive flood dS PF and the resulting wave reshaping (dS WR ) are infinitesimal compared to the sole wave-driven shoreline change. The latter is such that the difference between the (Y09 + PF)* and (Y09 + B) results is minimal (Figure 8a,b and Table 2). However, the difference between the two model results (Figure 8c) shows seasonal fluctuations, with larger magnitudes in the summer (low wave energy) and smaller magnitudes in the winters (high energy). In fact, in the (Y09 + PF)* model the wave reshaping (dS WR ) produced in response to SLR occurs at a slower rate during low energy periods (summer) and a higher rate during high energy periods (winters) (see Section 2.2, Figure 3c) compared to the average rate produced by the Bruun Rule (Appendix A). Such seasonal behavior highlights the role of wave energy trends (in this case seasonal) on the magnitude of the modeled wave reshaping component (dS WR ). Recall that this is a theoretical scenario, in which all Bruun's assumptions are satisfied. In reality, such conditions are not granted, and other processes, such as variability in wave climate, may dominate the evolution of the beach when combined with SLR [46]. Applications of the (Y09 + PF)* to these cases are expected to capture interactions between waves and SLR and produce results that diverge from the Bruun model application.

Bruun Rule Interpretation
The idealized model applications suggest that when the Bruun Rule's assumptions are satisfied, the results obtained with (Y09 + PF)* and the Bruun Rule are consistent. This indicates that under these specific conditions the physics of the Bruun model can be explained by our proposed interpretation (i.e., dS PF + dS WR ). The latter does not aim to validate or encourage the use of the Bruun model, but highlights that, in a geomorphology context where the Bruun's assumptions are tenable, our proposed mechanism provides a physical explanation for the magnitudes and time scales of the Bruun Rule's shoreline recession. When the conditions for the application of the Bruun model are not met (i.e., short time scale, long-term variability of wave climate, presence of sources or sinks of sediment, limited accommodation space, etc.), the Bruun Rule is expected to fail [27,[47][48][49], while, as long as the dS WR component is resolved at the time scale of wave events, the physical mechanism proposed here remains valid.
For instance, current rates of SLR are in the order of mm/year and can increase to cm/year by mid-century (SROCC projections [2]). Therefore, SLR can be considered as a quasi-static process when analyzing shoreline change by tide-and wave-induced processes. However, our concept can potentially be applied to any time scale, down to a single wave event (hours), as it captures the contribution to dS WR of each wave event in response to changes in MSL. Although on time scales shorter than years the recession due to SLR (i.e., dS PF and dS WR ) is infinitesimal compared to the wave-driven shoreline change, other irregular processes acting on these time scales (e.g., wave climate variability, longshore sediment transport, etc.) and/or short-term sea-level changes (e.g., storm surges) may affect the beach equilibrium state and, in turn, the instantaneous dS PF and cumulative dS WR .
Resolving the passive flooding effects on the time scale of wave events (hours) also allows us to model the shoreline response to the lowering sea level, accounting for the varying efficiency of erosion and accretion processes. In fact, lowering sea level would reduce the flooded area (inducing a first seaward shoreline migration), while the lowered MSL associates a reduced slope across the profile (that appears more eroded), which is more balanced with high energy waves and less balanced with low energy waves [31]. Thus, the lowering sea level drives the "beach-wave climate" system into a disequilibrium state where the wave-driven-erosion efficiency is reduced. Instead, although the Bruun model can be applied in such case, the underlying physical assumptions are not designed for lowering sea levels [20,21].

Integrating ESMs and SLR-Driven Shoreline Models
Based on the proposed physical interpretation of the SLR-driven shoreline response, our analytical discussion identifies multiple generalized approaches to integrate SLR effects into ESMs where passive flooding and the wave reshaping effect are explicitly resolved. The analysis shows that two approaches can be adopted to include the shoreline relaxation that follows passive flooding (Figure 3c, Sections 2.2 and 2.4) into ESMs. If the ESM disequilibrium condition is based on wave history (e.g., Equation (7)), the integration of dS PF does not alter the disequilibrium, so that the term dS WR should be computed separately and added. Instead, ESMs relying on the shoreline position to define the disequilibrium condition may require or not the additional dS WR term, depending on whether feedbacks between model equations are enabled or not. In the first two instances (ESMs based on wave history or shoreline position with no feedback) a specific function that quantifies dS WR in response to SLR is required. Instead, when feedback is allowed within ESMs based on previous shoreline positions (Equation (15), the model produces dS WR , provided that a relationship between dS PF (or SLR) and dE eq,SLR is assigned. Here, particular attention must be paid to the potential role of wave energy trends in the modeling of dS WR .
It must be noted that using the Bruun model with ESMs [17,19,36,38] assumes a constant active profile, which implicitly assumes a stationary wave climate. Therefore, applications to regions subject to statistically significant wave climate changes must be performed with care. The available global studies (0.5 • to 1.5 • spatial resolution) indicate significant changes in wave climate offshore (~200 m depth) of~50% of the world's coastline over the next century (COWCLIP, [50]), suggesting that significant nearshore wave changes may be expected in some localized areas. In addition, when addressing future shoreline variability (rather than overall trends), changes in the seasonal and interannual variability of the wave climate as well as storminess could modify the shoreline behavior [51,52], and should be accounted. Coupling SLR and ESM models with feedback resolves all SLR-induced shoreline change components on the same short time scale, reproducing the wave-event by wave-event response of the shoreline to passive flooding and capturing the joint effects of SLR with possible changes in the wave energy distribution.

Applications of Integrated Shoreline Models and Uncertainties
While ESMs are generally applicable to cross-shore-transport-dominated open beaches, the quantity, frequency and accuracy of the available shoreline data at a given site are a potential source of model uncertainty [18,19,52,53]. In real case applications of ESMs (with or without SLR-driven erosion component), such uncertainties can be empirically quantified and introduced in the frame of probabilistic shoreline predictions [18,19,[54][55][56].
Here, however, we performed an idealized model application to compare the behavior of different model combinations, and the evaluation of model uncertainty, as well as a full validation of the wave-driven models, is outside the scope of this work.
Our idealized model application is forced using a synthetic wave time series obtained by replicating 8 years of wave hindcast data. Although this is a simplistic approach, such wave series ensures that the Bruun's assumption of steady wave climate is satisfied. In real case applications, where the Bruun's assumptions are not necessarily satisfied, the model should be forced with realistic wave time series. Given the key role of short-and long-term variability of the incident wave energy in equilibrium-based models (integrated or not), uncertainties on the inherent variability of the wave climate should be considered [19,35]. The latter can be addressed by performing ensemble simulations using a large number of different plausible wave series generated, for instance, with statistical methods [35,57,58]. However, in this proof-of-concept study we aim to analyze the behavior of integrated models in an idealized context where Bruun's assumptions are satisfied, and the construction of a more realistic synthetic wave time series is beyond our scope.
The SLR-induced shoreline change component is very sensitive to the slope used into the SLR impact model (tan(α) and tan(β)). The Bruun model slope (tan(α)) is determined by the limits of the active beach profile. Local scale bathymetric surveys are scarce on most sandy coasts and do not allow for a reliable validation of the seaward limit of the active beach profile. The available global wave and topo-bathymetric merged datasets and local scale bathymetric surveys currently allow for estimates of tan(α) [59]. However, the accuracy of such data is limited owing to the spatiotemporal variability of the available merged data, merging procedures and estimation of the depth of closure [59,60]. Therefore, the Bruun slope tan(α) is typically associated to large uncertainties [24,28,61]. The foreshore beach slope, which is adopted in the current model, does not require the estimation of active beach boundaries and can be extracted from topographic profiles (regular enough to account for potential temporal and spatial variability) that can be more easily surveyed by various means at a low cost. Further, recently developed satellite-remote sensing techniques to derive foreshore beach slopes over large scales [62] may represent a promising complement to the model developed here. Thus, the ESMs combined with SLR-driven recession (including feedbacks) may result in reduced uncertainty in long-term shoreline change and a potential of application to a much wider range of coasts compared to the models based on the linear sum of the Bruun and ESM projections.

Beyond the Bruun Rule
The conceptualization of SLR-driven erosion provided in this contribution identifies the wave-reshaping process as a key component of this mechanism. The derivation of the wave-efficiency change requires an assumption on the complex interaction between the waves and the changing sea level, which amplifies the impact of the waves on the shoreline position. In this paper, the SLR-induced disequilibrium that drives wave-reshaping erosion is quantified based on the assumption that such disequilibrium corresponds to the effect of passive flooding on the shoreline. Future research efforts may aim at providing a generalized physics-based expression of this disequilibrium, for instance by investigating wave energy dissipation change in response to the changing bottom slope. For instance, Davidson [63] recently proposed a model where the influence of the changing sea level and bottom slope on wave energy dissipation is accounted.
It is to be noted that the aim of this paper is to present a proof of concept, and the next challenges include a validation of the concept against observations. Such validation would require a long-term (e.g., decadal) time series of accurate shoreline position (for calibration purposes) with persistent low-energy incident waves (e.g., <1m), a fast-changing relative sea level (of the order of centimeters per year) and a relatively stable sediment budget to depict SLR impacts at a given field site. To the best of the authors' knowledge, such dataset does not exist. Another avenue is to address the impact of sea-level anomalies of the order of tenths of centimeters on time scales of weeks to months, as observed off the coast of Japan [64] or in El Niño Southern Oscillation [65], associated with frequent shoreline measurements. Alternatively, laboratory data produced during reduced scale studies assessing the effects of SLR [6,7] may be used for the validation of the present framework.

Conclusions
The Bruun Rule is increasingly used in current coastal impact models, whether standalone (regional to global scale) or coupled to ESMs (local scale), often without full consideration of its underlying physics. This paper introduces a new interpretation of the SLR-driven shoreline response mechanism based on equilibrium beach theory, identifying two main contributing processes: passive flooding and wave reshaping. It is shown that, when Bruun's assumptions are satisfied, the proposed conceptualization provides a physics-based background to the Bruun model. Our concept also provides support to the integration of SLR impact into ESMs for future studies, and new ground for the development of more generalized shoreline models. We also showed that ESMs can be used to explicitly compute the portion of sea-level-driven shoreline change induced by the change in wave-erosion efficiency. In particular, using a disequilibrium condition based on the current shoreline position can be efficiently coupled (with feedback) with the SLR-induced passive flooding. Such modeling does not require the knowledge of the active profile (Bruun) slope but, instead, the foreshore slope, which can be easily measured and is typically associated with less uncertainties than the Bruun slope. This can ultimately avoid bias in long-term shoreline projection considering climate change and the sea-level rise.
We addressed the coupling of SLR-driven erosion with ESMs, for cross-shore-transportdominated coasts. However, our analysis also holds for the new generation of hybrid models that couple ESM with other transport processes (e.g., longshore transport), local geology and/or the backing coastal dunes (e.g., [17,36,38,66]), and to ESMs that continuously adapt the shoreline response to the non-stationary wave climate through time-varying model parameters [52]. In such models, not only SLR but also other processes can affect shoreline position and, in turn, the equilibrium profile conditions. We therefore anticipate that our findings will help future studies addressing the long-term evolution of wave-dominated sandy coasts as sea level rises.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/jmse9090974/s1, Table S1: List of symbols and notations used in this paper, Table S2: List of acronyms used in this paper. Data Availability Statement: This study includes the beach monitoring data study of Truc Vert labeled by the Service National d'Observation (SNO) Dynalit (https://www.dynalit.fr), with additional support from Observatoire Aquitaine de l'Univers (OASU) and Observatoire de la Côte Aquitaine (OCA).

Acknowledgments:
The authors thank their colleagues, including S. Bujan, S. Ferreira and V. Marieu, involved in the topographic data, and J. Rohmer and R. Thiéblemont for earlier discussions.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Passive Flooding, Wave Reshaping and the Bruun Rule
In this Appendix we show how expressing the Bruun Rule as the long-term cumulated effects of passive flooding and wave reshaping implies constant mean rates of these two contributions over time. The physical mechanism proposed in Section 2.2 indicates that SLR induces a shoreline change driven by the cumulated effects of passive flooding (S PF ) and wave reshaping (S WR ) over a time T: When the Bruun Rule assumptions (Section 2.1) are satisfied, namely the assumption of stationary long-term wave climate and constant SLR rate, over periods associating no trends in wave climate (e.g., dT ≥ O(year)), the mean rates of passive flooding and wave reshaping effects are constant over time: Hence, over a time ∆T = n dT (with n positive integer), Equation (A1) becomes: Therefore, under these specific conditions, on the time scale of dT (years or longer), the Bruun Rule can be expressed in terms of cumulated passive flooding and wave reshaping effects, as follows:

Appendix B. Integrating SLR-Driven Disequilibrium into ESMs with Feedback
In this Appendix, we provide a more detailed description of the approach introduced in Section 3.2, where SLR-driven erosion and an ESM based on the current shoreline position are integrated, enabling feedback between SLR effects and the disequilibrium condition. Then, we further analyze the assumption used to quantify the disequilibrium induced by sea-level rise (SLR), and we analytically show that such assumption is consistent with the Bruun model structure.
Given an ESM based on the current shoreline position, if the equilibrium energy in Equation (5) (E eq = E eq (S)) is expressed by a linear relationship between incident wave energy (E) and current shoreline cross-shore distance from a given reference (S) (i.e., E eq = aS + b, after Yates et al. [9]), the instantaneous disequilibrium state can be represented on the S-E Cartesian plane ( Figure A1). In Figure A1a, a given shoreline position is associated with an equilibrium wave energy (E eq ) that would not produce a shoreline change (black solid line), while positions located above and below the equilibrium line are associated with erosion (dS/dt < 0, red points) and accretion (dS/dt > 0, blue points), respectively. An increase in MSL (dMSL) induces passive flooding (dS PF ), corresponding to a landward shift of the reference shoreline ( Figure 3b of the main text), i.e., a smaller S value. Recall that this shift represents only a translation of the reference system with no morphologic changes to the beach profile. Therefore, the effect of passive flooding is projected on the S-E plane as a landward shift of the equilibrium shoreline position (red arrows and red point in Figure A1b). The second consequence of dMSL (observed in Figure 3b of the main text) is a general increase in the bottom slope relative to the local wave climate propagating over the new MSL. In morphologic terms, in respect to a given equilibrium profile, a generally steeper beach profile is closer to equilibrium with low energy wave conditions and farther to high energy wave conditions [31], resulting in an overall (erosive) disequilibrium of the beach-wave climate system ( Figure 3c of the main text). As the equilibrium condition is defined in terms of wave energy (E), the SLR-induced disequilibrium is projected on the S-E plane as a downward shift of the equilibrium line (dE eq,SLR and dashed black line in Figure A1b), indicating that for a given shoreline position (S) less wave energy is required to erode the beach and restore the equilibrium. Or, alternatively, for a given wave climate, the efficiency of erosion is increased, which can result from an increased depth and bottom slope across the beach profile due to SLR and wave breaking tending to occur closer to the shore. The color-bar indicates the magnitude of potential erosion (red) and accretion (blue). Red arrows indicate the geometric shift due to Passive Flooding (dS PF ); vertical black arrows indicate the disequilibrium change (dE eq,SLR ) corresponding to the change in the bottom slope relative to MSL; and the horizontal black arrow indicates the potential Wave Reshaping effect (dS WR ) that would re-establish the equilibrium state. The black dash line indicates the new equilibrium beach-wave condition post-SLR. The data used for this example refer to Truc Vert beach (France).
Quantifying dE eq,SLR is not straightforward. To our knowledge, the literature does not provide a physical expression of the increase in wave-driven erosion efficiency as a function of SLR. Therefore, to go one step further in the quantification of this complex interaction, we make the reasonable (though approximate) assumption that the increase in wave efficiency (dE eq,SLR ) has the same magnitude as the apparent wave-energy increase associated with the passive flooding shift on the S-E plane (dE PF = −a dS PF in Figure A1b). Hence, for a given wave energy, additional erosion must occur and, consequently, the magnitude of this wave-efficiency-driven erosion (dS WR ) is the same as dS PF . The latter implication is consistent with the Bruun Rule's structure, as demonstrated below.
Given that E eq = aS + b, the increase in disequilibrium (dE eq,SLR ) that results in dS WR is: The Bruun Rule can be expressed in terms of passive flooding and wave reshaping, as follows: S Bruun = SLR tan(α) = S PF + S WR (A6) As described in Section 2.2, the passive flooding contribution can be computed as: Where tan (β) is the slope of the beach foreshore.
In order for the combined recession to be consistent with the Bruun Rule recession (Equation (A6)), S WR can be expressed as: and, assuming c = 1 ((tan(β))/(tan(α)) = 2), we obtain dE eq,SLR = −a dS PF = dE PF (A13) Recall that, while realistic, the assumption described herein (Equation (A13)) is only used to provide an approximate quantification of dE eq,SLR . Here, we adopt Equation (A13) in order to analyze the behavior of an ESM where feedbacks between the modeled processes are enabled, but it may not be applicable in the generality of cases in practice.
In Section 3.2 of the main text, we proposed an approach for integrating SLR-driven recession into ESMs based on the current shoreline position, where feedbacks between modeled processes are allowed. In particular, we focused on ESMs where the equilibrium condition is expressed as E eq = aS + b.
For this case, if SLR is applied at a given time step, the consequent dS PF is added on the short-term to the ESM, accounting for the shoreline reference retreat induced by passive flooding. This first (geometric) effect of passive flooding is accounted in the equilibrium condition as a landward shift of the reference shoreline position ( Figure A1). The second effect of SLR (and thus passive flooding) is a disequilibrium increase due to the apparent increase in the beach profile slope, which ultimately results in dS WR . Therefore, the new equilibrium condition (E* eq ) can be rewritten as: where the term "ΣdS P " represents the cumulative landward shift of the reference shoreline position, while the term "ΣdE eq,SLR " quantifies the cumulative disequilibrium change induced by SLR.
Under the assumption described above, the disequilibrium increase (dE eq,SLR ) can be expressed by Equation (A5) and can be injected in Equation (A14), giving: Recall that although the additional SLR-induced disequilibrium is expressed as a function of dS PF , the dS WR resulting from it is of a different nature (Section 2.2). In response to the altered disequilibrium, the model captures the contribution of waves at each modeled time step to profile reshaping (dS WR ). In this approach, the effect of SLR is included correctly by adding the sole dS PF term to the ESM and using the equilibrium condition in Equation (A13): The presence of the dS PF term in Equation (A16) is so that the shoreline adjusts instantly to the passive flooding and evolves over time due to wave reshaping, consistently with the conceptual model.
We note that the magnitude of dS WR at each time step is quantified in Equation (A16) as a function of the sea-level-driven disequilibrium change and the instantaneous wave intensity, as follows: dS WR = k +/− F dE eq,SLR dt (A17)