Uncertainties in Sandy Shorelines Evolution under the Bruun Rule Assumption

In the current practice of sandy shoreline change assessments, the local sedimentary budget is evaluated using the sediment balance equation, that is, by summing the contributions of longshore and cross-shore processes. The contribution of future sea-level rise induced by climate change is usually obtained using the Bruun rule, which assumes that the shoreline retreat is equal to the change of sea-level divided by the slope of the upper shoreface. However, it remains unsure that this approach is appropriate to account for the impacts of future sea-level rise. This is due to the lack of relevant observations to validate the Bruun rule under the expected sea-level rise rates. To address this issue, this article estimates the coastal settings and period of time under which the use of the Bruun rule could be (in)validated, in the case of wave-exposed gently-sloping sandy beaches. Using the sedimentary budgets of Stive (2004) and probabilistic sea-level rise scenarios based on IPCC, we provide shoreline change projections that account for all uncertain hydrosedimentary processes affecting idealized low- and high-energy coasts. Hence, we incorporate uncertainties regarding the impacts of longshore processes, sea-level rise, storms, aeolian, and other cross-shore processes. We evaluate the relative importance of each source of uncertainties in the sediment balance equation using a global sensitivity analysis. For scenario RCP 6.0 and 8.5 and in the absence of coastal defenses, the model predicts a perceivable shift toward generalized beach erosion by the middle of the 21st century. In contrast, the model predictions are unlikely to differ from the current situation in case of scenario RCP 2.6. Finally, the contribution of sea-level rise and climate change scenarios to sandy shoreline change projections uncertainties increases with time during the 21st century. Our results have three primary implications for coastal settings similar to those provided described in Stive (2004) : ﬁrst, the validation of the Bruun rule will not necessarily be possible under scenario RCP 2.6. Second, even if the Bruun rule is assumed valid, the uncertainties around average values are large. Finally, despite these uncertainties, the Bruun rule predicts rapid shoreline retreat of sandy coasts during the second half of the 21st century, if greenhouse gas concentration in the atmosphere are not drastically reduced (scenarios RCP 4.5, 6.0, and 8.5).


