The Impact of Storm-Induced Breaches on Barrier Coast Systems Subject to Climate Change—A Stochastic Modelling Study

: Storms can have devastating impacts on barrier coasts causing coastal erosion, partial inundation, and possibly the breaching of barrier islands. The breaching of barrier islands provides a mechanism for the creation of new tidal inlets that connect the backbarrier basin (or lagoon) and the outer sea. As a new tidal inlet affects both the basin and the hydrodynamics of existing inlets, it is important to understand why an initial breach either closes or may evolve into a new tidal inlet. To this end, we performed a Monte Carlo analysis using an idealized model capable of simulating the long-term morphological evolution of multiple tidal inlets connected to a single backbarrier basin. To do so required the creation of a stochastic shell, as a new element around this existing barrier coast model. Our results demonstrate that barrier coast systems tend towards an equilibrium value for the number of inlets per kilometer of barrier coast and total inlet cross section. This even holds with the continuous stochastic forcing of storm-induced breaches. This ﬁnding implies that if a new breach opens in a coast that is already in equilibrium, existing inlets will shrink and may close if the new breach remains open. Furthermore, we ﬁnd that climate-driven changes in storm frequency will modify the timescales in which barrier coasts reach their equilibrium state. Finally, we ﬁnd that the distance between a new breach and its nearest neighbor is more important for its survival than the size of the breach or the degree of saturation of the barrier coast.


Introduction
Barrier coasts are important dynamic systems that cover around 10% of the coastlines worldwide [1,2]. Due to their ecologic, economic, and touristic worth, they are often densely populated [3,4]. Barrier coasts are highly dynamic systems whose morphology continuously changes under the influence of tides, waves, and storms [5]. The large-scale morphology of barrier coast systems is controlled by the relative importance of tides and waves [6]. This study focusses on meso-tidal barrier coasts, that is those in which tides and (fair-weather wind) waves are of equal importance.
A key aspect of the natural morphological evolution of barrier coast systems is the evolution of the tidal inlets, which form the connection that links the outer sea to the inner tidal basin. The impact of tides and waves on the evolution of a single tidal inlet system has already been studied by Escoffier [7] who showed that for an inlet to be in a morphological equilibrium, the sediment export due to the tide should be balanced by the wave-induced sediment import into the inlet. This equilibrium for multiple-inlet systems is also influenced by the interaction between various elements, such as between the inlets themselves and between the inlets and the backbarrier basin.
Evidence for this interaction comes from the observed residual circulation in backbarrier basins [8][9][10], the role of spatially varying backbarrier basin hydrodynamics in the stability of multiple tidal inlets [11,12], and the fact that the total tidal prism entering a backbarrier basin is primarily a system characteristic that is barely affected by the configuration of the tidal inlets [13].
Storms also significantly affect barrier coasts despite the fact that they are less frequent and more episodic than tides and waves. The effects of storms on barrier coasts range from coastal erosion and local inundation to the breaching of barrier islands [14,15]. This breaching of barrier islands is one of the mechanisms through which new tidal inlets may form. Regions that commonly experience storm-induced breaches include the east coast of the United States of America (USA) [16,17] and the Ria Formosa in Portugal [18,19]. These breaches can have a lasting impact on the entire barrier coast system, as a new inlet conveys a part of the total tidal prism entering the backbarrier basin, and alters tidal dynamics [20,21].
Examples of the breaching of barrier islands are plentiful. During Hurricane Sandy in 2012 the barrier chain south of Long Island (NY, USA) was breached at three locations, two of them at Fire Island and one at Westhampton [22]. However, at the beginning of 2020, only one of those breaches in Fire Island was still open, now called Old Inlet, whereas the other two breaches were closed artificially (see Figure 1). Another example is the Core Sound, the southern part of the North Carolina (USA) outer banks, where three inlets-Old Drum, New Drum, and Ophelia Inlet-have alternatingly opened and closed under the influence of storm-induced breaches (see Figure 2).
These examples illustrate that storm-induced breaches have a profound effect on nearby inlets and the tidal dynamics in the backbarrier basin. However, this effect is still poorly understood. Therefore, the goal of this study is to investigate the impact of storm-induced breaches on inlets connected to the same backbarrier basin and on the interaction among these inlets. Specifically, we aim to answer the following research questions: • How do storm-induced breaches affect existing inlets connecting to the same backbarrier basin and what are the interactions between them? • What breach characteristics determine whether a breach evolves into an open inlet or closes? • How are multiple-inlet systems affected by climate driven changes in the storm climate?
To answer our research questions, we need to overcome a methodical challenge as no standard approach exists to combine the two types of driving forces that we consider-storms, occurring on an irregular basis, and tides and waves that continuously affect barrier coast systems [4]. Our approach is to extend an existing idealized barrier coast model that simulates the morphological evolution of tidal inlets under the influence of tides and waves [12] by adding a stochastic shell around the model. This innovation allows the stochastic forcing of storm-induced breaches by determining whether or not breaches at each moment in the morphological evolution. We select an idealized barrier coast model as it allows us to study the effect of individual processes and its low computational cost allows us to systematically study the breach characteristics that determine whether a storm-induced breach remains open or not.
This work is organized as follows. First, the methods including the development of our stochastic shell are presented in Section 2, followed by the results in Section 3, the discussion in Section 4, and finally the conclusions in Section 5.

Methods
Previous studies on the long-term evolution of (multiple) inlet systems have used simulation models as long-term measurements are not available for a wide range of barrier coasts. Examples of models that have been used include semi-empirical models (e.g., ASMITA [24]), complex process-based models (e.g., Delft3D [25]), and idealized models (e.g., Reference [12]). However, none of these models have been combined with a stochastic forcing of storm-induced breaches. This can be partially attributed to the absence of a standard approach to combine stochastic forcing (storms) and deterministic forcing (tides and waves) in barrier coast modelling [4].
We build a stochastic shell around an existing idealized barrier coast model [12] to allow us to stochastically force storm-induced breaches. This model is outlined in Section 2.1, the stochastic shell is presented in Section 2.2, and the design of our model experiments in Section 2.3.

Idealized Barrier Coast Model
The idealized barrier coast model developed by Roos et al. [12] and described in Reef et al. [13] will be used in this study to simulate the morphological evolution of multiple tidal inlets as it combines the impact of tides and waves on the morphological evolution of inlets. The model domain is represented by a simplified geometry consisting of multiple tidal inlets that connect a rectangular inner tidal basin having a spatially uniform depth to the outer sea where the system is forced by a tidal wave. The j tidal inlets have a length l j , a cross-section A j , and constant shape defined by a shape factor γ 2 = h j /b j (e.g., Reference [5]). The rectangular tidal basin is characterized by a longshore length L, cross-shore width B, and depth h b . Finally, on the semi-infinite outer sea of depth h o , a tidal wave is forced with amplitude Z and tidal frequency ω.
The model consists of a hydrodynamic and a morphodynamic part. The morphodynamic part is based on the stability concept of Escoffier [7] and combines a constant sediment import M with a sediment export X j that is based on the velocity amplitude of the tide in an inlet. The balance of these sediment fluxes determines dA j /dt, that is an inlet accretes (M > X j ), erodes (M < X j ), or is in The hydrodynamic part is based on the linearized shallow water equations with linearized bottom friction in the basin and inlet channels according to Lorentz' linearization, and it simulates the hydrodynamics in the entire domain. It is forced by a tidal wave on the outer sea with amplitude Z and tidal frequency ω. The morphodynamic part of the model is solved numerically using a forward Euler discretization (with timestep ∆t), while the hydrodynamic part is analytically reduced to a linear system of equations that are solved using standard techniques.
In their simulations, Roos et al. [12] started with an oversaturated barrier coasts (i.e., with far more inlets than the equilibrium configuration) and simulate the morphological evolution of the tidal inlets over centuries towards an equilibrium (or near-equilibrium) configuration in which some inlets remained open and some have closed. They found that their model results agree with observed relationships between inlet characteristics and both tidal range and tidal basin width.