INTRODUCTION
One of the most important challenge for coastal adaptation is the rise of sea-level caused by anthropogenic climate change (Slangen et al., 2014b;Dangendorf et al., 2015). At present, there are already evidences that extreme water levels are becoming higher and more frequent (Marcos et al., 2009;Menéndez and Woodworth, 2010;Woodworth et al., 2011;Woodworth and Menéndez, 2015). At longer timescales, contributions from polar ice-sheets largely exceeding one meter are now recognized possible (Golledge et al., 2015;Hansen et al., 2015;Winkelmann et al., 2015). Meanwhile, increased rates of shoreline retreatand particularly of sandy beaches -are expected to take place (Nicholls and Cazenave, 2010).
A common approach to address sub-to multi-decadal shoreline variability on wave-exposed sandy coast with infinite sand availability is to use the sediment balance equation, which can be written as follows (e.g., Cowell et al., 2003b;Stive, 2004;Yates et al., 2011;Aagaard and Sørensen, 2013;Anderson et al., 2015): where: • S is the cross-shore shoreline displacement over a given period of time (typically a few decades), • ξ is the change of mean sea-level over the same period of time, • tan(β) is the average slope between the top of the beach and the closure depth, also refered to as active profile or upper shoreface, • f cross−shore and f longshore are the contributions of other processes causing losses or gains of sediments in the active beach profile.
Equation (1) applies to the upper shoreface of the coastal tract (Cowell et al., 2003a). It assumes that the upper shoreface keeps the same profile and translates seaward or landward depending on the sediment budget. In particular, the term ξ/tan(β) represents the impacts of sea-level change and corresponds to the Bruun rule (Bruun, 1962). While this term has been subject of debate over the last decades (Cooper and Pilkey, 2004;Ranasinghe and Stive, 2009;Woodroffe and Murray-Wallace, 2012;Passeri et al., 2015), there is no clear recommendation to leave it out. Recent results regarding the global impact of sealevel rise on shoreline change are largely based on the Bruun rule (Hinkel et al., 2013). Alternative approaches exist, but they are more complex and they require more data (Ranasinghe et al., 2012;Wainwright et al., 2015). Finally, even if the concepts of the Bruun rule are abandoned, the formula S sea−level−rise = ξ/tan(β) will remain. Indeed, other conceptual models have been developed and finally came up with the same formula (Davidson-Arnott, 2005). In this paper, we leave aside the key question of the relevance of Equation (1) as a modeling tool to evaluate future shoreline change. Instead, considering that the Bruun rule will continue to be widely used in the future, we evaluate the type of results and the uncertainties that can be expected by using it.
To evaluate the expected impacts of present day and future sea-level rise to sandy beaches erosion, Stive (2004) considered typical values for each term in Equation (1). These values are shown in Tables 1, 2. They are based on observations in the Netherlands and Australia. Stive (2004) showed that under present sea-level rise rates, the contribution of the Bruun effect to sediment losses and shoreline changes is of the same order of magnitude or lower than other effects. Because of the lack of longterm [O (10 years)] coastal data, firmly (in)validating the Bruun rule is difficult (e.g., Leatherman et al., 2000a,b;Sallenger et al., 2000), with the impacts of present-day sea-level rise on shoreline changes being challenging to observe (Stive, 2004;Le Cozannet et al., 2014). However, Stive (2004) showed that for higher rates of sea-level rise, the Bruun effect will not be neglectable any more, as it will significantly impact shoreline change. Hence, Stive (2004) not only helps understanding the behavior of Equation (1) in most general cases, but also provides a general background consistent with present-day observations. The periods of time by which one will be able to (in)validate the use of Bruun rule in Equation (1) is a critical unknown. The question is complex, because all terms in Equation (1) are uncertain, and because the duration by which sea-levelrise-induced erosion is decipherable likely depends on climate change scenario and regional coastal settings. Ultimately, this effect might never be observed in some regions, if the impacts of sea-level rise remain smaller than those of other sedimentary processes causing shoreline change. To investigate this question, we consider the case of idealized wave-exposed sandy beaches with infinite sand availability, and we adapt an approach developed to quantify uncertainties in future flooding occurrence (Le Cozannet et al., 2015): we first define realistic probability density function for each uncertain parameter in Equation (1), using values from Stive (2004) and the IPCC (Church et al., 2013) (Section 2). In Section 3, we propagate these uncertainties through Equation (1) to provide shoreline change projections for different idealized coastal settings and climate change scenarios. Using a global sensitivity analysis (Sobol', 2001;Saltelli et al., 2008), we evaluate the contribution of each uncertain parameter to the variance of future shoreline changes. We use these results to assess where and when the Bruun effect should become observable. Finally, in Section 4, we discuss the results with respect to previous work and question the representativeness of the idealized sandy shorelines considered.

MATERIALS AND METHODS
The method proceeds in three steps, which are detailed below.

Modeling Uncertainties of Input Parameters
We model the uncertainties of each input variable in Equation (1) using the probability density functions indicated in Table 3 and the values shown in Tables 1, 2. These tables consider different types of idealized coastal sites, based on the values provided by Stive (2004) for representative beaches in Australia and the Netherlands. For example, Table 2 indicates that if a groin extending far offshore is implemented on a high energy beach, it will act as a trap for the large longshore sediment transport and will contribute to shoreline change at a rate soundly exceeding 10 m/year. These orders of magnitudes can have a large site to site variability. This point is further discussed in the Section 5. Overall, in the absence of any other contribution similar to that of Stive (2004), we rely on the large uncertainties related to the contribution of each driving process shown in Tables 1, 2 to argue that they are applicable for a large number of sandy coasts. When defining probability distributions representing the uncertainties of input parameters (Table 3), the information available is often limited. As a general guidance, it should be avoided to introduce arbitrary constraints besides what is known already. In other words, the statistical entropy of the selected probability distribution functions should be maximized (Mishra, 2002). For example, a uniform probability density function should be selected when only boundaries are known. This case is met for most of the parameters given in Table 3, including the beach slopes and the contributions of cross-shore and longshore processes to the sedimentary budget. This choice reflects the fact that the data provided by Stive (2004) define boundary values for each uncertain input parameter (see Table 1).
More complex distribution can be elaborated when sufficient information is available. This is the case here for future sea-level rise. This source of uncertainties increases with time, and several components must be distinguished: (1) the selection of a climate change scenario, (2) the uncertainties of future sea-level rise for each climate change scenario, and (3) the regional variability of sea-level rise.
The IPCC provides likely range and median values for four climate change scenarios (RCP 2.6, 4.5, 6.0, and 8.5), corresponding to different trajectories of greenhouse gas emissions. To select one of these sets of value, we use a uniform discrete probability distribution, that is, a random application with a probability of 1/4 to select each climate change scenario.
This choice assumes that we do not know what will be future greenhouse gas emissions.
Even if the true climate change scenario was known, future sea-level rise would still remain uncertain, in particular because of insufficient knowledge on future ice-sheets melting and the Earth energy imbalance (Church et al., 2013;Kopp et al., 2014;von Schuckmann et al., 2016). Modeling results combined with expert judgements have suggested that this source of uncertainties can be adequately modeled by a skewed distribution with a finite support (Ben Abdallah et al., 2014;Jevrejeva et al., 2014). As in our previous article (Le Cozannet et al., 2015), we use a Beta distribution to represent this source of uncertainties. While many scientists estimate that the actual rates sea-level rise over the 21st century could exceed those anticipated by the IPCC , we strongly rely on their likely range and median values to elaborate our sea-level rise scenarios (Church et al., 2013). The upper bound of the distribution is a critical unknown, due to uncertainties in the dynamics of polar icesheets melting, which is currently accelerating (Rignot et al., 2011). Referring to recent modeling results of the Antarctic icesheet instability (Ritz et al., 2015), we use a conservative estimate of 1.5 m for the upper bound of sea-level rise by 2100. To our knowledge, there is presently no evidence that the lowprobability and high impact event of an ice-sheet instability can be avoided, even for the climate change scenario RCP 2.6. Therefore, we use the same upper bound of 1.5m by 2100 for all climate change scenarii. This results in probabilistic sea-level rise scenarios as indicated in Figure 1.
Due to spatially non-uniform warming of the ocean and the response of the Earth to ice-sheets melting, sea-level rise will display regional variability (Slangen et al., 2014a;Carson et al., 2016). This introduces an additional source of uncertainties. We assume that it includes a large part of natural randomness, which we evaluate using Church et al. (2013) and Carson et al. (2016). These references suggest that in addition to the boundaries of the distribution, its most likely value is known. Hence, the principle of maximum entropy leads to the choice of a triangular distribution (Mishra, 2002).
Extreme events are generated here assuming their occurrence follows a Poisson distribution. While this assumption is commonly used (Ranasinghe et al., 2012;Wainwright et al., 2015), it is challenged by some storm sequences, such as the winter 2013-2014 events in western Europe (Masselink et al., 2016). Moreover, as we use a Poisson distribution with fixed parameters, only the interannual variability of storms is captured, and the natural randomness is assumed unaffected by climate change (Plant et al., 2014). This hypothesis could be justified, as it is still difficult to assess where and to which extent climate change will affect extreme winds and waves (Planton et al., 2008;Mori et al., 2010;Hemer et al., 2013). However, we consider here that the return periods of storms are uncertain as well. The range of values considered for this uncertain input parameter is sufficiently large to exceed any plausible change of storminess due to climate change and any other limitation related to the Poisson assumption (see Table 1). Finally, we assume that the shoreline retreats after each extreme event ranges from 0 to 20 m. The first value can correspond to a storm occurring at  (2004); see Table 1 Aeolian processes Uniform distribution Same as in Stive (2004); see Table 1 Other cross-shore effects Uniform distribution Same as in Stive (2004); see Table 1 Longshore processes Uniform distribution Same as in Stive (2004); see Table 2 Exposure to marine forcing Uniform discrete distribution Low or high energy beaches Note that human impacts refer here to the presence of groins, which are coastal engineering structures acting as a trap for the longshore sediment transport to locally accrete beaches.
low tide. Field surveys shows that the second value can result from a collision regime after exceptional storms, and can even be exceeded in case of overwash or breaching, which are not considered in this study. Figure 2 shows examples of shoreline retreats observed in western France after the storm Xynthia in February 2010 (Pedreros et al., 2010;Garcin et al., 2011). The return period of water levels for this storm is estimated to be approximately 200 year in the harbor of La Rochelle (Bulteau et al., 2015). Such return periods are consistent with our working hypothesis, and we note that shoreline retreat ranged from a few meters to 20 m, except where breach occurred. Similar values can be found in many other studies (e.g., Forbes et al., 2004;Mendoza and Jiménez, 2006;Roelvink et al., 2009;Loureiro et al., 2014). As a summary, our approach relies on two major assumptions: (1) sea-level rise scenarios are based on IPCC; (2) we consider sites with gently sloping upper shorefaces, i.e., which are the most impacted by the Bruun effect. Other assumptions can be mentioned: some input parameters here are assumed independent, although they are not in practice. For example, for the same sediment size, beach slopes are expected to be steeper for low-energy beaches (Wright and Short, 1984). However, given the range of uncertainties provided in Table 3, such limitations can be left aside in a first attempt to propagate uncertainties through Equation (1).

Propagation of Uncertainties Through the Sediment Balance Equation
The second step of the approach consists in propagating these uncertainties through Equation (1).
The computational approach starts with defining auxiliary variables with uniform values in the interval [0,1]. Then we apply inverse cumulative distribution functions to these auxiliary variables, in order to generate distributions with the desired shapes (e.g., Beta distribution for sea-level rise).
This computational approach ensures the independence of the different input parameters considered in Table 3, so that the Sobol' parameters (see next section) can be computed using a classical procedure. To reduce the computation time, we use a quasi-Monte-Carlo approach, where we assign values following a quasi-random Sobol' sequence (Sobol', 1967) to the auxiliary variables with values in [0,1]. We empirically determine the number of computations required to converge, by comparing purely random against quasi-random simulations in a test case corresponding to a situation of the end of the 21st century.
This procedure allows to produce generic probabilistic shoreline change projections based on the coastal tract assumptions. Importantly, all input parameters vary simultaneously during these simulations. Therefore, the probability density function of future shoreline change projections integrates the complete variability.

Quantifying and Ranking the Sources of Uncertainties in Shoreline Change Projections
Once uncertainties of all input parameters have been propagated through Equation (1), we obtain a distribution for future shoreline change as a function of time. To understand what this distribution means, we analyse the contribution of each uncertain input parameter to the variance of the model outcome using a global sensitivity analysis (Sobol', 2001;Saltelli et al., 2008).
The principle of variance-based global sensitivity analysis is to decompose the variance of the model outcome Y into several terms that relate to each uncertain input parameter X i . These terms are called the Sobol' indices. Two subsets of terms are of particular interest for the evaluation of uncertainties: the first-oder and total order Sobol' indices.   (Bulteau et al., 2015), and is therefore compatible with the random sample of storms considered in this study.
Frontiers in Marine Science | www.frontiersin.org 5 April 2016 | Volume 3 | Article 49 The first-order Sobol' indices (S i ) correspond to the contribution of a given input parameter alone to the variance of the model outcome: They represent the expected proportion of the variance of the model outcome that would be removed if X i was known. The total order Sobol' index (S Ti ) represents the contribution of a given input parameter and all its interactions with other parameters to the variance of the model outcome. It is defined as: with X −i the set of all X j except X i . This index is used to identify which parameters can be set to any possible value without much impact to the variance of Y. All computations are done using R (R Core Team, 2014), using the same quasi Monte-Carlo approach as in the previous section. The values of the Sobol' indices are very similar for two consecutive years. Therefore, to reduce the computing time, we calculate these indices each 5 year and interpolate them over the 21st century. We use the codes developped by the Joint Research Center of the European Commission to compute the first and total order Sobol' indices (see Jansen, 1999;Saltelli et al., 2010, for the related algorithms and formulations). Figure 3 shows the probabilistic shoreline change projections S(t) at year t, once uncertainties from the input parameters (Table 3) have been propagated through Equation (1). This represents S(t) if no prior information is available, besides values in Tables 1-3. Figures 4-6 show S(t) for high and low energy beaches, with and without groins, and considering the different climate change scenarios separately.

RESULTS: WHEN WILL WE BE ABLE TO (IN)VALIDATE THE BRUUN RULE?
As expected, a slight shift toward erosion is observed over the 21st century (Figure 3). However, in practice, the response is very different depending on the climate change scenario and the type of coast considered. In case of RCP 2.6, we can not identify any clear separation between present and future probability density functions of shoreline change rates. Conversely, for other climate change scenarios and for beaches without groins, the likely range of future shoreline change rates separates or nearly separates (case of RCP 4.5) from present days values (Figures 4, 5). For example, in case of RCP 8.5 and coasts without groins, it is more than likely that observations will separate from their present values by the end of the 2060s for low-energy beaches, and in the early 2070s for high-energy beaches. Hence, using Figures 4, 5, it becomes possible to identify times of emergence for a more than likely observable shift toward erosion of sandy beaches under the assumptions listed above. As a practical implication, it should be possible to assess the Bruun rule validity by the mid-21st century on beaches without groins with gently sloping shoreface, if efforts to reduce greenhouse gas concentrations in the atmosphere fail.
If groins are present on the beach, the longshore processes have a larger impact on the sedimentary budget ( Table 2): on the one hand, large volumes of sediments can be trapped in protected areas, but on the other hand, these sediments are missing in adjacent locations. This does not prevent from observing a shift toward erosion in case of low-energy beaches and climate change scenarios RCP 4.5, 6.0, and 8.5 (Figure 4). In these cases, not enough sediments can be trapped to compensate the losses induced by sea-level rise. Conversely, if groins are built on high-energy beaches, the transport of sediment is severely altered. This results in high variability of shoreline advance and retreat, reaching ±10 m/year (see Stive, 2004, and Table 2), which may compensate the impacts of sea-level rise theoretically (Figure 6). However, it is widely acknowledged that groins are not a sustainable solution because they typically fix the erosion issue, but shift it nearby. Hence, it should not be concluded that groins are an efficient adaptation strategy for high-energy sandy coasts.
Drawing on these results, present day and future rates of shoreline change could remain indistinguishable if greenhouse gas emissions are reduced importantly (RCP 2.6 scenario). In addition, for coastal sites similar to those presented by Stive (2004), the times of emergence of a shift toward erosion induced FIGURE 3 | Shoreline retreat rate projections for the 21st century, with no prior information regarding the type of sandy coast considered, the impacts of man-made infrastructures and the climate change scenario. Positive values corresponds to shoreline retreat. For each date, the intensity of the color represents the probability density function, and the solid, dashed and dotted black lines represent the median, likely, and very likely ranges. The gray lines represent present days rates: they are extended toward the whole 21st century for comparison. They represent a virtual case in which sea-level remains at its present level, whereas the black lines take into account uncertain sea-level rise scenarios of Figure 1 Table 2). The meaning of gray and black lines is the same as in Figure 3.
by sea-level rise should occur during the second half of the 21st century, if sea-level rise follows the IPCC projections.  (Bird, 1985). This suggests that the uncertainties in shoreline change projections can be reduced importantly with more knowledge regarding local coastal processes. Hence, the ensuing question is to understand what drives the variability of shoreline change projections in Figure 3. In this section, we investigate this issue using the results of the global sensitivity analysis. In other words, we separate the variance of S(t) into several components corresponding to the contribution of the 12 uncertain input parameters, in order to gain further insight into the understanding of the upper shoreface sediment balance equation outcome under Table 3 constraints. Figures 7, 8 show the first and total order Sobol' indices obtained by the global sensitivity analysis. Using these figures, it becomes possible to classify input parameters according to their contribution to the variability of shoreline change rates: the first order Sobol' indices are used to rank input parameters, whereas the variability of parameters with total order Sobol' indices close to zero can be neglected without much impacts to the variability of the final results (Saltelli, 2004). In the following subsections, the uncertain input parameters of Equation (1) are classified according to their first and total order Sobol' indices, which reveal research priorities to better anticipate future shoreline change under the Bruun rule assumption.

Critical Unknowns for the Present Days
The random variables used to model the interannual variability of storms and the impacts of longshore sedimentary processes display large first and total order indices: Figure 7 shows that Frontiers in Marine Science | www.frontiersin.org 7 April 2016 | Volume 3 | Article 49 FIGURE 5 | Shoreline change rate projections for the 21st century for low-energy sandy coasts with no impacts of man-made infrastructures and for different IPCC climate change scenarios. The meaning of gray and black lines is the same as in Figure 3.
the value of the first order indices is of about 0.2. Figure 7 indicates that the total order indices decreases from about 0.6-0.3 for longshore processes (0.4-0.3 for the interannual variability of storms). Therefore, these parameters account for a large part of the uncertainties in shoreline change projections. As a consequence, knowing more about these processes should be a priority for reducing the uncertainties resulting from the application of the sediment balance equation for sandy upper shorefaces with gentle slopes. For example, our results implies that a perfect knowledge of longshore processes would reduce the variance of shoreline change projections by a factor of nearly 20%. Many studies in the field of coastal research have studied these sources of uncertainties (Allen, 1981;Inman and Dolan, 1989;Slott et al., 2006;Roelvink et al., 2009;Castelle et al., 2015). However, it remains difficult to quantify precisely all the phenomena involved (see Section 5). Here, the results of the global sensitivity analysis just highlight that sandy shoreline change projections are more accurate if the impacts of storms and longshore sedimentary processes are taken into account. While integrating these processes at regional to global scales is a challenge, it is increasingly addressed in recent assessments of climate change impacts (e.g., Anderson et al., 2015).

Uncertain Parameters with Minor Impacts
At the opposite, the Sobol' first and total order indices representing the aeolian processes, other cross-shore effects (e.g., due to waves asymmetry) and the regional variability of sealevel rise are small. Therefore, the variability of these parameters can be neglected without any large impact to the variability of shoreline change projections. This result can seem surprising. In fact, several studies have highlighted the importance of cross-shore sedimentary processes in explaining seasonal and interannual shoreline changes (e.g., Yates et al., 2009;Splinter et al., 2014). To understand how our results can be related to these findings, we need to point out that Stive (2004) distinguishes the rapid erosive impacts of highly energetic events from the slow accretion due to the net onshore contributions from mild waves ( Table 1).
The regional variability of sea-level rise also plays a minor role, whereas it can be much more important for strongly non-linear impact models. For example, urban coastal flooding generally becomes most damaging once a critical threshold has been exceeded (Hallegatte et al., 2013;Idier et al., 2013b;Miller et al., 2014). In this case, knowing more about the regional variability of sea-level rise would be one of the priorities for coastal managers interested in understanding how risk will evolve over the coming decades (Le Cozannet et al., 2015). According to Table 1, aeolian processes taking place on coastal beaches and dunes are expected to induce net sediment losses and shoreline retreat. However, these processes do not appear as a critical unknown either, in spite of their complexity (Arens, 1996;Hesp et al., 2013;Bauer et al., 2015). Again, our results highlight that when few data are available regarding other coastal processes, the values provided in Table 1 are sufficiently narrow to anticipate the impacts of the slow recovery of beaches satisfactorily. Finally, as this result comes directly from the order of magnitude provided in Tables 1-3, it can be invalidated if new observations demonstrate that these processes have much larger impacts. This question of the representativeness of the values of Tables 1, 2 is addressed in details in the discussion section.

Interactions between Processes: Why is a Global Sensitivity Analysis Recommended?
Figures 7, 8 also reveal that for two random variables, the Sobol' first order indices are close to zero, whereas the total order indices are large. This is the case for example for the random variables indicating if groins are built on the beach, or if the beach is exposed to low-or high-energetic conditions: their total order Sobol' indices ranges from 0.1 to 0.2 (Figure 8). In this case, the principles of the global sensitivity analysis implies that the variability of the input parameters cannot be neglected, because it interacts with other random variables (Saltelli, 2004). In other words, this result indicates that, for instance, groins will have little impacts on shoreline change if longshore drift is weak. Figure 7 shows that the sum of the first order Sobol' indices ranges from 0.45 to 0.6 depending on the period of time considered. This means that the interaction term is large. Hence, to a significant extent, the overall variability is driven by the combined variation of the parameters (joint effect). A very common way to perform a sensitivity analysis consists in evaluating the model response to each uncertain input factor independently. Our results show that despite the apparent simplicity of Equation (1), such an approach could fail to rank the impact of each uncertainty source: the total uncertainty on shoreline change can not be considered as the sum related to each uncertainty sources taken alone. On the contrary, the combined effect among them should be accounted for: here the interaction terms represent about the half of the total variance of shoreline change projections.
FIGURE 6 | Shoreline change rate projections for the 21st century for high-energy sandy coasts with groins and for different IPCC climate change scenarios. The meaning of gray and black lines is the same as in Figure 3.
Frontiers in Marine Science | www.frontiersin.org 9 April 2016 | Volume 3 | Article 49 FIGURE 7 | Evolution of the first order Sobol' indices over the time, indicating the contribution of each varying parameter alone to the final uncertainty. Parameters indicated in gray have the smallest 1st order Sobol' indices. The inflexion of the curves after 2060 is due to two phenomena: first, the different climate change scenarios (RCP 2.6-8.5) induce different sea-level rise projections by that time (see Figure 1); second, shoreline change projections uncertainties by 2100 is much larger by 2100 than now (see text). The colors (red, green, blue, and yellow) indicate uncertainties related to (1) climate change and sea-level rise, (2) local site characteristics, (3) impacts of storms, and (4) of coastal defenses (groins).

Critical Unknowns for the Next Decades
Figures 7, 8 shows two distinct periods, during which the relative importance of sea-level rise, climate change scenarios, and beach slopes increase successively. This result primarily comes from the propagation of uncertain sea-level rise scenarios (Figure 1) in the Bruun term of Equation (1). First, the climate change scenarios have small Sobol' indices during the 1st part of the 21st century. Therefore, as expected, the next decades do not appear as a relevant period of time to differentiate future shoreline change according to the future climate change scenarios. However, sea-level will be rising, and Equation (1) implies that this will induce a slight shift toward erosion. By 2050, uncertainties on global sea-level rise account for 20% of the variance of future shoreline change projections.
Conversely, there are large differences between e.g., RCP 2.6 and RCP 8.5 during the second part of the 21st century. Figure 7 shows that if a decision is taken to follow a given climate change scenario, the variance of shoreline change projections is reduced by 20%. This is almost twice the maximum reduction of variance that we could expect, if we were able to predict accurately future sea-level rise given greenhouse gas emissions by 2100. Hence, even if nothing is known regarding coastal processes besides Table 1, the benefits of climate change mitigation and greenhouse gaz emissions reductions appear significant, as suggested by Nicholls and Lowe (2004) and Pardaens et al. (2011). Here, Figures 7, 8 constitute an other way to communicate this result, which was already identified in Section 3 and Figure 3.  (1) will be straightforward by the end of the 21st century for RCP 6.0 and 8.5 scenarios (see Section 3). In a more general case, the results of the global sensitivity analysis (Figures 7,  8) enable to identify which types of coastal sites are potential candidates for testing the use of the Bruun term in Equation (1): beaches without human impacts, little exposed to storms and where longshore drift can be quantified are among such sites. This result is not a surprise: previous studies that tried to validate the Bruun rule and detect an observable impact of sea-level rise on shoreline change actually focused precisely on this type of coasts (Leatherman et al., 2000b;Zhang et al., 2004). In practice however, the selection of such sites is not straightforward, and there remains question regarding the possibility to remove the impacts of longshore drift in the analysis (Leatherman et al., 2000a;Sallenger et al., 2000). Nevertheless, these results provide a perspective for coastal observatories currently collecting coastal evolution observations with the aim of better understanding coastal impacts of climate change.

DISCUSSION: REPRESENTATIVENESS OF THE SELECTED IDEALIZED COASTAL SITES
Our study relies on the orders of magnitudes provided by Stive (2004) for typical coastal settings in the Netherlands and Australia. Though there is no study providing orders of magnitude at other representative sites, it is hypothesized that the respective contribution of each process in Tables 1, 2 vary considerably for other coastal settings. For instance, on open and straight sandy beaches exposed to high-energy oblique ocean waves, longshore sediment transport rate can be high [e.g., O(10 6 m 3 /year)] but without any alongshore gradients. In such a situation, longshore sedimentary processes do not contribute to shoreline change. This is the case of the southern Gironde Coast, SW France, where gradients in longshore sediment transport are negligible (Idier et al., 2013a) and where the observed shoreline variability on the timescales from hours to years can be explained by cross-shore processes only ). Yet, even in this simplified coastal setting, there is to date no detailed quantification of the contribution of aeolian processes and other cross-shore processes to the overall shoreline change. As far as extreme storms are concerned, it is impossible to quantify maximum storm wave erosion for given return period exceeding a few decades simply because there is no such long time series of storm-driven erosion. At the same beach, the cumulative impact of a series of outstanding storms during the winter 2013/2014 only drove a 15 m dune retreat over 3 months (Castelle et al., 2015). The same outstanding severe storms drove highly variable erosion patterns along Western Europe (Masselink et al., 2016), revealing the complexity of quantifying the erosion driven by multi-decadal or multi-centennial return period storms.
While providing numbers for this simple case of open beaches with no alongshore sediment transport gradients is arguably challenging, this is even worst for embayed beaches. Embayed beaches are ubiquitous worldwide given that nearly half of the world's coast consists of hilly or mountainous coastline. In addition, beaches artificially bounded by groins or jetty can behave as embayed beaches. Shoreline variability along embayed beaches is dominated by rotation signal (Ranasinghe et al., 2004). However, the respective contributions of cross-shore and longshore processes to the embayed beach rotation signal is still a subject of debate (Ranasinghe et al., 2004;Harley et al., 2011Harley et al., , 2015. Overall, although we do not pretend to address all sandy coasts worldwide, we advocate that the values given by Stive (2004) are sufficiently large to be representative of many coastal settings. It is obvious that further work is required (1) to decipher the respective contributions of the different driving processes to shoreline change in different coastal settings and (2) to develop alternative model framework in which the different processes can be accounted for explicitly. In this framework, reducedcomplexity modeling approaches (e.g., Splinter et al., 2014;Reeve et al., 2015) appear as a relevant avenue.
Noteworthy, Stive (2004) refers to coastal sites with gentle slopes, ranging from 0.001 to 0.002, which corresponds to very gentle slopes actually (Hinton and Nicholls, 1998;Marsh et al., 1998). Shoreface slopes are often larger. For instance in SW France, high-energy beaches display coastal slopes of about 0.015 , whereas they reach about 0.012 for low-energy Mediterranean beaches in Languedoc (Yates et al., 2011). In practice, however, the use of the Bruun rule in Equation (1) is more questionable for steep coasts. Other models would be required in such coastal sites, for example considering the sediment losses and gains during and after each storm (Larson et al., 2004;Ranasinghe et al., 2012). Nevertheless, we acknowledge that beach slopes values that can be calculated from Table 3 correspond to coastal sites where the Bruun effect will have the strongest impacts.
The approach presented here could be applied to other medium-complexity geomorphic evolution modeling approaches (French et al., 2015), in order to finally identify research priorities. However, in the case of complex coastal settings with limited sand availability, the constraints imposed by the geology, the continental sedimentary supply or the role of coastal lagoons should be considered as well.

CONCLUSION
In this article, we examined where and by which period of time the Bruun effect should become observable on waveexposed sandy coasts. To do this, we refered to Stive (2004) to evaluate the order of magnitude of each physical process contributing to shoreline change in the case of coastal uppershoreface with gentle slopes. We also considered probabilistic sea-level rise scenarios based on the IPCC (Church et al., 2013). We propagated these uncertainties through the sediment balance equation, that sums the Bruun effect with other drivers such as longshore and cross-shore processes. This allowed to identify a time of emergence of an observable Bruun effect by the middle of the 21st century in the case of beaches with gentle slopes without groins, and for climate change scenarios RCP 4.5, 6.0, and 8.5 only. Under the assumptions above, the absence of generalized sea-level-rise-driven beach erosion by the end of the 21st century is only achieved for RCP 2.6 climate change scenario.
Rather than predictions of future shoreline changes, our results can be viewed as an attempt to better understand where and when the Bruun rule can be (in)validated. Using a global sensitivity analysis, our results confirm that low-energy gently sloping beaches with little human impacts and small gradients in longshore drift and sheltered from storms are the most relevant to assess the validity of the Bruun rule in the sediment balance equation. These results could be further improved through a better integration of the primary driving processes of shoreline change in the sediment balance equation, and refined probability distribution functions. This also implies collecting site-specific and more accurate values than those listed in Tables 1, 2. Whatever the amount of data collected, part of the uncertainties in future shoreline change projections will remain. In addition, it should be noted that new knowledge may also illuminate some "not-yet envisaged" complex processes, of which our understanding may still be poor, and would require new research efforts in the future. Overall, and regardless the shortcomings discussed above, our approach is a relevant modeling framework to further improve communication to stakeholders and general public on uncertainties in sea-level rise impact assessments.

AUTHOR CONTRIBUTIONS
Overall design of the study: GC. Detailed design and implementation of the study: GC, CO, BC, MG, DI, RP, JR. Interpretation and discussion: All. Writing and correcting the paper: All.

ACKNOWLEDGMENTS
This project has been funded by BRGM through the "Adaptation to Climate Change Risks" project. BC was funded by the CHIPO project (ANR-14-ASTR-0004-01), financially supported by the Agence Nationale de la Recherche (ANR). The field mission following Xynthia (Figure 2) was funded by BRGM, and conducted by MG and RP, Daniel Monfort, Julie Mugica, Yann Krien, and Benjamin Francois. We thank Thomas Bulteau, Anny Cazenave, Franck Desmazes, Franck Lavigne, Cesca Ribas, and Marcel Stive for useful discussions and comments that finally resulted in this study. Finally, we thank the two reviewers for their very useful comments.