Stochastic Shell: Forcing Storm-Induced Breaches
To stochastically force storm-induced breaches in the model of Roos et al. [12], we develop a stochastic shell around this model. At every timestep, this shell determines whether or not storm-induced breaches are created-and, if so, what their properties are-see Figure 3. This is done in three steps: Step 1 determines the number of storms n s that occur during the timestep, Step 2 determines the number of breaches n b that occur (if n s > 0), Step 3 determines the initial inlet cross-section A init and inlet location y j of each newly created breach (if n b > 0).
In these three steps, we use three different Probability Density Functions (PDFs). In step one, to determine the number of storms n s that occur during a timestep we use a Poisson distribution as it shows the best fit with the Annual Hurricane Occurrence dataset [26], that is, Here we use a mean value λ s that is the inverse of the location specific recurrence times for hurricanes as determined by Keim et al. [27].  Because only little data is available on storm-induced breaches, and particularly about their occurrence, we assume that the occurrence and initial size can be modelled using the same PDFs as the occurrence and intensity of storms. Therefore, we also use a Poisson distribution in step two to determine the number of breaches n b per storm, that is, Here we assume a mean value λ b = 1 per storm due to a lack of data on the number of breaches per storm.
In step three, we use a Generalized Pareto distribution to determine the initial inlet cross-section A init as it was found to be the best for the intensity of storms in the Annual Hurricane Occurrence dataset [26] p(A init ) = 1 Here the scale σ = 243 m 2 , shape k = 7.84 × 10 −9 , and location θ = 38 m 2 were fitted using data from 19 historic breaches [22,23,[28][29][30][31][32][33][34][35][36]. Finally, also in the third step, we use a uniform distribution to determine the inlet location y j because we assume that the likelihood of a breach is considered the same for every part of the barrier coast (except for existing inlets). This implies p(y j ) = 0 for y j at an inlet, 1 L res for y j not at an inlet.
Here L res = L − ∑ J j=1 b j is the length of the barrier coast minus the cumulative width of all open inlets.

Outline of a Single Simulation
We start each simulation with a barrier coast without any inlets at t = 0 and subsequently force a storm event that induces at least one storm-induced breach. To do so, we start at step two in our stochastic shell and keep repeating this procedure until at least one (possibly more) storm-induced breach is created. This starting procedure is different from Roos et al. [12] where each simulation was initialized with an oversaturated barrier coast (i.e., with more inlets than in equilibrium). Next, the model simulates the evolution of the storm-induced breaches toward a final inlet configuration over a period of 1000 years.
At every timestep, two morphological changes are considered. First, the stochastic shell determines whether or not storm-induced breaches are generated (see Section 2.2). Second, similar to References [12,13] we simulate the morphological evolution of the tidal inlets based on the stability concept of Escoffier [7] (see Section 2.1).
The resulting final inlet configurations are analyzed using three metrics: the ratio of equilibrium inlet cross-section and the mean initial size A j,all /A µ (-), the number of open inlets per km barrier coast J total /L (km -1 ), and the relative total tidal prism P total /P ref,total (-) with Here, the tidal prism for a single inlet (i.e., 2|û j | ω A j ) is obtained by integrating the sinusoidal tidal velocity signal |û j | in an inlet over half a tidal cycle and multiplying the result by the inlet cross-section A j .
Two sets of parameter values are used in this work. The first represents the Great South Bay system (NY, USA) and the second the Core Sound system (NC, USA), both can be found in Table 1.

Monte Carlo Simulation
To systematically analyze the effect of storm-induced breaches we perform two Monte Carlo ensembles simulation consisting of n = 500 individual simulations, for both the Great South Bay and Core Sound sets of parameter values (see Table 1). In every simulation different storm-induced breaches are forced, therefore we expect different final configurations for every simulation as well.
To aggregate these results, we determine the median and 50% envelope over all simulations in the ensemble using the same metrics as above.
We also use these results to investigate which breach characteristics determine whether or not a storm-induced breach remains open. This is done by analyzing the effect of three breach characteristics: • the relative size of the breach A init /A µ , • the barrier coast saturation J/J µ , • the distance of a breach to the nearest neighboring inlet d nearest .
To do so, we compute 2D histograms showing the breach survival rate of all three combinations of taking two of the three breach characteristics. The chance of breach survival over the entire simulation is defined as where N open is the total number of breaches that remained open and N closed is the total number of breaches that closed, both determined per histogram bin.
To study the effect of changes in storm climate we vary the storm frequency λ s in the stochastic shell from −50% to +50% in 11 Monte Carlo ensemble simulations, consisting of 500 individual simulations each. The range of change in storm frequency from −50% to +50% is the span of the 67% confidence interval for tropical storm frequency (with a median just under 0), by the IPCC [37]. These results are analyzed by computing the same three metrics as before, as well as the timescale at which the median is at 90% of its final value (the value at t = 1000 yr).

Example Model Run
To illustrate individual simulations, Figure 4 shows two example runs, for sets of parameter values Great South Bay and Core Sound (found in Table 1) Furthermore, the example run for the Great South Bay (Figure 4 top row) shows an example of jump migration [38] where an inlet (the bottommost inlet) closes and a nearby breach remains open (effectively moving the inlet).

Monte Carlo Ensemble
To further analyze the impact of storm-induced breaches on barrier coast systems, we examine the results of two Monte Carlo ensembles for both the Great South Bay and Core Sound sets of parameter values (see Table 1).
Each Monte Carlo ensemble consists of 500 individual model runs, this number of individual model runs ensures that our results have converged, such that more runs would only change the result marginally. We find that even with the random forcing of breaches all three metrics (see Section 2.3.1) tend towards an equilibrium (see Figure 5). This means that the system reaches a dynamic equilibrium state, even with a continuous stochastic forcing of storm-induced breaches.
Comparison of the results for both sets of parameter values reveals that the equilibrium values and timescales are dependent on the set of parameter values. For all three metrics, the Core Sound system reaches its equilibrium values faster than the Great South Bay system. The former system approaches its equilibrium values very fast (the timescale in panel d & f coincides with the y-axis). The Core Sound system has less open inlets than the Great South Bay system (i.e., J total /L is lower), while the inlets are larger (i.e., A j,all /A µ is larger).
The results for the Great South Bay also show that the dimensionless total tidal prism P total /P ref,total approaches its dynamic equilibrium value first, followed by the ratio A j,all /A µ , and finally the number of inlets J total /L. Therefore, the configuration of the tidal inlets (i.e., size and spacing) takes longer to reach a dynamic equilibrium than the overall tidal prism of all inlets combined. This also explains why the ratio A j,all /A µ is initially relatively high; the dimensionless total tidal prism P total /P ref,total increases faster than the number of open inlets per km barrier coast J total /L, so each inlet conveys a larger part of the total tidal prism P total /P ref,total than in the final state.

Breach Survival Chance
Next, we investigate how the three breach characteristics already introduced in Section 2.3.2, control the chance of a breach to survive and evolve into a stable tidal inlet. This is done by analyzing all storm-induced breaches that were forced in the Monte Carlo ensemble presented in the previous section. Figure 6 shows three histograms for all possible combinations of two breach characteristics.
These histograms show that the distance to the nearest neighbor d nearest is the breach characteristic that best predicts whether a breach will remain open or not. If a barrier coast is undersaturated (i.e., J < 0.5J µ ), breaches closer to a nearest neighbor have an increased likelihood to remain open as well. Finally, the effect of the initial breach size A init is that larger breaches slightly increase the chance of survival. This pattern was observed for both the Greath South Bay (left part Figure 6) and the Core Sound (right part Figure 6) sets of parameter values.

The Effects of Climate Change
The effect of changes in storm frequency due to climate change [37] is that the timescales for the system to reach its dynamic equilibrium change, but the resulting values in dynamic equilibrium do not change. This change in timescale (here defined as time required to reach ±10% of the final value at t = 1000 yr) is shown in Figure 7. The results show that an increase (decrease) in storm frequency clearly leads to a shorter (longer) timescale for all three metrics for the Great South Bay system. The timescale for the inlet cross-section A j,all /A µ and the number of open inlets per km barrier coast J total /L to reach their equilibrium value decreased from around 700 yr for λ s = 0.5 to around 400 yr for λ s = 1.5. For the total tidal prism P total /P ref,total the timescale to reach its equilibrium decreased from over 300 yr for λ s = 0.5 to just over 100 yr for λ s = 1.5. For the Core Sound system, the timescale are very short already and barely changes due to variations in the storm frequency.  Table 1), in three histograms for varying values of relative initial breach size A init /A µ , barrier coast saturation J/J µ , and distance to nearest neighboring inlet d nearest .

Gr e a t S o u t h B a y
Co r e S o u n d  Table 1). The timescales required to reach 90% of the final value of the same three metrics as in Figure 5 are plotted while storm frequency λ s is changed between −50% and +50%. Shown are the timescales for: the ratio A j,all /A µ (a,d); the number of open inlets per km barrier coast J total /L (b,e); dimensionless total tidal prism P total /P ref,total (c,f). N.B. the timescales for λ s /λ s,0 = 1.0 corresponds to the dashed lines in Figure 5.

System-Wide Equilibrium
Our results show how a barrier coast system will reach a dynamic equilibrium in terms of inlet cross-section A j,all /A µ , the number of open inlets per km barrier coast J total /L, and total tidal prism P total /P ref,total while being stochastically forced by storm-induced breaches. Moreover, the systemwide equilibrium in total tidal prism P total /P ref,total implies that if a breach remains open once this equilibrium has been reached, it will take over part of the tidal prism from the existing inlets, potentially leading to inlets closing. The above implies that a storm-induced breach that remains open affects barrier coast systems in two ways. First, a new breach directly provides a new connection between the outer sea and basin resulting in changed hydrodynamics in both. Second, that as a result a new breach indirectly affects the existing tidal inlets and the impact that they have on the hydrodynamics in both the outer sea and basin. Importantly, these two aspects are nonlinearly coupled.
Furthermore, we found that variations in storm frequency do not lead to a different total tidal prism, only the equilibrium timescale might change. This result implies that if the storm frequency increases, barrier coast systems will reach their dynamic equilibrium faster. So, if a barrier coast is moved away from its equilibrium state (e.g., due to human interventions such as closing newly opened breaches) more effort will be required to maintain the non-equilibrium state if the storm frequency increases. These findings agree with previous findings indicating that the total tidal prism for a barrier coast system is a system characteristic [13] and thus independent of storm frequency.

Survival of the Farthest
Our analysis on which breach characteristics are adequate predictors for the survival of a storminduced breach showed that the distance to a neighboring inlet is the best predictor, followed by barrier coast saturation as second best and initial breach size as third.
The importance of a minimum distance between a new breach and a neighboring inlet indicates the need for a 'barrier' between different inlets for them both to remain stable. This finding is in agreement with previous findings of a minimum inlet spacing between multiple tidal inlets [12] and separations such as a topographic high (i.e., a tidal divide) leading to stable configurations for double inlet systems [39,40].
If a new breach close to an existing inlet drains the backbarrier basin more efficiently than the existing inlets, this could lead to the closure of the existing inlet. In this case, the new inlet succeeds in capturing part of the total tidal prism entering a backbarrier basin, which has been found to be a system characteristic [13]. Factors that were found to impact the minimum spacing of tidal inlets are the tidal range and backbarrier basin width [12,13]. Other factors that are likely to affect this minimum spacing are the depth of the basin h b , length of an inlet l j , and bottom friction in both the basin r b and inlets r j .

Model Validity & Limitations
In this work, we studied the impact of storm-induced breaches on barrier coast systems (i.e., not the actual breaching of barrier islands) and imitated the stochastic nature of storms in the idealized barrier coast model by Roos et al. [12]. Their earlier results showed a good qualitative agreement between the model results and observed relationships between inlet size and spacing and both tidal range and tidal basin width. Furthermore, Reef et al. [13] demonstrated that the model was capable of qualitatively reproducing the observed phenomena.
To assess the validity of how we stochastically simulate the storm capacity to induce breaches, we assess the results of the Monte Carlo ensemble. These results show that an equilibrium number of inlets is reached by the system, implying that most new breaches will close or cause a nearby inlet close. The example run (see Figure 4) also shows that most newly created breaches close as soon as the number of inlets has reached its equilibrium value. These results agree with observations from natural systems in which numerous tidal inlets, that started as storm-induced breaches, have closed [17]. Examples of these natural systems are the barrier coasts systems considered in this study: the Great South Bay and the Core Sound system.
As a result of incorporating the most important processes that affect the long-term morphological evolution of barrier coast systems, our model is capable of capturing the most important qualitative behavior of the natural system. To attain a better quantitative description of the natural system, future work could include more morphological features that affect inlet locations such as tidal divides [39,40], washover deposits [15,16], relict inlet features [16], channel networks [41,42], and paleographic river valleys [43] that control inlet locations. Furthermore, additional processes can be considered such as basin and outer sea morphodynamics [41,44], changing boundary conditions due to climate change [1,2], and nonlinear hydrodynamics [8].

Conclusions
We studied the impact of storm-induced breaches on the long-term evolution of multiple-inlet systems to an equilibrium configuration. To do so, we performed a Monte Carlo simulation in which storm-induced breaches were stochastically forced using a newly created stochastic shell around an idealized barrier coast model. Using parameters representative for the Great South Bay and Core Sound, we determined that most breaches close while some remain open and develop towards a new stable inlet. The formation of a new inlet has a profound impact on the barrier coast system dynamics as a whole, as a new inlet (nonlinearly) affects not only the basin hydrodynamics, but also the morphodynamics of the neighboring inlets. Despite this, the system tends towards an equilibrium in terms of number of inlets and total tidal prism, while stochastically forced by storm-induced breaches.
By analyzing all storm-induced breaches in our Monte Carlo ensemble, we identified the distance to a neighboring inlet as the breach characteristic that is the best predictor for a breach to remain open or not. Finally, we also determined that an increase in storm frequency due to climate change can lead to a decrease of the timescale over which these equilibria are reached.